Chasing Sequencing Perfection: Marching Toward Higher Accuracy and Lower Costs.

Next-generation sequencing (NGS), represented by Illumina platforms, has been an essential cornerstone of basic and applied research. However, the sequencing error rate of 1 per 1000 bp (10-3) represents a serious hurdle for research areas focusing on rare mutations, such as somatic mosaicism or microbe heterogeneity. By examining the high-fidelity sequencing methods developed in the past decade, we summarized three major factors underlying errors and the corresponding 12 strategies mitigating these errors. We then proposed a novel framework to classify 11 preexisting representative methods according to the corresponding combinatory strategies and identified three trends that emerged during methodological developments. We further extended this analysis to eight long-read sequencing methods, emphasizing error reduction strategies. Finally, we suggest two promising future directions that could achieve comparable or even higher accuracy with lower costs in both NGS and long-read sequencing.

Low-input PacBio sequencing generates high-quality individual fly genomes and characterizes mutational processes.

Long-read sequencing, exemplified by PacBio, revolutionizes genomics, overcoming challenges like repetitive sequences. However, the high DNA requirement ( > 1 µg) is prohibitive for small organisms. We develop a low-input (100 ng), low-cost, and amplification-free library-generation method for PacBio sequencing (LILAP) using Tn5-based tagmentation and DNA circularization within one tube. We test LILAP with two Drosophila melanogaster individuals, and generate near-complete genomes, surpassing preexisting single-fly genomes. By analyzing variations in these two genomes, we characterize mutational processes: complex transpositions (transposon insertions together with extra duplications and/or deletions) prefer regions characterized by non-B DNA structures, and gene conversion of transposons occurs on both DNA and RNA levels. Concurrently, we generate two complete assemblies for the endosymbiotic bacterium Wolbachia in these flies and similarly detect transposon conversion. Thus, LILAP promises a broad PacBio sequencing adoption for not only mutational studies of flies and their symbionts but also explorations of other small organisms or precious samples.

Genome editing in plants using the TnpB transposase system.

The widely used clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated nuclease (Cas) system is thought to have evolved from IS200/IS605 transposons. TnpB proteins, encoded by one type of IS200/IS605 transposon, are considered to be the evolutionary ancestors of Cas12 nucleases, which have been engineered to function as RNA-guided DNA endonucleases for genome editing in bacteria and human cells. TnpB nucleases, which are smaller than Cas nucleases, have been engineered for use in genome editing in animal systems, but the feasibility of this approach in plants remained unknown. Here, we obtained stably transformed genome-edited mutants in rice (Oryza sativa) by adapting three recently identified TnpB genome editing vectors, encoding distinct TnpB nucleases (ISAam1, ISDra2, and ISYmu1), for use in plants, demonstrating that the hypercompact TnpB proteins can effectively edit plant genomes. ISDra2 and ISYmu1 precisely edited their target sequences, with no off-target mutations detected, showing that TnpB transposon nucleases are suitable for development into a new genome editing tool for plants. Future modifications improving the genome-editing efficiency of the TnpB system will facilitate plant functional studies and breeding programs. Supplementary Information: The online version contains supplementary material available at 10.1007/s42994-024-00172-6.

Heterologous survey of 130 DNA transposons in human cells highlights their functional divergence and expands the genome engineering toolbox.

Experimental studies on DNA transposable elements (TEs) have been limited in scale, leading to a lack of understanding of the factors influencing transposition activity, evolutionary dynamics, and application potential as genome engineering tools. We predicted 130 active DNA TEs from 102 metazoan genomes and evaluated their activity in human cells. We identified 40 active (integration-competent) TEs, surpassing the cumulative number (20) of TEs found previously. With this unified comparative data, we found that the Tc1/mariner superfamily exhibits elevated activity, potentially explaining their pervasive horizontal transfers. Further functional characterization of TEs revealed additional divergence in features such as insertion bias. Remarkably, in CAR-T therapy for hematological and solid tumors, Mariner2_AG (MAG), the most active DNA TE identified, largely outperformed two widely used vectors, the lentiviral vector and the TE-based vector SB100X. Overall, this study highlights the varied transposition features and evolutionary dynamics of DNA TEs and increases the TE toolbox diversity.

Hagfish genome elucidates vertebrate whole-genome duplication events and their evolutionary consequences.

Polyploidy or whole-genome duplication (WGD) is a major event that drastically reshapes genome architecture and is often assumed to be causally associated with organismal innovations and radiations. The 2R hypothesis suggests that two WGD events (1R and 2R) occurred during early vertebrate evolution. However, the timing of the 2R event relative to the divergence of gnathostomes (jawed vertebrates) and cyclostomes (jawless hagfishes and lampreys) is unresolved and whether these WGD events underlie vertebrate phenotypic diversification remains elusive. Here we present the genome of the inshore hagfish, Eptatretus burgeri. Through comparative analysis with lamprey and gnathostome genomes, we reconstruct the early events in cyclostome genome evolution, leveraging insights into the ancestral vertebrate genome. Genome-wide synteny and phylogenetic analyses support a scenario in which 1R occurred in the vertebrate stem-lineage during the early Cambrian, and 2R occurred in the gnathostome stem-lineage, maximally in the late Cambrian-earliest Ordovician, after its divergence from cyclostomes. We find that the genome of stem-cyclostomes experienced an additional independent genome triplication. Functional genomic and morphospace analyses demonstrate that WGD events generally contribute to developmental evolution with similar changes in the regulatory genome of both vertebrate groups. However, appreciable morphological diversification occurred only in the gnathostome but not in the cyclostome lineage, calling into question the general expectation that WGDs lead to leaps of bodyplan complexity.

Forces driving transposable element load variation during Arabidopsis range expansion.

Genetic load refers to the accumulated and potentially life-threatening deleterious mutations in populations. Understanding the mechanisms underlying genetic load variation of transposable element (TE) insertion, a major large-effect mutation, during range expansion is an intriguing question in biology. Here, we used 1,115 global natural accessions of Arabidopsis (Arabidopsis thaliana) to study the driving forces of TE load variation during its range expansion. TE load increased with range expansion, especially in the recently established Yangtze River basin population. Effective population size, which explains 62.0% of the variance in TE load, high transposition rate, and selective sweeps contributed to TE accumulation in the expanded populations. We genetically mapped and identified multiple candidate causal genes and TEs, and revealed the genetic architecture of TE load variation. Overall, this study reveals the variation in TE genetic load during Arabidopsis expansion and highlights the causes of TE load variation from the perspectives of both population genetics and quantitative genetics.

Evolutionarily new genes in humans with disease phenotypes reveal functional enrichment patterns shaped by adaptive innovation and sexual selection.

New genes (or young genes) are structural novelties pivotal in mammalian evolution. Their phenotypic impacts on humans, however, remain elusive due to the technical and ethical complexities in functional studies. Through combining gene age dating with Mendelian disease phenotyping, our research reveals a steady integration of new genes with biomedical phenotypes into the human genome over macroevolutionary timescales (~0.07% per million years). Despite this stable pace, we observe distinct patterns in phenotypic enrichment, pleiotropy, and selective pressures shaped by different gene ages. Notably, young genes show significant enrichment in the male reproductive system, indicating strong sexual selection. Young genes also exhibit functions in tissues and systems potentially linked to human phenotypic innovations, such as increased brain size, musculoskeletal phenotypes, and color vision. Our findings further reveal increasing levels of pleiotropy over evolutionary time, which accompanies stronger selective constraints. We propose a "pleiotropy-barrier" model that delineates different potentials for phenotypic innovation between young and older genes subject to natural selection. Our study demonstrates that evolutionary new genes are critical in influencing human reproductive evolution and adaptive phenotypic innovations driven by sexual and natural selection, with low pleiotropy as a selective advantage.

The Morphological Transformation of the Thorax during the Eclosion of Drosophila melanogaster (Diptera: Drosophilidae).

The model organism Drosophila melanogaster, as a species of Holometabola, undergoes a series of transformations during metamorphosis. To deeply understand its development, it is crucial to study its anatomy during the key developmental stages. We describe the anatomical systems of the thorax, including the endoskeleton, musculature, nervous ganglion, and digestive system, from the late pupal stage to the adult stage, based on micro-CT and 3D visualizations. The development of the endoskeleton causes original and insertional changes in muscles. Several muscles change their shape during development in a non-uniform manner with respect to both absolute and relative size; some become longer and broader, while others shorten and become narrower. Muscular shape may vary during development. The number of muscular bundles also increases or decreases. Growing muscles are probably anchored by the tissues in the stroma. Some muscles and tendons are absent in the adult stage, possibly due to the hardened sclerites. Nearly all flight muscles are present by the third day of the pupal stage, which may be due to the presence of more myofibers with enough mitochondria to support flight power. There are sexual differences in the same developmental period. In contrast to the endodermal digestive system, the functions of most thoracic muscles change in the development from the larva to the adult in order to support more complex locomotion under the control of a more structured ventral nerve cord based on the serial homology proposed herein.

Adaptive evolution of the enigmatic Takakia now facing climate change in Tibet.

The most extreme environments are the most vulnerable to transformation under a rapidly changing climate. These ecosystems harbor some of the most specialized species, which will likely suffer the highest extinction rates. We document the steepest temperature increase (2010-2021) on record at altitudes of above 4,000 m, triggering a decline of the relictual and highly adapted moss Takakia lepidozioides. Its de-novo-sequenced genome with 27,467 protein-coding genes includes distinct adaptations to abiotic stresses and comprises the largest number of fast-evolving genes under positive selection. The uplift of the study site in the last 65 million years has resulted in life-threatening UV-B radiation and drastically reduced temperatures, and we detected several of the molecular adaptations of Takakia to these environmental changes. Surprisingly, specific morphological features likely occurred earlier than 165 mya in much warmer environments. Following nearly 400 million years of evolution and resilience, this species is now facing extinction.

Evolutionary mining and functional characterization of TnpB nucleases identify efficient miniature genome editors.

As the evolutionary ancestor of Cas12 nuclease, the transposon (IS200/IS605)-encoded TnpB proteins act as compact RNA-guided DNA endonucleases. To explore their evolutionary diversity and potential as genome editors, we screened TnpBs from 64 annotated IS605 members and identified 25 active in Escherichia coli, of which three are active in human cells. Further characterization of these 25 TnpBs enables prediction of the transposon-associated motif (TAM) and the right-end element RNA (reRNA) directly from genomic sequences. We established a framework for annotating TnpB systems in prokaryotic genomes and applied it to identify 14 additional candidates. Among these, ISAam1 (369 amino acids (aa)) and ISYmu1 (382 aa) TnpBs demonstrated robust editing activity across dozens of genomic loci in human cells. Both RNA-guided genome editors demonstrated similar editing efficiency as SaCas9 (1,053 aa) while being substantially smaller. The enormous diversity of TnpBs holds potential for the discovery of additional valuable genome editors.

The power of "controllers": Transposon-mediated duplicated genes evolve towards neofunctionalization.

Since the discovery of the first transposon by Dr. Barbara McClintock, the prevalence and diversity of transposable elements (TEs) have been gradually recognized. As fundamental genetic components, TEs drive organismal evolution not only by contributing functional sequences (e.g., regulatory elements or "controllers" as phrased by Dr. McClintock) but also by shuffling genomic sequences. In the latter respect, TE-mediated gene duplications have contributed to the origination of new genes and attracted extensive interest. In response to the development of this field, we herein attempt to provide an overview of TE-mediated duplication by focusing on common rules emerging across duplications generated by different TE types. Specifically, despite the huge divergence of transposition machinery across TEs, we identify three common features of various TE-mediated duplication mechanisms, including end bypass, template switching, and recurrent transposition. These three features lead to one common functional outcome, namely, TE-mediated duplicates tend to be subjected to exon shuffling and neofunctionalization. Therefore, the intrinsic properties of the mutational mechanism constrain the evolutionary trajectories of these duplicates. We finally discuss the future of this field including an in-depth characterization of both the duplication mechanisms and functions of TE-mediated duplicates.

Divergent contributions of coding and noncoding sequences to initial high-altitude adaptation in passerine birds endemic to the Qinghai-Tibet Plateau.

Early events in the evolution of an ancestral lineage can shape the adaptive patterns of descendant species, but the evolutionary mechanisms driving initial adaptation from an ancestor remain largely unexplored. High-altitude adaptations have been extensively explored from the viewpoint of protein-coding genes; however, the contribution of noncoding regions remains relatively neglected. Here, we integrate genomic and transcriptomic data to investigate adaptive evolution in the ancestor of three high-altitude snowfinch species endemic to the Qinghai-Tibet Plateau. Our genome-wide scan for adaptation in the snowfinch ancestor identifies strong adaptation signals in functions of development and metabolism for the coding genes, but in functions of the nervous system development for noncoding regions. This pattern is exclusive to the snowfinch ancestor compared to a control ancestral lineage subject to weak selection. Changes in noncoding regions in the snowfinch ancestor, especially those nearest to coding genes, may be disproportionately associated with the differential expression of genes in the brain tissue compared to other tissues. Extensive gene expression in the brain tissue can be further altered via genetic regulatory networks of transcription factors harbouring potential accelerated regulatory regions (e.g., the development-related transcription factor YEATS4). Altogether, our study provides new evidence concerning how coding and noncoding sequences work through decoupled pathways in initial adaptation to the selective pressure of high-altitude environments. The analysis highlights the idea that noncoding sequences may be promising elements in facilitating the rapid evolution and adaptation to high altitudes.

Resurrection of endogenous retroviruses during aging reinforces senescence.

Whether and how certain transposable elements with viral origins, such as endogenous retroviruses (ERVs) dormant in our genomes, can become awakened and contribute to the aging process is largely unknown. In human senescent cells, we found that HERVK (HML-2), the most recently integrated human ERVs, are unlocked to transcribe viral genes and produce retrovirus-like particles (RVLPs). These HERVK RVLPs constitute a transmissible message to elicit senescence phenotypes in young cells, which can be blocked by neutralizing antibodies. The activation of ERVs was also observed in organs of aged primates and mice as well as in human tissues and serum from the elderly. Their repression alleviates cellular senescence and tissue degeneration and, to some extent, organismal aging. These findings indicate that the resurrection of ERVs is a hallmark and driving force of cellular senescence and tissue aging.

Pan-cancer surveys indicate cell cycle-related roles of primate-specific genes in tumors and embryonic cerebrum.

BACKGROUND: Despite having been extensively studied, it remains largely unclear why humans bear a particularly high risk of cancer. The antagonistic pleiotropy hypothesis predicts that primate-specific genes (PSGs) tend to promote tumorigenesis, while the molecular atavism hypothesis predicts that PSGs involved in tumors may represent recently derived duplicates of unicellular genes. However, these predictions have not been tested. RESULTS: By taking advantage of pan-cancer genomic data, we find the upregulation of PSGs across 13 cancer types, which is facilitated by copy-number gain and promoter hypomethylation. Meta-analyses indicate that upregulated PSGs (uPSGs) tend to promote tumorigenesis and to play cell cycle-related roles. The cell cycle-related uPSGs predominantly represent derived duplicates of unicellular genes. We prioritize 15 uPSGs and perform an in-depth analysis of one unicellular gene-derived duplicate involved in the cell cycle, DDX11. Genome-wide screening data and knockdown experiments demonstrate that DDX11 is broadly essential across cancer cell lines. Importantly, non-neutral amino acid substitution patterns and increased expression indicate that DDX11 has been under positive selection. Finally, we find that cell cycle-related uPSGs are also preferentially upregulated in the highly proliferative embryonic cerebrum. CONCLUSIONS: Consistent with the predictions of the atavism and antagonistic pleiotropy hypotheses, primate-specific genes, especially those PSGs derived from cell cycle-related genes that emerged in unicellular ancestors, contribute to the early proliferation of the human cerebrum at the cost of hitchhiking by similarly highly proliferative cancer cells.

Ruminant-specific genes identified using high-quality genome data and their roles in rumen evolution.

Ruminants comprise a highly successful group of mammals with striking morphological innovations, including the presence of a rumen. Many studies have shown that species-specific or lineage-specific genes (referred to as new genes) play important roles in phenotypic evolution. In this study, we identified 1064 ruminant-specific genes based on the newly assembled high-quality genomes of representative members of two ruminant families and other publically available high-quality genomes. Ruminant-specific genes shared similar evolutionary and expression patterns with new genes found in other mammals, such as primates and rodents. Most new genes were derived from gene duplication and tended to be expressed in the testes or immune-related tissues, but were depleted in the adult brain. We also found that most genes expressed in the rumen were genes predating sheep-sperm whale split (referred to as old genes), but some new genes were also involved in the evolution of the rumen, and contributed more during rumen development than in the adult rumen. Notably, expression levels of members of the ruminant-specific PRD-SPRRII gene family, which are subject to positive selection, varied throughout rumen development and may thus play important roles in the development of the keratin-rich surface of the rumen. Overall, this study generated two novel ruminant genomes and also provided novel insights into the evolution of new mammalian organs.

Proteomic profiling of serum extracellular vesicles identifies diagnostic markers for echinococcosis.

Echinococcosis is a parasitic disease caused by the metacestodes of Echinococcus spp. The disease has a long latent period and is largely underdiagnosed, partially because of the lack of effective early diagnostic approaches. Using liquid chromatography-mass spectrometry, we profiled the serum-derived extracellular vesicles (EVs) of E. multilocularis-infected mice and identified three parasite-origin proteins, thioredoxin peroxidase 1 (TPx-1), transitional endoplasmic reticulum ATPase (TER ATPase), and 14-3-3, being continuously released by the parasites into the sera during the infection via EVs. Using ELISA, both TPx-1 and TER ATPase were shown to have a good performance in diagnosis of experimental murine echinococcosis as early as 10 days post infection and of human echinococcosis compared with that of control. Moreover, TER ATPase and TPx-1 were further demonstrated to be suitable for evaluation of the prognosis of patients with treatment. The present study discovers the potential of TER ATPase and TPx-1 as promising diagnostic candidates for echinococcosis.

Structural variations between small alarmone hydrolase dimers support different modes of regulation of the stringent response.

The bacterial stringent response involves wide-ranging metabolic reprogramming aimed at increasing long-term survivability during stress conditions. One of the hallmarks of the stringent response is the production of a set of modified nucleotides, known as alarmones, which affect a multitude of cellular pathways in diverse ways. Production and degradation of these molecules depend on the activity of enzymes from the RelA/SpoT homologous family, which come in both bifunctional (containing domains to both synthesize and hydrolyze alarmones) and monofunctional (consisting of only synthetase or hydrolase domain) variants, of which the structure, activity, and regulation of the bifunctional RelA/SpoT homologs have been studied most intensely. Despite playing an important role in guanosine nucleotide homeostasis in particular, mechanisms of regulation of the small alarmone hydrolases (SAHs) are still rather unclear. Here, we present crystal structures of SAH enzymes from Corynebacterium glutamicum (RelHCg) and Leptospira levettii (RelHLl) and show that while being highly similar, structural differences in substrate access and dimer conformations might be important for regulating their activity. We propose that a varied dimer form is a general property of the SAH family, based on current structural information as well as prediction models for this class of enzymes. Finally, subtle structural variations between monofunctional and bifunctional enzymes point to how these different classes of enzymes are regulated.

[Effects of CO2 concentration and soil water content on short-term water-use efficiency at whole-plant level]. / CO2浓度和土壤含水量对植物个体尺度短期水分利用效率的影响.

Uncovering the variations of short-term water-use efficiency (WUEp) at whole-plant level in response to CO2 concentration (Ca) and soil water content (SWC) can improve the understanding of plant survival strategies under climate change. In this study, Platycladus orientalis saplings were cultured in simulated climate chambers.There were totally 15 treatments, including Ca of 400 (C400), 600 (C600) and 800 (C800) µmol·mol-1 and SWC of 35%-45% field water holding capacity (FC), 50%-60%FC, 60%-70%FC, 70%-80%FC and 95%-100%FC. The WUEp was measured by mini-lysimeters, weighting method, and static assimilation chamber. The results showed that both daytime (0.12-1.87 mol·h-1) and nighttime transpiration rates (0.01-0.16 mol·h-1) at whole-plant level reached the maximum at C400×70%-80%FC, while the whole-plant daytime net photosynthetic rate (2.12-22.10 mmol·h-1) reached the maximum at C800×70%-80%FC. In contrast, nighttime respiration rate (0.84-4.41 mmol·h-1) increased with increasing SWC, but decreased with increasing of Ca, reaching the maximum at C400×95%-100%FC. For WUEp (5.37-24.35 mmol·mol-1), it reached the maximum at C800×50%-60%FC, indicating that plants could use less water and fixed more carbon by adjusting adaptation strategies under high Ca and drought conditions. In addition, leaf instantaneous water-use efficiency was a good predictor of WUEP when the canopy structure was similar.

Dosage sensitivity and exon shuffling shape the landscape of polymorphic duplicates in Drosophila and humans.

Despite polymorphic duplicate genes' importance for the early stages of duplicate gene evolution, they are less studied than old gene duplicates. Two essential questions thus remain poorly addressed: how does dosage sensitivity, imposed by stoichiometry in protein complexes or by X chromosome dosage compensation, affect the emergence of complete duplicate genes? Do introns facilitate intergenic and intragenic chimaerism as predicted by the theory of exon shuffling? Here, we analysed new data for Drosophila and public data for humans, to characterize polymorphic duplicate genes with respect to dosage, exon-intron structures and allele frequencies. We found that complete duplicate genes are under dosage constraint induced by protein stoichiometry but potentially tolerated by X chromosome dosage compensation. We also found that in the intron-rich human genome, gene fusions and intragenic duplications extensively use intronic breakpoints generating in-frame proteins, in accordance with the theory of exon shuffling. Finally, we found that only a small proportion of complete or partial duplicates are at high frequencies, indicating the deleterious nature of dosage or gene structural changes. Altogether, we demonstrate how mechanistic factors including dosage sensitivity and exon-intron structure shape the short-term functional consequences of gene duplication.