Full-length Transcriptome Sequencing Research Articles

Using long-read CAGE sequencing to profile cryptic-promoter-derived transcripts and their contribution to the immunopeptidome.

Recent studies have shown that the noncoding genome can produce unannotated proteins as antigens that induce immune response. One major source of this activity is the aberrant epigenetic reactivation of transposable elements (TEs). In tumors, TEs often provide cryptic or alternate promoters, which can generate transcripts that encode tumor-specific unannotated proteins. Thus, TE-derived transcripts (TE transcripts) have the potential to produce tumor-specific, but recurrent, antigens shared among many tumors. Identification of TE-derived tumor antigens holds the promise to improve cancer immunotherapy approaches; however, current genomics and computational tools are not optimized for their detection. Here we combined CAGE technology with full-length long-read transcriptome sequencing (long-read CAGE, or LRCAGE) and developed a suite of computational tools to significantly improve immunopeptidome detection by incorporating TE and other tumor transcripts into the proteome database. By applying our methods to human lung cancer cell line H1299 data, we show that long-read technology significantly improves mapping of promoters with low mappability scores and that LRCAGE guarantees accurate construction of uncharacterized 5' transcript structure. Augmenting a reference proteome database with newly characterized transcripts enabled us to detect noncanonical antigens from HLA-pulldown LC-MS/MS data. Lastly, we show that epigenetic treatment increased the number of noncanonical antigens, particularly those encoded by TE transcripts, which might expand the pool of targetable antigens for cancers with low mutational burden.

Comprehensive Long-Read Sequencing Analysis Discloses the Transcriptome Features of Papillary Thyroid Microcarcinoma.

Papillary thyroid microcarcinoma (PTMC) is the most common type of thyroid cancer. It has been shown that lymph node metastasis is associated with poor prognosis in patients with PTMC. We aim to characterize the PTMC transcriptome landscape and identify the candidate transcripts that are associated with lateral neck lymph node metastasis of PTMC. We performed full-length transcriptome sequencing in 64 PTMC samples. Standard bioinformatic pipelines were applied to characterize and annotate the full-length expression profiles of 2 PTMC subtypes. Functional open reading frame (ORF) annotation of the known and novel transcripts were predicted by HMMER, DeepLoc, and DeepTMHMM tools. Candidate transcripts associated with the pN1b subtype were identified after transcript quantification and differential gene expression analyses. We found that skipping exons accounted for the more than 27.82% of the alternative splicing events. At least 42.56% of the discovered transcripts were novel isoforms of annotated genes. A total of 39 193 ORFs in novel transcripts and 18 596 ORFs in known transcripts were identified. Distribution patterns of the characterized transcripts in functional domain, subcellular localization, and transmembrane structure were predicted. In total, 1033 and 1204 differentially expressed genes were identified in the pN0 and pN1b groups, respectively. Moreover, novel isoforms of FRMD3, NOD1, and SHROOM4 were highlighted for their association with pN1b subtype. Our data provided the global transcriptome landscape of PTMC and also revealed the novel isoforms that associated with PTMC aggressiveness.

PacBio Full-Length Transcriptome Sequencing Reveals the Mechanism of Salt Stress Response in Sonneratia apetala.

Sonneratia apetala is an essential mangrove wetland restoration tree species. Studying its molecular mechanism for salt tolerance could lay a foundation for further cultivating excellent resistant germplasm. This study used a combination of PacBio isoform sequencing (Iso-seq) and BGISEQ RNA sequencing (RNA-seq) to analyze the molecular mechanism to salt stress response of one-year-old S. apetala leaves. The growth and physiological analysis showed that physiological indexes such as growth rate, net photosynthetic rate and antioxidant enzyme activity all exhibit significant changes under salt stress. From Iso-seq, a total of 295,501 full-length transcripts, with an average length of 1418 bp, were obtained. RNA-seq produced 4712 differentially expressed genes (DEGs) as compared to a control group. Of these, 930 were identified to be co-expressed during the STEM time sequence analysis. Further, 715 and 444 co-expressed DEGs were annotated by GO and KEGG analyses, respectively. Moreover, 318 of the co-expressed DEGs were annotated as essential genes that were implicated in salt stress response of S. apetala, which were involved in transcription factors, signal transduction, hormone response, ROS homeostasis, osmotic balance, cell wall synthesis or modification. These results provide candidate targets for further characterization and offer insights into the salt-tolerant mechanism of S. apetala.

Full-Length Transcriptome Sequencing Analysis and Characterization of WRKY Transcription Factors Responsive to Cadmium Stress in Arabis paniculata.

Arabis paniculata is a newly discovered hyperaccumulator known for its ability to accumulate multiple metals. WRKY proteins play a significant role in plant responses to various stresses, including cadmium (Cd) stress. However, there is limited research on the molecular biology of Arabis paniculata, especially regarding the WRKY family. In this study, we conducted third-generation sequencing for functional annotation and structural analysis of Arabis paniculata. We obtained 41,196 high-quality isoforms from the full-length transcriptome, with an average length of 1043 bp. A total of 26,670 genes were predicted against NR, Swissprot, KOG, and KEGG databases. Functional comparison using the KOG database revealed excellent annotation in 25 functional categories, with general function prediction (1822 items) being the most predominant. MISA analysis identified 12,593 SSR loci, with single nucleotide repeats being the largest category (44.83% of the total). Moreover, our predictions provide insights into 20,022 coding sequences (CDS), 811 transcription factors, and 17,963 LncRNAs. In total, 34 WRKY gene sequences were identified in Arabis paniculata. Bioinformatics analysis revealed diverse numbers of amino acids in these WRKYs (113 to 545 aa), and a conserved WRKYGQK sequence within the N-terminus of the WRKY protein. Furthermore, all WRKYs were found to be localized in the nucleus. Phylogenetic analysis classified the WRKY genes into three categories: I (14 members), II (17 members), and III (3 members). Category II was subsequently divided into four sub-categories: II-a (8 members), II-b (1 member), II-c (1 member), and II-d (7 members). Our quantitative real-time polymerase chain reaction (qRT-PCR) experiments revealed that ApWRKY23 and ApWRKY34 exhibited the highest expression levels at the 24-h time point, suggesting their potential role as the candidate genes for Cd stress response. These findings contribute to our understanding of the genomic information of Arabis paniculata and provide a basis for the analysis of its genetic diversity. Additionally, this study paves the way for a comprehensive exploration of the molecular mechanisms underlying the WRKY genes in Arabis paniculata under Cd stress conditions.

Deciphering triterpenoid saponin biosynthesis by leveraging transcriptome response to methyl jasmonate elicitation in Saponaria vaccaria

Methyl jasmonate (MeJA) is a known elicitor of plant specialized metabolism, including triterpenoid saponins. Saponaria vaccaria is an annual herb used in traditional Chinese medicine, containing large quantities of oleanane-type triterpenoid saponins with anticancer properties and structural similarities to the vaccine adjuvant QS-21. Leveraging the MeJA-elicited saponin biosynthesis, we identify multiple enzymes catalyzing the oxidation and glycosylation of triterpenoids in S. vaccaria. This exploration is aided by Pacbio full-length transcriptome sequencing and gene expression analysis. A cellulose synthase-like enzyme can not only glucuronidate triterpenoid aglycones but also alter the product profile of a cytochrome P450 monooxygenase via preference for the aldehyde intermediate. Furthermore, the discovery of a UDP-glucose 4,6-dehydratase and a UDP-4-keto-6-deoxy-glucose reductase reveals the biosynthetic pathway for the rare nucleotide sugar UDP-d-fucose, a likely sugar donor for fucosylation of plant natural products. Our work enables the production and optimization of high-value saponins in microorganisms and plants through synthetic biology approaches.

Full-length transcriptome sequencing provides new insights into the complexity of flavonoid biosynthesis in Glechoma longituba.

Glechoma longituba has been frequently used in treating urolithiasis and cholelithiasis due to the presence of flavonoids, which are its major bioactive constituents. However, research on the molecular background of flavonoid biosynthesis in G. longituba is limited. In this study, we used single-molecule real-time combined with next-generation sequencing technologies to construct the complete transcriptome of G. longituba. We identified 404,648 non-redundant transcripts, including 249,697 coding sequences, 197,811 simple sequence repeats, 174,846 long noncoding RNA, and 176,554 coding RNA. Moreover, we functionally annotated 346,218 isoforms (85.56%) and identified 86,528 differentially expressed genes. We also identified 55 non-redundant full-length isoforms related to the flavonoid biosynthetic pathway. Pearson correlation analysis revealed that the expression levels of some key genes of the flavonoid biosynthesis pathway were significantly positively correlated with the flavonoid metabolites. Furthermore, we performed bioinformatics analysis (sequence and structural) of isoform_47029 (encoding flavanone 3-hydroxylase) and isoform_53692 (encoding flavonol synthase) to evaluate their potential biological functions. Finally, we validated gene expression levels of 12 flavonoid-related key enzyme genes using quantitative real-time PCR. Overall, this study provides full-length transcriptome information on G. longituba for the first time and valuable molecular resources for further research on the medicinal properties of this plant.

Genome-Wide Identification and Expression Analysis of the MADS Gene Family in Tulips (Tulipa gesneriana).

To investigate the cold response mechanism and low temperature regulation of flowering in tulips, this study identified 32 MADS-box transcription factor family members in tulips based on full-length transcriptome sequencing, named TgMADS1-TgMADS32. Phylogenetic analysis revealed that these genes can be divided into two classes: type I and type II. Structural analysis showed that TgMADS genes from different subfamilies have a similar distribution of conserved motifs. Quantitative real-time PCR results demonstrated that some TgMADS genes (e.g., TgMADS3, TgMADS15, TgMADS16, and TgMADS19) were significantly upregulated in buds and stems under cold conditions, implying their potential involvement in the cold response of tulips. In summary, this study systematically identified MADS family members in tulips and elucidated their evolutionary relationships, gene structures, and cold-responsive expression patterns, laying the foundation for further elucidating the roles of these transcription factors in flowering and the cold adaptability of tulips.

Full-Length Transcriptome Sequencing and Identification of Genes Related to Terpenoid Biosynthesis in Cinnamomum migao H. W. Li

Cinnamomum migao H. W. Li is an evergreen woody plant that is only distributed in southwest China. The volatile oils from the fruits of C. migao have long been used as a special authentic medicinal herb by local ethnic minorities. Due to its low seed germination rate, destructive lumbering and low rates of artificial planting, C. migao is facing the danger of extinction. Therefore, it is urgent to exploit and protect this species using molecular biological technology, especially to target the genes involved in the biosynthesis of terpenoids in the volatile oil. However, the genomic data for this plant are not available. In this study, the transcriptome sequencing of C. migao was performed to obtain the key genes involved in terpenoid biosynthesis via a combination of full-length transcriptome and next-generation sequencing (NGS). More than 39.9 Gb of raw data was obtained and 515,929 circular consensus sequences (CCSs) were extracted. After clustering similar 472,858 full-length non-chimeric (FLNC) reads and correction with NGS data, 139,871 consensus isoforms were obtained. Meanwhile, 73,575 non-redundant transcripts were generated by removing redundant transcripts. Then, 70,427 isoforms were successfully annotated using public databases. Moreover, differentially expressed transcripts (DETs) in four different developmental stages of the C. migao fruit were analyzed and 5764 transcripts showed period-specific expression. Finally, 15 transcripts, 6 transcripts, and only 1 transcript were identified as being involved in the biosynthesis of sesquiterpenoids, diterpenoids, and monoterpenoids, respectively. This study provides a basis for future research in the gene mining, genetic breeding, and metabolic engineering of C. migao.

The A1762T/G1764A mutations enhance HBV replication by alternating viral transcriptome.

The A1762T/G1764A mutations, one of the most common mutations in the hepatitis B virus basal core promoter, are associated with the progression of chronic HBV infection. However, effects of these mutations on HBV replication remains controversial. This study aimed to systematically investigate the effect of the mutations on HBV replication and its underlying mechanisms. Using the prcccDNA/pCMV-Cre recombinant plasmid system, a prcccDNA-A1762T/G1764A mutant plasmid was constructed. Compared with wild-type HBV, A1762T/G1764A mutant HBV showed enhanced replication ability with higher secreted HBV DNA and RNA levels, while Southern and Northern blot indicated higher intracellular levels of relaxed circular DNA, single-stranded DNA, and 3.5 kb RNA. Meanwhile, the mutations increased expression of intracellular core protein and decreased the production of HBeAg and HBsAg. In vitro infection based on HepG2-NTCP cells and mice hydrodynamic injection experiment also proved that these mutations promote HBV replication. 5'-RACE assays showed that these mutations upregulated transcription of pregenomic RNA (pgRNA) while downregulating that of preC RNA, which was further confirmed by full-length transcriptome sequencing. Moreover, a proportion of sub-pgRNAs with the potential to express polymerase were also upregulated by these mutations. The ChIP-qPCR assay showed that A1762T/G1764A mutations created a functional HNF1α binding site in the BCP region, and its overexpression enhanced the effect of A1762T/G1764A mutations on HBV. Our findings revealed the mechanism and importance of A1762T/G1764A mutations as an indicator for management of CHB patients, and provided HNF1α as a new target for curing HBV-infected patients.

Physiological and transcriptomic responses of Pinellia ternata to continuous cropping

Pinellia ternata (Thunb.) Breit. is one of the popular medicinal plants and industrial crops in China, but its production has been severely limited by continuous cropping obstacles (CCO). In order to elucidate the reaction of P. ternata to successive cultivation and the potentially uncover the possible causes of CCO, the research group examined P. ternata relevant indexes, including morphological, physio-biochemical, as well as its transcriptome, under continuous cropping (CC) conditions. Based on the results, the influences of CC were multilevel. Specifically, CC led to chlorophyll content in leaves and antioxidant enzymes (peroxidase and superoxide dismutase, i.e. POD & SOD) activities significant reduced, reactive oxygen species (ROS) and malondialdehyde (MDA) in cells increased significantly. And the obvious changes were also observed in osmoregulatory substances such as proline (Pro), soluble protein (SP) and soluble sugars (SS). As a result, the plants growth was hindered and had the appearance of a "sprout tumble" prematurely. By combining SMRT full-length and Illumina transcriptome sequencing, a database of the root and leaf transcriptome of P. ternata at different years of CC years was established, and 77,783 full-length non-redundant transcripts were assembled. 412 differentially expressed genes (DEGs) were identified in the root during the first planting year (FPR) vs. the root after 2 years of Continuous Cropping (CC2R), of which 138 were up-regulated, meanwhile 274 were down-regulated. Moreover, 5020 DEGs were identified in the leaves during the first planting year (FPL) vs. the leaves after 2 years of Continuous Cropping (CC2L), of which 1894 were up-regulated, and 3126 were down-regulated. The analysis of DEGs showed that photosynthesis and chlorophyll biosynthesis, fatty acid and wax biosynthesis, cell wall synthesis and modification, defense response and plant-pathogen interaction, as well as phytohormone signaling were the key genes differentially expressed. Weighted gene correlation network analysis (WGCNA) analysis revealed that the hub genes in the roots were enriched in pathways associated with sugar synthesis and metabolism, phosphate ion transport, jasmonate metabolism, and fungal defense. At the same time, phosphate ion transport, positive regulation of cellular response to phosphate starvation, chlorophyll biosynthetic process, photosynthesis, and stress response pathways were enriched in the leaves. As a result, it was inferred that phosphorus depletion and the alteration of the soil microbiota by the previous planting cycles, which resulted in nutrient element deficiencies, as well as the damage of roots by pathogenic fungi, were the main CCO causes of P. ternata, and the corresponding changes in the physiological and transcriptional homeostasis of the leaves due to a generalized stress response.

Targeting hnRNPC suppresses thyroid follicular epithelial cell apoptosis and necroptosis through m6A-modified ATF4 in autoimmune thyroid disease

Both environmental and genetic factors contribute to the etiology of autoimmune thyroid disease (AITD) including Graves’ disease (GD) and Hashimoto’s thyroiditis (HT). However, the exact pathogenesis and interactions that occur between environmental factors and genes remain unclear, and therapeutic targets require further investigation due to limited therapeutic options. To solve such problems, this study utilized single-cell transcriptome, whole transcriptome, full-length transcriptome (Oxford nanopore technology), and metabolome sequencing to examine thyroid lesion tissues from 2 HT patients and 2 GD patients as well as healthy thyroid tissue from 1 control subject. HT patients had increased ATF4-positive thyroid follicular epithelial (ThyFoEp) cells, which significantly increased endoplasmic reticulum stress. The enhanced sustained stress resulted in cell death mainly including apoptosis and necroptosis. The ATF4-based global gene regulatory network and experimental validation revealed that N6-methyladenosine (m6A) reader hnRNPC promoted the transcriptional activity, synthesis, and translation of ATF4 through mediating m6A modification of ATF4. Increased ATF4 expression initiated endoplasmic reticulum stress signaling, which when sustained, caused apoptosis and necroptosis in ThyFoEp cells, and mediated HT development. Targeting hnRNPC and ATF4 notably decreased ThyFoEp cell death, thus ameliorating disease progression. Collectively, this study reveals the mechanisms by which microenvironmental cells in HT and GD patients trigger and amplify the thyroid autoimmune cascade response. Furthermore, we identify new therapeutic targets for the treatment of autoimmune thyroid disease, hoping to provide a potential way for targeted therapy.

Differentially expressed transcripts study during pregnancy and postpartum anestrus of yak (Bos grunniens)

Background: Yak is the main livestock species in the plateau area, and its reproductive performance is low, usually two years or three years. A very few of yaks recover within a certain period of time after delivery and smoothly enter the next estrous cycle, while most of them enter the postpartum anestrus and show no estrus performance. However, the key biological factors and influencing mechanisms that cause postpartum anestrus in yaks are not clear. Objective: To study the expression of differential transcripts in ovaries of yak during pregnancy and postpartum anestrus. Methods: Each three yaks in pregnancy and anestrus under natural grazing conditions in Haiyan County, Qinghai Province were selected and slaughtered, and their ovaries were collected and sent to Biomarker Technologies. Oxford Nanopore Technologies single-molecule real-time electrical signal sequencing technology was used to perform full-length transcriptome sequencing. Astalavista software was used to identify the types of alternative splicing events in yak estrus and pregnancy, and TAPIS pipeline was used to identify alternative polyadenylation. Results: The results showed that there were 1751 differentially expressed transcripts (DETs) between pregnancy and anestrus in yak, of which 808 were upregulated and 943 were downregulated. GO analysis showed that the biological processes of DETs were mainly reproductive, reproductive and rhythmic processes. KEGG analysis showed that the DET cell junction-related adhesion junction protein (β-catenin) and amino terminal kinase (JNK) were involved in FAs (local adhesion). Phosphatidylinositol-3-kinase (PI3K) is involved in the PI3K/AKT/mTOR signaling pathway. Circadian rhythm output cycle failure (Clock) and brain and muscle tissue aromatic hydrocarbon receptor nuclear transporter-like protein 1 (Bmal1) are involved in circadian rhythm signaling pathway. Conclusion: This study found that β-catenin, JNK, PI3K, Clock and Bmal1 were closely related to postpartum anestrus in yak.

Integrating Full-Length Transcriptome and RNA Sequencing of Siberian Wildrye (Elymus sibiricus) to Reveal Molecular Mechanisms in Response to Drought Stress.

Drought is one of the most significant limiting factors affecting plant growth and development on the Qinghai-Tibet Plateau (QTP). Mining the drought-tolerant genes of the endemic perennial grass of the QTP, Siberian wildrye (Elymus sibiricus), is of great significance to creating new drought-resistant varieties which can be used in the development of grassland livestock and restoring natural grassland projects in the QTP. To investigate the transcriptomic responsiveness of E. sibiricus to drought stress, PEG-induced short- and long-term drought stress was applied to two Siberian wildrye genotypes (drought-tolerant and drought-sensitive accessions), followed by third- and second-generation transcriptome sequencing analysis. A total of 40,708 isoforms were detected, of which 10,659 differentially expressed genes (DEGs) were common to both genotypes. There were 2107 and 2498 unique DEGs in the drought-tolerant and drought-sensitive genotypes, respectively. Additionally, 2798 and 1850 DEGs were identified in the drought-tolerant genotype only under short- and long-term conditions, respectively. DEGs numbering 1641 and 1330 were identified in the drought-sensitive genotype only under short- and long-term conditions, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that all the DEGs responding to drought stress in E. sibiricus were mainly associated with the mitogen-activated protein kinase (MAKP) signaling pathway, plant hormone signal transduction, the linoleic acid metabolism pathway, the ribosome pathway, and plant circadian rhythms. In addition, Nitrate transporter 1/Peptide transporter family protein 3.1 (NPF3.1) and Auxin/Indole-3-Acetic Acid (Aux/IAA) family protein 31(IAA31) also played an important role in helping E. sibiricus resist drought. This study used transcriptomics to investigate how E. sibiricus responds to drought stress, and may provide genetic resources and references for research into the molecular mechanisms of drought resistance in native perennial grasses and for breeding drought-tolerant varieties.

Analysis on the salt tolerance of Nitraria sibirica Pall. based on Pacbio full-length transcriptome sequencing.

Nitraria sibirica Pall. regulates its tolerance to salt stress mainly by adjusting ion balance, modifying cell wall structure, and activating signal transduction pathways. N. sibirica, as a typical halophyte, can not only effectively restore saline-alkali land, but also has high economic value. However, studies on its salt tolerance at combining molecular and physiological levels were limited. In this study, the salt tolerance of N. sibirica was analyzed based on Pacbio full-length transcriptome sequencing, and the salt tolerance in the physiological level was verified by key genes. The results showed that 89,017 full-length transcripts were obtained, of which 84,632 sequences were annotated. A total of 86,482 coding sequences (CDS) were predicted and 6561 differentially expressed genes (DEGs) were identified. DEGs were significantly enriched in "sodium ion homeostasis", "response to osmotic stress", "reactive oxygen species metabolic process", "defense response by cell wall thickening", "signal transduction", etc. The expression levels for most of these DEGs increased under salt stress. A total of 69 key genes were screened based on weighted gene co-expression network analysis (WGCNA), of which 33 were first reported on salt tolerance. Moreover, NsRabE1c gene with the highest expression level was selected to verify its salt tolerance. Over-expression of NsRabE1c gene enhanced the germination potential and root length of transgenic Arabidopsis thaliana plants without salt treatment as compared to those of Col-0 and AtRabE1c mutant. The expression levels of NsRabE1c decreased in the growth stagnation phase, while significantly increased in the growth recovery phase under salt stress. We predicted that NsRabE1c gene help N. sibirica resist salt stress through the regulation of plant growth. The results of this study deepen the understanding of salinity resistance in N. sibirica.

Identifying and quantifying isoforms from accurate full-length transcriptome sequencing reads with Mandalorion

In this manuscript, we introduce and benchmark Mandalorion v4.1 for the identification and quantification of full-length transcriptome sequencing reads. It further improves upon the already strong performance of Mandalorion v3.6 used in the LRGASP consortium challenge. By processing real and simulated data, we show three main features of Mandalorion: first, Mandalorion-based isoform identification has very high precision and maintains high recall even in the absence of any genome annotation. Second, isoform read counts as quantified by Mandalorion show a high correlation with simulated read counts. Third, isoforms identified by Mandalorion closely reflect the full-length transcriptome sequencing data sets they are based on.

Transcriptomic analysis of Andrias davidianus meat and experimental validation for exploring its bioactive components as functional foods

Andrias davidianus (Chinese giant salamander, CGS) is the largest and oldest extant amphibian species in the world and is a source of prospective functional food in China. However, the progress of functional peptides mining was slow due to lack of reference genome and protein sequence data. In this study, we illustrated full-length transcriptome sequencing to interpret the proteome of CGS meat and obtain 10703 coding DNA sequences. By functional annotation and amino acid composition analysis, we have discovered various genes related to signal transduction, and 16 genes related to longevity. We have also found vast variety of functional peptides through protein coding sequence (CDS) analysis by comparing the data obtained with the functional peptide database. Val-Pro-Ile predicted by the CDS analysis was released from the CGS meat through enzymatic hydrolysis, suggesting that our approach is reliable. This study suggested that transcriptomic analysis can be used as a reference to guide polypeptide mining in CGS meat, thereby providing a powerful mining strategy for the bioresources with unknown genomic and proteomic sequences.

SMRT Sequencing Technology Was Used to Construct the Batocera horsfieldi (Hope) Transcriptome and Reveal Its Features.

Batocera horsfieldi (Hope) (Coleoptera: Cerambycidae) is an important forest pest in China that mainly infests timber and economic forests. This pest primarily causes plant tissue to necrotize, rot, and eventually die by feeding on the woody parts of tree trunks. To gain a deeper understanding of the genetic mechanism of B. horsfieldi, this study employed single-molecule real-time sequencing (SMRT) and Illumina RNA-seq technologies to conduct full-length transcriptome sequencing of the insect. Total RNA extracted from male and female adults was mixed and subjected to SMRT sequencing, generating a complete transcriptome. Transcriptome analysis, prediction of long non-coding RNA (lncRNA), coding sequences (CDs), analysis of simple sequence repeats (SSR), prediction of transcription factors, and functional annotation of transcripts were performed in this study. The collective 20,356,793 subreads (38.26 G, clean reads) were generated, including 432,091 circular consensus sequences and 395,851 full-length non-chimera reads. The full-length non-chimera reads (FLNC) were clustered and redundancies were removed, resulting in 39,912 consensus reads. SSR and ANGEL software v3.0 were used for predicting SSR and CDs. In addition, four tools were used for annotating 6058 lncRNAs, identifying 636 transcription factors. Furthermore, a total of 84,650 transcripts were functionally annotated in seven different databases. This is the first time that the full-length transcriptome of B. horsfieldi has been obtained using SMRT sequencing. This provides an important foundation for investigating the gene regulation underlying the interaction between B. horsfieldi and its host plants through gene editing in the future and provides a scientific basis for the prevention and control of B. horsfieldi.

Matrix metalloproteinases are key targets of acupuncture in the treatment of ulcerative colitis.

The aim of this study was to elucidate the key targets of acupuncture in the colon of ulcerative colitis (UC) mice model using full-length transcriptome sequencing. 2.5% dextran sodium sulfate (DSS)-induced colitis mice were treated with or without acupuncture. Intestinal pathology was observed, and full transcriptome sequencing and bioinformatic analysis were performed. The results demonstrated that acupuncture treatment reduced the UC symptoms, disease activity index score, and histological colitis score and increased body weight, colon length, and the number of intestinal goblet cells. In addition, acupuncture can also decrease the expression of necrotic biomarker phosphorylates mixed lineage kinase domain-like pseudo kinase (p-MLKL). Full-length transcriptome analysis indicated that acupuncture reversed the expression of 987 of the 1918 upregulated differentially expressed genes (DEGs), and 632 of the 1351 downregulated DEGs induced by DSS. DEGs regulated by acupuncture were mainly involved in inflammatory responses and intestinal barrier pathways. The protein-protein interaction network analysis revealed that matrix metalloproteinases (MMPs) are important genes regulated by acupuncture. Gene set enrichment analysis revealed that extracellular matrix (ECM)-receptor interaction was an important target of acupuncture. In addition, alternative splicing analysis suggested that acupuncture improved signaling pathways related to intestinal permeability, the biological processes of xenobiotics, sulfur compounds, and that monocarboxylic acids are closely associated with MMPs. Overall, our transcriptome analysis results indicate that acupuncture improves intestinal barrier function in UC through negative regulation of MMPs expression.

Full-Length Transcriptome Sequencing Reveals the Molecular Mechanism of Metasequoia glyptostroboides Seed Responding to Aging.

Metasequoia glyptostroboides, Hu and W. C. Cheng, as the only surviving relict species of the Taxodiaceae Metasequoia genus, is a critically endangered and protected species in China. There is a risk of extinction due to the low vigor of M. glyptostroboides seeds, and the physiological mechanism of seed aging in M. glyptostroboides is not yet clear. In order to investigate the physiological and molecular mechanisms underlying the aging process of M. glyptostroboides seeds, we analyzed the antioxidant system and transcriptome at 0, 2, 4, 6, and 8 days after artificial accelerated aging treatment at 40 °C and 100% relative humidity. It was found that the germination percentage of fresh dried M. glyptostroboides seeds was 54 ± 5.29%, and significantly declined to 9.33 ± 1.88% after 6 days of aging, and then gradually decreased until the seed died on day 8. Superoxide dismutase (SOD) activity, ascorbic acid (AsA), glutathione (GSH) content and superoxide anion (O2·-) content and production rate significantly decreased, while malondialdehyde (MDA) and hydrogen peroxide (H2O2) content and glutathione peroxidase (GPX) and catalase (CAT) activity gradually increased during the aging process. A total of 42,189 unigenes were identified in the whole transcriptome, and 40,446 (95.86%) unigenes were annotated in at least one protein database. A total of 15,376 differentially expressed genes (DEGs) were obtained; KEGG enrichment analysis results revealed that seed aging may be mainly involved in the protein-processing pathways in endoplasmic reticulum, oxidative phosphorylation, and ascorbate and aldarate metabolism. Weighted gene co-expression network analysis (WGCNA) revealed that the dark magenta, orange, and medium purple modules were highly correlated with physiological indicators such as SOD, CAT, and GSH and further identified 40 hub genes such as Rboh, ACO, HSF, and CML as playing important roles in the antioxidant network of M. glyptostroboides seeds. These findings provide a broader perspective for studying the regulatory mechanism of seed aging and a large number of potential target genes for the breeding of other endangered gymnosperms.

Full-Length Transcriptome Analysis of Soybean Cyst Nematode (Heterodera glycines) Reveals an Association of Behaviors in Response to Attractive pH and Salt Solutions with Activation of Transmembrane Receptors, Ion Channels, and Ca2+ Transporters.

Soybean cyst nematode (Heterodera glycines Ichinohe), a devastating pathogen in soybean, was chosen as a model system to investigate nematode behavior and gene expression changes in response to acidic and basic pH and salt signals (pH 4.5, 5.25, 8.6, and 10 and NaCl) through full-length transcriptome sequencing of 18 samples. An average of 4.36 Gbp of clean reads per sample were generated, and 3972 novel genes and 29,529 novel transcripts were identified. Sequence structural variation during or after transcription may be associated with the nematode's behavioral response. The functional analysis of 1817/4962 differentially expressed genes/transcripts showed that signal transduction pathways, including transmembrane receptors, ion channels, and Ca2+ transporters, were activated, but pathways involved in nematode development (e.g., ribosome) and energy production (e.g., oxidative phosphorylation) were inhibited. A corresponding model was established. Our findings suggest that these receptors and ion channels might be potential targets for nematicides or drug discovery.