Adapting to stress: The effects of hibernation and hibernacula temperature on the hepatic transcriptome of Rhinolophus pusillus.

Hibernation, a survival strategy in mammals for extreme climates, induces physiological phenomena such as ischemia-reperfusion and metabolic shifts that hold great potential for advancements in modern medicine. Despite this, the molecular mechanisms underpinning hibernation remain largely unclear. This study used RNA-seq and Iso-seq techniques to investigate the changes in liver transcriptome expression of Rhinolophus pusillus during hibernation and active periods, as well as under different microhabitat temperatures. We identified 11 457 differentially expressed genes during hibernation and active periods, of which 395 showed significant differential expression. Genes associated with fatty acid catabolism were significantly upregulated during hibernation, whereas genes related to carbohydrate metabolism and glycogen synthesis were downregulated. Conversely, immune-related genes displayed differential expression patterns: genes tied to innate immunity were significantly upregulated, while those linked to adaptive immunity and inflammatory response were downregulated. The analysis of transcriptomic data obtained from different microhabitat temperatures revealed that R. pusillus exhibited an upregulation of genes associated with lipid metabolism in lower microhabitat temperature. This upregulation facilitated an enhanced utilization rate of triglyceride, ultimately resulting in increased energy provision for the organism. Additionally, R. pusillus upregulated gluconeogenesis-related genes regardless of the microhabitat temperature, demonstrating the importance of maintaining blood glucose levels during hibernation. Our transcriptomic data reveal that these changes in liver gene expression optimize energy allocation during hibernation, suggesting that liver tissue adaptively responds to the inherent stress of its function during hibernation. This study sheds light on the role of differential gene expression in promoting more efficient energy allocation during hibernation. It contributes to our understanding of how liver tissue adapts to the stressors associated with this state.

A comparative full-length transcriptomic resource provides insight into the perennial monocarpic mass flowering.

Perennial monocarpic mass flowering represents as a key developmental innovation in flowering time diversity in several biological and economical essential families, such as the woody bamboos and the shrubby Strobilanthes. However, molecular and genetic mechanisms underlying this important biodiversity remain poorly investigated. Here, we generated a full-length transcriptome resource incorporated into the BlueOmics database (http://blueomics.iflora.cn) for two Strobilanthes species, which feature contrasting flowering time behaviors. Using about 112 and 104 Gb Iso-seq reads together with ~185 and ~75 Gb strand-specific RNA seq data, we annotated 80 971 and 79 985 non-redundant full-length transcripts for the perennial polycarpic Strobilanthes tetrasperma and the perennial monocarpic Strobilanthes biocullata, respectively. In S. tetrasperma, we identified 8794 transcripts showing spatiotemporal expression in nine tissues. In leaves and shoot apical meristems at two developmental stages, 977 and 1121 transcripts were differentially accumulated in S. tetrasperma and S. biocullata, respectively. Interestingly, among the 33 transcription factors showing differential expression in S. tetrasperma but without differential expression in S. biocullata, three were involved potentially in the photoperiod and circadian-clock pathway of flowering time regulation (FAR1 RELATED SEQUENCE 12, FRS12; NUCLEAR FACTOR Y A1, NFYA1; PSEUDO-RESPONSE REGULATOR 5, PRR5), hence provides an important clue in deciphering the flowering diversity mechanisms. Our data serve as a key resource for further dissection of molecular and genetic mechanisms underpinning key biological innovations, here, the perennial monocarpic mass flowering.

Genome-wide characterization of post-transcriptional processes related to wood formation in Dalbergia odorifera.

BACKGROUND: Alternative polyadenylation (APA), alternative splicing (AS), and long non-coding RNAs (lncRNAs) play regulatory roles in post-transcriptional processes in plants. However, little is known about their involvement in xylem development in Dalbergia odorifera, a valuable rosewood species with medicinal and commercial significance. We addressed this by conducting Isoform Sequencing (Iso-Seq) using PacBio's SMRT technology and combined it with RNA-seq analysis (RNA sequencing on Illumina platform) after collecting xylem samples from the transition zone and the sapwood of D. odorifera. RESULTS: We identified 14,938 full-length transcripts, including 9,830 novel isoforms, which has updated the D. odorifera genome annotation. Our analysis has revealed that 4,164 genes undergo APA, whereas 3,084 genes encounter AS. We have also annotated 118 lncRNAs. Furthermore, RNA-seq analysis identified 170 differential alternative splicing (DAS) events, 344 genes with differential APA site usage (DE-APA), and 6 differentially expressed lncRNAs in the transition zone when compared to the sapwood. AS, APA, and lncRNAs are differentially regulated during xylem development. Differentially expressed APA genes were enriched for terpenoid and flavonoid metabolism, indicating their role in the heartwood formation. Additionally, DE-APA genes were associated with cell wall biosynthesis and terpenoid metabolism, implying an APA's role in wood formation. A DAS gene (involved in chalcone accumulation) with a significantly greater inclusion of the last exon in the transition zone than in the sapwood was identified. We also found that differentially expressed lncRNAs targeted the genes related to terpene synthesis. CONCLUSIONS: This study enhances our understanding of the molecular regulatory mechanisms underlying wood formation in D. odorifera, and provides valuable genetic resources and insights for its molecular-assisted breeding.

Identification of key genes for triacylglycerol biosynthesis and storage in herbaceous peony (Paeonia lactifolra Pall.) seeds based on full-length transcriptome.

BACKGROUND: The herbaceous peony (Paeonia lactiflora Pall.) is extensively cultivated in China due to its root being used as a traditional Chinese medicine known as 'Radix Paeoniae Alba'. In recent years, it has been discovered that its seeds incorporate abundant unsaturated fatty acids, thereby presenting a potential new oilseed plant. Surprisingly, little is known about the full-length transcriptome sequencing of Paeonia lactiflora, limiting research into its gene function and molecular mechanisms. RESULTS: A total of 484,931 Reads of Inserts (ROI) sequences and 1,455,771 full-Length non-chimeric reads (FLNC) sequences were obtained for CDS prediction, TF analysis, SSR analysis and lncRNA identification. In addition, gene function annotation and gene structure analysis were performed. A total of 4905 transcripts were related to lipid metabolism biosynthesis pathway, belonging to 28 enzymes. We use these data to identify 10 oleosin (OLE) and 5 diacylglycerol acyltransferase (DGAT) gene members after de-redundancy. The analysis of physicochemical properties and secondary structure showed them similarity in gene family respectively. The phylogenetic analysis showed that the distribution of OLE and DGAT family members was roughly the same as that of Arabidopsis. Quantitative real-time polymerase chain reaction (qRT-PCR) analyses revealed expression changes in different seed development stages, and showed a trend of increasing and then decreasing. CONCLUSION: In summary, these results provide new insights into the molecular mechanism of triacylglycerol (TAG) biosynthesis and storage during the seedling stage in Paeonia lactiflora. It provides theoretical references for selecting and breeding oil varieties and understanding the functions of oil storage as well as lipid synthesis related genes in Paeonia lactiflora.

Annotated genome and transcriptome of the endangered Caribbean mountainous star coral (Orbicella faveolata) using PacBio long-read sequencing.

Long-read sequencing is revolutionizing de-novo genome assemblies, with continued advancements making it more readily available for previously understudied, non-model organisms. Stony corals are one such example, with long-read de-novo genome assemblies now starting to be publicly available, opening the door for a wide array of 'omics-based research. Here we present a new de-novo genome assembly for the endangered Caribbean star coral, Orbicella faveolata, using PacBio circular consensus reads. Our genome assembly improved the contiguity (51 versus 1,933 contigs) and complete and single copy BUSCO orthologs (93.6% versus 85.3%, database metazoa_odb10), compared to the currently available reference genome generated using short-read methodologies. Our new de-novo assembled genome also showed comparable quality metrics to other coral long-read genomes. Telomeric repeat analysis identified putative chromosomes in our scaffolded assembly, with these repeats at either one, or both ends, of scaffolded contigs. We identified 32,172 protein coding genes in our assembly through use of long-read RNA sequencing (ISO-seq) of additional O. faveolata fragments exposed to a range of abiotic and biotic treatments, and publicly available short-read RNA-seq data. With anthropogenic influences heavily affecting O. faveolata, as well as its increasing incorporation into reef restoration activities, this updated genome resource can be used for population genomics and other 'omics analyses to aid in the conservation of this species.

Multiple strategies, including 6mA methylation, affecting plant alternative splicing in allopolyploid peanut.

Alternative splicing (AS), an important post-transcriptional regulation mechanism in eukaryotes, can significantly increase transcript diversity and contribute to gene expression regulation and many other complicated developmental processes. While plant gene AS events are well described, few studies have investigated the comprehensive regulation machinery of plant AS. Here, we use multi-omics to analyse peanut AS events. Using long-read isoform sequencing, 146 464 full-length non-chimeric transcripts were obtained, resulting in annotation corrections for 1782 genes and the identification of 4653 new loci. Using Iso-Seq RNA sequences, 271 776 unique splice junctions were identified, 82.49% of which were supported by transcriptome data. We characterized 50 977 polyadenylation sites for 23 262 genes, 12 369 of which had alternative polyadenylation sites. AS allows differential regulation of the same gene by miRNAs at the isoform level coupled with polyadenylation. In addition, we identified many long non-coding RNAs and fusion transcripts. There is a suppressed effect of 6mA on AS and gene expression. By analysis of chromatin structures, the genes located in the boundaries of topologically associated domains, proximal chromosomal telomere regions, inter- or intra-chromosomal loops were found to have more unique splice isoforms, higher expression, lower 6mA and more transposable elements (TEs) in their gene bodies than the other genes, indicating that chromatin interaction, 6mA and TEs play important roles in AS and gene expression. These results greatly refine the peanut genome annotation and contribute to the study of gene expression and regulation in peanuts. This work also showed AS is associated with multiple strategies for gene regulation.

Aberrant splicing of PSEN2, but not PSEN1, in individuals with sporadic Alzheimer's disease.

Alzheimer's disease is the most common neurodegenerative disease, characterized by dementia and premature death. Early-onset familial Alzheimer's disease is caused in part by pathogenic variants in presenilin 1 (PSEN1) and presenilin 2 (PSEN2), and alternative splicing of these two genes has been implicated in both familial and sporadic Alzheimer's disease. Here, we leveraged targeted isoform-sequencing to characterize thousands of complete PSEN1 and PSEN2 transcripts in the prefrontal cortex of individuals with sporadic Alzheimer's disease, familial Alzheimer's disease (carrying PSEN1 and PSEN2 variants), and controls. Our results reveal alternative splicing patterns of PSEN2 specific to sporadic Alzheimer's disease, including a human-specific cryptic exon present in intron 9 of PSEN2 as well as a 77 bp intron retention product before exon 6 that are both significantly elevated in sporadic Alzheimer's disease samples, alongside a significantly lower percentage of canonical full-length PSEN2 transcripts versus familial Alzheimer's disease samples and controls. Both alternatively spliced products are predicted to generate a prematurely truncated PSEN2 protein and were corroborated in an independent cerebellum RNA-sequencing dataset. In addition, our data in PSEN variant carriers is consistent with the hypothesis that PSEN1 and PSEN2 variants need to produce full-length but variant proteins to contribute to the onset of Alzheimer's disease, although intriguingly there were far fewer full-length transcripts carrying pathogenic alleles versus wild-type alleles in PSEN2 variant carriers. Finally, we identify frequent RNA editing at Alu elements present in an extended 3' untranslated region in PSEN2. Overall, this work expands the understanding of PSEN1 and PSEN2 variants in Alzheimer's disease, shows that transcript differences in PSEN2 may play a role in sporadic Alzheimer's disease, and suggests novel mechanisms of Alzheimer's disease pathogenesis.

BaRTv2: a highly resolved barley reference transcriptome for accurate transcript-specific RNA-seq quantification.

Accurate characterisation of splice junctions (SJs) as well as transcription start and end sites in reference transcriptomes allows precise quantification of transcripts from RNA-seq data, and enables detailed investigations of transcriptional and post-transcriptional regulation. Using novel computational methods and a combination of PacBio Iso-seq and Illumina short-read sequences from 20 diverse tissues and conditions, we generated a comprehensive and highly resolved barley reference transcript dataset from the European 2-row spring barley cultivar Barke (BaRTv2.18). Stringent and thorough filtering was carried out to maintain the quality and accuracy of the SJs and transcript start and end sites. BaRTv2.18 shows increased transcript diversity and completeness compared with an earlier version, BaRTv1.0. The accuracy of transcript level quantification, SJs and transcript start and end sites have been validated extensively using parallel technologies and analysis, including high-resolution reverse transcriptase-polymerase chain reaction and 5'-RACE. BaRTv2.18 contains 39 434 genes and 148 260 transcripts, representing the most comprehensive and resolved reference transcriptome in barley to date. It provides an important and high-quality resource for advanced transcriptomic analyses, including both transcriptional and post-transcriptional regulation, with exceptional resolution and precision.

Long read isoform sequencing reveals hidden transcriptional complexity between cattle subspecies.

The Iso-Seq method of full-length cDNA sequencing is suitable to quantify differentially expressed genes (DEGs), transcripts (DETs) and transcript usage (DTU). However, the higher cost of Iso-Seq relative to RNA-seq has limited the comparison of both methods. Transcript abundance estimated by RNA-seq and deep Iso-Seq data for fetal liver from two cattle subspecies were compared to evaluate concordance. Inter-sample correlation of gene- and transcript-level abundance was higher within technology than between technologies. Identification of DEGs between the cattle subspecies depended on sequencing method with only 44 genes identified by both that included 6 novel genes annotated by Iso-Seq. There was a pronounced difference between Iso-Seq and RNA-seq results at transcript-level wherein Iso-Seq revealed several magnitudes more transcript abundance and usage differences between subspecies. Factors influencing DEG identification included size selection during Iso-Seq library preparation, average transcript abundance, multi-mapping of RNA-seq reads to the reference genome, and overlapping coordinates of genes. Some DEGs called by RNA-seq alone appear to be sequence duplication artifacts. Among the 44 DEGs identified by both technologies some play a role in immune system, thyroid function and cell growth. Iso-Seq revealed hidden transcriptional complexity in DEGs, DETs and DTU genes between cattle subspecies previously missed by RNA-seq.

Transcriptomic analysis of testis and epididymis tissues from Banna mini-pig inbred line boars with single-molecule long-read sequencing†.

In mammals, testis and epididymis are critical components of the male reproductive system for androgen production, spermatogenesis, sperm transportation, as well as sperm maturation. Here, we report single-molecule real-time sequencing data from the testis and epididymis of the Banna mini-pig inbred line (BMI), a promising laboratory animal for medical research. We obtained high-quality full-length transcriptomes and identified 9879 isoforms and 8761 isoforms in the BMI testis and epididymis, respectively. Most of the isoforms we identified have novel exon structures that will greatly improve the annotation of testis- and epididymis-expressed genes in pigs. We also found that 3055 genes (over 50%) were shared between BMI testis and epididymis, indicating widespread expression profiles of genes related to reproduction. We characterized extensive alternative splicing events in BMI testis and epididymis and showed that 96 testis-expressed genes and 79 epididymis-expressed genes have more than six isoforms, revealing the complexity of alternative splicing. We accurately defined the transcribed isoforms in BMI testis and epididymis by combining Pacific Biotechnology Isoform-sequencing (PacBio Iso-Seq) and Illumina RNA Sequencing (RNA-seq) techniques. The refined annotation of some key genes governing male reproduction will facilitate further understanding of the molecular mechanisms underlying BMI male sterility. In addition, the high-confident identification of 548 and 669 long noncoding RNAs (lncRNAs) in these two tissues has established a candidate gene set for future functional investigations. Overall, our study provides new insights into the role of the testis and epididymis during BMI reproduction, paving the path for further studies on BMI male infertility.

Insights on the stem elongation of spur-type bud sport mutant of 'Red Delicious' apple.

MAIN CONCLUSION: The decreased capacity of auxin-, CTK-, and BR-mediated cell division and cell enlargement pathways, combined with the enhanced capacity of GA and ETH-, JA-, ABA-, SA-mediated stress-resistant pathways were presumed to be the crucial reasons for the formation of spur-type 'Red Delicious' mutants. Vallee Spur', which exhibit short internodes and compact tree shape, is the fourth generation of the spur-type bud sport mutant of 'Red Delicious'. However, the underlying molecular mechanism of these properties remains unclear. Here, comparative phenotypic, full-length transcriptome and phytohormone analyses were performed between 'Red Delicious' (NSP) and 'Vallee Spur' (SP). The new shoot internode length of NSP was ˃ 1.53-fold higher than that of the SP mutant. Cytological analysis showed that the stem cells of the SP mutant were smaller and more tightly arranged relative to the NSP. By Iso-Seq, a total of 1426 differentially expressed genes (DEGs) were detected, including 808 upregulated and 618 downregulated genes in new shoot apex with 2 leaves of the SP mutant. Gene expressions involved in auxin, cytokinin (CTK), and brassinosteroid (BR) signal transduction were mostly downregulated in the SP mutant, whereas those involved in gibberellin (GA), ethylene (ETH), jasmonate (JA), ABA, and salicylic acid (SA) signal transduction were mostly upregulated. The overall thermogram analysis of hormone levels in the shoot apex carrying two leaves detected by LC-MS/MS absolute quantification showed that the levels of IAA-Asp, IAA, iP7G, OPDA, and 6-deoxyCS were significantly upregulated in the SP mutant, while the remaining 28 hormones were significantly downregulated. It is speculated that the decreased capacity of auxin, CTK, and BR-mediated cell division and cell enlargement pathways is crucial for the formation of the SP mutant. GA and stress-resistant pathways of ETH, JA, ABA, and SA also play vital roles in stem elongation. These results highlight the involvement of phytohormones in the formation of stem elongation occurring in 'Red Delicious' spur-type bud sport mutants and provide information for exploring its biological mechanism.

Polyadenylation-related isoform switching in human evolution revealed by full-length transcript structure.

Rhesus macaque is a unique nonhuman primate model for human evolutionary and translational study, but the error-prone gene models critically limit its applications. Here, we de novo defined full-length macaque gene models based on single molecule, long-read transcriptome sequencing in four macaque tissues (frontal cortex, cerebellum, heart and testis). Overall, 8 588 227 poly(A)-bearing complementary DNA reads with a mean length of 14 106 nt were generated to compile the backbone of macaque transcripts, with the fine-scale structures further refined by RNA sequencing and cap analysis gene expression sequencing data. In total, 51 605 macaque gene models were accurately defined, covering 89.7% of macaque or 75.7% of human orthologous genes. Based on the full-length gene models, we performed a human-macaque comparative analysis on polyadenylation (PA) regulation. Using macaque and mouse as outgroup species, we identified 79 distal PA events newly originated in humans and found that the strengthening of the distal PA sites, rather than the weakening of the proximal sites, predominantly contributes to the origination of these human-specific isoforms. Notably, these isoforms are selectively constrained in general and contribute to the temporospatially specific reduction of gene expression, through the tinkering of previously existed mechanisms of nuclear retention and microRNA (miRNA) regulation. Overall, the protocol and resource highlight the application of bioinformatics in integrating multilayer genomics data to provide an intact reference for model animal studies, and the isoform switching detected may constitute a hitherto underestimated regulatory layer in shaping the human-specific transcriptome and phenotypic changes.

Transcriptomic evidence for Brachyury expression in the caudal tip region of adult Ptychodera flava (Hemichordata).

Most metazoans have a single copy of the T-box transcription factor gene Brachyury. This gene is expressed in cells of the blastopore of late blastulae and the archenteron invagination region of gastrulae. It appears to be crucial for gastrulation and mesoderm differentiation of embryos. Although this expression pattern is shared by most deuterostomes, Brachyury expression has not been reported in adult stages. Here we show that Brachyury of an indirect developer, the hemichordate acorn worm Ptychodera flava, is expressed not only in embryonic cells, but also in cells of the caudal tip (anus) region of adults. This spatially restricted expression, shown by whole-mount in situ hybridization, was confirmed by Iso-Seq RNA sequencing and single-cell RNA-seq (scRNA-seq) analysis. Iso-Seq analysis showed that gene expression occurs only in the caudal region of adults, but not in anterior regions, including the stomochord. scRNA-seq analysis showed a cluster that contained Brachyury-expressing cells comprising epidermis- and mesoderm-related cells, but which is unlikely to be associated with the nervous system or muscle. Although further investigation is required to examine the roles of Brachyury in adults, this study provides important clues for extending studies on Brachyury expression involved in development of the most posterior region of deuterostomes.

Characterization of T-cell receptors and immunoglobulin heavy chains loci and identification of T/B cell clusters in teleost.

Bacterial disease is one of the important factors leading to economic losses in the turbot (Scophthalmus maximus) cultivation industry. T lymphocytes are major components of cellular immunity, whereas B lymphocytes produce immunoglobulins (Ig) that are key elements of humoral immune responses against infection. However, the genomic organization of genes encoding T-cell receptors (TCR) and immunoglobulin heavy chains (IgHs) in turbot remains largely unknown. In this study, abundant full-length transcripts of TCRs and IgHs were sequenced by Isoform-sequencing (Iso-seq), and we investigated and annotated the V, D, J and C gene loci of TCRα, TCRß, IgT, IgM and IgD in turbot. Furthermore, through single-cell RNA sequencing (scRNA-seq) of blood leukocytes, we confirmed that these identified TCRs and IgHs were highly expressed in T/B cell clusters, respectively. Meanwhile, we also identified the IgM+IgD+ B and IgT+ B cells with differential gene expression profiles and potential functions. Taken together, our results provide a comprehensive understanding of TCRs and IgHs loci in turbot, which will contribute to evolutionary and functional characterization of T and B lymphocytes in teleost.

Full-length transcriptome sequencing reveals extreme incomplete annotation of the goat genome.

Despite recent advances in generating high-quality reference genome assemblies, the genome sequences for most livestock species, including goats, are still poorly annotated. Single-molecule long-read sequencing has greatly facilitated gene annotation by obtaining full-length transcripts. In this study, we generated full-length transcriptome data for samples from abomasum (n = 2) and testicle (n = 1), using PacBio Iso-Seq technology. We further combined these data with published data from abomasum (5ZY, SRR8618141) to evaluate and improve the gene annotation of the goat genome. We identified 14.5-16.3% of novel genes per sample from the four Iso-Seq datasets. At the transcript level, 40.6% of them were novel, including 29.7% novel transcripts from known genes and 10.9% from novel genes. We further verified the expression of novel genes in four additional RNA-seq data and found that the expression level of novel genes was significantly lower than that of known genes, indicating that the lowly expressed genes tend to be missed in the current genome annotation. This study shows the superiority of full-length transcriptome data in gene annotation, and more such data are required to improve the gene annotation for goat genome and other species.

Aethionema arabicum genome annotation using PacBio full-length transcripts provides a valuable resource for seed dormancy and Brassicaceae evolution research.

Aethionema arabicum is an important model plant for Brassicaceae trait evolution, particularly of seed (development, regulation, germination, dormancy) and fruit (development, dehiscence mechanisms) characters. Its genome assembly was recently improved but the gene annotation was not updated. Here, we improved the Ae. arabicum gene annotation using 294 RNA-seq libraries and 136 307 full-length PacBio Iso-seq transcripts, increasing BUSCO completeness by 11.6% and featuring 5606 additional genes. Analysis of orthologs showed a lower number of genes in Ae. arabicum than in other Brassicaceae, which could be partially explained by loss of homeologs derived from the At-α polyploidization event and by a lower occurrence of tandem duplications after divergence of Aethionema from the other Brassicaceae. Benchmarking of MADS-box genes identified orthologs of FUL and AGL79 not found in previous versions. Analysis of full-length transcripts related to ABA-mediated seed dormancy discovered a conserved isoform of PIF6-ß and antisense transcripts in ABI3, ABI4 and DOG1, among other cases found of different alternative splicing between Turkey and Cyprus ecotypes. The presented data allow alternative splicing mining and proposition of numerous hypotheses to research evolution and functional genomics. Annotation data and sequences are available at the Ae. arabicum DB (https://plantcode.online.uni-marburg.de/aetar_db).

A survey of transcriptome complexity using full-length isoform sequencing in the tea plant Camellia sinensis.

Tea is one of the most popular beverages and its leaves are rich in catechins, contributing to the diverse flavor as well as beneficial for human health. However, the study of the post-transcriptional regulatory mechanism affecting the synthesis of catechins remains insufficient. Here, we sequenced the transcriptome using PacBio sequencing technology and obtained 63,111 full-length high-quality isoforms, including 1302 potential novel genes and 583 highly reliable fusion transcripts. We also identified 1204 lncRNAs with high quality, containing 188 known and 1016 novel lncRNAs. In addition, 311 mis-annotated genes were corrected based on the high-quality Isoseq reads. A large number of alternative splicing (AS) events (3784) and alternative polyadenylation (APA) genes (18,714) were analyzed, accounting for 8.84% and 43.7% of the total annotated genes, respectively. We also found that 2884 genes containing AS and APA features exhibited higher expression levels than other genes. These genes are mainly involved in amino acid biosynthesis, carbon fixation in photosynthetic organisms, phenylalanine, tyrosine, tryptophan biosynthesis, and pyruvate metabolism, suggesting that they play an essential role in the catechins content of tea polyphenols. Our results further improved the level of genome annotation and indicated that post-transcriptional regulation plays a crucial part in synthesizing catechins.

Landscape and regulation of alternative splicing and alternative polyadenylation in a plant pathogenic fungus.

Alternative splicing (AS) and alternative polyadenylation (APA) contribute significantly to the regulation of gene expression in higher eukaryotes. Their biological impact in filamentous fungi, however, is largely unknown. Here we combine PacBio Isoform-Sequencing and strand-specific RNA-sequencing of multiple tissues and mutant characterization to reveal the landscape and regulation of AS and APA in Fusarium graminearum. We generated a transcript annotation comprising 51 617 isoforms from 17 189 genes. In total, 4997 and 11 133 genes are alternatively spliced and polyadenylated, respectively. Majority of the AS events alter coding sequences. Unexpectedly, the AS transcripts containing premature-termination codons are not sensitive to nonsense-mediated messenger RNA decay. Unlike in yeasts and animals, distal APA sites have strong signals, but proximal APA isoforms are highly expressed in F. graminearum. The 3'-end processing factors FgRNA15, FgHRP1, and FgFIP1 play roles in promoting proximal APA site usage and intron splicing. A genome-wide increase in intron inclusion and distal APA site usage and downregulation of the spliceosomal and 3'-end processing factors were observed in older and quiescent tissues, indicating intron inclusion and 3'-untranslated region lengthening as novel mechanisms in regulating aging and dormancy in fungi. This study provides new insights into the complexity and regulation of AS and APA in filamentous fungi.

The importance of alternative splicing in adaptive evolution.

Although alternative splicing is a ubiquitous co-transcriptional gene regulatory mechanism in plants, animals and fungi, its contribution to evolutionary transitions is understudied. Alternative splicing enables different mRNA isoforms to be generated from the same gene, expanding transcriptomic and thus proteomic diversity. While the role of gene expression variation in adaptive evolution is widely accepted, biologists still debate the functional impact of alternative isoforms on phenotype. In light of recent empirical research linking splice variation to ecological adaptations, we propose that alternative splicing is an important substrate for adaptive evolution and speciation, particularly at short timescales. In this article we synthesise what is known about the role of alternative splicing in adaptive evolution. We discuss the contribution of standing splice variation to phenotypic plasticity and how hybridisation can produce novel splice forms. Going forwards, we propose that alternative splicing be included as a standard analysis alongside gene expression analysis so we can better understand of how alternative splicing contributes to adaptive divergence at the micro- and macroevolutionary levels.

Global identification of full-length cassava lncRNAs unveils the role of cold-responsive intergenic lncRNA 1 in cold stress response.

Long noncoding RNAs (lncRNAs) have been considered to be important regulators of gene expression in a range of biological processes in plants. A large number of lncRNAs have been identified in plants. However, most of their biological functions still remain to be determined. Here, we identified a total of 3004 lncRNAs in cassava under normal or cold-treated conditions from Iso-seq data. We further characterized a cold-responsive intergenic lncRNA 1 (CRIR1) as a novel positive regulator of the plant response to cold stress. CRIR1 can be significantly induced by cold treatment. Ectopic expression of CRIR1 in cassava enhanced the cold tolerance of transgenic plants. Transcriptome analysis demonstrated that CRIR1 regulated a range of cold stress-related genes in a CBF-independent pathway. We further found that CRIR1 RNA can interact with cassava cold shock protein 5 (MeCSP5), which acts as an RNA chaperone, indicating that CRIR1 may recruit MeCSP5 to improve the translation efficiency of messenger RNA. In summary, our study extends the repertoire of lncRNAs in plants as well as their role in cold stress responses. Moreover, it reveals a mechanism by which CRIR1 affected cold stress response by modulating the expression of stress-responsive genes and increasing their translational yield.