Navigating the pitfalls of mapping DNA and RNA modifications.

Chemical modifications to nucleic acids occur across the kingdoms of life and carry important regulatory information. Reliable high-resolution mapping of these modifications is the foundation of functional and mechanistic studies, and recent methodological advances based on next-generation sequencing and long-read sequencing platforms are critical to achieving this aim. However, mapping technologies may have limitations that sometimes lead to inconsistent results. Some of these limitations are technical in nature and specific to certain types of technology. Here, however, we focus on common (yet not always widely recognized) pitfalls that are shared among frequently used mapping technologies and discuss strategies to help technology developers and users mitigate their effects. Although the emphasis is primarily on DNA modifications, RNA modifications are also discussed.

mEnrich-seq: methylation-guided enrichment sequencing of bacterial taxa of interest from microbiome.

Metagenomics has enabled the comprehensive study of microbiomes. However, many applications would benefit from a method that sequences specific bacterial taxa of interest, but not most background taxa. We developed mEnrich-seq (in which 'm' stands for methylation and seq for sequencing) for enriching taxa of interest from metagenomic DNA before sequencing. The core idea is to exploit the self versus nonself differentiation by natural bacterial DNA methylation and rationally choose methylation-sensitive restriction enzymes, individually or in combination, to deplete host and background taxa while enriching targeted taxa. This idea is integrated with library preparation procedures and applied in several applications to enrich (up to 117-fold) pathogenic or beneficial bacteria from human urine and fecal samples, including species that are hard to culture or of low abundance. We assessed 4,601 bacterial strains with mapped methylomes so far and showed broad applicability of mEnrich-seq. mEnrich-seq provides microbiome researchers with a versatile and cost-effective approach for selective sequencing of diverse taxa of interest.

The Landscape of A-to-I RNA Editome Is Shaped by Both Positive and Purifying Selection.

The hydrolytic deamination of adenosine to inosine (A-to-I editing) in precursor mRNA induces variable gene products at the post-transcription level. How and to what extent A-to-I RNA editing diversifies transcriptome is not fully characterized in the evolution, and very little is known about the selective constraints that drive the evolution of RNA editing events. Here we present a study on A-to-I RNA editing, by generating a global profile of A-to-I editing for a phylogeny of seven Drosophila species, a model system spanning an evolutionary timeframe of approximately 45 million years. Of totally 9281 editing events identified, 5150 (55.5%) are located in the coding sequences (CDS) of 2734 genes. Phylogenetic analysis places these genes into 1,526 homologous families, about 5% of total gene families in the fly lineages. Based on conservation of the editing sites, the editing events in CDS are categorized into three distinct types, representing events on singleton genes (type I), and events not conserved (type II) or conserved (type III) within multi-gene families. While both type I and II events are subject to purifying selection, notably type III events are positively selected, and highly enriched in the components and functions of the nervous system. The tissue profiles are documented for three editing types, and their critical roles are further implicated by their shifting patterns during holometabolous development and in post-mating response. In conclusion, three A-to-I RNA editing types are found to have distinct evolutionary dynamics. It appears that nervous system functions are mainly tested to determine if an A-to-I editing is beneficial for an organism. The coding plasticity enabled by A-to-I editing creates a new class of binary variations, which is a superior alternative to maintain heterozygosity of expressed genes in a diploid mating system.

Decipher the ancestry of the plant-specific LBD gene family.

BACKGROUND: Lateral Organ Boundaries Domain (LBD) genes arise from charophyte algae and evolve essential functions in land plants in regulating organ development and secondary metabolism. Although diverse plant species have been investigated to construct the phylogeny of LBD gene family, a detailed and reliable ancestry that characterizes their evolutionary patterns has not been revealed. RESULTS: We develop an improved bioinformatic method that allows robust detection of 431 LBD genes in 11 high-quality land plant genomes. Phylogenetic analysis classifies the LBD genes into six subfamilies which support the existence of 7 ancient gene lineages. Phylogenetic relationship and gene collinearity are combined to retrace 11 ancestor genes for seed plants and 18 ancestor genes for angiosperms, which improves the resolution of LBD gene ancestry. The ancient gene lineages are strictly preserved in current plant genomes, including the previously controversial class IB gene in Selaginella moellendorphii, suggesting extreme reluctance of LBD genes to be lost during evolution. Meanwhile, whole-genome and dispersed gene duplications substantially expand LBD gene family in angiosperms, and elaborate functions of LBD genes through frequent expression pattern change and protein sequence variation. CONCLUSIONS: Through phylogenetic and gene collinearity analyses, we retrace the landscape of LBD gene ancestry which lays foundation for elucidating evolutionary diversification of LBD genes in land plants.

Conservation and functional influence of alternative splicing in wood formation of Populus and Eucalyptus.

BACKGROUND: Wood formation in tree species is regulated by multiple factors at various layers. Alternative splicing (AS) occurs within a large number of genes in wood formation. However, the functional implications and conservation of the AS occurrence are not well understood. RESULTS: In this study, we profiled AS events in wood-forming tissues of Populus and Eucalyptus, and analyzed their functional implications as well as inter-species conservation. 28.3% and 20.7% of highly expressed transcripts in the developing xylem of Populus and Eucalyptus respectively were affected by AS events. Around 42% of the AS events resulted in changes to the original reading frame. 25.0% (in Populus) and 26.8% (in Eucalyptus) of the AS events may cause protein domain modification. In the process of wood formation, about 28% of AS-occurring genes were putative orthologs and 71 conserved AS events were identified in the two species. CONCLUSION: Through analysis of AS events in developing xylem of two tree species, this study reveals an array of new information regarding AS occurrence and function in tree development.

Hygromycin-resistance vectors for gene expression in Pichia pastoris.

Pichia pastoris is a common host organism for heterologous protein expression and metabolic engineering. Zeocin-, G418-, nourseothricin- and blasticidin-resistance genes are the only dominant selectable markers currently available for selecting P. pastoris transformants. We describe here new P. pastoris expression vectors that confer a hygromycin resistance base on the Klebsiella pneumoniae hph gene. To demonstrate the application of the vectors for intracellular and secreted protein expression, green fluorescent protein (GFP) and human serum albumin (HSA) were cloned into the vectors and transformed into P. pastoris cells. The resulting strains expressed GFP and HSA constitutively or inducibly. The hygromycin resistance marker was also suitable for post-transformational vector amplication (PTVA) for obtaining strains with high plasmid copy numbers. A strain with multiple copies of the HSA expression cassette after PTVA had increased HSA expression compared with a strain with a single copy of the plasmid. To demonstrate compatibility of the new vectors with other vectors bearing antibiotic-resistance genes, P. pastoris was transformed with the Saccharomyces cerevisiae genes GSH1, GSH2 or SAM2 on plasmids containing genes for resistance to Zeocin, G418 or hygromycin. The resulting strain produced glutathione and S-adenosyl-L-methionine at levels approximately twice those of the parent strain. The new hygromycin-resistance vectors allow greater flexibility and potential applications in recombinant protein production and other research using P. pastoris.

Integrated multi-omic analyses uncover the effects of aging on cell-type regulation in glucose-responsive tissues.

Aging significantly influences cellular activity and metabolism in glucose-responsive tissues, yet a comprehensive evaluation of the impacts of aging and associated cell-type responses has been lacking. This study integrates transcriptomic, methylomic, single-cell RNA sequencing, and metabolomic data to investigate aging-related regulations in adipose and muscle tissues. Through coexpression network analysis of the adipose tissue, we identified aging-associated network modules specific to certain cell types, including adipocytes and immune cells. Aging upregulates the metabolic functions of lysosomes and downregulates the branched-chain amino acids (BCAAs) degradation pathway. Additionally, aging-associated changes in cell proportions, methylation profiles, and single-cell expressions were observed in the adipose. In the muscle tissue, aging was found to repress the metabolic processes of glycolysis and oxidative phosphorylation, along with reduced gene activity of fast-twitch type II muscle fibers. Metabolomic profiling linked aging-related alterations in plasma metabolites to gene expression in glucose-responsive tissues, particularly in tRNA modifications, BCAA metabolism, and sex hormone signaling. Together, our multi-omic analyses provide a comprehensive understanding of the impacts of aging on glucose-responsive tissues and identify potential plasma biomarkers for these effects.

Identification of molecular processes needed for vascular formation through transcriptome analysis of different vascular systems.

BACKGROUND: Vascular system formation has been studied through molecular and genetic approaches in Arabidopsis, a herbaceous dicot that is used as a model system. Different vascular systems have developed in other plants such as crops and trees. Uncovering shared mechanisms underlying vascular development by transcriptome analysis of different vascular systems may help to transfer knowledge acquired from Arabidopsis to other economically important species. RESULTS: Conserved vascular genes and biological processes fundamental to vascular development were explored across various plants. Through comparative transcriptome analysis, 226 genes from Arabidopsis, 217 genes from poplar and 281 genes from rice were identified as constituting 107 conserved vascular gene groups. These gene groups are expressed mainly in vascular tissues and form a complex coexpression network with multiple functional connections. To date, only half of the groups have been experimentally investigated. The conserved vascular gene groups were classified into 9 essential processes for vascular development. 18 groups (17%) lack of annotations were classified as having unknown functions. CONCLUSION: The study provides a map of fundamental biological processes conserved across different vascular systems. It identifies gaps in the experimental investigation of pathways active in vascular formation, which if explored, could lead to a more complete understanding of vascular development.

Global transcriptome profiles of Camellia sinensis during cold acclimation.

BACKGROUND: Tea is the most popular non-alcoholic health beverage in the world. The tea plant (Camellia sinensis (L.) O. Kuntze) needs to undergo a cold acclimation process to enhance its freezing tolerance in winter. Changes that occur at the molecular level in response to low temperatures are poorly understood in tea plants. To elucidate the molecular mechanisms of cold acclimation, we employed RNA-Seq and digital gene expression (DGE) technologies to the study of genome-wide expression profiles during cold acclimation in tea plants. RESULTS: Using the Illumina sequencing platform, we obtained approximately 57.35 million RNA-Seq reads. These reads were assembled into 216,831 transcripts, with an average length of 356 bp and an N50 of 529 bp. In total, 1,770 differentially expressed transcripts were identified, of which 1,168 were up-regulated and 602 down-regulated. These include a group of cold sensor or signal transduction genes, cold-responsive transcription factor genes, plasma membrane stabilization related genes, osmosensing-responsive genes, and detoxification enzyme genes. DGE and quantitative RT-PCR analysis further confirmed the results from RNA-Seq analysis. Pathway analysis indicated that the "carbohydrate metabolism pathway" and the "calcium signaling pathway" might play a vital role in tea plants' responses to cold stress. CONCLUSIONS: Our study presents a global survey of transcriptome profiles of tea plants in response to low, non-freezing temperatures and yields insights into the molecular mechanisms of tea plants during the cold acclimation process. It could also serve as a valuable resource for relevant research on cold-tolerance and help to explore the cold-related genes in improving the understanding of low-temperature tolerance and plant-environment interactions.

DeepEdit: single-molecule detection and phasing of A-to-I RNA editing events using nanopore direct RNA sequencing.

Single-molecule detection and phasing of A-to-I RNA editing events remain an unresolved problem. Long-read and PCR-free nanopore native RNA sequencing offers a great opportunity for direct RNA editing detection. Here, we develop a neural network model, DeepEdit, that not only recognizes A-to-I editing events in single reads of Oxford Nanopore direct RNA sequencing, but also resolves the phasing of RNA editing events on transcripts. We illustrate the robustness of DeepEdit by applying it to Schizosaccharomyces pombe and Homo sapiens transcriptome data. We anticipate DeepEdit to be a powerful tool for the study of RNA editing from a new perspective.

Critical assessment of DNA adenine methylation in eukaryotes using quantitative deconvolution.

The discovery of N6-methyldeoxyadenine (6mA) across eukaryotes led to a search for additional epigenetic mechanisms. However, some studies have highlighted confounding factors that challenge the prevalence of 6mA in eukaryotes. We developed a metagenomic method to quantitatively deconvolve 6mA events from a genomic DNA sample into species of interest, genomic regions, and sources of contamination. Applying this method, we observed high-resolution 6mA deposition in two protozoa. We found that commensal or soil bacteria explained the vast majority of 6mA in insect and plant samples. We found no evidence of high abundance of 6mA in Drosophila, Arabidopsis, or humans. Plasmids used for genetic manipulation, even those from Dam methyltransferase mutant Escherichia coli, could carry abundant 6mA, confounding the evaluation of candidate 6mA methyltransferases and demethylases. On the basis of this work, we advocate for a reassessment of 6mA in eukaryotes.

Optimizing de novo transcriptome assembly from short-read RNA-Seq data: a comparative study.

BACKGROUND: With the fast advances in nextgen sequencing technology, high-throughput RNA sequencing has emerged as a powerful and cost-effective way for transcriptome study. De novo assembly of transcripts provides an important solution to transcriptome analysis for organisms with no reference genome. However, there lacked understanding on how the different variables affected assembly outcomes, and there was no consensus on how to approach an optimal solution by selecting software tool and suitable strategy based on the properties of RNA-Seq data. RESULTS: To reveal the performance of different programs for transcriptome assembly, this work analyzed some important factors, including k-mer values, genome complexity, coverage depth, directional reads, etc. Seven program conditions, four single k-mer assemblers (SK: SOAPdenovo, ABySS, Oases and Trinity) and three multiple k-mer methods (MK: SOAPdenovo-MK, trans-ABySS and Oases-MK) were tested. While small and large k-mer values performed better for reconstructing lowly and highly expressed transcripts, respectively, MK strategy worked well for almost all ranges of expression quintiles. Among SK tools, Trinity performed well across various conditions but took the longest running time. Oases consumed the most memory whereas SOAPdenovo required the shortest runtime but worked poorly to reconstruct full-length CDS. ABySS showed some good balance between resource usage and quality of assemblies. CONCLUSIONS: Our work compared the performance of publicly available transcriptome assemblers, and analyzed important factors affecting de novo assembly. Some practical guidelines for transcript reconstruction from short-read RNA-Seq data were proposed. De novo assembly of C. sinensis transcriptome was greatly improved using some optimized methods.

Different adaptations of dopamine release in Nucleus Accumbens shell and core of individual alcohol drinking groups of mice.

Alcohol use disorder (AUD) places a tremendous burden on society, with approximately two billion alcohol users in the world. While most people drink alcohol recreationally, a subpopulation (3-5%) engages in reckless and compulsive drinking, leading to the development of AUD and alcohol dependence. The Ventral Tegmental Area (VTA)-Nucleus Accumbens (NAc) circuit has been shown to encode rewarding stimuli and drive individual alcohol drinking behavior. Our previous work successfully separated C57BL/6J isogenic mice into high or low alcohol drinking subgroups after a 12-day, two-bottle choice voluntary alcohol access paradigm. Electrophysiological studies revealed that low alcohol drinking mice exhibited elevated spontaneous and burst firing properties of their VTA dopamine (DA) neurons and specifically mimicking this pattern of activity in VTA-NAc neurons in high alcohol drinking mice using optogenetics decreased their alcohol preference. It is also known that VTA DA neurons encode the salience and rewarding properties of external stimuli while also regulating downstream dopamine concentrations. Here, as a follow-up to this study, we utilized Fast Scan Cyclic Voltammetry (FSCV) to examine dopamine release in the NAc shell and core between alcohol drinking groups. We observed dynamic changes of dopamine release in the core of high drinking mice, but failed to see widely significant differences of dopamine release in the shell of both groups, when compared with ethanol-naive controls. Overall, the present data suggest subregion-specific differences of evoked dopamine release in the NAc of low and high alcohol drinking mice, and may provide an anatomical substrate for individual alcohol drinking behavior. This article is part of the special issue on Stress, Addiction and Plasticity.

The evolutionary landscape of intergenic trans-splicing events in insects.

To explore the landscape of intergenic trans-splicing events and characterize their functions and evolutionary dynamics, we conduct a mega-data study of a phylogeny containing eight species across five orders of class Insecta, a model system spanning 400 million years of evolution. A total of 1,627 trans-splicing events involving 2,199 genes are identified, accounting for 1.58% of the total genes. Homology analysis reveals that mod(mdg4)-like trans-splicing is the only conserved event that is consistently observed in multiple species across two orders, which represents a unique case of functional diversification involving trans-splicing. Thus, evolutionarily its potential for generating proteins with novel function is not broadly utilized by insects. Furthermore, 146 non-mod trans-spliced transcripts are found to resemble canonical genes from different species. Trans-splicing preserving the function of 'breakup' genes may serve as a general mechanism for relaxing the constraints on gene structure, with profound implications for the evolution of genes and genomes.

The genome of Mesobuthus martensii reveals a unique adaptation model of arthropods.

Representing a basal branch of arachnids, scorpions are known as 'living fossils' that maintain an ancient anatomy and are adapted to have survived extreme climate changes. Here we report the genome sequence of Mesobuthus martensii, containing 32,016 protein-coding genes, the most among sequenced arthropods. Although M. martensii appears to evolve conservatively, it has a greater gene family turnover than the insects that have undergone diverse morphological and physiological changes, suggesting the decoupling of the molecular and morphological evolution in scorpions. Underlying the long-term adaptation of scorpions is the expansion of the gene families enriched in basic metabolic pathways, signalling pathways, neurotoxins and cytochrome P450, and the different dynamics of expansion between the shared and the scorpion lineage-specific gene families. Genomic and transcriptomic analyses further illustrate the important genetic features associated with prey, nocturnal behaviour, feeding and detoxification. The M. martensii genome reveals a unique adaptation model of arthropods, offering new insights into the genetic bases of the living fossils.