Transcriptional profiling of long noncoding RNAs associated with leaf-color mutation in Ginkgo biloba L

Effects of different light intensity on leaf color changes in a Chinese cabbage yellow cotyledon mutant

Leaf color is one of the most important phenotypic features in horticultural crops and directly related to the contents of photosynthetic pigments. Most leaf color mutants are determined by the altered chlorophyll or carotenoid, which can be affected by light quality and intensity. Our previous study obtained a Chinese cabbage yellow cotyledon mutant that exhibited obvious yellow phenotypes in the cotyledons and the new leaves. However, the underlying mechanisms in the formation of yellow cotyledons and leaves remain unclear. In this study, the Chinese cabbage yellow cotyledon mutant 19YC-2 exhibited obvious difference in leaf color and abnormal chloroplast ultrastructure compared to the normal green cotyledon line 19GC-2. Remarkably, low-intensity light treatment caused turn-green leaves and a significant decrease in carotenoid content in 19YC-2. RNA-seq analysis revealed that the pathways of photosynthesis antenna proteins and carotenoid biosynthesis were significantly enriched during the process of leaf color changes, and many differentially expressed genes related to the two pathways were identified to respond to different light intensities. Remarkably, BrPDS and BrLCYE genes related to carotenoid biosynthesis showed significantly higher expression in 19YC-2 than that in 19GC-2, which was positively related to the higher carotenoid content in 19YC-2. In addition, several differentially expressed transcription factors were also identified and highly correlated to the changes in carotenoid content, suggesting that they may participate in the regulatory pathway of carotenoid biosynthesis. These findings provide insights into the molecular mechanisms of leaf color changes in yellow cotyledon mutant 19YC-2 of Chinese cabbage.

Molecular insights into the mechanisms of a leaf color mutant in Anoectochilus roxburghii by gene mapping and transcriptome profiling based on PacBio Sequel II

Plants with partial or complete loss of chlorophylls and other pigments are frequently occurring in nature but not commonly found. In the present study, we characterize a leaf color mutant 'arly01' with an albino stripe in the middle of the leaf, which is an uncommon ornamental trait in Anoectochilus roxburghii. The albino "mutant" middle portion and green "normal" leaf parts were observed by transmission electron microscopy (TEM), and their pigment contents were determined. The mutant portion exhibited underdevelopment of plastids and had reduced chlorophyll and other pigment (carotenoid, anthocyanin, and flavonoid) content compared to the normal portion. Meanwhile, comparative transcript analysis and metabolic pathways mapping showed that a total of 599 differentially expressed genes were mapped to 78 KEGG pathways, most of which were down-regulated in the mutant portion. The five most affected metabolic pathways were determined to be oxidative phosphorylation, photosynthesis system, carbon fixation & starch and sucrose metabolism, porphyrin and chlorophyll metabolism, and flavonoid biosynthesis. Our findings suggested that the mutant 'arly01' was a partial albinism of A. roxburghii, characterized by the underdevelopment of chloroplasts, low contents of photosynthetic and other color pigments, and a number of down-regulated genes and metabolites. With the emergence of ornamental A. roxburghii in southern China, 'arly01' could become a popular cultivar due to its unique aesthetics.

Construction and analysis of the tapping panel dryness-related lncRNA/circRNA-miRNA-mRNA ceRNA network in latex of Hevea brasiliensis

Tapping panel dryness (TPD) results in a severe reduction in latex yield in Hevea brasiliensis. However, the molecular regulatory mechanisms of TPD occurrence are still largely unclear. In this study, whole-transcriptome sequencing was carried out on latex from TPD and healthy trees. In total, 7078 long noncoding RNAs (lncRNAs), 3077 circular RNAs (circRNAs), 4956 miRNAs, and 25041 mRNAs were identified in latex, among which 435 lncRNAs, 68 circRNAs, 320 miRNAs, and 1574 mRNAs were differentially expressed in the latex of TPD trees. GO and KEGG analyses indicated that plant hormone signal transduction, MAPK signaling pathway, and ubiquitin-mediated proteolysis were the key pathways associated with TPD onset. Phytohormone profiling revealed significant changes in the contents of 28 hormonal compounds, among which ACC, ABA, IAA, GA, and JA contents were increased, while SA content was reduced in TPD latex, suggesting that hormone homeostasis is disrupted in TPD trees. Furthermore, we constructed a TPD-related competitive endogenous RNA (ceRNA) regulatory network of lncRNA/circRNA-miRNA-mRNA with 561 edges and 434 nodes (188 lncRNAs, 5 circRNAs, 191 miRNAs, and 50 mRNAs) and identified two hub lncRNAs (MSTRG.11908.1 and MSTRG.8791.1) and four hub miRNAs (hbr-miR156, miR156-x, miRf10477-y, and novel-m0452-3p). Notably, the lncRNA-miR156/157-SPL module containing three hubs probably plays a crucial role in TPD onset. The expression of network hubs and the lncRNA-miR156/157-SPL module were further validated by qRT-PCR. Our results reveal the TPD-associated ceRNA regulatory network of lncRNA/circRNA-miRNA-mRNA in latex and lay a foundation for further investigation of molecular regulatory mechanisms for TPD onset in H. brasiliensis.

A long-read and short-read transcriptomics approach provides the first high-quality reference transcriptome and genome annotation for Pseudotsuga menziesii (Douglas-fir)

Douglas-fir (Pseudotsuga menziesii) is native to western North America. It grows in a wide range of environmental conditions and is an important timber tree. Although there are several studies on the gene expression responses of Douglas-fir to abiotic cues, the absence of high-quality transcriptome and genome data is a barrier to further investigation. Like for most conifers, the available transcriptome and genome reference dataset for Douglas-fir remains fragmented and requires refinement. We aimed to generate a highly accurate, and complete reference transcriptome and genome annotation. We deep-sequenced the transcriptome of Douglas-fir needles from seedlings that were grown under nonstress control conditions or a combination of heat and drought stress conditions using long-read (LR) and short-read (SR) sequencing platforms. We used 2 computational approaches, namely de novo and genome-guided LR transcriptome assembly. Using the LR de novo assembly, we identified 1.3X more high-quality transcripts, 1.85X more "complete" genes, and 2.7X more functionally annotated genes compared to the genome-guided assembly approach. We predicted 666 long noncoding RNAs and 12,778 unique protein-coding transcripts including 2,016 putative transcription factors. We leveraged the LR de novo assembled transcriptome with paired-end SR and a published single-end SR transcriptome to generate an improved genome annotation. This was conducted with BRAKER2 and refined based on functional annotation, repetitive content, and transcriptome alignment. This high-quality genome annotation has 51,419 unique gene models derived from 322,631 initial predictions. Overall, our informatics approach provides a new reference Douglas-fir transcriptome assembly and genome annotation with considerably improved completeness and functional annotation.

Transcriptional profiling of long non-coding RNAs regulating fruit cracking in Punica granatum L. under bagging

Fruit cracking tremendously damages the appearance of fruit, easily leads to pathogen invasion, greatly reduces the marketability and causes immense economic losses. The pivotal role of long non-coding RNAs (lncRNAs) in diverse biological processes has been confirmed, while the roles of lncRNAs underlying fruit cracking remain poorly understood. In this study, the incidence of fruit cracking was 7.26% under the bagging treatment, the control group was 38.11%, indicating that bagging considerably diminished the fruit cracking rate. LncRNA libraries for fruit cracking (FC), fruit non-cracking (FNC) and fruit non-cracking under bagging (FB) in pomegranate (Punica granatum L.) were performed and analysed via high-throughput transcriptome sequencing. A total of 3194 lncRNAs were obtained with a total length of 4898846 nt and an average length of 1533.77 nt in pomegranate. We identified 42 differentially expressed lncRNAs (DELs) and 137 differentially expressed mRNAs (DEGs) in FC vs FNC and 35 DELs and 160 DEGs in FB vs FC that formed co-expression networks respectively, suggesting that there are involved in phytohormone signaling pathway, lignin catabolic process, lipid transport/binding, cutin biosynthetic process and cell wall organization. We also found that 18 cis-acting DELs regulated 18 target genes, and 10 trans-acting DELs regulated 24 target genes in FC vs FNC, 23 DELs regulate 23 target genes for the cis-acting lncRNAs and 12 DELs regulated 36 target genes in FB vs FC, which provides an understanding for the regulation of the fruit cracking. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis results demonstrated that DELs participated in calcium ion binding, glycerophospholipid metabolism, flavonoid biosynthetic process, cell wall biogenesis, xyloglucan metabolic process, hormone signal transduction and starch and sucrose metabolism. Our findings provide new insights into the roles of lncRNAs in regulating the fruit cracking and lay the foundation for further improvement of pomegranate quality.

Dynamic transcriptome and network-based analysis of yellow leaf mutant Ginkgo biloba

BACKGROUND

Golden leaf in autumn is a prominent feature of deciduous tree species like Ginkgo biloba L., a landscape tree widely cultivated worldwide. However, little was known about the molecular mechanisms of leaf yellowing, especially its dynamic regulatory network. Here, we performed a suite of comparative physiological and dynamic transcriptional analyses on the golden-leaf cultivar and the wild type (WT) ginkgo to investigate the underlying mechanisms of leaf yellowing across different seasons.

RESULTS

In the present study, we used the natural bud mutant cultivar with yellow leaves "Wannianjin" (YL) as materials. Physiological analysis revealed that higher ratios of chlorophyll a to chlorophyll b and carotenoid to chlorophyll b caused the leaf yellowing of YL. On the other hand, dynamic transcriptome analyses showed that genes related to chlorophyll metabolism played key a role in leaf coloration. Genes encoding non-yellow coloring 1 (NYC1), NYC1-like (NOL), and chlorophyllase (CLH) involved in the degradation of chlorophyll were up-regulated in spring. At the summer stage, down-regulated HEMA encoding glutamyl-tRNA reductase functioned in chlorophyll biosynthesis, while CLH involved in chlorophyll degradation was up-regulated, causing a lower chlorophyll accumulation. In carotenoid metabolism, genes encoding zeaxanthin epoxidase (ZEP) and 9-cis-epoxy carotenoid dioxygenase (NCED) showed significantly different expression levels in the WT and YL. Moreover, the weighted gene co-expression network analysis (WGCNA) suggested that the most associated transcriptional factor, which belongs to the AP2/ERF-ERF family, was engaged in regulating pigment metabolism. Furthermore, quantitative experiments validated the above results.

CONCLUSIONS

By comparing the golden-leaf cultivar and the wide type of ginkgo across three seasons, this study not only confirm the vital role of chlorophyll in leaf coloration of YL but also provided new insights into the seasonal transcriptome landscape and co-expression network. Our novel results pinpoint candidate genes for further wet-bench experiments in tree species.

Population transcriptomic analysis identifies the comprehensive lncRNAs landscape of spike in wheat (Triticum aestivum L.)

BACKGROUND

Long noncoding RNAs (lncRNAs) are emerging as the important regulators involving in growth and development as well as stress response in plants. However, current lncRNA studies were mainly performed at the individual level and the significance of it is not well understood in wheat.

RESULTS

In this study, the lncRNA landscape of wheat spike was characterized through analysing a total of 186 spike RNA-seq datasets from 93 wheat genotypes. A total of 35,913 lncRNAs as well as 1,619 lncRNA-mRNA pairs comprised of 443 lncRNAs and 464 mRNAs were obtained. Compared to coding genes, these lncRNAs displayed rather low conservation among wheat and other gramineous species. Based on re-sequencing data, the genetic variations of these lncRNA were investigated and obvious genetic bottleneck were found on them during wheat domestication process. Furthermore, 122 lncRNAs were found to act as ceRNA to regulate endogenous competition. Finally, association and co-localization analysis of the candidate lncRNA-mRNA pairs identified 170 lncRNAs and 167 target mRNAs significantly associated with spike-related traits, including lncRNA.127690.1/TraesCS2A02G518500.1 (PMEI) and lncRNA.104854.1/TraesCS6A02G050300.1 (ATG5) associated with heading date and spike length, respectively.

CONCLUSIONS

This study reported the lncRNA landscape of wheat spike through the population transcriptome analysis, which not only contribute to better understand the wheat evolution from the perspective of lncRNA, but also lay the foundation for revealing roles of lncRNA playing in spike development.

Genome-wide identification of lncRNAs, miRNAs, mRNAs and their regulatory networks involved in tapping panel dryness in rubber tree (Hevea brasiliensis)

Noncoding RNAs (ncRNAs) play pivotal roles in various biological processes in plants. However, the role of ncRNAs in tapping panel dryness (TPD) of rubber tree (Hevea brasiliensis Muell. Arg.) is largely unknown. Here, the whole transcriptome analyses of bark tissues from healthy and TPD trees were performed to identify differentially expressed long ncRNAs (DELs), microRNAs/miRNAs (DEMs), genes (DEGs) and their regulatory networks involved in TPD. A total of 263 DELs, 174 DEMs and 1574 DEGs were identified in the bark of TPD tree compared with that of healthy tree. Kyoto Encyclopedia of Genes and Genomes analysis revealed that most of the DEGs and targets of DELs and DEMs were mainly enriched in metabolic pathways, biosynthesis of secondary metabolites and plant hormone signal transduction. Additionally, the majority of DEGs and DELs related to rubber biosynthesis were downregulated in TPD trees. Furthermore, 98 DEGs and 44 DELs were targeted by 54 DEMs, 190 DEGs were identified as putative targets of 56 DELs, and 2 and 44 DELs were predicted as precursors and endogenous target mimics of 2 and 6 DEMs, respectively. Based on these, the DEL-DEM-DEG regulatory network involved in TPD was constructed, and 13 hub DELs, 3 hub DEMs and 2 hub DEGs were identified. The results provide novel insights into the regulatory roles of ncRNAs underlying TPD and lay a foundation for future functional characterization of long ncRNAs, miRNAs and genes involved in TPD in rubber tree.

Differential Expression of lncRNAs and miRNAs Between Self-Rooting Juvenile and Donor Clones Unveils Novel Insight Into the Molecular Regulation of Rubber Biosynthesis in Hevea brasiliensis

The rubber tree (Hevea brasiliensis Muell. Arg.) is a tropical tree species that produce natural rubber. Self-rooted juvenile clones (SRJCs) are novel rubber tree planting materials developed through primary somatic embryogenesis. SRJCs have a higher rubber yield compared with donor clones (DCs). The molecular basis underlying increased rubber yield in SRJCs remains largely unknown. Here, the latex from SRJCs and DCs were collected for strand-specific and small RNA-seq methods. A total of 196 differentially expressed long noncoding RNAs (DELs), and 11 differentially expressed microRNAs were identified in latex between SRJCs and DCs. Targeted genes of DELs were markedly enriched for various biological pathways related to plant hormone signal transduction, photosynthesis, glutathione metabolism, and amino acids biosynthesis. DELs probably acted as cis-acting regulation was calculated, and these DELs relevant to potentially regulate rubber biosynthesis, reactive oxygen species metabolism, and epigenetic modification. Furthermore, the DELs acting as microRNA targets were studied. The interaction of microRNA and DELs might involve in the regulation of natural rubber biosynthesis.

Bidirectional lncRNA Transfer between Cuscuta Parasites and Their Host Plant

Dodder species (Cuscuta spp.) are holoparasites that have extensive material exchange with their host plants through vascular connections. Recent studies on cross-species transfer have provided breakthrough insights, but little is known about the interaction mechanisms of the inter-plant mobile substances in parasitic systems. We sequenced the transcriptomes of dodder growing on soybean hosts to characterize the long non-coding RNA (lncRNA) transfer between the two species, and found that lncRNAs can move in high numbers (365 dodder lncRNAs and 14 soybean lncRNAs) in a bidirectional manner. Reverse transcription-polymerase chain reaction further confirmed that individual lncRNAs were trafficked in the dodder-soybean parasitic system. To reveal the potential functions of mobile transcripts, the Gene Ontology terms of mobile lncRNA target genes were predicted, and mobile dodder target genes were found to be mainly enriched in "metabolic process", "catalytic activity", "signaling", and "response to stimulus" categories, whereas mobile soybean target genes were enriched in organelle-related categories, indicating that specific mobile lncRNAs may be important in regulating dodder parasitism. Our findings reveal that lncRNAs are transferred between dodder and its host soybean plants, which may act as critical regulators to coordinate the host-dodder interaction at the whole parasitic level.

Fine dissection of limber pine resistance to Cronartium ribicola using targeted sequencing of the NLR family

BACKGROUND

Proteins with nucleotide binding site (NBS) and leucine-rich repeat (LRR) domains (NLR) make up one of most important resistance (R) families for plants to resist attacks from various pathogens and pests. The available transcriptomes of limber pine (Pinus flexilis) allow us to characterize NLR genes and related resistance gene analogs (RGAs) in host resistance against Cronartium ribicola, the causal fungal pathogen of white pine blister rust (WPBR) on five-needle pines throughout the world. We previously mapped a limber pine major gene locus (Cr4) that confers complete resistance to C. ribicola on the Pinus consensus linkage group 8 (LG-8). However, genetic distribution of NLR genes as well as their divergence between resistant and susceptible alleles are still unknown.

RESULTS

To identify NLR genes at the Cr4 locus, the present study re-sequenced a total of 480 RGAs using targeted sequencing in a Cr4-segregated seed family. Following a call of single nucleotide polymorphisms (SNPs) and genetic mapping, a total of 541 SNPs from 155 genes were mapped across 12 LGs. Three putative NLR genes were newly mapped in the Cr4 region, including one that co-segregated with Cr4. The tight linkage of NLRs with Cr4-controlled phenotypes was further confirmed by bulked segregation analysis (BSA) using extreme-phenotype genome-wide association study (XP-GWAS) for significance test. Local tandem duplication in the Cr4 region was further supported by syntenic analysis using the sugar pine genome sequence. Significant gene divergences have been observed in the NLR family, revealing that diversifying selection pressures are relatively higher in local duplicated genes. Most genes showed similar expression patterns at low levels, but some were affected by genetic background related to disease resistance. Evidence from fine genetic dissection, evolutionary analysis, and expression profiling suggests that two NLR genes are the most promising candidates for Cr4 against WPBR.

CONCLUSION

This study provides fundamental insights into genetic architecture of the Cr4 locus as well as a set of NLR variants for marker-assisted selection in limber pine breeding. Novel NLR genes were identified at the Cr4 locus and the Cr4 candidates will aid deployment of this R gene in combination with other major/minor genes in the limber pine breeding program.

Regulatory role of non-coding RNA in ginseng rusty root symptom tissue

Ginseng rusty root symptom (GRS) is one of the primary diseases of ginseng. It leads to a severe decline in the quality of ginseng and significantly affects the ginseng industry. The regulatory mechanism of non-coding RNA (ncRNA) remains unclear in the course of disease. This study explored the long ncRNAs (lncRNAs), circular RNAs (circRNAs), and microRNAs (miRNAs) in GRS tissues and healthy ginseng (HG) tissues and performed functional enrichment analysis of the screened differentially expressed ncRNAs. Considering the predictive and regulatory effects of ncRNAs on mRNAs, we integrated ncRNA and mRNA data to analyze and construct relevant regulatory networks. A total of 17,645 lncRNAs, 245 circRNAs, and 299 miRNAs were obtained from HG and GRS samples, and the obtained ncRNAs were characterized, including the classification of lncRNAs, length and distribution of circRNA, and the length and family affiliations of miRNAs. In the analysis of differentially expressed ncRNA target genes, we found that lncRNAs may be involved in the homeostatic process of ginseng tissues and that lncRNAs, circRNAs, and miRNAs are involved in fatty acid-related regulation, suggesting that alterations in fatty acid-related pathways may play a key role in GRS. Besides, differentially expressed ncRNAs play an essential role in regulating transcriptional translation processes, primary metabolism such as starch and sucrose, and secondary metabolism such as alkaloids in ginseng tissues. Finally, we integrated the correlations between ncRNAs and mRNAs, constructed corresponding interaction networks, and identified ncRNAs that may play critical roles in GRS. These results provide a basis for revealing GRS's molecular mechanism and enrich our understanding of ncRNAs in ginseng.

Comparative Analyses of Full-Length Transcriptomes Reveal Gnetum luofuense Stem Developmental Dynamics

Genus Gnetum, of which the majority species are pantropical liana, have broad industrial uses including for string, nets, and paper production. Although numerous studies have investigated anatomical structures during stem development, the underlying molecular mechanisms that regulate this developmental trajectory in Gnetum species remain poorly understood. A total of 12 full-length transcriptomes were generated from four stem developmental stages of an arborescent representative of this genus, Gnetum luofuense, using Oxford Nanopore Technologies. The results of this analysis reveal a total of 24,151 alternative splicing (AS) and 134,391 alternative polyadenylation events. A remarkably dynamic pattern of AS events, especially in the case of intron retentions, was found across the four developmental stages while no dynamic pattern was found among transcript numbers with varied poly(A) sites. A total of 728 long non-coding RNAs were also detected; the number of cis-regulated target genes dramatically increased while no changes were found among trans-regulated target genes. In addition, a K-means clustering analysis of all full-length transcripts revealed that primary growth is associated with carbohydrate metabolism and fungi defense, while secondary growth is closely linked with photosynthesis, nitrogen transportation, and leaf ontogenesis. The use of weighted gene co-expression network analysis as well as differentially expressed transcripts reveals that bHLH, GRF, and MYB-related transcription factors are involved in primary growth, while AP2/ERF, MYB, NAC, PLAZ, and bZIP participate in G. luofuense stem secondary growth. The results of this study provide further evidence that Nanopore sequencing technology provides a cost-effective method for generating full-length transcriptome data as well as for investigating seed plant organ development.

A full-length transcriptome and gene expression analysis reveal genes and molecular elements expressed during seed development in Gnetum luofuense

BACKGROUND

Gnetum is an economically important tropical and subtropical gymnosperm genus with various dietary, industrial and medicinal uses. Many carbohydrates, proteins and fibers accumulate during the ripening of Gnetum seeds. However, the molecular mechanisms related to this process remain unknown.

RESULTS

We therefore assembled a full-length transcriptome from immature and mature G. luofuense seeds 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 of G. luofuense.

CONCLUSIONS

These findings provide a valuable molecular resource for understanding seed development of gymnosperms.