Genome-wide identification of microRNAs responding to early stages of phosphate deficiency in maize

Integration of small RNA, degradome, and transcriptome sequencing data illustrates the mechanism of low phosphorus adaptation in Camellia oleifera

Phosphorus (P) is an indispensable macronutrient for plant growth and development, and it is involved in various cellular biological activities in plants. Camellia oleifera is a unique high-quality woody oil plant that grows in the hills and mountains of southern China. However, the available P content is deficient in southern woodland soil. Until now, few studies focused on the regulatory functions of microRNAs (miRNAs) and their target genes under low inorganic phosphate (Pi) stress. In this study, we integrated small RNA, degradome, and transcriptome sequencing data to investigate the mechanism of low Pi adaptation in C. oleifera. We identified 40,689 unigenes and 386 miRNAs by the deep sequencing technology and divided the miRNAs into four different groups. We found 32 miRNAs which were differentially expressed under low Pi treatment. A total of 414 target genes of 108 miRNAs were verified by degradome sequencing. Gene ontology (GO) functional analysis of target genes found that they were related to the signal response to the stimulus and transporter activity, indicating that they may respond to low Pi stress. The integrated analysis revealed that 31 miRNA-target pairs had negatively correlated expression patterns. A co-expression regulatory network was established based on the profiles of differentially expressed genes. In total, three hub genes (ARF22, WRKY53, and SCL6), which were the targets of differentially expressed miRNAs, were discovered. Our results showed that integrated analyses of the small RNA, degradome, and transcriptome sequencing data provided a valuable basis for investigating low Pi in C. oleifera and offer new perspectives on the mechanism of low Pi tolerance in woody oil plants.

Combined Transcriptome and Proteome Analysis of Maize (Zea mays L.) Reveals A Complementary Profile in Response to Phosphate Deficiency

A deficiency in the macronutrient phosphate (Pi) brings about various changes in plants at the morphological, physiological and molecular levels. However, the molecular mechanism for regulating Pi homeostasis in response to low-Pi remains poorly understood, particularly in maize (Zea mays L.), which is a staple crop and requires massive amounts of Pi. Therefore, in this study, we performed expression profiling of the shoots and roots of maize seedlings with Pi-tolerant genotype at both the transcriptomic and proteomic levels using RNA sequencing and isobaric tags for relative and absolute quantitation (iTRAQ). We identified 1944 differentially expressed transcripts and 340 differentially expressed proteins under low-Pi conditions. Most of the differentially expressed genes were clustered as regulators, such as transcription factors involved in the Pi signaling pathway at the transcript level. However, the more functional and metabolism-related genes showed expression changes at the protein level. Moreover, under low-Pi conditions, Pi transporters and phosphatases were specifically induced in the roots at both the transcript and protein levels, and increased amounts of mRNA and protein of two purple acid phosphatases (PAPs) and one UDP-sulfoquinovose synthase (SQD) were specifically detected in the roots. The new insights provided by this study will help to improve the P-utilization efficiency of maize.

Identification and characterization of circRNAs in maize seedlings under deficient nitrogen

Here, deep sequencing results of the maize transcriptome in leaves and roots were compared under high-nitrogen (HN) and low-nitrogen (LN) conditions to identify differentially expressed circRNAs (DECs). Circular RNAs (circRNAs) are covalently closed non-coding RNA with widely regulatory potency that has been identified in animals and plants. However, the understanding of circRNAs involved in responsive nitrogen deficiency remains to be elucidated. A total of 24 and 22 DECs were obtained from the leaves and roots, respectively. Ten circRNAs were validated by divergent and convergent primers, and 6 DECs showed the same expression tendency validated by reverse transcriptase-quantitative PCR. Integrating the identified differentially expressed miRNAs, 34 circRNAs could act as miRNA decoys, which might play important roles in multiple biological processes, including organonitrogen compound biosynthesis and regulation of the metabolic process. A total of 51 circRNA-parent genes located in the genome-wide association study identified loci were assessed between HN and LN conditions and were associated with root growth and development. In summary, our results provide valuable information regarding further study of maize circRNAs under nitrogen deficiency and provide new insights into screening of candidate genes as well as the improvement of maize regarding nitrogen deficiency resistance. CircRNA-miRNA-mRNA co-expression networks were constructed to explore the circRNAs that participated in biological development and nitrogen metabolism.

MicroRNAs: Potential Targets for Developing Stress-Tolerant Crops

Crop yield is challenged every year worldwide by changing climatic conditions. The forecasted climatic scenario urgently demands stress-tolerant crop varieties to feed the ever-increasing global population. Molecular breeding and genetic engineering approaches have been frequently exploited for developing crops with desired agronomic traits. Recently, microRNAs (miRNAs) have emerged as powerful molecules, which potentially serve as expression markers during stress conditions. The miRNAs are small non-coding endogenous RNAs, usually 20-24 nucleotides long, which mediate post-transcriptional gene silencing and fine-tune the regulation of many abiotic- and biotic-stress responsive genes in plants. The miRNAs usually function by specifically pairing with the target mRNAs, inducing their cleavage or repressing their translation. This review focuses on the exploration of the functional role of miRNAs in regulating plant responses to abiotic and biotic stresses. Moreover, a methodology is also discussed to mine stress-responsive miRNAs from the enormous amount of transcriptome data available in the public domain generated using next-generation sequencing (NGS). Considering the functional role of miRNAs in mediating stress responses, these molecules may be explored as novel targets for engineering stress-tolerant crop varieties.

Transcriptomic and genome-wide association study reveal long noncoding RNAs responding to nitrogen deficiency in maize

BACKGROUND

Long noncoding RNAs (lncRNAs) play important roles in essential biological processes. However, our understanding of lncRNAs as competing endogenous RNAs (ceRNAs) and their responses to nitrogen stress is still limited.

RESULTS

Here, we surveyed the lncRNAs and miRNAs in maize inbred line P178 leaves and roots at the seedling stage under high-nitrogen (HN) and low-nitrogen (LN) conditions using lncRNA-Seq and small RNA-Seq. A total of 894 differentially expressed lncRNAs and 38 different miRNAs were identified. Co-expression analysis found that two lncRNAs and four lncRNA-targets could competitively combine with ZmmiR159 and ZmmiR164, respectively. To dissect the genetic regulatory by which lncRNAs might enable adaptation to limited nitrogen availability, an association mapping panel containing a high-density single-nucleotide polymorphism (SNP) array (56,110 SNPs) combined with variable LN tolerant-related phenotypes obtained from hydroponics was used for a genome-wide association study (GWAS). By combining GWAS and RNA-Seq, 170 differently expressed lncRNAs within the range of significant markers were screened. Moreover, 40 consistently LN-responsive genes including those involved in glutamine biosynthesis and nitrogen acquisition in root were identified. Transient expression assays in Nicotiana benthamiana demonstrated that LNC_002923 could inhabit ZmmiR159-guided cleavage of Zm00001d015521.

CONCLUSIONS

These lncRNAs containing trait-associated significant SNPs could consider to be related to root development and nutrient utilization. Taken together, the results of our study can provide new insights into the potential regulatory roles of lncRNAs in response to LN stress, and give valuable information for further screening of candidates as well as the improvement of maize resistance to LN stress.

Transcriptional, metabolic, physiological and developmental responses of switchgrass to phosphorus limitation

Knowing how switchgrass (Panicum virgatum L.) responds and adapts to phosphorus (P)-limitation will aid efforts to optimize P acquisition and use in this species for sustainable biomass production. This integrative study investigated the impacts of mild, moderate, and severe P-stress on genome transcription and whole-plant metabolism, physiology and development in switchgrass. P-limitation reduced overall plant growth, increased root/shoot ratio, increased root branching at moderate P-stress, and decreased root diameter with increased density and length of root hairs at severe P-stress. RNA-seq analysis revealed thousands of genes that were differentially expressed under moderate and severe P-stress in roots and/or shoots compared to P-replete plants, with many stress-induced genes involved in transcriptional and other forms of regulation, primary and secondary metabolism, transport, and other processes involved in P-acquisition and homeostasis. Amongst the latter were multiple miRNA399 genes and putative targets of these. Metabolite profiling showed that levels of most sugars and sugar alcohols decreased with increasing P stress, while organic and amino acids increased under mild and moderate P-stress in shoots and roots, although this trend reversed under severe P-stress, especially in shoots.

Genome-wide identification of low phosphorus responsive microRNAs in two soybean genotypes by high-throughput sequencing

MicroRNAs (miRNAs) have been reported to be correlated with various stress responses in soybean, but only a few miRNAs have been demonstrated to respond to low phosphorus (LP) stress. To unravel the response mechanisms of miRNAs to low-P stress, the roots of two representative soybean genotypes with different P efficiency, Nannong94-156 (a LP-tolerant genotype) and Bogao (a LP-sensitive genotype), were used for the construction of RNA sequencing (RNA-seq) libraries under low/normal-P treatment by high-throughput sequencing. In total, 603 existing miRNAs and 1699 novel miRNAs belonging to 248 and 1582 families in all samples were identified, respectively. Among these miRNAs, 777 miRNAs were differentially expressed (DE) across different P levels and genotypes. Furthermore, putative targets of DE miRNAs were predicted, and these miRNAs mainly targeted ERF (ethylene responsive factor), auxin response factors (ARF), zinc finger protein, MYB, and NAC domain transcription factors. Gene ontology (GO) analysis showed that targets of DE miRNAs were significantly enriched in binding, metabolic processes, biological regulation, response to stress, and phosphorus metabolic processes. In addition, the expression profiles of chosen P-responsive miRNAs and target genes were validated by quantitative real-time PCR (qRT-PCR). Our study focused on genome-wide miRNA identification in two representative soybean genotypes under low-P stress. Overall, the DE miRNAs across different P levels and genotypes and their putative target genes will provide useful information for further study of miRNAs mediating low-P response and facilitate improvements in soybean breeding.

Epigenomic landscape and epigenetic regulation in maize

Epigenetic regulation has been implicated in the control of multiple agronomic traits in maize. Here, we review current advances in our understanding of epigenetic regulation, which has great potential for improving agronomic traits and the environmental adaptability of crops. Epigenetic regulation plays vital role in the control of complex agronomic traits. Epigenetic variation could contribute to phenotypic diversity and can be used to improve the quality and productivity of crops. Maize (Zea mays L.), one of the most widely cultivated crops for human food, animal feed, and ethanol biofuel, is a model plant for genetic studies. Recent advances in high-throughput sequencing technology have made possible the study of epigenetic regulation in maize on a genome-wide scale. In this review, we discuss recent epigenetic studies in maize many achieved by Chinese research groups. These studies have explored the roles of DNA methylation, posttranslational modifications of histones, chromatin remodeling, and noncoding RNAs in the regulation of gene expression in plant development and environment response. We also provide our future prospects for manipulating epigenetic regulation to improve crops.

Genome-Wide Association Study of 13 Traits in Maize Seedlings under Low Phosphorus Stress

Low P stress is a global issue for grain production. Significant phenotypic differences were detected among 13 traits in 356 maize lines under P-sufficient and P-deficient conditions. Significant single nucleotide polymorphisms (SNPs) and low-P stress-responsive genes were identified for 13 maize root traits based on a genome-wide association study. Hap5, harboring 12 favorable SNPs, could enhance strong root systems and P absorption under low-P stress. Phosphorus is an essential macronutrient required for normal plant growth and development. Determining the genetic basis of root traits will enhance our understanding of maize's (Zea mays L.) tolerance to low-P stress. Here, we identified significant phenotypic differences for 13 traits in maize seedlings subjected to P-sufficient and P-deficient conditions. Six extremely sensitive and seven low-P stress tolerant inbreds were selected from 356 inbred lines of maize. No significant differences were observed between temperate and tropical-subtropical groups with respect to trait ratios associated with the adaptation to low-P stress. The broad-sense heritability of these traits ranged from relatively moderate (0.59) to high (0.90). Through genome-wide association mapping with 541,575 informative single nucleotide polymorphisms (SNPs), 551, 1140 and 1157 significant SNPs were detected for the 13 traits in 2012, 2016 and both years combined, respectively, along with 23 shared candidate genes, seven of which overlapped with reported quantitative trait loci and genes for low-P stress. Five haplotypes located in candidate gene GRMZM2G009544 were identified; among these, Hap5, harboring 12 favorable SNP alleles, showed significantly greater values for the root traits studied than the other four haplotypes under both experimental conditions. The candidate genes and favorable haplotypes and alleles identified here provide promising resources for genetic studies and molecular breeding for improving tolerance to abiotic stress in maize.

Comparative profile analysis reveals differentially expressed microRNAs regulate anther and pollen development in kenaf cytoplasmic male sterility line

Cytoplasmic male sterility (CMS) is advantageous in extensive crop breeding and represents a perfect model for understanding anther and pollen development research. MicroRNAs (miRNAs) play key roles in regulating various biological processes. However, the miRNA-mediated regulatory network in kenaf CMS occurrence remains largely unknown. In the present study, a comparative deep sequencing approach was used to investigate the miRNAs and their roles in regulating anther and pollen development during CMS occurrence. We identified 283 known and 46 new candidate miRNAs in kenaf anther. A total of 67 differentially expressed miRNAs (DEMs) were discovered between CMS and its maintainer line. Among them, 40 and 27 miRNAs were up- and downregulated, respectively. These 67 DEMs were predicted to target 189 genes. Validation of DEMs and putative target genes were confirmed by using real-time quantitative PCR. In addition, a potential miRNA-mediated regulatory network, which mainly involves the auxin signaling pathway, signal transduction, glycolysis and energy metabolism, gene expression, transmembrane transport, protein modification and metabolism, and floral development, that mediates anther development during CMS occurrence was proposed. Taken together, our findings provide a better understanding of the molecular mechanism of miRNA regulation in pollen development and CMS occurrence in kenaf.

Identification of Arbuscular Mycorrhiza Fungi Responsive microRNAs and Their Regulatory Network in Maize

Maize can form symbiotic relationships with arbuscular mycorrhiza (AM) fungus to increase productivity and resistance, but the miRNAs in maize responsible for this process have not been discovered. In this study, 155 known and 28 novel miRNAs were identified by performing high-throughput sequencing of sRNA in maize roots colonized by AM fungi. Similar to the profiles in other AM-capable plants, a large proportion of identified maize miRNAs were 24 nt in length. Fourteen and two miRNAs were significantly down- and up-regulated in response to AM fungus Glomus intraradices inoculation, respectively, suggesting potential roles of these miRNAs in AM symbiosis. Interestingly, 12 of 14 significantly down-regulated known maize miRNAs belong to the miR399 family, which was previously reported to be involved in the interaction between Medicago truncatula and AM fungi. This result indicated that the miR399 family should regulate AM symbiosis conservatively across different plant lineages. Pathway and network analyses showed that the differentially expressed miRNAs might regulate lipid metabolism and phosphate starvation response in maize during the symbiosis process via their target genes. Several members of the miR399 family and the miR397 family should be involved in controlling the fatty acid metabolism and promoting lipid delivering from plants to AM fungi. To the best of our knowledge, this is the first report on miRNAs mediating fatty acids from plant to AM fungi. This study provides insight into the regulatory roles of miRNAs in the symbiosis between plants and AM fungi.

Genome-wide analysis of miRNAs and Tasi-RNAs in Zea mays in response to phosphate deficiency

Globally important cereal crop maize provides important nutritions and starch in dietary foods. Low phosphate (LPi) availability in the soil frequently limits the maize quality and yield across the world. Small non-coding RNAs (Snc-RNAs) play crucial roles in growth and adaptation of plants to the environment. Snc-RNAs like microRNAs (miRs) and trans-acting small interfering RNAs (Tasi-Rs) play important functions in posttranscriptional regulation of gene expression, which controls plant development, reproduction, and biotic/abiotic stress responses. In order to identify the miR and Tasi-R alterations in leaf and root of maize in response to sufficient phosphate and LPi at 3LS and 4LS, the snc-RNA population libraries for 0th, 1st, 2nd, 4th, and 8th day were constructed. These libraries were used for genome-wide alignment and RNA-fold analysis for possible prediction of potential miRs and Tasi-Rs. This study reported 174 known and conserved differentially expressed miRs of 27 miR families of maize plant. In addition, leaf and root specific potential novel miRs representing 155 new families were also discovered. Differentially expressed conserved as well as novel miR functions in root and leaf during early stage of Pi starvation were extensively discussed. Leaf and root specific miRs as well as common miRs with their target genes, participating in different biological, cellular, and metabolic processes were explored. Further, four miR390-directed Tasi-Rs which belong to TAS3 gene family along with other orthologs of Tasi-Rs were also identified. Finally, the study provides an insight into the composite regulatory mechanism of miRs in maize in response to Pi deficiency.

Identification and characterization of microRNAs in maize endosperm response to exogenous sucrose using small RNA sequencing

Sucrose acts as a signaling molecule for genes critical to starch biosynthesis in maize endosperm. Previously, we showed that sucrose could regulate starch biosynthesis in maize via transcription factors. To better understand the complex regulation of starch biosynthesis, the 10days after pollination endosperm from Zea mays L. B73 inbred line was collected and treated with sucrose for small RNA sequencing. The sequencing results revealed that 24 known miRNAs and 190 novel miRNAs were significantly differentially expressed in response to sucrose. In addition, most of target mRNAs were characterized as transcription factors, mainly including, MYB, ARF, NAC, AP2/ERF, WRKY, and GRAS, which play important roles in starch biosynthesis and seed development in maize endosperm. The expression profiles of miR398a/b and miR159b/j/k followed opposite expression trends to their target genes when analyzed by qPCR. In conclusion, these results show that sucrose regulates the expression of starch synthetic genes through miRNAs.