Genome-wide identification and functional characterization of natural antisense transcripts in Salvia miltiorrhiza

Genome-Wide Identification and Characterization of miRNAs and Natural Antisense Transcripts Show the Complexity of Gene Regulatory Networks for Secondary Metabolism in Aristolochia contorta

Aristolochia contorta Bunge is an academically and medicinally important plant species. It belongs to the magnoliids, with an uncertain phylogenetic position, and is one of the few plant species lacking a whole-genome duplication (WGD) event after the angiosperm-wide WGD. A. contorta has been an important traditional Chinese medicine material. Since it contains aristolochic acids (AAs), chemical compounds with nephrotoxity and carcinogenicity, the utilization of this plant has attracted widespread attention. Great efforts are being made to increase its bioactive compounds and reduce or completely remove toxic compounds. MicroRNAs (miRNAs) and natural antisense transcripts (NATs) are two classes of regulators potentially involved in metabolism regulation. Here, we report the identification and characterization of 223 miRNAs and 363 miRNA targets. The identified miRNAs include 51 known miRNAs belonging to 20 families and 172 novel miRNAs belonging to 107 families. A negative correlation between the expression of miRNAs and their targets was observed. In addition, we identified 441 A. contorta NATs and 560 NAT-sense transcript (ST) pairs, of which 12 NATs were targets of 13 miRNAs, forming 18 miRNA-NAT-ST modules. Various miRNAs and NATs potentially regulated secondary metabolism through the modes of miRNA-target gene-enzyme genes, NAT-STs, and NAT-miRNA-target gene-enzyme genes, suggesting the complexity of gene regulatory networks in A. contorta. The results lay a solid foundation for further manipulating the production of its bioactive and toxic compounds.

Genome-wide view and characterization of natural antisense transcripts in Cannabis Sativa L

Natural Antisense Transcripts (NATs) are a kind of complex regulatory RNAs that play crucial roles in gene expression and regulation. However, the NATs in Cannabis Sativa L., a widely economic and medicinal plant rich in cannabinoids remain unknown. In this study, we comprehensively predicted C. sativa NATs genome-wide using strand-specific RNA sequencing (ssRNA-Seq) data, and validated the expression profiles by strand-specific quantitative reverse transcription PCR (ssRT-qPCR). Consequently, a total of 307 NATs were predicted in C. sativa, including 104 cis- and 203 trans- NATs. Functional enrichment analysis demonstrated the potential involvement of the C. sativa NATs in DNA polymerase activity, RNA-DNA hybrid ribonuclease activity, and nucleic acid binding. Finally, 18 cis- and 376 trans- NAT-ST pairs were predicted to produce 621 cis- and 5,679 trans- small interfering RNA (nat-siRNAs), respectively. These nat-siRNAs were potentially involved in the biosynthesis of cannabinoids and cellulose. All these results will shed light on the regulation of NATs and nat-siRNAs in C. sativa.

Global identification of natural antisense transcripts in Gossypium hirsutum and Gossypium barbadense under chilling stress

Natural antisense transcripts (NATs) in model plants have been recognized as important regulators of gene expression under abiotic stresses. However, the functional roles of NATs in crops under low temperature are still unclear. Here, we identified 815 and 689 NATs from leaves of Gossypium hirsutum and G. barbadense under chilling stress. Among those, 224 NATs were identified as interspecific homologs between the two species. The correlation coefficients for expression of NATs and their cognate sense genes (CSG) were 0.43 and 0.37 in G. hirsutum and G. barbadense, respectively. Furthermore, expression of interspecific NATs and CSGs alike was highly consistent under chilling stress with correlation coefficients of 0.90-0.91. Four cold-associated NATs were selected for functional validation using virus-induced gene silencing (VIGS). Our results suggest that CAN1 engage in the molecular regulation of chilling stress by regulating SnRK2.8 expression. This highly conserved NAT have valuable potential for applications in breeding cold-tolerant cotton.

Integrated analysis of lncRNAs, mRNAs, and TFs to identify network modules underlying diterpenoid biosynthesis in Salvia miltiorrhiza

Long non-coding RNAs (lncRNAs) are transcripts of more than 200 nucleotides (nt) in length, with minimal or no protein-coding capacity. Increasing evidence indicates that lncRNAs play important roles in the regulation of gene expression including in the biosynthesis of secondary metabolites. Salvia miltiorrhiza Bunge is an important medicinal plant in China. Diterpenoid tanshinones are one of the main active components of S. miltiorrhiza. To better understand the role of lncRNAs in regulating diterpenoid biosynthesis in S. miltiorrhiza, we integrated analysis of lncRNAs, mRNAs, and transcription factors (TFs) to identify network modules underlying diterpenoid biosynthesis based on transcriptomic data. In transcriptomic data, we obtained 6,651 candidate lncRNAs, 46 diterpenoid biosynthetic pathway genes, and 11 TFs involved in diterpenoid biosynthesis. Combining the co-expression and genomic location analysis, we obtained 23 candidate lncRNA-mRNA/TF pairs that were both co-expressed and co-located. To further observe the expression patterns of these 23 candidate gene pairs, we analyzed the time-series expression of S. miltiorrhiza induced by methyl jasmonate (MeJA). The results showed that 19 genes were differentially expressed at least a time-point, and four lncRNAs, two mRNAs, and two TFs formed three lncRNA-mRNA and/or TF network modules. This study revealed the relationship among lncRNAs, mRNAs, and TFs and provided new insight into the regulation of the biosynthetic pathway of S. miltiorrhiza diterpenoids.

Genome-wide identification of antisense lncRNAs and their association with susceptibility to Flavobacterium psychrophilum in rainbow trout

Eukaryotic genomes encode long noncoding natural antisense transcripts (lncNATs) that have been increasingly recognized as regulatory members of gene expression. Recently, we identified a few antisense transcripts correlating in expression with immune-related genes. However, a systematic genome-wide analysis of lncNATs in rainbow trout is lacking. This study used 134 RNA-Seq datasets from five different projects to identify antisense transcripts. A total of 13,503 lncNATs were identified genome-wide. About 75% of lncNATs showed multiple exons compared to 36.5% of the intergenic lncRNAs. RNA-Seq datasets from resistant, control, and susceptible rainbow trout genetic lines with significant differences in survival rate following Flavobacterium psychrophilum (Fp) infection were analyzed to investigate the potential role of the lncNATs during infection. Twenty-four pairwise comparisons between the different genetic lines, infectious status, and time points revealed 581 differentially expressed (DE) lncNATs and 179 differentially used exons (DUEs). Most of the DE lncNATs strongly and positively correlated in expression with their corresponding sense transcripts across 24 RNA-Seq datasets. LncNATs complementary to genes related to immunity, muscle contraction, proteolysis, and iron/heme metabolism were DE following infection. LncNATs complementary to hemolysis-related genes were DE in the resistant fish compared to susceptible fish on day 5 post-infection, suggesting enhanced clearance of free hemoglobin (Hb) and heme and increased erythropoiesis. LncNATs complementary to hepcidin, a master negative regulator of the plasma iron concentration, were the most downregulated lncNATs on day 5 of bacterial infection in the resistant fish. Ninety-four DE lncNAT, including five complementary to hepcidin, are located within 26 QTL regions previously identified in association with bacterial cold water disease (BCWD) in rainbow trout. Collectively, lncNATs are involved in the molecular architecture of fish immunity and should be further investigated for potential applications in genomic selection and genetic manipulation in aquaculture.

Genome-Wide Prediction and Analysis of Oryza Species NRP Genes in Rice Blast Resistance

Members of the N-rich proteins (NRPs) gene family play important roles in the plant endoplasmic reticulum stress in response, which can be triggered by plant pathogens’ infection. Previous studies of the NRP gene family have been limited to only a few plants, such as soybean and Arabidopsis thaliana. Thus, their evolutionary characteristics in the Oryza species and biological functions in rice defense against the pathogenic fungus Magnaporthe oryzae have remained unexplored. In the present study, we demonstrated that the NRP genes family may have originated in the early stages of plant evolution, and that they have been strongly conserved during the evolution of the Oryza species. Domain organization of NRPs was found to be highly conserved within but not between subgroups. OsNRP1, an NRP gene in the Oryza sativa japonica group, was specifically up-regulated during the early stages of rice-M. oryzae interactions-inhibited M. oryzae infection. Predicted protein-protein interaction networks and transcription-factor binding sites revealed a candidate interactor, bZIP50, which may be involved in OsNRP1-mediated rice resistance against M. oryzae infection. Taken together, our results established a basis for future studies of the NRP gene family and provided molecular insights into rice immune responses to M. oryzae.

Genome-wide identification and multiple abiotic stress transcript profiling of potassium transport gene homologs in Sorghum bicolor

Potassium (K⁺) is the most abundant cation that plays a crucial role in various cellular processes in plants. Plants have developed an efficient mechanism for the acquisition of K⁺ when grown in K⁺ deficient or saline soils. A total of 47 K⁺ transport gene homologs (27 HAKs, 4 HKTs, 2 KEAs, 9 AKTs, 2 KATs, 2 TPCs, and 1 VDPC) have been identified in Sorghum bicolor. Of 47 homologs, 33 were identified as K⁺ transporters and the remaining 14 as K⁺ channels. Chromosome 2 has been found as the hotspot of K⁺ transporters with 9 genes. Phylogenetic analysis revealed the conservation of sorghum K⁺ transport genes akin to Oryza sativa. Analysis of regulatory elements indicates the key roles that K⁺ transport genes play under different biotic and abiotic stress conditions. Digital expression data of different developmental stages disclosed that expressions were higher in milk, flowering, and tillering stages. Expression levels of the genes SbHAK27 and SbKEA2 were higher during milk, SbHAK17, SbHAK11, SbHAK18, and SbHAK7 during flowering, SbHAK18, SbHAK10, and 23 other gene expressions were elevated during tillering inferring the important role that K⁺ transport genes play during plant growth and development. Differential transcript expression was observed in different tissues like root, stem, and leaf under abiotic stresses such as salt, drought, heat, and cold stresses. Collectively, the in-depth genome-wide analysis and differential transcript profiling of K⁺ transport genes elucidate their role in ion homeostasis and stress tolerance mechanisms.

Genome-Wide Identification of Eucalyptus Heat Shock Transcription Factor Family and Their Transcriptional Analysis under Salt and Temperature Stresses

Heat shock transcription factors (HSFs) activate heat shock protein gene expression by binding their promoters in response to heat stress and are considered to be pivotal transcription factors in plants. Eucalyptus is a superior source of fuel and commercial wood. During its growth, high temperature or other abiotic stresses could impact its defense capability and growth. Hsf genes have been cloned and sequenced in many plants, but rarely in Eucalyptus. In this study, we used bioinformatics methods to analyze and identify Eucalyptus Hsf genes, their chromosomal localization and structure. The phylogenetic relationship and conserved domains of their encoded proteins were further analyzed. A total of 36 Hsf genes were identified and authenticated from Eucalyptus, which were scattered across 11 chromosomes. They could be classified into three classes (A, B and C). Additionally, a large number of stress-related cis-regulatory elements were identified in the upstream promoter sequence of HSF, and cis-acting element analysis indicated that the expression of EgHsf may be regulated by plant growth and development, metabolism, hormones and stress responses. The expression profiles of five representative Hsf genes, EgHsf4, EgHsf9, EgHsf13, EgHsf24 and EgHsf32, under salt and temperature stresses were examined by qRT-PCR. The results show that the expression pattern of class B genes (EgHsf4, EgHsf24 and EgHsf32) was more tolerant to abiotic stresses than that of class A genes (EgHsf9 and EgHsf13). However, the expressions of all tested Hsf genes in six tissues were at different levels. Finally, we investigated the network of interplay between genes, and the results suggest that there may be synergistic effects between different Hsf genes in response to abiotic stresses. We conclude that the Hsf gene family played an important role in the growth and developmental processes of Eucalyptus and could be vital for maintaining cell homeostasis against external stresses. This study provides basic information on the members of the Hsf gene family in Eucalyptus and lays the foundation for the functional identification of related genes and the further investigation of their biological functions in plant stress regulation.

miR398b and AtC2GnT form a negative feedback loop to regulate Arabidopsis thaliana resistance against Phytophthora parasitica

Oomycetes are diploid eukaryotic microorganisms that seriously threaten sustainable crop production. MicroRNAs (miRNAs) and corresponding natural antisense transcripts (NATs) are important regulators of multiple biological processes. However, little is known about their roles in plant immunity against oomycete pathogens. In this study, we report the identification and functional characterization of miR398b and its cis-NAT, the core-2/I-branching beta-1,6-N-acetylglucosaminyltransferase gene (AtC2GnT), in plant immunity. Gain- and loss-of-function assays revealed that miR398b mediates Arabidopsis thaliana susceptibility to Phytophthora parasitica by targeting Cu/Zn-Superoxidase Dismutase1 (CSD1) and CSD2, leading to suppressed expression of CSD1 and CSD2 and decreased plant disease resistance. We further showed that AtC2GnT transcripts could inhibit the miR398b-CSDs module via inhibition of pri-miR398b expression, leading to elevated plant resistance to P. parasitica. Furthermore, quantitative reverse transcription PCR, RNA ligase-mediated 5'-amplification of cDNA ends (RLM-5' RACE), and transient expression assays indicated that miR398b suppresses the expression of AtC2GnT. We generated AtC2GnT-silenced A. thaliana plants by CRISPR/Cas9 or RNA interference methods, and the Nicotiana benthamiana NbC2GnT-silenced plants by virus-induced gene silencing. Pathogenicity assays showed that the C2GnT-silenced plants were more susceptible, while AtC2GnT-overexpressing plants exhibited elevated resistance to P. parasitica. AtC2GnT encodes a Golgi-localized protein, and transient expression of AtC2GnT enhanced N. benthamiana resistance to Phytophthora pathogens. Taken together, our results revealed a positive role of AtC2GnT and a negative regulatory loop formed by miR398b and AtC2GnT in regulating plant resistance to P. parasitica.

The PAP Gene Family in Tomato: Comprehensive Comparative Analysis, Phylogenetic Relationships and Expression Profiles

Purple acid phosphatase (PAP) plays a vital role in plant phosphate acquisition and utilization, as well as cell wall synthesis and redox reactions. In this study, comprehensive comparative analyses of PAP genes were carried out using the integration of phylogeny, chromosomal localization, intron/exon structural characteristics, and expression profiling. It was shown that the number of introns of the PAP genes, which were distributed unevenly on 12 chromosomes, ranged from 1 to 12. These findings pointed to the existence of complex structures. Phylogenetic analyses revealed that PAPs from tomato, rice, and Arabidopsis could be divided into three groups (Groups I, II, and III). It was assumed that the diversity of these PAP genes occurred before the monocot–dicot split. RNA-seq analysis revealed that most of the genes were expressed in all of the tissues analyzed, with the exception of SlPAP02, SlPAP11, and SlPAP14, which were not detected. It was also found that expression levels of most of the SlPAP gene family of members were changed under phosphorus stress conditions, suggesting potential functional diversification. The findings of this work will help us to achieve a better insight into the function of SlPAP genes in the future, as well as enhance our understanding of their evolutionary relationships in plants.

Comprehensive Analysis of Carotenoid Cleavage Dioxygenases Gene Family and Its Expression in Response to Abiotic Stress in Poplar

Carotenoid cleavage dioxygenases (CCDs) catalyzes the cleavage of various carotenoids into smaller apocarotenoids which are essential for plant growth and development and response to abiotic stresses. CCD family is divided into two subfamilies: 9-cis epoxycarotenoid dioxygenases (NCED) family and CCD family. A better knowledge of carotenoid biosynthesis and degradation could be useful for regulating carotenoid contents. Here, 23 CCD genes were identified from the Populus trichocarpa genome, and their characterizations and expression profiling were validated. The PtCCD members were divided into PtCCD and PtNCED subfamilies. The PtCCD family contained the PtCCD1, 4, 7, and 8 classes. The PtCCDs clustered in the same clade shared similar intron/exon structures and motif compositions and distributions. In addition, the tandem and segmental duplications resulted in the PtCCD gene expansion based on the collinearity analysis. An additional integrated collinearity analysis among poplar, Arabidopsis, rice, and willow revealed the gene pairs between poplar and willow more than that between poplar and rice. Identifying tissue-special expression patterns indicated that PtCCD genes display different expression patterns in leaves, stems, and roots. Abscisic acid (ABA) treatment and abiotic stress suggested that many PtCCD genes are responsive to osmotic stress regarding the comprehensive regulation networks. The genome-wide identification of PtCCD genes may provide the foundation for further exploring the putative regulation mechanism on osmotic stress and benefit poplar molecular breeding.

Genome-Wide Characterization and Expression Analysis of Pathogenesis-Related 1 (PR-1) Gene Family in Tea Plant (Camellia sinensis (L.) O. Kuntze) in Response to Blister-Blight Disease Stress

Pathogenesis-related 1 (PR-1) proteins, which are defense proteins in plant–pathogen interactions, play an important role in the resistance and defense of plants against diseases. Blister blight disease is caused by Exobasidium vexans Massee and a major leaf disease of tea plants (Camellia sinensis (L.) O. Kuntze). However, the systematic characterization and analysis of the PR-1 gene family in tea plants is still lacking, and the defense mechanism of this family remains unknown. In this study, 17 CsPR-1 genes were identified from the tea plant genome and classified into five groups based on their signal peptide, isoelectric point, and C-terminus extension. Most of the CsPR-1 proteins contained an N-terminal signal peptide and a conserved PR-1 like domain. CsPR-1 genes comprised multiple cis-acting elements and were closely related to the signal-transduction pathways involving TCA, NPR1, EDS16, BGL2, PR4, and HCHIB. These characteristics imply an important role of the genes in the defense of the tea plant. In addition, the RNA-seq data and real-time PCR analysis demonstrated that the CsPR-1-2, -4, -6, -7, -8, -9, -10, -14, -15, and -17 genes were significantly upregulated under tea blister-blight stress. This study could help to increase understanding of CsPR-1 genes and their defense mechanism in response to tea blister blight.

Transcriptome Analysis Reveals Candidate Genes Involved in Light-Induced Primordium Differentiation in Pleurotus eryngii

Pleurotus eryngii, a highly valued edible fungus, is one of the major commercially cultivated mushrooms in China. The development of P. eryngii, especially during the stage of primordium differentiation, is easily affected by light. However, the molecular mechanism underlying the response of primordium differentiation to light remains unknown. In the present study, primordium expression profiles under blue-light stimulation, red-light stimulation, and exposure to darkness were compared using high-throughput sequencing. A total of 16,321 differentially expressed genes (DEGs) were identified from three comparisons. GO enrichment analysis showed that a large number of DEGs were related to light stimulation and amino acid biosynthesis. KEGG analyses demonstrated that the MAPK signaling pathway, oxidative phosphorylation pathway, and RNA transport were most active during primordium differentiation. Furthermore, it was predicted that the blue-light photoreceptor WC-1 and Deoxyribodipyrimidine photolyase PHR play important roles in the primordium differentiation of P. eryngii. Taken together, the results of this study provide a speculative mechanism that light induces primordium differentiation and a foundation for further research on fruiting body development in P. eryngii.

AvNAC030, a NAC Domain Transcription Factor, Enhances Salt Stress Tolerance in Kiwifruit

Kiwifruit (Actinidia chinensis Planch) is suitable for neutral acid soil. However, soil salinization is increasing in kiwifruit production areas, which has adverse effects on the growth and development of plants, leading to declining yields and quality. Therefore, analyzing the salt tolerance regulation mechanism can provide a theoretical basis for the industrial application and germplasm improvement of kiwifruit. We identified 120 NAC members and divided them into 13 subfamilies according to phylogenetic analysis. Subsequently, we conducted a comprehensive and systematic analysis based on the conserved motifs, key amino acid residues in the NAC domain, expression patterns, and protein interaction network predictions and screened the candidate gene AvNAC030. In order to study its function, we adopted the method of heterologous expression in Arabidopsis. Compared with the control, the overexpression plants had higher osmotic adjustment ability and improved antioxidant defense mechanism. These results suggest that AvNAC030 plays a positive role in the salt tolerance regulation mechanism in kiwifruit.

Genome-Wide Identification, Genomic Organization, and Characterization of Potassium Transport-Related Genes in Cajanus cajan and Their Role in Abiotic Stress

Potassium is the most important and abundant inorganic cation in plants and it can comprise up to 10% of a plant's dry weight. Plants possess complex systems of transporters and channels for the transport of K+ from soil to numerous parts of plants. Cajanus cajan is cultivated in different regions of the world as an economical source of carbohydrates, fiber, proteins, and fodder for animals. In the current study, 39 K+ transport genes were identified in C. cajan, including 25 K+ transporters (17 carrier-like K+ transporters (KUP/HAK/KTs), 2 high-affinity potassium transporters (HKTs), and 6 K+ efflux transporters (KEAs) and 14 K+ channels (9 shakers and 5 tandem-pore K+ channels (TPKs). Chromosomal mapping indicated that these genes were randomly distributed among 10 chromosomes. A comparative phylogenetic analysis including protein sequences from Glycine max, Arabidopsis thaliana, Oryza sativa, Medicago truncatula Cicer arietinum, and C. cajan suggested vital conservation of K+ transport genes. Gene structure analysis showed that the intron/exon organization of K+ transporter and channel genes is highly conserved in a family-specific manner. In the promoter region, many cis-regulatory elements were identified related to abiotic stress, suggesting their role in abiotic stress response. Abiotic stresses (salt, heat, and drought) adversely affect chlorophyll, carotenoids contents, and total soluble proteins. Furthermore, the activities of catalase, superoxide, and peroxidase were altered in C. cajan leaves under applied stresses. Expression analysis (RNA-seq data and quantitative real-time PCR) revealed that several K+ transport genes were expressed in abiotic stress-responsive manners. The present study provides an in-depth understanding of K+ transport system genes in C. cajan and serves as a basis for further characterization of these genes.

Genome-Wide Analysis and the Expression Pattern of the MADS-Box Gene Family in Bletilla striata

Bletilla striata (Thunb. ex A. Murray) Rchb. f., a species of the perennial herb Orchidaceae, has potent anti-inflammatory and antiviral biological activities. MADS-box transcription factors play critical roles in the various developmental processes of plants. Although this gene family has been extensively investigated in many species, it has not been analyzed for B. striata. In total, 45 MADS-box genes were identified from B. striata in this study, which were classified into five subfamilies (Mδ, MIKC, Mα, Mβ, and Mγ). Meanwhile, the highly correlated protein domains, motif compositions, and exon-intron structures of BsMADSs were investigated according to local B. striata databases. Chromosome distribution and synteny analyses revealed that segmental duplication and homologous exchange were the main BsMADSs expansion mechanisms. Further, RT-qPCR analysis revealed that BsMADSs had different expression patterns in response to various stress treatments. Our results provide a potential theoretical basis for further investigation of the functions of MADS genes during the growth of B. striata.

Genomic-Wide Identification and Characterization of the Uridine Diphosphate Glycosyltransferase Family in Eucommia ulmoides Oliver

Eucommia ulmoides Oliver is a woody plant with great economic and medicinal value. Its dried bark has a long history of use as a traditional medicinal material in East Asia, which led to many glycosides, such as aucubin, geniposide, hyperoside, astragalin, and pinoresinol diglucoside, being recognized as pharmacologically active ingredients. Uridine diphosphate glycosyltransferases (UGTs) catalyze a glycosyl-transferring reaction from the donor molecule uridine-5'-diphosphate-glucose (UDPG) to the substrate, which plays an important role in many biological processes, such as plant growth and development, secondary metabolism, and environmental adaptation. In order to explore the biosynthetic pathways of glycosides in E. ulmoides, 91 putative EuUGT genes were identified throughout the complete genome of E. ulmoides through function annotation and an UDPGT domain search. Phylogenetic analysis categorized them into 14 groups. We also performed GO annotations on all the EuUGTs to gain insights into their functions in E. ulmoides. In addition, transcriptomic analysis indicated that most EuUGTs showed different expression patterns across diverse organs and various growing seasons. By protein-protein interaction predication, a biosynthetic routine of flavonoids and their glycosides was also proposed. Undoubtedly, these results will help in future research into the biosynthetic pathways of glycoside compounds in E. ulmoides.