An integrated analysis of mRNA and sRNA transcriptional profiles in Coffea arabica L. roots: insights on nitrogen starvation responses

Transcriptomic Data Analysis Using the Galaxy Platform: Coffee (Coffea arabica L.) Flowers as Example

Coffee, an important agricultural product for tropical producing countries, is facing challenges due to climate change, including periods of drought, irregular rain distribution, and high temperatures. These changes result in plant water stress, leading to significant losses in coffee productivity and quality. Understanding the processes that affect coffee flowering is crucial for improving productivity and quality. In this chapter, we describe a protocol for transcriptome analysis using available Internet software, mainly in the Galaxy Platform, using RNA-Seq data from flowers collected from different parts of the coffee tree. The methods presented in this chapter provide a comprehensive protocol for transcriptome analysis of differentially expressed genes from flowers of coffee plant. This knowledge can be utilized in coffee genetic improvement programs, particularly in the selection of cultivars that are tolerant to water deficit.

Using soil nitrogen amendments in mitigating ozone stress in agricultural crops: a case study of cluster beans

The climate change scenario in the coming years is liable to have serious negative consequences on agricultural productivity. Increasing tropospheric ozone concentration is an important aspect of climate change, which, due to its oxidative nature, is injurious to the plants. Due to the multifarious nature and continuously increasing concentration of tropospheric ozone, it is prerequisite to develop strategies to manage ozone stress in plants. Present study not only evaluates the potential of soil nitrogen amendments in ameliorating ozone stress in plants, but also focuses upon the mechanistic approaches adopted by the different plant cultivars to combat ozone stress. Three doses of nitrogen amendments, recommended (N1), 1.5× recommended (N2) and 2× recommended (N3), were given to two cultivars (S-151 and PUSA-N) of Cymopsis tetragonoloba exposed to ambient ozone stress. Control plants were also maintained in which no nitrogen treatment was given. Nitrogen supplementation reduced the root nodulation frequency and leghaemoglobin content, which subsequently increased the cellular nitrogen metabolism as evident through increase in the activities of nitrate reductase and nitrite reductase in both the test cultivars. The positive effects of nitrogen amendments are clearly evident in the 1D protein profile studies which showed a greater accumulation of larger sub-units of RuBisCO in nitrogen amended plants. The results clearly indicate that N2 treatment effectively enhanced the yield of both the cultivars (84.8% and 76.37%, in S-151 and PUSA-N, respectively); however, the mechanistic approach adopted by the two cultivars was different. Whereas the yield quantity showed higher increments in S-151, the yield quality parameters (carbohydrates and nitrogen contents) responded more positively in PUSA-N.

Small RNAs: Promising Molecules to Tackle Climate Change Impacts in Coffee Production

Over the centuries, human society has evolved based on the ability to select and use more adapted species for food supply, which means making plant species tastier and more productive in particular environmental conditions. However, nowadays, this scenario is highly threatened by climate change, especially by the changes in temperature and greenhouse gasses that directly affect photosynthesis, which highlights the need for strategic studies aiming at crop breeding and guaranteeing food security. This is especially worrying for crops with complex phenology, genomes with low variability, and the ones that support a large production chain, such as Coffea sp. L. In this context, recent advances shed some light on the genome function and transcriptional control, revealing small RNAs (sRNAs) that are responsible for environmental cues and could provide variability through gene expression regulation. Basically, sRNAs are responsive to environmental changes and act on the transcriptional and post-transcriptional gene silencing pathways that regulate gene expression and, consequently, biological processes. Here, we first discuss the predicted impact of climate changes on coffee plants and coffee chain production and then the role of sRNAs in response to environmental changes, especially temperature, in different species, together with their potential as tools for genetic improvement. Very few studies in coffee explored the relationship between sRNAs and environmental cues; thus, this review contributes to understanding coffee development in the face of climate change and towards new strategies of crop breeding.

MicroRNAs modulating nutrient homeostasis: a sustainable approach for developing biofortified crops

During their lifespan, sessile plants have to cope with bioavailability of the suboptimal nutrient concentration and have to constantly sense/evolve the connecting web of signal cascades for efficient nutrient uptake, storage, and translocation for proper growth and metabolism. However, environmental fluctuations and escalating anthropogenic activities are making it a formidable challenge for plants. This is adding to (micro)nutrient-deficient crops and nutritional insecurity. Biofortification is emerging as a sustainable and efficacious approach which can be utilized to combat the micronutrient malnutrition. A biofortified crop has an enriched level of desired nutrients developed using conventional breeding, agronomic practices, or advanced biotechnological tools. Nutrient homeostasis gets hampered under nutrient stress, which involves disturbance in short-distance and long-distance cell-cell/cell-organ communications involving multiple cellular and molecular components. Advanced sequencing platforms coupled with bioinformatics pipelines and databases have suggested the potential roles of tiny signaling molecules and post-transcriptional regulators, the microRNAs (miRNAs) in key plant phenomena including nutrient homeostasis. miRNAs are seen as emerging targets for biotechnology-based biofortification programs. Thus, understanding the mechanistic insights and regulatory role of miRNAs could open new windows for exploring them in developing nutrient-efficient biofortified crops. This review discusses significance and roles of miRNAs in plant nutrition and nutrient homeostasis and how they play key roles in plant responses to nutrient imbalances/deficiencies/toxicities covering major nutrients-nitrogen (N), phosphorus (P), sulfur (S), magnesium (Mg), iron (Fe), and zinc (Zn). A perspective view has been given on developing miRNA-engineered biofortified crops with recent success stories. Current challenges and future strategies have also been discussed.

Identification of Nitrogen Starvation-Responsive miRNAs to Reveal the miRNA-Mediated Regulatory Network in Betula luminifera

Because of the immobility, plants encounter a series of stresses, such as varied nutrient concentrations in soil, which regulate plant growth, development, and phase transitions. Nitrogen (N) is one of the most limiting factors for plants, which was exemplified by the fact that low nitrogen (LN) has a great adverse effect on plant growth and development. In the present study, we explored the potential role of microRNAs (miRNAs) in response to LN stress in Betula luminifera. We identified 198 miRNAs using sRNA sequencing, including 155 known and 43 novel miRNAs. Among them, 98 known miRNAs and 31 novel miRNAs were differentially expressed after 0.5 h or 24 h of LN stress. Based on degradome data, 122 differential expressed miRNAs (DEmiRNAs) including 102 known miRNAs and 20 novel miRNAs targeted 203 genes, comprising 321 miRNA–target pairs. A big proportion of target genes were transcription factors and functional proteins, and most of the Gene Ontology terms were enriched in biological processes; moreover, one Kyoto Encyclopedia of Genes and Genomes term “ascorbate and aldarate metabolism” was significantly enriched. The expression patterns of six miRNAs and their corresponding target genes under LN stress were monitored. According to the potential function for targets of DEmiRNAs, a proposed regulatory network mediated by miRNA–target pairs under LN stress in B. luminifera was constructed. Taken together, these findings provide useful information to elucidate miRNA functions and establish a framework for exploring N signaling networks mediated by miRNAs in B. luminifera. It may provide new insights into the genetic engineering of the high use efficiency of N in forestry trees.

Genome-Wide Identification and Characterization of Long Noncoding RNAs in Populus × canescens Roots Treated With Different Nitrogen Fertilizers

Nitrate (NO₃ ^-) and ammonium (NH₄ ⁺) are the primary forms of inorganic nitrogen acquired by plant roots. LncRNAs, as key regulators of gene expression, are a class of non-coding RNAs larger than 200 bp. However, knowledge about the regulatory role of lncRNAs in response to different nitrogen forms remains limited, particularly in woody plants. Here, we performed strand-specific RNA-sequencing of P. × canescens roots under three different nitrogen fertilization treatments. In total, 324 lncRNAs and 6,112 mRNAs were identified as showing significantly differential expression between the NO₃ ^- and NH₄NO₃ treatments. Moreover, 333 lncRNAs and 6,007 mRNAs showed significantly differential expression between the NH₄ ⁺ and NH₄NO₃ treatments. Further analysis suggested that these lncRNAs and mRNAs have different response mechanisms for different nitrogen forms. In addition, functional annotation of cis and trans target mRNAs of differentially expressed lncRNAs indicated that 60 lncRNAs corresponding to 49 differentially expressed cis and trans target mRNAs were involved in plant nitrogen metabolism and amino acid biosynthesis and metabolism. Furthermore, 42 lncRNAs were identified as putative precursors of 63 miRNAs, and 28 differentially expressed lncRNAs were potential endogenous target mimics targeted by 96 miRNAs. Moreover, ceRNA regulation networks were constructed. MSTRG.6097.1, MSTRG.13550.1, MSTRG.2693.1, and MSTRG.12899.1, as hub lncRNAs in the ceRNA networks, are potential candidate lncRNAs for studying the regulatory mechanism in poplar roots under different nitrogen fertilization treatments. The results provide a basis for obtaining insight into the molecular mechanisms of lncRNA responses to different nitrogen forms in woody plants.

Identification of Genes Associated with Nitrogen Stress Responses in Apple Leaves

Nitrogen (N) is an essential macronutrient that regulates diverse physiological processes for plant survival and development. In apple orchards, inappropriate N conditions can cause imbalanced growth and subsequent physiological disorders in trees. In order to investigate the molecular basis underlying the physiological signals for N stress responses, we examined the metabolic signals responsive to contrasting N stress conditions (deficient/excessive) in apple leaves using transcriptome approaches. The clustering of differentially expressed genes (DEGs) showed the expression dynamics of genes associated with each N stress group. Functional analyses of gene ontology and pathway enrichments revealed the potential candidates of metabolic signals responsible for N-deficient/excessive stress responses. The functional interactions of DEGs in each cluster were further explored by protein-protein interaction network analysis. Our results provided a comprehensive insight into molecular signals responsive to N stress conditions, and will be useful in future research to enhance the nutrition tolerance of tree crops.

miRNAomic Approach to Plant Nitrogen Starvation

Nitrogen (N) is one of the indispensable nutrients required by plants for their growth, development, and survival. Being a limited nutrient, it is mostly supplied exogenously to the plants, to maintain quality and productivity. The increased use of N fertilizers is associated with high-cost inputs and negative environmental consequences, which necessitates the development of nitrogen-use-efficient plants for sustainable agriculture. Understanding the regulatory mechanisms underlying N metabolism in plants under low N is one of the prerequisites for the development of nitrogen-use-efficient plants. One of the important and recently discovered groups of regulatory molecules acting at the posttranscriptional and translational levels are microRNAs (miRNAs). miRNAs are known to play critical roles in the regulation of gene expression in plants under different stress conditions including N stress. Several classes of miRNAs associated with N metabolism have been identified so far. These nitrogen-responsive miRNAs may provide a platform for a better understanding of the regulation of N metabolism and pave a way for the development of genotypes for better N utilization. The current review presents a brief outline of miRNAs and their regulatory role in N metabolism.

MicroRNAs in Woody Plants

MicroRNAs (miRNAs) are small (∼21-nucleotides) non-coding RNAs found in plant and animals. MiRNAs function as critical post-transcriptional regulators of gene expression by binding to complementary sequences in their target mRNAs, leading to mRNA destabilization and translational inhibition. Plant miRNAs have some distinct characteristics compared to their animal counterparts, including greater evolutionary conservation and unique miRNA processing methods. The lifecycle of a plant begins with embryogenesis and progresses through seed germination, vegetative growth, reproductive growth, flowering and fruiting, and finally senescence and death. MiRNAs participate in the transformation of plant growth and development and directly monitor progression of these processes and the expression of certain morphological characteristics by regulating transcription factor genes involved in cell growth and differentiation. In woody plants, a large and rapidly increasing number of miRNAs have been identified, but their biological functions are largely unknown. In this review, we summarize the progress of miRNA research in woody plants to date. In particular, we discuss the potential roles of these miRNAs in growth, development, and biotic and abiotic stresses responses in woody plants.

Physiological characteristics and miRNA sequencing of two root zones with contrasting ammonium assimilation patterns in Populus

BACKGROUND

The net ammonium fluxes differ among the different root zones of Populus, but the physiological and microRNA regulatory mechanisms are unclear.

OBJECTIVE

To elucidate the physiological and miRNA regulatory mechanisms, we investigated the two root zones displaying significant differences in net NH₄⁺ effluxes of P. × canescens.

METHODS

Populus plantlets were cultivated with 500 μM NH₄Cl for 10 days. Six plants were randomly selected to determine the net NH₄⁺ fluxes using a noninvasive microtest technique. High-throughput sequencing were used to determine the dynamic expression profile of miRNA among the different root zones of Populus.

RESULTS

Net NH₄⁺ efflux in zone I (from 0 to 40 mm from the root apex) was - 19.64 pmol cm^-2 s^-1 and in zone II (from 40 to 80 mm) it was - 43.96 pmol cm^-2 s^-1. The expression of eleven miRNAs was significantly upregulated, whereas fifteen miRNAs were downregulated. Moreover, eighty-eight target genes of the significantly differentially expressed miRNAs were identified in root zone II compared with zone I. Particularly, ptc-miR171a/b/e and their target, SCL6, were found to be important for the difference in net NH₄⁺ effluxes in the two root zones. Moreover, the expression of the target of ptc-miR169d, NFYA3 was upregulated in root zone II compared with root zone I, contributing to increased NH₄⁺ efflux and decreased NH₄⁺ assimilation in root zone II.

CONCLUSION

These results indicate that miRNAs regulate the expression levels of their target genes and thus play key roles in net NH₄⁺ fluxes and NH₄⁺ assimilation in different poplar root zones.

The urea transporter DUR3 is differentially regulated by abiotic and biotic stresses in coffee plants

The high costs of N fertilizers in the coffee production emphasizes the need to optimize fertilization practices and improve nitrogen use efficiency. Urea is widespread in nature, characterizing itself as a significant source of nitrogen for the growth and development of several organisms. Thus, the characterization of genes involved in urea transport in coffee plants is an important research topic for the sustainable production of this valuable cash crop. In the current study, we evaluated the expression of the DUR3 gene under abiotic and biotic stresses in coffee plants. Here, we show that the expression of a high-affinity urea transporter gene (CaDUR3) was up-regulated by N starvation in leaves and roots of two out of three C. arabica cultivars examined. Moreover, the CaDUR3 gene was differentially expressed in coffee plants under different abiotic and biotic stresses. In plants of cv. IAPAR59, CaDUR3 showed an increased expression in leaves after exposure to water deficit and heat stress, while it was downregulated in plants under salinity. Upon infection with H. vastatrix (coffee rust), the CaDUR3 was markedly up-regulated at the beginning of the infection process in the disease susceptible Catuaí Vermelho 99 in comparison with the resistant cultivar. These results indicate that besides urea acquisition and N-remobilization, CaDUR3 gene may be closely involved in the response to various stresses.

Genome-Wide Screening and Characterization of Non-Coding RNAs in Coffea canephora

Coffea canephora grains are highly traded commodities worldwide. Non-coding RNAs (ncRNAs) are transcriptional products involved in genome regulation, environmental responses, and plant development. There is not an extensive genome-wide analysis that uncovers the ncRNA portion of the C. canephora genome. This study aimed to provide a curated characterization of six ncRNA classes in the Coffea canephora genome. For this purpose, we employed a combination of similarity-based and structural-based computational approaches with stringent curation. Candidate ncRNA loci had expression evidence analyzed using sRNA-seq libraries. We identified 7455 ncRNA loci (6976 with transcriptional evidence) in the C. canephora genome. This comprised of total 115 snRNAs, 1031 snoRNAs, 92 miRNA precursors, 602 tRNAs, 72 rRNAs, and 5064 lncRNAs. For miRNAs, we identified 159 putative high-confidence targets. This study was the most extensive genomic catalog of curated ncRNAs in the Coffea genus. This data might help elaborating more robust hypotheses in future comparative genomic studies as well as gene regulation and genome dynamics, helping to understand the molecular basis of domestication, environmental adaptation, resistance to pests and diseases, and coffee productivity.

Transcriptional patterns of Coffea arabica L. nitrate reductase, glutamine and asparagine synthetase genes are modulated under nitrogen suppression and coffee leaf rust

This study evaluated the transcriptional profile of genes related to nitrogen (N) assimilation in coffee plants susceptible and resistant to rust fungi under N sufficiency and N suppression. For this purpose, we inoculated young coffee leaves with Hemileia vastatrix uredospores and collected them at 0, 12, 24 and 48 hours post-inoculation (HPI) to evaluate the relative expressions of genes encoding cytosolic glutamine synthetase (CaGS₁), plastid glutamine synthetase (CaGS₂), nitrate reductase (CaNR), and asparagine synthetase (CaAS). The genes exhibited distinct patterns of transcriptional modulation for the different genotypes and N nutritional regimes. The resistant genotype (I59) presented high levels of transcription in response to pathogen inoculation for CaNR and CaGS₁ genes, evaluated under N sufficiency in the initial moments of infection (12 HPI). The gene CaGS₁ also showed a peak at 48 HPI. The susceptible genotype (CV99) showed increased transcript rates of CaNR at 12 and 24 HPI in response to rust inoculation. The transcriptional patterns observed for CV99, under N suppression, were high levels for CaAS and CaGS₂ at all post-inoculation times in response to coffee leaf rust disease. In addition, CaGS₁ was up-regulated at 48 HPI for CV99. Cultivar I59 showed high transcript levels at 12 HPI for CaAS and peaks at 24 and 48 HPI for CaGS₂ in inoculated samples. Consequently, total chlorophyl concentration was influenced by N suppression and by rust infection. Regarding enzyme activities in vitro for glutamine synthetase and CaNR, there was an increase in infected coffee leaves (I59) and under N sufficiency. Moreover, CV99 was modulated in both N nutritional regimes for GS activity in response to rust. Our results indicate that N transport genes trigger a differential modulation between genotypes through the action of rust disease.

A Small RNA-Mediated Regulatory Network in Arabidopsis thaliana Demonstrates Connectivity Between phasiRNA Regulatory Modules and Extensive Co-Regulation of Transcription by miRNAs and phasiRNAs

Gene regulation involves the orchestrated action of multiple regulators to fine-tune the expression of genes. Hierarchical interactions and co-regulation among regulators are commonly observed in biological systems, leading to complex regulatory networks. Small RNA (sRNAs) have been shown to be important regulators of gene expression due to their involvement in multiple cellular processes. In plants, microRNA (miRNAs) and phased small interfering RNAs (phasiRNAs) correspond to two well-characterized types of sRNAs involved in the regulation of posttranscriptional gene expression, although information about their targets and interactions with other gene expression regulators is limited. We describe an extended sRNA-mediated regulatory network in Arabidopsis thaliana that provides a reference frame to understand sRNA biogenesis and activity at the genome-wide level. This regulatory network combines a comprehensive evaluation of phasiRNA production and sRNA targets supported by degradome data. The network includes ~17% of genes in the A. thaliana genome, representing ~50% annotated gene ontology (GO) functional categories. Approximately 14% of genes with GO annotations corresponding to regulation of gene expression were found to be under sRNA control. The unbiased bioinformatic approach used to produce the network was able to detect 107 PHAS loci (regions of phasiRNA production), 5,047 active phasiRNAs (~70% of which were non-canonical), and reconstruct 17 regulatory modules resulting from complex regulatory interactions between different sRNA-regulatory pathways. Known regulatory modules like miR173-TAS-PPR/TPR and miR390-TAS3-ARF/F-box were faithfully reconstructed and expanded, illustrating the accuracy and sensitivity of the methods and providing confidence for the validity of findings of previously unrecognized modules. The network presented here includes a 2X increase in the number of identified PHAS loci, a large complement (~70%) of non-canonical phasiRNAs, and the most comprehensive evaluation of sRNA cleavage activity in A. thaliana to date. Structural analysis showed similarities to networks of other biological systems and demonstrated connectivity between phasiRNA regulatory modules with extensive co-regulation of transcripts by miRNAs and phasiRNAs. The described regulatory network provides a reference that will facilitate global analyses of individual plant regulatory programs such as those that control homeostasis, development, and responses to biotic and abiotic environmental changes.