Bioinformatic identification and experimental validation of miRNAs from foxtail millet (Setaria italica)

Realizing visionary goals for the International Year of Millet (IYoM): accelerating interventions through advances in molecular breeding and multiomics resources

MAIN CONCLUSION

Leveraging advanced breeding and multi-omics resources is vital to position millet as an essential "nutricereal resource," aligning with IYoM goals, alleviating strain on global cereal production, boosting resilience to climate change, and advancing sustainable crop improvement and biodiversity. The global challenges of food security, nutrition, climate change, and agrarian sustainability demand the adoption of climate-resilient, nutrient-rich crops to support a growing population amidst shifting environmental conditions. Millets, also referred to as "Shree Anna," emerge as a promising solution to address these issues by bolstering food production, improving nutrient security, and fostering biodiversity conservation. Their resilience to harsh environments, nutritional density, cultural significance, and potential to enhance dietary quality index made them valuable assets in global agriculture. Recognizing their pivotal role, the United Nations designated 2023 as the "International Year of Millets (IYoM 2023)," emphasizing their contribution to climate-resilient agriculture and nutritional enhancement. Scientific progress has invigorated efforts to enhance millet production through genetic and genomic interventions, yielding a wealth of advanced molecular breeding technologies and multi-omics resources. These advancements offer opportunities to tackle prevailing challenges in millet, such as anti-nutritional factors, sensory acceptability issues, toxin contamination, and ancillary crop improvements. This review provides a comprehensive overview of molecular breeding and multi-omics resources for nine major millet species, focusing on their potential impact within the framework of IYoM. These resources include whole and pan-genome, elucidating adaptive responses to abiotic stressors, organelle-based studies revealing evolutionary resilience, markers linked to desirable traits for efficient breeding, QTL analysis facilitating trait selection, functional gene discovery for biotechnological interventions, regulatory ncRNAs for trait modulation, web-based platforms for stakeholder communication, tissue culture techniques for genetic modification, and integrated omics approaches enabled by precise application of CRISPR/Cas9 technology. Aligning these resources with the seven thematic areas outlined by IYoM catalyzes transformative changes in millet production and utilization, thereby contributing to global food security, sustainable agriculture, and enhanced nutritional consequences.

Integrative genomics analysis of the ever-shrinking pectin methylesterase (PME) gene family in foxtail millet (Setaria italica)

Pectin methylesterase (PME) plays a vital role in the growth and development of plants. Their genes can be classified into two types, with Type-1 having an extra domain, PMEI. PME genes in foxtail millet (Setaria italica L.) have not been identified, and their sequence features and evolution have not been explored. Here, we identified 41 foxtail millet PME genes. Decoding the pro-region, containing the PMEI domain, revealed its more active nature than the DNA encoding PME domain, easier to be lost to produce Type-2 PME genes. We inferred that the active nature of the pro-region could be related to its harbouring more repetitive DNA sequences. Further, we revealed that though whole-genome duplication and tandem duplication contributed to producing new copies of PME genes, phylogenetic analysis provided clear evidence of ever-shrinking gene family size in foxtail millet and the other grasses in the past 100 million years. Phylogenetic analysis also supports the existence of two gene groups, Group I and Group II, with genes in Group II being more conservative. Our research contributes to understanding how DNA sequence structure affects the functional innovation and evolution of PME genes.

The Construction and Exploration of a Comprehensive MicroRNA Centered Regulatory Network in Foxtail Millet (Setaria italica L.)

MicroRNA (miRNA) is an essential endogenous post-transcriptional regulatory factor, and foxtail millet (Setaria italica L.) is an ideal C4 model cereal that is a highly valuable crop in semiarid and arid areas. The Research on comprehensive and high confidence identification and annotation of foxtail millet miRNAs needs to be strengthened, and to our knowledge, there is no information on the regulatory network of foxtail millet miRNA. In this study, 136 high confidence miRNAs were identified through high-throughput sequencing of the small RNAs in seven tissues at the shooting and grain filling stages of foxtail millet. A total of 2,417 target genes were obtained by combining computational biology software and degradome sequencing methods. Furthermore, an analysis using transcriptome sequencing revealed the relationships between miRNAs and their target genes and simultaneously explored key regulatory modules in panicles during the grain filling stage. An miRNA regulatory network was constructed to explore the functions of miRNA in more detail. This network, centered on miRNAs and combining upstream transcriptional factors and downstream target genes, is primarily composed of feed forward loop motifs, which greatly enhances our knowledge of the potential functions of miRNAs and uncovers numerous previously unknown regulatory links. This study provides a solid foundation for research on the function and regulatory network of miRNAs in foxtail millet.

Identification, Molecular Characteristics, and Evolution of GRF Gene Family in Foxtail Millet (Setaria italica L.)

Growth-regulating factor (GRF) is a multigene family that plays a vital role in the growth and development of plants. In the past, the GRF family of many plants has been studied. However, there is not a report about identification and evolution of GRF in foxtail millet (Setaria italia). Here, we identified 10 GRF genes in foxtail millet. Seven (70.00%) were regulated by Sit-miR396, and there were 19 optimal codons in GRFs of foxtail millet. Additionally, we found that WGD or segmental duplication have affected GRFs in foxtail millet between 15.07 and 45.97 million years ago. Regarding the GRF gene family of land plants, we identified a total of 157 GRF genes in 15 representative land plants. We found that GRF gene family originated from Group E, and the GRF gene family in monocots was gradually shrinking. Also, more loss resulted from the small number of GRF genes in lower plants. Exploring the evolution of GRF and functional analysis in the foxtail millet help us to understand GRF better and make a further study about the mechanism of GRF. These results provide a basis for the genetic improvement of foxtail millet and indicate an improvement of the yield.

Azadirachta indica MicroRNAs: Genome-Wide Identification, Target Transcript Prediction, and Expression Analyses

MicroRNAs are short, endogenous, non-coding RNAs, liable for essential regulatory function. Numerous miRNAs have been identified and studied in plants with known genomic or small RNA resources. Despite the availability of genomic and transcriptomic resources, the miRNAs have not been reported in the medicinal tree Azadirachta indica (Neem) till date. Here for the first time, we report extensive identification of miRNAs and their possible targets in A. indica which might help to unravel their therapeutic potential. A comprehensive search of miRNAs in the A. indica genome by C-mii tool was performed. Overall, 123 miRNAs classified into 63 families and their stem-loop hairpin structures were predicted. The size of the A. indica (ain)-miRNAs ranged between 19 and 23 nt in length, and their corresponding ain-miRNA precursor sequence MFEI value averaged as -1.147 kcal/mol. The targets of ain-miRNAs were predicted in A. indica as well as Arabidopsis thaliana plant. The gene ontology (GO) annotation revealed the involvement of ain-miRNA targets in developmental processes, transport, stress, and metabolic processes including secondary metabolism. Stem-loop qRT-PCR was carried out for 25 randomly selected ain-miRNAs and differential expression patterns were observed in different A. indica tissues. Expression of miRNAs and its targets shows negative correlation in a dependent manner.

Design of novel viral attachment inhibitors of the spike glycoprotein (S) of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) through virtual screening and dynamics

To date, the global COVID-19 pandemic has been associated with 11.8 million cases and over 545481 deaths. In this study, we have employed virtual screening approaches and selected 415 lead-like compounds from 103 million chemical substances, based on the existing drugs, from PubChem databases as potential candidates for the S protein-mediated viral attachment inhibition. Thereafter, based on drug-likeness and Lipinski's rules, 44 lead-like compounds were docked within the active side pocket of the viral-host attachment site of the S protein. Corresponding ligand properties and absorption, distribution, metabolism, excretion, and toxicity (ADMET) profile were measured. Furthermore, four novel inhibitors were designed and assessed computationally for efficacy. Comparative analysis showed the screened compounds in this study maintain better results than the proposed mother compounds, VE607 and SSAA09E2. The four designed novel lead compounds possessed more fascinating output without deviating from any of Lipinski's rules. They also showed higher bioavailability and the drug-likeness score was 0.56 and 1.81 compared with VE607 and SSAA09E2, respectively. All the screened compounds and novel compounds showed promising ADMET properties. Among them, the compound AMTM-02 was the best candidate, with a docking score of -7.5 kcal/mol. Furthermore, the binding study was verified by molecular dynamics simulation over 100 ns by assessing the stability of the complex. The proposed screened compounds and the novel compounds may give some breakthroughs for the development of a therapeutic drug to treat SARS-CoV-2 proficiently in vitro and in vivo.

Identification and characterization of Populus microRNAs in response to plant growth-promoting endophytic Streptomyces sp. SSD49

Endophytic Streptomyces sp. SSD49 inhibited eight pathogens, including the human opportunistic pathogenic microorganisms, the plant pathogenic fungi and bacteria. The growth of soybeans, tomatoes, peppers and Populus tomentosa seedings inoculated with SSD49 are remarkably promoted. Here, we constructed two P. tomentosa seedling microRNA (miRNA) libraries inoculated with (PS30d) and without SSD49 (PC30d) to explore the molecular regulatory roles in the plant response to the beneficial bacteria. Totals of 314 known and 144 novel miRNAs were identified, among which 27 known and 11 novel miRNA had significantly different expression. The targets of up-regulated miR160, miR156, ptc114 and down-regulated miR319 and other differential expressed miRNAs primarily regulated genes encoding transcription factors (auxin response factor, small auxin-up RNA, and GRAS proteins), disease resistance proteins, phytohormone oxidase, and response regulators, which could promote plant growth, influence disease resistance and miRNA biosynthesis in P. tomentosa. This is the first report on the genome-wide identification of biocontrol endophytic Streptomyces inoculation-responsive miRNAs using small RNA sequencing in P. tomentosa and these findings provide new insight into understanding the biocontrol effects of endophytic Streptomyces.

Grasses suppress shoot-borne roots to conserve water during drought

Many important crops are members of the Poaceae family, which develop root systems characterized by a high degree of root initiation from the belowground basal nodes of the shoot, termed the crown. Although this postembryonic shoot-borne root system represents the major conduit for water uptake, little is known about the effect of water availability on its development. Here we demonstrate that in the model C4 grass Setaria viridis, the crown locally senses water availability and suppresses postemergence crown root growth under a water deficit. This response was observed in field and growth room environments and in all grass species tested. Luminescence-based imaging of root systems grown in soil-like media revealed a shift in root growth from crown-derived to primary root-derived branches, suggesting that primary root-dominated architecture can be induced in S. viridis under certain stress conditions. Crown roots of Zea mays and Setaria italica, domesticated relatives of teosinte and S. viridis, respectively, show reduced sensitivity to water deficit, suggesting that this response might have been influenced by human selection. Enhanced water status of maize mutants lacking crown roots suggests that under a water deficit, stronger suppression of crown roots actually may benefit crop productivity.

Suppression Subtractive Hybridization Versus Next-Generation Sequencing in Plant Genetic Engineering: Challenges and Perspectives

Suppression subtractive hybridization (SSH) is an effective method to identify different genes with different expression levels involved in a variety of biological processes. This method has often been used to study molecular mechanisms of plants in complex relationships with different pathogens and a variety of biotic stresses. Compared to other techniques used in gene expression profiling, SSH needs relatively smaller amounts of the initial materials, with lower costs, and fewer false positives present within the results. Extraction of total RNA from plant species rich in phenolic compounds, carbohydrates, and polysaccharides that easily bind to nucleic acids through cellular mechanisms is difficult and needs to be considered. Remarkable advancement has been achieved in the next-generation sequencing (NGS) field. As a result of progress within fields related to molecular chemistry and biology as well as specialized engineering, parallelization in the sequencing reaction has exceptionally enhanced the overall read number of generated sequences per run. Currently available sequencing platforms support an earlier unparalleled view directly into complex mixes associated with RNA in addition to DNA samples. NGS technology has demonstrated the ability to sequence DNA with remarkable swiftness, therefore allowing previously unthinkable scientific accomplishments along with novel biological purposes. However, the massive amounts of data generated by NGS impose a substantial challenge with regard to data safe-keeping and analysis. This review examines some simple but vital points involved in preparing the initial material for SSH and introduces this method as well as its associated applications to detect different novel genes from different plant species. This review evaluates general concepts, basic applications, plus the probable results of NGS technology in genomics, with unique mention of feasible potential tools as well as bioinformatics.

Combined small RNA and degradome sequencing to identify miRNAs and their targets in response to drought in foxtail millet

Background

Foxtail millet (Setaria italica) is a diploid C₄ panicoid species. Because of its prominent drought resistance, small genome size, self-pollination, and short life cycle, foxtail millet has become an ideal model system for studying drought tolerance of crops. MicroRNAs (miRNAs) are endogenous, small RNAs that play important regulatory roles in the development and stress response in plants.

Results

In this study, we applied Illumina sequencing to systematically investigate the drought-responsive miRNAs derived from S. italica inbred An04-4783 seedlings grown under control and drought conditions. Degradome sequencing was applied to confirm the targets of these miRNAs at a global level. A total of 81 known miRNAs belonging to 28 families were identified, among which 14 miRNAs were upregulated and four were downregulated in response to drought. In addition, 72 potential novel miRNAs were identified, three of which were differentially expressed under drought conditions. Degradome sequencing analysis showed that 56 and 26 genes were identified as targets of known and novel miRNAs, respectively.

Conclusions

Our analysis revealed post-transcriptional remodeling of cell development, transcription factors, ABA signaling, and cellar homeostasis in S.italica in response to drought. This preliminary characterization provided useful information for further studies on the regulatory networks of drought-responsive miRNAs in foxtail millet.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0364-7) contains supplementary material, which is available to authorized users.

Dehydration-responsive miRNAs in foxtail millet: genome-wide identification, characterization and expression profiling

A set of novel and known dehydration-responsive miRNAs have been identified in foxtail millet. These findings provide new insights into understanding the functional role of miRNAs and their respective targets in regulating plant response to dehydration stress. MicroRNAs perform significant regulatory roles in growth, development and stress response of plants. Though the miRNA-mediated gene regulatory networks under dehydration stress remain largely unexplored in plant including foxtail millet (Setaria italica), which is a natural abiotic stress tolerant crop. To find out the dehydration-responsive miRNAs at the global level, four small RNA libraries were constructed from control and dehydration stress treated seedlings of two foxtail millet cultivars showing contrasting tolerance behavior towards dehydration stress. Using Illumina sequencing technology, 55 known and 136 novel miRNAs were identified, representing 22 and 48 miRNA families, respectively. Eighteen known and 33 novel miRNAs were differentially expressed during dehydration stress. After the stress treatment, 32 dehydration-responsive miRNAs were up-regulated in tolerant cultivar and 22 miRNAs were down-regulated in sensitive cultivar, suggesting that miRNA-mediated molecular regulation might play important roles in providing contrasting characteristics to these cultivars. Predicted targets of identified miRNAs were found to encode various transcription factors and functional enzymes, indicating their involvement in broad spectrum regulatory functions and biological processes. Further, differential expression patterns of seven known miRNAs were validated by northern blot and expression of ten novel dehydration-responsive miRNAs were confirmed by SL-qRT PCR. Differential expression behavior of five miRNA-target genes was verified under dehydration stress treatment and two of them also validated by RLM RACE. Overall, the present study highlights the importance of dehydration stress-associated post-transcriptional regulation governed by miRNAs and their targets in a naturally stress-tolerant model crop.

Global analysis of uncapped mRNA changes under drought stress and microRNA-dependent endonucleolytic cleavages in foxtail millet

BACKGROUND

mRNA degradation plays an important role in the determination of mRNA abundance and can quickly regulate gene expression. The production of uncapped mRNAs, an important mechanism of mRNA degradation, can be initiated by decapping enzymes, endonucleases or small RNAs such as microRNAs (miRNAs). Little is known, however, about the role of uncapped mRNAs in plants under environmental stress.

RESULTS

Using a novel approach called parallel analysis of RNA ends (PARE), we performed a global study of uncapped mRNAs under drought stress in foxtail millet (Setaria italica [L.] P. Beauv.). When both gene degradation (PARE) and gene transcription (RNA-sequencing) data were considered, four types of mRNA decay patterns were identified under drought stress. In addition, 385 miRNA-target interactions were identified in the PARE data using PAREsnip. The PARE analysis also suggested that two miRNA hairpin processing mechanisms--loop-last and loop-first processing--operate in foxtail millet, with both miR319 and miR156 gene families undergoing precise processing via the unusual loop-first mechanism. Finally, we found 11 C4 photosynthesis-related enzymes encoded by drought-responsive genes.

CONCLUSIONS

We performed a global analysis of mRNA degradation under drought stress and uncovered diverse drought-response mechanisms in foxtail millet. This information will deepen our understanding of mRNA expression under stressful environmental conditions in gramineous plants. In addition, PARE analysis identified many miRNA targets and revealed miRNA-precursor processing modes in foxtail millet.

Brachypodium distachyon and Setaria viridis: Model Genetic Systems for the Grasses

The family of grasses encompasses the world's most important food, feed, and bioenergy crops, yet we are only now beginning to develop the genetic resources to explore the diversity of form and function that underlies economically important traits. Two emerging model systems, Brachypodium distachyon and Setaria viridis, promise to greatly accelerate the process of gene discovery in the grasses and to serve as bridges in the exploration of panicoid and pooid grasses, arguably two of the most important clades of plants from a food security perspective. We provide both a historical view of the development of plant model systems and highlight several recent reports that are providing these developing communities with the tools for gene discovery and pathway engineering.