A genome-wide analysis of the auxin/indole-3-acetic acid gene family in hexaploid bread wheat (Triticum aestivum L.)

Genome-Wide Analysis of Aux/IAA Gene Family in Artemisia argyi: Identification, Phylogenetic Analysis, and Determination of Response to Various Phytohormones

Artemisia argyi is a traditional herbal medicine plant, and its folium artemisia argyi is widely in demand due to moxibustion applications globally. The Auxin/indole-3-acetic acid (Aux/IAA, or IAA) gene family has critical roles in the primary auxin-response process, with extensive involvement in plant development and stresses, controlling various essential traits of plants. However, the systematic investigation of the Aux/IAA gene family in A. argyi remains limited. In this study, a total of 61 Aux/IAA genes were comprehensively identified and characterized. Gene structural analysis indicated that 46 Aux/IAA proteins contain the four typical domains, and 15 Aux/IAA proteins belong to non-canonical IAA proteins. Collinear prediction and phylogenetic relationship analyses suggested that Aux/IAA proteins were grouped into 13 distinct categories, and most Aux/IAA genes might experience gene loss during the tandem duplication process. Promoter cis-element investigation indicated that Aux/IAA promoters contain a variety of plant hormone response and stress response cis-elements. Protein interaction prediction analysis demonstrated that AaIAA26/29/7/34 proteins are possibly core members of the Aux/IAA family interaction. Expression analysis in roots and leaves via RNA-seq data indicated that the expression of some AaIAAs exhibited tissue-specific expression patterns, and some AaIAAs were involved in the regulation of salt and saline-alkali stresses. In addition, RT-qPCR results indicated that AaIAA genes have differential responses to auxin, with complex response patterns in response to other hormones, indicating that Aux/IAA may play a role in connecting auxin and other hormone signaling pathways. Overall, these findings shed more light on AaIAA genes and offer critical foundational knowledge toward the elucidation of their function during plant growth, stress response, and hormone networking of Aux/IAA family genes in A. argyi.

Genome-Wide Identification of the Paulownia fortunei Aux/IAA Gene Family and Its Response to Witches' Broom Caused by Phytoplasma

The typical symptom of Paulownia witches' broom (PaWB), caused by phytoplasma infection, is excessive branching, which is mainly triggered by auxin metabolism disorder. Aux/IAA is the early auxin-responsive gene that participates in regulating plant morphogenesis such as apical dominance, stem elongation, lateral branch development, and lateral root formation. However, no studies have investigated the response of the Aux/IAA gene family to phytoplasma infection in Paulownia fortunei. In this study, a total of 62 Aux/IAA genes were found in the genome. Phylogenetic analysis showed that PfAux/IAA genes could be divided into eight subgroups, which were formed by tandem duplication and fragment replication. Most of them had a simple gene structure, and several members lacked one or two conserved domains. By combining the expression of PfAux/IAA genes under phytoplasma stress and SA-treated phytoplasma-infected seedlings, we found that PfAux/IAA13/33/45 may play a vital role in the occurrence of PaWB. Functional analysis based on homologous relationships showed a strong correlation between PfAux/IAA45 and branching. Protein-protein interaction prediction showed that PfARF might be the binding partner of PfAux/IAA, and the yeast two-hybrid assay and bimolecular fluorescent complementary assay confirmed the interaction of PfAux/IAA45 and PfARF13. This study provides a theoretical basis for further understanding the function of the PfAux/IAA gene family and exploring the regulatory mechanism of branching symptoms caused by PaWB.

Genome-wide identification, expression analysis, and transcriptome analysis of the IAA gene family in Zoysia japonica

BACKGROUND

AUX/IAA is an essential signaling molecule and has great physiological importance in various plants, but its function in Zoysia japonica remains unknown.

METHODS AND RESULTS

Genome-wide identification and analysis of AUX/IAA genes used bioinformatics methods to investigate the ZjIAA genes' expression of exogenous IAA hydroponics treatment for 2 h by qRT-PCR, control and exogenous IAA treated zoysia were subjected to transcriptome sequencing. ZjIAAs were distributed across the 13 subfamilies by phylogenetic analysis with Oryza sativa and Arabidopsis thaliana. Multiple sequence alignment revealed that the majority of genes were non-canonical ZjIAAs with incomplete domain. The optimal growth concentration of the IAA hormone was 0.05 mM, and the qRT-PCR analysis revealed that eight ZjIAAs were differentially expressed, with seven genes considerably upregulating and one gene significantly downregulating. The result of transcriptome sequencing revealed that 515 differentially expressed genes (DEGs) were identified, with 344 upregulated genes and 171 downregulated genes. A total of 18 genes were annotated as involved in the plant hormone signal transduction pathway. And 8 ZjIAAs exhibited distinct expressions, 7 upregulated, and only one downregulated, according to the qRT-PCR study.

CONCLUSIONS

Genome-wide identification and analysis increased the understanding of the evolution and function of the IAA family in zoysia. DEGs of control and treatment with 0.05 mM exogenous IAA hormone were investigated by transcriptome sequencing. ZjIAAs had substantial variations in the expression of associated genes, with the majority of genes upregulated and 18 genes implicated in plant hormone signal transduction.

Genome-wide identification of Aux/IAA and ARF gene families in bread wheat (Triticum aestivum L.)

Wheat (Triticum aestivum L.) is one of the most important food crops in the world. Somatic embryogenesis is an event that is triggered by the presence of auxin hormone for the induction of somatic cells to get converted to embryonic cells. Somatic embryogenesis represents the most important process of totipotency of plants. The role of auxins is widely understood during various stages of embryogenesis including polarity establishment, de-differentiation, re-differentiations, and morphogenesis. Many of the Aux/IAAs and ARFs which are part of auxin signaling components have been identified to play various roles during embryogenesis. In this analysis, the Aux/IAAs and ARFs of T. aestivum have been analyzed at the genome-scale; their structure, function, and evolutionary relatedness were determined. Several Aux/IAAs and ARFs components of T. aestivum have been found to exclusively regulate axis formation, meristem commitment, and other re-differentiation processes by differential expression studies.

Genome-wide identification of the mitogen-activated protein kinase kinase kinase (MAPKKK) in pear (Pyrus bretschneideri) and their functional analysis in response to black spot

Identification of MAPKKK genes in pear and functional characterization of PbrMAPKKK82 in response to pear black spot. Mitogen-activated protein kinase kinase kinase (MAPKKK) is located upstream of the MAPK cascade pathway. This region senses extracellular stimuli via the signaling molecule or by themselves and is activated by phosphorylation. In this study, we identified 108 PbrMAPKKK genes from the pear genome. The genes were divided into three subfamilies and contained the conserved domain. Except for chromosome 7, there were 93 PbrMAPKKK genes randomly distributed on 16 out of the 17 chromosomes, while 15 PbrMAPKKK genes were detected on unknown chromosomes. They largely originated from whole-genome duplication (WGD) and dispersed events. In the expression analysis of PbrMAPKKK genes in seven pear tissue types by using a database, 20 PbrMAPKKK genes were selected to verify the expression associated with different resistance in two varieties by quantitative real-time PCR (qRT-PCR). The results showed that PbrMAPKKK12, PbrMAPKKK13, PbrMAPKKK53, PbrMAPKKK60, PbrMAPKKK65, PbrMAPKKK82, PbrMAPKKK83, and PbrMAPKKK96 were correlated with black spot resistance. PbrMAPKKK3, PbrMAPKKK9, PbrMAPKKK11, PbrMAPKKK34, PbrMAPKKK80, PbrMAPKKK81, PbrMAPKKK99, and PbrMAPKKK100 were correlated with black spot susceptibility, while the PbrMAPKKK gene positively responded to the life process of pear resistance to black spot. Furthermore, virus-induced gene silencing (VIGS) indicated that the PbrMAPKKK82 gene enhanced resistance to pear black spot disease.

Characterization of the Key Bibenzyl Synthase in Dendrobium sinense

Dendrobium sinense, an endemic medicinal herb in Hainan Island, is rich in bibenzyls. However, the key rate-limited enzyme involved in bibenzyl biosynthesis has yet to be identified in D. sinense. In this study, to explore whether there is a significant difference between the D. sinense tissues, the total contents of bibenzyls were determined in roots, pseudobulbs, and leaves. The results indicated that roots had higher bibenzyl content than pseudobulbs and leaves. Subsequently, transcriptomic sequencings were conducted to excavate the genes encoding type III polyketide synthase (PKS). A total of six D. sinense PKS (DsPKS) genes were identified according to gene function annotation. Phylogenetic analysis classified the type III DsPKS genes into three groups. Importantly, the c93636.graph_c0 was clustered into bibenzyl synthase (BBS) group, named as D. sinense BBS (DsBBS). The expression analysis by FPKM and RT-qPCR indicated that DsBBS showed the highest expression levels in roots, displaying a positive correlation with bibenzyl contents in different tissues. Thus, the recombinant DsBBS-HisTag protein was constructed and expressed to study its catalytic activity. The molecular weight of the recombinant protein was verified to be approximately 45 kDa. Enzyme activity analysis indicated that the recombinant DsBBS-HisTag protein could use 4-coumaryol-CoA and malonyl-CoA as substrates for resveratrol production in vitro. The Vmax of the recombinant protein for the resveratrol production was 0.88 ± 0.07 pmol s⁻¹ mg⁻¹. These results improve our understanding with respect to the process of bibenzyl biosynthesis in D. sinense.

Genome-wide identification of CCT genes in wheat (Triticum aestivum L.) and their expression analysis during vernalization

Numerous CCT genes are known to regulate various biological processes, such as circadian rhythm regulation, flowering, light signaling, plant development, and stress resistance. The CCT gene family has been characterized in many plants but remains unknown in the major cereal wheat (Triticum aestivum L.). Extended exposure to low temperature (vernalization) is necessary for winter wheat to flower successfully. VERNALIZATION2 (VRN2), a specific CCT-containing gene, has been proved to be strongly associated with vernalization in winter wheat. Mutation of all VRN2 copies in three subgenomes results in the eliminated demands of low temperature in flowering. However, no other CCT genes have been reported to be associated with vernalization to date. The present study screened CCT genes in the whole wheat genome, and preliminarily identified the vernalization related CCT genes through expression analysis. 127 CCT genes were identified in three subgenomes of common wheat through a hidden Markov model-based method. Based on multiple alignment, these genes were grouped into 40 gene clusters, including the duplicated gene clusters TaCMF6 and TaCMF8, each tandemly arranged near the telomere. The phylogenetic analysis classified these genes into eight groups. The transcriptome analysis using leaf tissues collected before, during, and after vernalization revealed 49 upregulated and 31 downregulated CCT genes during vernalization, further validated by quantitative real-time PCR. Among the differentially expressed and well-investigated CCT gene clusters analyzed in this study, TaCMF11, TaCO18, TaPRR95, TaCMF6, and TaCO16 were induced during vernalization but decreased immediately after vernalization, while TaCO1, TaCO15, TaCO2, TaCMF8, and TaPPD1 were stably suppressed during and after vernalization. These data imply that some vernalization related CCT genes other than VRN2 may exist in wheat. This study improves our understanding of CCT genes and provides a foundation for further research on CCT genes related to vernalization in wheat.

Toward synthetic plant development

The ability to engineer plant form will enable the production of novel agricultural products designed to tolerate extreme stresses, boost yield, reduce waste, and improve manufacturing practices. While historically, plants were altered through breeding to change their size or shape, advances in our understanding of plant development and our ability to genetically engineer complex eukaryotes are leading to the direct engineering of plant structure. In this review, I highlight the central role of auxin in plant development and the synthetic biology approaches that could be used to turn auxin-response regulators into powerful tools for modifying plant form. I hypothesize that recoded, gain-of-function auxin response proteins combined with synthetic regulation could be used to override endogenous auxin signaling and control plant structure. I also argue that auxin-response regulators are key to engineering development in nonmodel plants and that single-cell -omics techniques will be essential for characterizing and modifying auxin response in these plants. Collectively, advances in synthetic biology, single-cell -omics, and our understanding of the molecular mechanisms underpinning development have set the stage for a new era in the engineering of plant structure.

Silencing GhIAA43, a member of cotton AUX/IAA genes, enhances wilt resistance via activation of salicylic acid-mediated defenses

Auxin-mediated degradation of Aux/IAA proteins is a crucial step in auxin signaling. Recent researches indicate that Aux/IAA members also play a role in biotic and abiotic stresses. For example, Pseudomonas syringae infection causes Arabidopsis Aux/IAA protein (AXR2, AXR3) turnover. Here, by analyzing RNA-seq data we found that several cotton Aux/IAA genes are responsive to Verticillium dahliae infection, one of these named GhIAA43, was investigated for its role in cotton defense against V. dahliae infection. We demonstrate that the transcript levels of GhIAA43 were responsive to both V. dahliae infection and exogenous IAA application. By producing transgenic Arabidopsis plants overexpressing GhIAA43-GUS fusion, we show that IAA treatment and V. dahliae infection promoted GhIAA43 protein turnover. Silencing GhIAA43 in cotton enhanced wilt resistance, suggesting that GhIAA43 is a negative regulator in cotton defense against V. dahliae attack. By monitoring SA marker gene expression and measurement of SA content in GhIAA43-silenced cotton plants, we found that the enhanced resistance in GhIAA43-silenced cotton plants is due to the activation of SA-related defenses, and the activated defenses specifically occurred in the presence of V. dahliae. Furthermore, exogenous IAA application improve wilt resistance in cotton plants tested. Our results provide novel connection between auxin signaling and SA-related defenses in cotton upon V. dahliae attack.

Identification of GhLOG gene family revealed that GhLOG3 is involved in regulating salinity tolerance in cotton (Gossypium hirsutum L.)

Cytokinin (CK) is an important plant hormone that promotes plant cell division and differentiation, and participates in salt response under osmotic stress. LOGs (LONELY GUY) are CK-activating enzymes involved in CK synthesis. The LOG gene family has not been comprehensively characterized in cotton. In this study we identified 151 LOG genes from nine plant species, including 28 LOG genes in Gossypium hirsutum. Phylogenetic analysis divided LOG genes into three groups. Exon/intron structures and protein motifs of GhLOG genes were highly conserved. Synteny analysis revealed that several gene loci were highly conserved between the A and D sub-genomes of G. hirsutum with purifying selection pressure during evolution. Expression profiles showed that most LOG genes were constitutively expressed in eight different tissues. Furthermore, LOG genes can be regulated by abiotic stresses and phytohormone treatments. Moreover, subcellular localization revealed that GhLOG3_At resides inside the cell membrane. Overexpression of GhLOG3 enhanced salt tolerance in Arabidopsis. Virus-induced gene silencing (VIGS) of GhLOG3_At in cotton enhanced sensitivity of plants to salt stress with increased H₂O₂ contents and decreased chlorophyll and proline (PRO) activity. Our results suggested that GhLOG3_At induces salt stress tolerance in cotton, and provides a basis for the use of CK synthesis genes to regulate cotton growth and stress resistance.

Genome-wide identification and expression analysis of the NCED family in cotton (Gossypium hirsutum L.)

Abscisic acid (ABA) is an important plant hormone that plays multiple roles in regulating growth and development as well as in stress responses in plants. The NCED gene family includes key genes involved in the process of ABA synthesis. This gene family has been found in many species; however, the function of the NCED gene family in cotton is unclear. Here, a total of 23 NCED genes (designated as GhNCED1 to GhNCED23) were identified in cotton. Phylogenetic analysis indicated that the identified NCED proteins from cotton and Arabidopsis could be classified into 4 subgroups. Conserved motif analysis revealed that the gene structure and motif distribution of proteins within each subgroup were highly conserved. qRT-PCR and ABA content analyses indicated that NCED genes exhibited stage-specific expression patterns at tissue development stages. GhNCED5, GhNCED6 and GhNCED13 expression was similar to the change in ABA content, suggesting that this gene family plays a role in ABA synthesis. These results provide a better understanding of the potential functions of GhNCED genes.

Genome-Wide Identification, Evolution, and Expression Analysis of LBD Transcription Factor Family in Bread Wheat (Triticum aestivum L.)

The lateral organ boundaries domain (LBD) genes, as the plant-specific transcription factor family, play a crucial role in controlling plant architecture and stress tolerance. Although it has been thoroughly characterized in many species, the LBD family was not well studied in wheat. Here, the wheat LBD family was systematically investigated through an in silico genome-wide search method. A total of 90 wheat LBD genes (TaLBDs) were identified, which were classified into class I containing seven subfamilies, and class II containing two subfamilies. Exon-intron structure, conserved protein motif, and cis-regulatory elements analysis showed that the members in the same subfamily shared similar gene structure organizations, supporting the classification. Furthermore, the expression patterns of these TaLBDs in different types of tissues and under diverse stresses were identified through public RNA-seq data analysis, and the regulation networks of TaLBDs involved were predicted. Finally, the expression levels of 12 TaLBDs were validated by quantitative PCR (qPCR) analysis and the homoeologous genes showed differential expression. Additionally, the genetic diversity of TaLBDs in the landrace population showed slightly higher than that of the genetically improved germplasm population while obvious asymmetry at the subgenome level. This study not only provided the potential targets for further functional analysis but also contributed to better understand the roles of LBD genes in regulating development and stress tolerance in wheat and beyond.

Genome-wide characterization and expression analyses of the auxin/indole-3-acetic acid (Aux/IAA) gene family in apple (Malus domestica)

Auxin is a necessary phytohormone for fruit development, accompanying the whole process of fruit growth and development. The Aux/IAA gene family is one of the early auxin-responsive gene families. At present, there were few reports involved in Aux/IAA genes in the fruit, especially in apple. In our study, we identified 42 MdAux/IAAs, phylogenetic analysis showed that Aux/IAA proteins from apple, tomato, and strawberry were clustered into 5 groups, 42 MdAux/IAAs randomly distributed on 13 chromosomes. Additionally, a comprehensive analysis of Aux/IAA gene family was completed, including gene structures, conserved motifs, phylogenetic analysis, chromosome mapping, orthologous identification, selection pressure analyses, synteny analysis, and protein interaction. We also tested the expression of MdAux/IAAs in different tissues and fruit development stages using quantitative reverse transcription-polymerase chain reaction (qRT-PCR), we found that the most members of Aux/IAA showed higher expression in seeds compared within stem and leaves, indicating they may play a role in regulating fruit development. We also declared that the expression of Aux/IAA gene was not consistent in the pericarp and seeds at the same developmental stage, 3 MdAux/IAAs of the pericarp were upregulated over 20-fold at 90 d and 5 MdAux/IAAs of the seeds were upregulated over 40-fold at 90 d. It was MdAux/IAA23 that showed extreme up-regulated expression in both pericarp and seeds. This study proved that the Aux/IAA gene families may perform a different function in apple fruit development and ripening, more importantly, it provided a foundation for further exploring the biological function of these MdAux/IAAs.

Genome-wide characterization and expression analysis of geranylgeranyl diphosphate synthase genes in cotton (Gossypium spp.) in plant development and abiotic stresses

Background

GGPP (geranylgeranyl diphosphate) is produced in the isoprenoid pathway and mediates the function of various plant metabolites, which is synthesized by GGPPS (GGPP synthases) in plants. GGPPS characterization has not been performed in any plant species except Arabidopsis thaliana. Here, we performed a complete computational and bioinformatics analysis of GGPPS and detected their transcription expression pattern in Gossypium hirsutum for the first time so that to explore their evolutionary relationship and potential functions. Finally, we unravelled evolutionary relationship, conserved sequence logos, gene duplication and potential involvement in plant development and abiotic stresses tolerance of GGPPS genes in G. hirsutum and other plant species.

Results

A total of 159 GGPPS genes from 18 plant species were identified and evolutionary analysis divided these GGPPS genes into five groups to indicate their divergence from a common ancestor. Further, GGPPS family genes were conserved during evolution and underwent segmental duplication. The identified 25 GhGGPPS genes showed diverse expression pattern particularly in ovule and fiber development indicating their vital and divers roles in the fiber development. Additionally, GhGGPPS genes exhibited wide range of responses when subjected to abiotic (heat, cold, NaCl and PEG) stresses and hormonal (BL, GA, IAA, SA and MeJA) treatments, indicating their potential roles in various biotic and abiotic stresses tolerance.

Conclusions

The GGPPS genes are evolutionary conserved and might be involve in different developmental stages and stress response. Some potential key genes (e.g. GhGGPP4, GhGGPP9, and GhGGPP15) were suggested for further study and provided valuable source for cotton breeding to improve fiber quality and resistant to various stresses.

Genome-wide characterization and expression analyses of the auxin/indole-3-acetic acid (Aux/IAA) gene family in barley (Hordeum vulgare L.)

Aux/IAA genes are early auxin-responsive genes and essential for auxin signaling transduction. There is little information about Aux/IAAs in the agriculturally important cereal, barley. Using in silico method, we identified and subsequently characterized 36 Aux/IAAs from the barley genome. Based on their genomic sequences and the phylogenic relationship with Arabidopsis and rice Aux/IAA, the 36 HvIAAs were categorized into two major groups and 14 subgroups. The indication of the presence or absence of these domains for the biological functions and acting mechanisms was discussed. The cis-element distributions in HvIAA promoters suggests that the HvIAAs expressions may not only regulated by auxin (the presence of AuxREs and TGA-element) but also by other hormones and developmental and environmental cues. We then studied the HvIAAs expression in response to NAA (1-Naphthaleneacetic acid) using quantitative real-time PCR (qRT-PCR). Like the promoter analysis, only 14 HvIAAs were upregulated by NAA over two-fold at 4 h. HvIAAs were clustered into three groups based on the spatiotemporal expression data. We confirmed by qRT-PCR that most HvIAAs, especially HvIAA3, HvIAA7, HvIAA8, HvIAA18, HvIAA24 and HvIAA34, are expressed in the developing barley spike compared within seedling, suggesting their roles in regulating spike development. Taken together, our data provide a foundation for further revealing the biological function of these HvIAAs.

Identification and evolutionary analysis of polycistronic miRNA clusters in domesticated and wild wheat

MicroRNAs are ~22 nucleotide long non-coding RNAs that regulate gene expression at posttranscriptional level. Genome-wide analysis was performed to identify polycistronic miRNAs from wheat. Total 89 polycistronic miRNAs were identified in bread wheat which were distributed on three component sub-genomes (A = 26, B = 33 and D = 30). Except some, most of the identified polycistronic miRNAs were also present in other cultivated and wild wheat species. Expression of 11 identified polycistronic miRNAs could be validated using previously assembled transcriptomes, RNA-seq/s-RNA seq data of cultivated and wild wheats and RT-PCR. Polycistronic miRNAs orthologs were also localized on rice and Brachypodium genomes. As a case study, we also analyzed molecular evolution of miR395 family polycistrons in wheat. Both tandem and segmental duplications contributed to expansion of miR395 family polycistrons. Our findings provide a comprehensive view on wheat polycitronic miRNAs that will enable their in-depth functional analysis in the future.

Molecular functional analysis of auxin/indole-3-acetic acid proteins (Aux/IAAs) in plant disease resistance in cassava

Auxin/indole-3-acetic acid proteins (Aux/IAAs) play important roles in auxin signaling pathways, with extensive involvement in plant development and plant response to abiotic and biotic stresses. Manihot esculenta (Cassava) is one of the most important biomass energy crops in tropical regions; however, the information about Aux/IAA proteins remain limited in cassava. In this study, 37 MeAux/IAA gene family members were identified in cassava and a phylogenetic analysis was performed. The transcript levels of MeAux/IAAs were commonly regulated by the pathogen Xanthomonas axonopodis pv manihotis (Xam), and some of them were specifically localized to the nucleus. Moreover, the overexpression of MeAux/IAAs confers an improved disease resistance against Xam in Nicotiana benthamiana, while MeAux/IAAs-silenced plants show disease sensitivity against Xam in cassava, as evidenced by the leaf phenotype and leaf bacterial population. Consistent with the disease resistance, MeAux/IAAs regulated the transcript levels of PATHOGENESIS-RELATED GENES (MePRs), reactive oxygen species accumulation and callose development in the plants' defense response. Taken together, gene profile and functional analysis identified several MeAux/IAAs as novel members in plant disease resistance, providing important information for further utilization of MeAux/IAAs.

Genome-Wide Identification and Expression Profiles of Late Embryogenesis-Abundant (LEA) Genes during Grain Maturation in Wheat (Triticum aestivum L.)

Late embryogenesis-abundant (LEA) genes play important roles in plant growth and development, especially the cellular dehydration tolerance during seed maturation. In order to comprehensively understand the roles of LEA family members in wheat, we carried out a series of analyses based on the latest genome sequence of the bread wheat Chinese Spring. 121 Triticum aestivum L. LEA (TaLEA) genes, classified as 8 groups, were identified and characterized. TaLEA genes are distributed in all chromosomes, most of them with a low number of introns (≤3). Expression profiles showed that most TaLEA genes expressed specifically in grains. By qRT-PCR analysis, we confirmed that 12 genes among them showed high expression levels during late stage grain maturation in two spring wheat cultivars, Yangmai16 and Yangmai15. For most genes, the peak of expression appeared earlier in Yangmai16. Statistical analysis indicated that expression level of 8 genes in Yangmai 16 were significantly higher than Yangmai 15 at 25 days after anthesis. Taken together, our results provide more knowledge for future functional analysis and potential utilization of TaLEA genes in wheat breeding.

The roles of Aux/IAA gene family in development of Dendrocalamus sinicus (Poaceae: Bambusoideae) inferred by comprehensive analysis and expression profiling

Auxin is an important hormone in many plant developmental processes. In this study, the auxin/indole acetic acid (Aux/IAA) gene family was comprehensively identified using Dendrocalamus sinicus transcriptome data. A total of 26 Aux/IAA genes (DsIAA1-DsIAA26) were mined using four conserved Aux/IAA family motifs (PF02309). They encoded hydrophilic proteins, including one or two nuclear localisation signals. The D. sinicus Aux/IAA proteins were classified into two groups, including seven sister-gene pairs based on their phylogenetic relationships. A phylogenetic tree generated by aligning 108 predicted protein sequences of 26 DsIAAs, 43 PhIAAs (Phyllostachys heterocycla), 29 AtIAAs (Arabidopsis), 31 OsIAAs (Oryza sativa) and 22 PtIAAs (Populus) revealed nine major groups. Among them, four groups, including 96 IAA proteins of all five species, suggested that the genes originated before divergence of monocots and dicots. The expression profiling in different tissues showed that most of the DsIAAs preferentially expressed in leaves and shoots, suggesting their important roles in the development of leaves and shoots in D. sinicus. Continuously high expression of DsIAA3, DsIAA4, DsIAA15, and DsIAA20 may be important for regulating shoot development in D. sinicus. These results provide useful information for further research into the function of Aux/IAA genes in woody sympodial bamboos.

Genome-wide characterization and phylogenetic analysis of GSK gene family in three species of cotton: evidence for a role of some GSKs in fiber development and responses to stress

BACKGROUND

The glycogen synthase kinase 3/shaggy kinase (GSK3) is a serine/threonine kinase with important roles in animals. Although GSK3 genes have been studied for more than 30 years, plant GSK genes have been studied only since the last decade. Previous research has confirmed that plant GSK genes are involved in diverse processes, including floral development, brassinosteroid signaling, and responses to abiotic stresses.

RESULT

In this study, 20, 15 (including 5 different transcripts) and 10 GSK genes were identified in G. hirsutum, G. raimondii and G. arboreum, respectively. A total of 65 genes from Arabidopsis, rice, and cotton were classified into 4 clades. High similarities were found in GSK3 protein sequences, conserved motifs, and gene structures, as well as good concordance in gene pairwise comparisons (G. hirsutum vs. G. arboreum, G. hirsutum vs. G. raimondii, and G. arboreum vs. G. raimondii) were observed. Whole genome duplication (WGD) within At and Dt sub-genomes has been central to the expansion of the GSK gene family. Furthermore, GhSK genes showed diverse expression patterns in various tissues. Additionally, the expression profiles of GhSKs under different stress treatments demonstrated that many are stress-responsive genes. However, none were induced by brassinolide treatment. Finally, nine co-expression sub-networks were observed for GhSKs and the functional annotations of these genes suggested that some GhSKs might be involved in cotton fiber development.

CONCLUSION

In this present work, we identified 45 GSK genes from three cotton species, which were divided into four clades. The gene features, muti-alignment, conversed motifs, and syntenic blocks indicate that they have been highly conserved during evolution. Whole genome duplication was determined to be the dominant factor for GSK gene family expansion. The analysis of co-expressed sub-networks and tissue-specific expression profiles suggested functions of GhSKs during fiber development. Moreover, their different responses to various abiotic stresses indicated great functional diversity amongst the GhSKs. Briefly, data presented herein may serve as the basis for future functional studies of GhSKs.

Genome-Wide Identification and Transcriptional Regulation of Aquaporin Genes in Bread Wheat (Triticum aestivum L.) under Water Stress

Aquaporins (AQPs) are transmembrane proteins essential for controlling the flow of water and other molecules required for development and stress tolerance in plants, including important crop species such as wheat (Triticum aestivum). In this study, we utilized a genomic approach for analyzing the information about AQPs available in public databases to characterize their structure and function. Furthermore, we validated the expression of a suite of AQP genes, at the transcriptional level, including accessions with contrasting responses to drought, different organs and water stress levels. We found 65 new AQP genes, from which 60% are copies expanded by polyploidization. Sequence analysis of the AQP genes showed that the purifying selection pressure acted on duplicate genes, which was related to a high conservation of the functions. This situation contrasted with the expression patterns observed for different organs, developmental stages or genotypes under water deficit conditions, which indicated functional divergence at transcription. Expression analyses on contrasting genotypes showed high gene transcription from Tonoplast Intrinsic Protein 1 (TIP1) and 2 (TIP2), and Plasma Membrane Intrinsic Protein 1 (PIP1) and 2 (PIP2) subfamilies in roots and from TIP1 and PIP1 subfamilies in leaves. Interestingly, during severe drought stress, 4 TIP genes analyzed in leaves of the tolerant accession reached up to 15-fold the level observed at the susceptible genotype, suggesting a positive relationship with drought tolerance. The obtained results extend our understanding of the structure and function of AQPs, particularly under water stress conditions.

Genome-wide sequence and expressional analysis of autophagy Gene family in bread wheat (Triticum aestivum L.)

Autophagy, a highly conserved intracellular degradation system, is regarded to be responsible for self-defense and protect cells from abiotic stress. Extensive studies have demonstrated that autophagy plays a crucial role in regulating plant growth and development as well as in response to diverse stresses. However, little is known about autophagy-associated genes (ATGs) in wheat, especially those involved in the regulatory network of stress processes. In this study, a total of 108 putative wheat ATGs (TaATG) were obtained based on a genome-wide search approach. Phylogenetic analysis classified them into 13 subfamilies, of which the TaAtg16 subfamily consisted of 29 members, ranking it the largest subfamily. The conserved motif compositions as well as their exon-intron structures were systematically analyzed and strongly supported the classification. The homoeologous genes tended to have similar gene features during wheat polyploidization. Furthermore, a total of 114 putative cis-elements were found, and those related to hormone, stress, and light responsiveness were abundantly presented in the promoter regions. Co-expression network analysis revealed that orthologous VAMP727 was the hub node of the whole network, and complex interactions were also found. Finally, the expression profiles of TaATGs among different tissues and under abiotic stresses were investigated to identify tissue-specific or stress-responsive candidates, and then 14 were validated by wet-lab analysis. Results showed that the TaAtg8 subfamily played a crucial role in tissue autophagy and stress defense, which could be considered as processes that are candidates for further functional study. This was the first study to comprehensively investigate the ATG family in wheat, which ultimately provided important clues for further functional analysis and also took a step toward uncovering the evolutionary mechanism of ATG genes in wheat and beyond.

Genome-Wide Identification and Expression Analysis of the HD-Zip Gene Family in Wheat (Triticum aestivum L.)

The homeodomain-leucine zipper (HD-Zip) gene family, as plant-specific transcription factors, plays an important role in plant development and growth as well as in the response to diverse stresses. Although HD-Zip genes have been extensively studied in many plants, they had not yet been studied in wheat, especially those involved in response to abiotic stresses. In this study, 46 wheat HD-Zip genes were identified using a genome-wide search method. Phylogenetic analysis classified these genes into four groups, numbered 4, 5, 17 and 20 respectively. In total, only three genes with A, B and D homoeologous copies were identified. Furthermore, the gene interaction networks found that the TaHDZ genes played a critical role in the regulatory pathway of organ development and osmotic stress. Finally, the expression profiles of the wheat HD-Zips in different tissues and under various abiotic stresses were investigated using the available RNA sequencing (RNA-Seq) data and then validated by quantitative real-time polymerase chain reaction (qRT-PCR) to obtain the tissue-specific and stress-responsive candidates. This study systematically identifies the HD-Zip gene family in wheat at the genome-wide level, providing important candidates for further functional analysis and contributing to the better understanding of the molecular basis of development and stress tolerance in wheat.

Aux/IAA Gene Family in Plants: Molecular Structure, Regulation, and Function

Auxin plays a crucial role in the diverse cellular and developmental responses of plants across their lifespan. Plants can quickly sense and respond to changes in auxin levels, and these responses involve several major classes of auxin-responsive genes, including the Auxin/Indole-3-Acetic Acid (Aux/IAA) family, the auxin response factor (ARF) family, small auxin upregulated RNA (SAUR), and the auxin-responsive Gretchen Hagen3 (GH3) family. Aux/IAA proteins are short-lived nuclear proteins comprising several highly conserved domains that are encoded by the auxin early response gene family. These proteins have specific domains that interact with ARFs and inhibit the transcription of genes activated by ARFs. Molecular studies have revealed that Aux/IAA family members can form diverse dimers with ARFs to regulate genes in various ways. Functional analyses of Aux/IAA family members have indicated that they have various roles in plant development, such as root development, shoot growth, and fruit ripening. In this review, recently discovered details regarding the molecular characteristics, regulation, and protein-protein interactions of the Aux/IAA proteins are discussed. These details provide new insights into the molecular basis of the Aux/IAA protein functions in plant developmental processes.

Characterization and Expression Patterns of Auxin Response Factors in Wheat

Auxin response factors (ARFs) are important transcription factors involved in both the auxin signaling pathway and the regulatory development of various plant organs. In this study, 23 TaARF members encoded by a total of 68 homeoalleles were isolated from 18 wheat chromosomes (excluding chromosome 4). The TaARFs, including their conserved domains, exon/intron structures, related microRNAs, and alternative splicing (AS) variants, were then characterized. Phylogenetic analysis revealed that members of the TaARF family share close homology with ARFs in other grass species. qRT-PCR analyses revealed that 20 TaARF members were expressed in different organs and tissues and that the expression of some members significantly differed in the roots, stems, and leaves of wheat seedlings in response to exogenous auxin treatment. Moreover, protein network analyses and co-expression results showed that TaTIR1-TaARF15/18/19-TaIAA13 may interact at both the protein and genetic levels. The results of subsequent evolutionary analyses showed that three transcripts of TaARF15 in the A subgenome of wheat exhibited high evolutionary rate and underwent positive selection. Transgenic analyses indicated that TaARF15-A.1 promoted the growth of roots and leaves of Arabidopsis thaliana and was upregulated in the overexpression plants after auxin treatment. Our results will provide reference information for subsequent research and utilization of the TaARF gene family.

Genome-wide analysis of the auxin/indoleacetic acid (Aux/IAA) gene family in allotetraploid rapeseed (Brassica napus L.)

BACKGROUND

Auxin/Indoleacetic acid (Aux/IAA) genes participate in the auxin signaling pathway and play key roles in plant growth and development. Although the Aux/IAA gene family has been identified in many plants, within allotetraploid Brassica napus little is known.

RESULTS

In this study, a total of 119 Aux/IAA genes were found in the genome of B. napus. They were distributed non-randomly across all 19 chromosomes and other non-anchored random scaffolds, with a symmetric distribution in the A and C subgenomes. Evolutionary and comparative analysis revealed that 111 (94.1%) B. napus Aux/IAA genes were multiplied due to ancestral Brassica genome triplication and recent allotetraploidy from B. rapa and B. oleracea. Phylogenetic analysis indicated seven subgroups containing 29 orthologous gene sets and two Brassica-specific gene sets. Structures of genes and proteins varied across different genes but were conserved among homologous genes in B. napus. Furthermore, analysis of transcriptional profiles revealed that the expression patterns of Aux/IAA genes in B. napus were tissue dependent. Auxin-responsive elements tend to be distributed in the proximal region of promoters, and are significantly associated with early exogenous auxin up-regulation.

CONCLUSIONS

Members of the Aux/IAA gene family were identified and analyzed comprehensively in the allotetraploid B. napus genome. This analysis provides a deeper understanding of diversification of the Aux/IAA gene family and will facilitate further dissection of Aux/IAA gene function in B. napus.

Evolution of the Aux/IAA Gene Family in Hexaploid Wheat

The Aux/IAA (IAA) gene family, involved in the auxin signalling pathway, acts as an important regulator in plant growth and development. In this study, we explored the evolutionary trajectory of the IAA family in common wheat. The results showed ten pairs of paralogs among 34 TaIAA family members. Seven of the pairs might have undergone segmental duplication, and the other three pairs appear to have experienced tandem duplication. Except for TaIAA15-16, these duplication events occurred in the ancestral genomes before the divergence of Triticeae. After that point, two polyploidization events shaped the current TaIAA family consisting of three subgenomic copies. The structure or expression pattern of the TaIAA family begins to differentiate in the hexaploid genome, where TaIAAs in the D genome lost more genes (eight) and protein secondary structures (α1, α3 and β5) than did the other two genomes. Expression analysis showed that six members of the TaIAA family were not expressed, and members such as TaIAA8, 15, 16, 28 and 33 exhibited tissue-specific expression patterns. In addition, three of the ten pairs of paralogs (TaIAA5-12, TaIAA15-16 and TaIAA29-30) showed similar expression patterns, and another five paralog pairs displayed differential expression patterns. Phylogenetic analysis showed that paralog pairs with high rates of evolution (ω > ω ₀), particularly TaIAA15-16 and TaIAA29-30, experienced greater motif loss, with only zero to two interacting IAA proteins. In contrast, most paralogous genes with low ω, such as TaIAA5-12, had more complete motifs and higher degrees of interaction with other family members.

Genome-wide identification and analysis of the MADS-box gene family in bread wheat (Triticum aestivum L.)

The MADS-box genes encode transcription factors with key roles in plant growth and development. A comprehensive analysis of the MADS-box gene family in bread wheat (Triticum aestivum) has not yet been conducted, and our understanding of their roles in stress is rather limited. Here, we report the identification and characterization of the MADS-box gene family in wheat. A total of 180 MADS-box genes classified as 32 Mα, 5 Mγ, 5 Mδ, and 138 MIKC types were identified. Evolutionary analysis of the orthologs among T. urartu, Aegilops tauschii and wheat as well as homeologous sequences analysis among the three sub-genomes in wheat revealed that gene loss and chromosomal rearrangements occurred during and/or after the origin of bread wheat. Forty wheat MADS-box genes that were expressed throughout the investigated tissues and development stages were identified. The genes that were regulated in response to both abiotic stresses (i.e., phosphorus deficiency, drought, heat, and combined drought and heat) and biotic stresses (i.e., Fusarium graminearum, Septoria tritici, stripe rust and powdery mildew) were detected as well. A few notable MADS-box genes were specifically expressed in a single tissue and those showed relatively higher expression differences between the stress and control treatment. The expression patterns of considerable MADS-box genes differed from those of their orthologs in Brachypodium, rice, and Arabidopsis. Collectively, the present study provides new insights into the possible roles of MADS-box genes in response to stresses and will be valuable for further functional studies of important candidate MADS-box genes.

Genome-wide characterization of JASMONATE-ZIM DOMAIN transcription repressors in wheat (Triticum aestivum L.)

BACKGROUND

The JASMONATE-ZIM DOMAIN (JAZ) repressor family proteins are jasmonate co-receptors and transcriptional repressor in jasmonic acid (JA) signaling pathway, and they play important roles in regulating the growth and development of plants. Recently, more and more researches on JAZ gene family are reported in many plants. Although the genome sequencing of common wheat (Triticum aestivum L.) and its relatives is complete, our knowledge about this gene family remains vacant.

RESULTS

Fourteen JAZ genes were identified in the wheat genome. Structural analysis revealed that the TaJAZ proteins in wheat were as conserved as those in other plants, but had structural characteristics. By phylogenetic analysis, all JAZ proteins from wheat and other plants were clustered into 11 sub-groups (G1-G11), and TaJAZ proteins shared a high degree of similarity with some JAZ proteins from Aegliops tauschii, Brachypodium distachyon and Oryza sativa. The Ka/Ks ratios of TaJAZ genes ranged from 0.0016 to 0.6973, suggesting that the TaJAZ family had undergone purifying selection in wheat. Gene expression patterns obtained by quantitative real-time PCR (qRT-PCR) revealed differential temporal and spatial regulation of TaJAZ genes under multifarious abiotic stress treatments of high salinity, drought, cold and phytohormone. Among these, TaJAZ7, 8 and 12 were specifically expressed in the anther tissues of the thermosensitive genic male sterile (TGMS) wheat line BS366 and normal control wheat line Jing411. Compared with the gene expression patterns in the normal wheat line Jing411, TaJAZ7, 8 and 12 had different expression patterns in abnormally dehiscent anthers of BS366 at the heading stage 6, suggesting that specific up- or down-regulation of these genes might be associated with the abnormal anther dehiscence in TGMS wheat line.

CONCLUSION

This study analyzed the size and composition of the JAZ gene family in wheat, and investigated stress responsive and differential tissue-specific expression profiles of each TaJAZ gene in TGMS wheat line BS366. In addition, we isolated 3 TaJAZ genes that would be more likely to be involved in the regulation of abnormal anther dehiscence in TGMS wheat line. In conclusion, the results of this study contributed some novel and detailed information about JAZ gene family in wheat, and also provided 3 potential candidate genes for improving the TGMS wheat line.

Genome-Wide Investigation of WRKY Transcription Factors Involved in Terminal Drought Stress Response in Common Bean

WRKY transcription factor plays a key role in drought stress. However, the characteristics of the WRKY gene family in the common bean (Phaseolus vulgaris L.) are unknown. In this study, we identified 88 complete WRKY proteins from the draft genome sequence of the "G19833" common bean. The predicted genes were non-randomly distributed in all chromosomes. Basic information, amino acid motifs, phylogenetic tree and the expression patterns of PvWRKY genes were analyzed, and the proteins were classified into groups 1, 2, and 3. Group 2 was further divided into five subgroups: 2a, 2b, 2c, 2d, and 2e. Finally, we detected 19 WRKY genes that were responsive to drought stress using qRT-PCR; 11 were down-regulated, and 8 were up-regulated under drought stress. This study comprehensively examines WRKY proteins in the common bean, a model food legume, and it provides a foundation for the functional characterization of the WRKY family and opportunities for understanding the mechanisms of drought stress tolerance in this plant.

Genome-wide identification, phylogeny and expressional profiles of mitogen activated protein kinase kinase kinase (MAPKKK) gene family in bread wheat (Triticum aestivum L.)

BACKGROUND

Mitogen-activated protein kinase kinase kinases (MAPKKKs) are the important components of MAPK cascades, which play the crucial role in plant growth and development as well as in response to diverse stresses. Although this family has been systematically studied in many plant species, little is known about MAPKKK genes in wheat (Triticum aestivum L.), especially those involved in the regulatory network of stress processes.

RESULTS

In this study, we identified 155 wheat MAPKKK genes through a genome-wide search method based on the latest available wheat genome information, of which 29 belonged to MEKK, 11 to ZIK and 115 to Raf subfamily, respectively. Then, chromosome localization, gene structure and conserved protein motifs and phylogenetic relationship as well as regulatory network of these TaMAPKKKs were systematically investigated and results supported the prediction. Furthermore, a total of 11 homologous groups between A, B and D sub-genome and 24 duplication pairs among them were detected, which contributed to the expansion of wheat MAPKKK gene family. Finally, the expression profiles of these MAPKKKs during development and under different abiotic stresses were investigated using the RNA-seq data. Additionally, 10 tissue-specific and 4 salt-responsive TaMAPKKK genes were selected to validate their expression level through qRT-PCR analysis.

CONCLUSIONS

This study for the first time reported the genome organization, evolutionary features and expression profiles of the wheat MAPKKK gene family, which laid the foundation for further functional analysis of wheat MAPKKK genes, and contributed to better understanding the roles and regulatory mechanism of MAPKKKs in wheat.

Glutathione S-Transferase Gene Family in Gossypium raimondii and G. arboreum: Comparative Genomic Study and their Expression under Salt Stress

Glutathione S-transferases (GSTs) play versatile functions in multiple aspects of plant growth and development. A comprehensive genome-wide survey of this gene family in the genomes of G. raimondii and G. arboreum was carried out in this study. Based on phylogenetic analyses, the GST gene family of both two diploid cotton species could be divided into eight classes, and approximately all the GST genes within the same subfamily shared similar gene structure. Additionally, the gene structures between the orthologs were highly conserved. The chromosomal localization analyses revealed that GST genes were unevenly distributed across the genome in both G. raimondii and G. arboreum. Tandem duplication could be the major driver for the expansion of GST gene families. Meanwhile, the expression analysis for the selected 40 GST genes showed that they exhibited tissue-specific expression patterns and their expression were induced or repressed by salt stress. Those findings shed lights on the function and evolution of the GST gene family in Gossypium species.