Search

Under specific cultivation conditions, the sepal color of Hydrangea macrophylla (H. macrophylla) changes from red to blue due to the complexation of aluminum ions (Al3+), delphinidin 3-glucoside, and copigments. However, this phenomenon cannot occur in all cultivars despite the presence of sufficient Al3+ and copigments. To explore the mechanism of sepal bluing in H. macrophylla, there is an urgent need to study the molecular regulation of the anthocyanin biosynthesis pathway. However, the key structural genes, other than CHS, regulating anthocyanin biosynthesis in the sepals of H. macrophylla have not been identified. In this study, based on full-length transcriptome data from H.macrophylla ‘Bailmer’, the key structural genes regulating anthocyanin biosynthesis in the sepals of H. macrophylla were isolated and investigated. Ultimately, seven key structural genes, HmCHS1, HmCHI, HmF3H1, HmF3′H1, HmF3′5′H, HmDFR2, and HmANS3, were demonstrated to show high expression levels in colored sepals. The expression levels of these seven genes increased gradually with the development of sepals and were highest in the full-bloom stage. The trend of gene expression was consistent with the trend of anthocyanin contents. It was concluded that the seven selected genes were involved in anthocyanin biosynthesis in the sepals of H. macrophylla. The full-length sequence data have been deposited into the NCBI Sequence Read Archive (SRA) with accession number PRJNA849710. This study lays a good foundation for the further elucidation of the molecular mechanism of sepal coloration in H. macrophylla.

The Great Himalayan Leaf-nosed bat (Hipposideros armiger) is one of the most representative species of all echolocating bats and is an ideal model for studying the echolocation system of bats. An incomplete reference genome and limited availability of full-length cDNAs have hindered the identification of alternatively spliced transcripts, which slowed down related basic studies on bats’ echolocation and evolution. In this study, we analyzed five organs from H. armiger for the first time using PacBio single-molecule real-time sequencing (SMRT). There were 120 GB of subreads generated, including 1,472,058 full-length non-chimeric (FLNC) sequences. A total of 34,611 alternative splicing (AS) events and 66,010 Alternative Polyadenylation (APA) sites were detected by transcriptome structural analysis. Moreover, a total of 110,611 isoforms were identified, consisting of 52% new isoforms of known genes and 5% of novel gene loci, as well as 2112 novel genes that have not been annotated before in the current reference genome of H. armiger. Furthermore, several key novel genes, including Pol, RAS, NFKB1, and CAMK4, were identified as being associated with nervous, signal transduction, and immune system processes, which may be involved in regulating the auditory nervous perception and immune system that helps bats to regulate in echolocation. In conclusion, the full-length transcriptome results optimized and replenished existing H. armiger genome annotation in multiple ways and offer advantages for newly discovered or previously unrecognized protein-coding genes and isoforms, which can be used as a reference resource.

Background Light switches on the photomorphogenic development of young plant seedlings, allowing young seedlings to acquire photosynthetic capacities and gain survival fitness. Light regulates gene expression at all levels of the central dogma, including alternative splicing (AS) during the photomorphogenic development. However, accurate determination of full-length (FL) splicing variants has been greatly hampered by short-read RNA sequencing technologies. Result In this study, we adopt PacBio isoform sequencing (Iso-seq) to overcome the limitation of the short-read RNA-seq technologies. Normalized cDNA libraries used for Iso-seq allows for comprehensive and effective identification of FL AS variants. Our analyses reveal more than 30,000 splicing variant models from approximately 16,500 gene loci and additionally identify approximately 700 previously unannotated genes. Among the variants, approximately 12,000 represent new gene models. Intron retention (IR) is the most frequently observed form of variants, and many IR-containing AS variants show evidence of engagement in translation. Our study reveals the formation of heterodimers of transcription factors composed of annotated and IR-containing AS variants. Moreover, transgenic plants overexpressing the IR forms of two B-BOX DOMAIN PROTEINs exhibits light-hypersensitive phenotypes, suggesting their regulatory roles in modulating optimal light responses. Conclusions This study provides an accurate and comprehensive portrait of full-length transcript isoforms and experimentally confirms the presence of de novo synthesized AS variants that impose regulatory functions in photomorphogenic development in Arabidopsis.

Alsophila spinulosa, a relict tree fern, is a valuable plant for investigating environmental adaptations. Its genetic resources, however, are scarce. We used the PacBio and Illumina platforms to sequence the polyadenylated RNA of A. spinulosa root, rachis, and pinna, yielding 125,758, 89,107, and 89,332 unigenes, respectively. Combining the unigenes from three organs yielded a non-redundant reference transcriptome with 278,357 unigenes and N50 of 4141 bp, which were further reconstructed into 38,470 UniTransModels. According to functional annotation, pentatricopeptide repeat genes and retrotransposon-encoded polyprotein genes are the most abundant unigenes. Clean reads mapping to the full-length transcriptome is used to assess the expression of unigenes. The stress-induced ASR genes are highly expressed in all three organs, which is validated by qRT-PCR. The organ-specific upregulated genes are enriched for pathways involved in stress response, secondary metabolites, and photosynthesis. Genes for five types of photoreceptors, CRY signaling pathway, ABA biosynthesis and transduction pathway, and stomatal movement-related ion channel/transporter are profiled using the high-quality unigenes. The gene expression pattern coincides with the previously identified stomatal characteristics of fern. This study is the first multi-organ full-length transcriptome report of a tree fern species, the abundant genetic resources and comprehensive analysis of A. spinulosa, which provides the groundwork for future tree fern research.

Goldfish Carassius auratus is an ideal model for exploring fish morphology evolution. Although genes underlying several ornamental traits have been identified, little is known about the effects of artificial selection on embryo gene expression. In the present study, hybrid transcriptome sequencing was conducted to reveal gene expression profiles of Celestial-Eye (CE) and Ryukin (RK) goldfish embryos. Full-length transcriptome sequencing on the PacBio platform identified 54,218 and 54,106 transcript isoforms in CE and RK goldfish, respectively. Of particular note was that thousands of alternative splicing (AS) and alternative polyadenylation (APA) events were identified in both goldfish breeds, and most of them were inter-breed specific. RT-PCR and Sanger sequencing showed that most of the predicted AS and APA were correct. Moreover, abundant long non-coding RNA and fusion genes were detected, and again most of them were inter-breed specific. Through RNA-seq, we detected thousands of differentially expressed genes (DEGs) in each embryonic stage between the two goldfish breeds. KEGG enrichment analysis on DEGs showed extensive differences between CE and RK goldfish in gene expression. Taken together, our results demonstrated that artificial selection has led to far-reaching influences on goldfish gene expression, which probably laid the genetic basis for hundreds of goldfish variations.

Tree peony (Paoenia ostii) is a famous Chinese traditional flower well-known in many countries of the world. However, the short and concentrated flowering period of tree peony greatly affects the ornamental and economic value of the flowers. Exogenous brassinosteroid (BR) treatment can delay the flowering period of ostii T. Hong et J. X. Zhang var. lishizhenenii B. A. Shen for 3 days, but the underlying regulatory mechanism remains elusive. Here, full-length transcriptome and transcriptome sequencing were used to mine key genes related to BR-induced delayed flowering in tree peony. The transcriptome sequencing of the petals yielded 21.27 G clean data and 62,229 isoforms. Among them, 58,218 isoforms were annotated in NR, NT, SwissProt, KEGG, KOG, InterPro and GO databases. GO and KEGG analyses showed that 2460 DEGs were related to delayed flowering in response to BR. Additionally, a total of seven genes affecting flowering were annotated from 11 isoforms, which responded to BR through three pathways to delay the flowering of P. ostii var. lishizhenii. BR treatment increased the expression of BRASSINOSTEROID-SIGNALING KINASE3 (BSK3), potentially by promoting BRI1 ASSOCIATED KINASE RECEPTOR 1 (BAK1). Moreover, BR treatment suppressed the expression of SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE 1 (SPL1), possibly to inhibit the expressions of downstream WRKY genes and APETALA 2 (AP2). Meanwhile, BR treatment promoted the expression of DELLA, which might inhibit the expressions of CONSTANS (CO) and SPL. These results suggest a theoretical basis for further analyses of the molecular mechanism of flowering regulation in tree peony.

Adipose tissues have a central role in organisms, and adipose content is a crucial economic trait of poultry. Pekin duck is an ideal model to study the mechanism of abdominal and subcutaneous adipose deposition for its high ability of adipose synthesis and deposition. Alternative splicing contributes to functional diversity in abdominal and subcutaneous adipose. However, there has been no systematic analysis of the dynamics of differential alternative splicing of abdominal and subcutaneous adipose in Pekin duck. In our study, the Pacific Biosciences (PacBio) Iso-Seq technology was applied to explore the transcriptional complexity of abdominal and subcutaneous adipose in Pekin ducks. In total, 143,931 and 111,337 full-length non-chimeric transcriptome sequences of abdominal and subcutaneous adipocytes were obtained from 41.78 GB raw data, respectively. These data led us to identify 19,212 long non-coding RNAs (lncRNAs) and 74,571 alternative splicing events. In addition, combined with the next-generation sequencing technology, we correlated the structure and function annotation with the differential expression profiles of abdominal and subcutaneous adipose transcripts. This study identified lots of novel alternative splicing events and major transcripts of transcription factors related to adipose synthesis. STAT3 was reported as a vital gene for adipogenesis, and we found that its major transcript is STAT3-1, which may play a considerable role in the process of adipose synthesis in Pekin duck. This study greatly increases our understanding of the gene models, genome annotations, genome structures, and the complexity and diversity of abdominal and subcutaneous adipose in Pekin duck. These data provide insights into the regulation of alternative splicing events, which form an essential part of transcript diversity during adipogenesis in poultry. The results of this study provide an invaluable resource for studying alternative splicing and tissue-specific expression.

Background and Objectives: Castor (Ricinus communis L.) is an important non-edible oilseed crop. Lm-type female strains and normal amphiprotic strains are important castor cultivars, and are mainly different in their inflorescence structures and leaf shapes. To better understand the mechanisms underlying these differences at the molecular level, we performed a comparative transcriptional analysis.Materials and Methods: Full-length transcriptome sequencing and short-read RNA sequencing were employed.Results: A total of 76,068 and 44,223 non-redundant transcripts were obtained from high-quality transcripts of Lm-type female strains and normal amphiprotic strains, respectively. In Lm-type female strains and normal amphiprotic strains, 51,613 and 20,152 alternative splicing events were found, respectively. There were 13,239 transcription factors identified from the full-length transcriptomes. Comparative analysis showed a great variety of gene expression of common and unique transcription factors between the two cultivars. Meanwhile, a functional analysis of the isoforms was conducted. The full-length sequences were used as a reference genome, and a short-read RNA sequencing analysis was performed to conduct differential gene analysis. Furthermore, the function of DEGs were performed to annotation analysis.Conclusion: The results revealed considerable differences and expression diversity between the two cultivars, well beyond what was reported in previous studies and likely reflecting the differences in architecture between these two cultivars.

Seed-borne endophyte Epichloë gansuensis enhance NaCl tolerance in Achnatherum inebrians and increase its biomass. However, the molecular mechanism by which E. gansuensis increases the tolerance of host grasses to NaCl stress is unclear. Hence, we firstly explored the full-length transcriptome information of A. inebrians by PacBio RS II. In this work, we obtained 738,588 full-length non-chimeric reads, 36,105 transcript sequences and 27,202 complete CDSs from A. inebrians. We identified 3558 transcription factors (TFs), 15,945 simple sequence repeats and 963 long non-coding RNAs of A. inebrians. The present results show that 2464 and 1817 genes were differentially expressed by E. gansuensis in the leaves of E+ and E− plants at 0 mM and 200 mM NaCl concentrations, respectively. In addition, NaCl stress significantly regulated 4919 DEGs and 502 DEGs in the leaves of E+ and E− plants, respectively. Transcripts associated with photosynthesis, plant hormone signal transduction, amino acids metabolism, flavonoid biosynthetic process and WRKY TFs were differentially expressed by E. gansuensis; importantly, E. gansuensis up-regulated biology processes (brassinosteroid biosynthesis, oxidation–reduction, cellular calcium ion homeostasis, carotene biosynthesis, positive regulation of proteasomal ubiquitin-dependent protein catabolism and proanthocyanidin biosynthesis) of host grass under NaCl stress, which indicated an increase in the ability of host grasses’ adaptation to NaCl stress. In conclusion, our study demonstrates the molecular mechanism for E. gansuensis to increase the tolerance to salt stress in the host, which provides a theoretical basis for the molecular breed to create salt-tolerant forage with endophytes.

Tanshinones are the bioactive constituents of Danshen (Salvia miltiorrhiza Bunge), which is used in Traditional Chinese Medicine to treat cardiovascular and other diseases, and they synthesize and accumulate in the root periderm of S. miltiorrhiza. However, there is no relevant report on the initial stage of tanshinone synthesis, as well as the root structure and gene expression characteristics. The present study aims to provide new insights into how these bioactive principles begin to synthesize by characterizing possible differences in their biosynthesis and accumulation during early root development from both spatial and temporal aspects. The morphological characteristics and the content of tanshinones in roots of S. miltiorrhiza were investigated in detail by monitoring the seedlings within 65 days after germination (DAGs). The ONT transcriptome sequencing was applied to investigate gene expression patterns. The periderm of the S. miltiorrhiza storage taproot initially synthesized tanshinone on about 30 DAGs. Three critical stages of tanshinone synthesis were preliminarily determined: preparation, the initial synthesis, and the continuous rapid synthesis. The difference of taproots in the first two stages was the smallest, and the differentially expressed genes (DEGs) were mainly enriched in terpene synthesis. Most genes involved in tanshinone synthesis were up regulated during the gradual formation of the red taproot. Plant hormone signal transduction and ABC transport pathways were widely involved in S. miltiorrhiza taproot development. Five candidate genes that may participate in or regulate tanshinone synthesis were screened according to the co-expression pattern. Moreover, photosynthetic ferredoxin (FD), cytochrome P450 reductase (CPR), and CCAAT binding transcription factor (CBF) were predicted to interact with the known downstream essential enzyme genes directly. The above results provide a necessary basis for analyzing the initial synthesis and regulation mechanism of Tanshinones.

Air exposure is a common stress for Corbicula fluminea, an economically important freshwater shellfish consumed in China, during aquaculture and transportation. However, little is known about its molecular responses to air exposure. Therefore, this study used a combination of PacBio full-length and Illumina transcriptomes to investigate its molecular responses to air exposure. A total of 36,772 transcripts were obtained using PacBio sequencing. Structural analysis identified 32,069 coding sequences, 1906 transcription factors, 8873 simple sequence repeats, and 17,815 long non-coding RNAs. Subcellular localization analysis showed that most transcripts were located in the cytoplasm and nucleus. After 96-h of air exposure, 210 differentially expressed genes (DEGs) in the gill were obtained via Illumina sequencing. Among these DEGs, most of the genes related to glycolysis, tricarboxylic acid cycle, lipid metabolism, and amino acid metabolism were upregulated. Additionally, many DEGs associated with immunity, cytoskeleton reorganization, autophagy, and ferroptosis were identified. These findings indicated that metabolic strategy change, immune response, cytoskeleton reconstruction, autophagy, and ferroptosis might be the important mechanisms that C. fluminea use to cope with air exposure. This study will enrich the gene resources of C. fluminea and provide valuable data for studying the molecular mechanisms coping with air exposure in C. fluminea and other freshwater mollusks.

Background Kiwifruit (Actinidia Lindl.) is considered an important fruit species worldwide. Due to its temperate origin, this species is highly vulnerable to freezing injury while under low-temperature stress. To obtain further knowledge of the mechanism underlying freezing tolerance, we carried out a hybrid transcriptome analysis of two A. arguta (Actinidi arguta) genotypes, KL and RB, whose freezing tolerance is high and low, respectively. Both genotypes were subjected to − 25 °C for 0 h, 1 h, and 4 h. Results SMRT (single-molecule real-time) RNA-seq data were assembled using the de novo method, producing 24,306 unigenes with an N50 value of 1834 bp. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of DEGs showed that they were involved in the ‘starch and sucrose metabolism’, the ‘mitogen-activated protein kinase (MAPK) signaling pathway’, the ‘phosphatidylinositol signaling system’, the ‘inositol phosphate metabolism’, and the ‘plant hormone signal transduction’. In particular, for ‘starch and sucrose metabolism’, we identified 3 key genes involved in cellulose degradation, trehalose synthesis, and starch degradation processes. Moreover, the activities of beta-GC (beta-glucosidase), TPS (trehalose-6-phosphate synthase), and BAM (beta-amylase), encoded by the abovementioned 3 key genes, were enhanced by cold stress. Three transcription factors (TFs) belonging to the AP2/ERF, bHLH (basic helix-loop-helix), and MYB families were involved in the low-temperature response. Furthermore, weighted gene coexpression network analysis (WGCNA) indicated that beta-GC, TPS5, and BAM3.1 were the key genes involved in the cold response and were highly coexpressed together with the CBF3, MYC2, and MYB44 genes. Conclusions Cold stress led various changes in kiwifruit, the ‘phosphatidylinositol signaling system’, ‘inositol phosphate metabolism’, ‘MAPK signaling pathway’, ‘plant hormone signal transduction’, and ‘starch and sucrose metabolism’ processes were significantly affected by low temperature. Moreover, starch and sucrose metabolism may be the key pathway for tolerant kiwifruit to resist low temperature damages. These results increase our understanding of the complex mechanisms involved in the freezing tolerance of kiwifruit under cold stress and reveal a series of candidate genes for use in breeding new cultivars with enhanced freezing tolerance.

Coelomactra antiquata is an important aquatic economic shellfish with high medicinal value. However, because C. antiquata has no reference genome, a lot of molecular biology research cannot be carried out, so the analysis of its transcripts is an important step to study the regulatory genes of various substances in C. antiquata. In the present study, we conducted the first full-length transcriptome analysis of C. antiquata by using PacBio single-molecule real-time (SMRT) sequencing technology. The results identified a total of 39,209 unigenes with an average length of 2,732 bp, 23,338 CDSs, 251 AS events, 9,881 lncRNAs, 20,106 SSRs, and 2,316 TFs. Subsequently, 59.22% (23,220) of the unigenes were successfully annotated, of which 23,164, 18,711, 15,840, 13,534, and 13,474 unigenes could be annotated using NR, Swiss-prot, KOG, GO, and KEGG databases, respectively. This study lays the foundation for the follow-up research of molecular biology and provides a reference for studying the more medicinal value of C. antiquata.

Background The pathological process of myocardial ischemia (MI) is very complicated. Acupuncture at PC6 has been proved to be effective against MI injury, but the mechanism remains unclear. This study investigated the mechanism that underlies the effect of acupuncture on MI through full-length transcriptome. Methods Adult male C57/BL6 mice were randomly divided into control, MI, and PC6 groups. Mice in MI and PC6 group generated MI model by ligating the left anterior descending (LAD) coronary artery. The samples were collected 5 days after acupuncture treatment. Results The results showed that treatment by acupuncture improved cardiac function, decreased myocardial infraction area, and reduced the levels of cTnT and cTnI. Based on full-length transcriptome sequencing, 5083 differential expression genes (DEGs) and 324 DEGs were identified in the MI group and PC6 group, respectively. These genes regulated by acupuncture were mainly enriched in the inflammatory response pathway. Alternative splicing (AS) is a post-transcriptional action that contributes to the diversity of protein. In all samples, 8237 AS events associated with 1994 genes were found. Some differential AS-involved genes were enriched in the pathway related to heart disease. We also identified 602 new genes, 4 of which may the novel targets of acupuncture in MI. Conclusions Our findings suggest that the effect of acupuncture on MI may be based on the multi-level regulation of the transcriptome.

Several main families of Ranunculales are rich in alkaloids and other medicinal compounds; many species of these families are used in traditional and folk medicine. Dichocarpum is a representative medicinal genus of Ranunculaceae, but the genetic basis of its metabolic phenotype has not been investigated, which hinders its sustainable conservation and utilization. We use the third-generation high-throughput sequencing and metabolomic techniques to decipher the full-length transcriptomes and metabolomes of five Dichocarpum species endemic in China, and 71,598 non-redundant full-length transcripts were obtained, many of which are involved in defense, stress response and immunity, especially those participating in the biosynthesis of specialized metabolites such as benzylisoquinoline alkaloids (BIAs). Twenty-seven orthologs extracted from trancriptome datasets were concatenated to reconstruct the phylogenetic tree, which was verified by the clustering analysis based on the metabolomic profile and agreed with the Pearson correlation between gene expression patterns of Dichocarpum species. The phylogenomic analysis of phytometabolite biosynthesis genes, e.g., (S)-norcoclaurine synthase, methyltransferases, cytochrome p450 monooxygenases, berberine bridge enzyme and (S)-tetrahydroprotoberberine oxidase, revealed the evolutionary trajectories leading to the chemodiversity, especially that of protoberberine type, aporphine type and bis-BIA abundant in Dichocarpum and related genera. The biosynthesis pathways of these BIAs are proposed based on full-length transcriptomes and metabolomes of Dichocarpum. Within Ranunculales, the gene duplications are common, and a unique whole genome duplication is possible in Dichocarpum. The extensive correlations between metabolite content and gene expression support the co-evolution of various genes essential for the production of different specialized metabolites. Our study provides insights into the transcriptomic and metabolomic landscapes of Dichocarpum, which will assist further studies on genomics and application of Ranunculales plants.

As the core of heterosis utilization, cytoplasmic male sterility (CMS) has been widely used in hybrid seed production. Previous studies have shown that CMS is always closely related to the altered programming of mitochondrial genes. To explore candidate CMS genes in cotton (Gossypium hirsutum), sequencing and de novo assembly were performed on the mitochondrial genome of the G. hirsutum CMS line SI3A, with G. harknessii CMS-D2 cytoplasm, and the corresponding G. hirsutum restorer line 0-613-2R. Remarkable variations in genome structure and gene transcripts were detected. The mitochondrial genome of SI3A has three circle molecules, including one main circle and two sub-circles, while 0-613-2R only has one. RNA-seq and RT-qPCR analysis proved that orf606a and orf109a, which have a chimeric structure and transmembrane domain, were highly expressed in abortive anthers of SI3A. In addition, comparative analysis of RNA-seq and full-length transcripts revealed the complex I gene nad4 to be expressed at a lower level in SI3A than in its restorer and that it featured an intron retention splicing pattern. These two novel chimeric ORFs and nad4 are potential candidates that confer CMS character in SI3A. This study provides new insight into the molecular basis of the nuclear–cytoplasmic interaction mechanism, and that putative CMS genes might be important sources for future precise design cross-breeding of cotton.

Background Spodoptera frugiperda (J. E. Smith), commonly known as fall armyworm (FAW), is one of the most destructive agricultural pests in the world and has posed a great threat to crops. The improper use of insecticides has led to rapid development of resistance. However, the genetic data available for uncovering the insecticide resistance mechanisms are scarce. Methods In this study, we used PacBio single-molecule real-time (SMRT) sequencing aimed at revealing the full-length transcriptome profiling of the FAW larval brain to obtain detoxification genes. Results A total of 18,642 high-quality transcripts were obtained with an average length of 2,371 bp, and 11,230 of which were successfully annotated in six public databases. Among these, 5,692 alternative splicing events were identified.

The Tibetan Schizothoracinae fish Gymnocypris przewalskii has the ability to adapt to the extreme plateau environment, making it an ideal biological material for evolutionary biology research. However, the lack of well-annotated reference genomes has limited the study of the molecular genetics of G. przewalskii. To characterize its transcriptome features, we first used long-read sequencing technology in combination with RNA-seq for transcriptomic analysis. A total of 159,053 full-length (FL) transcripts were captured by Iso-Seq, having a mean length of 3,445 bp with N50 value of 4,348. Of all FL transcripts, 145,169 were well-annotated in the public database and 134,537 contained complete open reading frames. There were 4,149 pairs of alternative splicing events, of which three randomly selected were defined by RT–PCR and sequencing, and 13,293 long non-coding RNAs detected, based on all-vs.-all BLAST. A total of 118,185 perfect simple sequence repeats were identified from FL transcripts. The FL transcriptome might provide basis for further research of G. przewalskii.

Kiwifruit an important horticultural crop that is widely cultivated and is known as the king of fruits. Recently, a new seedless kiwifruit cultivar, ‘Chengxiang’ (Actinidia arguta), was discovered by field transplantation. It exhibited distinguishable characteristics such as parthenocarpy, and a unique flavor and appearance when compared to other cultivated type. Flavonoids are known to play an important role in fertility and parthenocarpy in plants. However, the genes responsible for flavonoid biosynthesis in seedless kiwifruit remain largely unknown. Especially, chalcone synthase (CHS), as a key enzyme catalyzing the first committed step in the flavonoid pathway, remains a mystery. In this study, we combined a full-length transcriptome survey by PacBio single-molecule real-time (SMRT) sequencing, CHS gene family analysis, and analysis of the gene expression involved in flavonoid pathways to further enhance the understanding of parthenocarpy. Based on SMRT, we obtained 80,615 high-quality full-length consensus transcripts. In total, 52,406 (90.79%) transcripts were functionally annotated, and more than 80% of the transcripts were longer than 1Kb. Among them, 39,117 (74.64%) transcripts were assigned to GO terms, the majority of which were associated with the cell (19,089, 48.80%) and metabolic process (19,859, 50.77%). Furthermore, 25,289 (48.26%) transcripts were mapped into 129 KEGG pathways. We identified the majority of putative genes as being involved in the flavonoid biosynthesis pathway, including 14 key enzyme gene families, such as CHS, chalcone isomerase (CHI), flavonol synthase (FLS), and so on. Moreover, we also identified 13 CHS genes and characterized the CHS gene family in seedless kiwifruit. We further evaluated the expression pattern of 10 flavonoid-related key enzyme genes in flowers using quantitative real-time PCR. This is the first time that the full-length transcriptome have been studied in seedless kiwifruit, and the findings enhance our understanding the molecular mechanisms of parthenocarpy.

Actinidia latifolia is one of the very few kiwifruit genotypes with extremely high ascorbic acid (AsA) content. However, a transcriptome atlas of this species is lacking. The accumulation of AsA during fruit development and ripening and the associated molecular mechanisms are still poorly understood. Herein, dynamic changes in AsA content at six different stages of A. latifolia fruit development and ripening were determined. AsA content of A. latifolia fruit reached 1108.76 ± 35.26 mg 100 g−1 FW at full maturity. A high-quality, full-length (FL) transcriptome of A. latifolia was successfully constructed for the first time using third-generation sequencing technology. The transcriptome comprises 326,926 FL non-chimeric reads, 15,505 coding sequences, 2882 transcription factors, 18,797 simple sequence repeats, 3328 long noncoding RNAs, and 231 alternative splicing events. The genes involved in AsA biosynthesis and recycling pathways were identified and compared with those in different kiwifruit genotypes. The correlation between the AsA content and expression levels of key genes in AsA biosynthesis and recycling pathways was revealed. LncRNAs that participate in AsA-related gene expression regulation were also identified. Gene expression patterns in AsA biosynthesis and metabolism exhibited a trend similar to that of AsA accumulation. Overall, this study paves the way for genetic engineering to develop kiwifruits with super-high AsA content.

Sirex noctilio is a major international quarantine pest that recently emerged in northeast China to specifically invade conifers. During female oviposition, venom is injected into the host together with its symbiotic fungus to alter the normal Pinus physiology and weaken or even kill the tree. In China, the Mongolian pine (Pinus sylvestris var. mongolica), an important wind-proof and sand-fixing species, is the unique host of S. noctilio. To explore the interplay between S. noctilio venom and Mongolian pine, we performed a transcriptome comparative analysis of a 10-year-old Mongolian pine after wounding and inoculation with S. noctilio venom. The analysis was performed at 12 h, 24 h and 72 h. PacBio ISO-seq was used and integrated with RNA-seq to construct an accurate full-length transcriptomic database. We obtained 52,963 high-precision unigenes, consisting of 48,654 (91.86%) unigenes that were BLASTed to known sequences in the public database and 4309 unigenes without any annotation information, which were presumed to be new genes. The number of differentially expressed genes (DEGs) increased with the treatment time, and the DEGs were most abundant at 72 h. A total of 706 inoculation-specific DEGs (475 upregulated and 231 downregulated) and 387 wounding-specific DEGs (183 upregulated and 204 downregulated) were identified compared with the control. Under venom stress, we identified 6 DEGs associated with reactive oxygen species (ROS) and 20 resistance genes in Mongolian pine. Overall, 52 transcription factors (TFs) were found under venom stress, 45 of which belonged to the AP2/ERF TF family and were upregulated. A total of 13 genes related to the photosystem, 3 genes related photo-regulation, and 9 TFs were identified under wounding stress. In conclusion, several novel putative genes were found in Mongolian pine by PacBio ISO seq. Meanwhile, we also identified various genes that were resistant to S. noctilio venom, such as GAPDH, GPX, CAT, FL2, CERK1, and HSP83A, etc.

Inbred animals were historically chosen for genome analysis to circumvent assembly issues caused by haplotype variation but this resulted in a composite of the two genomes. Here we report a haplotype-aware scaffolding and polishing pipeline which was used to create haplotype-resolved, chromosome-level genome assemblies of Angus (taurine) and Brahman (indicine) cattle subspecies from contigs generated by the trio binning method. These assemblies reveal structural and copy number variants that differentiate the subspecies and that variant detection is sensitive to the specific reference genome chosen. Six genes with immune related functions have additional copies in the indicine compared with taurine lineage and an indicus-specific extra copy of fatty acid desaturase is under positive selection. The haplotyped genomes also enable transcripts to be phased to detect allele-specific expression. This work exemplifies the value of haplotype-resolved genomes to better explore evolutionary and functional variations., Taurine and indicine cattle have different desirable traits making them better adapted to different climates across the world. Here, Low et al. describe a pipeline to produce haplotype-resolved, chromosome-level genomes of Angus and Brahman cattle breeds from a crossbred individual and report on comparisons of the two genomes.

Background Caucasian clover (Trifolium ambiguum M. Bieb.) is a strongly rhizomatous, low-crowned perennial leguminous and ground-covering grass. The species may be used as an ornamental plant and is resistant to cold, arid temperatures and grazing due to a well-developed underground rhizome system and a strong clonal reproduction capacity. However, the posttranscriptional mechanism of the development of the rhizome system in caucasian clover has not been comprehensively studied. Additionally, a reference genome for this species has not yet been published, which limits further exploration of many important biological processes in this plant. Result We adopted PacBio sequencing and Illumina sequencing to identify differentially expressed genes (DEGs) in five tissues, including taproot (T1), horizontal rhizome (T2), swelling of taproot (T3), rhizome bud (T4) and rhizome bud tip (T5) tissues, in the caucasian clover rhizome. In total, we obtained 19.82 GB clean data and 80,654 nonredundant transcripts were analysed. Additionally, we identified 78,209 open reading frames (ORFs), 65,227 coding sequences (CDSs), 58,276 simple sequence repeats (SSRs), 6821 alternative splicing (AS) events, 2429 long noncoding RNAs (lncRNAs) and 4501 putative transcription factors (TFs) from 64 different families. Compared with other tissues, T5 exhibited more DEGs, and co-upregulated genes in T5 are mainly annotated as involved in phenylpropanoid biosynthesis. We also identified betaine aldehyde dehydrogenase (BADH) as a highly expressed gene-specific to T5. A weighted gene co-expression network analysis (WGCNA) of transcription factors and physiological indicators were combined to reveal 11 hub genes (MEgreen-GA3), three of which belong to the HB-KNOX family, that are up-regulated in T3. We analysed 276 DEGs involved in hormone signalling and transduction, and the largest number of genes are associated with the auxin (IAA) signalling pathway, with significant up-regulation in T2 and T5. Conclusions This study contributes to our understanding of gene expression across five different tissues and provides preliminary insight into rhizome growth and development in caucasian clover.

Background Mepiquat chloride (DPC) is a chemical that is extensively used to control internode growth and create compact canopies in cultured plants. Previous studies have suggested that DPC could also inhibit gibberellin biosynthesis in sugarcane. Unfortunately, the molecular mechanism underlying the suppressive effects of DPC on plant growth is still largely unknown. Results In the present study, we first obtained high-quality long transcripts from the internodes of sugarcane using the PacBio Sequel System. A total of 72,671 isoforms, with N50 at 3073, were generated. These long isoforms were used as a reference for the subsequent RNA-seq. Afterwards, short reads generated from the Illumina HiSeq 4000 platform were used to compare the differentially expressed genes in both the DPC and the control groups. Transcriptome profiling showed that most significant gene changes occurred after six days post DPC treatment. These genes were related to plant hormone signal transduction and biosynthesis of several metabolites, indicating that DPC affected multiple pathways, in addition to suppressing gibberellin biosynthesis. The network of DPC on the key stage was illustrated by weighted gene co-expression network analysis (WGCNA). Among the 36 constructed modules, the top positive correlated module, at the stage of six days post spraying DPC, was sienna3. Notably, Stf0 sulfotransferase, cyclin-like F-box, and HOX12 were the hub genes in sienna3 that had high correlation with other genes in this module. Furthermore, the qPCR validated the high accuracy of the RNA-seq results. Conclusion Taken together, we have demonstrated the key role of these genes in DPC-induced growth inhibition in sugarcane.

AlthoughPsathyrostachys huashanicahas excellent potential for resistance gene mining and molecular genetic breeding, no reference genome is available. To date, most studies ofP. huashanicahave been focused on the creation of translocation lines and additional lines, as well as the development of molecular markers. Therefore, research at the transcriptional level is lacking. In this study, the full-length transcriptome ofP. huashanicawas sequenced using PacBio isoform sequencing (Iso-Seq) of a pooled RNA sample to explore the potential full-length transcript isoforms. We obtained 112,596 unique transcript isoforms with a total length of 114,957,868 base pairs (bp). Subsequently, Illumina sequencing reads were used to correct and trim the PacBio isoforms. We annotated 103,875 unigenes in at least one functional database, and identified a plethora of differentially-expressed genes (DEGs) that are involved in the defense responses ofP. huashanicaagainst barley yellow dwarf virus-GAV (BYDV-GAV). Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that these DEGs were mostly involved in plant-pathogen interaction, plant hormone signal transduction, and the mitogen-activated protein kinase (MAPK) signaling pathway. Additionally, we selected twenty of the RNA-seq identified resistance-related up-regulated genes, including MAPKs, cysteine-rich receptor-like protein kinases (CRPKs), calcium-dependent protein kinases (CDPKs), pathogenesis-related protein (PR) proteins, WRKYs, and disease resistance proteins, and validated their up-regulation in response to BYDV-GAV by quantitative real-time PCR. Our results indicate that a series of defense-related genes were induced inP. huashanicaduring BYDV-GAV infection. The full-length transcriptome dataset will contribute to improved use of stress-resistance genes ofP. huashanica, and serves as a reference database for the analysis of transcript expression inP. huashanica.

Whole genome duplication events have occurred frequently during the course of vertebrate evolution. To better understand the influence of polyploidization on the fish genome, we herein used the autotetraploid Carassius auratus (4n = 200, RRRR) (4nRR) resulting from the whole genome duplication of Carassius auratus (2n = 100, RR) (RCC) to explore the genomic and epigenetic alterations after polyploidization. We subsequently performed analyses of full-length transcriptome dataset, amplified fragment length polymorphism (AFLP) and methylation sensitive amplification polymorphism (MSAP) on 4nRR and RCC. By matching the results of 4nRR and RCC isoforms with reference genome in full-length transcriptome dataset, 649 and 1,971 novel genes were found in the RCC and 4nRR full-length geneset, respectively. Compared to Carassius auratus and Megalobrama amblycephala, 4nRR presented 3,661 unexpressed genes and 2,743 expressed genes. Furthermore, GO enrichment analysis of expressed genes in 4nRR revealed that they were enriched in meiosis I, whereas KEGG enrichment analysis displayed that they were mainly enriched in proteasome. Using AFLP analysis, we noted that 32.61% of RCC fragments had disappeared, while 32.79% of new bands were uncovered in 4nRR. Concerning DNA methylation, 4nRR exhibited a lower level of global DNA methylation than RCC. Additionally, 60.31% of methylation patterns in 4nRR were altered compared to RCC. These observations indicated that transcriptome alterations, genomic changes and regulation of DNA methylation levels and patterns had occurred in the newly established autotetraploid genomes, suggesting that genetic and epigenetic alterations were influenced by autotetraploidization. In summary, this study provides valuable novel insights into vertebrate genome evolution and generates relevant information for fish breeding.

BackgroundGnetumis an economically important tropical and subtropical gymnosperm genus with various dietary, industrial and medicinal uses. Many carbohydrates, proteins and fibers accumulate during the ripening ofGnetumseeds. However, the molecular mechanisms related to this process remain unknown.ResultsWe therefore assembled a full-length transcriptome from immature and matureG. luofuenseseeds using PacBio sequencing reads. We identified a total of 5726 novel genes, 9061 alternative splicing events, 3551 lncRNAs, 2160 transcription factors, and we found that 8512 genes possessed at least one poly(A) site. In addition, gene expression comparisons of six transcriptomes generated by Illumina sequencing showed that 14,323 genes were differentially expressed from an immature stage to a mature stage with 7891 genes upregulated and 6432 genes downregulated. The expression of 14 differentially expressed transcription factors from the MADS-box, Aux/IAA and bHLH families was validated by qRT-PCR, suggesting that they may have important roles in seed ripening ofG. luofuense.ConclusionsThese findings provide a valuable molecular resource for understanding seed development of gymnosperms.

Background Red palm weevil Rhynchophorus ferrugineus (Coleoptera: Curculionidae) is one of the most destructive insects for palm trees in the world. However, its genome resources are still in the blank stage, which limits the study of molecular and growth development analysis. Methods In this study, we used PacBio Iso-Seq and Illumina RNA-seq to first generate transcriptome from three developmental stages of R. ferrugineus (pupa, 7th larva, female and male) to increase our understanding of the life cycle and molecular characteristics of R. ferrugineus. Results A total of 63,801 nonredundant full-length transcripts were generated with an average length of 2,964 bp from three developmental stages, including the 7th instar larva, pupa, female adult and male adult. These transcripts showed a high annotation rate in seven public databases, with 54,999 (86.20%) successfully annotated. Meanwhile, 2,184 alternative splicing (AS) events, 2,084 transcription factors (TFs), 66,230 simple sequence repeats (SSR) and 9,618 Long noncoding RNAs (lncRNAs) were identified. In summary, our results provide a new source of full-length transcriptional data and information for the further study of gene expression and genetics in R. ferrugineus.

Seed-borne endophyte Epichloë gansuensis enhance NaCl tolerance in Achnatherum inebrians and increase its biomass. However, the molecular mechanism by which E. gansuensis increases the tolerance of host grasses to NaCl stress is unclear. Hence, we firstly explored the full-length transcriptome information of A. inebrians by PacBio RS II. In this work, we obtained 738,588 full-length non-chimeric reads, 36,105 transcript sequences and 27,202 complete CDSs from A. inebrians. We identified 3558 transcription factors (TFs), 15,945 simple sequence repeats and 963 long non-coding RNAs of A. inebrians. The present results show that 2464 and 1817 genes were differentially expressed by E. gansuensis in the leaves of E+ and E− plants at 0 mM and 200 mM NaCl concentrations, respectively. In addition, NaCl stress significantly regulated 4919 DEGs and 502 DEGs in the leaves of E+ and E− plants, respectively. Transcripts associated with photosynthesis, plant hormone signal transduction, amino acids metabolism, flavonoid biosynthetic process and WRKY TFs were differentially expressed by E. gansuensis; importantly, E. gansuensis up-regulated biology processes (brassinosteroid biosynthesis, oxidation–reduction, cellular calcium ion homeostasis, carotene biosynthesis, positive regulation of proteasomal ubiquitin-dependent protein catabolism and proanthocyanidin biosynthesis) of host grass under NaCl stress, which indicated an increase in the ability of host grasses’ adaptation to NaCl stress. In conclusion, our study demonstrates the molecular mechanism for E. gansuensis to increase the tolerance to salt stress in the host, which provides a theoretical basis for the molecular breed to create salt-tolerant forage with endophytes.

Several main families of Ranunculales are rich in alkaloids and other medicinal compounds; many species of these families are used in traditional and folk medicine. Dichocarpum is a representative medicinal genus of Ranunculaceae, but the genetic basis of its metabolic phenotype has not been investigated, which hinders its sustainable conservation and utilization. We use the third-generation high-throughput sequencing and metabolomic techniques to decipher the full-length transcriptomes and metabolomes of five Dichocarpum species endemic in China, and 71,598 non-redundant full-length transcripts were obtained, many of which are involved in defense, stress response and immunity, especially those participating in the biosynthesis of specialized metabolites such as benzylisoquinoline alkaloids (BIAs). Twenty-seven orthologs extracted from trancriptome datasets were concatenated to reconstruct the phylogenetic tree, which was verified by the clustering analysis based on the metabolomic profile and agreed with the Pearson correlation between gene expression patterns of Dichocarpum species. The phylogenomic analysis of phytometabolite biosynthesis genes, e.g., (S)-norcoclaurine synthase, methyltransferases, cytochrome p450 monooxygenases, berberine bridge enzyme and (S)-tetrahydroprotoberberine oxidase, revealed the evolutionary trajectories leading to the chemodiversity, especially that of protoberberine type, aporphine type and bis-BIA abundant in Dichocarpum and related genera. The biosynthesis pathways of these BIAs are proposed based on full-length transcriptomes and metabolomes of Dichocarpum. Within Ranunculales, the gene duplications are common, and a unique whole genome duplication is possible in Dichocarpum. The extensive correlations between metabolite content and gene expression support the co-evolution of various genes essential for the production of different specialized metabolites. Our study provides insights into the transcriptomic and metabolomic landscapes of Dichocarpum, which will assist further studies on genomics and application of Ranunculales plants.

Since the last century, episodes of coral reef bleaching and mortality have occurred almost annually in tropical or subtropical seas. When the temperature exceeds the tolerant limit of a coral–zooxanthellae holobiont, it induces physiological stress and disrupts the vulnerable fine-tuned balance between the two partners, leading to bleaching. The gene expression profiles of a scleractinian coral and its symbiotic zooxanthellae can offer important information with which to decipher this balanced relationship at the functional level of genes. Here, we sequence a full-length transcriptome of a well-known, common and frequently dominant reef-building coral, Pocillopora damicornis, to acquire gene expression information for the coral–zooxanthellae holobiont. To this end, we identify 21,926 and 465 unique genes in the coral and algal symbiont, respectively, and examine the functional enrichment among these genes based on GO (gene ontology) terms and KEGG (the Kyoto Encyclopedia of Genes and Genomes) pathways. The results show that the zooxanthellae provide for their coral host through energy and nutrition metabolism by photosynthesis, and that both the coral host and zooxanthellae have an anti-stress molecular mechanism, though the two parties have independent abilities to survive in the short term. This work sheds light on the valuable gene expression profile of a coral–zooxanthellae holobiont and provides grounds for further molecular biological research to support ecological protection work.

Background Spodoptera frugiperda (J. E. Smith), commonly known as fall armyworm (FAW), is one of the most destructive agricultural pests in the world and has posed a great threat to crops. The improper use of insecticides has led to rapid development of resistance. However, the genetic data available for uncovering the insecticide resistance mechanisms are scarce. Methods In this study, we used PacBio single-molecule real-time (SMRT) sequencing aimed at revealing the full-length transcriptome profiling of the FAW larval brain to obtain detoxification genes. Results A total of 18,642 high-quality transcripts were obtained with an average length of 2,371 bp, and 11,230 of which were successfully annotated in six public databases. Among these, 5,692 alternative splicing events were identified.

Leaves of C. porrectum are rich in essential oils containing monoterpenes, sesquiterpenes and aromatic compounds, but the molecular mechanism of terpenoid biosynthesis in C. porrectum is still unclear. In this paper, the differences in the contents and compositions of terpenoids among three chemotypes were analyzed using gas chromatography mass spectrometry (GC/MS). Furthermore, the differential expression of gene transcripts in the leaf tissues of the three C. porrectum chemotypes were analyzed through a comparison of full-length transcriptomes and expression profiles. The essential oil of the three C. porrectum chemotypes leaves was mainly composed of monoterpenes. In the full-length transcriptome of C. porrectum, 104,062 transcripts with 306,337,921 total bp, an average length of 2944 bp, and an N50 length of 5449 bp, were obtained and 94025 transcripts were annotated. In the eucalyptol and linalool chemotype, the camphor and eucalyptol chemotype, and the camphor and linalool chemotype comparison groups, 21, 22 and 18 terpene synthase (TPS) unigenes were identified respectively. Three monoterpene synthase genes, CpTPS3, CpTPS5 and CpTPS9, were upregulated in the eucalyptol chemotype compared to the linalool chemotype and camphor chemotype. CpTPS1 was upregulated in the camphor chemotype compared to the linalool chemotype and the eucalyptol chemotype. CpTPS4 was upregulated in the linalool chemotype compared to the camphor chemotype and the eucalyptol chemotype. Different unigenes had different expression levels among the three chemotypes, but the unigene expression levels of the 2-C-methyl-D-erythritol 4phosphate (MEP) pathway were generally higher than those of the mevalonate acid (MVA) pathway. Quantitative reverse transcription PCR(qRT-PCR) further validated these expression levels. The present study provides new clues for the functional exploration of the terpenoid synthesis mechanism and key genes in different chemotypes of C. porrectum.

Fenneropenaeus chinensis (F. chinensis) is one of the most commercially important cultured shrimps in China. The adult F. chinensis exhibit sexual dimorphism in growth and body color. In this research, we profiled the whole transcriptome of F. chinensis by using single molecule real-time-based full-length transcriptome sequencing. We further performed Illumina-based short reads RNA-seq on muscle and gonad of two sexes to detect the sex-biased expression genes. In muscle, we observed significantly more female-biased transcripts. With the differentially expressed transcripts (DETs) in muscle, some pathways related to the energy metabolism were enriched, which may be responsible for the difference of growth. We also digged out a pathway named porphyrin and chlorophyll metabolism. It was speculated to relevant to the difference of body color between the two sexes of shrimp. Interestingly, almost all DETs in these pathways were female-biased expression in muscle, which could explain the phenomenon of better growth performance and darker body color in female. In gonad, several pathways involved in reproduction were enriched. For instance, some female-biased DETs participated in the arachidonic acid metabolism, which was reported crucial in female reproduction. In conclusion, our studies identified abundant sex-biased expression transcripts and important pathways involved in sexual dimorphism by using the RNA-seq method. It provided a basis for future researches on the sexual dimorphism of F. chinensis.

Nitraria tangutorum Bobrov is a halophyte that is resistant to salt and alkali and is widely distributed in northwestern China. However, its genome has not been sequenced, thereby limiting studies on this particular species. For species without a reference genome, the full-length transcriptome is a convenient and rapid way to obtain reference gene information. To better study N. tangutorum, we used PacBio single-molecule real-time technology to perform full-length transcriptome analysis of this halophyte. In this study, a total of 21.83 Gb of data were obtained, and 198,300 transcripts, 51,875 SSRs (simple sequence repeats), 55,574 CDS (coding sequence), and 74,913 lncRNAs (long non-coding RNA) were identified. In addition, using this full-length transcriptome, we identified the key Na+/H+ antiporter (NHX) genes that maintain ion balance in plants and found that these are induced to express under salt stress. The results indicate that the full-length transcriptome of N. tangutorum can be used as a database and be utilized in elucidating the salt tolerance mechanism of N. tangutorum.

The white-backed planthopper Sogatella furcifera is an economically important rice pest distributed throughout Asia. It damages rice crops by sucking phloem sap, resulting in stunted growth and plant virus transmission. We aimed to obtain the full-length transcriptome data of S. furcifera using PacBio single-molecule real-time (SMRT) sequencing. Total RNA extracted from S. furcifera at various developmental stages (egg, larval, and adult stages) was mixed and used to generate a full-length transcriptome for SMRT sequencing. Long non-coding RNA (lncRNA) identification, full-length coding sequence prediction, full-length non-chimeric (FLNC) read detection, simple sequence repeat (SSR) analysis, transcription factor detection, and transcript functional annotation were performed. A total of 12,514,449 subreads (15.64 Gbp, clean reads) were generated, including 630,447 circular consensus sequences and 388,348 FLNC reads. Transcript cluster analysis of the FLNC reads revealed 251,109 consensus reads including 29,700 high-quality reads. Additionally, 100,360 SSRs and 121,395 coding sequences were identified using SSR analysis and ANGEL software, respectively. Furthermore, 44,324 lncRNAs were annotated using four tools and 1,288 transcription factors were identified. In total, 95,495 transcripts were functionally annotated based on searches of seven different databases. To the best of our knowledge, this is the first study of the full-length transcriptome of the white-backed planthopper obtained using SMRT sequencing. The acquired transcriptome data can facilitate further studies on the ecological and viral-host interactions of this agricultural pest.

The amphipod Gammarus lacustris has been distributing in the Tibetan region with well-known uplifts of the Tibetan plateau. It is hence considered as a good model for investigating stress adaptations of the plateau. Here, we sequenced the whole-genome and full-length transcriptome of G. lacustris, and compared the transcriptome results with its counterpart Gammarus pisinnus from a nearby plain. Our main goal was to provide a genomic resource for investigation of genetic mechanisms, by which G. lacustris adapted to living on the plateau. The final draft genome assembly of G. lacustris was 5.07 gigabases (Gb), and it contained 443,304 scaffolds (>2 kb) with an N50 of 2,578 bp. A total of 8,858 unigenes were predicted in the full-length transcriptome of G. lacustris, with an average gene length of 1,811 bp. Compared with the G. pisinnus transcriptome, 2,672 differentially expressed genes (DEGs) were up-regulated and 2,881 DEGs were down-regulated in the G. lacustris transcriptome. Along with these critical DEGs, several enriched metabolic pathways, such as oxidative phosphorylation, ribosome, cell energy homeostasis, glycolysis and gluconeogenesis, were predicted to play essential roles in the plateau adaptation. In summary, the present study provides a genomic basis for understanding the plateau adaption of G. lacustris, which lays a fundamental basis for further biological and ecological studies on other resident aquatic species in the Tibetan plateau.

Background Poxviruses are large DNA viruses that infect humans and animals. Vaccinia virus (VACV) has been applied as a live vaccine for immunization against smallpox, which was eradicated by 1980 as a result of worldwide vaccination. VACV is the prototype of poxviruses in the investigation of the molecular pathogenesis of the virus. Short-read sequencing methods have revolutionized transcriptomics; however, they are not efficient in distinguishing between the RNA isoforms and transcript overlaps. Long-read sequencing (LRS) is much better suited to solve these problems and also allow direct RNA sequencing. Despite the scientific relevance of VACV, no LRS data have been generated for the viral transcriptome to date. Findings For the deep characterization of the VACV RNA profile, various LRS platforms and library preparation approaches were applied. The raw reads were mapped to the VACV reference genome and also to the host (Chlorocebus sabaeus) genome. In this study, we applied the Pacific Biosciences RSII and Sequel platforms, which altogether resulted in 937,531 mapped reads of inserts (1.42 Gb), while we obtained 2,160,348 aligned reads (1.75 Gb) from the different library preparation methods using the MinION device from Oxford Nanopore Technologies. Conclusions By applying cutting-edge technologies, we were able to generate a large dataset that can serve as a valuable resource for the investigation of the dynamic VACV transcriptome, the virus-host interactions, and RNA base modifications. These data can provide useful information for novel gene annotations in the VACV genome. Our dataset can also be used to analyze the currently available LRS platforms, library preparation methods, and bioinformatics pipelines.