Genetic manipulation of gut microbes enables single-gene interrogation in a complex microbiome.

Hundreds of microbiota genes are associated with host biology/disease. Unraveling the causal contribution of a microbiota gene to host biology remains difficult because many are encoded by nonmodel gut commensals and not genetically targetable. A general approach to identify their gene transfer methodology and build their gene manipulation tools would enable mechanistic dissections of their impact on host physiology. We developed a pipeline that identifies the gene transfer methods for multiple nonmodel microbes spanning five phyla, and we demonstrated the utility of their genetic tools by modulating microbiome-derived short-chain fatty acids and bile acids in vitro and in the host. In a proof-of-principle study, by deleting a commensal gene for bile acid synthesis in a complex microbiome, we discovered an intriguing role of this gene in regulating colon inflammation. This technology will enable genetically engineering the nonmodel gut microbiome and facilitate mechanistic dissection of microbiota-host interactions.

De novo deletions and duplications at recombination hotspots in mouse germlines.

Numerous DNA double-strand breaks (DSBs) arise during meiosis to initiate homologous recombination. These DSBs are usually repaired faithfully, but here, we uncover a distinct type of mutational event in which deletions form via joining of ends from two closely spaced DSBs (double cuts) within a single hotspot or at adjacent hotspots on the same or different chromatids. Deletions occur in normal meiosis but are much more frequent when DSB formation is dysregulated in the absence of the ATM kinase. Events between chromosome homologs point to multi-chromatid damage and aborted gap repair. Some deletions contain DNA from other hotspots, indicating that double cutting at distant sites creates substrates for insertional mutagenesis. End joining at double cuts can also yield tandem duplications or extrachromosomal circles. Our findings highlight the importance of DSB regulation and reveal a previously hidden potential for meiotic mutagenesis that is likely to affect human health and genome evolution.

Genome-wide de novo L1 Retrotransposition Connects Endonuclease Activity with Replication.

L1 retrotransposon-derived sequences comprise approximately 17% of the human genome. Darwinian selective pressures alter L1 genomic distributions during evolution, confounding the ability to determine initial L1 integration preferences. Here, we generated high-confidence datasets of greater than 88,000 engineered L1 insertions in human cell lines that act as proxies for cells that accommodate retrotransposition in vivo. Comparing these insertions to a null model, in which L1 endonuclease activity is the sole determinant dictating L1 integration preferences, demonstrated that L1 insertions are not significantly enriched in genes, transcribed regions, or open chromatin. By comparison, we provide compelling evidence that the L1 endonuclease disproportionately cleaves predominant lagging strand DNA replication templates, while lagging strand 3'-hydroxyl groups may prime endonuclease-independent L1 retrotransposition in a Fanconi anemia cell line. Thus, acquisition of an endonuclease domain, in conjunction with the ability to integrate into replicating DNA, allowed L1 to become an autonomous, interspersed retrotransposon.

A New Way to Diversify Antibodies by DNA Transposition.

While searching for new therapeutics against malaria, Lanzavecchia and colleagues discovered that antibodies can be diversified by DNA sequences encoded outside of antibody genes.

Bridge RNAs direct programmable recombination of target and donor DNA.

Genomic rearrangements, encompassing mutational changes in the genome such as insertions, deletions or inversions, are essential for genetic diversity. These rearrangements are typically orchestrated by enzymes that are involved in fundamental DNA repair processes, such as homologous recombination, or in the transposition of foreign genetic material by viruses and mobile genetic elements1,2. Here we report that IS110 insertion sequences, a family of minimal and autonomous mobile genetic elements, express a structured non-coding RNA that binds specifically to their encoded recombinase. This bridge RNA contains two internal loops encoding nucleotide stretches that base-pair with the target DNA and the donor DNA, which is the IS110 element itself. We demonstrate that the target-binding and donor-binding loops can be independently reprogrammed to direct sequence-specific recombination between two DNA molecules. This modularity enables the insertion of DNA into genomic target sites, as well as programmable DNA excision and inversion. The IS110 bridge recombination system expands the diversity of nucleic-acid-guided systems beyond CRISPR and RNA interference, offering a unified mechanism for the three fundamental DNA rearrangements-insertion, excision and inversion-that are required for genome design.

Transposase-assisted target-site integration for efficient plant genome engineering.

The current technologies to place new DNA into specific locations in plant genomes are low frequency and error-prone, and this inefficiency hampers genome-editing approaches to develop improved crops1,2. Often considered to be genome 'parasites', transposable elements (TEs) evolved to insert their DNA seamlessly into genomes3-5. Eukaryotic TEs select their site of insertion based on preferences for chromatin contexts, which differ for each TE type6-9. Here we developed a genome engineering tool that controls the TE insertion site and cargo delivered, taking advantage of the natural ability of the TE to precisely excise and insert into the genome. Inspired by CRISPR-associated transposases that target transposition in a programmable manner in bacteria10-12, we fused the rice Pong transposase protein to the Cas9 or Cas12a programmable nucleases. We demonstrated sequence-specific targeted insertion (guided by the CRISPR gRNA) of enhancer elements, an open reading frame and a gene expression cassette into the genome of the model plant Arabidopsis. We then translated this system into soybean-a major global crop in need of targeted insertion technology. We have engineered a TE 'parasite' into a usable and accessible toolkit that enables the sequence-specific targeting of custom DNA into plant genomes.

Chemical-genetic interrogation of RNA polymerase mutants reveals structure-function relationships and physiological tradeoffs.

The multi-subunit bacterial RNA polymerase (RNAP) and its associated regulators carry out transcription and integrate myriad regulatory signals. Numerous studies have interrogated RNAP mechanism, and RNAP mutations drive Escherichia coli adaptation to many health- and industry-relevant environments, yet a paucity of systematic analyses hampers our understanding of the fitness trade-offs from altering RNAP function. Here, we conduct a chemical-genetic analysis of a library of RNAP mutants. We discover phenotypes for non-essential insertions, show that clustering mutant phenotypes increases their predictive power for drawing functional inferences, and demonstrate that some RNA polymerase mutants both decrease average cell length and prevent killing by cell-wall targeting antibiotics. Our findings demonstrate that RNAP chemical-genetic interactions provide a general platform for interrogating structure-function relationships in vivo and for identifying physiological trade-offs of mutations, including those relevant for disease and biotechnology. This strategy should have broad utility for illuminating the role of other important protein complexes.

Alu insertion-mediated dsRNA structure formation with pre-existing Alu elements as a disease-causing mechanism.

We previously identified a homozygous Alu insertion variant (Alu_Ins) in the 3'-untranslated region (3'-UTR) of SPINK1 as the cause of severe infantile isolated exocrine pancreatic insufficiency. Although we established that Alu_Ins leads to the complete loss of SPINK1 mRNA expression, the precise mechanisms remained elusive. Here, we aimed to elucidate these mechanisms through a hypothesis-driven approach. Initially, we speculated that, owing to its particular location, Alu_Ins could independently disrupt mRNA 3' end formation and/or affect other post-transcriptional processes such as nuclear export and translation. However, employing a 3'-UTR luciferase reporter assay, Alu_Ins was found to result in only an ∼50% reduction in luciferase activity compared to wild type, which is insufficient to account for the severe pancreatic deficiency in the Alu_Ins homozygote. We then postulated that double-stranded RNA (dsRNA) structures formed between Alu elements, an upstream mechanism regulating gene expression, might be responsible. Using RepeatMasker, we identified two Alu elements within SPINK1's third intron, both oriented oppositely to Alu_Ins. Through RNAfold predictions and full-length gene expression assays, we investigated orientation-dependent interactions between these Alu repeats. We provide compelling evidence to link the detrimental effect of Alu_Ins to extensive dsRNA structures formed between Alu_Ins and pre-existing intronic Alu sequences, including the restoration of SPINK1 mRNA expression by aligning all three Alu elements in the same orientation. Given the widespread presence of Alu elements in the human genome and the potential for new Alu insertions at almost any locus, our findings have important implications for detecting and interpreting Alu insertions in disease genes.

Common single-base insertions in the VNTR of the carboxyl ester lipase (CEL) gene are benign and also likely to arise somatically in the exocrine pancreas.

The CEL gene encodes carboxyl ester lipase, a pancreatic digestive enzyme. CEL is extremely polymorphic due to a variable number tandem repeat (VNTR) located in the last exon. Single-base deletions within this VNTR cause the inherited disorder MODY8, whereas little is known about VNTR single-base insertions in pancreatic disease. We therefore mapped CEL insertion variants (CEL-INS) in 200 Norwegian patients with pancreatic neoplastic disorders. Twenty-eight samples (14.0%) carried CEL-INS alleles. Most common were insertions in repeat 9 (9.5%), which always associated with a VNTR length of 13 repeats. The combined INS allele frequency (0.078) was similar to that observed in a control material of 416 subjects (0.075). We performed functional testing in HEK293T cells of a set of CEL-INS variants, in which the insertion site varied from the first to the 12th VNTR repeat. Lipase activity showed little difference among the variants. However, CEL-INS variants with insertions occurring in the most proximal repeats led to protein aggregation and endoplasmic reticulum stress, which upregulated the unfolded protein response. Moreover, by using a CEL-INS-specific antibody, we observed patchy signals in pancreatic tissue from humans without any CEL-INS variant in the germline. Similar pancreatic staining was seen in knock-in mice expressing the most common human CEL VNTR with 16 repeats. CEL-INS proteins may therefore be constantly produced from somatic events in the normal pancreatic parenchyma. This observation along with the high population frequency of CEL-INS alleles strongly suggests that these variants are benign, with a possible exception for insertions in VNTR repeats 1-4.

The evolutionary history of 2,658 cancers.

Cancer develops through a process of somatic evolution1,2. Sequencing data from a single biopsy represent a snapshot of this process that can reveal the timing of specific genomic aberrations and the changing influence of mutational processes3. Here, by whole-genome sequencing analysis of 2,658 cancers as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA)4, we reconstruct the life history and evolution of mutational processes and driver mutation sequences of 38 types of cancer. Early oncogenesis is characterized by mutations in a constrained set of driver genes, and specific copy number gains, such as trisomy 7 in glioblastoma and isochromosome 17q in medulloblastoma. The mutational spectrum changes significantly throughout tumour evolution in 40% of samples. A nearly fourfold diversification of driver genes and increased genomic instability are features of later stages. Copy number alterations often occur in mitotic crises, and lead to simultaneous gains of chromosomal segments. Timing analyses suggest that driver mutations often precede diagnosis by many years, if not decades. Together, these results determine the evolutionary trajectories of cancer, and highlight opportunities for early cancer detection.

An insertion in the promoter of a malate dehydrogenase gene regulates malic acid content in apple fruit.

Malic acid is an important flavor determinant in apple (Malus × domestica Borkh.) fruit. One known variation controlling malic acid is the A/G single nucleotide polymorphism in an aluminum-activated malate transporter gene (MdMa1). Nevertheless, there are still differences in malic acid content in apple varieties with the same Ma1 genotype (Ma1/Ma1 homozygous), such as 'Honeycrisp' (high malic acid content) and 'Qinguan' (low malic acid content), indicating that other loci may influence malic acid and fruit acidity. Here, the F1 (Filial 1) hybrid generation of 'Honeycrisp' × 'Qinguan' was used to analyze quantitative trait loci for malic acid content. A major locus (Ma7) was identified on chromosome 13. Within this locus, a malate dehydrogenase gene, MDH1 (MdMa7), was the best candidate for further study. Subcellular localization suggested that MdMa7 encodes a cytosolic protein. Overexpression and RNA interference of MdMa7 in apple fruit increased and decreased malic acid content, respectively. An insertion/deletion (indel) in the MdMa7 promoter was found to affect MdMa7 expression and malic acid content in both hybrids and other cultivated varieties. The insertion and deletion genotypes were designated as MA7 and ma7, respectively. The transcription factor MdbHLH74 was found to stimulate MdMa7 expression in the MA7 genotype but not in the ma7 genotype. Transient transformation of fruit showed that MdbHLH74 affected MdMa7 expression and malic acid content in 'Gala' (MA7/MA7) but not in 'Fuji' (ma7/ma7). Our results indicated that genetic variation in the MdMa7 (MDH1) promoter alters the binding ability of the transcription factor MdbHLH74, which alters MdMa7 (MDH1) transcription and the malic acid content in apple fruit, especially in Ma1/Ma1 homozygous accessions.

Long-read sequencing identifies an SVA_D retrotransposon insertion deep within the intron of ATP7A as a novel cause of occipital horn syndrome.

BACKGROUND: SINE-VNTR-Alu (SVA) retrotransposons move from one genomic location to another in a 'copy-and-paste' manner. They continue to move actively and cause monogenic diseases through various mechanisms. Currently, disease-causing SVA retrotransposons are classified into human-specific young SVA_E or SVA_F subfamilies. In this study, we identified an evolutionarily old SVA_D retrotransposon as a novel cause of occipital horn syndrome (OHS). OHS is an X-linked, copper metabolism disorder caused by dysfunction of the copper transporter, ATP7A. METHODS: We investigated a 16-year-old boy with OHS whose pathogenic variant could not be detected via routine molecular genetic analyses. RESULTS: A 2.8 kb insertion was detected deep within the intron of the patient's ATP7A gene. This insertion caused aberrant mRNA splicing activated by a new donor splice site located within it. Long-read circular consensus sequencing enabled us to accurately read the entire insertion sequence, which contained highly repetitive and GC-rich segments. Consequently, the insertion was identified as an SVA_D retrotransposon. Antisense oligonucleotides (AOs) targeting the new splice site restored the expression of normal transcripts and functional ATP7A proteins. AO treatment alleviated excessive accumulation of copper in patient fibroblasts in a dose-dependent manner. Pedigree analysis revealed that the retrotransposon had moved into the OHS-causing position two generations ago. CONCLUSION: This is the first report of a human monogenic disease caused by the SVA_D retrotransposon. The fact that the evolutionarily old SVA_D is still actively transposed, leading to increased copy numbers may make a notable impact on rare genetic disease research.

DNA Repair Pathway Choices in CRISPR-Cas9-Mediated Genome Editing.

Many clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated protein 9 (Cas9)-based genome editing technologies take advantage of Cas nucleases to induce DNA double-strand breaks (DSBs) at desired locations within a genome. Further processing of the DSBs by the cellular DSB repair machinery is then necessary to introduce desired mutations, sequence insertions, or gene deletions. Thus, the accuracy and efficiency of genome editing are influenced by the cellular DSB repair pathways. DSBs are themselves highly genotoxic lesions and as such cells have evolved multiple mechanisms for their repair. These repair pathways include homologous recombination (HR), classical nonhomologous end joining (cNHEJ), microhomology-mediated end joining (MMEJ) and single-strand annealing (SSA). In this review, we briefly highlight CRISPR-Cas9 and then describe the mechanisms of DSB repair. Finally, we summarize recent findings of factors that can influence the choice of DNA repair pathway in response to Cas9-induced DSBs.

Deletion of the s2m RNA Structure in the Avian Coronavirus Infectious Bronchitis Virus and Human Astrovirus Results in Sequence Insertions.

Coronaviruses infect a wide variety of host species, resulting in a range of diseases in both humans and animals. The coronavirus genome consists of a large positive-sense single-stranded molecule of RNA containing many RNA structures. One structure, denoted s2m and consisting of 41 nucleotides, is located within the 3' untranslated region (3' UTR) and is shared between some coronavirus species, including infectious bronchitis virus (IBV), severe acute respiratory syndrome coronavirus (SARS-CoV), and SARS-CoV-2, as well as other pathogens, including human astrovirus. Using a reverse genetic system to generate recombinant viruses, we investigated the requirement of the s2m structure in the replication of IBV, a globally distributed economically important Gammacoronavirus that infects poultry causing respiratory disease. Deletion of three nucleotides predicted to destabilize the canonical structure of the s2m or the deletion of the nucleotides corresponding to s2m impacted viral replication in vitro. In vitro passaging of the recombinant IBV with the s2m sequence deleted resulted in a 36-nucleotide insertion in place of the deletion, which was identified to be composed of a duplication of flanking sequences. A similar result was observed following serial passage of human astrovirus with a deleted s2m sequence. RNA modeling indicated that deletion of the nucleotides corresponding to the s2m impacted other RNA structures present in the IBV 3' UTR. Our results indicated for both IBV and human astrovirus a preference for nucleotide occupation in the genome location corresponding to the s2m, which is independent of the specific s2m sequence. IMPORTANCE Coronaviruses infect many species, including humans and animals, with substantial effects on livestock, particularly with respect to poultry. The coronavirus RNA genome consists of structural elements involved in viral replication whose roles are poorly understood. We investigated the requirement of the RNA structural element s2m in the replication of the Gammacoronavirus infectious bronchitis virus, an economically important viral pathogen of poultry. Using reverse genetics to generate recombinant IBVs with either a disrupted or deleted s2m, we showed that the s2m is not required for viral replication in cell culture; however, replication is decreased in tracheal tissue, suggesting a role for the s2m in the natural host. Passaging of these viruses as well as human astrovirus lacking the s2m sequence demonstrated a preference for nucleotide occupation, independent of the s2m sequence. RNA modeling suggested deletion of the s2m may negatively impact other essential RNA structures.

A self-propagating, barcoded transposon system for the dynamic rewiring of genomic networks.

In bacteria, natural transposon mobilization can drive adaptive genomic rearrangements. Here, we build on this capability and develop an inducible, self-propagating transposon platform for continuous genome-wide mutagenesis and the dynamic rewiring of gene networks in bacteria. We first use the platform to study the impact of transposon functionalization on the evolution of parallel Escherichia coli populations toward diverse carbon source utilization and antibiotic resistance phenotypes. We then develop a modular, combinatorial assembly pipeline for the functionalization of transposons with synthetic or endogenous gene regulatory elements (e.g., inducible promoters) as well as DNA barcodes. We compare parallel evolutions across alternating carbon sources and demonstrate the emergence of inducible, multigenic phenotypes and the ease with which barcoded transposons can be tracked longitudinally to identify the causative rewiring of gene networks. This work establishes a synthetic transposon platform that can be used to optimize strains for industrial and therapeutic applications, for example, by rewiring gene networks to improve growth on diverse feedstocks, as well as help address fundamental questions about the dynamic processes that have sculpted extant gene networks.

Insertion of YFP at P5CS1 and AFL1 shows the potential, and potential complications, of gene tagging for functional analyses of stress-related proteins.

Crispr/CAS9-enabled homologous recombination to insert a tag in frame with an endogenous gene can circumvent difficulties such as context-dependent promoter activity that complicate analysis of gene expression and protein accumulation patterns. However, there have been few reports examining whether such gene targeting/gene tagging (GT) can alter expression of the target gene. The enzyme encoded by Δ1-pyrroline-5-carboxylate synthetase 1 (P5CS1) is key for stress-induced proline synthesis and drought resistance, yet its expression pattern and protein localisation have been difficult to assay. We used GT to insert YFP in frame with the 5' or 3' ends of the endogenous P5CS1 and At14a-Like 1 (AFL1) coding regions. Insertion at the 3' end of either gene generated homozygous lines with expression of the gene-YFP fusion indistinguishable from the wild type allele. However, for P5CS1 this occurred only after selfing and advancement to the T5 generation allowed initial homozygous lethality of the insertion to be overcome. Once this was done, the GT-generated P5CS1-YFP plants revealed new information about P5CS1 localisation and tissue-specific expression. In contrast, insertion of YFP at the 5' end of either gene blocked expression. The results demonstrate that GT can be useful for functional analyses of genes that are problematic to properly express by other means but also show that, in some cases, GT can disrupt expression of the target gene.

Novel Alu insertion in the ZEB2 gene causing Mowat-Wilson syndrome.

Alu elements are short, interspersed elements located throughout the genome, playing a role in human diversity, and occasionally causing genetic diseases. Here, we report a novel Alu insertion causing Mowat-Wilson syndrome, a rare neurodevelopmental disorder, in an 8-year-old boy displaying the typical clinical features for Mowat-Wilson syndrome. The variant was not initially detected in genome sequencing data, but through deep phenotyping, which pointed to only one plausible candidate gene, manual inspection of genome sequencing alignment data enabled us to identify a de novo heterozygous Alu insertion in exon 8 of the ZEB2 gene. Nanopore long-read sequencing confirmed the Alu insertion, leading to the formation of a premature stop codon and likely haploinsufficiency of ZEB2. This underscores the importance of deep phenotyping and mobile element insertion analysis in uncovering genetic causes of monogenic disorders as these elements might be overlooked in standard next-generation sequencing protocols.

Dna2 nuclease deficiency results in large and complex DNA insertions at chromosomal breaks.

Insertions of mobile elements1-4, mitochondrial DNA5 and fragments of nuclear chromosomes6 at DNA double-strand breaks (DSBs) threaten genome integrity and are common in cancer7-9. Insertions of chromosome fragments at V(D)J recombination loci can stimulate antibody diversification10. The origin of insertions of chromosomal fragments and the mechanisms that prevent such insertions remain unknown. Here we reveal a yeast mutant, lacking evolutionarily conserved Dna2 nuclease, that shows frequent insertions of sequences between approximately 0.1 and 1.5 kb in length into DSBs, with many insertions involving multiple joined DNA fragments. Sequencing of around 500 DNA inserts reveals that they originate from Ty retrotransposons (8%), ribosomal DNA (rDNA) (15%) and from throughout the genome, with preference for fragile regions such as origins of replication, R-loops, centromeres, telomeres or replication fork barriers. Inserted fragments are not lost from their original loci and therefore represent duplications. These duplications depend on nonhomologous end-joining (NHEJ) and Pol4. We propose a model in which alternative processing of DNA structures arising in Dna2-deficient cells can result in the release of DNA fragments and their capture at DSBs. Similar DNA insertions at DSBs are expected to occur in any cells with linear extrachromosomal DNA fragments.

Virus-like insertions with sequence signatures similar to those of endogenous nonretroviral RNA viruses in the human genome.

Understanding the genetics and taxonomy of ancient viruses will give us great insights into not only the origin and evolution of viruses but also how viral infections played roles in our evolution. Endogenous viruses are remnants of ancient viral infections and are thought to retain the genetic characteristics of viruses from ancient times. In this study, we used machine learning of endogenous RNA virus sequence signatures to identify viruses in the human genome that have not been detected or are already extinct. Here, we show that the k-mer occurrence of ancient RNA viral sequences remains similar to that of extant RNA viral sequences and can be differentiated from that of other human genome sequences. Furthermore, using this characteristic, we screened RNA viral insertions in the human reference genome and found virus-like insertions with phylogenetic and evolutionary features indicative of an exogenous origin but lacking homology to previously identified sequences. Our analysis indicates that animal genomes still contain unknown virus-derived sequences and provides a glimpse into the diversity of the ancient virosphere.

Conjugative plasmids interact with insertion sequences to shape the horizontal transfer of antimicrobial resistance genes.

It is well established that plasmids play an important role in the dissemination of antimicrobial resistance (AMR) genes; however, little is known about the role of the underlying interactions between different plasmid categories and other mobile genetic elements (MGEs) in shaping the promiscuous spread of AMR genes. Here, we developed a tool designed for plasmid classification, AMR gene annotation, and plasmid visualization and found that most plasmid-borne AMR genes, including those localized on class 1 integrons, are enriched in conjugative plasmids. Notably, we report the discovery and characterization of a massive insertion sequence (IS)-associated AMR gene transfer network (245 combinations covering 59 AMR gene subtypes and 53 ISs) linking conjugative plasmids and phylogenetically distant pathogens, suggesting a general evolutionary mechanism for the horizontal transfer of AMR genes mediated by the interaction between conjugative plasmids and ISs. Moreover, our experimental results confirmed the importance of the observed interactions in aiding the horizontal transfer and expanding the genetic range of AMR genes within complex microbial communities.