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

Comparative genome-wide characterization of salt responsive micro RNA and their targets through integrated small RNA and de novo transcriptome profiling in sugarcane and its wild relative Erianthus arundinaceus

Soil salinity and saline irrigation water are major constraints in sugarcane affecting the production of cane and sugar yield. To understand the salinity induced responses and to identify novel genomic resources, integrated de novo transcriptome and small RNA sequencing in sugarcane wild relative, Erianthus arundinaceus salt tolerant accession IND 99-907 and salt-sensitive sugarcane genotype Co 97010 were performed. A total of 362 known miRNAs belonging to 62 families and 353 miRNAs belonging to 63 families were abundant in IND 99-907 and Co 97010 respectively. The miRNA families such as miR156, miR160, miR166, miR167, miR169, miR171, miR395, miR399, miR437 and miR5568 were the most abundant with more than ten members in both genotypes. The differential expression analysis of miRNA reveals that 221 known miRNAs belonging to 48 families and 130 known miRNAs belonging to 42 families were differentially expressed in IND 99-907 and Co 97010 respectively. A total of 12,693 and 7982 miRNA targets against the monoploid mosaic genome and a total of 15,031 and 12,152 miRNA targets against the de novo transcriptome were identified for differentially expressed known miRNAs of IND 99-907 and Co 97010 respectively. The gene ontology (GO) enrichment analysis of the miRNA targets revealed that 24, 12 and 14 enriched GO terms (FDR < 0.05) for biological process, molecular function and cellular component respectively. These miRNAs have many targets that associated in regulation of biotic and abiotic stresses. Thus, the genomic resources generated through this study are useful for sugarcane crop improvement through biotechnological and advanced breeding approaches.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03867-7.

Analysis of miRNAs responsive to long-term calcium deficiency in tef (Eragrostis tef (Zucc.) Trotter)

MicroRNAs (miRNAs) play an important role in growth, development, stress resilience, and epigenetic modifications of plants. However, the effect of calcium (Ca²⁺) deficiency on miRNA expression in the orphan crop tef (Eragrostis tef) remains unknown. In this study, we analyzed expression of miRNAs in roots and shoots of tef in response to Ca²⁺ treatment. miRNA-seq followed by bioinformatic analysis allowed us to identify a large number of small RNAs (sRNAs) ranging from 17 to 35 nt in length. A total of 1380 miRNAs were identified in tef experiencing long-term Ca²⁺ deficiency while 1495 miRNAs were detected in control plants. Among the miRNAs identified in this study, 161 miRNAs were similar with those previously characterized in other plant species and 348 miRNAs were novel, while the remaining miRNAs were uncharacterized. Putative target genes and their functions were predicted for all the known and novel miRNAs that we identified. Based on gene ontology (GO) analysis, the predicted target genes are known to have various biological and molecular functions including calcium uptake and transport. Pairwise comparison of differentially expressed miRNAs revealed that some miRNAs were specifically enriched in roots or shoots of low Ca²⁺-treated plants. Further characterization of the miRNAs and their targets identified in this study may help in understanding Ca²⁺ deficiency responses in tef and related orphan crops.

Genome-wide miRNAs profiles of pearl millet under contrasting high vapor pressure deficit reveal their functional roles in drought stress adaptations

Pearl millet (Pennisetum glaucum [L.] R. Br.) is an important crop capable of growing in harsh and marginal environments, with the highest degree of tolerance to drought and heat stresses among cereals. Diverse germplasm of pearl millet shows a significant phenotypic variation in response to abiotic stresses, making it a unique model to study the mechanisms responsible for stress mitigation. The present study focuses on identifying the physiological response of two pearl millet high-resolution cross (HRC) genotypes, ICMR 1122 and ICMR 1152, in response to low and high vapor pressure deficit (VPD). Under high VPD conditions, ICMR 1152 exhibited a lower transpiration rate (Tr), higher transpiration efficiency, and lower root sap exudation than ICMR 1122. Further, Pg-miRNAs expressed in the contrasting genotypes under low and high VPD conditions were identified by deep sequencing analysis. A total of 116 known and 61 novel Pg-miRNAs were identified from ICMR 1152, while 26 known and six novel Pg-miRNAs were identified from ICMR 1122 genotypes, respectively. While Pg-miR165, 168, 170, and 319 families exhibited significant differential expression under low and high VPD conditions in both genotypes, ICMR 1152 showed abundant expression of Pg-miR167, Pg-miR172, Pg-miR396 Pg-miR399, Pg-miR862, Pg-miR868, Pg-miR950, Pg-miR5054, and Pg-miR7527 indicating their direct and indirect role in root physiology and abiotic stress responses. Drought responsive Pg-miRNA targets showed upregulation in response to high VPD stress, further narrowing down the miRNAs involved in regulation of drought tolerance in pearl millet.

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.

Noncoding-RNA-Mediated Regulation in Response to Macronutrient Stress in Plants

Macronutrient elements including nitrogen (N), phosphorus (P), potassium (K), calcium (Ca), magnesium (Mg), and sulfur (S) are required in relatively large and steady amounts for plant growth and development. Deficient or excessive supply of macronutrients from external environments may trigger a series of plant responses at phenotypic and molecular levels during the entire life cycle. Among the intertwined molecular networks underlying plant responses to macronutrient stress, noncoding RNAs (ncRNAs), mainly microRNAs (miRNAs) and long ncRNAs (lncRNAs), may serve as pivotal regulators for the coordination between nutrient supply and plant demand, while the responsive ncRNA-target module and the interactive mechanism vary among elements and species. Towards a comprehensive identification and functional characterization of nutrient-responsive ncRNAs and their downstream molecules, high-throughput sequencing has produced massive omics data for comparative expression profiling as a first step. In this review, we highlight the recent findings of ncRNA-mediated regulation in response to macronutrient stress, with special emphasis on the large-scale sequencing efforts for screening out candidate nutrient-responsive ncRNAs in plants, and discuss potential improvements in theoretical study to provide better guidance for crop breeding practices.

Genetics of Germination and Seedling Traits under Drought Stress in a MAGIC Population of Maize

Drought is one of the most detrimental abiotic stresses hampering seed germination, development, and productivity. Maize is more sensitive to drought than other cereals, especially at seedling stage. Our objective was to study genetic regulation of drought tolerance at germination and during seedling growth in maize. We evaluated 420 RIL with their parents from a multi-parent advanced generation inter-cross (MAGIC) population with PEG-induced drought at germination and seedling establishment. A genome-wide association study (GWAS) was carried out to identify genomic regions associated with drought tolerance. GWAS identified 28 and 16 SNPs significantly associated with germination and seedling traits under stress and well-watered conditions, respectively. Among the SNPs detected, two SNPs had significant associations with several traits with high positive correlations, suggesting a pleiotropic genetic control. Other SNPs were located in regions that harbored major QTLs in previous studies, and co-located with QTLs for cold tolerance previously published for this MAGIC population. The genomic regions comprised several candidate genes related to stresses and plant development. These included numerous drought-responsive genes and transcription factors implicated in germination, seedling traits, and drought tolerance. The current analyses provide information and tools for subsequent studies and breeding programs for improving drought tolerance.

Deciphering the transcriptomic insight during organogenesis in Castor (Ricinus communis L.), Jatropha (Jatropha curcas L.) and Sunflower (Helianthus annuus L.)

Cultivation of the castor crop is hindered by various factors and one of the approaches for genetic improvement of the crop is through exploitation of biotechnological tools. Response of castor tissues to in vitro culture is poor which necessitated this study on understanding the molecular basis of organogenesis in cultured tissues of castor, through de novo transcriptome analysis and by comparing with jatropha and sunflower having good regeneration ability. Transcriptome profiling analysis was carried out with hypocotyl explants from castor, jatropha and cotyledons from sunflower cultured on MS media supplemented with different concentrations of hormones. Differentially expressed genes during dedifferentiation and organogenic differentiation stages of callus included components of auxin and cytokinin signaling, secondary metabolite synthesis, genes encoding transcription factors, receptor kinases and protein kinases. In castor, many genes involved in auxin biosynthesis and homeostasis like WAT1, vacuolar transporter genes, transcription factors like short root like protein were down-regulated while genes like DELLA were up-regulated accounting for regeneration recalcitrance. Validation of 62 DEGs through qRT-PCR showed a consensus of 77.4% of the genes expressed. Overall study provides set of genes involved in the process of organogenesis in three oilseed crops which forms a basis for understanding and improving the efficiency of plant regeneration and genetic transformation in castor.

Whole Genome Sequencing and Comparative Genomic Analysis Reveal Allelic Variations Unique to a Purple Colored Rice Landrace (Oryza sativa ssp. indica cv. Purpleputtu)

Purpleputtu (Oryza sativa ssp. indica cv. Purpleputtu) is a unique rice landrace from southern India that exhibits predominantly purple color. This study reports the underlying genetic complexity of the trait, associated domestication and de-domestication processes during its coevolution with present day cultivars. Along-with genome level allelic variations in the entire gene repertoire associated with the purple, red coloration of grain and other plant parts. Comparative genomic analysis using 'a panel of 108 rice lines' revealed a total of 3,200,951 variants including 67,774 unique variations in Purpleputtu (PP) genome. Multiple sequence alignment uncovered a 14 bp deletion in Rc (Red colored, a transcription factor of bHLH class) locus of PP, a key regulatory gene of anthocyanin biosynthetic pathway. Interestingly, this deletion in Rc gene is a characteristic feature of the present-day white pericarped rice cultivars. Phylogenetic analysis of Rc locus revealed a distinct clade showing proximity to the progenitor species Oryza rufipogon and O. nivara. In addition, PP genome exhibits a well conserved 4.5 Mbp region on chromosome 5 that harbors several loci associated with domestication of rice. Further, PP showed 1,387 unique when SNPs compared to 3,023 lines of rice (SNP-Seek database). The results indicate that PP genome is rich in allelic diversity and can serve as an excellent resource for rice breeding for a variety of agronomically important traits such as disease resistance, enhanced nutritional values, stress tolerance, and protection from harmful UV-B rays.

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.

Deciphering genome-wide WRKY gene family of Triticum aestivum L. and their functional role in response to Abiotic stress

WRKY transcription factors (TFs) act in regulating plant growth and development as well as in response to different stress. Some earlier studies done by individual researchers reported different wheat WRKY TFs. Although, the recently released wheat genome has opened an avenue to investigate wheat WRKYs (TaWRKY) TFs. Prime objective of this study to performed genome-wide classifications of TaWRKYs and their functional annotation. The classification of 107 individual identified characterized sequences of TaWRKY (IICS-TaWRKY) and 160 uncharacterized draft sequences of TaWRKY (UDS-TaWRKY), along with their gene structures and motifs analysis was performed. Along with comparative sequence analysis and microarray analysis was performed to mimic out TaWRKYs functions in response to different abiotic stresses, accompanied by in-vitro validation. The comparative phylogenetic analysis and estimation of Ka/Ks ratio with Triticum urartu, illustrate group based clasifications of TaWRKYs and evolutionary divergences. Furthermore, motif-based and protein-DNA interaction analysis of TaWRKYs helps to identify, their putative function in target DNA recognition sites. Subsequently, results of microarray and comparative sequence analysis provides the evidence of TaWRKYs involved in heat and/or drought stress. Further, in-vitro results validates that TaWRKY014, TaWRKY090 are found to participate in response of drought stress, whereas TaWRKY008, TaWRKY122, and WRKY45 are involved in response of heat and drought stress. These findings can be utilized in developing novel heat and drought-tolerant wheat cultivars using marker-assisted breeding and transgenic development.

Extrapolation of significant genes and transcriptional regulatory networks involved in Zea mays in response in UV-B stress

A wide range of plant species growth influenced when they exposed to solar UV-B radiation. Leaves of the plant are highly affected by UV-B radiation lead to the reduction in the growth of the plant. Current work demonstrates the comparative transcriptional changes and visible symptoms occurred in the maize leaf growth zone (GZ). Primary objective of this study was to identify differentially expressed genes (DEGs) responsible for leaf growth and their association in the transcriptional regulatory network under UV-B stress. Whole transcriptomic data was analysed and the quality check was tested for each sample and further genome-wide mapping and DEGs were performed. Gene Ontology (GO) based functional annotation, associated transcriptional networks and molecular pathways were annotated. Reduction in cell production due to UV-B stress causes a decrease in leaf's length and size was observed. Further, the specific role of the DEGs, in UV-B signalling pathways and other molecular functions responsible for leaf cell death was discovered. Results also infer that the major changes occurred in the cell cycle, transcriptional regulation, post-transcriptional modification, phytohormones, flavonoids biosynthesis, and chromatin remodeling. UV-B signalling pathways and the transcriptional regulatory networks infer the different molecular steps along with downstream transcriptional and post-transcriptional control of metabolic enzymes used in long-term memory adoption and attainment resistance to UV-B stress identified. Effects of UV-B radiation on leaf growth was noted in this study. UV-B stress response genes and associated transcriptional regulatory networks were identified, can be used in developing the marker assist UB-B stress tolerant genotypes of the maize.

Analysis of physiological and miRNA responses to Pi deficiency in alfalfa (Medicago sativa L.)

The induction of miR399 and miR398 and the inhibition of miR156, miR159, miR160, miR171, miR2111, and miR2643 were observed under Pi deficiency in alfalfa. The miRNA-mediated genes involved in basic metabolic process, root and shoot development, stress response and Pi uptake. Inorganic phosphate (Pi) deficiency is known to be a limiting factor in plant development and growth. However, the underlying miRNAs associated with the Pi deficiency-responsive mechanism in alfalfa are unclear. To elucidate the molecular mechanism at the miRNA level, we constructed four small RNA (sRNA) libraries from the roots and shoots of alfalfa grown under normal or Pi-deficient conditions. In the present study, alfalfa plants showed reductions in biomass, photosynthesis, and Pi content and increases in their root-to-shoot ratio and citric, malic, and succinic acid contents under Pi limitation. Sequencing results identified 47 and 44 differentially expressed miRNAs in the roots and shoots, respectively. Furthermore, 909 potential target genes were predicted, and some targets were validated by RLM-RACE assays. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed prominent enrichment in signal transducer activity, binding and basic metabolic pathways for carbohydrates, fatty acids and amino acids; cellular response to hormone stimulus and response to auxin pathways were also enriched. qPCR results verified that the differentially expressed miRNA profile was consistent with sequencing results, and putative target genes exhibited opposite expression patterns. In this study, the miRNAs associated with the response to Pi limitation in alfalfa were identified. In addition, there was an enrichment of miRNA-targeted genes involved in biological regulatory processes such as basic metabolic pathways, root and shoot development, stress response, Pi transportation and citric acid secretion.

High-Throughput Sequencing of Small RNA Transcriptomes in Maize Kernel Identifies miRNAs Involved in Embryo and Endosperm Development

Maize kernel development is a complex biological process that involves the temporal and spatial expression of many genes and fine gene regulation at a transcriptional and post-transcriptional level, and microRNAs (miRNAs) play vital roles during this process. To gain insight into miRNA-mediated regulation of maize kernel development, a deep-sequencing technique was used to investigate the dynamic expression of miRNAs in the embryo and endosperm at three developmental stages in B73. By miRNA transcriptomic analysis, we characterized 132 known miRNAs and six novel miRNAs in developing maize kernel, among which, 15 and 14 miRNAs were commonly differentially expressed between the embryo and endosperm at 9 days after pollination (DAP), 15 DAP and 20 DAP respectively. Conserved miRNA families such as miR159, miR160, miR166, miR390, miR319, miR528 and miR529 were highly expressed in developing embryos; miR164, miR171, miR393 and miR2118 were highly expressed in developing endosperm. Genes targeted by those highly expressed miRNAs were found to be largely related to a regulation category, including the transcription, macromolecule biosynthetic and metabolic process in the embryo as well as the vitamin biosynthetic and metabolic process in the endosperm. Quantitative reverse transcription-PCR (qRT-PCR) analysis showed that these miRNAs displayed a negative correlation with the levels of their corresponding target genes. Importantly, our findings revealed that members of the miR169 family were highly and dynamically expressed in the developing kernel, which will help to exploit new players functioning in maize kernel development.

Comparative analysis of miRNAs of two rapeseed genotypes in response to acetohydroxyacid synthase-inhibiting herbicides by high-throughput sequencing

Acetohydroxyacid synthase (AHAS), also called acetolactate synthase, is a key enzyme involved in the first step of the biosynthesis of the branched-chain amino acids valine, isoleucine and leucine. Acetohydroxyacid synthase-inhibiting herbicides (AHAS herbicides) are five chemical families of herbicides that inhibit AHAS enzymes, including imidazolinones (IMI), sulfonylureas (SU), pyrimidinylthiobenzoates, triazolinones and triazolopyrimidines. Five AHAS genes have been identified in rapeseed, but little information is available regarding the role of miRNAs in response to AHAS herbicides. In this study, an AHAS herbicides tolerant genotype and a sensitive genotype were used for miRNA comparative analysis. A total of 20 small RNA libraries were obtained of these two genotypes at three time points (0h, 24 h and 48 h) after spraying SU and IMI herbicides with two replicates. We identified 940 conserved miRNAs and 1515 novel candidate miRNAs in Brassica napus using high-throughput sequencing methods combined with computing analysis. A total of 3284 genes were predicted to be targets of these miRNAs, and their functions were shown using GO, KOG and KEGG annotations. The differentiation expression results of miRNAs showed almost twice as many differentiated miRNAs were found in tolerant genotype M342 (309 miRNAs) after SU herbicide application than in sensitive genotype N131 (164 miRNAs). In additiond 177 and 296 miRNAs defined as differentiated in sensitive genotype and tolerant genotype in response to SU herbicides. The miR398 family was observed to be associated with AHAS herbicide tolerance because their expression increased in the tolerant genotype but decreased in the sensitive genotype. Moreover, 50 novel miRNAs from 39 precursors were predicted. There were 8 conserved miRNAs, 4 novel miRNAs and 3 target genes were validated by quantitative real-time PCR experiment. This study not only provides novel insights into the miRNA content of AHAS herbicides tolerant rapeseed in response to AHAS herbicides, but also demonstrates that miRNAs may be involved in AHAS herbicides tolerance.

MicroRNAs in model and complex organisms

Non-coding RNAs such as microRNAs (miRNAs) are very tiny ribonucleotides having an essential role in gene regulation at both post-transcriptional and translational levels. They are very conserved and expressed in worms, flies, plants, and mammals in a sequence-specific manner. Furthermore, it is now possible to clone miRNAs using the new genome editing tool CRISPR/cas9, which shows benefit in control of untargeted effect. In this special issue, we tried to cover researches associated with functional roles of miRNAs accross model and complex organisms.

Identification and characterization of a grain micronutrient-related OsFRO2 rice gene ortholog from micronutrient-rich little millet (Panicum sumatrense)

Minor millets are considered as nutrient-rich cereals having significant effect in improving human health. In this study, a rice ortholog of Ferric Chelate Reductase (FRO2) gene involved in plant metal uptake has been identified in iron-rich Little millet (LM) using PCR and next generation sequencing-based strategy. FRO2 gene-specific primers designed from rice genome amplified 2.7 Kb fragment in LM genotype RLM-37. Computational genomics analyses of the sequenced amplicon showed high level sequence similarity with rice OsFRO2 gene. The predicted gene structure showed the presence of 6 exons and 5 introns and its protein sequence was found to contain ferric reductase and NOX_Duox_Like_FAD_NADP domains. Further, 3D structure analysis of FCR-LM model protein (494 amino acids) shows that it has 18 helices, 10 beta sheets, 10 strands, 41 beta turn and 5 gamma turn with slight deviation from the FCR-Os structure. Besides, the structures of FCR-LM and FCR-Os were modelled followed by molecular dynamics simulations. The overall study revealed both sequence and structural similarity between the identified gene and OsFRO2. Thus, a putative ferric chelate reductase gene has been identified in LM paving the way for using this approach for identification of orthologs of other metal genes from millets. This also facilitates mining of effective alleles of known genes for improvement of staple crops like rice.

Transcriptomic Analysis of Soil Grown T. aestivum cv. Root to Reveal the Changes in Expression of Genes in Response to Multiple Nutrients Deficiency

Deficiency of necessary macronutrients, i.e., Potassium (K), Magnesium (Mg), Nitrogen (N), Phosphorus (P), and Sulfate (S) in the soil leads to a reduction in plant growth and yield, which is a result of changes in expression level of various genes. This study was performed to identify the differentially expressed genes and its associated metabolic pathways occurred in soil grown wheat root samples excavated from the control and treated fields. To identify the difference in gene expression levels due to deficiency of the said nutrients, a transcriptomic, meta-analysis was performed on array expression profile data. A set of 435 statistically significant probes encoding 398 Nutrient Deficiency Response Genes (NRGs) responding at-least one nutrients deficiency (ND) were identified. Out of them 55 NRGs were found to response to minimum two ND. Singular Enrichment Analysis (SEA) predicts ontological based classifications and functional analysis of NRGs in different cellular/molecular pathways involved in root development and growth. Functional annotation and reaction mechanism of differentially expressed genes, proteins/enzymes in the different metabolic pathway through MapMan analysis were explored. Further the meta-analysis was performed to revels the active involvement each NRGs in distinct tissues and their comparative potential expression analysis in different stress conditions. The study results in exploring the role of major acting candidate genes such as Non-specific serine/threonine protein kinase, Xyloglucan endotransglucosylase/hydrolase, Peroxides, Glycerophosphoryl diester phosphodiesterase, S-adenosylmethionine decarboxylase proenzyme, Dehydrin family proteins, Transcription factors, Membrane Proteins, Metal binding proteins, Photosystem proteins, Transporter and Transferase associated in different metabolic pathways. Finally, the differences of transcriptional responses in the soil-grown root of T. aestivum cv. and in-vitro grown model plants under nutrients deficiency were summarized.