Evolution of Arabidopsis microRNA families through duplication events

Analysis of bZIP transcription factors in Rhododendron simsii and functional study of RsbZIP6 in regulating anthocyanin biosynthesis

The basic leucine zipper (bZIP) transcription factors play a critical role in various plant biological processes, including anthocyanin biosynthesis. This study focuses on Rhododendron simsii, a notable ornamental species with insufficiently explored bZIP transcription factors. We identified 66 bZIP transcription factors in the R. simsii genome and conducted comprehensive bioinformatics analyses to determine their gene localization, phylogenetic relationships, grouping, gene/protein structure, duplication events, synteny, and expression profiles. Our analysis identified RsbZIP6, a homolog of HY5 known to influence anthocyanin biosynthesis in many plants, as a potential regulator of this pathway. We cloned the complete coding sequence of RsbZIP6, which encodes a 170-amino acid protein spanning 510 bp. Subcellular localization analysis verified the nuclear presence of the RsbZIP6 protein. RT-qPCR analysis revealed the highest expression of RsbZIP6 in petals, which correlated with anthocyanin accumulation. Transgenic experiments indicated that overexpressing RsbZIP6 in Arabidopsis enhanced anthocyanin accumulation by upregulating genes involved in anthocyanin biosynthesis (4CL, CHS, CHI, DFR, F3H, F3'H, ANS and UF3GT). Our findings enhance understanding of the bZIP transcription factor family in R. simsii and underscore the vital role of RsbZIP6 in anthocyanin biosynthesis, providing insights for future genetic enhancement strategies.

Genome-wide survey of KT/HAK/KUP genes in the genus Citrullus and analysis of their involvement in K+-deficiency and drought stress responses in between C. lanatus and C. amarus

BACKGROUND

The KT/HAK/KUP is the largest K+ transporter family in plants, playing crucial roles in K+ absorption, transport, and defense against environmental stress. Sweet watermelon is an economically significant horticultural crop belonging to the genus Citrullus, with a high demand for K+ during its growth process. However, a comprehensive analysis of the KT/HAK/KUP gene family in watermelon has not been reported.

RESULTS

14 KT/HAK/KUP genes were identified in the genomes of each of seven Citrullus species. These KT/HAK/KUPs in watermelon were unevenly distributed across seven chromosomes. Segmental duplication is the primary driving force behind the expansion of the KT/HAK/KUP family, subjected to purifying selection during domestication (Ka/Ks < 1), and all KT/HAK/KUPs exhibit conserved motifs and could be phylogenetically classified into four groups. The promoters of KT/HAK/KUPs contain numerous cis-regulatory elements related to plant growth and development, phytohormone response, and stress response. Under K+ deficiency, the growth of watermelon seedlings was significantly inhibited, with cultivated watermelon experiencing greater impacts (canopy width, redox enzyme activity) compared to the wild type. All KT/HAK/KUPs in C. lanatus and C. amarus exhibit specific expression responses to K+-deficiency and drought stress by qRT-PCR. Notably, ClG42_07g0120700/CaPI482276_07g014010 were predominantly expressed in roots and were further induced by K+-deficiency and drought stress. Additionally, the K+ transport capacity of ClG42_07g0120700 under low K+ stress was confirmed by yeast functional complementation assay.

CONCLUSIONS

KT/HAK/KUP genes in watermelon were systematically identified and analyzed at the pangenome level and provide a foundation for understanding the classification and functions of the KT/HAK/KUPs in watermelon plants.

Genome-wide analysis of the bZIP gene family in Cinnamomum camphora ('Gantong 1') reveals the putative function in anthocyanin biosynthesis

Basic leucine zipper (bZIP) transcription factors (TFs) regulate plant development, growth, and secondary metabolism. The formation of red bark of new ornamental cultivar 'Gantong 1' is regulated mainly by anthocyanin anabolism. However, it is unclear whether and which bZIP TFs are involved in this process. We identified 89 genes encoding CcbZIP TFs in Cinnamomum camphora genome that could be divided into 14 subfamilies with similar gene structures and conserved motifs. CcbZIP38 and CcbZIP57 were highly conserved compared to HY5 in Arabidopsis thaliana and they were highly expressed in the bark and leaves of 'Gantong 1' at different stages. The target gene enrichment analysis showed that indicating indirect involvement of CcbZIP38 and CcbZIP57 in the regulation of anthocyanin synthesis. Our study contributes to understanding the molecular mechanism of anthocyanin synthesis regulation by CcbZIP TFs and provides a theoretical basis for genetic improvement of ornamental traits in C. camphora.

MicroRNA166: Old Players and New Insights into Crop Agronomic Traits Improvement

MicroRNA (miRNA), a type of non-coding RNA, is crucial for controlling gene expression. Among the various miRNA families, miR166 stands out as a highly conserved group found in both model and crop plants. It plays a key role in regulating a wide range of developmental and environmental responses. In this review, we explore the diverse sequences of MIR166s in major crops and discuss the important regulatory functions of miR166 in plant growth and stress responses. Additionally, we summarize how miR166 interacts with other miRNAs and highlight the potential for enhancing agronomic traits by manipulating the expression of miR166 and its targeted HD-ZIP III genes.

Identification and expression of the Di19 gene family in response to abiotic stress in common bean (Phaseolus vulgaris L.)

Drought-induced 19 (Di19) protein plays critical biological functions in response to adversity as well as in plant growth and development. Exploring the role and mechanism of Di19 in abiotic stress responses is of great significance for improving plant tolerance. In this study, six Di19 genes were identified in the common bean (Phaseolus vulgaris L.), which were mainly derived from segmental duplications. These genes share conserved exon/intron structures and were classified into three subfamilies based on their phylogenetic relationships. The composition and arrangement of conserved motifs were consistent with their phylogenetic relationships. Many hormone- and stress-responsive elements were distributed in the promoters region of PvDi19 genes. Variations in histidine residues in the Cys2/His2 (C2H2) zinc-finger domains resulted in an atypical tertiary structure of PvDi19-5. Gene expression analysis showed rapid induction of PvDi19-1 in roots by 10% PEG treatment, and PvDi19-2 in leaves by 20% PEG treatment, respectively. Most PvDi19s exhibited insensitivity to saline-alkali stress, except for PvDi19-6, which was notably induced during later stages of treatment. The most common bean Di19 genes were inhibited or not regulated by cadmium stress, but the expression of PvDi19-6 in roots was significantly upregulated when subjected to lower concentrations of cadmium (5 mmol). Moreover, Di19s exhibited greater sensitivity to severe cold stress (6°C). These findings enhance our understanding of the role of PvDi19s in common bean abiotic stress responses and provide a basis for future genetic enhancements in common bean stress tolerance.

Comprehensive in silico characterization of NAC transcription factor family of Pinellia ternata and functional analysis of PtNAC66 under high-temperature tolerance in transgenic Arabidopsis thaliana

Pinellia ternata, a valuable Chinese herb, suffers yield reduction due to "sprout tumble" under high temperatures. However, the mechanisms underlying its high-temperature stress remain poorly understood. NAM, ATAF1/2, and CUC2 (NAC) transcription factors regulate plant tissue growth and abiotic stress. Hence, there has been no comprehensive research conducted on NAC transcription factors in P. ternata. We identified 98 PtNAC genes unevenly distributed across 13 chromosomes, grouped into 15 families via phylogenetic analysis. Gene expression analysis revealed diverse expression patterns of PtNAC genes in different tissue types. Further studies revealed that PtNAC5/7/17/35/43/47/57/66/86 genes were highly expressed in various tissues of P. ternata and induced by heat stress, among which PtNAC66 was up-regulated at the highest folds induced by heat temperature. PtNAC66 is a nuclear protein that can selectively bind to the cis-responsive region NACRS but lacks the ability to activate transcription in yeast. For further research, PtNAC66 was cloned and transgenic Arabidopsis was obtained. PtNAC66 overexpression increased high-temperature tolerance compared to wild-type plants. Transcriptome profiling demonstrated that overexpression of PtNAC66 led to significant modification of genes responsible for regulating binding, catalytic activity, transcription regulator activity and transporter activity response genes. Additionally, PtNAC66 was found to bind the promoters of CYP707A3, MYB102 and NAC055, respectively, and inhibited their expression, affecting the high-temperature stress response in Arabidopsis. Our research established the foundation for functional studies of PtNAC genes in response to high-temperature forcing by characterizing the P. ternata NAC gene family and examining the biological role of PtNAC66 in plant high-temperature tolerance.

PbRbohH/J mediates ROS generation to regulate the growth of pollen tube in pear

Respiratory burst oxidase homolog (Rboh) family genes play crucial functions in development and growth. However, comprehensive and systematic investigation of Rboh family members in Rosaceae and their specific functions during pear pollen development are still limited. In the study, 63 Rboh genes were identified from eight Rosaceae genomes (Malus domestica, Pyrus bretschneideri, Pyrus communis, Prunus persica, Rubus occidentalis, Fragaria vesca, Prunus mume and Prunus avium) and divided into seven main subfamilies (I-VII) according to phylogenetic and structural features. Different modes of gene duplication led to the expansion of Rboh family, with purifying selection playing a vital role in the evolution of Rboh genes. In addition, RNA sequencing and qRT-PCR results indicated that PbRbohH and PbRbohJ were specifically high-expressed in pear pollen. Subsequently, subcellular localization revealed that PbRbohH/J distributed at the plasma membrane. Furthermore, by pharmacological analysis and antisense oligodeoxynucleotide assay, PbRbohH/J were demonstrated to mediate the formation of reactive oxygen species (ROS) to manage pollen tube growth. In conclusion, our results provide useful insights into the functions, expression patterns, evolutionary history of the Rboh genes in pear and other Rosaceae species.

Genome-wide identification and comprehensive analysis of the phytochrome-interacting factor (PIF) gene family in wheat

Phytochrome-interacting factors (PIFs) are essential transcription factors for plant growth, development, and stress responses. Although PIF genes have been extensively studied in many plant species, they have not been thoroughly investigated in wheat. Here, we identified 18 PIF genes in cultivated hexaploid wheat (Triticum aestivum L). Phylogenetic analysis, exon-intron structures, and motif compositions revealed the presence of four distinct groups of TaPIFs. Genome-wide collinearity analysis of PIF genes revealed the evolutionary history of PIFs in wheat, Oryza sativa, and Brachypodium distachyon. Cis-regulatory element analysis suggested that TaPIF genes indicated participated in plant development and stress responses. Subcellular localization assays indicated that TaPIF2-1B and TaPIF4-5B were transcriptionally active. Both were found to be localized to the nucleus. Gene expression analyses demonstrated that TaPIFs were primarily expressed in the leaves and were induced by various biotic and abiotic stresses and phytohormone treatments. This study provides new insights into PIF-mediated stress responses and lays a strong foundation for future investigation of PIF genes in wheat.

Identification and analysis of BAHD superfamily related to malonyl ginsenoside biosynthesis in Panax ginseng

Introduction

The BAHD (benzylalcohol O-acetyl transferase, anthocyanin O-hydroxycinnamoyl transferase, N-hydroxycinnamoyl anthranilate benzoyl transferase and deacetylvindoline 4-O-acetyltransferase), has various biological functions in plants, including catalyzing the biosynthesis of terpenes, phenolics and esters, participating in plant stress response, affecting cell stability, and regulating fruit quality.

Methods

Bioinformatics methods, real-time fluorescence quantitative PCR technology, and ultra-high-performance liquid chromatography combined with an Orbitrap mass spectrometer were used to explore the relationship between the BAHD gene family and malonyl ginsenosides in Panax ginseng.

Results

In this study, 103 BAHD genes were identified in P. ginseng, mainly distributed in three major clades. Most PgBAHDs contain cis-acting elements associated with abiotic stress response and plant hormone response. Among the 103 genes, 68 PgBAHDs are WGD (whole-genome duplication) genes. The significance of malonylation in biosynthesis has garnered considerable attention in the study of malonyltransferases. The phylogenetic tree results showed 34 PgBAHDs were clustered with genes that have malonyl characterization. Among them, seven PgBAHDs (PgBAHD4, 45, 65, 74, 90, 97, and 99) showed correlations > 0.9 with crucial enzyme genes involved in ginsenoside biosynthesis and > 0.8 with malonyl ginsenosides. These seven genes were considered potential candidates involved in the biosynthesis of malonyl ginsenosides.

Discussion

These results help elucidate the structure, evolution, and functions of the P. ginseng BAHD gene family, and establish the foundation for further research on the mechanism of BAHD genes in ginsenoside biosynthesis.

Systematic analysis of the Rboh gene family in seven gramineous plants and its roles in response to arbuscular mycorrhizal fungi in maize

BACKGROUND

Plant respiratory burst oxidase homolog (Rboh) gene family produces reactive oxygen species (ROS), and it plays key roles in plant-microbe interaction. Most Rboh gene family-related studies mainly focused on dicotyledonous plants; however, little is known about the roles of Rboh genes in gramineae.

RESULTS

A total of 106 Rboh genes were identified in seven gramineae species, including Zea mays, Sorghum bicolor, Brachypodium distachyon, Oryza sativa, Setaria italica, Hordeum vulgare, and Triticum aestivum. The Rboh protein sequences showed high similarities, suggesting that they may have conserved functions across different species. Duplication mode analysis detected whole-genome/segmental duplication (WGD)/(SD) and dispersed in the seven species. Interestingly, two local duplication (LD, including tandem and proximal duplication) modes were found in Z. mays, S. italica and H. vulgare, while four LD were detected in T. aestivum, indicating that these genes may have similar functions. Collinearity analysis indicated that Rboh genes are at a stable evolution state in all the seven species. Besides, Rboh genes from Z. mays were closely related to those from S. bicolor, consistent with the current understanding of plant evolutionary history. Phylogenetic analysis showed that the genes in the subgroups I and II may participate in plant-AM fungus symbiosis. Cis-element analysis showed that different numbers of elements are related to fungal induction in the promoter region. Expression profiles of Rboh genes in Z. mays suggested that Rboh genes had distinct spatial expression patterns. By inoculation with AM fungi, our transcriptome analysis showed that the expression of Rboh genes varies upon AM fungal inoculation. In particularly, ZmRbohF was significantly upregulated after inoculation with AM fungi. pZmRbohF::GUS expression analyses indicated that ZmRbohF was induced by arbuscular mycorrhizal fungi in maize. By comparing WT and ZmRbohF mutant, we found ZmRbohF had limited impact on the establishment of maize-AM fungi symbiosis, but play critical roles in regulating the proper development of arbuscules.

CONCLUSIONS

This study provides a comprehensive analysis of the evolution relationship of Rboh genes in seven gramineae species. Results showed that several Rboh genes regulate maize-AM fungal symbiosis process. This study provides valuable information for further studies of Rboh genes in gramineae.

Identification of SRS transcription factor family in Solanum lycopersicum, and functional characterization of their responses to hormones and abiotic stresses

The SHI RELATED SEQUENCE (SRS) family plays a vital role in the development of multiple plant organs such as floral meristem determinacy, organ morphogenesis, and signal transduction. Nevertheless, there is little understanding of the biological significance of tomato SRS family at this point. Our research identified eight SlSRS family members and classified them into three subfamilies based on phylogenetics, conserved motifs, and characteristic domain analysis. The intraspecies and interspecies collinearity analysis revealed clues of SRS family evolution. Many cis-elements related to hormones, stresses, and plant development can be found in the promoter region of SlSRS genes. All of eight SlSRS proteins were located in the nucleus and possessed transcriptional activity, half of which were transcriptional activators, and the other half were transcriptional repressors. Except for SlSRS1, which showed high transcript accumulation in vegetative organs, most SlSRS genes expressed ubiquitously in all flower organs. In addition, all SlSRS genes could significantly respond to at least four different plant hormones. Further, expression of SlSRS genes were regulated by various abiotic stress conditions. In summary, we systematically analyzed and characterized the SlSRS family, reviewed the expression patterns and preliminarily investigated the protein function, and provided essential information for further functional research of the tomato SRS genes in the determination of reproductive floral organs and the development of plants, and possibly other plants.

Genome-wide identification and bioinformatics analysis of the WD40 transcription factor family and candidate gene screening for anthocyanin biosynthesis in Rhododendron simsii

WD40 transcription factors (TFs) constitute a large gene family in eukaryotes, playing diverse roles in cellular processes. However, their functions in the major ornamental plant, Rhododendron simsii, remain poorly understood. In this study, we identified 258 WD40 proteins in the R. simsii genome, which exhibited an uneven distribution across chromosomes. Based on domain compositions and phylogenetic analysis, we classified these 258 RsWD40 proteins into 42 subfamilies and 47 clusters. Comparative genomic analysis suggested that the expansion of the WD40 gene family predates the divergence of green algae and higher plants, indicating an ancient origin. Furthermore, by analyzing the duplication patterns of RsWD40 genes, we found that transposed duplication played a major role in their expansion. Notably, the majority of RsWD40 gene duplication pairs underwent purifying selection during evolution. Synteny analysis identified significant orthologous gene pairs between R. simsii and Arabidopsis thaliana, Oryza sativa, Vitis vinifera, and Malus domestica. We also investigated potential candidate genes involved in anthocyanin biosynthesis during different flower development stages in R. simsii using RNA-seq data. Specifically, we identified 10 candidate genes during the bud stage and 7 candidate genes during the full bloom stage. GO enrichment analysis of these candidate genes revealed the potential involvement of the ubiquitination process in anthocyanin biosynthesis. Overall, our findings provide a valuable foundation for further investigation and functional analysis of WD40 genes, as well as research on the molecular mechanisms underlying anthocyanin biosynthesis in Rhododendron species.

Genome-wide identification and bioinformatics analysis of the WRKY transcription factors and screening of candidate genes for anthocyanin biosynthesis in azalea (Rhododendron simsii)

WRKY transcription factors have been demonstrated to influence the anthocyanin biosynthesis in many plant species. However, there is limited knowledge about the structure and function of WRKY genes in the major ornamental plant azalea (Rhododendron simsii). In this study, we identified 57 RsWRKY genes in the R. simsii genome and classified them into three main groups and several subgroups based on their structural and phylogenetic characteristics. Comparative genomic analysis suggested WRKY gene family has significantly expanded during plant evolution from lower to higher species. Gene duplication analysis indicated that the expansion of the RsWRKY gene family was primarily due to whole-genome duplication (WGD). Additionally, selective pressure analysis (Ka/Ks) suggested that all RsWRKY duplication gene pairs underwent purifying selection. Synteny analysis indicated that 63 and 24 pairs of RsWRKY genes were orthologous to Arabidopsis thaliana and Oryza sativa, respectively. Furthermore, RNA-seq data was used to investigate the expression patterns of RsWRKYs, revealing that 17 and 9 candidate genes may be associated with anthocyanin synthesis at the bud and full bloom stages, respectively. These findings provide valuable insights into the molecular mechanisms underlying anthocyanin biosynthesis in Rhododendron species and lay the foundation for future functional studies of WRKY genes.

Genome-Wide Characterization and Expression of the Hsf Gene Family in Salvia miltiorrhiza (Danshen) and the Potential Thermotolerance of SmHsf1 and SmHsf7 in Yeast

Salvia miltiorrhiza Bunge (Danshen) is a traditional Chinese herb with significant medicinal value. The yield and quality of Danshen are greatly affected by climatic conditions, in particular high temperatures. Heat shock factors (Hsfs) play important regulatory roles in plant response to heat and other environmental stresses. However, little is currently known about the role played by the Hsf gene family in S. miltiorrhiza. Here, we identified 35 SmHsf genes and classified them into three major groups: SmHsfA (n = 22), SmHsfB (n = 11), and SmHsfC (n = 2) using phylogenetic analysis. The gene structure and protein motifs were relatively conserved within subgroups but diverged among the different groups. The expansion of the SmHsf gene family was mainly driven by whole-genome/segmental and dispersed gene duplications. The expression profile of SmHsfs in four distinct organs revealed its members (23/35) are predominantly expressed in the root. The expression of a large number of SmHsfs was regulated by drought, ultraviolet, heat and exogenous hormones. Notably, the SmHsf1 and SmHsf7 genes in SmHsfB2 were the most responsive to heat and are conserved between dicots and monocots. Finally, heterologous expression analysis showed that SmHsf1 and SmHsf7 enhance thermotolerance in yeast. Our results provide a solid foundation for further functional investigation of SmHsfs in Danshen plants as a response to abiotic stresses.

Genome-Wide Characterization and Analysis of bHLH Transcription Factors Related to Anthocyanin Biosynthesis in Cinnamomum camphora ('Gantong 1')

Cinnamomum camphora is one of the most commonly used tree species in landscaping. Improving its ornamental traits, particularly bark and leaf colors, is one of the key breeding goals. The basic helix-loop-helix (bHLH) transcription factors (TFs) are crucial in controlling anthocyanin biosynthesis in many plants. However, their role in C. camphora remains largely unknown. In this study, we identified 150 bHLH TFs (CcbHLHs) using natural mutant C. camphora 'Gantong 1', which has unusual bark and leaf colors. Phylogenetic analysis revealed that 150 CcbHLHs were divided into 26 subfamilies which shared similar gene structures and conserved motifs. According to the protein homology analysis, we identified four candidate CcbHLHs that were highly conserved compared to the TT8 protein in A. thaliana. These TFs are potentially involved in anthocyanin biosynthesis in C. camphora. RNA-seq analysis revealed specific expression patterns of CcbHLHs in different tissue types. Furthermore, we verified expression patterns of seven CcbHLHs (CcbHLH001, CcbHLH015, CcbHLH017, CcbHLH022, CcbHLH101, CcbHLH118, and CcbHLH134) in various tissue types at different growth stages using qRT-PCR. This study opens a new avenue for subsequent research on anthocyanin biosynthesis regulated by CcbHLH TFs in C. camphora.

Genome-Wide Identification and Phylogenetic Analysis of TRP Gene Family Members in Saurian

The transient receptor potential plays a critical role in the sensory nervous systems of vertebrates in response to various mechanisms and stimuli, such as environmental temperature. We studied the physiological adaptive evolution of the TRP gene in the saurian family and performed a comprehensive analysis to identify the evolution of the thermo-TRPs channels. All 251 putative TRPs were divided into 6 subfamilies, except TRPN, from the 8 saurian genomes. Multiple characteristics of these genes were analyzed. The results showed that the most conserved proteins of TRP box 1 were located in motif 1, and those of TRP box 2 were located in motif 10. The TRPA and TRPV in saurian tend to be one cluster, as a sister cluster with TRPC, and the TRPM is the root of group I. The TRPM, TRPV, and TRPP were clustered into two clades, and TRPP were organized into TRP PKD1-like and PKD2-like. Segmental duplications mainly occurred in the TRPM subfamily, and tandem duplications only occurred in the TRPV subfamily. There were 15 sites to be under positive selection for TRPA1 and TRPV2 genes. In summary, gene structure, chromosomal location, gene duplication, synteny analysis, and selective pressure at the molecular level provided some new evidence for genetic adaptation to the environment. This result provides a basis for identifying and classifying TRP genes and contributes to further elucidating their potential function in thermal sensors.

Evolution and function of ubiquitin-specific proteases (UBPs): Insight into seed development roles in tung tree (Vernicia fordii)

The tung oil produced by the tung tree (Vernicia fordii) provides resources for the manufacture of biodiesel. Ubiquitin-specific proteases (UBPs) are the largest group of deubiquitinases and play key roles in regulating development and stress responses. Here, 21 UBPs were identified in V. fordii, roughly one-half the number found in Manihot esculenta and Hevea brasiliensis. Most UBP duplications are produced from whole-genome duplication (WGD), and significant differences in gene retention existed among Euphorbiaceae. The great majority of UBP-containing blocks in V. fordii, V. montana, Ricinus communis, and Jatropha curcas exhibited extensive conservation with the duplicated regions of M. esculenta and H. brasiliensis. These blocks formed 14 orthologous groups, indicating they shared WGD with UBPs in M. esculenta and H. brasiliensis, but most of these UBPs copies were lost. The UBP orthologs contained significant functional divergence which explained the susceptibility of V. fordii to Fusarium wilt and the resistance of V. montana to Fusarium wilt. The expression patterns and experiments suggested that Vf03G1417 could affect the seed-related traits and positively regulate the seed oil accumulation. This study provided important insights into the evolution of UBPs in Euphorbiaceae and identified important candidate VfUBPs for marker-assisted breeding in V. fordii.

Chromosome-level genome assembly of Mentha longifolia L. reveals gene organization underlying disease resistance and essential oil traits

Mentha longifolia (L.) Huds., a wild, diploid mint species, has been developed as a model for mint genetic and genomic research to aid breeding efforts that target Verticillium wilt disease resistance and essential oil monoterpene composition. Here, we present a near-complete, chromosome-scale mint genome assembly for M. longifolia USDA accession CMEN 585. This new assembly is an update of a previously published genome draft, with dramatic improvements. A total of 42,107 protein-coding genes were annotated and placed on 12 chromosomal scaffolds. One hundred fifty-three genes contained conserved sequence domains consistent with nucleotide binding site-leucine-rich-repeat plant disease resistance genes. Homologs of genes implicated in Verticillium wilt resistance in other plant species were also identified. Multiple paralogs of genes putatively involved in p-menthane monoterpenoid biosynthesis were identified and several cases of gene clustering documented. Heterologous expression of candidate genes, purification of recombinant target proteins, and subsequent enzyme assays allowed us to identify the genes underlying the pathway that leads to the most abundant monoterpenoid volatiles. The bioinformatic and functional analyses presented here are laying the groundwork for using marker-assisted selection in improving disease resistance and essential oil traits in mints.

Deciphering the role of miRNA in reprogramming plant responses to drought stress

Drought is the most prevalent environmental stress that affects plants' growth, development, and crop productivity. However, plants have evolved adaptive mechanisms to respond to the harmful effects of drought. They reprogram their: transcriptome, proteome, and metabolome that alter their cellular and physiological processes and establish cellular homeostasis. One of the crucial regulatory processes that govern this reprogramming is post-transcriptional regulation by microRNAs (miRNAs). miRNAs are small non-coding RNAs, involved in the downregulation of the target mRNA via translation inhibition/mRNA degradation/miRNA-mediated mRNA decay/ribosome drop off/DNA methylation. Many drought-inducible miRNAs have been identified and characterized in plants. Their main targets are regulatory genes that influence growth, development, osmotic stress tolerance, antioxidant defense, phytohormone-mediated signaling, and delayed senescence during drought stress. Overexpression of drought-responsive miRNAs (Osa-miR535, miR160, miR408, Osa-miR393, Osa-miR319, and Gma-miR394) in certain plants has led to tolerance against drought stress indicating their vital role in stress mitigation. Similarly, knock down (miR166/miR398c) or deletion (miR169 and miR827) of miRNAs has also resulted in tolerance to drought stress. Likewise, engineered Arabidopsis plants with miR165, miR166 using short tandem target mimic strategy, exhibited drought tolerance. Since miRNAs regulate the expression of an array of drought-responsive genes, they can act as prospective targets for genetic manipulations to enhance drought tolerance in crops and achieve sustainable agriculture. Further investigations toward functional characterization of diverse miRNAs, and understanding stress-responses regulated by these miRNAs and their utilization in biotechnological applications is highly recommended.

Genome-wide identification and comparative analysis of the PYL gene family in eight Rosaceae species and expression analysis of seeds germination in pear

Abscisic acid (ABA) is a phytohormone that plays important roles in the regulation of plant growth, seed germination, and stress responses. The pyrabactin resistance 1-like (PYR/PYL) protein, an ABA receptor, was involved in the initial step in ABA signal transduction. However, the evolutionary history and characteristics of PYL genes expression remain unclear in Chinese white pear (Pyrus bretschneideri) or other Rosaceae species. In this study, 67 PYL genes were identified in eight Rosaceae species, and have been classified into three subgroups based on specific motifs and phylogenetic analysis. Intriguingly, we observed that whole-genome duplication (WGD) and dispersed duplication (DSD) have a major contribution to PYL family expansion. Purifying selection was the major force in PYL genes evolution. Expression analysis finds that PYL genes may function in multiple pear tissues. qRT-PCR validation of 11 PbrPYL genes indicates their roles in seed germination and abiotic stress responses. Our study provides a basis for further elucidation of the function of PYL genes and analysis of their expansion, evolution and expression patterns, which helps to understand the molecular mechanism of pear response to seed germination and seedling abiotic stress.

Genome-wide identification and expression analysis of HAK/KUP/KT potassium transporter provides insights into genes involved in responding to potassium deficiency and salt stress in pepper (Capsicum annuum L.)

In plants, the HAK/KUP/KT family is the largest group of potassium transporters, and it plays an important role in mineral element absorption, plant growth, environmental stress adaptation, and symbiosis. Although these important genes have been investigated in many plant species, limited information is currently available on the HAK/KUP/KT genes for Pepper (Capsicum annuum L.). In the present study, a total of 20 CaHAK genes were identified from the pepper genome and the CaHAK genes were numbered 1 - 20 based on phylogenetic analysis. For the genes and their corresponding proteins, the physicochemical properties, phylogenetic relationship, chromosomal distribution, gene structure, conserved motifs, gene duplication events, and expression patterns were analyzed. Phylogenetic analysis divided CaHAK genes into four cluster (I-IV) based on their structural features and the topology of the phylogenetic tree. Purifying selection played a crucial role in the evolution of CaHAK genes, while whole-genome triplication contributed to the expansion of the CaHAK gene family. The expression patterns showed that CaHAK proteins exhibited functional divergence in terms of plant K⁺ uptake and salt stress response. In particular, transcript abundance of CaHAK3 and CaHAK7 was strongly and specifically up-regulated in pepper roots under low K⁺ or high salinity conditions, suggesting that these genes are candidates for high-affinity K⁺ uptake transporters and may play crucial roles in the maintenance of the Na⁺/K⁺ balance during salt stress in pepper. In summary, the results not only provided the important information on the characteristics and evolutionary relationships of CaHAKs, but also provided potential genes responding to potassium deficiency and salt stress.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-022-03136-z.

Impact of polyploidy on plant tolerance to abiotic and biotic stresses

Polyploidy, defined as the coexistence of three or more complete sets of chromosomes in an organism's cells, is considered as a pivotal moving force in the evolutionary history of vascular plants and has played a major role in the domestication of several crops. In the last decades, improved cultivars of economically important species have been developed artificially by inducing autopolyploidy with chemical agents. Studies on diverse species have shown that the anatomical and physiological changes generated by either natural or artificial polyploidization can increase tolerance to abiotic and biotic stresses as well as disease resistance, which may positively impact on plant growth and net production. The aim of this work is to review the current literature regarding the link between plant ploidy level and tolerance to abiotic and biotic stressors, with an emphasis on the physiological and molecular mechanisms responsible for these effects, as well as their impact on the growth and development of both natural and artificially generated polyploids, during exposure to adverse environmental conditions. We focused on the analysis of those types of stressors in which more progress has been made in the knowledge of the putative morpho-physiological and/or molecular mechanisms involved, revealing both the factors in common, as well as those that need to be addressed in future research.

Plant RNA-mediated gene regulatory network

Not all transcribed RNAs are protein-coding RNAs. Many of them are non-protein-coding RNAs in diverse eukaryotes. However, some of them seem to be non-functional and are resulted from spurious transcription. A lot of non-protein-coding transcripts have a significant function in the translation process. Gene expressions depend on complex networks of diverse gene regulatory pathways. Several non-protein-coding RNAs regulate gene expression in a sequence-specific system either at the transcriptional level or post-transcriptional level. They include a significant part of the gene expression regulatory network. RNA-mediated gene regulation machinery is evolutionarily ancient. They well-evolved during the evolutionary time and are becoming much more complex than had been expected. In this review, we are trying to summarizing the current knowledge in the field of RNA-mediated gene silencing.

Genome-wide identification and expression analysis of the bZIP transcription factors, and functional analysis in response to drought and cold stresses in pear (Pyrus breschneideri)

BACKGROUND

Transcription factors (TFs) are involved in many important biological processes, including cell stretching, histological differentiation, metabolic activity, seed storage, gene regulation, and response to abiotic and biotic stresses. Little is known about the functions, evolutionary history, and expression patterns of basic region-leucine zipper TF family genes in pear, despite the release of the genome of Chinese white pears ("Dangshansuli").

RESULTS

Overall, 92 bZIP genes were identified in the pear genome (Pyrus breschneideri). Of these, 83 were randomly distributed on all 17 chromosomes except chromosome 4, and the other 9 genes were located on loose scaffolding. The genes were divided into 14 subgroups. Whole-genome duplications, dispersed duplication, and purifying selection for whole-genome duplications are the main reasons for the expansion of the PbrbZIP gene family. The analysis of functional annotation enrichment indicated that most of the functions of PbrbZIP genes were enriched in Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathways involved in the abiotic stress response. Next, expression analysis and virus-induced gene silencing results indicated that PbrbZIP genes might play critical roles in response to drought and cold stresses, especially for the genes from subgroups A, C, G, I, and S.

CONCLUSIONS

Ninety-two PbrbZIP genes were identified from the pear genome and classified into 14 subgroups. PbrbZIP genes were mainly expanded from whole-genome duplications and dispersed duplications and retained by purifying selection. PbrbZIP genes were induced by cold and drought stresses and played important roles in drought and cold tolerance. These results provided useful information for further increasing the tolerance of pears to stresses and a foundation to study the cold and drought tolerance mechanism of PbrbZIP genes.

Deciphering Evolutionary Dynamics of WRKY I Genes in Rosaceae Species

WRKY transcription factors participate in various regulation processes at different developmental stages in higher plants. Here, 98 WRKY I genes were identified in seven Rosaceae species. The WRKY I genes are highly enriched in some subgroups and are selectively expanded in Chinese pear [Pyrus bretschneideri (P. bretschneideri)] and apple [Malus domestica (M. domestica)]. By searching for intra-species gene microsynteny, we found the majority of chromosomal segments for WRKY I-containing segments in both P. bretschneideri and M. domestica genomes, while paired segments were hardly identified in the other five genomes. Furthermore, we analyzed the environmental selection pressure of duplicated WRKY I gene pairs, which indicated that the strong purifying selection for WRKY domains may contribute to the stability of its structure and function. The expression patterns of duplication PbWRKY genes revealed that functional redundancy for some of these genes was derived from common ancestry and neo-functionalization or sub-functionalization for some of them. This study traces the evolution of WRKY I genes in Rosaceae genomes and lays the foundation for functional studies of these genes in the future. Our results also show that the rates of gene loss and gain in different Rosaceae genomes are far from equilibrium.

Exploring diverse roles of micro RNAs in banana: Current status and future prospective

Micro RNAs (miRNAs) are 20-24 nucleotides long non-coding RNA sequences identified and characterized in multiple plant and animal systems. miRNAs play multifarious roles ranging from plant development to stress tolerance by synchronizing physiological processes at the level of transcription and translation. Banana is a major horticultural crop with colossal production worldwide. Despite the recent encouraging developments, the information on functions of miRNAs in banana plants is still in its infancy. The available literature pertaining to miRNAs in banana plants hints towards their contribution as master regulators in crucial physiological processes for instance abiotic stress responses, pathogenic defence response, fruit ripening and so on. This review is focused on biogenesis of miRNAs, their identification and deciphering their respective roles in banana plants with special emphasis on abiotic stress responses, plant immune responses, fruit ripening and storage. Based on the prior reports, we identified a few miRNAs with prospective roles in stress tolerance and illustrated the potential applications of miRNAs in banana crop improvement utilizing recent biotechnological tools such as CRISPR cas9, RNAi and the nano particle based foliar spray of miRNAs. The review briefly explained the future directions in banana research with a special emphasis on miRNA regulatory networks and agronomic traits improvement. Finally, future domains in miRNA research in plants and their possible applications towards crop improvement in agriculture are described briefly.

The APETALA2 homolog CaFFN regulates flowering time in pepper

Flowering time is an important agronomic trait that contributes to fitness in plants. However, the genetic basis of flowering time has not been extensively studied in pepper. To understand the genetics underlying flowering time, we constructed an F₂ population by crossing a spontaneous early flowering mutant and a late-flowering pepper line. Using bulked segregant RNA-seq, a major locus controlling flowering time in this population was mapped to the end of chromosome 2. An APETALA2 (AP2) homolog (CaFFN) cosegregated with flowering time in 297 individuals of the F₂ population. A comparison between the parents revealed a naturally occurring rare SNP (SNP2T > C) that resulted in the loss of a start codon in CaFFN in the early flowering mutant. Transgenic Nicotiana benthamiana plants with high CaFFN expression exhibited a delay in flowering time and floral patterning defects. On the other hand, pepper plants with CaFFN silencing flowered early. Therefore, the CaFFN gene acts as a flowering repressor in pepper. CaFFN may function as a transcriptional activator to activate the expression of CaAGL15 and miR156e and as a transcriptional repressor to repress the expression of CaAG, CaAP1, CaSEP3, CaSOC1, and miR172b based on a qRT-PCR assay. Direct activation of CaAGL15 by CaFFN was detected using yeast one-hybrid and dual-luciferase reporter assays, consistent with the hypothesis that CaFFN regulates flowering time. Moreover, the CaFFN gene association analysis revealed a significant association with flowering time in a natural pepper population, indicating that the CaFFN gene has a broad effect on flowering time in pepper. Finally, the phylogeny, evolutionary expansion and expression patterns of CaFFN/AP2 homologs were analyzed to provide valuable insight into CaFFN. This study increases our understanding of the involvement of CaFFN in controlling flowering time in pepper, thus making CaFFN a target gene for breeding early maturing pepper.

Genome-wide identification of the BURP domain-containing genes in Phaseolus vulgaris

Plant-specific BURP domain-containing proteins have an essential role in the plant's development and stress responses. Although BURP domain-containing proteins have been identified in several plant species, genome-wide analysis of the BURP gene family has not been investigated in the common bean. In the present study, we identified 11 BURP family members in the common bean (Phaseolus vulgaris) genome with a comprehensive in silico analysis. Pairwise alignment and phylogenetic analyses grouped PvBURP members into four subfamilies [RD-22 like (3), PG1β-like (4), BNM2-like (3), and USP-like (1)] according to their amino acid motifs, protein domains and intron-exon structure. The physical and biochemical characteristics of amino acids, motif and intron-exon structure, and cis-regulatory elements of BURPs members were determined. Promoter regions of BURP members included stress, light, and hormone response-related cis-elements. Therefore, expression profiles of PvBURP genes were identified with in silico tools and qRT-PCR analyses under stress (salt and drought) and hormone treatment (ABA, IAA) in the current study. While significant activity changes were not observed in BURP genes in RNA-seq data sets related to salt stress, it was determined that some BURP genes were expressed differently in those with drought stress. We identified 12 different miRNA, including miRNA395, miRNA156, miRNA169, miRNA171, miRNA319, and miRNA390, targeting the nine PvBURP genes using two different in silico tools based on perfect or near-perfect complementarity to their targets. Here we present the first study to identify and characterize the BURP genes in common bean using whole-genome analysis, and the findings may serve as a reference for future functional research in common bean.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01052-9.

The BBX gene family in Moso bamboo (Phyllostachys edulis): identification, characterization and expression profiles

BACKGROUND

The BBX (B-box) family are zinc finger protein (ZFP) transcription factors that play an essential role in plant growth, development and response to abiotic stresses. Although BBX genes have been characterized in many model organisms, genome-wide identification of the BBX family genes have not yet been reported in Moso bamboo (Phyllostachys edulis), and the biological functions of this family remain unknown.

RESULT

In the present study, we identified 27 BBX genes in the genome of Moso bamboo, and analysis of their conserved motifs and multiple sequence alignments revealed that they all shared highly similar structures. Additionally, phylogenetic and homology analyses indicated that PeBBX genes were divided into three clusters, with whole-genome duplication (WGD) events having facilitated the expansion of this gene family. Light-responsive and stress-related cis-elements were identified by analyzing cis-elements in the promoters of all PeBBX genes. Short time-series expression miner (STEM) analysis revealed that the PeBBX genes had spatiotemporal-specific expression patterns and were likely involved in the growth and development of bamboo shoots. We further explored the downstream target genes of PeBBXs, and GO/KEGG enrichment analysis predicted multiple functions of BBX target genes, including those encoding enzymes involved in plant photosynthesis, pyruvate metabolism and glycolysis/gluconeogenesis.

CONCLUSIONS

In conclusion, we analyzed the PeBBX genes at multiple different levels, which will contribute to further studies of the BBX family and provide valuable information for the functional validation of this family.

Analysis of Evolution, Expression and Genetic Transformation of TCP Transcription Factors in Blueberry Reveal That VcTCP18 Negatively Regulates the Release of Flower Bud Dormancy

Plant-specific TEOSINTE BRANCHED 1, CYCLOIDEA, PROLIFERATING CELL FACTORS (TCP) transcription factors have versatile functions in plant growth, development and response to environmental stress. Despite blueberry's value as an important fruit crop, the TCP gene family has not been systematically studied in this plant. The current study identified blueberry TCP genes (VcTCPs) using genomic data from the tetraploid blueberry variety 'Draper'; a total of 62 genes were obtained. Using multiple sequence alignment, conserved motif, and gene structure analyses, family members were divided into two subfamilies, of which class II was further divided into two subclasses, CIN and TB1. Synteny analysis showed that genome-wide or segment-based replication played an important role in the expansion of the blueberry TCP gene family. The expression patterns of VcTCP genes during fruit development, flower bud dormancy release, hormone treatment, and tissue-specific expression were analyzed using RNA-seq and qRT-PCR. The results showed that the TB1 subclass members exhibited a certain level of expression in the shoot, leaf, and bud; these genes were not expressed during fruit development, but transcript levels decreased uniformly during the release of flower bud dormancy by low-temperature accumulation. The further transgenic experiments showed the overexpression of VcTCP18 in Arabidopsis significantly decreased the seed germination rate in contrast to the wild type. The bud dormancy phenomena as late-flowering, fewer rosettes and main branches were also observed in transgenic plants. Overall, this study provides the first insight into the evolution, expression, and function of VcTCP genes, including the discovery that VcTCP18 negatively regulated bud dormancy release in blueberry. The results will deepen our understanding of the function of TCPs in plant growth and development.

Conservation and divergence of the TaSOS1 gene family in salt stress response in wheat (Triticum aestivum L.)

Salinity is one of the most important problems that adversely affect crops growth, productivity and quality worldwide. Salt Overly Sensitive 1 (SOS1) gene family plays vital roles in plant response to salt stress. Herein, we report the identification of the SOS family in wheat and the exploration of the expression profiles of SOSs under salt stress. Complete genome sequences of T. aestivum were downloaded from Ensembl plant database. Conservation and divergence of TaSOS1 family were conducted by using phylogenetic tree, gene structure and synteny distribution analysis. Expression profiles of TaSOS1s were obtained based on transcriptome and qRT-PCR analysis. Totally, 119 TaSOS1 proteins in wheat were identified at the genome-wide level and classified into three groups. Six motifs were conserved in TaSOS1 gene family. Moreover, 25 TaSOS1 genes had three copies distributing in three sub-genomes (A, B and D). A total of 32, 28 and 29 TaSOS1 genes were located on the sub-genomes A, B and D, respectively. Moreover, there were 19, 12, 6, 7, 28, 5 and 12 genes located on the three homologous of chromosomes 1, 2, 3, 4, 5, 6 and 7, respectively. Two genes were mapped to unattributed scaffolds. The duplication events analysis indicated that tandem repeats contributed to the expansion of the SOS1 family in wheat. Collinearity analysis demonstrated that segmental duplications play an important role in the expansion of SOS1 members. Chromosome 7, 5, 3, and 2 showed collinear relationship. Tissue specific expression pattern analysis revealed that 41 TaSOS1 genes expressed in various tissues, such as root, shoot, leaf, spike and grain. Transcriptomic analysis revealed that 28 and 26 genes were up- and down-regulated under salinity stress, respectively, of which 18 genes were further confirmed by RT-qPCR. The plants with high expression level of these genes displayed higher tolerance to salinity stress, stronger root system, higher Fv/Fm value and water potential. The results could be helpful for further elucidating the molecular mechanism of TaSOS1 related to salt tolerance in wheat and provide a toolkit for improving the salinity tolerance of wheat.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01009-y.

Genome-wide identification, evolutionary estimation and functional characterization of two cotton CKI gene types

BACKGROUND

Casein kinase I (CKI) is a kind of serine/threonine protein kinase highly conserved in plants and animals. Although molecular function of individual member of CKI family has been investigated in Arabidopsis, little is known about their evolution and functions in Gossypium.

RESULTS

In this study, five cotton species were applied to study CKI gene family in cotton, twenty-two species were applied to trace the origin and divergence of CKI genes. Four important insights were gained: (i) the cotton CKI genes were classified into two types based on their structural characteristics; (ii) two types of CKI genes expanded with tetraploid event in cotton; (iii) two types of CKI genes likely diverged about 1.5 billion years ago when red and green algae diverged; (iv) two types of cotton CKI genes which highly expressed in leaves showed stronger response to photoperiod (circadian clock) and light signal, and most two types of CKI genes highly expressed in anther showed identical heat inducible expression during anther development in tetraploid cotton (Gossypium hirsutum).

CONCLUSION

This study provides genome-wide insights into the evolutionary history of cotton CKI genes and lays a foundation for further investigation of the functional differentiation of two types of CKI genes in specific developmental processes and environmental stress conditions.

Identification and Expression Analysis of SLAC/SLAH Gene Family in Brassica napus L

Slow type anion channels (SLAC/SLAHs) play important roles during anion transport, growth and development, abiotic stress responses and hormone responses in plants. However, there is few report on SLAC/SLAHs in rapeseed (Brassica napus). Genome-wide identification and expression analysis of SLAC/SLAH gene family members were performed in B. napus. A total of 23 SLAC/SLAH genes were identified in B. napus. Based on the structural characteristics and phylogenetic analysis of these members, the SLAC/SLAHs could be classified into three main groups. Transcriptome data demonstrated that BnSLAH3 genes were detected in various tissues of the rapeseed and could be up-regulated by low nitrate treatment in roots. BnSLAC/SLAHs were exclusively localized on the plasma membrane in transient expression of tobacco leaves. These results will increase our understanding of the evolution and expression of the SLAC/SLAHs and provide evidence for further research of biological functions of candidates in B. napus.

MicroRNA annotation in plants: current status and challenges

Last two decades, the studies on microRNAs (miRNAs) and the numbers of annotated miRNAs in plants and animals have surged. Herein, we reviewed the current progress and challenges of miRNA annotation in plants. Via the comparison of plant and animal miRNAs, we pinpointed out the difficulties on plant miRNA annotation and proposed potential solutions. In terms of recalling the history of methods and criteria in plant miRNA annotation, we detailed how the major progresses made and evolved. By collecting and categorizing bioinformatics tools for plant miRNA annotation, we surveyed their advantages and disadvantages, especially for ones with the principle of mimicking the miRNA biogenesis pathway by parsing deeply sequenced small RNA (sRNA) libraries. In addition, we summarized all available databases hosting plant miRNAs, and posted the potential optimization solutions such as how to increase the signal-to-noise ratio (SNR) in these databases. Finally, we discussed the challenges and perspectives of plant miRNA annotations, and indicated the possibilities offered by an all-in-one tool and platform according to the integration of artificial intelligence.

Genome-wide identification of PbrbHLH family genes, and expression analysis in response to drought and cold stresses in pear (Pyrus bretschneideri)

BACKGROUND

The basic helix-loop-helix (bHLH) transcription factors play important roles in many processes in plant growth, metabolism and responses to abiotic stresses. Although, the sequence of Chinese white pear genome (cv. 'Dangshansuli') has already been reported, there is still a lack of clarity regarding the bHLH family genes and their evolutionary history.

RESULTS

In this work, a genome-wide identification of the bHLH genes in Chinese white pear was performed, and we characterized the functional roles of these PbrbHLH genes in response to abiotic stresses. Based on the phylogenetic analysis and structural characteristics, 197 identified bHLH genes could be well classified into 21 groups. Expansion of PbrbHLH gene family was mainly driven by WGD and dispersed duplication with the purifying selection from the recent WGD. The functional annotation enrichment showed that the majority of PbrbHLHs were enriched in the GO terms and KEGG pathways involved in responds to stress conditions as TFs. Transcriptomic profiles and qRT-PCR revealed that PbrbHLH7, PbrbHLH8, PbrbHLH128, PbrbHLH160, PbrbHLH161 and PbrbHLH195 were significantly up-regulated under cold and drought treatments. In addition, PbrbHLH195-silenced pear seedlings display significant reduced cold tolerance, exhibiting reduced chlorophyll content, as well as increased electrolyte leakage and concentrations of malondialdehyde and H₂O₂.

CONCLUSION

For the first time, a comprehensive analysis identified the bHLH genes in Chinese white pear and demonstrated that PbrbHLH195 is involved in the production of ROS in response to cold stress, suggesting that members of the PbrbHLH family play an essential role in the stress tolerance of pear.

Redundant and specific roles of individual MIR172 genes in plant development

Evolutionarily conserved microRNAs (miRNAs) usually have high copy numbers in the genome. The redundant and specific roles of each member of a multimember miRNA gene family are poorly understood. Previous studies have shown that the miR156-SPL-miR172 axis constitutes a signaling cascade in regulating plant developmental transitions. Here, we report the feasibility and utility of CRISPR-Cas9 technology to investigate the functions of all 5 MIR172 family members in Arabidopsis. We show that an Arabidopsis plant devoid of miR172 is viable, although it displays pleiotropic morphological defects. MIR172 family members exhibit distinct expression pattern and exert functional specificity in regulating meristem size, trichome initiation, stem elongation, shoot branching, and floral competence. In particular, we find that the miR156-SPL-miR172 cascade is bifurcated into specific flowering responses by matching pairs of coexpressed SPL and MIR172 genes in different tissues. Our results thus highlight the spatiotemporal changes in gene expression that underlie evolutionary novelties of a miRNA gene family in nature. The expansion of MIR172 genes in the Arabidopsis genome provides molecular substrates for the integration of diverse floral inductive cues, which ensures that plants flower at the optimal time to maximize seed yields.

This study uses CRISPR-Cas9 technology to investigate the functions of all five miR172 genes in Arabidopsis, finding that miRNA172 family members exhibit distinct expression pattern and exert functional specificity in regulating meristem size, trichome initiation, stem elongation, shoot branching and floral competence.

Maize microRNA166 Inactivation Confers Plant Development and Abiotic Stress Resistance

MicroRNAs are important regulators in plant developmental processes and stress responses. In this study, we generated a series of maize STTM166 transgenic plants. Knock-down of miR166 resulted in various morphological changes, including rolled leaves, enhanced abiotic stress resistance, inferior yield-related traits, vascular pattern and epidermis structures, tassel architecture, as well as abscisic acid (ABA) level elevation and indole acetic acid (IAA) level reduction in maize. To profile miR166 regulated genes, we performed RNA-seq and qRT-PCR analysis. A total of 178 differentially expressed genes (DEGs) were identified, including 118 up-regulated and 60 down-regulated genes. These DEGs were strongly enriched in cell and intercellular components, cell membrane system components, oxidoreductase activity, single organism metabolic process, carbohydrate metabolic process, and oxidation reduction process. These results indicated that miR166 plays important roles in auxin and ABA interaction in monocots, yet the specific mechanism may differ from dicots. The enhanced abiotic stress resistance is partly caused via rolling leaves, high ABA content, modulated vascular structure, and the potential changes of cell membrane structure. The inferior yield-related traits and late flowering are partly controlled by the decreased IAA content, the interplay of miR166 with other miRNAs and AGOs. Taken together, the present study uncovered novel functions of miR166 in maize, and provide insights on applying short tandem target mimics (STTM) technology in plant breeding.

Identification of genuine and novel miRNAs in Amaranthus hypochondriacus from high-throughput sequencing data

Amaranth has been proposed as an exceptional alternative for food security and climate change mitigation. Information about the distribution, abundance, or specificity of miRNAs in amaranth species is scare. Here, small RNAs from seedlings under control, drought, heat, and cold stress conditions of the Amaranthus hypocondriacus variety "Gabriela" were sequenced and miRNA loci identified in the amaranth genome using the ShortStack software. Fifty-three genuine miRNA clustersthirty-nine belonging to conserved families, and fourteen novel, were identified. Identification of their target genes suggests that conserved amaranth miRNAs are involved in growth and developmental processes, as well as stress responses. MiR0005, an amaranth-specific miRNA, exhibited an unusual high level of expression, akin to that of conserved miRNAs. Overall, our results broaden our knowledge regarding the distribution, abundance and expression of miRNAs in amaranth, providing the basis for future research on miRNAs and their functions in this important species.

Independent recurrent evolution of MICRORNA genes converging onto similar non-canonical organisation across green plant lineages is driven by local and segmental duplication events in species, family and lineages

The relationship between evolutionary history, organisation and transcriptional regulation of genes are intrinsically linked. These have been well studied in canonically organised protein-coding genes but not of MIRNAs. In the present study, we investigated the non-canonical arrangement of MIRNAs across taxonomic boundaries from algae to angiosperms employing a combination of genome organization, phylogeny and synteny. We retrieved the complete dataset of MIRNA from twenty-five species to identify and classify based on organisational patterns. The median size of cluster was between 2-5 kb and between 1-20 % of all MIRNAs are organized in head-to-head (with bidirectional promoter), head-to-tail (tandem), and overlapping manner. Although majority of the clusters are composed of MIRNA homologs, 25% of all clusters comprises of non-homologous genes with a potential of generating functional and regulatory complexity. A comparison of phylogeny and organizational patterns revealed that multiple independent events, some of which are species-specific, and some ancient, in different lineages, are responsible for non-canonical organization. Detailed investigation of MIR395 family across the plants revealed a complex origin of non-canonical arrangement through ancient and recent, segmental and local duplications; analysis of MIR399 family revealed major expansion occurred prior to monocot-dicot split, with few lineage-specific events. Evolution of "convergent" organization pattern of non-canonical arrangement originating from independent loci through recurrent event highlights our poor understanding of evolutionary process of MIRNA genes. The present investigation thus paves way for comparative functional genomics to understand the role of non-canonical organization on transcriptional regulation and regulatory diversity in MIRNA gene families.

Genome-wide survey of sucrose non-fermenting 1-related protein kinase 2 in Rosaceae and expression analysis of PbrSnRK2 in response to ABA stress

Background

The members of the sucrose non-fermenting 1-related protein kinase 2 (SnRK2) family are specific serine/threonine protein kinases in plants that play important roles in stress signal transduction and adaptation. Because of their positive regulatory roles in response to adverse conditions, the genes encoding thes proteins are considered potential candidates for breeding of plants for disease resistance and genetic improvement. However, there is far less information about this kinase family, and the function of these genes has not been explored in Rosaceae.

Results

A genome-wide survey and analysis of the genes encoding members of the SnRK2 family were performed in pear (Pyrus bretschneideri) and seven other Rosaceae species. A total of 71 SnRK2 genes were identified from the eight Rosaceae species and classified into three subgroups based on phylogenetic analysis and structural characteristics. Purifying selection played a crucial role in the evolution of SnRK2 genes, and whole-genome duplication and dispersed duplication were the primary forces underlying the characteristics of the SnRK2 gene family in Rosaceae. Transcriptome data and qRT-PCR assay results revealed that the distribution of PbrSnRK2s was very extensive, including across the roots, leaves, pollen, styles, and flowers, although most of them were mainly expressed in leaves. In addition, under stress conditions, the transcript levels of some of the genes were upregulated in leaves in response to ABA treatment.

Conclusions

This study provides useful information and a theoretical introduction for the study of the evolution, expression, and functions of the SnRK2 gene family in plants.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07201-w.

MIR2111-5 locus and shoot-accumulated mature miR2111 systemically enhance nodulation depending on HAR1 in Lotus japonicus

Legumes utilize a shoot-mediated signaling system to maintain a mutualistic relationship with nitrogen-fixing bacteria in root nodules. In Lotus japonicus, shoot-to-root transfer of microRNA miR2111 that targets TOO MUCH LOVE, a nodulation suppressor in roots, has been proposed to explain the mechanism underlying nodulation control from shoots. However, the role of shoot-accumulating miR2111s for the systemic regulation of nodulation was not clearly shown. Here, we find L. japonicus has seven miR2111 loci, including those mapped through RNA-seq. MIR2111-5 expression in leaves is the highest among miR2111 loci and repressed after rhizobial infection depending on a shoot-acting HYPERNODULATION ABERRANT ROOT FORMATION1 (HAR1) receptor. MIR2111-5 knockout mutants show significantly decreased nodule numbers and miR2111 levels. Furthermore, grafting experiments using transformants demonstrate scions with altered miR2111 levels influence nodule numbers in rootstocks in a dose-dependent manner. Therefore, miR2111 accumulation in leaves through MIR2111-5 expression is required for HAR1-dependent systemic optimization of nodule number.

Mining and evolution analysis of lateral organ boundaries domain (LBD) genes in Chinese white pear (Pyrus bretschneideri)

BACKGROUND

The lateral organ boundaries domain (LBD) gene is a plant-specific transcription factor that plays a critical role in diverse biological processes. However, the evolution and functional divergence of the LBD gene family has not yet been characterized for the Chinese White Pear.

RESULTS

In our study, a total of 60 PbrLBDs were identified in the pear genome. The PbrLBD gene family was divided into two classes based on gene structure and phylogenetic analysis: class I (53) and class II (7). Cis-acting element analysis results suggested that PbrLBDs may participate in various biological processes, such as flavonoid biosynthetic and stress response. Synteny analysis results indicated that segmental duplication played a key role in the expansion of the PbrLBD gene family. The mean Ks and 4DTv values showed that the PbrLBD gene family had undergone only one recent whole-genome duplication event occurring at 30-45 MYA. Purifying selection was a primary force during the PbrLBD gene family evolution process. Transcriptome data analysis revealed that 10 PbrLBDs were expressed in all six examined tissues, and 73.33% of members in the PbrLBD gene family were expressed in pear sepal. qRT-PCR was conducted to verify the expression levels of 11 PbrLBDs in these six tissues. Specifically, PbrLBD20, PbrLBD35 and PbrLBD53 genes were down-regulated when anthocyanin concentrations were high, whereas PbrLBD33 was significantly up-regulated in pear when anthocyanin concentrations were high. Furthermore, PbrLBD20, one of the candidate genes related to anthocyanins was localized in the nucleus.

CONCLUSIONS

Our analysis provides valuable information for understanding the evolution of the PbrLBD gene family, and provides new insights into the regulation of pear pigment metabolism and lays a foundation for the future disclosure of the molecular mechanism of LBD gene regulating flavonoid metabolism.

Genome-wide identification of lysin motif containing protein family genes in eight rosaceae species, and expression analysis in response to pathogenic fungus Botryosphaeria dothidea in Chinese white pear

BACKGROUND

Lysin motif-containing proteins (LYP), which act as pattern-recognition receptors, play central roles in growth, node formation, and responses to biotic stresses. The sequence of Chinese white pear genome (cv. 'Dangshansuli') along with the seven other species of Rosaceae has already been reported. Although, in these fruit crops, there is still a lack of clarity regarding the LYP family genes and their evolutionary history.

RESULTS

In the existing study, eight Rosaceae species i.e., Pyrus communis, Prunus persica, Fragaria vesca, Pyrus bretschneideri, Prunus avium, Prunus mume, Rubus occidentalis, and Malus × domestica were evaluated. Here, we determined a total of 124 LYP genes from the underlined Rosaceae species. While eighteen of the genes were from Chinese white pear, named as PbrLYPs. According to the LYPs structural characteristics and their phylogenetic analysis, those genes were classified into eight groups (group LYK1, LYK2, LYK3, LYK4/5, LYM1/3, LYM2, NFP, and WAKL). Dispersed duplication and whole-genome duplication (WGD) were found to be the most contributing factors of LYP family expansion in the Rosaceae species. More than half of the duplicated PbrLYP gene pairs were dated back to the ancient WGD (~ 140 million years ago (MYA)), and PbrLYP genes have experienced long-term purifying selection. The transcriptomic results indicated that the PbrLYP genes expression was tissue-specific. Most PbrLYP genes showed differential expression in leaves under fungal pathogen infection with two of them located in the plasmalemma.

CONCLUSION

A comprehensive analysis identified 124 LYP genes in eight Rosaceae species. Our findings have provided insights into the functions and characteristics of the Rosaceae LYP genes and a guide for the identification of other candidate LYPs for further genetic improvements for pathogen-resistance in higher plants.

Expansion of MIR482/2118 by a class-II transposable element in cotton

Some plant microRNA (miRNA) families contain multiple members generating identical or highly similar mature miRNA variants. Mechanisms underlying the expansion of miRNA families remain elusive, although tandem and/or segmental duplications have been proposed. In this study of two tetraploid cottons, Gossypium hirsutum and Gossypium barbadense, and their extant diploid progenitors, Gossypium arboreum and Gossypium raimondii, we investigated the gain and loss of members of the miR482/2118 superfamily, which modulates the expression of nucleotide-binding site leucine-rich repeat (NBS-LRR) disease resistance genes. We found significant expansion of MIR482/2118d in G. barbadense, G. hirsutum and G. raimondii, but not in G. arboreum. Several newly expanded MIR482/2118d loci have mutated to produce different miR482/2118 variants with altered target-gene specificity. Based on detailed analysis of sequences flanking these MIR482/2118 loci, we found that this expansion of MIR482/2118d and its derivatives resulted from an initial capture of an MIR482/2118d by a class-II DNA transposable element (TE) in G. raimondii prior to the tetraploidization event, followed by transposition to new genomic locations in G. barbadense, G. hirsutum and G. raimondii. The 'GosTE' involved in the capture and proliferation of MIR482/2118d and its derivatives belongs to the PIF/Harbinger superfamily, generating a 3-bp target site duplication upon insertion at new locations. All orthologous MIR482/2118 loci in the two diploids were retained in the two tetraploids, but mutation(s) in miR482/2118 were observed across all four species as well as in different cultivars of both G. barbadense and G. hirsutum, suggesting a dynamic co-evolution of miR482/2118 and its NBS-LRR targets. Our results provide fresh insights into the mechanisms contributing to MIRNA proliferation and enrich our knowledge on TEs.

The alcohol dehydrogenase gene family in sugarcane and its involvement in cold stress regulation

BACKGROUND

Alcohol dehydrogenases (ADHs) in plants are encoded by a multigene family. ADHs participate in growth, development, and adaptation in many plant species, but the evolution and function of the ADH gene family in sugarcane is still unclear.

RESULTS

In the present study, 151 ADH genes from 17 species including 32 ADH genes in Saccharum spontaneum and 6 ADH genes in modern sugarcane cultivar R570 were identified. Phylogenetic analysis demonstrated two groups of ADH genes and suggested that these genes underwent duplication during angiosperm evolution. Whole-genome duplication (WGD)/segmental and dispersed duplications played critical roles in the expansion of ADH family in S. spontaneum and R570, respectively. ScADH3 was cloned and preferentially expressed in response to cold stress. ScADH3 conferred improved cold tolerance in E. coli cells. Ectopic expression showed that ScADH3 can also enhance cold tolerance in transgenic tobacco. The accumulation of reactive oxygen species (ROS) in leaves of transgenic tobacco was significantly lower than in wild-type tobacco. The transcript levels of ROS-related genes in transgenic tobacco increased significantly. ScADH3 seems to affect cold tolerance by regulating the ROS-related genes to maintain the ROS homeostasis.

CONCLUSIONS

This study depicted the size and composition of the ADH gene family in 17 species, and investigated their evolution pattern. Comparative genomics analysis among the ADH gene families of S. bicolor, R570 and S. spontaneum revealed their close evolutionary relationship. Functional analysis suggested that ScADH3, which maintained the steady state of ROS by regulating ROS-related genes, was related to cold tolerance. These findings will facilitate research on evolutionary and functional aspects of the ADH genes in sugarcane, especially for the understanding of ScADH3 under cold stress.

Genome-Wide Analysis of the GRAS Gene Family in Barley (Hordeum vulgare L.)

The GRAS (named after first three identified proteins within this family, GAI, RGA, and SCR) family contains plant-specific genes encoding transcriptional regulators that play a key role in gibberellin (GA) signaling, which regulates plant growth and development. Even though GRAS genes have been characterized in some plant species, little research is known about the GRAS genes in barley (Hordeum vulgare L.). In this study, we observed 62 GRAS members from barley genome, which were grouped into 12 subgroups by using phylogenomic analysis together with the GRAS genes from Arabidopsis (Arabidopsis thaliana), maize (Zea mays), and rice (Oryza sativa). Chromosome localization and gene structure analysis suggested that duplication events and abundant presence of intronless genes might account for the massive expansion of GRAS gene family in barley. The analysis of RNA-seq data indicates the expression pattern of GRAS genes in various tissues at different stages in barley. Noteworthy, our qRT-PCR analysis showed the expression of 18 candidate GRAS genes abundantly in the developing inflorescence, indicating their potential roles in the barley inflorescence development and reproduction. Collectively, our evolutionary and expression analysis of GRAS family are useful for future functional characterization of GA signaling in barley and agricultural improvement.

Characterization of the Nitrate Transporter gene family and functional identification of HvNRT2.1 in barley (Hordeum vulgare L.)

Nitrogen use efficiency (NUE) is the efficiency with which plants acquire and use nitrogen. Plants have high-affinity nitrate transport systems, which involve certain nitrate transporter (NRT) genes. However, limited data are available on the contribution of the NRT2/3 gene family in barley nitrate transport. In the present study, ten putative NRT2 and three putative NRT3 genes were identified using bioinformatics methods. All the HvNRT2/3 genes were located on chromosomes 3H, 5H, 6H or 7H. Remarkably, the presence of tandem repeats indicated that duplication events contributed to the expansion of the NRT2 gene family in barley. In addition, the HvNRT2/3 genes displayed various expression patterns at selected developmental stages and were induced in the roots by both low and high nitrogen levels. Furthermore, the overexpression of HvNRT2.1 improved the yield related traits in Arabidopsis. Taken together, the data generated in the present study will be useful for genome-wide analyses to determine the precise role of the HvNRT2/3 genes during barley development, with the ultimate goal of improving NUE and crop production.

The Cytochrome P450 Monooxygenase Inventory of Grapevine (Vitis vinifera L.): Genome-Wide Identification, Evolutionary Characterization and Expression Analysis

The cytochrome P450 (CYP) monooxygenase superfamily, belonging to heme-thiolate protein products, plays a vital role in metabolizing physiologically valuable compounds in plants. To date, CYP superfamily genes have not yet been characterized in grapevine (V. vinifera L.), and their functions remain unclear. In this study, a sum of 236 VvCYPs, divided into 46 families and clustered into nine clans, have been identified based on bioinformatics analyses in grapevine genome. The characteristics of both exon-intron organizations and motif structures further supported the close evolutionary relationships of VvCYP superfamily as well as the reliability of phylogenetic analysis. The gene number-based hierarchical cluster of CYP subfamilies of different plants demonstrated that the loss of CYP families seems to be limited to single species or single taxa. Promoter analysis elucidated various cis-regulatory elements related to phytohormone signaling, plant growth and development, as well as abiotic/biotic stress responses. The tandem duplication mainly contributed to the expansion of the VvCYP superfamily, followed by singleton duplication in grapevine. Global RNA-sequencing data of grapevine showed functional divergence of VvCYPs as diverse expression patterns of VvCYPs in various organs, tissues, and developmental phases, which were confirmed by quantitative real-time reverse transcription PCR (qRT-PCR). Taken together, our results provided valuable inventory for understanding the classification and biological functions of the VvCYPs and paved the way for further functional verification of these VvCYPs and are helpful to grapevine molecular breeding.

Integrated Analysis of Large-Scale Omics Data Revealed Relationship Between Tissue Specificity and Evolutionary Dynamics of Small RNAs in Maize (Zea mays)

The evolutionary dynamics and tissue specificity of protein-coding genes are well documented in plants. However, the evolutionary consequences of small RNAs (sRNAs) on tissue-specific functions remain poorly understood. Here, we performed integrated analysis of 195 deeply sequenced sRNA libraries of maize B73, representing more than 10 tissues, and identified a comprehensive list of 419 maize microRNA (miRNA) genes, 271 of which were newly discovered in this study. We further characterized the evolutionary dynamics and tissue specificity of miRNA genes and corresponding miRNA isoforms (isomiRs). Our analysis revealed that tissue specificity of isomiR events tends to be associated with miRNA gene abundance and suggested that the frequencies of isomiR types are affected by the local genomic regions. Moreover, genome duplication (GD) events have dramatic effect on evolutionary dynamics of maize miRNA genes, and the abundance divergence for tissue-specific miRNA genes is associated with GD events. Further study indicated that duplicate miRNA genes with tissue-specific expression patterns, such as miR2275a, a phased siRNA (phasiRNA) trigger, contribute to phenotypic traits in maize. Additionally, our study revealed the expression preference of 21- and 24-nt phasiRNAs in relation to tissue specificity. This large-scale sRNAomic study depicted evolutionary implications of tissue-specific maize sRNAs, which coordinate genome duplication, isomiR modification, phenotypic traits and phasiRNAs differentiation.