Genome-wide identification and expression analysis of auxin response factor gene family in Medicago truncatula

Large-scale analysis of the ARF and Aux/IAA gene families in 406 horticultural and other plants

The auxin response factor (ARF) and auxin/indole-3-acetic acid (Aux/IAA) family of genes are central components of the auxin signaling pathway and play essential roles in plant growth and development. Their large-scale analysis and evolutionary trajectory of origin are currently not known. Here, we identified the corresponding ARF and Aux/IAA family members and performed a large-scale analysis by scanning 406 plant genomes. The results showed that the ARF and Aux/IAA gene families originated from charophytes. The ARF family sequences were more conserved than the Aux/IAA family sequences. Dispersed duplications were the common expansion mode of ARF and Aux/IAA families in bryophytes, ferns, and gymnosperms; however, whole-genome duplication was the common expansion mode of the ARF and Aux/IAA families in basal angiosperms, magnoliids, monocots, and dicots. Expression and regulatory network analyses revealed that the Arabidopsis thaliana ARF and Aux/IAA families responded to multiple hormone, biotic, and abiotic stresses. The APETALA2 and serum response factor-transcription factor gene families were commonly enriched in the upstream and downstream genes of the ARF and Aux/IAA gene families. Our study provides a comprehensive overview of the evolutionary trajectories, structural functions, expansion mechanisms, expression patterns, and regulatory networks of these two gene families.

Genome-wide survey, molecular evolution and expression analysis of Auxin Response Factor (ARF) gene family indicating their key role in seed number per pod in pigeonpea (C. cajan L. Millsp.)

Auxin Response Factors (ARF) are a family of transcription factors that mediate auxin signalling and regulate multiple biological processes. Their crucial role in increasing plant biomass/yield influenced this study, where a systematic analysis of ARF gene family was carried out to identify the key proteins controlling embryo/seed developmental pathways in pigeonpea. A genome-wide scan revealed the presence of 12 ARF genes in pigeonpea, distributed across the chromosomes 1, 3, 4, 8 and 11. Domain analysis of ARF proteins showed the presence of B3 DNA binding, AUX response, and IAA domains. Majority of them are of nuclear origin, and do not exhibit the level of genomic expansion as observed in Glycine max (51 members). The duplication events seem to range from 31.6 to 42.3 million years ago (mya). Promoter analysis revealed the presence of multiple cis-acting elements related to stress responses, hormone signalling and other development processes. The expression atlas data highlighted the expression of CcARF8 in hypocotyl, bud and flower whereas, CcARF7 expression was significantly high in pod. The real-time expression of CcARF2, CcARF3 and CcARF18 was highest in genotypes with high seed number indicating their key role in regulating embryo development and determining seed set in pigeonpea.

Analysis of Auxin-Encoding Gene Family in Vigna radiata and It's Cross-Species Expression Modulating Waterlogging Tolerance in Wild Vigna umbellata

Mungbean is known to be susceptible to waterlogging (WL) stress. Some of the wild species have the potential to tolerate this through various physiological and molecular mechanisms. Auxin Response Factor (ARF) and Auxin/Indole Acetic Acid (AUX/IAA), an early responsive gene family, has multiple functions in growth, development, and stress tolerance. Here, we report the first comprehensive analysis of the ARF and AUX/IAA gene family in mungbean. A total of 26 ARF and 19 AUX/IAA genes were identified from the mungbean genome. The ARF and AUX/IAA candidates were clearly grouped into two major clades. Further, the subgrouping within the major clades indicated the presence of significant diversity. The gene structure, motif analysis, and protein characterization provided the clue for further fundamental research. Out of the10 selected candidate genes, VrARF-5, VrARF-11, VrARF-25, and VrAUX/IAA-9 were found to significantly multiple-fold gene expression in the hypocotyl region of WL-tolerant wild relatives (PRR 2008-2) provides new insight into a role in the induction of lateral root formation under WL stress. The analysis provides an insight into the structural diversity of ARF and AUX/IAA genes in mungbean. These results increase our understanding of ARF and AUX/IAA genes and therefore offer robust information for functional investigations, which can be taken up in the future and will form a foundation for improving tolerance against waterlogging stress.

Identification of potential Auxin Response Candidate genes for soybean rapid canopy coverage through comparative evolution and expression analysis

Glycine max, soybean, is an abundantly cultivated crop worldwide. Efforts have been made over the past decades to improve soybean production in traditional and organic agriculture, driven by growing demand for soybean-based products. Rapid canopy cover development (RCC) increases soybean yields and suppresses early-season weeds. Genome-wide association studies have found natural variants associated with RCC, however causal mechanisms are unclear. Auxin modulates plant growth and development and has been implicated in RCC traits. Therefore, modulation of auxin regulatory genes may enhance RCC. Here, we focus on the use of genomic tools and existing datasets to identify auxin signaling pathway RCC candidate genes, using a comparative phylogenetics and expression analysis approach. We identified genes encoding 14 TIR1/AFB auxin receptors, 61 Aux/IAA auxin co-receptors and transcriptional co-repressors, and 55 ARF auxin response factors in the soybean genome. We used Bayesian phylogenetic inference to identify soybean orthologs of Arabidopsis thaliana genes, and defined an ortholog naming system for these genes. To further define potential auxin signaling candidate genes for RCC, we examined tissue-level expression of these genes in existing datasets and identified highly expressed auxin signaling genes in apical tissues early in development. We identified at least 4 TIR1/AFB, 8 Aux/IAA, and 8 ARF genes with highly specific expression in one or more RCC-associated tissues. We hypothesize that modulating the function of these genes through gene editing or traditional breeding will have the highest likelihood of affecting RCC while minimizing pleiotropic effects.

Identification and Evolutionary Analysis of the Auxin Response Factor (ARF) Family Based on Transcriptome Data from Caucasian Clover and Analysis of Expression Responses to Hormones

Caucasian clover (Trifolium ambiguum M. Bieb.) is an excellent perennial plant in the legume family Fabaceae, with a well-developed rhizome and strong clonal growth. Auxin is one of the most important phytohormones in plants and plays an important role in plant growth and development. Auxin response factor (ARF) can regulate the expression of auxin-responsive genes, thus participating in multiple pathways of auxin transduction signaling in a synergistic manner. No genomic database has been established for Caucasian clover. In this study, 71 TaARF genes were identified through a transcriptomic database of Caucasian clover rhizome development. Phylogenetic analysis grouped the TaARFs into six (1-6) clades. Thirty TaARFs contained a complete ARF structure, including three relatively conserved regions. Physical and chemical property analysis revealed that TaARFs are unstable and hydrophilic proteins. We also analyzed the expression pattern of TaARFs in different tissues (taproot, horizontal rhizome, swelling of taproot, rhizome bud and rhizome bud tip). Quantitative real-time RT-PCR revealed that all TaARFs were responsive to phytohormones (indole-3-acetic acid, gibberellic acid, abscisic acid and methyl jasmonate) in roots, stems and leaves. These results helped elucidate the role of ARFs in responses to different hormone treatments in Caucasian clover.

Genome-wide identification and expression analysis of the Auxin-Response factor (ARF) gene family in Medicago sativa under abiotic stress

BACKGROUND

Alfalfa (Medicago sativa) is the most widely planted legume forage and one of the most economically valuable crops in the world. The periodic changes in its growth and development and abiotic stress determine its yield and economic benefits. Auxin controls many aspects of alfalfa growth by regulating gene expression, including organ differentiation and stress response. Auxin response factors (ARF) are transcription factors that play an essential role in auxin signal transduction and regulate the expression of auxin-responsive genes. However, the function of ARF transcription factors is unclear in autotetraploid-cultivated alfalfa.

RESULT

A total of 81 ARF were identified in the alfalfa genome in this study. Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were analyzed, identifying that ARF genes are mainly involved in transcriptional regulation and plant hormone signal transduction pathways. Phylogenetic analysis showed that MsARF was divided into four clades: I, II, III, and IV, each containing 52, 13, 7, and 9 genes, respectively. The promoter region of the MsARF gene contained stress-related elements, such as ABRE, TC-rich repeats, MBS, LTR. Proteins encoded by 50 ARF genes were localized in the nucleus without guide peptides, signal peptides, or transmembrane structures, indicating that most MsARF genes are not secreted or transported but only function in the nucleus. Protein structure analysis revealed that the secondary and tertiary structures of the 81 MsARF genes varied. Chromosomal localization analysis showed 81 MsARF genes were unevenly distributed on 25 chromosomes, with the highest distribution on chromosome 5. Furthermore, 14 segmental duplications and two sets of tandem repeats were identified. Expression analysis indicated that the MsARF was differentially expressed in different tissues and under various abiotic stressors. The quantitative reverse transcription polymerase chain reaction (qRT-PCR) analysis showed that the expression profiles of 23 MsARF genes were specific to abiotic stresses such as drought, salt, high temperature, and low temperature, as well as tissue-specific and closely related to the duration of stress.

CONCLUSION

This study identified MsARF in the cultivated alfalfa genome based on the autotetraploid level, which GO, KEGG analysis, phylogenetic analysis, sequence characteristics, and expression pattern analysis further confirmed. Together, these findings provide clues for further investigation of MsARF functional verification and molecular breeding of alfalfa. This study provides a novel approach to systematically identify and characterize ARF transcription factors in autotetraploid cultivated alfalfa, revealing 23 MsARF genes significantly involved in response to various stresses.

Combined genomic and transcriptomic analysis reveals the contribution of tandem duplication genes to low-temperature adaptation in perennial ryegrass

Perennial ryegrass (Lolium perenne L.) is an agronomically important cool-season grass species that is widely used as forage for ruminant animal production and cultivated in temperate regions for the establishment of lawns. However, the underlying genetic mechanism of the response of L. perenne to low temperature is still unclear. In the present study, we performed a comprehensive study and identified 3,770 tandem duplication genes (TDGs) in L. perenne, and evolutionary analysis revealed that L. perenne might have undergone a duplication event approximately 7.69 Mya. GO and KEGG pathway functional analyses revealed that these TDGs were mainly enriched in photosynthesis, hormone-mediated signaling pathways and responses to various stresses, suggesting that TDGs contribute to the environmental adaptability of L. perenne. In addition, the expression profile analysis revealed that the expression levels of TDGs were highly conserved and significantly lower than those of all genes in different tissues, while the frequency of differentially expressed genes (DEGs) from TDGs was much higher than that of DEGs from all genes in response to low-temperature stress. Finally, in-depth analysis of the important and expanded gene family indicated that the members of the ELIP subfamily could rapidly respond to low temperature and persistently maintain higher expression levels during all low temperature stress time points, suggesting that ELIPs most likely mediate low temperature responses and help to facilitate adaptation to low temperature in L. perenne. Our results provide evidence for the genetic underpinning of low-temperature adaptation and valuable resources for practical application and genetic improvement for stress resistance in L. perenne.

Identification and transcriptome data analysis of ARF family genes in five Orchidaceae species

The Orchidaceae is a large family of perennial herbs especially noted for the exceptional diversity of specialized flowers. Elucidating the genetic regulation of flowering and seed development of orchids is an important research goal with potential utility in orchid breeding programs. Auxin Response Factor (ARF) genes encode auxin-responsive transcription factors, which are involved in the regulation of diverse morphogenetic processes, including flowering and seed development. However, limited information on the ARF gene family in the Orchidaceae is available. In this study, 112 ARF genes were identified in the genomes of 5 orchid species (Apostasia shenzhenica, Dendrobium catenatum, Phalaenopsis aphrodite, Phalaenopsis equestris and Vanilla planifolia,). These genes were grouped into 7 subfamilies based on their phylogenetic relationships. Compared with the ARF family in model plants, such as Arabidopsis thaliana and Oryza sativa, one group of ARF genes involved in pollen wall synthesis has been lost during evolution of the Orchidaceae. This loss corresponds with absence of the exine in the pollinia. Through mining of the published genomic and transcriptomic data for the 5 orchid species: the ARF genes of subfamily 4 may play an important role in flower formation and plant growth, whereas those of subfamily 3 are potentially involved in pollen wall development. the study results provide novel insights into the genetic regulation of unique morphogenetic phenomena of orchids, which lay a foundation for further analysis of the regulatory mechanisms and functions of sexual reproduction-related genes in orchids.

Histone demethylase MaJMJ15 is involved in the regulation of postharvest banana fruit ripening

Histone methylation plays important roles in plant development. However, the role of histone methylation in fruit ripening remains unclear. Here, a total of 16 Jumonji domain-containing proteins (JMJs) were identified from banana genome. During fruit ripening, expression of MaJMJ15 was significantly upregulated. Exogenous ethylene accelerated the upregulation whereas 1-methylcyclopropene delayed the process, suggesting that MaJMJ15 positively regulates banana fruit ripening. MaJMJ15 is an H3K27me3 site-specific demethylase. Transient overexpression of MaJMJ15 promoted banana fruit ripening. Moreover, the global H3K27me3 was decreased by MaJMJ15. Furthermore, MaJMJ15 directly targeted several key ripening-related genes (RRGs) in banana including NAC transcription factor 1/2 (MaNAC1/2), 1-aminocyclopropane-1-carboxylate synthase 1 (MaACS1), 1-aminocyclopropane-1-carboxylate oxidase 1 (MaACO1) and expansin 2 (MaEXP2), removed H3K27me3 from their chromatin, and activated their expression. Our data suggest that MaJMJ15 is an H3K27me3 demethylase, which is involved in the regulation of banana fruit ripening by activating expression of key RRGs via removal of H3K27me3.

Phylogeny, gene structures, and expression patterns of the auxin response factor (GhARF2) in upland cotton (Gossypium hirsutum L.)

BACKGROUND

Auxin response factors (ARFs) are a class of transcription factors that regulate the expression of auxin-responsive genes and play important functions in plant growth and development. To understand the biological functions of the auxin response factor GhARF2 gene in upland cotton, the coding sequence (CDS) of GhARF2 gene was cloned, and its protein sequence, evolutionary relationship, subcellular localization and expression pattern were analysed.

METHODS

The CDS sequence of GhARF2 gene was cloned from upland cotton variety Baimian No.1, and its protein sequence was analyzed by bioinformatics method. The subcellular localization of GhARF2 protein was detected by tobacco epidermal transient transformation system, and the tissue expression and stress expression pattern of GhARF2 were analyzed by quantitative Real‑Time PCR (qRT-PCR).

RESULTS

The full-length CDS of GhARF2 gene was 2583 bp, encoded 860 amino acids, and had a molecular weight and an isoelectric point of 95.46 KDa and 6.02, respectively. The GhARF2 protein had multiple phosphorylation sites, no transmembrane domain, and secondary structures dominated by random coils and alpha helix. The GhARF2 protein had 3 conserved typical domains of ARF gene family members, including the B3 DNA binding domain, the Auxin_resp domain, and the Aux/IAA domain. Phylogenetic analysis revealed that ARF2 proteins in different species were clustered in the Group A subgroup, in which GhARF2 was closely related to TcARF2 of Theobroma cacao L. (Malvaceae). The subcellular localization results showed that the GhARF2 protein was localized in the nucleus. Analysis of tissue expression pattern showed that the GhARF2 gene was expressed in all tested tissues, with the highest expression levels in sepal, followed by leaf, and the lowest expression levels in fiber. Further stress expression analysis showed that the GhARF2 gene was induced by drought, high-temperature, low-temperature and salt stress, and had different expression patterns under different stress conditions.

CONCLUSION

These results established a foundation for understanding the functions of GhARF2 and breeding varieties with high-stress tolerance in cotton.

Genome-Wide Identification and Characterization of Auxin Response Factor (ARF) Gene Family Involved in Wood Formation and Response to Exogenous Hormone Treatment in Populus trichocarpa

Auxin is a key regulator that virtually controls almost every aspect of plant growth and development throughout its life cycle. As the major components of auxin signaling, auxin response factors (ARFs) play crucial roles in various processes of plant growth and development. In this study, a total of 35 PtrARF genes were identified, and their phylogenetic relationships, chromosomal locations, synteny relationships, exon/intron structures, cis-elements, conserved motifs, and protein characteristics were systemically investigated. We also analyzed the expression patterns of these PtrARF genes and revealed that 16 of them, including PtrARF1, 3, 7, 11, 13-17, 21, 23, 26, 27, 29, 31, and 33, were preferentially expressed in primary stems, while 15 of them, including PtrARF2, 4, 6, 9, 10, 12, 18-20, 22, 24, 25, 28, 32, and 35, participated in different phases of wood formation. In addition, some PtrARF genes, with at least one cis-element related to indole-3-acetic acid (IAA) or abscisic acid (ABA) response, responded differently to exogenous IAA and ABA treatment, respectively. Three PtrARF proteins, namely PtrARF18, PtrARF23, and PtrARF29, selected from three classes, were characterized, and only PtrARF18 was a transcriptional self-activator localized in the nucleus. Moreover, Y2H and bimolecular fluorescence complementation (BiFC) assay demonstrated that PtrARF23 interacted with PtrIAA10 and PtrIAA28 in the nucleus, while PtrARF29 interacted with PtrIAA28 in the nucleus. Our results provided comprehensive information regarding the PtrARF gene family, which will lay some foundation for future research about PtrARF genes in tree development and growth, especially the wood formation, in response to cellular signaling and environmental cues.

Cell-specific pathways recruited for symbiotic nodulation in the Medicago truncatula legume

Medicago truncatula is a model legume species that has been studied for decades to understand the symbiotic relationship between legumes and soil bacteria collectively named rhizobia. This symbiosis called nodulation is initiated in roots with the infection of root hair cells by the bacteria, as well as the initiation of nodule primordia from root cortical, endodermal, and pericycle cells, leading to the development of a new root organ, the nodule, where bacteria fix and assimilate the atmospheric dinitrogen for the benefit of the plant. Here, we report the isolation and use of the nuclei from mock and rhizobia-inoculated roots for the single nuclei RNA-seq (sNucRNA-seq) profiling to gain a deeper understanding of early responses to rhizobial infection in Medicago roots. A gene expression map of the Medicago root was generated, comprising 25 clusters, which were annotated as specific cell types using 119 Medicago marker genes and orthologs to Arabidopsis cell-type marker genes. A focus on root hair, cortex, endodermis, and pericycle cell types, showing the strongest differential regulation in response to a short-term (48 h) rhizobium inoculation, revealed not only known genes and functional pathways, validating the sNucRNA-seq approach, but also numerous novel genes and pathways, allowing a comprehensive analysis of early root symbiotic responses at a cell type-specific level.

Genome-wide identification of the auxin response factor (ARF) gene family in Magnolia sieboldii and functional analysis of MsARF5

Auxin plays an essential role in flowering, embryonic development, seed dormancy, and germination. Auxin response factors (ARFs) are plant-specific key transcriptional factors in mediating the gene expression network of auxin signaling. Although ARFs in model plants such as Arabidopsis had been well characterized, their identities and potential roles in non-model plants are less studied. Here, we performed genome-wide identification of ARFs in Magnolia sieboldii K. Koch, a primitive species with high taxonomic importance and medicinal values. We found 25 ARF genes in M. sieboldii, which were widely distributed across multiple chromosomes. Based on sequence similarity, the encoded proteins could be either transcriptional repressors or activators. Gene expression analysis showed a dynamic pattern for many ARFs including MsARF5 during seed germination. In addition, overexpressing of MsARF5 showed that it restores many developmental defects in the Arabidopsis mutant. Moreover, two phenotypically distinct transgenic Arabidopsis lines were obtained, indicating a link between gene expression levels and developmental phenotypes. Taken together, we provided a systematic investigation of the ARF gene family in M. sieboldii and revealed an important role of MsARF5 in mediating auxin signaling.

Bioinformatics Study of Aux/IAA Family Genes and Their Expression in Response to Different Hormones Treatments during Japanese Apricot Fruit Development and Ripening

Auxin/indole-3-acetic acid (Aux/IAA) is a transcriptional repressor in the auxin signaling pathway that plays a role in several plant growth and development as well as fruit and embryo development. However, it is unclear what role they play in Japanese apricot (Prunus mume) fruit development and maturity. To investigate the role of Aux/IAA genes in fruit texture, development, and maturity, we comprehensively identified and expressed 19 PmIAA genes, and demonstrated their conserved domains and homology across species. The majority of PmIAA genes are highly responsive and expressed in different hormone treatments. PmIAA2, PmIAA5, PmIAA7, PmIAA10, PmIAA13, PmIAA18, and PmIAA19 showed a substantial increase in expression, suggesting that these genes are involved in fruit growth and maturity. During fruit maturation, alteration in the expression of PmIAA genes in response to 1-Methylcyclopropene (1-MCP) treatment revealed an interaction between auxin and ethylene. The current study investigated the response of Aux/IAA development regulators to auxin during fruit ripening, with the goal of better understanding their potential application in functional genomics.

Systematic Identification and Expression Analysis of the Auxin Response Factor (ARF) Gene Family in Ginkgo biloba L

Auxin participates in various physiological and molecular response-related developmental processes and is a pivotal hormone that regulates phenotypic formation in plants. Auxin response factors (ARFs) are vital transcription factors that mediate downstream auxin signaling by explicitly binding to auxin-responsive genes' promoters. Here, to investigate the possible developmental regulatory functions of ARFs in Ginkgo biloba, through employing comprehensive bioinformatics, we recognized 15 putative GbARF members. Conserved domains and motifs, gene and protein structure, gene duplication, GO enrichment, transcriptome expression profiles, and qRT-PCR all showed that Group I and III members were highly conserved. Among them, GbARF10b and GbARF10a were revealed as transcriptional activators in the auxin response for the development of Ginkgo male flowers through sequences alignment, cis-elements analysis and GO annotation; the results were corroborated for the treatment of exogenous SA. Moreover, the GbARFs expansion occurred predominantly by segmental duplication, and most GbARFs have undergone purifying selection. The Ka/Ks ratio test identified the functional consistence of GbARF2a and GbARF2c, GbARF10b, and GbARF10a in tissue expression profiles and male flower development. In summary, our study established a new research basis for exploring Ginkgo GbARF members' roles in floral organ development and hormone response.

Physio-Morphological, Biochemical and Transcriptomic Analyses Provide Insights Into Drought Stress Responses in Mesona chinensis Benth

Drought stress affects the normal growth and development of Mesona chinensis Benth (MCB), which is an important medicinal and edible plant in China. To investigate the physiological and molecular mechanisms of drought resistance in MCB, different concentrations of polyethylene glycol 6000 (PEG6000) (0, 5, 10, and 15%) were used to simulate drought conditions in this study. Results showed that the growth of MCB was significantly limited under drought stress conditions. Drought stress induced the increases in the contents of Chla, Chlb, Chla + b, soluble protein, soluble sugar, and soluble pectin and the activities of superoxide dismutase (SOD), catalase (CAT), total antioxidant capacity (TAC), hydrogen peroxide (H2O2), and malondialdehyde (MDA). Transcriptome analysis revealed 3,494 differentially expressed genes (DEGs) (1,961 up-regulated and 1,533 down-regulated) between the control and 15% PEG6000 treatments. These DEGs were identified to be involved in the 10 metabolic pathways, including "plant hormone signal transduction," "brassinosteroid biosynthesis," "plant-pathogen interaction," "MAPK signaling pathway-plant," "starch and sucrose metabolism," "pentose and glucuronate interconversions," "phenylpropanoid biosynthesis," "galactose metabolism," "monoterpenoid biosynthesis," and "ribosome." In addition, transcription factors (TFs) analysis showed 8 out of 204 TFs, TRINITY_DN3232_c0_g1 [ABA-responsive element (ABRE)-binding transcription factor1, AREB1], TRINITY_DN4161_c0_g1 (auxin response factor, ARF), TRINITY_DN3183_c0_g2 (abscisic acid-insensitive 5-like protein, ABI5), TRINITY_DN28414_c0_g2 (ethylene-responsive transcription factor ERF1b, ERF1b), TRINITY_DN9557_c0_g1 (phytochrome-interacting factor, PIF3), TRINITY_DN11435_c1_g1, TRINITY_DN2608_c0_g1, and TRINITY_DN6742_c0_g1, were closely related to the "plant hormone signal transduction" pathway. Taken together, it was inferred that these pathways and TFs might play important roles in response to drought stress in MCB. The current study provided important information for MCB drought resistance breeding in the future.

Genome-Wide cis-Regulatory Element Based Discovery of Auxin-Responsive Genes in Higher Plant

Auxin has a profound impact on plant physiology and participates in almost all aspects of plant development processes. Auxin exerts profound pleiotropic effects on plant growth and differentiation by regulating the auxin response genes’ expressions. The classical auxin reaction is usually mediated by auxin response factors (ARFs), which bind to the auxin response element (AuxRE) in the promoter region of the target gene. Experiments have generated only a limited number of plant genes with well-characterized functions. It is still unknown how many genes respond to exogenous auxin treatment. An economical and effective method was proposed for the genome-wide discovery of genes responsive to auxin in a model plant, Arabidopsis thaliana (A. thaliana). Our method relies on cis-regulatory-element-based targeted gene finding across different promoters in a genome. We first exploit and analyze auxin-specific cis-regulatory elements for the transcription of the target genes, and then identify putative auxin responsive genes whose promoters contain the elements in the collection of over 25,800 promoters in the A. thaliana genome. Evaluating our result by comparing with a published database and the literature, we found that this method has an accuracy rate of 65.2% (309/474) for predicting candidate genes responsive to auxin. Chromosome distribution and annotation of the putative auxin-responsive genes predicted here were also mined. The results can markedly decrease the number of identified but merely potential auxin target genes and also provide useful clues for improving the annotation of gene that lack functional information.

Genome-Wide Analysis and Expression Profile of Nuclear Factor Y (NF-Y) Gene Family in Z. jujuba

Nuclear factor-Y (NF-Y) is an important transcription factor in the plant species, which potentially provides a higher level of functional diversity including for abiotic stress tolerance. The genome-wide study and expression analysis of NF-Y gene family in Ziziphus, an elite abiotic stress-tolerant species, assist bioprospecting of genes. Here, a total of 32 NF-Y (8 NF-YA, 15 NF-YB, and 9 NF-YC) genes were identified in genome-wide search of Z. jujuba genome. Physicochemical properties, cellular localization, gene structure, chromosomal location, and protein motifs were analyzed for structural and functional understanding. Identified 12 NF-Ys were responsible for the expansion of NF-Y gene family by tandem duplication in Z. jujuba. Phylogenetic and comparative physical mapping of Z. jujuba NF-Ys with its orthologs illustrated evolutionary and functional insights into NF-Y gene family. A total of 45 perfect microsatellites (20bp to 40bp) were extracted across the ZjNF-Y genes. The promoter and gene ontology study suggested that Z. jujuba NF-Y gene family is functionally diverse and could play a wide-ranging role in plant abiotic stress, development, and cellular processes. An expression study revealed that large numbers of the NF-Ys are differentially expressed in response to drought and salinity. The total 15 and 18 ZjNF-Y genes that are upregulated under drought and salinity stress, respectively, are the potential candidates for further functional analysis for development of climate-resilient crops. The present study established a base for understanding the role of NF-Ys in Z. jujuba under abiotic stress conditions and paved a way for further research.

Genome-Wide Identification and Expression Analysis of the Aux/IAA and Auxin Response Factor Gene Family in Medicago truncatula

Aux/IAA and auxin response transcription factor (ARF) genes are key regulators of auxin responses in plants. A total of 25 MtIAA and 40 MtARF genes were identified based on the latest updated Medicago truncatula reference genome sequence. They were clustered into 10 and 8 major groups, respectively. The homologs among M. truncatula, soybean, and Arabidopsis thaliana shared close relationships based on phylogenetic analysis. Gene structure analysis revealed that MtIAA and MtARF genes contained one to four concern motifs and they are localized to eight chromosomes, except chromosome 6 without MtARFs. In addition, some MtIAA and MtARF genes were expressed in all tissues, while others were specifically expressed in specific tissues. Analysis of cis-acting elements in promoter region and expression profiles revealed the potential response of MtIAA and MtARF genes to hormones and abiotic stresses. The prediction protein–protein interaction network showed that some ARF proteins could interact with multiple Aux/IAA proteins, and the reverse is also true. The investigation provides valuable, basic information for further studies on the biological functions of MtIAA and MtARF genes in the regulation of auxin-related pathways in M. truncatula.

Plant-Microbe Interaction: Aboveground to Belowground, from the Good to the Bad

Soil health and fertility issues are constantly addressed in the agricultural industry. Through the continuous and prolonged use of chemical heavy agricultural systems, most agricultural lands have been impacted, resulting in plateaued or reduced productivity. As such, to invigorate the agricultural industry, we would have to resort to alternative practices that will restore soil health and fertility. Therefore, in recent decades, studies have been directed towards taking a Magellan voyage of the soil rhizosphere region, to identify the diversity, density, and microbial population structure of the soil, and predict possible ways to restore soil health. Microbes that inhabit this region possess niche functions, such as the stimulation or promotion of plant growth, disease suppression, management of toxicity, and the cycling and utilization of nutrients. Therefore, studies should be conducted to identify microbes or groups of organisms that have assigned niche functions. Based on the above, this article reviews the aboveground and below-ground microbiomes, their roles in plant immunity, physiological functions, and challenges and tools available in studying these organisms. The information collected over the years may contribute toward future applications, and in designing sustainable agriculture.

Genome-Wide Identification of ARF Transcription Factor Gene Family and Their Expression Analysis in Sweet Potato

Auxin response factors (ARFs) are a family of transcription factors that play an important role of auxin regulation through their binding with auxin response elements. ARF genes are represented by a large multigene family in plants; however, to our knowledge, the ARF gene family has not been well studied and characterized in sweet potatoes. In this study, a total of 25 ARF genes were identified in Ipomea trifida. The identified ItrARF genes’ conserved motifs, chromosomal locations, phylogenetic relationships, and their protein characteristics were systemically investigated using different bioinformatics tools. The expression patterns of ItfARF genes were analyzed within the storage roots and normal roots at an early stage of development. ItfARF16b and ItfARF16c were both highly expressed in the storage root, with minimal to no expression in the normal root. ItfARF6a and ItfARF10a exhibited higher expression in the normal root but not in the storage root. Subsequently, ItfARF1a, ItfARF2b, ItfARF3a, ItfARF6b, ItfARF8a, ItfARF8b, and ItfARF10b were expressed in both root types with moderate to high expression for each. All ten of these ARF genes and their prominent expression signify their importance within the development of each respective root type. This study provides comprehensive information regarding the ARF family in sweet potatoes, which will be useful for future research to discover further functional verification of these ItfARF genes.

Small RNAs, Degradome, and Transcriptome Sequencing Provide Insights into Papaya Fruit Ripening Regulated by 1-MCP

As an inhibitor of ethylene receptors, 1-methylcyclopropene (1-MCP) can delay the ripening of papaya. However, improper 1-MCP treatment will cause a rubbery texture in papaya. Understanding of the underlying mechanism is still lacking. In the present work, a comparative sRNA analysis was conducted after different 1-MCP treatments and identified a total of 213 miRNAs, of which 44 were known miRNAs and 169 were novel miRNAs in papaya. Comprehensive functional enrichment analysis indicated that plant hormone signal pathways play an important role in fruit ripening. Through the comparative analysis of sRNAs and transcriptome sequencing, a total of 11 miRNAs and 12 target genes were associated with the ethylene and auxin signaling pathways. A total of 1741 target genes of miRNAs were identified by degradome sequencing, and nine miRNAs and eight miRNAs were differentially expressed under the ethylene and auxin signaling pathways, respectively. The network regulation diagram of miRNAs and target genes during fruit ripening was drawn. The expression of 11 miRNAs and 12 target genes was verified by RT-qPCR. The target gene verification showed that cpa-miR390a and cpa-miR396 target CpARF19-like and CpERF RAP2-12-like, respectively, affecting the ethylene and auxin signaling pathways and, therefore, papaya ripening.

Genome-wide analysis of ARF transcription factors reveals HcARF5 expression profile associated with the biosynthesis of β-ocimene synthase in Hedychium coronarium

Herein, 37 ARF genes were identified and analyzed in Hedychium coronarium and HcARF5 showed a potential role in the regulation of HcTPS3. Auxin is an important plant hormone, implicated in various aspects of plant growth and development processes especially in the biosynthesis of various secondary metabolites. Auxin response factors (ARF) belong to the transcription factors (TFs) gene family and play a crucial role in transcriptional activation/repression of auxin-responsive genes by directly binding to their promoter region. Nevertheless, whether ARF genes are involved in the regulatory mechanism of volatile compounds in flowering plants is largely unknown. β-ocimene is a key floral volatile compound synthesized by terpene synthase 3 (HcTPS3) in Hedychium coronarium. A comprehensive analysis of H. coronarium genome reveals 37 candidate ARF genes in the whole genome. Tissue-specific expression patterns of HcARFs family members were assessed using available transcriptome data. Among them, HcARF5 showed a higher expression level in flowers, and significantly correlated with the key structural β-ocimene synthesis gene (HcTPS3). Furthermore, transcript levels of both genes were associated with the flower development. Under hormone treatments, the response of HcARF5 and HcTPS3, and the emission level of β-ocimene contents were evaluated. Subcellular and transcriptional activity assay showed that HcARF5 localizes to the nucleus and possesses transcriptional activity. Yeast one-hybrid (Y1H) and dual-luciferase assays revealed that HcARF5 directly regulates the transcriptional activity of HcTPS3. Yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) assays showed that HcARF5 interacts with scent-related HcIAA4, HcIAA6, and HcMYB1 in vivo. Overall, these results indicate that HcARF5 is potentially involved in the regulation of β-ocimene synthesis in H. coronarium.

Genome-wide analysis of ARF transcription factors reveals HcARF5 expression profile associated with the biosynthesis of β-ocimene synthase in Hedychium coronarium

Herein, 37 ARF genes were identified and analyzed in Hedychium coronarium and HcARF5 showed a potential role in the regulation of HcTPS3. Auxin is an important plant hormone, implicated in various aspects of plant growth and development processes especially in the biosynthesis of various secondary metabolites. Auxin response factors (ARF) belong to the transcription factors (TFs) gene family and play a crucial role in transcriptional activation/repression of auxin-responsive genes by directly binding to their promoter region. Nevertheless, whether ARF genes are involved in the regulatory mechanism of volatile compounds in flowering plants is largely unknown. β-ocimene is a key floral volatile compound synthesized by terpene synthase 3 (HcTPS3) in Hedychium coronarium. A comprehensive analysis of H. coronarium genome reveals 37 candidate ARF genes in the whole genome. Tissue-specific expression patterns of HcARFs family members were assessed using available transcriptome data. Among them, HcARF5 showed a higher expression level in flowers, and significantly correlated with the key structural β-ocimene synthesis gene (HcTPS3). Furthermore, transcript levels of both genes were associated with the flower development. Under hormone treatments, the response of HcARF5 and HcTPS3, and the emission level of β-ocimene contents were evaluated. Subcellular and transcriptional activity assay showed that HcARF5 localizes to the nucleus and possesses transcriptional activity. Yeast one-hybrid (Y1H) and dual-luciferase assays revealed that HcARF5 directly regulates the transcriptional activity of HcTPS3. Yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) assays showed that HcARF5 interacts with scent-related HcIAA4, HcIAA6, and HcMYB1 in vivo. Overall, these results indicate that HcARF5 is potentially involved in the regulation of β-ocimene synthesis in H. coronarium.

AUXIN RESPONSE FACTOR 18-HISTONE DEACETYLASE 6 module regulates floral organ identity in rose (Rosa hybrida)

The phytohormone auxin plays a pivotal role in floral meristem initiation and gynoecium development, but whether and how auxin controls floral organ identity remain largely unknown. Here, we found that auxin levels influence organ specification, and changes in auxin levels influence homeotic transformation between petals and stamens in rose (Rosa hybrida). The PIN-FORMED-LIKES (PILS) gene RhPILS1 governs auxin levels in floral buds during floral organogenesis. RhAUXIN RESPONSE FACTOR 18 (RhARF18), whose expression decreases with increasing auxin content, encodes a transcriptional repressor of the C-class gene RhAGAMOUS (RhAG), and controls stamen-petal organ specification in an auxin-dependent manner. Moreover, RhARF18 physically interacts with the histone deacetylase (HDA) RhHDA6. Silencing of RhHDA6 increases H3K9/K14 acetylation levels at the site adjacent to the RhARF18-binding site in the RhAG promoter and reduces petal number, indicating that RhARF18 might recruit RhHDA6 to the RhAG promoter to reinforce the repression of RhAG transcription. We propose a model for how auxin homeostasis controls floral organ identity via regulating transcription of RhAG.

Auxin-induced AUXIN RESPONSE FACTOR4 activates APETALA1 and FRUITFULL to promote flowering in woodland strawberry

Flowering time is known to be regulated by numerous pathways, such as the autonomous, gibberellin, aging, photoperiod-mediated, and vernalization pathways. These regulatory mechanisms involve both environmental triggers and endogenous hormonal cues. Additional flowering control mechanisms mediated by other phytohormones, such as auxin, are less well understood. We found that in cultivated strawberry (Fragaria × ananassa), the expression of auxin response factor4 (FaARF4) was higher in the flowering stage than in the vegetative stage. Overexpression of FaARF4 in Arabidopsis thaliana and woodland strawberry (Fragaria vesca) resulted in transgenic plants flowering earlier than control plants. In addition, FveARF4-silenced strawberry plants showed delayed flowering compared to control plants, indicating that FaARF4 and FveARF4 function similarly in regulating flowering. Further studies showed that ARF4 can bind to the promoters of the floral meristem identity genes APETALA1 (AP1) and FRUITFULL (FUL), inducing their expression and, consequently, flowering in woodland strawberry. Our studies reveal an auxin-mediated flowering pathway in strawberry involving the induction of ARF4 expression.

Auxin Response Factor 2 (ARF2), ARF3, and ARF4 Mediate Both Lateral Root and Nitrogen Fixing Nodule Development in Medicago truncatula

Auxin Response Factors (ARFs) constitute a large family of transcription factors that mediate auxin-regulated developmental programs in plants. ARF2, ARF3, and ARF4 are post-transcriptionally regulated by the microRNA390 (miR390)/trans-acting small interference RNA 3 (TAS3) module through the action of TAS3-derived trans - acting small interfering RNAs (ta-siRNA). We have previously reported that constitutive activation of the miR390/TAS3 pathway promotes elongation of lateral roots but impairs nodule organogenesis and infection by rhizobia during the nitrogen-fixing symbiosis established between Medicago truncatula and its partner Sinorhizobium meliloti. However, the involvement of the targets of the miR390/TAS3 pathway, i.e., MtARF2, MtARF3, MtARF4a, and MtARF4b, in root development and establishment of the nitrogen-fixing symbiosis remained unexplored. Here, promoter:reporter fusions showed that expression of both MtARF3 and MtARF4a was associated with lateral root development; however, only the MtARF4a promoter was active in developing nodules. In addition, up-regulation of MtARF2, MtARF3, and MtARF4a/b in response to rhizobia depends on Nod Factor perception. We provide evidence that simultaneous knockdown of MtARF2, MtARF3, MtARF4a, and MtARF4b or mutation in MtARF4a impaired nodule formation, and reduced initiation and progression of infection events. Silencing of MtARF2, MtARF3, MtARF4a, and MtARF4b altered mRNA levels of the early nodulation gene nodulation signaling pathway 2 (MtNSP2). In addition, roots with reduced levels of MtARF2, MtARF3, MtARF4a, and MtARF4b, as well as arf4a mutant plants exhibited altered root architecture, causing a reduction in primary and lateral root length, but increasing lateral root density. Taken together, our results suggest that these ARF members are common key players of the morphogenetic programs that control root development and the formation of nitrogen-fixing nodules.

Alternative Splicing Generates a MONOPTEROS Isoform Required for Ovule Development

The plant hormone auxin is a fundamental regulator of organ patterning and development that regulates gene expression via the canonical AUXIN RESPONSE FACTOR (ARF) and AUXIN/INDOLE-3-ACETIC ACID (Aux/IAA) combinatorial system. ARF and Aux/IAA factors interact, but at high auxin concentrations, the Aux/IAA transcriptional repressor is degraded, allowing ARF-containing complexes to activate gene expression. ARF5/MONOPTEROS (MP) is an important integrator of auxin signaling in Arabidopsis development and activates gene transcription in cells with elevated auxin levels. Here, we show that in ovules, MP is expressed in cells with low levels of auxin and can activate the expression of direct target genes. We identified and characterized a splice variant of MP that encodes a biologically functional isoform that lacks the Aux/IAA interaction domain. This MP11ir isoform was able to complement inflorescence, floral, and ovule developmental defects in mp mutants, suggesting that it was fully functional. Our findings describe a novel scenario in which ARF post-transcriptional regulation controls the formation of an isoform that can function as a transcriptional activator in regions of subthreshold auxin concentration.

Auxin Profiling and GmPIN Expression in Phytophthora sojae-Soybean Root Interactions

Auxin (indole-3-acetic acid, IAA) has been implicated as a susceptibility factor in both beneficial and pathogenic molecular plant-microbe interactions. Previous studies have identified a large number of auxin-related genes underlying quantitative disease resistance loci (QDRLs) for Phytophthora sojae. Thus, we hypothesized that auxin may be involved the P. sojae-soybean interaction. The levels of IAA and related metabolites were measured in mycelia and media supernatant as well as in mock and inoculated soybean roots in a time course assay. The expression of 11 soybean Pin-formed (GmPIN) auxin efflux transporter genes was also examined. Tryptophan, an auxin precursor, was detected in the P. sojae mycelia and media supernatant. During colonization of roots, levels of IAA and related metabolites were significantly higher in both moderately resistant Conrad and moderately susceptible Sloan inoculated roots compared with mock controls at 48 h postinoculation (hpi) in one experiment and at 72 hpi in a second, with Sloan accumulating higher levels of the auxin catabolite IAA-Ala than Conrad. Additionally, one GmPIN at 24 hpi, one at 48 hpi, and three at 72 hpi had higher expression in inoculated compared with the mock control roots in Conrad. The ability of resistant cultivars to cope with auxin accumulation may play an important role in quantitative disease resistance. Levels of jasmonic acid (JA), another plant hormone associated with defense responses, were also higher in inoculated roots at these same time points, suggesting that JA also plays a role during the later stages of infection.

Identification of Candidate Auxin Response Factors Involved in Pomegranate Seed Coat Development

Auxin response factors (ARFs) are transcription factors, regulating the auxin signaling pathways involved in plant development and related processes. In this study, we performed the genome-wide identification and characterization of ARFs in pomegranate and compared them with ARFs from three other species. Seventeen PgrARFs were identified and clustered into four groups, according to their phylogenetic relationship with the remaining 59 ARFs. A recent whole-genome duplication event in pomegranate may have contributed to the expansion and diversification of PgrARFs. Genomic truncation and variant splicing mechanisms contributed to the divergence of PgrARFs, a conclusion that was supported by different exon-intron structures of genes and incomplete conserved domains of PgrARFs in a specific phylogenetic group (group III). Interestingly, the absence of motifs from certain PgrARF genes corresponded to their low transcription levels, which contrasted to the highly expressed PgrARFs with intact motifs. Specifically, PgrARF1 and PgrARF2 highly expressed in both inner and outer seed coat, and phylogenetically related to Arabidopsis orthologs which mediates cell divisions in seed coat. We infer these two PgrARFs might involve in seed coat development through cell divisions in response to auxin regulation. These findings provided information on the characteristics and evolutionary relationships of PgrARFs, but also shed lights on their potential roles during seed coat development in pomegranate.

Auxin-mediated responses under salt stress: from developmental regulation to biotechnological applications

As sessile organisms, plants are exposed to multiple abiotic stresses commonly found in nature. To survive, plants have developed complex responses that involve genetic, epigenetic, cellular, and morphological modifications. Among different environmental cues, salt stress has emerged as a critical problem contributing to yield losses and marked reductions in crop production. Moreover, as the climate changes, it is expected that salt stress will have a significant impact on crop production in the agroindustry. On a mechanistic level, salt stress is known to be regulated by the crosstalk of many signaling molecules such as phytohormones, with auxin having been described as a key mediator of the process. Auxin plays an important role in plant developmental responses and stress, modulating a complex balance of biosynthesis, transport, and signaling that among other things, finely tune physiological changes in plant architecture and Na+ accumulation. In this review, we describe current knowledge on auxin's role in modulating the salt stress response. We also discuss recent and potential biotechnological approaches to tackling salt stress.

Genome-wide sequence identification and expression analysis of ARF family in sugar beet (Beta vulgaris L.) under salinity stresses

Auxin response factor (ARF) proteins respond to biological and abiotic stresses and play important roles in regulating plant growth and development. In this study, based on the genome-wide database of sugar beet, 16 BvARF proteins were identified. A detailed investigation into the BvARF family is performed, including analysis of the conserved domains, chromosomal locations, phylogeny, exon-intron structure, conserved motifs, subcellular localization, gene ontology (GO) annotations and expression profiles of BvARF under salt-tolerant condition. The majority of BvARF proteins contain B3 domain, AUX_RESP domain and AUX/IAA domain and a few lacked of AUX/IAA domain. Phylogenetic analysis suggests that the 16 BvARF proteins are clustered into six groups. Expression profile analysis shows that most of these BvARF genes in sugar beet under salinity stress were up-regulated or down-regulated to varying degrees and nine of the BvARF genes changed significantly. They were thought to have a significant response to salinity stress. The current study provides basic information for the BvARF genes and will pave the way for further studies on the roles of BvARF genes in regulating sugar beet's growth, development and responses to salinity stress.

Genome-wide characterization and expression analysis of the auxin response factor (ARF) gene family during melon (Cucumis melo L.) fruit development

ARFs in plants mediate auxin signaling transduction and regulate growth process. To determine genome-wide characterization of ARFs family in melon (Cucumis melo L.), ARFs were identified via analysis of information within the melon genomic database, and bioinformatic analyses were performed using various types of software. Based on different treatment methods involving dipping with the growth regulator Fengchanji No. 2 and artificial pollination, Jingmi No. 11 melon was used as the test material, and melon plants with unpollinated ovaries served as controls. The expression of ARFs during the early development of melon was analyzed via qRT-PCR. Seventeen genes that encode ARF proteins were identified in the melon genome for the first time. The expression of these ARFs differed in different tissues. The expression levels of CmARF2, CmARF16-like, CmARF18-like2, and CmARF19-like were especially high in melon fruits. The expression of ARFs during the early development of melon fruits differed in response to the different treatments, which suggested that CmARF9, CmARF16-like, CmARF19-like, CmARF19, CmARF1, CmARF2, CmARF3, and CmARF5 may be associated with melon fruit growth during early development. Interestingly, the increase in the transverse diameter of fruits treated with growth regulators was significantly greater than that of fruits resulting from artificial pollination, while the increase in the longitudinal diameter of the fruits resulting from artificial pollination was significantly greater.

Time-course RNA-seq analysis provides an improved understanding of gene regulation during the formation of nodule-like structures in rice

Using a time-course RNA-seq analysis we identified transcriptomic changes during formation of nodule-like structures (NLS) in rice and compared rice RNA-seq dataset with a nodule transcriptome dataset in Medicago truncatula. Plant hormones can induce the formation of nodule-like structures (NLS) in plant roots even in the absence of bacteria. These structures can be induced in roots of both legumes and non-legumes. Moreover, nitrogen-fixing bacteria can recognize and colonize these root structures. Therefore, identifying the genetic switches controlling the NLS organogenesis program in crops, especially cereals, can have important agricultural implications. Our recent study evaluated the transcriptomic response occurring in rice roots during NLS formation, 7 days post-treatment (dpt) with auxin, 2,4-D. In this current study, we investigated the regulation of gene expression occurring in rice roots at different stages of NLS formation: early (1-dpt) and late (14-dpt). At 1-dpt and 14-dpt, we identified 1662 and 1986 differentially expressed genes (DEGs), respectively. Gene ontology enrichment analysis revealed that the dataset was enriched with genes involved in auxin response and signaling; and in anatomical structure development and morphogenesis. Next, we compared the gene expression profiles across the three time points (1-, 7-, and 14-dpt) and identified genes that were uniquely or commonly differentially expressed at all three time points. We compared our rice RNA-seq dataset with a nodule transcriptome dataset in Medicago truncatula. This analysis revealed there is some amount of overlap between the molecular mechanisms governing nodulation and NLS formation. We also identified that some key nodulation genes were not expressed in rice roots during NLS formation. We validated the expression pattern of several genes via reverse transcriptase polymerase chain reaction (RT-PCR). The DEGs identified in this dataset may serve as a useful resource for future studies to characterize the genetic pathways controlling NLS formation in cereals.

Effects of Trichoderma asperellum and its siderophores on endogenous auxin in Arabidopsis thaliana under iron-deficiency stress

Iron (Fe) deficiency is one of the major limiting factors affecting crop yields. Trichoderma asperellum Q1, a biocontrol and plant growth promoting fungus, can produce the siderophore which has a high affinity to Fe³⁺ in the absence of iron. In this study, Trichoderma asperellum Q1 was found to be able to promote growth of Arabidopsis thaliana in an iron-deficient or insoluble iron-containing (Fe₂O₃) medium. It also can produce more siderophore and indole-3-acetic acid (IAA) as the concentration of iron ions decreased. However, it is unclear that the relationship between siderophore and IAA in promoting plant growth. Both Trichoderma asperellum Q1 and siderophore promotes not only the DR5::GFP transgenic Arabidopsis thaliana seedlings, in which the root IAA is labeled by green fluorescent protein gene, but also increases the content of endogenous IAA in the roots, which was shown by the fluorescence study. The strongest fluorescence was observed in the treated group inoculated with Trichoderma asperellum Q1 under the condition of insoluble iron. In the case of iron-free medium, adding siderophore also increased the observed fluorescence intensity. These results suggest that the siderophores produced by Trichoderma asperellum Q1 increased the content of IAA in Arabidopsis roots by enhancing the conversion of poorly soluble iron or by the siderophore itself.

Analysis of Aldo-Keto Reductase Gene Family and Their Responses to Salt, Drought, and Abscisic Acid Stresses in Medicago truncatula

Salt and drought stresses are two primary abiotic stresses that inhibit growth and reduce the activity of photosynthetic apparatus in plants. Abscisic acid (ABA) plays a key role in abiotic stress regulation in plants. Some aldo-keto reductases (AKRs) can enhance various abiotic stresses resistance by scavenging cytotoxic aldehydes in some plants. However, there are few comprehensive reports of plant AKR genes and their expression patterns in response to abiotic stresses. In this study, we identified 30 putative AKR genes from Medicago truncatula. The gene characteristics, coding protein motifs, and expression patterns of these MtAKRs were analyzed to explore and identify candidate genes in regulation of salt, drought, and ABA stresses. The phylogenetic analysis result indicated that the 52 AKRs in Medicago truncatula and Arabidopsis thaliana can be divided into three groups and six subgroups. Fifteen AKR genes in M. truncatula were randomly selected from each group or subgroup, to investigate their response to salt (200 mM of NaCl), drought (50 g·L^-1 of PEG 6000), and ABA (100 µM) stresses in both leaves and roots. The results suggest that MtAKR1, MtAKR5, MtAKR11, MtAKR14, MtAKR20, and MtAKR29 may play important roles in response to these stresses.

Comprehensive Atlas of Wheat (Triticum aestivum L.) AUXIN RESPONSE FACTOR Expression During Male Reproductive Development and Abiotic Stress

AUXIN RESPONSE FACTOR (ARF) proteins regulate a wide range of signaling pathways, from general plant growth to abiotic stress responses. Here, we performed a genome-wide survey in wheat (Triticum aestivum) and identified 69 TaARF members that formed 24 homoeologous groups. Phylogenetic analysis clustered TaARF genes into three clades, similar to ARF genes in Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa). Structural characterization suggested that ARF gene structure and domain composition are well conserved between plant species. Expression profiling revealed diverse patterns of TaARF transcript levels across a range of developmental stages, tissues, and abiotic stresses. A number of TaARF genes shared similar expression patterns and were preferentially expressed in anthers. Moreover, our systematic analysis identified three anther-specific TaARF genes (TaARF8, TaARF9, and TaARF21) whose expression was significantly altered by low temperature in thermosensitive genic male-sterile (TGMS) wheat; these TaARF genes are candidates to participate in the cold-induced male sterility pathway, and offer potential applications in TGMS wheat breeding and hybrid seed production. Moreover, we identified putative functions for a set of TaARFs involved in responses to abscisic acid and abiotic stress. Overall, this study characterized the wheat ARF gene family and generated several hypotheses for future investigation of ARF function during anther development and abiotic stress.

Auxin Response Factor 2A Is Part of the Regulatory Network Mediating Fruit Ripening Through Auxin-Ethylene Crosstalk in Durian

Fruit ripening is a highly coordinated developmental process driven by a complex hormonal network. Ethylene is the main regulator of climacteric fruit ripening. However, a putative role of other key phytohormones in this process cannot be excluded. We previously observed an increasing level of auxin during the post-harvest ripening of the durian fruit, which occurred concomitantly with the rise in the climacteric ethylene biosynthesis. Herein, we connect the key auxin signaling component, auxin response factors (ARFs), with the regulatory network that controls fruit ripening in durian through the identification and functional characterization of a candidate ripening-associated ARF. Our transcriptome-wide analysis identified 15 ARF members in durian (DzARFs), out of which 12 were expressed in the fruit pulp. Most of these DzARFs showed a differential expression, but DzARF2A had a marked ripening-associated expression pattern during post-harvest ripening in Monthong, a commercial durian cultivar from Thailand. Phylogenetic analysis of DzARF2A based on its tomato orthologue predicted a role in ripening through the regulation of ethylene biosynthesis. Transient expression of DzARF2A in Nicotiana benthamiana leaves significantly upregulated the expression levels of ethylene biosynthetic genes, pointing to a ripening-associated role of DzARF2A through the transcriptional regulation of ethylene biosynthesis. Dual-luciferase reporter assay determined that DzARF2A trans-activates durian ethylene biosynthetic genes. We previously reported significantly higher auxin level during post-harvest ripening in a fast-ripening cultivar (Chanee) compared to a slow-ripening one (Monthong). DzARF2A expression was significantly higher during post-harvest ripening in the fast-ripening cultivars (Chanee and Phuangmanee) compared to that of the slow-ripening ones (Monthong and Kanyao). Thus, higher auxin level could upregulate the expression of DzARF2A during ripening of a fast-ripening cultivar. The auxin-induced expression of DzARF2A confirmed its responsiveness to exogenous auxin treatment in a dose-dependent manner, suggesting an auxin-mediated role of DzARF2A in fruit ripening. We suggest that high DzARF2A expression would activate ARF2A-mediated transcription of ethylene biosynthetic genes, leading to increased climacteric ethylene biosynthesis (auxin-ethylene crosstalk) and faster ripening. Hence, we demonstrated DzARF2A as a new component of the regulatory network possibly mediating durian fruit ripening through transcriptional regulation of ethylene biosynthetic genes.

A comprehensive analysis of the B3 superfamily identifies tissue-specific and stress-responsive genes in chickpea (Cicer arietinum L.)

The aim of this study was to provide a comprehensive analysis of the plant-specific B3 domain-containing transcription factors (TFs) in chickpea. Scanning of the chickpea genome resulted in the identification of 51 B3 domain-containing TFs that were located on seven out of eight chickpea chromosomes. Based on the presence of additional domains other than the B3 domain, the candidates were classified into four subfamilies, i.e., ARF (24), REM (19), LAV (6) and RAV (2). Phylogenetic analysis classified them into four groups in which members of the same group had similar intron-exon organization and motif composition. Genome duplication analysis of the candidate B3 genes revealed an event of segmental duplication that was instrumental in the expansion of the B3 gene family. Ka/Ks analysis showed that the B3 gene family was under purifying selection. Further, chickpea B3 genes showed maximum orthology with Medicago followed by soybean and Arabidopsis. Promoter analyses of the B3 genes led to the identification of several tissue-specific and stress-responsive cis-regulatory elements. Expression profiling of the candidate B3 genes using publicly available RNA-seq data of several chickpea tissues indicated their putative role in plant development and abiotic stress response. These findings were further validated by real-time expression analysis. Overall, this study provides a comprehensive analysis of the B3 domain-containing proteins in chickpea that would aid in devising strategies for crop manipulation in chickpea.

geneHummus: an R package to define gene families and their expression in legumes and beyond

BACKGROUND

During the last decade, plant biotechnological laboratories have sparked a monumental revolution with the rapid development of next sequencing technologies at affordable prices. Soon, these sequencing technologies and assembling of whole genomes will extend beyond the plant computational biologists and become commonplace within the plant biology disciplines. The current availability of large-scale genomic resources for non-traditional plant model systems (the so-called 'orphan crops') is enabling the construction of high-density integrated physical and genetic linkage maps with potential applications in plant breeding. The newly available fully sequenced plant genomes represent an incredible opportunity for comparative analyses that may reveal new aspects of genome biology and evolution. The analysis of the expansion and evolution of gene families across species is a common approach to infer biological functions. To date, the extent and role of gene families in plants has only been partially addressed and many gene families remain to be investigated. Manual identification of gene families is highly time-consuming and laborious, requiring an iterative process of manual and computational analysis to identify members of a given family, typically combining numerous BLAST searches and manually cleaning data. Due to the increasing abundance of genome sequences and the agronomical interest in plant gene families, the field needs a clear, automated annotation tool.

RESULTS

Here, we present the geneHummus package, an R-based pipeline for the identification and characterization of plant gene families. The impact of this pipeline comes from a reduction in hands-on annotation time combined with high specificity and sensitivity in extracting only proteins from the RefSeq database and providing the conserved domain architectures based on SPARCLE. As a case study we focused on the auxin receptor factors gene (ARF) family in Cicer arietinum (chickpea) and other legumes.

CONCLUSION

We anticipate that our pipeline should be suitable for any taxonomic plant family, and likely other gene families, vastly improving the speed and ease of genomic data processing.

RNA-seq reveals differentially expressed genes in rice (Oryza sativa) roots during interactions with plant-growth promoting bacteria, Azospirillum brasilense

Major non-legume crops can form beneficial associations with nitrogen-fixing bacteria like Azospirillum brasilense. Our current understanding of the molecular aspects and signaling that occur between important crops like rice and these nitrogen-fixing bacteria is limited. In this study, we used an experimental system where the bacteria could colonize the plant roots and promote plant growth in wild type rice and symbiotic mutants (dmi3 and pollux) in rice. Our data suggest that plant growth promotion and root penetration is not dependent on these genes. We then used this colonization model to identify regulation of gene expression at two different time points during this interaction: at 1day post inoculation (dpi), we identified 1622 differentially expressed genes (DEGs) in rice roots, and at 14dpi, we identified 1995 DEGs. We performed a comprehensive data mining to classify the DEGs into the categories of transcription factors (TFs), protein kinases (PKs), and transporters (TRs). Several of these DEGs encode proteins that are involved in the flavonoid biosynthetic pathway, defense, and hormone signaling pathways. We identified genes that are involved in nitrate and sugar transport and are also implicated to play a role in other plant-microbe interactions. Overall, findings from this study will serve as an excellent resource to characterize the host genetic pathway controlling the interactions between non-legumes and beneficial bacteria which can have long-term implications towards sustainably improving agriculture.

New insight into the molecular basis of cadmium stress responses of wild paper mulberry plant by transcriptome analysis

BACKGROUND

Heavy metal contamination is becoming a limitation to the utilization of soil and the distribution of vegetation. In particular, cadmium (Cd) pollution has had a serious impact on the food chain. Broussonetia papyrifera is a widely distributed pioneer tree species of heavy metal contaminated areas with important economic value. However, little is known about the genomic background of the Cd-tolerance mechanism in B. papyrifera.

RESULTS

The CdCl₂ responsive physiology was evaluated and proved to be involved in antioxidase activity and active oxygen species (ROS) accumulation. The leaf and root transcriptomes derived from B. papyrifera grown under normal and CdCl₂ stress conditions were systematically investigated using the Illumina HiSeq method. A total of 180,678,660 bp (27.1 GB) clean reads were assembled into 589,487 high-quality unigenes, of which 256,025 (43.43% of the total) and 250,251 (42.45% of the total) were aligned in Gene Ontology (GO) and Protein family (Pfam), respectively. A total of 24,414 differentially expressed genes (DEGs) were GO-annotated into 53, 23, 55, and 60 terms from the transcriptomes of root and leaf tissues under Cd stress and control conditions. A total of 117,547 Kyoto Encyclopedia of Genes and Genomes (KEGG) Orthology (KO)-annotated DEGs were enriched in at least 47 KEGG pathway terms among the four comparisons. Many genes encoding important transcription factors (e.g., auxin/indole-3-acetic acid (AUX/IAA), basic helix-loop-helix (bHLH), DNA-binding one zinc finger (Dof), and MYB) and proteins involved in plant-pathogen interactions, phenylpropanoid biosynthesis, plant hormone signal transduction, oxidative phosphorylation, carbon fixation, peroxisomes, flavonoid biosynthesis, and glutathione metabolism, among others, were substantially upregulated under CdCl₂ stress.

CONCLUSIONS

These genes represent important candidates for studying Cd-response mechanisms and molecular biology of B. papyrifera and related species. Our findings provide a genomic sequence resource for functional genetic assignments in B. papyrifera, which will help elucidate the molecular mechanisms of its Cd-stress responses and facilitate the bioremediation of heavy metal contaminated areas via breeding of new stress-tolerant cultivars.

Genome-wide identification of the GhARF gene family reveals that GhARF2 and GhARF18 are involved in cotton fibre cell initiation

Auxin signalling plays an essential role in regulating plant development. Auxin response factors (ARFs), which are critical components of auxin signalling, modulate the expression of early auxin-responsive genes by binding to auxin response factor elements (AuxREs). However, there has been no comprehensive characterization of this gene family in cotton. Here, we identified 56 GhARF genes in the assembled Gossypium hirsutum genome. This gene family was divided into 17 subfamilies, and 44 members of them were distributed across 21 chromosomes. GhARF6 and GhARF11 subfamily genes were predominantly expressed in vegetative tissues, whereas GhARF2 and GhARF18 subfamily genes were highly expressed during seed fibre cell initiation. GhARF2-1 and GhARF18-1 were exclusively expressed in trichomes, organs similar to cotton seed fibre cells, and overexpression of these genes in Arabidopsis enhances trichome initiation. Comparative transcriptome analysis combined with AuxRE prediction revealed 11 transcription factors as potential target genes of GhARF2 and GhARF18. Six of these genes were significantly expressed during seed fibre cell initiation and were bound by GhARF2-1 and GhARF18-1 in yeast one-hybrid assays. Our results suggest that GhARF2 and GhARF18 genes may be key regulators of cotton seed fibre initiation by regulating the expression of several transcription factor genes. This study deepens our understanding of auxin-mediated initiation of cotton seed fibre cells and helps us in breeding better cotton varieties in the future.

Hormone modulation of legume-rhizobial symbiosis

Leguminous plants can establish symbiotic associations with diazotropic rhizobia to form nitrogen-fixating nodules, which are classified as determinate or indeterminate based on the persistence of nodule meristem. The formation of nitrogen-fixing nodules requires coordinating rhizobial infection and root nodule organogenesis. The formation of an infection thread and the extent of nodule formation are largely under plant control, but vary with environmental conditions and the physiological state of the host plants. Many achievements in these two areas have been made in recent decades. Phytohormone signaling pathways have gradually emerged as important regulators of root nodule symbiosis. Cytokinin, strigolactones (SLs) and local accumulation of auxin can promote nodule development. Ethylene, jasmonic acid (JA), abscisic acid (ABA) and gibberellic acid (GA) all negatively regulate infection thread formation and nodule development. However, salicylic acid (SA) and brassinosteroids (BRs) have different effects on the formation of these two nodule types. Some peptide hormones are also involved in nodulation. This review summarizes recent findings on the roles of these plant hormones in legume-rhizobial symbiosis, and we propose that DELLA proteins may function as a node to integrate plant hormones to regulate nodulation.

Identification, characterization, and expression analysis of auxin response factor (ARF) gene family in Brachypodium distachyon

Auxin response factors (ARFs) are one type of essential family of transcription factors that bind with auxin response elements (AuxRE), and play vital roles in variety of plant development and physiological processes. Brachypodium distachyon, related to the major cereal grain species, were recently developed to be a good model organism for functional genomics research. So far, genome-wide overview of the ARF gene family in B. distachyon was not available. Here, a systemic analysis of ARF gene family members in B. distachyon was performed. A comprehensive overview of the characterization of the BdARFs was obtained by multiple bioinformatics analyses, including the gene and protein structure, chromosome locations, conserved motifs of proteins, phylogenetic analysis, and cis-elements in promoters of BdARF. Results showed that all BdARFs contained conserved DBD, MR, and CTD could be divided into four classes, Ia, IIa, IIb, and III. Expression profiles of BdARF genes indicated that they were expressed across various tissues and organs, which could be clustered into three main expression groups, and most of BdARF genes were involved in phytohormone signal transduction pathways and regulated physiological process in responding to multiple environmental stresses. And predicted regulatory network between B. distachyon ARFs and IAAs was also discussed. Our genomics analysis of BdARFs could yield new insights into the complexity of the control of BdARF genes and lead to potential applications in the investigation of the accurate regulatory mechanisms of ARFs in herbaceous plants.

Genome-Wide Identification, Phylogeny, and Expression Analysis of ARF Genes Involved in Vegetative Organs Development in Switchgrass

Auxin response factors (ARFs) have been reported to play vital roles during plant growth and development. In order to reveal specific functions related to vegetative organs in grasses, an in-depth study of the ARF gene family was carried out in switchgrass (Panicum virgatum L.), a warm-season C4 perennial grass that is mostly used as bioenergy and animal feedstock. A total of 47 putative ARF genes (PvARFs) were identified in the switchgrass genome (2n = 4x = 36), 42 of which were anchored to the seven pairs of chromosomes and found to be unevenly distributed. Sixteen PvARFs were predicted to be potential targets of small RNAs (microRNA160 and 167). Phylogenetically speaking, PvARFs were divided into seven distinct subgroups based on the phylogeny, exon/intron arrangement, and conserved motif distribution. Moreover, 15 pairs of PvARFs have different temporal-spatial expression profiles in vegetative organs (2nd, 3rd, and 4th internode and leaves), which implies that different PvARFs have specific functions in switchgrass growth and development. In addition, at least 14 pairs of PvARFs respond to naphthylacetic acid (NAA) treatment, which might be helpful for us to study on auxin response in switchgrass. The comprehensive analysis, described here, will facilitate the future functional analysis of ARF genes in grasses.

Identification and expression profiling of the Aux/IAA gene family in Chinese hickory (Carya cathayensis Sarg.) during the grafting process

Auxin is an essential regulator in various aspects of organism growth and development. Members of the Aux/IAA family of genes encode short-lived nuclear proteins and mediate the responses of auxin-regulated gene expression. Here, the first identification and characterization of 22 cDNAs encoding the open reading frame of the Aux/IAA family in Chinese hickory (named as CcIAA) has been performed. The proteins encoded by these genes contain four whole or partially conserved domains of the Aux/IAA family. Phylogenetic analysis indicated that CcIAAs were unevenly distributed among eight different subgroups. The spatio-specific expression profiles showed that most of the CcIAAs preferentially expressed in specific tissues. Three CcIAA genes, including CcIAA11, CcIAA27a2 and CcIAAx, were predominantly expressed in stem. The predominant expression of CcIAA genes in stems might play important roles in vascular reconnection during the graft process. Furthermore, expression profiles of Aux/IAA genes during the grafting process of Chinese hickory have been analysed. Our data suggested that 19 CcIAAs were down-regulated and 3 CcIAAs (including CcIAA28, CcIAA8a and CcIAA27b) were induced, indicating their specializations during the grafting process. The involvement of CcIAA genes at the early stage after grafting gives us an opportunity to understand the role of auxin signalling in the grafting process.

Genome-wide identification of the auxin response factor gene family in Cicer arietinum

Background

Auxin Response Factors act as critical components of the auxin-signaling pathway by regulating the transcription of auxin-responsive genes. The release of the chickpea reference genome provides an opportunity to identify and characterize the ARF gene family in this important legume by a data mining coupled by comparative genomics approaches.

Results

We performed a comprehensive characterization and analysis of 24 ARF genes in the chickpea reference genome. Comparative phylogenetic analysis of the ARF from chickpea, Medicago and Arabidopsis suggests that recent duplications have played a very limited role in the expansion of the ARF chickpea family. Gene structure analysis based on exon-intron organization provides additional evidence to support the evolutionary relationship among the ARF members. Conserved motif analysis shows that most of the proteins fit into the canonical ARF structure model, but 9 proteins lack or have a truncated dimerization domain. The mechanisms underlying the diversification of the ARF gene family are based on duplications, variations in domain organization and alternative splicing. Concerning duplications, segmental, but not tandem duplications, have contributed to the expansion of the gene family. Moreover, the duplicated pair genes have evolved mainly under the influence of purifying selection pressure with restricted functional divergence. Expression profiles responding to various environmental stimuli show a close relationship between tissue and expression patterns. Promoter sequence analysis reveals an enrichment of several cis-regulatory elements related to symbiosis, and modulation of plant gene expression during the interaction with microbes.

Conclusions

In conclusion, this study provides a comprehensive overview of the ARF gene family in chickpea. Globally, our data supports that auxin signaling pathway regulates a wide range of physiological processes and stress responses. Our findings could further provide new insights into the complexity of the regulation of ARF at the transcription level that may be useful to develop rational chickpea breeding strategies to improve development or stress responses. Our study also provides a foundation for comparative genomic analyses and a framework to trace the dynamic evolution of ARF genes on a large time-scale within the legume family.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4695-9) contains supplementary material, which is available to authorized users.

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.