Genome wide identification and characterization of abiotic stress responsive lncRNAs in Capsicum annuum

Systematic profiling and analysis of growth and development responsive DE-lncRNAs in cluster bean (Cyamopsis tetragonoloba)

Long non-coding RNAs (lncRNAs) play crucial role in regulating genes involved in various processes including growth & development, flowering, and stress response in plants. The study aims to identify and characterize tissue-specific, growth & development and floral responsive differentially expressed lncRNAs (DE-lncRNAs) in cluster bean from a high-throughput RNA sequencing data. We have identified 3309 DE-lncRNAs, with an average length of 818 bp. Merely, around 4 % of DE-lncRNAs across the tissues were found to be conserved as rate of evolution of lncRNAs is high. Among the identified DE-lncRNAs, 204 were common in leaf vs. shoot, leaf vs. flower and flower vs. shoot. A total of 60 DE-lncRNAs targeted 10 protein-coding genes involved in flower development and initiation processes. We investigated 179 tissue-specific DE-lncRNAs based on tissue specificity index. Three DE-lncRNAs: Cb_lnc_0820, Cb_lnc_0430, Cb_lnc_0260 and their target genes show their involvement in floral development and stress mechanisms, which were validated by Quantitative real-time PCR (qRT-PCR). The identified DE-lncRNAs were expressed higher in flower bud than in leaf and similar expression pattern was observed in both RNA-seq data and qRT-PCR analyses. Notably, 362 DE-lncRNAs were predicted as eTM-lncRNAs with the participation of 84 miRNAs. Whereas 46 DE-lncRNAs were predicted to possess the internal ribosomal entry sites (IRES) and can encode for small peptides. The regulatory networks established between DE-lncRNAs, mRNAs and miRNAs have provided an insight into their association with plant growth & development, flowering, and stress mechanisms. Comprehensively, the characterization of DE-lncRNAs in various tissues of cluster bean shed a light on interactions among lncRNAs, miRNAs and mRNAs and help understand their involvement in growth & development and floral initiation processes. The information retrieved from the analyses was shared in the public domain in the form of a database: Cb-DElncRNAdb, and made available at http://backlin.cabgrid.res.in/Cb-DElncRNA/index.php, which may be useful for the scientific community engaged cluster bean research.

Predicted roles of long non-coding RNAs in abiotic stress tolerance responses of plants

The plant genome exhibits a significant amount of transcriptional activity, with most of the resulting transcripts lacking protein-coding potential. Non-coding RNAs play a pivotal role in the development and regulatory processes in plants. Long non-coding RNAs (lncRNAs), which exceed 200 nucleotides, may play a significant role in enhancing plant resilience to various abiotic stresses, such as excessive heat, drought, cold, and salinity. In addition, the exogenous application of chemicals, such as abscisic acid and salicylic acid, can augment plant defense responses against abiotic stress. While how lncRNAs play a role in abiotic stress tolerance is relatively well-studied in model plants, this review provides a comprehensive overview of the current understanding of this function in horticultural crop plants. It also delves into the potential role of lncRNAs in chemical priming of plants in order to acquire abiotic stress tolerance, although many limitations exist in proving lncRNA functionality under such conditions.

A lncRNA bra-miR156HG regulates flowering time and leaf morphology as a precursor of miR156 in Brassica campestris and Arabidopsis thaliana

Long non-coding RNAs (lncRNAs) are important regulators in plant growth and development. Here the function of a lncRNA fragment was studied, which was predicted as an endogenous target mimic (eTM) of miR156 in Brassica campesrtis. Unexpectedly, the transformation of this lncRNA into Arabidopsis thaliana neither inhibited the expression of miR156a nor resulted in any phenotypes that differed from the control plants (CK). The full-length sequence of the lncRNA (named bra-miR156HG) was then obtained using RACE and transferred into A. thaliana. The transgenic plants displayed a delay in flowering time, an increasing number of rosette leaves, and a changed morphology of cauline leaves, which was similar to the plants that expressed bra-miR156a. In contrast, the overexpression of bra-miR156HG in B. campestris resulted in an increased tip angle of leaves and changed the length-width ratio of leaves at different nodes, suggesting that bra-miR156HG may be involved in regulating the leaf morphology. Collectively, our study showed that bra-miR156HG functions as a precursor of bra-miR156a involved in regulating plant flowering time and leaf development under different biological backgrounds. The secondary structure of lncRNA is essential not only for the normal roles that it plays but also for expanding the functional diversities.

Novel transcriptome networks are associated with adaptation of capsicum fruit development to a light-blocking glasshouse film

Light-blocking films (LBFs) can contribute to significant energy savings for protected cropping via altering light transmitting, such as UVA, photosynthetically active radiation, blue and red spectra affecting photosynthesis, and capsicum yield. Here, we investigated the effects of LBF on orange color capsicum (O06614, Capsicum annuum L.) fruit transcriptome at 35 (mature green) and 65 (mature ripe) days after pollination (DAP) relative to untreated control in a high-technology glasshouse. The results of targeted metabolites showed that LBF significantly promotes the percentage of lutein but decreased the percentage of zeaxanthin and neoxanthin only at 35 DAP. At 35 DAP, fruits were less impacted by LBF treatment (versus control) with a total of 1,192 differentially expressed genes (DEGs) compared with that at 65 DAP with 2,654 DEGs. Response to stress and response to light stimulus in biological process of Gene Ontology were found in 65-DAP fruits under LBF vs. control, and clustering analysis revealed a predominant role of light receptors and phytohormone signaling transduction as well as starch and sucrose metabolism in LBF adaptation. The light-signaling DEGs, UV light receptor UVR8, transcription factors phytochrome-interacting factor 4 (PIF4), and an E3 ubiquitin ligase (COP1) were significantly downregulated at 65 DAP. Moreover, key DEGs in starch and sucrose metabolism (SUS, SUC, and INV), carotenoid synthesis (PSY2 and BCH1), ascorbic acid biosynthesis (VTC2, AAO, and GME), abscisic acid (ABA) signaling (NCED3, ABA2, AO4, and PYL2/4), and phenylpropanoid biosynthesis (PAL and DFR) are important for the adaptation of 65-DAP fruits to LBF. Our results provide new candidate genes for improving quality traits of low-light adaptation of capsicum in protected cropping.

CircRNAs responsive to winter dormancy and spring flushing conditions of tea leaf buds

Circular RNAs (circRNAs) are important regulators of diverse biological processes of plants. However, the evolution and potential functions of circRNAs during winter dormancy and spring bud flushing of tea plant is largely unknown. Using RNA-seq data, a total of 1184 circRNAs were identified in the winter dormant and spring bud flushing leaf samples of tea plants in two different cultivars exhibiting different duration of winter dormancy. A total of 156 circRNAs are found to be differentially expressed and the weighted gene co-expression network (WGCNA) analysis revealed that 22 and 20 differentially expressed-circRNAs (DE-circRNAs) positively correlated with the flushing and dormant leaf traits, respectively, in both the tea cultivars used. Some transcription factors (TFs) viz. MYB, WRKY, ERF, bHLH and several genes related to secondary metabolite biosynthetic pathways are found to co-express with circRNAs. DE-circRNAs also predicted to interact with miRNAs and can regulate phytohormone biosynthesis and various signalling pathways in tea plant. This study uncovers the potential roles of circRNAs to determine winter dormancy and spring bud flushing conditions in tea plants.

The plant noncoding transcriptome: a versatile environmental sensor

Plant noncoding RNA transcripts have gained increasing attention in recent years due to growing evidence that they can regulate developmental plasticity. In this review article, we comprehensively analyze the relationship between noncoding RNA transcripts in plants and their response to environmental cues. We first provide an overview of the various noncoding transcript types, including long and small RNAs, and how the environment modulates their performance. We then highlight the importance of noncoding RNA secondary structure for their molecular and biological functions. Finally, we discuss recent studies that have unveiled the functional significance of specific long noncoding transcripts and their molecular partners within ribonucleoprotein complexes during development and in response to biotic and abiotic stress. Overall, this review sheds light on the fascinating and complex relationship between dynamic noncoding transcription and plant environmental responses, and highlights the need for further research to uncover the underlying molecular mechanisms and exploit the potential of noncoding transcripts for crop resilience in the context of global warming.

In silico and computational analysis of zinc finger motif-associated homeodomain (ZF-HD) family genes in chilli (Capsicum annuum L)

Zinc finger-homeodomain (ZHD) proteins are mostly expressed in plants and are involved in proper growth and development and minimizing biotic and abiotic stress. A recent study identified and characterized the ZHD gene family in chilli (Capsicum annuum L.) to determine their probable molecular function. ZHD genes with various physicochemical characteristics were discovered on twelve chromosomes in chilli. We separated ZHD proteins into two major groups using sequence alignment and phylogenetic analysis. These groups differ in gene structure, motif distribution, and a conserved ZHD and micro-zinc finger ZF domain. The majority of the CaZHDs genes are preserved, early duplication occurred recently, and significant pure selection took place throughout evolution, according to evolutionary study. According to expression profiling, the genes were found to be equally expressed in tissues above the ground, contribute to plant growth and development and provide tolerance to biotic and abiotic stress. This in silico analysis, taken as a whole, hypothesized that these genes perform distinct roles in molecular and phytohormone signaling processes, which may serve as a foundation for subsequent research into the roles of these genes in other crops.

A comprehensive review on genomic resources in medicinally and industrially important major spices for future breeding programs: Status, utility and challenges

In the global market, spices possess a high-value but low-volume commodities of commerce. The food industry depends largely on spices for taste, flavor, and therapeutic properties in replacement of cheap synthetic ones. The estimated growth rate for spices demand in the world is ∼3.19%. Since spices grow in limited geographical regions, India is one of the leading producer of spices, contributing 25-30 percent of total world trade. Hitherto, there has been no comprehensive review of the genomic resources of industrially important major medicinal spices to overcome major impediments in varietal improvement and management. This review focuses on currently available genomic resources of 24 commercially significant spices, namely, Ajwain, Allspice, Asafoetida, Black pepper, Cardamom large, Cardamom small, Celery, Chillies, Cinnamon, Clove, Coriander, Cumin, Curry leaf, Dill seed, Fennel, Fenugreek, Garlic, Ginger, Mint, Nutmeg, Saffron, Tamarind, Turmeric and Vanilla. The advent of low-cost sequencing machines has contributed immensely to the voluminous data generation of these spices, cracking the complex genomic architecture, marker discovery, and understanding comparative and functional genomics. This review of spice genomics resources concludes the perspective and way forward to provide footprints by uncovering genome assemblies, sequencing and re-sequencing projects, transcriptome-based studies, non-coding RNA-mediated regulation, organelles-based resources, developed molecular markers, web resources, databases and AI-directed resources in candidate spices for enhanced breeding potential in them. Further, their integration with molecular breeding could be of immense use in formulating a strategy to protect and expand the production of the spices due to increased global demand.

miRNAs for crop improvement

Climate change significantly impacts crop production by inducing several abiotic and biotic stresses. The increasing world population, and their food and industrial demands require focused efforts to improve crop plants to ensure sustainable food production. Among various modern biotechnological tools, microRNAs (miRNAs) are one of the fascinating tools available for crop improvement. miRNAs belong to a class of small non-coding RNAs playing crucial roles in numerous biological processes. miRNAs regulate gene expression by post-transcriptional target mRNA degradation or by translation repression. Plant miRNAs have essential roles in plant development and various biotic and abiotic stress tolerance. In this review, we provide propelling evidence from previous studies conducted around miRNAs and provide a one-stop review of progress made for breeding stress-smart future crop plants. Specifically, we provide a summary of reported miRNAs and their target genes for improvement of plant growth and development, and abiotic and biotic stress tolerance. We also highlight miRNA-mediated engineering for crop improvement and sequence-based technologies available for the identification of miRNAs associated with stress tolerance and plant developmental events.

Role of bZIP transcription factors in the regulation of plant secondary metabolism

MAIN CONCLUSION

This study provides an overview of the structure, classification, regulatory mechanisms, and biological functions of the basic (region) leucine zipper transcription factors and their molecular mechanisms in flavonoid, terpenoid, alkaloid, phenolic acid, and lignin biosynthesis. Basic (region) leucine zippers (bZIPs) are evolutionarily conserved transcription factors (TFs) in eukaryotic organisms. The bZIP TFs are widely distributed in plants and play important roles in plant growth and development, photomorphogenesis, signal transduction, resistance to pathogenic microbes, biotic and abiotic stress, and secondary metabolism. Moreover, the expression of bZIP TFs not only promotes or inhibits the accumulation of secondary metabolites in medicinal plants, but also affects the stress response of plants to the external adverse environment. This paper describes the structure, classification, biological function, and regulatory mechanisms of bZIP TFs. In addition, the molecular mechanism of bZIP TFs regulating the biosynthesis of flavonoids, terpenoids, alkaloids, phenolic acids, and lignin are also elaborated. This review provides a summary for in-depth study of the molecular mechanism of bZIP TFs regulating the synthesis pathway of secondary metabolites and plant molecular breeding, which is of significance for the generation of beneficial secondary metabolites and the improvement of plant varieties.

Genome-wide identification and characterization of tissue-specific non-coding RNAs in black pepper (Piper nigrum L.)

Long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) are the two classes of non-coding RNAs (ncRNAs) present predominantly in plant cells and have various gene regulatory functions at pre- and post-transcriptional levels. Previously deemed as "junk", these ncRNAs have now been reported to be an important player in gene expression regulation, especially in stress conditions in many plant species. Black pepper, scientifically known as Piper nigrum L., despite being one of the most economically important spice crops, lacks studies related to these ncRNAs. From a panel of 53 RNA-Seq datasets of black pepper from six tissues, namely, flower, fruit, leaf, panicle, root, and stem of six black pepper cultivars, covering eight BioProjects across four countries, we identified and characterized a total of 6406 lncRNAs. Further downstream analysis inferred that these lncRNAs regulated 781 black pepper genes/gene products via miRNA-lncRNA-mRNA network interactions, thus working as competitive endogenous RNAs (ceRNAs). The interactions may be various mechanisms like miRNA-mediated gene silencing or lncRNAs acting as endogenous target mimics (eTMs) of the miRNAs. A total of 35 lncRNAs were also identified to be potential precursors of 94 miRNAs after being acted upon by endonucleases like Drosha and Dicer. Tissue-wise transcriptome analysis revealed 4621 circRNAs. Further, miRNA-circRNA-mRNA network analysis showed 432 circRNAs combining with 619 miRNAs and competing for the binding sites on 744 mRNAs in different black pepper tissues. These findings can help researchers to get a better insight to the yield regulation and responses to stress in black pepper in endeavor of higher production and improved breeding programs in black pepper varieties.

Long Non-Coding RNAs of Plants in Response to Abiotic Stresses and Their Regulating Roles in Promoting Environmental Adaption

Abiotic stresses triggered by climate change and human activity cause substantial agricultural and environmental problems which hamper plant growth. Plants have evolved sophisticated mechanisms in response to abiotic stresses, such as stress perception, epigenetic modification, and regulation of transcription and translation. Over the past decade, a large body of literature has revealed the various regulatory roles of long non-coding RNAs (lncRNAs) in the plant response to abiotic stresses and their irreplaceable functions in environmental adaptation. LncRNAs are recognized as a class of ncRNAs that are longer than 200 nucleotides, influencing a variety of biological processes. In this review, we mainly focused on the recent progress of plant lncRNAs, outlining their features, evolution, and functions of plant lncRNAs in response to drought, low or high temperature, salt, and heavy metal stress. The approaches to characterize the function of lncRNAs and the mechanisms of how they regulate plant responses to abiotic stresses were further reviewed. Moreover, we discuss the accumulating discoveries regarding the biological functions of lncRNAs on plant stress memory as well. The present review provides updated information and directions for us to characterize the potential functions of lncRNAs in abiotic stresses in the future.

Identification of YABBY Transcription Factors and Their Function in ABA and Salinity Response in Nelumbo nucifera

The plant-specific transcription factor family YABBY plays important roles in plant responses to biotic and abiotic stresses. Although the function of YABBY has been identified in many species, systematic analysis in lotus (Nelumbo nucifera) is still relatively lacking. The present study aimed to characterize all of the YABBY genes in lotus and obtain better insights into NnYABBYs in response to salt stress by depending on ABA signaling. Here, we identified nine YABBY genes by searching the whole lotus genome based on the conserved YABBY domain. Further analysis showed that these members were distributed on six different chromosomes and named from YABBY1 to YABBY9, which were divided into five subgroups, including YAB1, YAB2, YAB5, INO, and CRC. The analysis of cis-elements in promotors revealed that NnYABBYs could be involved in plant hormone signaling and plant responses to abiotic stresses. Quantitative real-time PCR (qRT-PCR) showed that NnYABBYs could be up-regulated or down-regulated by ABA, fluridone, and salt treatment. Subcellular localization indicated that NnYABBY4, NnYABBY5, and NnYABBY6 were mainly localized in the cell membrane and cytoplasm. In addition, the intrinsic trans-activity of NnYABBY was tested by a Y2H assay, which revealed that NnYABBY4, NnYABBY5, and NnYABBY6 are deprived of such a property. This study provided a theoretical basis and reference for the functional research of YABBY for the molecular breeding of lotus.

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 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.

Comprehensive transcriptomic analysis of two RIL parents with contrasting salt responsiveness identifies polyadenylated and non-polyadenylated flower lncRNAs in chickpea

Salinity severely affects the yield of chickpea. Understanding the role of lncRNAs can shed light on chickpea salt tolerance mechanisms. However, because lncRNAs are encoded by multiple sites within the genome, their classification to reveal functional versatility at the transcriptional and the post-transcriptional levels is challenging. To address this, we deep sequenced 24 salt-challenged flower transcriptomes from two parental genotypes of a RIL population that significantly differ in salt tolerance ability. The transcriptomes for the first time included 12 polyadenylated and 12 non-polyadenylated RNA libraries to a sequencing depth of ~50 million reads. The ab initio transcriptome assembly comprised ~34 082 transcripts from three biological replicates of salt-tolerant (JG11) and salt-sensitive (ICCV2) flowers. A total of 9419 lncRNAs responding to salt stress were identified, 2345 of which were novel lncRNAs specific to chickpea. The expression of poly(A+) lncRNAs and naturally antisense transcribed RNAs suggest their role in post-transcriptional modification and gene silencing. Notably, 178 differentially expressed lncRNAs were induced in the tolerant genotype but repressed in the sensitive genotype. Co-expression network analysis revealed that the induced lncRNAs interacted with the FLOWERING LOCUS (FLC), chromatin remodelling and DNA methylation genes, thus inducing flowering during salt stress. Furthermore, 26 lncRNAs showed homology with reported lncRNAs such as COOLAIR, IPS1 and AT4, thus confirming the role of chickpea lncRNAs in controlling flowering time as a crucial salt tolerance mechanism in tolerant chickpea genotype. These robust set of differentially expressed lncRNAs provide a deeper insight into the regulatory mechanisms controlled by lncRNAs under salt stress.

Genome-Wide Identification of Powdery Mildew Responsive Long Non-Coding RNAs in Cucurbita pepo

Cucurbita pepo L. is an essential economic vegetable crop worldwide, and its production is severely affected by powdery mildew (PM). However, our understanding of the molecular mechanism of PM resistance in C. pepo is very limited. Long non-coding RNAs (lncRNAs) play an important role in regulating plant responses to biotic stress. Here, we systematically identified 2,363 reliably expressed lncRNAs from the leaves of PM-susceptible (PS) and PM-resistant (PR) C. pepo. The C. pepo lncRNAs are shorter in length and expressed at a lower level than the protein-coding transcripts. Among the 2,363 lncRNAs, a total of 113 and 146 PM-responsive lncRNAs were identified in PS and PR, respectively. Six PM-responsive lncRNAs were predicted as potential precursors of microRNAs (miRNAs). In addition, 58 PM-responsive lncRNAs were predicted as targets of miRNAs and one PM-responsive lncRNA was predicted as an endogenous target mimic (eTM). Furthermore, a total of 5,200 potential cis target genes and 5,625 potential trans target genes were predicted for PM-responsive lncRNAs. Functional enrichment analysis showed that these potential target genes are involved in different biological processes, such as the plant-pathogen interaction pathway, MAPK signaling pathway, and plant hormone signal transduction pathway. Taken together, this study provides a comprehensive view of C. pepo lncRNAs and explores the putative functions of PM-responsive lncRNAs, thus laying the foundation for further study of the regulatory mechanisms of lncRNAs responding to PM.

Identification and characterization of heat-responsive lncRNAs in maize inbred line CM1

BACKGROUND

Frequent occurrence of extreme high temperature is a major threat to crop production. Increasing evidence demonstrates that long non-coding RNAs (lncRNAs) have important biological functions in the regulation of the response to heat stress. However, the regulatory mechanism of lncRNAs involved in heat response requires further exploration and the regulatory network remains poorly understood in maize.

RESULTS

In this research, high-throughput sequencing was adopted to systematically identify lncRNAs in maize inbred line CM1. In total, 53,249 lncRNAs (259 known lncRNAs and 52,990 novel lncRNAs) were detected, of which 993 lncRNAs showed significantly differential expression (DElncRNAs) under heat stress. By predicting the target genes, 953 common targets shared by cis- and trans-regulation of the DElncRNAs were identified, which exhibited differential expression between the control and the heat stress treatments. Functional annotation indicated that a number of important biological processes and pathways, including photosynthesis, metabolism, translation, stress response, hormone signal transduction, and spliceosome, were enriched for the common targets, suggesting that they play important roles in heat response. A lncRNA-mediated regulatory network was constructed to visualize the molecular response mechanism in response to heat stress, which represented the direct regulatory relationships of DElncRNAs, differentially expressed miRNAs, target genes, and functional annotations.

CONCLUSIONS

This study lays a foundation for further elucidation of the regulatory mechanism for the response to heat stress in the maize inbred line CM1. The findings provide important information for identification of heat-responsive genes, which will be beneficial for the molecular breeding in the cultivation of heat-tolerant maize germplasm.

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.

Identification and functional deciphering suggested the regulatory roles of long intergenic ncRNAs (lincRNAs) in increasing grafting pepper resistance to Phytophthora capsici

BACKGROUND

As a popular and valuable technique, grafting is widely used to protect against soil-borne diseases and nematodes in vegetable production. Growing evidences have revealed that long intergenic ncRNAs (lincRNAs) are strictly regulated and play essential roles in plants development and stress responses. Nevertheless, genome-wide identification and function deciphering of pepper lincRNAs, especially for their roles in improving grafting pepper resistance to Phytophthora capsici is largely unknown.

RESULTS

In this study, RNA-seq data of grafting and control pepper plants with or without P. capsici inoculation were used to identify lincRNAs. In total, 2,388 reliable lincRNAs were identified. They were relatively longer and contained few exons than protein-coding genes. Similar to coding genes, lincRNAs had higher densities in euchromatin regions; and longer chromosome transcribed more lincRNAs. Expression pattern profiling suggested that lincRNAs commonly had lower expression than mRNAs. Totally, 607 differentially expressed lincRNAs (DE-lincRANs) were identified, of which 172 were found between P. capsici resistance grafting pepper sample GR and susceptible sample LDS. The neighboring genes of DE-lincRNAs and miRNAs competitively sponged by DE-lincRNAs were identified. Subsequently, the expression level of DE-lincRNAs was further confirmed by qRT-PCR and regulation patterns between DE-lincRNAs and neighboring mRNAs were also validated. Function annotation revealed that DE-lincRNAs increased the resistance of grafting prepper to P. capsici by modulating the expression of disease-defense related genes through cis-regulating and/or lincRNA-miRNA-mRNA interaction networks.

CONCLUSIONS

This study identified pepper lincRNAs and suggested their potential roles in increasing the resistance level of grafting pepper to P. capsici.

Genome-Wide Identification of ARF Gene Family Suggests a Functional Expression Pattern during Fruitlet Abscission in Prunus avium L

Auxin response factors (ARFs) play a vital role in plant growth and development. In the current study, 16 ARF members have been identified in the sweet cherry (Prunus avium L.) genome. These genes are all located in the nucleus. Sequence analysis showed that genes in the same subgroup have similar exon-intron structures. A phylogenetic tree has been divided into five groups. The promoter sequence includes six kinds of plant hormone-related elements, as well as abiotic stress response elements such as low temperature or drought. The expression patterns of PavARF in different tissues, fruitlet abscission, cold and drought treatment were comprehensively analyzed. PavARF10/13 was up-regulated and PavARF4/7/11/12/15 was down-regulated in fruitlet abscising. These genes may be involved in the regulation of fruit drop in sweet cherry fruits. This study comprehensively analyzed the bioinformatics and expression pattern of PavARF, which can lay the foundation for further understanding the PavARF family in plant growth development and fruit abscission.

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 and Functional Analysis of the Basic Helix-Loop-Helix (bHLH) Transcription Family Reveals Candidate PtFBH Genes Involved in the Flowering Process of Populus trichocarpa

As one of the largest TF families+ in plants, the basic helix-loop-helix (bHLH) family plays an important part in the growth and development of many plants. FLOWERING BHLH (FBH) encodes a bHLH-type transcriptional factor related to the flowering process. Poplar is a model woody plant as well as an important economic tree species with a small genome. However, the characteristics of the bHLHs and FBHs gene family in the newest version of Populustrichocarpa genome have not been analyzed yet. We identified 233 PtbHLHs and 10 PtFBHs in the newest version genome, and PtbHLHs were classified into 21 groups with FBH subfamily occupying one, supported by phylogenetic analysis, exon–intron patterns, and conserved protein motifs. These PtHLHs were distributed on 19 chromosomes unevenly and expressed in nucleus mainly. Gene duplication and synteny analysis have indicated that the PtbHLHs gene family has undergone strong purification selection during the evolution process. The cis-elements analysis has suggested that PtbHLHs may be related to the growth and development. Conserved residues of FBHs among Arabidopsis and poplar were also identified. Expression of 227 PtHLH genes (6 unmatched, 13 no expressed) showed diverse patterns in different tissues, implying their multiple functions. Protein–protein interaction network prediction and expression patterns in three states of the flowering process (Flowers-Dormant, Flowers-Expanding and Flowers-Expanded) suggested that some members of PtbHLH and PtFBH family may be involved in the flowering process. Our comprehensive and systematic analysis can provide some valuable clues and basic reference toward further investigations on physiological and molecular functions of PtbHLHs.

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.

Global identification of long non-coding RNAs involved in the induction of spinach flowering

Background

Spinach is a beneficial annual vegetable species and sensitive to the bolting or early flowering, which causes a large reduction in quality and productivity. Indeed, bolting is an event induced by the coordinated effects of various environmental factors and endogenous genetic components. Although some key flowering responsive genes have been identified in spinach, non-coding RNA molecules like long non-coding RNAs (lncRNAs) were not investigated yet. Herein, we used bioinformatic approaches to analyze the transcriptome datasets from two different accessions Viroflay and Kashan at two vegetative and reproductive stages to reveal novel lncRNAs and the construction of the lncRNA-mRNA co-expression network. Additionally, correlations among gene expression modules and phenotypic traits were investigated; day to flowering was chosen as our interesting trait.

Results

In the present study, we identified a total of 1141 lncRNAs, of which 111 were differentially expressed between vegetative and reproductive stages. The GO and KEGG analyses carried out on the cis target gene of lncRNAs showed that the lncRNAs play an important role in the regulation of flowering spinach. Network analysis pinpointed several well-known flowering-related genes such as ELF, COL1, FLT, and FPF1 and also some putative TFs like MYB, WRKY, GATA, and MADS-box that are important regulators of flowering in spinach and could be potential targets for lncRNAs.

Conclusions

This study is the first report on identifying bolting and flowering-related lncRNAs based on transcriptome sequencing in spinach, which provides a useful resource for future functional genomics studies, genes expression researches, evaluating genes regulatory networks and molecular breeding programs in the regulation of the genetic mechanisms related to bolting in spinach.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07989-1.

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.