Transcriptome profiling of Gossypium arboreum during fiber initiation and the genome-wide identification of trihelix transcription factors

Functional divergence of two arginine decarboxylase genes in tropane alkaloid biosynthesis and root growth in Atropa belladonna

Putrescine, produced via the arginine decarboxylase (ADC)/ornithine decarboxylase (ODC)-mediated pathway, is an initial precursor for polyamines metabolism and the root-specific biosynthesis of medicinal tropane alkaloids (TAs). These alkaloids are widely used as muscarinic acetylcholine antagonists in clinics. Although the functions of ODC in biosynthesis of polyamines and TAs have been well investigated, the role of ADC is still poorly understood. In this study, enzyme inhibitor treatment showed that ADC was involved in the biosynthesis of putrescine-derived metabolites and root growth in Atropa belladonna. Further analysis found that there were six ADC unigenes in the A. belladonna transcriptome, with two of them, AbADC1 and AbADC2, exhibiting high expression in the roots. To investigate their roles in TAs/polyamines metabolism and root growth, RNA interference (RNAi) was used to suppress either AbADC1 or AbADC2 expression in A. belladonna hairy roots. Suppression of the AbADC1 expression resulted in a significant reduction in the putrescine content and hairy root biomass. However, it had no noticeable effect on the levels of N-methylputrescine and the TAs hyoscyamine, anisodamine, and scopolamine. On the other hand, suppression of AbADC2 expression markedly reduced the levels of putrescine, N-methylputrescine, and TAs, but had no significant effect on hairy root biomass. According to β-glucuronidase (GUS) staining assays, AbADC1 was mainly expressed in the root elongation and division region while AbADC2 was mainly expressed in the cylinder of the root maturation region. These differences in expression led to functional divergence, with AbADC1 primarily regulating root growth and AbADC2 contributing to TA biosynthesis.

Large-scale analysis of trihelix transcription factors reveals their expansion and evolutionary footprint in plants

The trihelix transcription factor (TTF) gene family is an important class of transcription factors that play key roles in regulating developmental processes and responding to various stresses. To date, no comprehensive analysis of the TTF gene family in large-scale species has been performed. A cross-genome exploration of its origin, copy number variation, and expression pattern in plants is also unavailable. Here, we identified and characterized the TTF gene family in 110 species representing typical plant phylogenetic taxa. Interestingly, we found that the number of TTF genes was significantly expanded in Chara braunii compared to other species. Based on the available plant genomic datasets, our comparative analysis suggested that the TTF gene family likely originated from the GT-1-1 group and then expanded to form other groups through duplication or deletion of some domains. We found evidence that whole-genome duplication/triplication contributed most to the expansion of the TTF gene family in dicots, monocots and basal angiosperms. In contrast, dispersed and proximal duplications contributed to the expansion of the TTF gene family in algae and bryophyta. The expression patterns of TTF genes and their upstream and downstream genes in different treatments showed a functional divergence of TTF-related genes. Furthermore, we constructed the interaction network between TTF genes and the corresponding upstream and downstream genes, providing a blueprint for their regulatory pathways. This study provided a cross-genome comparative analysis of TTF genes in 110 species, which contributed to understanding their copy number expansion and evolutionary footprint in plants.

Identification and analysis of differentially expressed trihelix genes in maize (Zea mays) under abiotic stresses

Background

Trihelix transcription factors play important roles in triggering plant growth and imparting tolerance against biotic and abiotic stresses. However, a systematical analysis of the trihelix transcription factor family under heat and drought stresses in maize has not been reported.

Methods

PlantTFDB and TBtools were employed to identify the trihelix domain-containing genes in the maize genome. The heat-regulated transcriptome data for maize were obtained from NCBI to screen differentially expressed ZmTHs genes through statistical analysis. The basic protein sequences, chromosomal localization, and subcellular localization were analyzed using Maize GDB, Expasy, SOMPA, TBtools, and Plant-mPLoc. The conserved motifs, evolutionary relationships, and cis-elements, were analyzed by MEME, MEGA7.0 and PlantCARE software, respectively. The tissue expression patterns of ZmTHs and their expression profiles under heat and drought stress were detected using quantitative real-time PCR (qRT-PCR).

Results

A total of 44 trihelix family members were discovered, and members were distributed over 10 chromosomes in the maize genome. A total of 11 genes were identified that were regulated by heat stress; these were unevenly distributed on chromosomes 1, 2, 4, 5, and 10. ZmTHs encoded a total of 16 proteins, all of which were located in the nucleus; however, ZmTH04.1 was also distributed in the chloroplast. The protein length varied from 206 to 725 amino acids; the molecular weight ranged from 22.63 to 76.40 kD; and the theoretical isoelectric point (pI) ranged from 5.24 to 11.2. The protein's secondary structures were mainly found to be random coils and α-helices, with fewer instances of elongation chains and β-rotations. Phylogenetic relationship analysis showed that these can be divided into five sub-groups. The conserved domain of ZmTHs was GT1 or MyB_DNA-Bind_4. The protein and gene structure of ZmTHs differed greatly among the subfamilies, while the structures within the subfamilies were similar. The promoter of ZmTHs contained abundant tissue-specific expression cis-acting elements and abiotic stress response elements. qRT-PCR analysis showed that ZmTHs expression levels were significantly different in different tissues. Furthermore, the expression of ZmTH08 was dramatically up-regulated by heat stress, while the expression of ZmTH03, ZmTH04, ZmTH05, ZmTH06, ZmTH07, ZmTH09, ZmTH10, and ZmTH11 were down-regulated by heat stress. Upon PEG-simulated drought stress, ZmTH06 was significantly up-regulated, while ZmTH01 and ZmTH07 were down-regulated.

Conclusions

We performed a genome-wide, systematic identification and analysis of differentially expressed trihelix genes under heat and drought stresses in maize.

Conservation and Divergence of the Trihelix Genes in Brassica and Expression Profiles of BnaTH Genes in Brassica napus under Abiotic Stresses

Trihelix (TH) proteins are a family of plant-specific transcription factors that play a role in light response and are extensively involved in plant growth and development, as well as in various stress responses. However, the function of TH genes in Brassica napus (B. napus) remains unclear, as does the evolution and differentiation pattern of TH genes in Brassica plants. Here, we identified a total of 455 TH genes in seven species, including six Brassica species and Arabidopsis, which were grouped into five clades, GT-1, GT-2, GTγ, SH4, and SIP1, each with 69, 142, 44, 55, and 145 members, respectively. The types and distributions of motifs of the TH proteins and the structures of the TH genes are conserved in the same subgroup, and some variations in certain amino acid residues occur in B. napus when inheriting motifs from Brassica rapa (B. rapa) and Brassica oleracea (B. oleracea). Collinearity analysis revealed that the massive expansion of TH genes in tetraploid species was attributed to the hetero-tetraploidization of diploid ancestors and gene duplication events within the tetraploid species. Comparative analysis of the membership numbers of five subgroups in different species revealed that the GT-2 and SIP1 genes underwent significant expansion during evolution, possibly to support the better adaptation of plants to their environments. The differential expression of the BnaTH genes under five stresses indicates that the BnaTH genes are involved in plant responses to stresses such as drought, cold, and heat. The presence of different stress-responsive cis-elements in the upstream promoter region of the genes indicated that BnaTH genes have the potential to cope with variable environments. Meanwhile, qRT-PCR analyses also confirmed that five TH genes respond to different abiotic stresses. Our results provide information and candidates for further studies on the role of TH genes in stress resistance of B. napus.

Arginine Decarboxylase Gene ADC2 Regulates Fiber Elongation in Cotton

Cotton is an important agro-industrial crop providing raw material for the textile industry. Fiber length is the key factor that directly affects fiber quality. ADC, arginine decarboxylase, is the key rate-limiting enzyme in the polyamine synthesis pathway; whereas, there is no experimental evidence that ADC is involved in fiber development in cotton yet. Our transcriptome analysis of the fiber initiation material of Gossypium arboreum L. showed that the expression profile of GaADC2 was induced significantly. Here, GhADC2, the allele of GaADC2 in tetraploid upland cotton Gossypium hirsutum L., exhibited up-regulated expression pattern during fiber elongation in cotton. Levels of polyamine are correlated with fiber elongation; especially, the amount of putrescine regulated by ADC was increased. Scanning electron microscopy showed that the fiber length was increased with exogenous addition of an ADC substrate or product putrescine; whereas, the fiber density was decreased with exogenous addition of an ADC specific inhibitor. Next, genome-wide transcriptome profiling of fiber elongation with exogenous putrescine addition was performed to determine the molecular basis in Gossypium hirsutum. A total of 3163 differentially expressed genes were detected, which mainly participated in phenylpropanoid biosynthesis, fatty acid elongation, and sesquiterpenoid and triterpenoid biosynthesis pathways. Genes encoding transcription factors MYB109, WRKY1, and TCP14 were enriched. Therefore, these results suggested the ADC2 and putrescine involvement in the development and fiber elongation of G. hirsutum, and provides a basis for cotton fiber development research in future.

Protoplast Dissociation and Transcriptome Analysis Provides Insights to Salt Stress Response in Cotton

As one of the pioneer crops widely planted in saline-alkaline areas, Gossypium provides daily necessities, including natural fiber, vegetable proteins, and edible oils. However, cotton fiber yield and quality are highly influenced by salt stress. Therefore, elucidating the molecular mechanisms of cotton in response to salinity stress is importance to breed new cultivars with high tolerance. In this study, we first developed a method for single-cell RNA-seq based on isolating protoplast from cotton root tips; then, we studied the impact of salinity stress on gene expression profiling and their dynamic changes using the developed high-efficiency method for protoplast dissociation suitable for single-cell RNA-seq. A total of 3391 and 2826 differentially expressed genes (DEGs) were identified in salt-treated samples before and after protoplast dissociation, respectively, which were enriched into several molecular components, including response to stimulus, response to stress, and cellular macromolecule metabolic process by gene ontology (GO) analysis. Plant hormone signal transduction, phenylpropanoid biosynthesis, and MAPK signaling pathway were found to be enriched via Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Twenty-two and nine salinity-responsive DEGs participated in plant hormone signaling and MAPK signaling in roots, before and after protoplast dissociation, respectively; six upregulated DEGs were involved in ABA signaling transduction, namely, Ga04G2111, Ga07G0142, Ga09G2061, Ga10G0262, Ga01G0063, and Ga08G1915 which indicates their potential functions on plants adapting to salt stress. Additionally, 384 and 257 transcription factors (TFs) were differentially expressed in salt-stress roots before and after protoplast dissociation, respectively, of which significantly up-regulated TFs mainly belonged to the AP2/ERF-ERF family, which implied their potential roles responding to salt stress. These results not only provide novel insights to reveal the regulatory networks in plant’s root response to salt stress, but also lay the solid foundation for further exploration on cellular heterogeneity by single-cell transcriptome sequencing.

Genome-Wide Identification and Expression Profiling Analysis of the Trihelix Gene Family Under Abiotic Stresses in Medicago truncatula

The trihelix transcription factor (GT) family is widely involved in regulating plant growth and development, and most importantly, responding to various abiotic stresses. Our study first reported the genome-wide identification and analysis of GT family genes in Medicago truncatula. Overall, 38 trihelix genes were identified in the M. truncatula genome and were classified into five subfamilies (GT-1, GT-2, SH4, GTγ and SIP1). We systematically analyzed the phylogenetic relationship, chromosomal distribution, tandem and segmental duplication events, gene structures and conserved motifs of MtGTs. Syntenic analysis revealed that trihelix family genes in M. truncatula had the most collinearity relationship with those in soybean followed by alfalfa, but very little collinearity with those in the maize and rice. Additionally, tissue-specific expression analysis of trihelix family genes suggested that they played various roles in the growth and development of specific tissues in M. truncatula. Moreover, the expression of some MtGT genes, such as MtGT19, MtGT20, MtGT22, and MtGT33, was dramatically induced by drought, salt, and ABA treatments, illustrating their vital roles in response to abiotic stresses. These findings are helpful for improving the comprehensive understanding of trihelix family; additionally, the study provides candidate genes for achieving the genetic improvement of stress resistance in legumes.

Transcriptome Analysis Reveals a Gene Expression Pattern Associated with Fuzz Fiber Initiation Induced by High Temperature in Gossypium barbadense

Gossypium barbadense is an important source of natural textile fibers, as is Gossypium hirsutum. Cotton fiber development is often affected by various environmental factors, such as abnormal temperature. However, little is known about the underlying mechanisms of temperature regulating the fuzz fiber initiation. In this study, we reveal that high temperatures (HT) accelerate fiber development, improve fiber quality, and induced fuzz initiation of a thermo-sensitive G. barbadense variety L7009. It was proved that fuzz initiation was inhibited by low temperature (LT), and 4 dpa was the stage most susceptible to temperature stress during the fuzz initiation period. A total of 43,826 differentially expressed genes (DEGs) were identified through comparative transcriptome analysis. Of these, 9667 were involved in fiber development and temperature response with 901 transcription factor genes and 189 genes related to plant hormone signal transduction. Further analysis of gene expression patterns revealed that 240 genes were potentially involved in fuzz initiation induced by high temperature. Functional annotation revealed that the candidate genes related to fuzz initiation were significantly involved in the asparagine biosynthetic process, cell wall biosynthesis, and stress response. The expression trends of sixteen genes randomly selected from the RNA-seq data were almost consistent with the results of qRT-PCR. Our study revealed several potential candidate genes and pathways related to fuzz initiation induced by high temperature. This provides a new view of temperature-induced tissue and organ development in Gossypium barbadense.