Genome-wide analysis of auxin response factor gene family members in medicinal model plant Salvia miltiorrhiza

A non-canonical Aux/IAA gene MsIAA32 regulates peltate glandular trichome development in spearmint

Phytohormone auxin controls various aspects of plant growth and development. The typical auxin signalling involves the degradation of canonical Aux/IAA proteins upon auxin perception releasing the auxin response factors (ARF) to activate auxin-regulated gene expression. Extensive research has been pursued in deciphering the role of canonical Aux/IAAs, however, the function of non-canonical Aux/IAA genes remains elusive. Here we identified a non-canonical Aux/IAA gene, MsIAA32 from spearmint (Mentha spicata), which lacks the TIR1-binding domain and shows its involvement in the development of peltate glandular trichomes (PGT), which are the sites for production and storage of commercially important essential oils. Using yeast two-hybrid studies, two canonical Aux/IAAs, MsIAA3, MsIAA4 and an ARF, MsARF3 were identified as the preferred binding partners of MsIAA32. Expression of a R2R3-MYB gene MsMYB36 and a cyclin gene MsCycB2-4 was altered in MsIAA32 suppressed plants indicating that these genes are possible downstream targets of MsIAA32 mediated signalling. Ectopic expression of MsIAA32 in Arabidopsis affected non-glandular trichome formation along with other auxin related developmental traits. Our findings establish the role of non-canonical Aux/IAA mediated auxin signalling in PGT development and reveal species-specific functionalization of Aux/IAAs.

Comparative Microscopic, Transcriptome and IAA Content Analyses Reveal the Stem Growth Variations in Two Cultivars Ilex verticillata

Ilex verticillata is not only an excellent ornamental tree species for courtyards, but it is also a popular bonsai tree. 'Oosterwijk' and 'Red sprite' are two varieties of Ilex verticillata. The former has a long stem with few branches, while the latter has a short stem. In order to explain the stem growth differences between the two cultivars 'Oosterwijk' and 'Red sprite', determination of the microstructure, transcriptome sequence and IAA content was carried out. The results showed that the xylem thickness, vessel area and vessel number of 'Oosterwijk' were larger than in 'Red sprite'. In addition, our analysis revealed that the differentially expressed genes which were enriched in phenylpropanoid biosynthesis; phenylalanine metabolism and phenylalanine, tyrosine and tryptophan biosynthesis in the black and tan modules of the two varieties. We found that AST, HCT and bHLH 94 may be key genes in the formation of shoot difference. Moreover, we found that the IAA content and auxin-related DEGs GH3.6, GH3, ATRP5, IAA27, SAUR36-like, GH3.6-like and AIP 10A5-like may play important roles in the formation of shoot differences. In summary, these results indicated that stem growth variations of 'Oosterwijk' and 'Red sprite' were associated with DEGs related to phenylpropanoid biosynthesis, phenylalanine metabolism and phenylalanine, tyrosine and tryptophan biosynthesis, as well as auxin content and DEGs related to the auxin signaling pathway.

LkARF7 and LkARF19 overexpression promote adventitious root formation in a heterologous poplar model by positively regulating LkBBM1

Cuttage propagation involves adventitious root formation induced by auxin. In our previous study, Larix kaempferi BABY BOOM 1 (LkBBM1), which is known to regulate adventitious root formation, was affected by auxin. However, the relationship between LkBBM1 and auxin remains unclear. Auxin response factors (ARFs) are a class of important transcription factors in the auxin signaling pathway and modulate the expression of early auxin-responsive genes by binding to auxin response elements. In the present study, we identified 14 L. kaempferi ARFs (LkARFs), and found LkARF7 and LkARF19 bound to LkBBM1 promoter and enhanced its transcription using yeast one-hybrid, ChIP-qPCR, and dual-luciferase assays. In addition, the treatment with naphthalene acetic acid promoted the expression of LkARF7 and LkARF19. We also found that overexpression of these two genes in poplar promoted adventitious root formation. Furthermore, LkARF19 interacted with the DEAD-box ATP-dependent RNA helicase 53-like protein to form a heterodimer to regulate adventitious root formation. Altogether, our results reveal an additional regulatory mechanism underlying the control of adventitious root formation by auxin.

Tubby-like proteins (TLPs) transcription factor in different regulatory mechanism in plants: a review

Tubby-like proteins (TLPs) transcription factors are found in single-celled to multi-cellular eukaryotes in the form of large multigene families. TLPs are identified through a specific signature of carboxyl terminal tubby domain, required for plasma membrane tethering and amino terminal F-box domain communicate as functional SCF-type E3 ligases. The comprehensive distribution of TLP gene family members in diverse species indicates some conserved functions of TLPs in multicellular organisms. Plant TLPs have higher gene members than animals and these members reported important role in multiple physiological and developmental processes and various environmental stress responses. Although the TLPs are suggested to be a putative transcription factors but their functional mechanism is not much clear. This review provides significant recent updates on TLP-mediated regulation with an insight into its functional roles, origin and evolution and also phytohormones related regulation to combat with various stresses and its involvement in adaptive stress response in crop plants.

The Expression Profiles of the Salvia miltiorrhiza 3-Hydroxy-3-methylglutaryl-coenzyme A Reductase 4 Gene and Its Influence on the Biosynthesis of Tanshinones

Salvia miltiorrhiza is a medicinal plant that synthesises biologically-active tanshinones with numerous therapeutic properties. An important rate-limiting enzyme in the biosynthesis of their precursors is 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMGR). This study presents the organ-specific expression profile of the S. miltiorrhiza HMGR4 gene and its sensitivity to potential regulators, viz. gibberellic acid (GA₃), indole-3-acetic acid (IAA) and salicylic acid (SA). In addition, it demonstrates the importance of the HMGR4 gene, the hormone used, the plant organ, and the culture environment for the biosynthesis of tanshinones. HMGR4 overexpression was found to significantly boost the accumulation of dihydrotanshinone I (DHTI), cryptotanshinone (CT), tanshinone I (TI) and tanshinone IIA (TIIA) in roots by 0.44 to 5.39 mg/g dry weight (DW), as well as TIIA in stems and leaves. S. miltiorrhiza roots cultivated in soil demonstrated higher concentrations of the examined metabolites than those grown in vitro. GA₃ caused a considerable increase in the quantity of CT (by 794.2 µg/g DW) and TIIA (by 88.1 µg/g DW) in roots. In turn, IAA significantly inhibited the biosynthesis of the studied tanshinones in root material.

Genome-wide identification of the Tubby-Like Protein (TLPs) family in medicinal model plant Salvia miltiorrhiza

Tubby-Like Proteins (TLPs) are important transcription factors with many functions and are found in both animals and plants. In plants, TLPs are thought to be involved in the abiotic stress response. To reveal the potential function of TLPs in the medicinal model plant Salvia miltiorrhiza, we identified 12 S. miltiorrhiza TLPs (SmTLPs) and conducted a comprehensive analysis. We examined SmTLP gene structure, protein structure, phylogenetics, and expression analysis. Our results show that all SmTLPs, except SmTLP11, have a complete typical Tub domain. Promoter analysis revealed that most SmTLPs are involved in hormone and abiotic stress responses. Expression analysis revealed that the 12 SmTLPs could be divided into three categories: those specifically expressed in roots, those specifically expressed in stems, and those specifically expressed in leaves. Additional studies have shown that SmTLP10 may play an important role in the plant cold resistance, while SmTLP12 may be involved in the S. miltiorrhiza ABA metabolic pathway. Our study represents the first comprehensive investigation of TLPs in S. miltiorrhiza. These data may provide useful clues for future studies and may support the hypotheses regarding the role of TLPs in plant abiotic stress process. All in all, we may provide a reference for improving S. miltiorrhiza quality using genetic engineering technology.

Genome-Wide Characterization and Analysis of bHLH Transcription Factors Related to Crocin Biosynthesis in Gardenia jasminoides Ellis (Rubiaceae)

Crocins, enriched in Gardenia jasminoides fruits, have a pharmacological activity against central nervous system diseases, cardiovascular diseases, and cancer cell growth. The biosynthesis of crocins has been widely explored, but its regulatory mechanism remains unknown. Here, the basic helix-loop-helix (bHLH) transcription factors related to crocin biosynthesis were systematically identified on the basis of the genome of G. jasminoides. A total of 95 GjbHLH transcription factor genes were identified, and their phylogenetic analysis indicated that they could be classified into 23 subfamilies. The combination of gene-specific bHLH expression patterns, the coexpression analysis of biosynthesis genes, and the analysis of promoter sequences in crocin biosynthesis pathways suggested that nine bHLHs in G. jasminoides might negatively regulate crocin biosynthesis. This study laid a foundation for understanding the regulatory mechanism of crocin biosynthesis and the improvement and breeding of G. jasminoides varieties.

ARF family identification in Tamarix chinensis reveals the salt responsive expression of TcARF6 targeted by miR167

Auxin response factors (ARFs) are important transcription factors (TFs) that are differentially expressed in response to various abiotic stresses. The important roles of ARFs and small RNA-ARF pathways in mediating plant growth and stress responses have emerged in several recent studies. However, no studies on the involvement of ARFs in tamarisk trees, which are resistant to salinity, have been conducted. In this study, systematic analysis revealed 12 TcARF genes belonging to five different groups in Tamarix chinensis. The microRNA response elements of miR160, which belongs to group I and miR167, which belongs to group III, were conserved in terms of their location and sequence. Moreover, digital gene expression profiles suggested that a potential miR167 target gene, TcARF6, was rapidly expressed in response to salt stress. Cloning of TcARF6 revealed that TcARF6 could be an activation TF with a glutamine-rich region and expression pattern analysis revealed that the expression of TcARF6 was significantly downregulated specifically in the roots. A significant negative correlation in the expression pattern of tch-miR167/TcARF6 indicated that this module may play a key role in the response to salt stress. Overall, these results provide basic information on the posttranscriptional regulation of TcARF6 for future investigations of the T. chinensis salt-stress response.

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.

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

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

The auxin response factor gene family in allopolyploid Brassica napus

Auxin response factor (ARF) is a member of the plant-specific B3 DNA binding superfamily. Here, we report the results of a comprehensive analysis of ARF genes in allotetraploid Brassica napus (2n = 38, AACC). Sixty-seven ARF genes were identified in B. napus (BnARFs) and divided into four subfamilies (I–IV). Sixty-one BnARFs were distributed on all chromosomes except C02; the remaining were on Ann and Cnn. The full length of the BnARF proteins was highly conserved especially within each subfamily with all members sharing the N-terminal DNA binding domain (DBD) and the middle region (MR), and most contained the C-terminal dimerization domain (PBI). Twenty-one members had a glutamine-rich MR that may be an activator and the remaining were repressors. Accordingly, the intron patterns are highly conserved in each subfamily or clade, especially in DBD and PBI domains. Several members in subfamily III are potential targets for miR167. Many putative cis-elements involved in phytohormones, light signaling responses, and biotic and abiotic stress were identified in BnARF promoters, implying their possible roles. Most ARF proteins are likely to interact with auxin/indole-3-acetic acid (Aux/IAA) -related proteins, and members from different subfamilies generally shared many common interaction proteins. Whole genome-wide duplication (WGD) by hybridization between Brassica rapa and Brassica oleracea and segmental duplication led to gene expansion. Gene loss following WGD is biased with the A_n-subgenome retaining more ancestral genes than the C_n-subgenome. BnARFs have wide expression profiles across vegetative and reproductive organs during different developmental stages. No obvious expression bias was observed between A_n- and C_n-subgenomes. Most synteny-pair genes had similar expression patterns, indicating their functional redundancy. BnARFs were sensitive to exogenous IAA and 6-BA treatments especially subfamily III. The present study provides insights into the distribution, phylogeny, and evolution of ARF gene family.

Differential Expression of MicroRNAs During Root Formation in Taxus Chinensis Var. mairei Cultivars

MicroRNAs (miRNAs) have been shown to play key roles in the regulation of plant growth and development by modifying the expression of their target genes. However, the influence of miRNAs on root formation and development in woody plants, such as Taxus chinensis, remains largely unknown. In the current study, we explored the phytohormone-response and nutrition-response miRNA expression profiles during T. chinensis rooting by quantitative real-time PCR (qPCR). We identified six phytohormone-response miRNAs, namely, miR164a, miR165, miR167a, miR171b, miR319, and miR391, and eight nutrition-response miRNAs, namely, miR169b, miR395a, miR399c, miR408, miR826, miR827, miR857, and miR2111a, that were differentially expressed at different rooting phases of T. chinensis. Using northern blot analysis of the putative target genes of these miRNAs, we detected the relative gene expression changes of the target genes. Taken together, our results suggest that miRNAs are involved in root formation of T. chinensis and that miRNAs may play important regulatory roles in primary root, crown root, and root hair formation by targeting phytohormone and/or nutrition response genes in T. chinensis. For the first time, these results expand our understanding of the molecular mechanisms of plant root formation and development in a conifer species.

Trends in herbgenomics

From Shen Nong's Herbal Classic (Shennong Bencao Jing) to the Compendium of Materia Medica (Bencao Gangmu) and the first scientific Nobel Prize for the mainland of China, each milestone in the historical process of the development of traditional Chinese medicine (TCM) involves screening, testing and integrating. After thousands of years of inheritance and development, herbgenomics (bencaogenomics) has bridged the gap between TCM and international advanced omics studies, promoting the application of frontier technologies in TCM. It is a discipline that uncovers the genetic information and regulatory networks of herbs to clarify their molecular mechanism in the prevention and treatment of human diseases. The main theoretical system includes genomics, functional genomics, proteomics, transcriptomics, metabolomics, epigenomics, metagenomics, synthetic biology, pharmacogenomics of TCM, and bioinformatics, among other fields. Herbgenomics is mainly applicable to the study of medicinal model plants, genomic-assisted breeding, herbal synthetic biology, protection and utilization of gene resources, TCM quality evaluation and control, and TCM drug development. Such studies will accelerate the application of cutting-edge technologies, revitalize herbal research, and strongly promote the development and modernization of TCM.

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.

Genome-wide characterization and analysis of bHLH transcription factors in Panax ginseng

Ginseng (Panax ginseng C.A. Meyer) is one of the best-selling herbal medicines, with ginsenosides as its main pharmacologically active constituents. Although extensive chemical and pharmaceutical studies of these compounds have been performed, genome-wide studies of the basic helix-loop-helix (bHLH) transcription factors of ginseng are still limited. The bHLH transcription factor family is one of the largest transcription factor families found in eukaryotic organisms, and these proteins are involved in a myriad of regulatory processes. In our study, 169 bHLH transcription factor genes were identified in the genome of P. ginseng, and phylogenetic analysis indicated that these PGbHLHs could be classified into 24 subfamilies. A total of 21 RNA-seq data sets, including two sequencing libraries for jasmonate (JA)-responsive and 19 reported libraries for organ-specific expression analyses were constructed. Through a combination of gene-specific expression patterns and chemical contents, 6 PGbHLH genes from 4 subfamilies were revealed to be potentially involved in the regulation of ginsenoside biosynthesis. These 6 PGbHLHs, which had distinct target genes, were further divided into two groups depending on the absence of MYC-N structure. Our results would provide a foundation for understanding the molecular basis and regulatory mechanisms of bHLH transcription factor action in P. ginseng.

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.

Genomic survey of bZIP transcription factor genes related to tanshinone biosynthesis in Salvia miltiorrhiza

Tanshinones are a class of bioactive components in the traditional Chinese medicine Salvia miltiorrhiza, and their biosynthesis and regulation have been widely studied. Current studies show that basic leucine zipper (bZIP) proteins regulate plant secondary metabolism, growth and developmental processes. However, the bZIP transcription factors involved in tanshinone biosynthesis are unknown. Here, we conducted the first genome-wide survey of the bZIP gene family and analyzed the phylogeny, gene structure, additional conserved motifs and alternative splicing events in S. miltiorrhiza. A total of 70 SmbZIP transcription factors were identified and categorized into 11 subgroups based on their phylogenetic relationships with those in Arabidopsis. Moreover, seventeen SmbZIP genes underwent alternative splicing events. According to the transcriptomic data, the SmbZIP genes that were highly expressed in the Danshen root and periderm were selected. Based on the prediction of bZIP binding sites in the promoters and the co-expression analysis and co-induction patterns in response to Ag⁺ treatment via quantitative real-time polymerase chain reaction (qRT-PCR), we concluded that SmbZIP7 and SmbZIP20 potentially participate in the regulation of tanshinone biosynthesis. These results provide a foundation for further functional characterization of the candidate SmbZIP genes, which have the potential to increase tanshinone production.

Identification and Expression Profiling of the Auxin Response Factors in Capsicum annuum L. under Abiotic Stress and Hormone Treatments

Auxin response factors (ARFs) play important roles in regulating plant growth and development and response to environmental stress. An exhaustive analysis of the CaARF family was performed using the latest publicly available genome for pepper (Capsicum annuum L.). In total, 22 non-redundant CaARF gene family members in six classes were analyzed, including chromosome locations, gene structures, conserved motifs of proteins, phylogenetic relationships and Subcellular localization. Phylogenetic analysis of the ARFs from pepper (Capsicum annuum L.), tomato (Solanum lycopersicum L.), Arabidopsis and rice (Oryza sativa L.) revealed both similarity and divergence between the four ARF families, and aided in predicting biological functions of the CaARFs. Furthermore, expression profiling of CaARFs was obtained in various organs and tissues using quantitative real-time RT-PCR (qRT-PCR). Expression analysis of these genes was also conducted with various hormones and abiotic treatments using qRT-PCR. Most CaARF genes were regulated by exogenous hormone treatments at the transcriptional level, and many CaARF genes were altered by abiotic stress. Systematic analysis of CaARF genes is imperative to elucidate the roles of CaARF family members in mediating auxin signaling in the adaptation of pepper to a challenging environment.

Functional Characterization of Selected Universal Stress Protein from Salvia miltiorrhiza (SmUSP) in Escherichia coli

The multigene universal stress protein (USP) family is evolutionarily conserved. Members play indispensable roles in plant tolerance to abiotic stresses. Although relatively well-characterized in model plants, such as Arabidopsis thaliana and Oryzasativa, this family has not been investigated in Salvia miltiorrhiza, an important herbal plant for which yields can be limited by various abiotic stresses. Here, we identified 32 USP family members in the S. miltiorrhiza genome, and used phylogenetic analysis to sort these SmUSPs into four groups. Groups A and B belong to the ATP-binding class whereas Groups C and D are in the non-ATP-binding class. Motif analysis and multiple sequence alignment hinted that members of group A and B were able to bind ATP. Our qRT-PCR data from different tissues/organs and under salt and heat stresses provided an overall expression pattern for those genes. Three SmUSPs (SmUSP1, SmUSP8, and SmUSP27) were cloned from S. miltiorrhiza and functionally characterized in Escherichiacoli. Compared with the control cells, those that expressed SmUSPs exhibited enhanced tolerance to salt, heat, and a combination of the two. This suggested that the protein has a protective role in cells when exposed to single-stress and multiple-stress conditions. Our findings provide valuable information that helps improve our understanding of the evolutionary and functional conservation and diversity associated with the USP gene family in S. miltiorrhiza.