All-RNA-mediated targeted gene integration in mammalian cells with rationally engineered R2 retrotransposons.

All-RNA-mediated targeted gene integration methods, rendering reduced immunogenicity, effective deliverability with non-viral vehicles, and a low risk of random mutagenesis, are urgently needed for next-generation gene addition technologies. Naturally occurring R2 retrotransposons hold promise in this context due to their site-specific integration profile. Here, we systematically analyzed the biodiversity of R2 elements and screened several R2 orthologs capable of full-length gene insertion in mammalian cells. Robust R2 system gene integration efficiency was attained using combined donor RNA and protein engineering. Importantly, the all-RNA-delivered engineered R2 system showed effective integration activity, with efficiency over 60% in mouse embryos. Unbiased high-throughput sequencing demonstrated that the engineered R2 system exhibited high on-target integration specificity (99%). In conclusion, our study provides engineered R2 tools for applications based on hit-and-run targeted DNA integration and insights for further optimization of retrotransposon systems.

The Ancient Gamete Fusogen HAP2 Is a Eukaryotic Class II Fusion Protein.

Sexual reproduction is almost universal in eukaryotic life and involves the fusion of male and female haploid gametes into a diploid cell. The sperm-restricted single-pass transmembrane protein HAP2-GCS1 has been postulated to function in membrane merger. Its presence in the major eukaryotic taxa-animals, plants, and protists (including important human pathogens like Plasmodium)-suggests that many eukaryotic organisms share a common gamete fusion mechanism. Here, we report combined bioinformatic, biochemical, mutational, and X-ray crystallographic studies on the unicellular alga Chlamydomonas reinhardtii HAP2 that reveal homology to class II viral membrane fusion proteins. We further show that targeting the segment corresponding to the fusion loop by mutagenesis or by antibodies blocks gamete fusion. These results demonstrate that HAP2 is the gamete fusogen and suggest a mechanism of action akin to viral fusion, indicating a way to block Plasmodium transmission and highlighting the impact of virus-cell genetic exchanges on the evolution of eukaryotic life.

Triaging of α-helical proteins to the mitochondrial outer membrane by distinct chaperone machinery based on substrate topology.

Mitochondrial outer membrane ⍺-helical proteins play critical roles in mitochondrial-cytoplasmic communication, but the rules governing the targeting and insertion of these biophysically diverse proteins remain unknown. Here, we first defined the complement of required mammalian biogenesis machinery through genome-wide CRISPRi screens using topologically distinct membrane proteins. Systematic analysis of nine identified factors across 21 diverse ⍺-helical substrates reveals that these components are organized into distinct targeting pathways that act on substrates based on their topology. NAC is required for the efficient targeting of polytopic proteins, whereas signal-anchored proteins require TTC1, a cytosolic chaperone that physically engages substrates. Biochemical and mutational studies reveal that TTC1 employs a conserved TPR domain and a hydrophobic groove in its C-terminal domain to support substrate solubilization and insertion into mitochondria. Thus, the targeting of diverse mitochondrial membrane proteins is achieved through topological triaging in the cytosol using principles with similarities to ER membrane protein biogenesis systems.

Mechanism of signal-anchor triage during early steps of membrane protein insertion.

Most membrane proteins use their first transmembrane domain, known as a signal anchor (SA), for co-translational targeting to the endoplasmic reticulum (ER) via the signal recognition particle (SRP). The SA then inserts into the membrane using either the Sec61 translocation channel or the ER membrane protein complex (EMC) insertase. How EMC and Sec61 collaborate to ensure SA insertion in the correct topology is not understood. Using site-specific crosslinking, we detect a pre-insertion SA intermediate adjacent to EMC. This intermediate forms after SA release from SRP but before ribosome transfer to Sec61. The polypeptide's N-terminal tail samples a cytosolic vestibule bordered by EMC3, from where it can translocate across the membrane concomitant with SA insertion. The ribosome then docks on Sec61, which has an opportunity to insert those SAs skipped by EMC. These results suggest that EMC acts between SRP and Sec61 to triage SAs for insertion during membrane protein biogenesis.

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.

Transposon Insertion Sequencing, a Global Measure of Gene Function.

The goal of genomics and systems biology is to understand how complex systems of factors assemble into pathways and structures that combine to form living organisms. Great advances in understanding biological processes result from determining the function of individual genes, a process that has classically relied on characterizing single mutations. Advances in DNA sequencing has made available the complete set of genetic instructions for an astonishing and growing number of species. To understand the function of this ever-increasing number of genes, a high-throughput method was developed that in a single experiment can measure the function of genes across the genome of an organism. This occurred approximately 10 years ago, when high-throughput DNA sequencing was combined with advances in transposon-mediated mutagenesis in a method termed transposon insertion sequencing (TIS). In the subsequent years, TIS succeeded in addressing fundamental questions regarding the genes of bacteria, many of which have been shown to play central roles in bacterial infections that result in major human diseases. The field of TIS has matured and resulted in studies of hundreds of species that include significant innovations with a number of transposons. Here, we summarize a number of TIS experiments to provide an understanding of the method and explanation of approaches that are instructive when designing a study. Importantly, we emphasize critical aspects of a TIS experiment and highlight the extension and applicability of TIS into nonbacterial species such as yeast.

Function of the Omp85 Superfamily of Outer Membrane Protein Assembly Factors and Polypeptide Transporters.

The Omp85 protein superfamily is found in the outer membrane (OM) of all gram-negative bacteria and eukaryotic organelles of bacterial origin. Members of the family catalyze both the membrane insertion of ß-barrel proteins and the translocation of proteins across the OM. Although the mechanism(s) by which these proteins function is unclear, striking new insights have emerged from recent biochemical and structural studies. In this review we discuss the entire Omp85 superfamily but focus on the function of the best-studied member, BamA, which is an essential and highly conserved component of the bacterial barrel assembly machinery (BAM). Because BamA has multiple functions that overlap with those of other Omp85 proteins, it is likely the prototypical member of the Omp85 superfamily. Furthermore, BamA has become a protein of great interest because of the recent discovery of small-molecule inhibitors that potentially represent an important new class of antibiotics.

The Landscape of L1 Retrotransposons in the Human Genome Is Shaped by Pre-insertion Sequence Biases and Post-insertion Selection.

L1 retrotransposons are transposable elements and major contributors of genetic variation in humans. Where L1 integrates into the genome can directly impact human evolution and disease. Here, we experimentally induced L1 retrotransposition in cells and mapped integration sites at nucleotide resolution. At local scales, L1 integration is mostly restricted by genome sequence biases and the specificity of the L1 machinery. At regional scales, L1 shows a broad capacity for integration into all chromatin states, in contrast to other known mobile genetic elements. However, integration is influenced by the replication timing of target regions, suggesting a link to host DNA replication. The distribution of new L1 integrations differs from those of preexisting L1 copies, which are significantly reshaped by natural selection. Our findings reveal that the L1 machinery has evolved to efficiently target all genomic regions and underline a predominant role for post-integrative processes on the distribution of endogenous L1 elements.

Scleraxis-lineage cells are required for correct muscle patterning.

Movement of the vertebrate body is supported by the connection of muscle, tendon and bone. Each skeletal muscle in the vertebrate body has a unique shape and attachment site; however, the mechanism that ensures reproducible muscle patterning is incompletely understood. In this study, we conducted targeted cell ablation using scleraxis (Scx)-Cre to examine the role of Scx-lineage cells in muscle morphogenesis and attachment in mouse embryos. We found that muscle bundle shapes and attachment sites were significantly altered in embryos with Scx-lineage cell ablation. Muscles in the forelimb showed impaired bundle separation and limb girdle muscles distally dislocated from their insertion sites. Scx-lineage cells were required for post-fusion myofiber morphology, but not for the initial segregation of myoblasts in the limb bud. Furthermore, muscles could change their attachment site, even after formation of the insertion. Lineage tracing suggested that the muscle patterning defect was primarily attributed to the reduction of tendon/ligament cells. Our study demonstrates an essential role of Scx-lineage cells in the reproducibility of skeletal muscle attachment, in turn revealing a previously unappreciated tissue-tissue interaction in musculoskeletal morphogenesis.

Mechanistic analysis of a novel membrane-interacting variable loop in the pleckstrin-homology domain critical for dynamin function.

Classical dynamins are best understood for their ability to generate vesicles by membrane fission. During clathrin-mediated endocytosis (CME), dynamin is recruited to the membrane through multivalent protein and lipid interactions between its proline-rich domain (PRD) with SRC Homology 3 (SH3) domains in endocytic proteins and its pleckstrin-homology domain (PHD) with membrane lipids. Variable loops (VL) in the PHD bind lipids and partially insert into the membrane thereby anchoring the PHD to the membrane. Recent molecular dynamics (MD) simulations reveal a novel VL4 that interacts with the membrane. Importantly, a missense mutation that reduces VL4 hydrophobicity is linked to an autosomal dominant form of Charcot-Marie-Tooth (CMT) neuropathy. We analyzed the orientation and function of the VL4 to mechanistically link data from simulations with the CMT neuropathy. Structural modeling of PHDs in the cryo-electron microscopy (cryo-EM) cryoEM map of the membrane-bound dynamin polymer confirms VL4 as a membrane-interacting loop. In assays that rely solely on lipid-based membrane recruitment, VL4 mutants with reduced hydrophobicity showed an acute membrane curvature-dependent binding and a catalytic defect in fission. Remarkably, in assays that mimic a physiological multivalent lipid- and protein-based recruitment, VL4 mutants were completely defective in fission across a range of membrane curvatures. Importantly, expression of these mutants in cells inhibited CME, consistent with the autosomal dominant phenotype associated with the CMT neuropathy. Together, our results emphasize the significance of finely tuned lipid and protein interactions for efficient dynamin function.

Digital microfluidics-based digital counting of single-cell copy number variation (dd-scCNV Seq).

Single-cell copy number variations (CNVs), major dynamic changes in humans, result in differential levels of gene expression and account for adaptive traits or underlying disease. Single-cell sequencing is needed to reveal these CNVs but has been hindered by single-cell whole-genome amplification (scWGA) bias, leading to inaccurate gene copy number counting. In addition, most of the current scWGA methods are labor intensive, time-consuming, and expensive with limited wide application. Here, we report a unique single-cell whole-genome library preparation approach based on digital microfluidics for digital counting of single-cell Copy Number Variation (dd-scCNV Seq). dd-scCNV Seq directly fragments the original single-cell DNA and uses these fragments as templates for amplification. These reduplicative fragments can be filtered computationally to generate the original partitioned unique identified fragments, thereby enabling digital counting of copy number variation. dd-scCNV Seq showed an increase in uniformity in the single-molecule data, leading to more accurate CNV patterns compared to other methods with low-depth sequencing. Benefiting from digital microfluidics, dd-scCNV Seq allows automated liquid handling, precise single-cell isolation, and high-efficiency and low-cost genome library preparation. dd-scCNV Seq will accelerate biological discovery by enabling accurate profiling of copy number variations at single-cell resolution.

Insertions and Deletions: Computational Methods, Evolutionary Dynamics, and Biological Applications.

Insertions and deletions constitute the second most important source of natural genomic variation. Insertions and deletions make up to 25% of genomic variants in humans and are involved in complex evolutionary processes including genomic rearrangements, adaptation, and speciation. Recent advances in long-read sequencing technologies allow detailed inference of insertions and deletion variation in species and populations. Yet, despite their importance, evolutionary studies have traditionally ignored or mishandled insertions and deletions due to a lack of comprehensive methodologies and statistical models of insertions and deletion dynamics. Here, we discuss methods for describing insertions and deletion variation and modeling insertions and deletions over evolutionary time. We provide practical advice for tackling insertions and deletions in genomic sequences and illustrate our discussion with examples of insertions and deletion-induced effects in human and other natural populations and their contribution to evolutionary processes. We outline promising directions for future developments in statistical methodologies that would allow researchers to analyze insertions and deletion variation and their effects in large genomic data sets and to incorporate insertions and deletions in evolutionary inference.

Insertion sequence excision is enhanced by a protein that catalyzes branch migration and promotes microhomology-mediated end joining.

Insertion sequence (IS)-excision enhancer (IEE) promotes the excision of ISs in the genome of enterohemorrhagic Escherichia coli O157. Because IEE-dependent IS excision occurs in the presence of transposase, the process of IS transposition may be involved in IS excision; however, little is understood about the molecular mechanisms of IS excision. Our in vitro analysis revealed that IEE exhibits DNA-dependent ATPase activity, which is activated by branched DNA. IEE also catalyzes the branch migration of fork-structured DNA. These results suggest that IEE remodels branched structures of the IS transposition intermediate. Sequence analysis of recombination sites in IS-excision products suggested that microhomologous sequences near the ends of the IS are involved in IS excision. IEE promoted microhomology-mediated end joining (MMEJ), in which base pairing between 6-nucleotides complementary ends of two 3'-protruding DNAs and subsequent elongation of the paired DNA strand occurred. IS-excision frequencies were significantly decreased in cells producing IEE mutants that had lost either branch migration or MMEJ activity, which suggests that these activities of IEE are required for IS excision. Based on our results, we propose a model for IS excision triggered by IEE and transposase.

Experimental and computational approaches for membrane protein insertion and topology determination.

Membrane proteins play pivotal roles in a wide array of cellular processes and constitute approximately a quarter of the protein-coding genes across all organisms. Despite their ubiquity and biological significance, our understanding of these proteins remains notably less comprehensive compared to their soluble counterparts. This disparity in knowledge can be attributed, in part, to the inherent challenges associated with employing specialized techniques for the investigation of membrane protein insertion and topology. This review will center on a discussion of molecular biology methodologies and computational prediction tools designed to elucidate the insertion and topology of helical membrane proteins.

Correlations between alignment gaps and nucleotide substitution or amino acid replacement.

To assess the conventional treatment in evolutionary inference of alignment gaps as missing data, we propose a simple nonparametric test of the null hypothesis that the locations of alignment gaps are independent of the nucleotide substitution or amino acid replacement process. When we apply the test to 1,390 protein alignments that are informed by protein tertiary structure and use a 5% significance level, the null hypothesis of independence between amino acid replacement and gap location is rejected for ∼65% of datasets. Via simulations that include substitution and insertion-deletion, we show that the test performs well with true alignments. When we simulate according to the null hypothesis and then apply the test to optimal alignments that are inferred by each of four widely used software packages, the null hypothesis is rejected too frequently. Via further simulations and analyses, we show that the overly frequent rejections of the null hypothesis are not solely due to weaknesses of widely used software for finding optimal alignments. Instead, our evidence suggests that optimal alignments are unrepresentative of true alignments and that biased evolutionary inferences may result from relying upon individual optimal alignments.

Precision targeting tumor cells using cancer-specific InDel mutations with CRISPR-Cas9.

An ideal cancer therapeutic strategy involves the selective killing of cancer cells without affecting the surrounding normal cells. However, researchers have failed to develop such methods for achieving selective cancer cell death because of shared features between cancerous and normal cells. In this study, we have developed a therapeutic strategy called the cancer-specific insertions-deletions (InDels) attacker (CINDELA) to selectively induce cancer cell death using the CRISPR-Cas system. CINDELA utilizes a previously unexplored idea of introducing CRISPR-mediated DNA double-strand breaks (DSBs) in a cancer-specific fashion to facilitate specific cell death. In particular, CINDELA targets multiple InDels with CRISPR-Cas9 to produce many DNA DSBs that result in cancer-specific cell death. As a proof of concept, we demonstrate here that CINDELA selectively kills human cancer cell lines, xenograft human tumors in mice, patient-derived glioblastoma, and lung patient-driven xenograft tumors without affecting healthy human cells or altering mouse growth.

A 163-bp insertion in the Capana10g000198 encoding a MYB transcription factor causes male sterility in pepper (Capsicum annuum L.).

Male sterility provides an efficient approach for commercial exploitation of heterosis. Despite more than 20 genic male sterile (GMS) mutants documented in pepper (Capsicum annuum L.), only two causal genes have been successfully identified. Here, a novel spontaneous recessive GMS mutant, designated msc-3, is identified and characterized at both phenotypic and histological levels. Pollen abortion of msc-3 mutant may be due to the delayed tapetum degradation, leading to the non-degeneration of tetrads callosic wall. Then, a modified MutMap method and molecular marker linkage analysis were employed to fine mapping the msc-3 locus, which was delimited to the ~139.91-kb region harboring 10 annotated genes. Gene expression and structure variation analyses indicate the Capana10g000198, encoding a R2R3-MYB transcription factor, is the best candidate gene for the msc-3 locus. Expression profiling analysis shows the Capana10g000198 is an anther-specific gene, and a 163-bp insertion in the Capana10g000198 is highly correlated with the male sterile (MS) phenotype. Additionally, downregulation of Capana10g000198 in male fertile plants through virus-induced gene silencing resulted in male sterility. Finally, possible regulatory relationships of the msc-3 gene with the other two reported pepper GMS genes, msc-1 and msc-2, have been studied, and comparative transcriptome analysis reveals the expression of 16 GMS homologs are significantly downregulated in the MS anthers. Overall, our results reveal that Capana10g000198 is the causal gene underlying the msc-3 locus, providing important theoretical clues and basis for further in-depth study on the regulatory mechanisms of pollen development in pepper.

The role of leucine and isoleucine in tuning the hydropathy of class A GPCRs.

Leucine and Isoleucine are two amino acids that differ only by the positioning of one methyl group. This small difference can have important consequences in α-helices, as the ß-branching of Ile results in helix destabilization. We set out to investigate whether there are general trends for the occurrences of Leu and Ile residues in the structures and sequences of class A GPCRs (G protein-coupled receptors). GPCRs are integral membrane proteins in which α-helices span the plasma membrane seven times and which play a crucial role in signal transmission. We found that Leu side chains are generally more exposed at the protein surface than Ile side chains. We explored whether this difference might be attributed to different functions of the two amino acids and tested if Leu tunes the hydrophobicity of the transmembrane domain based on the Wimley-White whole-residue hydrophobicity scales. Leu content decreases the variation in hydropathy between receptors and correlates with the non-Leu receptor hydropathy. Both measures indicate that hydropathy is tuned by Leu. To test this idea further, we generated protein sequences with random amino acid compositions using a simple numerical model, in which hydropathy was tuned by adjusting the number of Leu residues. The model was able to replicate the observations made with class A GPCR sequences. We speculate that the hydropathy of transmembrane domains of class A GPCRs is tuned by Leu (and to some lesser degree by Lys and Val) to facilitate correct insertion into membranes and/or to stably anchor the receptors within membranes.

Characterization of active transposable elements and their new insertions in tuber propagated greater yam (Dioscorea alata).

BACKGROUND: Greater yam is a key staple crop grown in tropical and subtropical regions, while its asexual propagation mode had led to non-flowering mutations. How transposable elements contribute to its genetic variations is rarely analyzed. We used transcriptome and whole genome sequencing data to identify active transposable elements (TEs) and genetic variation caused by these active TEs. Our aim was to shed light on which TEs would lead to its intraspecies variation. RESULTS: Annotation of de novo assembly transcripts indicated that 0.8 - 0.9% of transcripts were TE related, with LTR retrotransposons (LTR-RTs) accounted for 65% TE transcripts. A large portion of these transcripts were non-autonomous TEs, which had incomplete functional domains. The majority of mapped transcripts were distributed in genic deficient regions, with 9% of TEs overlapping with genic regions. Moreover, over 90% TE transcripts exhibited low expression levels and insufficient reads coverage to support full-length structure assembly. Subfamily analysis of Copia and Gypsy, the two LTR-RTs revealed that a small number of subfamilies contained a significantly larger number of members, which play a key role in generating TE transcript. Based on resequencing data, 15,002 L-RT insertion loci were detected for active LTR-RT members. The insertion loci of LTR-RTs were highly divergent among greater yam accessions. CONCLUSIONS: This study showed the ongoing transcription and transpositions of TEs in greater yam, despite low transcription levels and incomplete proteins insufficient for autonomous transposition. While our research did not directly link these TEs to specific yam traits such as tuber yield and propagation mode, it lays a crucial foundation for further research on how these TE insertion polymorphisms (TIPs) might be related to variations in greater yam traits and its tuber propagation mode. Future research may explore the potential roles of TEs in trait variations, such as tuber yield and stress resistance, in greater yam.

Developing an SNP dataset for efficiently evaluating soybean germplasm resources using the genome sequencing data of 3,661 soybean accessions.

BACKGROUND: Single nucleotide polymorphism (SNP) markers play significant roles in accelerating breeding and basic crop research. Several soybean SNP panels have been developed. However, there is still a lack of SNP panels for differentiating between wild and cultivated populations, as well as for detecting polymorphisms within both wild and cultivated populations. RESULTS: This study utilized publicly available resequencing data from over 3,000 soybean accessions to identify differentiating and highly conserved SNP and insertion/deletion (InDel) markers between wild and cultivated soybean populations. Additionally, a naturally occurring mutant gene library was constructed by analyzing large-effect SNPs and InDels in the population. CONCLUSION: The markers obtained in this study are associated with numerous genes governing agronomic traits, thus facilitating the evaluation of soybean germplasms and the efficient differentiation between wild and cultivated soybeans. The natural mutant gene library permits the quick identification of individuals with natural mutations in functional genes, providing convenience for accelerating soybean breeding using reverse genetics.