DNA replication is a highly coordinated cell cycle process that can become dysregulated in cancer, increasing both proliferation and mutation rates. Single-cell whole genome sequencing holds potential for studying replication dynamics of cancer cells; however, computational methods for identifying S-phase cells and inferring single-cell replication timing profiles remain immature for samples with heterogeneous copy number. Here we report a new method, PERT, which jointly infers replication and somatic copy number states of S-phase cells. This method enabled us to analyze the replication dynamics of >10,000 S-phase single-cell genomes across various triple negative breast cancers and cell lines with subclonal copy number heterogeneity. We show that PERT robustly predicts cell cycle phase, quantifies replication timing variability, and approximates relative proliferation rates between tumor subclones. Our results illuminate how aberrant DNA replication processes can both drive and result from evolution of human tumors.

DNA strand exchange, the central step of homologous recombination, is considered to occur approximately independently of DNA sequence content. However, certain prokaryotic and eukaryotic genomic loci display either an enhanced or reduced frequency of genetic exchange. Here we show that the Homo sapiens DNA strand exchange protein, HsRad51, shows a preference for binding to single-stranded DNA sequences primarily rich in G-residues and poor in A- and C-residues, and that these DNA sequences manifest enhanced HsRad51 protein-dependent homologous pairing. Both of these properties are common to all DNA strand exchange proteins examined thus far. These preferred DNA pairing sequences resemble those found at genetic loci in human cells that cause genomic instability and lead to genetic diseases.

Double strand break (DSB) repair mechanisms guard genome integrity and their deterioration causes genomic instability. Common and rare fragile sites (CFS and RFS, respectively) are particularly vulnerable to instability, and there is an inverse correlation between fragile site (FS) expression and DSB repair protein levels. Upon DSB repair dysfunction, genes residing at these sites are at greater risk of deregulation compared to genes located at non-FS. In this regard, it remains enigmatic why the incidence of miRNA genes at FS is higher compared to non-FS. Herein, using bioinformatics, we examined whether miRNA genes localized at FS inhibit components of DSB repair pathways and assessed their effects on cancer. We show that such miRNAs over-accumulate in RFS, and that FRAXA, which is expressed in Fragile X syndrome, is a conserved hotspot for miRNAs inhibiting DSB repair. Axes of FRAXA-residing miRNAs/DSB repair targets affect survival in a cancer type-specific manner. Moreover, copy number variations in the region encompassing these miRNA genes discriminate survival between male and female patients. Given that, thus far, only CFS have been considered relevant for carcinogenesis, our data are the first to associate RFS with cancer, through the impairment of DSB repair by the FRAXA-residing miRNAs.

Hypermutable simple sequence repeats (SSRs) are major drivers of phase variation in Campylobacter jejuni. The presence of multiple SSR-mediated phase-variable genes encoding enzymes that modify surface structures, including capsular polysaccharide (CPS) and lipooligosaccharide (LOS), generates extreme cell surface diversity within bacterial populations, thereby promoting adaptation to selective pressures in host environments. Therefore, genetically controlling SSR-mediated phase variation can be important for achieving stable and reproducible research on C. jejuni. Here, we show that natural "cotransformation" is an effective method for C. jejuni genome editing. Cotransformation is a trait of naturally competent bacteria that causes uptake/integration of multiple different DNA molecules, which has been recently adapted to multiplex genome editing by natural transformation (MuGENT), a method for introducing multiple mutations into the genomes of these bacteria. We found that cotransformation efficiently occurred in C. jejuni. To examine the feasibility of MuGENT in C. jejuni, we "locked" different polyG SSR tracts in strain NCTC11168 (which are located in the biosynthetic CPS/LOS gene clusters) into either the ON or OFF configurations. This approach, termed "MuGENT-SSR," enabled the generation of all eight edits within 2 weeks and the identification of a phase-locked strain with a highly stable type of Penner serotyping, a CPS-based serotyping scheme. Furthermore, extensive genome editing of this strain by MuGENT-SSR identified a phase-variable gene that determines the Penner serotype of NCTC11168. Thus, MuGENT-SSR provides a platform for genetic and phenotypic engineering of genetically unstable C. jejuni, making it a reliable approach for elucidating the mechanisms underlying phase-variable expression of specific phenotypes. IMPORTANCE Campylobacter jejuni is the leading bacterial cause of foodborne gastroenteritis in developed countries and occasionally progresses to the autoimmune disease Guillain-Barré syndrome. A relatively large number of hypermutable simple sequence repeat (SSR) tracts in the C. jejuni genome markedly decreases its phenotypic stability through reversible changes in the ON or OFF expression states of the genes in which they reside, a phenomenon called phase variation. Thus, controlling SSR-mediated phase variation can be important for achieving stable and reproducible research on C. jejuni. In this study, we developed a feasible and effective approach for genetically manipulate multiple SSR tracts in the C. jejuni genome using natural cotransformation, a trait of naturally transformable bacterial species that causes the uptake and integration of multiple different DNA molecules. This approach will greatly help to improve the genetic and phenotypic stability of C. jejuni to enable diverse applications in research and development.

Accurate and efficient microbial diagnosis is crucial for effective molecular diagnostics, especially in the field of human healthcare. The gold standard equipment widely employed for detecting specific microorganisms in molecular diagnosis is quantitative real-time polymerase chain reaction (qRT-PCR). However, its limitations in low metagenomic DNA yield samples necessitate exploring alternative approaches. Digital PCR, by quantifying the number of copies of the target sequence, provides absolute quantification results for the bacterial strain. In this study, we compared the diagnostic efficiency of qRT-PCR and digital PCR in detecting a particular bacterial strain (Staphylococcus aureus), focusing on skin-derived DNA samples. Experimentally, specific primer for S. aureus were designed at transcription elongation factor (greA) gene and the target amplicon were cloned and sequenced to validate efficiency of specificity to the greA gene of S. aureus. To quantify the absolute amount of microorganisms present on the skin, the variable region 5 (V5) of the 16S rRNA gene was used, and primers for S. aureus identification were used to relative their amount in the subject's skin. The findings demonstrate the absolute convenience and efficiency of digital PCR in microbial diagnostics. We suggest that the high sensitivity and precise quantification provided by digital PCR could be a promising tool for detecting specific microorganisms, especially in skin-derived DNA samples with low metagenomic DNA yields, and that further research and implementation is needed to improve medical practice and diagnosis.

Cryptic unstable transcripts (CUTs) are a largely unexplored class of nuclear exosome degraded, non-coding RNAs in budding yeast. It is highly debated whether CUT transcription has a functional role in the cell or whether CUTs represent noise in the yeast transcriptome. We sought to ascertain the extent of conserved CUT expression across a variety of Saccharomyces yeast strains to further understand and characterize the nature of CUT expression.

Protein filaments play important roles in many biological processes. We discovered an actin homolog in halophilic archaea, which we call Salactin. Just like the filaments that segregate DNA in eukaryotes, Salactin grows out of the cell poles towards the middle, and then quickly depolymerizes, a behavior known as dynamic instability. Furthermore, we see that Salactin affects the distribution of DNA in daughter cells when cells are grown in low-phosphate media, suggesting Salactin filaments might be involved in segregating DNA when the cell has only a few copies of the chromosome.

Canine Lafora disease is a recessively inherited, rapidly progressing neurodegenerative disease caused by the accumulation of abnormally constructed insoluble glycogen Lafora bodies in the brain and other tissues due to the loss of NHL repeat containing E3 ubiquitin protein ligase 1 (NHLRC1). Dogs have a dodecamer repeat sequence within the NHLRC1 gene, which is prone to unstable (dynamic) expansion and loss of function. Progressive signs of Lafora disease include hypnic jerks, reflex and spontaneous myoclonus, seizures, vision loss, ataxia and decreased cognitive function. We studied five dogs (one Chihuahua, two French Bulldogs, one Griffon Bruxellois, one mixed breed) with clinical signs associated with canine Lafora disease. Identification of polyglucosan bodies (Lafora bodies) in myocytes supported diagnosis in the French Bulldogs; muscle areas close to the myotendinous junction and the myofascial union segment had the highest yield of inclusions. Postmortem examination of one of the French Bulldogs revealed brain Lafora bodies. Genetic testing for the known canine NHLRC1 mutation confirmed the presence of a homozygous mutation associated with canine Lafora disease. Our results show that Lafora disease extends beyond previous known breeds to the French Bulldog, Griffon Bruxellois and even mixed-breed dogs, emphasizing the likely species-wide nature of this genetic problem. It also establishes these breeds as animal models for the devastating human disease. Genetic testing should be used when designing breeding strategies to determine the frequency of the NHLRC1 mutation in affected breeds. Lafora diseases should be suspected in any older dog presenting with myoclonus, hypnic jerks or photoconvulsions.

An A-to-G transition at position 16247 in the human mtDNA genome denotes haplogroup B4a1a1a and its sublineages. Informally known as the "Polynesian motif," this haplogroup has been widely used as a marker in Oceania of genetic affiliation with the Austronesian expansion. The 16247G allele has arisen only once in the human mtDNA phylogeny, about 7,000 thousand years ago, and is nearly fixed in Remote Oceania. We analyzed 536 complete mtDNA genome sequences from the Solomon Islands from haplogroup B4a1a1 and associated subhaplogroups and found multiple independent back mutations from 16247G to 16247A. We also find elevated levels of heteroplasmy at this position in samples with the 16247G allele, suggesting the ongoing occurrence of somatic back-mutations and/or transmission of heteroplasmy. Moreover, the G allele is predicted to introduce a novel stem-loop structure in the DNA sequence that may be structurally unfavorable, thereby accounting for the remarkable number of back-mutations observed at the 16247G allele in this short evolutionary time span. More generally, haplogroup-calling scripts result in inaccurate haplogroup calls involving the back-mutation and need to be supplemented with other types of analyses; this may be true for other mtDNA lineages because no other lineage has been investigated to the same extent (over 500 complete mtDNA sequences).

Genome instability, associated with chromosome breakage syndromes and most human cancers, is still poorly understood. In the yeast Saccharomyces cerevisiae, numerous genes with roles in the preservation of genome integrity have been identified. DNA-damage-checkpoint-deficient yeast cells that lack Sgs1, a RecQ-like DNA helicase related to the human Bloom's-syndrome-associated helicase BLM, show an increased rate of genome instability, and we have previously shown that they accumulate recurring chromosomal translocations between three similar genes, CAN1, LYP1 and ALP1. Here, the chromosomal location, copy number and sequence similarity of the translocation targets ALP1 and LYP1 were altered to gain insight into the formation of complex translocations. Among 844 clones with chromosomal rearrangements, 93 with various types of simple and complex translocations involving CAN1, LYP1 and ALP1 were identified. Breakpoint sequencing and mapping showed that the formation of complex translocation types is strictly dependent on the location of the initiating DNA break and revealed that complex translocations arise via a combination of interchromosomal translocation and template-switching, as well as from unstable dicentric intermediates. Template-switching occurred between sequences on the same chromosome, but was inhibited if the genes were transferred to different chromosomes. Unstable dicentric translocations continuously gave rise to clones with multiple translocations in various combinations, reminiscent of intratumor heterogeneity in human cancers. Base substitutions and evidence of DNA slippage near rearrangement breakpoints revealed that translocation formation can be accompanied by point mutations, and their presence in different translocation types within the same clone provides evidence that some of the different translocation types are derived from each other rather than being formed de novo. These findings provide insight into eukaryotic genome instability, especially the formation of translocations and the sources of intraclonal heterogeneity, both of which are often associated with human cancers.

The considerable genome size variation in Arabidopsis thaliana has been shown largely to be due to copy number variation (CNV) in 45S ribosomal RNA (rRNA) genes. Surprisingly, attempts to map this variation by means of genome-wide association studies (GWAS) failed to identify either of the two likely sources, namely the nucleolus organizer regions (NORs). Instead, GWAS implicated a trans-acting locus, as if rRNA gene CNV was a phenotype rather than a genotype. To explain these results, we investigated the inheritance and stability of rRNA gene copy number using the variety of genetic resources available in A. thaliana - F2 crosses, recombinant inbred lines, the multiparent advanced-generation inter-cross population, and mutation accumulation lines. Our results clearly show that rRNA gene CNV can be mapped to the NORs themselves, with both loci contributing equally to the variation. However, NOR size is unstably inherited, and dramatic copy number changes are visible already within tens of generations, which explains why it is not possible to map the NORs using GWAS. We did not find any evidence of trans-acting loci in crosses, which is also expected since changes due to such loci would take very many generations to manifest themselves. rRNA gene copy number is thus an interesting example of "missing heritability"-a trait that is heritable in pedigrees, but not in the general population.

Animals have stable dominant gut microbiomes under similar diets. Similar diets can also lead to similar gut microbial communities within host species levels. Giant pandas (Ailuropoda melanoleuca) and red pandas (Ailurus fulgens) have had long-term and stable bamboo diets, and seem well adapted to this highly fibrous diet. When compared to the gut microbiomes of Père David's deer (Elaphurus davidianus), humans, cheetah (Acinonyx jubatus), black-backed jackal (Canis-mesomelas), and black bear (Ursus thibetanus), giant panda gut microbiomes have high variation in the abundance of Pseudomonadaceae and Clostridiaceae, and are somewhat unstable. This high instability and dissimilarity may reflect an unstable gut environment, perturbation or selective pressure because of their carnivorous gastrointestinal system. A short digestive tract, brief digestion time and fast intestinal peristalsis may result in higher oxygen concentrations that select for the growth of aerobes and facultative anaerobes in giant pandas. Potential selection of high proportion of Pseudomonadaceae in giant panda (GP-HP) and red panda gut microbiomes may arise because of their postulated ability to degrade secondary compounds (e.g., cyanide compounds and aromatic compounds). However, high proportion of Clostridiaceae (GP-HF) may focus on cellulose and hemicellulose digestion. Thus, GP-HP and GP-HF groups have high dissimilarity on the functional level. These findings show that long-term similarities in diet do not always lead to similar or stable gut microbial system within the same host species and that other factors can drive the selection of gut taxa.

Colorectal cancer develops through two main genetic instability pathways characterized by distinct pathologic features and clinical outcome.

More than 30 human disorders are caused by the expansion of simple sequence DNA repeats, among which triplet repeats remain the most frequent. Most trinucleotide repeat expansion disorders affect primarily the nervous system, through mechanisms of neurodysfunction and/or neurodegeneration. While trinucleotide repeat tracts are short and stably transmitted in unaffected individuals, disease-associated expansions are highly dynamic in the germline and in somatic cells, with a tendency toward further expansion. Since longer repeats are associated with increasing disease severity and earlier onset of symptoms, intergenerational repeat size gains account for the phenomenon of anticipation. In turn, higher levels of age-dependent somatic expansion have been linked with increased disease severity and earlier age of onset, implicating somatic instability in the onset and progression of disease symptoms. Hence, tackling the root cause of symptoms through the control of repeat dynamics may provide therapeutic modulation of clinical manifestations. DNA repair pathways have been firmly implicated in the molecular mechanism of repeat length mutation. The demonstration that repeat expansion depends on functional DNA mismatch repair (MMR) proteins, points to MMR as a potential therapeutic target. Similarly, a role of DNA base excision repair (BER) in repeat expansion has also been suggested, particularly during the removal of oxidative lesions. Using a well-characterized mouse cell model system of an unstable CAG•CTG trinucleotide repeat, we tested if expanded repeat tracts can be stabilized by small molecules with reported roles in both pathways: cadmium (an inhibitor of MMR activity) and a variety of antioxidants (capable of neutralizing oxidative species). We found that chronic exposure to sublethal doses of cadmium and antioxidants did not result in significant reduction of the rate of trinucleotide repeat expansion. Surprisingly, manganese yielded a significant stabilization of the triplet repeat tract. We conclude that treatment with cadmium and antioxidants, at doses that do not interfere with cell survival and cell culture dynamics, is not sufficient to modify trinucleotide repeat dynamics in cell culture.

Microsatellite instability-high (MSI-H) tumors are characterized by high tumor mutation burden and responsiveness to checkpoint blockade. We identified tumor-specific frameshifts encoding multiple epitopes that originated from indel mutations shared among patients with MSI-H endometrial, colorectal, and stomach cancers. Epitopes derived from these shared frameshifts have high population occurrence rates, wide presence in many tumor subclones, and are predicted to bind to the most frequent MHC alleles in MSI-H patient cohorts. Neoantigens arising from these mutations are distinctly unlike self and viral antigens, signifying novel groups of potentially highly immunogenic tumor antigens. We further confirmed the immunogenicity of frameshift peptides in T cell stimulation experiments using blood mononuclear cells isolated from both healthy donors and MSI-H cancer patients. Our study uncovers the widespread occurrence and strong immunogenicity of tumor-specific antigens derived from shared frameshift mutations in MSI-H cancer and Lynch syndrome patients, suitable for the design of common "off-the-shelf" cancer vaccines.

Purification of proteins that participate in large transient complexes is impeded by low amounts, heterogeneity, instability and poor solubility. To circumvent these difficulties we set up a methodology that enables the production of stable complexes for structural and functional studies. This procedure is benchmarked and applied to two challenging protein families: the human steroid nuclear receptors (SNR) and the HIV-1 pre-integration complex. In the context of transcriptional regulation studies, we produce and characterize the ligand-binding domains of the glucocorticoid nuclear receptor and the oestrogen receptor beta in complex with a TIF2 (transcriptional intermediary factor 2) domain containing the three SNR-binding motifs. In the context of retroviral integration, we demonstrate the stabilization of the HIV-1 integrase by formation of complexes with partner proteins and DNA. This procedure provides a powerful research tool for structural and functional studies of proteins participating in non-covalent macromolecular complexes.

Metastatic cancer affects millions of people worldwide annually and is the leading cause of cancer-related deaths. Most patients with metastatic disease are not eligible for surgical resection, and current therapeutic regimens have varying success rates, some with 5-year survival rates below 5%. Here we test the hypothesis that metastatic cancer can be genetically targeted by exploiting single base substitution mutations unique to individual cells that occur as part of normal aging prior to transformation. These mutations are targetable because ~10% of them form novel tumor-specific "NGG" protospacer adjacent motif (PAM) sites targetable by CRISPR-Cas9.

Many eukaryotic membrane proteins are prone to misfolding, which compromises their functional expression at the plasma membrane. This is particularly true for mammalian gonadotropin-releasing hormone receptors (GnRHRs), which are G protein-coupled receptors involved in reproductive steroidogenesis. We recently demonstrated that evolutionary modifications within mammalian GnRHRs appear to have coincided with adaptive changes in cotranslational folding efficiency. Though changes in protein stability are known to shape evolutionary interactions, it is unclear how the energetic drivers of cotranslational folding in the membrane may modify epistatic interactions. We therefore surveyed the pairwise epistatic interactions that modify the expression of two destabilized GnRHR variants bearing mutations that selectively compromise either its membrane topology (V276T) or its native tertiary structure (W107A). Using deep mutational scanning (DMS), we evaluated how the effects of these mutations on the expression of the mature form of the protein at the plasma membrane are modified by hundreds of secondary mutations. A focused analysis of 251 mutants with high-quality measurements in three genetic backgrounds reveals that V276T and W107A form distinct epistatic interactions that depend on both the degree to which they destabilize the protein and the mechanism of their destabilization. An unsupervised learning analysis shows that V276T forms predominantly negative epistatic interactions that are most pronounced among destabilizing mutations within soluble loop regions. In contrast, W107A forms interactions with mutations in both loops and transmembrane domains that skew positive as a result of the diminishing impact of the destabilizing mutations in the context of an already unstable variant. These findings provide general insights into how pairwise epistasis is remodeled by conformational defects in membrane proteins and, more generally, in unstable proteins.

Tumor progression and evolution are frequently associated with chromosomal instability (CIN). Tumor cells often express high levels of the mitotic checkpoint protein MAD2, leading to mitotic arrest and cell death. However, some tumor cells are capable of exiting mitosis and consequently increasing CIN. How cells escape the mitotic arrest induced by MAD2 and proliferate with CIN is not well understood. Here, we explored loss-of-function screens and drug sensitivity tests associated with MAD2 levels in aneuploid cells and identified that aneuploid cells with high MAD2 levels are more sensitive to FOXM1 depletion. Inhibition of FOXM1 promotes MAD2-mediated mitotic arrest and exacerbates CIN. Conversely, elevating FOXM1 expression in MAD2-overexpressing human cell lines reverts prolonged mitosis and rescues mitotic errors, cell death and proliferative disadvantages. Mechanistically, we found that FOXM1 facilitates mitotic exit by inhibiting the spindle assembly checkpoint (SAC) and the expression of Cyclin B. Notably, we observed that FOXM1 is upregulated upon aneuploid induction in cells with dysfunctional SAC and error-prone mitosis, and these cells are sensitive to FOXM1 knockdown, indicating a novel vulnerability of aneuploid cells.

Host-virus association data underpin research into the distribution and eco-evolutionary correlates of viral diversity and zoonotic risk across host species. However, current knowledge of the wildlife virome is inherently constrained by historical discovery effort, and there are concerns that the reliability of ecological inference from host-virus data may be undermined by taxonomic and geographical sampling biases. Here, we evaluate whether current estimates of host-level viral diversity in wild mammals are stable enough to be considered biologically meaningful, by analysing a comprehensive dataset of discovery dates of 6571 unique mammal host-virus associations between 1930 and 2018. We show that virus discovery rates in mammal hosts are either constant or accelerating, with little evidence of declines towards viral richness asymptotes, even in highly sampled hosts. Consequently, inference of relative viral richness across host species has been unstable over time, particularly in bats, where intensified surveillance since the early 2000s caused a rapid rearrangement of species' ranked viral richness. Our results illustrate that comparative inference of host-level virus diversity across mammals is highly sensitive to even short-term changes in sampling effort. We advise caution to avoid overinterpreting patterns in current data, since it is feasible that an analysis conducted today could draw quite different conclusions than one conducted only a decade ago.