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.

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.

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.

Expansion and improvement of ChinaMu by MuT-seq and chromosome-level assembly of the Mu-starter genome.

ChinaMu is the largest sequence-indexed Mutator (Mu) transposon insertional library in maize (Zea mays). In this study, we made significant improvements to the size and quality of the ChinaMu library. We developed a new Mu-tag isolation method Mu-Tn5-seq (MuT-seq). Compared to the previous method used by ChinaMu, MuT-seq recovered 1/3 more germinal insertions, while requiring only about 1/14 of the sequencing volume and 1/5 of the experimental time. Using MuT-seq, we identified 113,879 germinal insertions from 3,168 Mu-active F1 families. We also assembled a high-quality genome for the Mu-active line Mu-starter, which harbors the initial active MuDR element and was used as the pollen donor for the mutation population. Using the Mu-starter genome, we recovered 33,662 (15.6%) additional germinal insertions in 3,244 (7.4%) genes in the Mu-starter line. The Mu-starter genome also improved the assignment of 117,689 (54.5%) germinal insertions. The newly upgraded ChinaMu dataset currently contains 215,889 high-quality germinal insertions. These insertions cover 32,224 pan-genes in the Mu-starter and B73Ref5 genomes, including 23,006 (80.4%) core genes shared by the two genomes. As a test model, we investigated Mu insertions in the pentatricopeptide repeat (PPR) superfamily, discovering insertions for 92% (449/487) of PPR genes in ChinaMu, demonstrating the usefulness of ChinaMu as a functional genomics resource for maize.

DREB1A/CBF3 Is Repressed by Transgene-Induced DNA Methylation in the Arabidopsis ice1 -1 Mutant.

DREB1/CBFs are key transcription factors involved in plant cold stress adaptation. The expression of DREB1/CBFs triggers a cold-responsive transcriptional cascade, after which many stress tolerance genes are expressed. Thus, elucidating the mechanisms of cold stress-inducible DREB1/CBF expression is important to understand the molecular mechanisms of plant cold stress responses and tolerance. We analyzed the roles of a transcription factor, INDUCER OF CBF EXPRESSION1 (ICE1), that is well known as an important transcriptional activator in the cold-inducible expression of DREB1A/CBF3 in Arabidopsis (Arabidopsis thaliana). ice1-1 is a widely accepted mutant allele known to abolish cold-inducible DREB1A expression, and this evidence has strongly supported ICE1-DREB1A regulation for many years. However, in ice1-1 outcross descendants, we unexpectedly discovered that ice1-1 DREB1A repression was genetically independent of the ice1-1 allele ICE1(R236H). Moreover, neither ICE1 overexpression nor double loss-of-function mutation of ICE1 and its homolog SCRM2 altered DREB1A expression. Instead, a transgene locus harboring a reporter gene in the ice1-1 genome was responsible for altering DREB1A expression. The DREB1A promoter was hypermethylated due to the transgene. We showed that DREB1A repression in ice1-1 results from transgene-induced silencing and not genetic regulation by ICE1. The ICE1(R236H) mutation has also been reported as scrm-D, which confers constitutive stomatal differentiation. The scrm-D phenotype and the expression of a stomatal differentiation marker gene were confirmed to be linked to the ICE1(R236H) mutation. We propose that the current ICE1-DREB1 regulatory model should be revalidated without the previous assumptions.

Cognate restriction of transposition by piggyBac-like proteins.

Mobile genetic elements have been harnessed for gene transfer for a wide variety of applications including generation of stable cell lines, recombinant protein production, creation of transgenic animals, and engineering cell and gene therapy products. The piggyBac transposon family includes transposase or transposase-like proteins from a variety of species including insect, bat and human. Recently, human piggyBac transposable element derived 5 (PGBD5) protein was reported to be able to transpose piggyBac transposons in human cells raising possible safety concerns for piggyBac-mediated gene transfer applications. We evaluated three piggyBac-like proteins across species including piggyBac (insect), piggyBat (bat) and PGBD5 (human) for their ability to mobilize piggyBac transposons in human cells. We observed a lack of cross-species transposition activity. piggyBac and piggyBat activity was restricted to their cognate transposons. PGBD5 was unable to mobilize piggyBac transposons based on excision, colony count and plasmid rescue analysis, and it was unable to bind piggyBac terminal repeats. Within the piggyBac family, we observed a lack of cross-species activity and found that PGBD5 was unable to bind, excise or integrate piggyBac transposons in human cells. Transposition activity appears restricted within species within the piggyBac family of mobile genetic elements.