Single-molecule long-read sequencing of the full-length transcriptome of Rhododendron lapponicum L

Exploring the innate immune system of Urechis unicinctus: Insights from full-length transcriptome analysis

The Echiura worm Urechis unicinctus refers to a common benthic invertebrate found in the intertidal zone of Huanghai as well as Bohai Bay. U. unicinctus is known to contain various physiologically active substances, making it highly valuable in terms of its edibility, medicinal properties, and economic potential. Nonetheless, the limited study on the immune system of U. unicinctus poses difficulties for its aquaculture and artificial reproduction. Marine invertebrates, including shellfish and U. unicinctus, are thought to primarily depend on their innate immune system for disease protection, owing to the severalinnate immune molecules they possess. Herein, we employed PacBio single-molecule real-time (SMRT) sequencing technology to perform the full-length transcriptome analysis of U. unicinctus individuals under five different conditions (room temperature (RT), low temperature (LT), high temperature (HT), without water (DRY), ultraviolet irradiation (UV)). Concequently, we identified 59,371 unigenes that had a 2,779 bp average length, 2,613 long non-coding RNAs (lncRNAs), 59,190 coding sequences (CDSs), 35,166 simple sequence repeats (SSRs), and 1,733 transcription factors (TFs), successfully annotating 90.58 % (53,778) of the unigenes. Subsequently, key factors associated with immune-related processes, such as non-self-recognition, cellular immune defenses, and humoral immune defenses, were searched. Our study also identified pattern recognition receptors (PRRs) that included 17 peptidoglycan recognition proteins (PGRPs), 13 Gram-negative binding proteins (GNBPs), 18 scavenger receptors (SRs), 74 toll-like receptors (TLRs), and 89 C-type lectins (CLTs). Altogether, the high-quality transcriptome obtained data will offer valuable insights for further investigations into U. unicinctus innate immune response, laying the foundation for subsequent molecular biology studies and aquaculture.

Analysis of the Spatholobus suberectus full-length transcriptome identified an R2R3-MYB transcription factor-encoding gene SsMYB158 that regulates flavonoid biosynthesis

Spatholobus suberectus Dunn (Leguminosae) has been used for medicinal purposes for a long period. Flavonoids are the major bioactive components of S. suberectus. However, there is still limited knowledge of the exact method via which transcription factors (TFs) regulate flavonoid biosynthesis. The full-length transcriptome of S. suberectus was analyzed using SMRT sequencing; 61,548 transcripts were identified, including 12,311 new gene loci, 53,336 novel transcripts, 44,636 simple sequence repeats, 36,414 complete coding sequences, 871 long non-coding RNAs and 6781 TFs. The SsMYB158 TF, which is associated with flavonoid biosynthesis, belongs to the R2R3-MYB class and is localized subcellularly to the nucleus. The overexpression of SsMYB158 in Nicotiana benthamiana and the transient overexpression of SsMYB158 in S. suberectus resulted in a substantial enhancement in both flavonoids and catechin levels. In addition, there was a remarkable upregulation in the expression of essential enzyme-coding genes associated with the flavonoid biosynthesis pathways. Our study revealed SsMYB158 as a critical regulator of flavonoid biosynthesis in S. suberectus and laying the foundation for its molecular breeding.

SSR Marker Acquisition and Application from Transcriptome of Captive Chinese Forest Musk Deer (Moschus berezovskii)

Forest musk deer (Moschus berezovskii) is one of the most endangered medicinally important wild animals in the world. Forest musk deer farming is the main way of production of musk. However, the single provenance and lack of genetic information lead to reduced genetic diversity of forest musk deer. Therefore, more SSR markers need to be developed to identify forest musk deer germplasm. In this study, bone marrow derived mesenchymal cells were used to construct cDNA library for transcriptome sequencing. The datasets were de novo assembled and annotated. 9 polymorphic simple sequence repeat (SSR) markers were finally identified and used to detect population genetic diversity. 6.07 Gb clean data were generated using Illumina sequencing technology, and de novo assembled into 138,591 transcripts and 81,553 unigenes. 5,777 simple sequence repeats (SSRs) were identified, in which there were 578 repeating motif types, with mono-nucleotide and tri-nucleotides comprising 55.88% and 25.60%, respectively. 100 primer pairs were designed to validate amplification and polymorphism using DNA from fecal samples. 9 polymorphic SSRs were developed and used to detect population genetic diversity of 122 forest musk deer in 2 farms. The average number of alleles per locus varied from 4 to 15 (average = 8.3). The observed heterozygosity (HO) per locus ranged from 0.102 to 0.941, while the expected heterozygosity (HE) per locus was from 0.111 to 0.651. All loci deviated significantly from the Hardy-Weinberg equilibrium (p < 0.001). The polymorphism information content (PIC) of these loci varied from 0.108 to 0.619. 9 polymorphic SSR markers were developed in this research. These sites can be used for breeding planning and conservation of germplasm resources.

Recent advancements in the physiological, genetic, and genomic research on Rhododendrons for trait improvement

High species diversity, hybridization potential, broad geographical dispersal range and ornamental characteristics (i.e., attractive size, shape, structure, flowers, and evergreen) have fetched a good international market for Rhododendron. However, most species are restricted to specific geographic areas due to their habitat specificity in acidic soil and cold climates, resulting many species being classified under threat categories of the IUCN. In this review, advances in research on Rhododendron for improvement to floral display quality and stress resistance have been described. The low genetic barrier among species has created opportunities for extensive hybridization and ploidy alteration for introducing quality and adaptive traits during the development of new varieties. Recent technological advances have supported investigations into the mechanism of flower development, as well as cold tolerance and pathogen resistance mechanisms in the Rhododendron. However, most of the species have limited adaptability to drought, line-tolerance, pathogen resistance, and high-temperature conditions and this resistance ability present in few species largely remains unexplored. Additionally, the available genetic diversity and genomic information on species, and possibilities for their application in molecular breeding have been summarized. Overall, genomic resource data are scarce in the majority of the members of this genus. Finally, various research gaps such as genetic mapping of quality traits, understanding the molecular mechanism of quality-related traits and genomic assortment in Rhododendron members have been discussed in the future perspective section.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-04006-6.

Characterization and analysis of multi-organ full-length transcriptomes in Sphaeropteris brunoniana and Alsophila latebrosa highlight secondary metabolism and chloroplast RNA editing pattern of tree ferns

BACKGROUND

Sphaeropteris brunoniana and Alsophila latebrosa are both old relict and rare tree ferns, which have experienced the constant changes of climate and environment. However, little is known about their high-quality genetic information and related research on environmental adaptation mechanisms of them. In this study, combined with PacBio and Illumina platforms, transcriptomic analysis was conducted on the roots, rachis, and pinna of S. brunoniana and A. latebrosa to identify genes and pathways involved in environmental adaptation. Additionally, based on the transcriptomic data of tree ferns, chloroplast genes were mined to analyze their gene expression levels and RNA editing events.

RESULTS

In the study, we obtained 11,625, 14,391 and 10,099 unigenes of S. brunoniana root, rachis, and pinna, respectively. Similarly, a total of 13,028, 11,431 and 12,144 unigenes were obtained of A. latebrosa root, rachis, and pinna, respectively. According to the enrichment results of differentially expressed genes, a large number of differentially expressed genes were enriched in photosynthesis and secondary metabolic pathways of S. brunoniana and A. latebrosa. Based on gene annotation results and phenylpropanoid synthesis pathways, two lignin synthesis pathways (H-lignin and G-lignin) were characterized of S. brunoniana. Among secondary metabolic pathways of A. latebrosa, three types of WRKY transcription factors were identified. Additionally, based on transcriptome data obtained in this study, reported transcriptome data, and laboratory available transcriptome data, positive selection sites were identified from 18 chloroplast protein-coding genes of four tree ferns. Among them, RNA editing was found in positive selection sites of four tree ferns. RNA editing affected the protein secondary structure of the rbcL gene. Furthermore, the expression level of chloroplast genes indicated high expression of genes related to the chloroplast photosynthetic system in all four species.

CONCLUSIONS

Overall, this work provides a comprehensive transcriptome resource of S. brunoniana and A. latebrosa, laying the foundation for future tree fern research.

Full-length transcriptome characterization of Platycladus orientalis based on the PacBio platform

As a unique and native conifer in China, Platycladus orientalis is widely used in soil erosion control, garden landscapes, timber, and traditional Chinese medicine. However, due to the lack of reference genome and transcriptome, it is limited to the further molecular mechanism research and gene function mining. To develop a full-length reference transcriptome, tissues from five different parts of P. orientalis and four cone developmental stages were sequenced and analyzed by single-molecule real-time (SMRT) sequencing through the PacBio platform in this study. Overall, 37,111 isoforms were detected by PacBio with an N50 length of 2,317 nt, an average length of 1,999 bp, and the GC content of 41.81%. Meanwhile, 36,120 coding sequences, 5,645 simple sequence repeats (SSRs), 1,201 non-coding RNAs (lncRNAs), and 182 alternative splicing (AS) events with five types were identified using the results obtained from the PacBio transcript isoforms. Furthermore, 1,659 transcription factors (TFs) were detected and belonged to 51 TF families. A total of 35,689 transcripts (96.17%) were annotated through the NCBI nr, KOG, Swiss-Prot and KEGG databases, and 385 transcript isoforms related to 8 types of hormones were identified incorporated into plant hormone signal transduction pathways. The assembly and revelation of the full-length transcriptome of P. orientalis offer a pioneering insight for future investigations into gene function and genetic breeding within Platycladus species.

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.

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 (O₂^·-) content and production rate significantly decreased, while malondialdehyde (MDA) and hydrogen peroxide (H₂O₂) 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.

Integrative Analysis of Metabolome and Transcriptome Reveals Molecular Insight into Metabolomic Variations during Hawthorn Fruit Development

Hawthorn (Crataegus pinnatifida var. major), a cultivated fruit tree, is native and unique to China. Its fruits have high nutritional, health, and medicinal values. However, the development and ripening process of hawthorns is accompanied by dramatic changes in flavor, aroma, and bioactive phytochemicals, which are the fundamental factors that contribute to the potential health benefits and establishment of fruit quality. Therefore, an exploration of the dynamic changes in metabolites and their regulatory networks during the development of hawthorn fruits can elucidate the formation mechanisms of active substances in hawthorn fruits. In this study, we used a broad targeted metabolomics approach to identify and analyze the dynamics of metabolites in hawthorn fruits at five developmental stages. The results revealed 998 primary and secondary metabolites that were classified into 15 categories. The accumulation levels of most sugars increased during fruit development and then accelerated at the fruit ripening stage. The accumulation levels of a few organic acids (e.g., citric acid, isocitric acid, and quinic acid) continuously increased. Many organic acids exhibited significant decreasing trends. Among the 561 secondary metabolites detected, 189 were phenolic acids and 199 were flavonoids. The levels of many flavonoids were significantly reduced at later stages of fruit development; in contrast, the levels of two anthocyanins significantly increased during fruit ripening. Correlation analysis revealed that there is a certain correlation within and between primary as well as secondary metabolites during fruit development. Furthermore, the integration of metabolomic and transcriptomic data in this study revealed that changes in the expression of some differentially expressed genes (DEGs) were associated with the accumulation of metabolites such as sugars, organic acids, and flavonoids, e.g., the upregulated expression levels of CS (citrate synthase) genes were consistent with the continued accumulation of citric acid. Overall, this study demonstrates the metabolic changes that occur during the development of hawthorn fruit, explores the molecular mechanisms that underlie metabolite changes during fruit development, and lays a strong theoretical foundation for the improvement of hawthorn fruit quality and the development of functional components.

Screening and identification of multiple abiotic stress responsive candidate genes based on hybrid-sequencing in Vicia sativa

Common vetch is an important leguminous forage for both livestock fodder and green manure and has a tremendous latent capacity in a sustainable agroecosystem. In the present study, a comprehensive transcriptome analysis of the aboveground leaves and underground roots of common vetch under multiple abiotic stress treatments, including NaCl, drought, cold, and cold drought, was performed using hybrid-sequencing technology, i. e. single-molecule real-time sequencing technology (SMRT) and supplemented by next-generation sequencing (NGS) technology. A total of 485,038 reads of insert (ROIs) with a mean length of 2606 bp and 228,261 full-length nonchimeric (FLNC) reads were generated. After deduplication, 39,709 transcripts were generated. Of these transcripts, we identified 1059 alternative splicing (AS) events, 17,227 simple sequence repeats (SSRs), and 1647 putative transcription factors (TFs). Furthermore, 640 candidates long noncoding RNAs (lncRNAs) and 28,256 complete coding sequences (CDSs) were identified. In gene annotation analyses, a total of 38,826 transcripts (97.78%) were annotated in eight public databases. Finally, seven multiple abiotic stress-responsive candidate genes were obtained through gene expression, annotation information, and protein-protein interaction (PPI) networks. Our research not only enriched the structural information of FL transcripts in common vetch, but also provided useful information for exploring the molecular mechanism of multiple abiotic stress tolerance between aboveground and underground tissues in common vetch and related legumes.

Full-length transcriptome analysis of Zanthoxylum nitidum (Roxb.) DC

Zanthoxylum nitidum (Roxb.) DC. (Z. nitidum) is a type of Chinese Dao-di herb, also called Liangmianzhen, which is widely used to treat arthralgia, rheumatic arthralgia, and stomach pain. However, genomic resources for Z. nitidum are still scarce. This study provides transcriptomic resources for Z. nitidum by applying single-molecule real-time (SMRT) sequencing technology. In total, 456,109 circular consensus sequencing (CCS) reads were generated with a mean length of 2,216 bp from Z. nitidum roots, old stems, young branches, leaves, flowers, and fruits. Of these total reads, 353,932 were full-length nonchimeric (FLNC) reads with an average length of 1,996 bp. A total of 16,163 transcripts with a mean length of 1,171 bp were acquired. Of these transcripts, 14,231 (88%) were successfully annotated using public databases. Across all the 16,163 transcripts, we identified 6,255 long non-coding RNAs (lncRNAs) and 22,780 simple sequence repeats (SSRs). Furthermore, 3,482 transcription factors were identified. Among the SSR loci, 1-3 nucleotide repeats were dominant, occupying 99.36% of the total SSR loci, with mono-, di-, and tri-nucleotide repeats accounting for 61.80%, 19.89%, and 5.02% of the total SSR loci, respectively. A total of 36 out of 100 randomly selected primer pairs were verified to be positive, 20 of which showed polymorphism. These findings enrich the genetic resources available for facilitating future studies and research on relevant topics such as population genetics in Z. nitidum.

Microsatellite analysis and polymorphic marker development based on the full-length transcriptome of Camellia chekiangoleosa

Camellia chekiangoleosa is a popular variety of Oil-camellia that has high oil production and ornamental value. Microsatellite (SSR) markers are the preferred tool for the molecular marker-assisted breeding of C. chekiangoleosa. By focusing on the problems of the low development efficiency of polymorphic SSR markers and the lack of available functional markers in Oil-camellia, we identified 97,510 SSR loci based on the full-length transcriptome sequence of C. chekiangoleosa. An analysis of SSR characteristics showed that mononucleotide (51.29%) and dinucleotide (34.36%) SSRs were the main repeat types. The main SSR distribution areas based on proportion covered were ordered as follows: 5'UTR > 3'UTR > CDS. By comparing our data with those in databases such as GO and KEGG, we obtained functional annotations of unigene sequences containing SSR sites. The data showed that the amplification efficiency of the SSR primers was 51.72%, and the development efficiency of polymorphic SSR primers was 26.72%. Experiments verified that dinucleotide and pentanucleotide SSRs located in UTR regions could produce more polymorphic markers. An investigation into the genetic diversity of several C. chekiangoleosa populations also suggested that the developed SSR markers had higher levels of polymorphism. This study will provide a reference and high-quality markers for the large-scale development of functional SSR markers and genetic research in Oil-camellia.

Single-molecule real-time sequencing of the full-length transcriptome of Halophila beccarii

Ecologically, Halophila beccarii Asch. is considered as a colonizing or a pioneer seagrass species and a “tiny but mighty” seagrass species, since it may recover quickly from disturbance generally. The use of transcriptome technology can provide a better understanding of the physiological processes of seagrasses. To date, little is known about the genome and transcriptome information of H. beccarii. In this study, we used single molecule real-time (SMRT) sequencing to obtain full-length transcriptome data and characterize the transcriptome structure. A total of 11,773 of the 15,348 transcripts were successfully annotated in seven databases. In addition, 1573 long non-coding RNAs, 8402 simple sequence repeats and 2567 transcription factors were predicted in all the transcripts. A GO analysis showed that 5843 transcripts were divided into three categories, including biological process (BP), cellular component (CC) and molecular function (MF). In these three categories, metabolic process (1603 transcripts), protein-containing complex (515 transcripts) and binding (3233 transcripts) were the primary terms in BP, CC, and MF, respectively. The major types of transcription factors were involved in MYB-related and NF-YB families. To the best of our knowledge, this is the first report of the transcriptome of H. beccarii using SMRT sequencing technology.

Phytophthora cinnamomi causing root rot on Rhododendron lapponicum and control it using potential biocontrol agents

Rhododendron lapponicum (R. lapponicum) is a dwarf Rhododendron species, which is severely infected with root rot and wilt in Yunnan province, China. However, the causal agent causing these symptoms was unknown. An isolate, Pci-1 was identified as Phytophthora cinnamomi, based on its morphology and the sequences of β-tubulin, internal transcribed spacer, and Ypt1 genes and verified according to the Koch's postulate. We found that this pathogen could infect 14 species of plants, including Althaea rosea, Viburnum cylindricum, and Brassica napus. Strain Pci-1 could cause R. lapponicum to wither and die; and it grows best in an oat medium with pH 7.0 - 8.0 and an optimum temperature of 27°C. We suggest that the rhizosphere of R. lapponicum treated with biocontrol strains Paenibacillus polymyxoides P2-5 and Trichoderma asperellum Tv-1 showed a significant inhibitory effect on pathogen Pci-1. The inhibitory effect of 70% dimethomorph + cymoxanil was significantly higher with EC₅₀ and EC₉₀ values of 0.1894 and 0.3618 a.i. µg/ml, respectively. Greenhouse experiments revealed that the pathogen load is decreased in the presence of potential antagonists. This study provides fundamentals on risk assessment and theoretical support for the management of P. cinnamomi pathogen and contributes significantly to the planting of forest and horticultural crops in a disease-free environment.

Full-Length Transcriptome Characterization and Comparative Analysis of Chosenia arbutifolia

As a unique tree species in the Salicaceae family, Chosenia arbutifolia is used primarily for construction materials and landscape planting in China. Compared with other Salicaceae species members, the genomic resources of C. arbutifolia are extremely scarce. Thus, in the present study, the full-length transcriptome of C. arbutifolia was sequenced by single-molecular real-time sequencing (SMRT) technology based on the PacBio platform. Then, it was compared against those of other Salicaceae species. We generated 17,397,064 subreads and 95,940 polished reads with an average length of 1812 bp, which were acquired through calibration, clustering, and polishing. In total, 50,073 genes were reconstructed, of which 48,174 open reading frames, 4281 long non-coding RNAs, and 3121 transcription factors were discovered. Functional annotation revealed that 47,717 genes had a hit in at least one of five reference databases. Moreover, a set of 12,332 putative SSR markers were screened among the reconstructed genes. Single-copy and special orthogroups, and divergent and conserved genes, were identified and analyzed to find divergence among C. arbutifolia and the five Salicaceae species. To reveal genes involved in a specific function and pathway, enrichment analyses for GO and KEGG were also performed. In conclusion, the present study empirically confirmed that SMRT sequencing realistically depicted the C. arbutifolia transcriptome and provided a comprehensive reference for functional genomic research on Salicaceae species.

Full-length transcriptome analysis of Adiantum flabellulatum gametophyte

Ferns are important components of plant communities on earth, but their genomes are generally very large, with many redundant genes, making whole genome sequencing of ferns prohibitively expensive and time-consuming. This means there is a significant lack of fern reference genomes, making molecular biology research difficult. The gametophytes of ferns can survive independently, are responsible for sexual reproduction and the feeding of young sporophytes, and play an important role in the alternation of generations. For this study, we selected Adiantum flabellulatum as it has both ornamental and medicinal value and is also an indicator plant of acidic soil. The full-length transcriptome sequencing of its gametophytes was carried out using PacBio three-generation sequencing technology. A total of 354,228 transcripts were obtained, and 231,705 coding sequences (CDSs) were predicted, including 5,749 transcription factors (TFs), 2,214 transcription regulators (TRs) and 4,950 protein kinases (PKs). The transcripts annotated by non-redundant protein sequence database (NR), Kyoto encyclopedia of genes and genomes (KEGG), eukaryotic ortholog groups (KOG), Swissprot, protein family (Pfma), nucleotide sequence database (NT) and gene ontology (GO) were 251,501, 197,474, 193,630, 194,639, 195,956, 113,069 and 197,883, respectively. In addition, 138,995 simple sequence repeats (SSRs) and 111,793 long non-coding RNAs (lncRNAs) were obtained. We selected nine chlorophyll synthase genes for qRT-PCR, and the results showed that the full-length transcript sequences and the annotation information were reliable. This study can provide a reference gene set for subsequent gene expression quantification.

SMRT Sequencing of the Full-Length Transcriptome of the Coelomactra antiquata

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.

The chromosome-scale genome assembly, annotation and evolution of Rhododendron henanense subsp. lingbaoense

Rhododendron henanense subsp. lingbaoense (hereafter referred to as R. henanense) is an endemic species naturally distributed in the Henan province, China, with high horticultural, ornamental and medicinal value. Herein, we report a de novo genome assembly for R. henanense using a combination of PacBio long read and Illumina short read sequencing technologies. In total, we assembled 634.07 Mb with a contig N50 of 2.5 Mb, representing ~96.93% of the estimated genome size. By applying Hi-C data, 13 pseudochromosomes of R. henanense genome were assembled, covering ~98.21% of the genome assembly. The genome was composed of ~65.76% repetitive sequences and 31,098 protein-coding genes, 88.77% of which could be functionally annotated. Rhododendron henanense displayed a high level of synteny with other Rhododendron species from the Hymenanthes subgenus. Our data also suggests that R. henanense genes related to stress responses have undergone expansion, which may underly the unique abiotic and biotic stress resistance of the species. This alpine Rhododendron chromosome-scale genome assembly provides fundamental molecular resources for germplasm conservation, breeding efforts, evolutionary studies, and elucidating the unique biological characteristics of R. henanense.

Transcriptome Profile Analysis of Arabidopsis Reveals the Drought Stress-Induced Long Non-coding RNAs Associated With Photosynthesis, Chlorophyll Synthesis, Fatty Acid Synthesis and Degradation

Long non-coding RNAs (lncRNAs) play an important role in the response of plants to drought stress. The previous studies have reported that overexpression of LEA3 and VOC could enhance drought tolerance and improve the oil content in Brassica napus and Arabidopsis thaliana, and most of the efforts have been invested in the gene function analysis, there is little understanding of how genes that involved in these important pathways are regulated. In the present study, the transcriptomic results of LEA3 and VOC over-expressed (OE) lines were compared with the RNAi lines, mutant lines and control lines under long-term and short-term drought treatment, a series of differentially expressed lncRNAs were identified, and their regulation patterns in mRNA were also investigated in above mentioned materials. The regulation of the target genes of differentially expressed lncRNAs on plant biological functions was studied. It was revealed that the mutant lines had less drought-response related lncRNAs than that of the OE lines. Functional analysis demonstrated that multiple genes were involved in the carbon-fixing and chlorophyll metabolism, such as CDR1, CHLM, and CH1, were regulated by the upregulated lncRNA in OE lines. In LEA-OE, AT4G13180 that promotes the fatty acid synthesis was regulated by five lncRNAs that were upregulated under both long-term and short-term drought treatments. The key genes, including of SHM1, GOX2, and GS2, in the methylglyoxal synthesis pathway were all regulated by a number of down-regulated lncRNAs in OE lines, thereby reducing the content of such harmful compounds produced under stress in plants. This study identified a series of lncRNAs related to the pathways that affect photosynthesis, chlorophyll synthesis, fatty acid synthesis, degradation, and other important effects on drought resistance and oil content. The present study provided a series of lncRNAs for further improvement of crop varieties, especially drought resistant and oil content traits.

The Rhododendron Plant Genome Database (RPGD): a comprehensive online omics database for Rhododendron

Background

The genus Rhododendron L. has been widely cultivated for hundreds of years around the world. Members of this genus are known for great ornamental and medicinal value. Owing to advances in sequencing technology, genomes and transcriptomes of members of the Rhododendron genus have been sequenced and published by various laboratories. With increasing amounts of omics data available, a centralized platform is necessary for effective storage, analysis, and integration of these large-scale datasets to ensure consistency, independence, and maintainability.

Results

Here, we report our development of the Rhododendron Plant Genome Database (RPGD; http://bioinfor.kib.ac.cn/RPGD/), which represents the first comprehensive database of Rhododendron genomics information. It includes large amounts of omics data, including genome sequence assemblies for R. delavayi, R. williamsianum, and R. simsii, gene expression profiles derived from public RNA-Seq data, functional annotations, gene families, transcription factor identification, gene homology, simple sequence repeats, and chloroplast genome. Additionally, many useful tools, including BLAST, JBrowse, Orthologous Groups, Genome Synteny Browser, Flanking Sequence Finder, Expression Heatmap, and Batch Download were integrated into the platform.

Conclusions

RPGD is designed to be a comprehensive and helpful platform for all Rhododendron researchers. Believe that RPGD will be an indispensable hub for Rhododendron studies.

Full-length SMRT transcriptome sequencing and microsatellite characterization in Paulownia catalpifolia

Paulownia catalpifolia is an important, fast-growing timber species known for its high density, color and texture. However, few transcriptomic and genetic studies have been conducted in P. catalpifolia. In this study, single-molecule real-time sequencing technology was applied to obtain the full-length transcriptome of P. catalpifolia leaves treated with varying degrees of drought stress. The sequencing data were then used to search for microsatellites, or simple sequence repeats (SSRs). A total of 28.83 Gb data were generated, 25,969 high-quality (HQ) transcripts with an average length of 1624 bp were acquired after removing the redundant reads, and 25,602 HQ transcripts (98.59%) were annotated using public databases. Among the HQ transcripts, 16,722 intact coding sequences, 149 long non-coding RNAs and 179 alternative splicing events were predicted, respectively. A total of 7367 SSR loci were distributed throughout 6293 HQ transcripts, of which 763 complex SSRs and 6604 complete SSRs. The SSR appearance frequency was 28.37%, and the average distribution distance was 5.59 kb. Among the 6604 complete SSR loci, 1-3 nucleotide repeats were dominant, occupying 97.85% of the total SSR loci, of which mono-, di- and tri-nucleotide repeats were 44.68%, 33.86% and 19.31%, respectively. We detected 112 repeat motifs, of which A/T (42.64%), AG/CT (12.22%), GA/TC (9.63%), GAA/TTC (1.57%) and CCA/TGG (1.54%) were most common in mono-, di- and tri-nucleotide repeats, respectively. The length of the repeat SSR motifs was 10-88 bp, and 4997 (75.67%) were ≤ 20 bp. This study provides a novel full-length transcriptome reference for P. catalpifolia and will facilitate the identification of germplasm resources and breeding of new drought-resistant P. catalpifolia varieties.

Impact of different numbers of microsatellite markers on population genetic results using SLAF-seq data for Rhododendron species

Microsatellites (simple sequence repeats, SSRs) are co-dominant nuclear markers that are widely used in population genetic studies. Population genetic parameters from different studies might be significantly influenced by differences in marker number. In our study, 265 sequences with polymorphic microsatellites were obtained from SLAF-seq data. Then, subpopulations containing different numbers (5, 6, 7,…, 15, 20, 25, 30, 35, 40) of markers were genotyped 10 times to investigate the impact of marker numbers on population genetic diversity results. Our results show that genotyping with less than 11 or 12 microsatellite markers lead to significant deviations in the population genetic diversity or genetic structure results. In order to provide markers for population genetic and conservation studies for Rhododendron, 26 SSR primers were designed and validated in three species.

Single-Molecule Long-Read Sequencing of Purslane (Portulaca oleracea) and Differential Gene Expression Related with Biosynthesis of Unsaturated Fatty Acids

This study aimed to obtain the full-length transcriptome of purslane (Portulaca oleracea); assorted plant samples were used for single-molecule real-time (SMRT) sequencing. Based on SMRT, functional annotation of transcripts, transcript factors (TFs) analysis, simple sequence repeat analysis and long non-coding RNAs (LncRNAs) prediction were accomplished. Total 15.33-GB reads were produced; with 9,350,222 subreads and the average length of subreads, 1640 bp was counted. With 99.99% accuracy, after clustering, 132,536 transcripts and 78,559 genes were detected. All unique SMART transcripts were annotated in seven functional databases. 4180 TFs (including transcript regulators) and 7289 LncRNAs were predicted. The results of RNA-seq were confirmed with qRT–PCR analysis. Illumina sequencing of leaves and roots of two purslane genotypes was carried out. Amounts of differential expression genes and related KEGG pathways were found. The expression profiles of related genes in the biosynthesis of unsaturated fatty acids pathway in leaves and roots of two genotypes of purslane were analyzed. Differential expression of genes in this pathway built the foundation of ω-3 fatty acid accumulation in different organs and genotypes of purslane. The aforementioned results provide sequence information and may be a valuable resource for whole-genome sequencing of purslane in the future.