Genome-wide characterization of the GRF family and their roles in response to salt stress in Gossypium

AaMYC3 bridges the regulation of glandular trichome density and artemisinin biosynthesis in Artemisia annua

Artemisinin, the well-known natural product for treating malaria, is biosynthesised and stored in the glandular-secreting trichomes (GSTs) of Artemisia annua. While numerous efforts have clarified artemisinin metabolism and regulation, the molecular association between artemisinin biosynthesis and GST development remains elusive. Here, we identified AaMYC3, a bHLH transcription factor of A. annua, induced by jasmonic acid (JA), which simultaneously regulates GST density and artemisinin biosynthesis. Overexpressing AaMYC3 led to a substantial increase in GST density and artemisinin accumulation. Conversely, in the RNAi-AaMYC3 lines, both GST density and artemisinin content were markedly reduced. Through RNA-seq and analyses conducted both in vivo and in vitro, AaMYC3 not only directly activates AaHD1 transcription, initiating GST development, but also up-regulates the expression of artemisinin biosynthetic genes, including CYP71AV1 and ALDH1, thereby promoting artemisinin production. Furthermore, AaMYC3 acts as a co-activator, interacting with AabHLH1 and AabHLH113, to trigger the transcription of two crucial enzymes in the artemisinin biosynthesis pathway, ADS and DBR2, ultimately boosting yield. Our findings highlight a critical connection between GST initiation and artemisinin biosynthesis in A. annua, providing a new target for molecular design breeding of traditional Chinese medicine.

Time-Course Transcriptomics Analysis Reveals Molecular Mechanisms of Salt-Tolerant and Salt-Sensitive Cotton Cultivars in Response to Salt Stress

Salt stress is an environmental factor that limits plant seed germination, growth, and survival. We performed a comparative RNA sequencing transcriptome analysis during germination of the seeds from two cultivars with contrasting salt tolerance responses. A transcriptomic comparison between salt-tolerant cotton cv Jin-mian 25 and salt-sensitive cotton cv Su-mian 3 revealed both similar and differential expression patterns between the two genotypes during salt stress. The expression of genes related to aquaporins, kinases, reactive oxygen species (ROS) scavenging, trehalose biosynthesis, and phytohormone biosynthesis and signaling that include ethylene (ET), gibberellin (GA), abscisic acid (ABA), jasmonic acid (JA), and brassinosteroid (BR) were systematically investigated between the cultivars. Despite the involvement of these genes in cotton's response to salt stress in positive or negative ways, their expression levels were mostly similar in both genotypes. Interestingly, a PXC2 gene (Ghir_D08G025150) was identified, which encodes a leucine-rich repeat receptor-like protein kinase (LRR-RLK). This gene showed an induced expression pattern after salt stress treatment in salt-tolerant cv Jin-mian 25 but not salt-sensitive cv Su-mian 3. Our multifaceted transcriptome approach illustrated a differential response to salt stress between salt-tolerant and salt-sensitive cotton.

In silico genome-wide analysis of the growth-regulating factor gene family and their expression profiling in Vitis vinifera under biotic stress

Growth regulatory factors (GRFs) are transcription factors that encode the proteins involved in plant growth and development. However, no comprehensive analysis of Vitis vinifera GRF genes has yet been conducted. In the current study, we performed a genome-wide analysis of the GRF gene family to explore the VvGRF gene's role in Vitis vinifera. We identified 30 VvGRF genes in the Vitis vinifera genome, localized over 20 chromosomes. Based on evolutionary analysis, 49 GRF genes (nine AtGRF, ten FvGRF, and 30 VvGRF) were clustered into six groups. Many cis-elements involved in light control, defense, and plant growth have been identified in the promoter region of VvGRF genes, and multiple miRNAs have been predicted to be involved in regulating VvGRF gene expression. Protein-protein interaction analysis showed that nine VvGRF proteins formed a complex protein interaction network. Furthermore, the gene expression analysis of VvGRF revealed that VvGRF-5 and VvGRF-6 were highly upregulated suggesting that these genes are involved in biotic responses. This study provides comprehensive insights into the functional characteristics and occurrence of the VvGRF gene family in Vitis vinifera, which may be applied in breeding programs to enhance the growth of Vitis vinifera varieties under stress and growth changes.

Genome-wide identification and expression profiling of growth‑regulating factor (GRF) and GRF‑interacting factor (GIF) gene families in chickpea and pigeonpea

The growth-regulating factor (GRF) and GRF-interacting factor (GIF) families encode plant-specific transcription factors and play vital roles in plant development and stress response processes. Although GRF and GIF genes have been identified in various plant species, there have been no reports of the analysis and identification of the GRF and GIF transcription factor families in chickpea (Cicer arietinum) and pigeonpea (Cajanus cajan). The present study identified seven CaGRFs, eleven CcGRFs, four CaGIFs, and four CcGIFs. The identified proteins were grouped into eight and three clades for GRFs and GIFs, respectively based on their phylogenetic relationships. A comprehensive in-silico analysis was performed to determine chromosomal location, sub-cellular localization, and types of regulatory elements present in the putative promoter region. Synteny analysis revealed that GRF and GIF genes showed diploid-polyploid topology in pigeonpea, but not in chickpea. Tissue-specific expression data at the vegetative and reproductive stages of the plant showed that GRFs and GIFs were strongly expressed in tissues like embryos, pods, and seeds, indicating that GRFs and GIFs play vital roles in plant growth and development. This research characterized GRF and GIF families and hints at their primary roles in the chickpea and pigeonpea growth and developmental process. Our findings provide potential gene resources and vital information on GRF and GIF gene families in chickpea and pigeonpea, which will help further understand the regulatory role of these gene families in plant growth and development.

Expression profile analysis of cotton fiber secondary cell wall thickening stage

To determine the genes associated with the fiber strength trait in cotton, three different cotton cultivars were selected: Sea Island cotton (Xinhai 32, with hyper-long fibers labeled as HL), and upland cotton (17-24, with long fibers labeled as L, and 62-33, with short fibers labeled as S). These cultivars were chosen to assess fiber samples with varying qualities. RNA-seq technology was used to analyze the expression profiles of cotton fibers at the secondary cell wall (SCW) thickening stage (20, 25, and 30 days post-anthesis (DPA)). The results showed that a large number of differentially expressed genes (DEGs) were obtained from the three assessed cotton cultivars at different stages of SCW development. For instance, at 20 DPA, Sea Island cotton (HL) had 6,215 and 5,364 DEGs compared to upland cotton 17-24 (L) and 62-33 (S), respectively. Meanwhile, there were 1,236 DEGs between two upland cotton cultivars, 17-24 (L) and 62-33 (S). Gene Ontology (GO) term enrichment identified 42 functions, including 20 biological processes, 11 cellular components, and 11 molecular functions. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis identified several pathways involved in SCW synthesis and thickening, such as glycolysis/gluconeogenesis, galactose metabolism, propanoate metabolism, biosynthesis of unsaturated fatty acids pathway, valine, leucine and isoleucine degradation, fatty acid elongation pathways, and plant hormone signal transduction. Through the identification of shared DEGs, 46 DEGs were found to exhibit considerable expressional differences at different fiber stages from the three cotton cultivars. These shared DEGs have functions including REDOX enzymes, binding proteins, hydrolases (such as GDSL thioesterase), transferases, metalloproteins (cytochromatin-like genes), kinases, carbohydrates, and transcription factors (MYB and WRKY). Therefore, RT-qPCR was performed to verify the expression levels of nine of the 46 identified DEGs, an approach which demonstrated the reliability of RNA-seq data. Our results provided valuable molecular resources for clarifying the cell biology of SCW biosynthesis during fiber development in cotton.

Genome-wide identification of the GRF family in sweet orange (Citrus sinensis) and functional analysis of the CsGRF04 in response to multiple abiotic stresses

BACKGROUND

Citrus is one of the most valuable fruits worldwide and an economic pillar industry in southern China. Nevertheless, it frequently suffers from undesirable environmental stresses during the growth cycle, which severely restricts the growth, development and yield of citrus. In plants, the growth-regulating factor (GRF) family of transcription factors (TF) is extensively distributed and plays an vital part in plant growth and development, hormone response, as well as stress adaptation. However, the systematic identification and functional analysis of GRF TFs in citrus have not been reported.

RESULTS

Here, a genome-wide identification of GRF TFs was performed in Citrus sinensis, 9 members of CsGRFs were systematically identified and discovered to be scattered throughout 5 chromosomes. Subsequently, physical and chemical properties, phylogenetic relationships, structural characteristics, gene duplication events, collinearity and cis-elements of promoter were elaborately analyzed. In particular, the expression patterns of the CsGRF genes in response to multiple phytohormone and abiotic stress treatments were investigated. Predicated on this result, CsGRF04, which exhibited the most differential expression pattern under multiple phytohormone and abiotic stress treatments was screened out. Virus-induced gene silencing (VIGS) technology was utilized to obtain gene silenced plants for CsGRF04 successfully. After the three stress treatments of high salinity, low temperature and drought, the CsGRF04-VIGS lines showed significantly reduced resistance to high salinity and low temperature stresses, but extremely increased resistance to drought stress.

CONCLUSIONS

Taken together, our findings systematically analyzed the genomic characterization of GRF family in Citrus sinensis, and excavated a CsGRF04 with potential functions under multiple abiotic stresses. Our study lay a foundation for further study on the function of CsGRFs in abiotic stress and hormone signaling response.

Genome-wide identification of GROWTH-REGULATING FACTORs in Liriodendron chinense and functional characterization of LcGRF2 in leaf size regulation

GROWTH-REGULATING FACTORs (GRFs) play a pivotal role in the regulation of leaf size in plants and have been widely reported in plants. However, their specific functions in leaf size regulation in Liriodendron chinense remains unclear. Therefore, in this study, we identified GRF genes on a genome-wide scale in L. chinense to characterize the roles of LcGRFs in regulating leaf size. A total of nine LcGRF genes were identified, and these genes exhibited weak expression in mature leaves but strong expression in shoot apex. Notably, LcGRF2 exhibited the highest expression level in the shoot apex of L. chinense. Further RT-qPCR assay revealed that the expression level of LcGRF2 gradually decreased along with the leaf development process, and also displayed a gradient along the leaf proximo-distal and medio-lateral axes. Furthermore, overexpression of LcGRF2 in Arabidopsis thaliana resulted in increased leaf size, and significantly up-regulated the expression of genes involved in cell division like AtCYCD3;1, AtKNOLLE, and AtCYCB1;1, indicating that LcGRF2 may influence leaf size by promoting cell proliferation. This work contributes to a better understanding of the roles and molecular mechanisms of LcGRFs in the regulation of leaf size in L. chinense.

Genome-wide identification and expression analysis of growth-regulating factors in Dendrobium officinale and Dendrobium chrysotoxum

Background

Dendrobium, one of the largest genera in Orchidaceae, is popular not only for its aesthetic appeal but for its significant medicinal value. Growth-regulating factors (GRFs) play an essential role in plant growth and development. However, there is still a lack of information about the evolution and biological function analysis of the GRF gene family among Dendrobiumspecies.

Methods

Growth-regulating factors from Dendrobium officinale Kimura et Migo and Dendrobium chrysotoxum Lindl. were identified by HMMER and BLAST. Detailed bioinformatics analysis was conducted to explore the evolution and function of GRF gene family in D. officinale and D. chrysotoxum using genomic data, transcriptome data and qRT-PCR technology.

Results

Here, we evaluated the evolution of the GRF gene family based on the genome sequences of D. officinale and D. chrysotoxum. Inferred from phylogenetic trees, the GRF genes were classified into two clades, and each clade contains three subclades. Sequence comparison analysis revealed relatively conserved gene structures and motifs among members of the same subfamily, indicating a conserved evolution of GRF genes within Dendrobiumspecies. However, considering the distribution of orthologous DoGRFs and DcGRFs, and the differences in the number of GRFs among species, we suggest that the GRF gene family has undergone different evolutionary processes. A total of 361 cis-elements were detected, with 33, 141, and 187 related to plant growth and development, stress, and hormones, respectively. The tissue-specific expression of GRFs showed that DoGRF8 may have a significant function in the stem elongation of D. officinale. Moreover, four genes were up-regulated under Methyl-jasmonic acid/methyl jasmonate (MeJA) treatment, showing that DoGRFs and DcGRFs play a crucial role in stress response. These findings provide valuable information for further investigations into the evolution and function of GRF genes in D. officinale and D. chrysotoxum.

Genome-Wide Analysis of Flax (Linum usitatissimum L.) Growth-Regulating Factor (GRF) Transcription Factors

Flax is an important cash crop globally with a variety of commercial uses. It has been widely used for fiber, oil, nutrition, feed and in composite materials. Growth regulatory factor (GRF) is a transcription factor family unique to plants, and is involved in regulating many processes of growth and development. Bioinformatics analysis of the GRF family in flax predicted 17 LuGRF genes, which all contained the characteristic QLQ and WRC domains. Equally, 15 of 17 LuGRFs (88%) are predicted to be regulated by lus-miR396 miRNA. Phylogenetic analysis of GRFs from flax and several other well-characterized species defined five clades; LuGRF genes were found in four clades. Most LuGRF gene promoters contained cis-regulatory elements known to be responsive to hormones and stress. The chromosomal locations and collinearity of LuGRF genes were also analyzed. The three-dimensional structure of LuGRF proteins was predicted using homology modeling. The transcript expression data indicated that most LuGRF family members were highly expressed in flax fruit and embryos, whereas LuGRF3, LuGRF12 and LuGRF16 were enriched in response to salt stress. Real-time quantitative fluorescent PCR (qRT-PCR) showed that both LuGRF1 and LuGRF11 were up-regulated under ABA and MeJA stimuli, indicating that these genes were involved in defense. LuGRF1 was demonstrated to be localized to the nucleus as expected for a transcription factor. These results provide a basis for further exploration of the molecular mechanism of LuGRF gene function and obtaining improved flax breeding lines.

Growth-regulating factors: conserved and divergent roles in plant growth and development and potential value for crop improvement

High yield and stress resistance are the major prerequisites for successful crop cultivation, and can be achieved by modifying plant architecture. Evolutionarily conserved growth-regulating factors (GRFs) control the growth of different tissues and organs of plants. Here, we provide a systematic overview of the expression patterns of GRF genes and the structural features of GRF proteins in different plant species. Moreover, we illustrate the conserved and divergent roles of GRFs, microRNA396 (miR396), and GRF-interacting factors (GIFs) in leaf, root, and flower development. We also describe the molecular networks involving the miR396-GRF-GIF module, and illustrate how this module coordinates with different signaling molecules and transcriptional regulators to control development of different plant species. GRFs promote leaf growth, accelerate grain filling, and increase grain size and weight. We also provide some molecular insight into how coordination between GRFs and other signaling modules enhances crop productivity; for instance, how the GRF-DELLA interaction confers yield-enhancing dwarfism while increasing grain yield. Finally, we discuss how the GRF-GIF chimera substantially improves plant transformation efficiency by accelerating shoot formation. Overall, we systematically review the conserved and divergent roles of GRFs and the miR396-GRF-GIF module in growth regulation, and also provide insights into how GRFs can be utilized to improve the productivity and nutrient content of crop plants.

Genome-wide identification of GRF gene family and their contribution to abiotic stress response in pitaya (Hylocereus polyrhizus)

Growth-regulating factors (GRFs) are plant-specific transcription factors identified in many land plants. Recently, their indispensable roles in stress response are highlighted. In present work, 11 HpGRFs were cloned in pitaya. Segmental duplication is considered essential for the expansion of HpGRFs. A phylogenetic tree suggested that GRFs could be divided into eight categories, among which G-I was a Caryophyllales-specific one. The categorization was further evidenced by differences in the gene structure, collinearity, protein domain of HpGRFs. Five miR396 hairpins giving rise to two types of matured miR396s were identified in pitaya via sRNA-Seq in combination with bioinformatic analysis. Parallel analysis of RNA ends proved that HpGRFs except HpGRF5 were degraded by miR396-directed cleavages at the regions which code the conserved WRC motifs of HpGRFs. Multiple cis-regulatory elements were discovered in the promoters of HpGRFs. Among the elements, most are involved in stress and phytohormone response as well as plant growth, indicating a crosstalk between them. Expression analysis showed the responsive patterns of the miR396-GRF module under abiotic stresses. To conclude, our work systematically identified the miR396-targeted HpGRFs in pitaya and confirmed their involvement in stress response, providing novel insights into the comprehensive understanding of the stress resistance of pitaya.

Comparative Genomics and Functional Studies of Putative m6A Methyltransferase (METTL) Genes in Cotton

N6-methyladenosine (m6A) RNA modification plays important regulatory roles in plant development and adapting to the environment, which requires methyltransferases to achieve the methylation process. However, there has been no research regarding m6A RNA methyltransferases in cotton. Here, a systematic analysis of the m6A methyltransferase (METTL) gene family was performed on twelve cotton species, resulting in six METTLs identified in five allotetraploid cottons, respectively, and three to four METTLs in the seven diploid species. Phylogenetic analysis of protein-coding sequences revealed that METTL genes from cottons, Arabidopsis thaliana, and Homo sapiens could be classified into three clades (METTL3, METTL14, and METTL-like clades). Cis-element analysis predicated the possible functions of METTL genes in G. hirsutum. RNA-seq data revealed that GhMETTL14 (GH_A07G0817/GH_D07G0819) and GhMETTL3 (GH_A12G2586/GH_D12G2605) had high expressions in root, stem, leaf, torus, petal, stamen, pistil, and calycle tissues. GhMETTL14 also had the highest expression in 20 and 25 dpa fiber cells, implying a potential role at the cell wall thickening stage. Suppressing GhMETTL3 and GhMETTL14 by VIGS caused growth arrest and even death in G. hirsutum, along with decreased m6A abundance from the leaf tissues of VIGS plants. Overexpression of GhMETTL3 and GhMETTL14 produced distinct differentially expressed genes (DEGs) in A. thaliana, indicating their possible divergent functions after gene duplication. Overall, GhMETTLs play indispensable but divergent roles during the growth of cotton plants, which provides the basis for the systematic investigation of m6A in subsequent studies to improve the agronomic traits in cotton.

Genome-wide characterization and expression analysis of the growth-regulating factor family in Saccharum

BACKGROUND

Growth regulating factors (GRFs) are transcription factors that regulate diverse biological and physiological processes in plants, including growth, development, and abiotic stress. Although GRF family genes have been studied in a variety of plant species, knowledge about the identification and expression patterns of GRFs in sugarcane (Saccharum spp.) is still lacking.

RESULTS

In the present study, a comprehensive analysis was conducted in the genome of wild sugarcane (Saccharum spontaneum) and 10 SsGRF genes were identified and characterized. The phylogenetic relationship, gene structure, and expression profiling of these genes were analyzed entirely under both regular growth and low-nitrogen stress conditions. Phylogenetic analysis suggested that the 10 SsGRF members were categorized into six clusters. Gene structure analysis indicated that the SsGRF members in the same group were greatly conserved. Expression profiling demonstrated that most SsGRF genes were extremely expressed in immature tissues, implying their critical roles in sugarcane growth and development. Expression analysis based on transcriptome data and real-time quantitative PCR verification revealed that GRF1 and GRF3 were distinctly differentially expressed in response to low-nitrogen stress, which meant that they were additional participated in sugarcane stress tolerance.

CONCLUSION

Our study provides a scientific basis for the potential functional prediction of SsGRF and will be further scrutinized by examining their regulatory network in sugarcane development and abiotic stress response, and ultimately facilitating their application in cultivated sugarcane breeding.

Genome-Wide Analysis of the Growth-Regulating Factor (GRF) Family in Aquatic Plants and Their Roles in the ABA-Induced Turion Formation of Spirodela polyrhiza

Growth-regulating factors (GRFs) are plant-specific transcription factors that play essential roles in regulating plant growth and stress response. The GRF gene families have been described in several terrestrial plants, but a comprehensive analysis of these genes in diverse aquatic species has not been reported yet. In this study, we identified 130 GRF genes in 13 aquatic plants, including floating plants (Azolla filiculoides, Wolffia australiana, Lemna minuta, Spirodela intermedia, and Spirodela polyrhiza), floating-leaved plants (Nymphaea colorata and Euryale ferox), submersed plants (Zostera marina, Ceratophyllum demersum, Aldrovanda vesiculosa, and Utricularia gibba), an emergent plant (Nelumbo nucifera), and an amphibious plant (Cladopus chinensis). The gene structures, motifs, and cis-acting regulatory elements of these genes were analyzed. Phylogenetic analysis divided these GRFs into five clusters, and ABRE cis-elements were highly enriched in the promoter region of the GRFs in floating plants. We found that abscisic acid (ABA) is efficient at inducing the turion of Spirodela polyrhiza (giant duckweed), accompanied by the fluctuated expression of SpGRF genes in their fronds. Our results provide information about the GRF gene family in aquatic species and lay the foundation for future studies on the functions of these genes.

Genome-wide identification of growth-regulating factor transcription factor family related to leaf and stem development in alfalfa

Growth-regulating factors (GRFs) play crucial roles in plant growth and stress response. To date, there have been no reports of the analysis and identification of the GRF transcription factor family in alfalfa. In this study, we identified 27 GRF family members from alfalfa (Medicago sativa L.) "Xinjiang Daye", and analyzed their physicochemical properties. Based on phylogenetic analysis, these MsGRFs were divided into five subgroups, each with a similar gene structure and conserved motifs. MsGRFs genes are distributed on 23 chromosomes, and all contain QLQ and WRC conserved domains. The results of the collinearity analysis showed that all MsGRFs are involved in gene duplication, including multiple whole-genome duplication or segmental duplication and a set of tandem duplication, indicating that large-scale duplication is important for the expansion of the GRF family in alfalfa. Several hormone-related and stress-related cis-acting elements have been found in the promoter regions of MsGRFs. Some MsGRFs were highly expressed in young leaves and stems, and their expression decreased during development. In addition, the leaf size of different varieties was found to vary, and MsGRF1 to 4, MsGRF18 to 20, and MsGRF22 to 23 were differentially expressed in large and small leaf alfalfa varieties, suggesting that they are critical in the regulation of leaf size. The results of this study can benefit further exploration of the regulatory functions of MsGRFs in growth and development, and can identify candidate genes that control leaf size development.

TRY intron2 determined its expression in inflorescence activated by SPL9 and MADS-box genes in Arabidopsis

Plant trichomes are specialized epidermal cells that protect plants from insects and pathogens. In Arabidopsis, epidermal hairs decrease as internodes increase in height, with only few epidermal hairs produced on the sepals abaxial surface of the early flowers. TRIPTYCHON (TRY) is known to be a negative regulator of epidermal hair development in Arabidopsis, suppressing the formation of ectopic epidermal hairs in the inflorescence. Here, we reported that the second intron of TRY gene plays a critical role in trichome spatial distribution in Arabidopsis. The expression of TRY rises with the increasing stem nodes and reaches the peak in the inflorescence, while the trichomes distribution decrease. The transgenic plants showed that TRY promoter could only drive the genomic instead of coding sequences combined with GUS reporter gene, which indicates that the regulatory elements of TRY expression in inflorescence could be located in the intron regions. Multiple SPLs and MADS-box binding sites were found in the TRY intron2 sequence. Further genetic and biochemistry assays revealed that the flowering-related genes such as SPL9 could bind to these cis-elements directly, contributing to the TRY spatial expression. Since cotton fiber and Arabidopsis trichomes share similar regulatory mechanism, extended analysis showed that the intron2 of cotton TRY genes also contain the cis-elements. Thus, the introns harboring the transcription element may be the general way to regulate the gene expression in different plants and provides molecular clues for the related crops' traits design.

Comparative transcriptome analysis unveiling reactive oxygen species scavenging system of Sonneratia caseolaris under salinity stress

Many mangrove forests have undergone major changes as a result of human activity and global climate change. Sonneratia caseolaris is a common tree located in inner mangroves, and its range extends inland along tidal creeks, as far as the influence of salinity extends. This study investigated the physiological and molecular response mechanisms of S. caseolaris by analyzing its antioxidant defense capacity, including its differentially expressed genes (DEGs) under similar salt stress conditions. Salt treatment significantly affected the osmoprotectants and lipid peroxidation in S. caseolaris seedlings, which increased proline (Pro) content by 31.01-54.90% during all sample periods and decreased malonaldehyde (MDA) content by 12.81 and 18.17% at 25 and 40 days under 3.0% NaCl treatment. Antioxidant enzyme activities increased significantly following 3.0% NaCl treatment. Transcriptome analysis following De novo assembly showed 26,498 matched unigenes. The results showed that 1,263 DEGs responded to transcription factors (TFs) and plant phytohormones and mediated oxidoreductase activity to scavenge reactive oxygen species (ROS) in the control vs. 3.0% NaCl comparison. In addition, the transcription levels of genes associated with auxin and ethylene signal transduction also changed. Under salt stress, ROS scavenging genes (POD, CAT, and APX) and part of AP2, MYB, NAC, C2C2, bHLH, and WRKY TFs were upregulated. This study identified important pathways and candidate genes involved in S. caseolaris salinity tolerance and provided suggestions for further research into the mechanisms of salt tolerance in S. caseolaris.

Systematical characterization of GRF gene family in sorghum, and their potential functions in aphid resistance

Sorghum (Sorghum bicolor) is the fifth important cereal and an industrial energy crop in the world. Growth Regulation Factors (GRFs) play an important role in response to environmental stress, however, the knowledge of GRFs relating to the pest resistance is lacking. Here, we identified 8 GRF genes harboring the typical QLQ (glutamine, leucine, glutamine) and WRC (tryptophan, arginine, cysteine) domains in Sorghum, which could be classified into 4 clades through phylogenetic analysis. The SbGRF genes express in most tissues, while more than half of them express at the highest level in inflorescence. To further investigate their possible role in stress response, we analyzed the transcriptomics data. The results showed that SbGRFs could respond to the abiotic stresses including heat, salt and drought stress. Furthermore, combined the data with qRT-PCR, SbGRF1, 2, 4 and 7 were identified as dominant genes response to the aphid-induced stress. SSR markers close to these genes were also searched. Above all, we summarized the SbGRFs and provided their potential roles in aphid response.

Comprehensive identification and expression analysis of CRY gene family in Gossypium

BACKGROUND

The cryptochromes (CRY) are specific blue light receptors of plants and animals, which play crucial roles in physiological processes of plant growth, development, and stress tolerance.

RESULTS

In the present work, a systematic analysis of the CRY gene family was performed on twelve cotton species, resulting in 18, 17, 17, 17, and 17 CRYs identified in five alloteraploid cottons (Gossypium hirsutum, G. barbadense, G. tomentosum, G. mustelinum and G. darwinii), respectively, and five to nine CRY genes in the seven diploid species. Phylogenetic analysis of protein-coding sequences revealed that CRY genes from cottons and Arabidopsis thaliana could be classified into seven clades. Synteny analysis suggested that the homoeolog of G. hirsutum Gh_A02G0384 has undergone an evolutionary loss event in the other four allotetraploid cotton species. Cis-element analysis predicated the possible functions of CRY genes in G. hirsutum. RNA-seq data revealed that Gh_D09G2225, Gh_A09G2012 and Gh_A11G1040 had high expressions in fiber cells of different developmental states. In addition, the expression levels of one (Gh_A03G0120), 15 and nine GhCRY genes were down-regulated following the PEG, NaCl and high-temperature treatments, respectively. For the low-temperature treatment, five GhCRY genes were induced, and five were repressed. These results indicated that most GhCRY genes negatively regulate the abiotic stress treatments.

CONCLUSION

We report the structures, domains, divergence, synteny, and cis-elements analyses systematically of G. hirsutum CRY genes. Possible biological functions of GhCRY genes in differential tissues as well as in response to abiotic stress during the cotton plant life cycle were predicted.

Genome-Wide Identification and Expression Analysis of the 14-3-3 Gene Family in Mango (Mangifera indica L.)

Members of the Mi14-3-3 gene family interact with target proteins that are widely involved in plant hormone signal transduction and physiology-related metabolism and play important roles in plant growth, development and stress responses. In this study, 14-3-3s family members are identified by the bioinformatic analysis of the mango (Mangifera indica L.) genome. The gene structures, chromosomal distributions, genetic evolution, and expression patterns of these genes and the physical and chemical properties and conserved motifs of their proteins are analysed systematically. The results identified 16 members of the 14-3-3 genes family in the mango genome. The members were not evenly distributed across the chromosomes, and the gene structure analysis showed that the gene sequence length and intron number varied greatly among the different members. Protein sequence analysis showed that the Mi14-3-3 proteins had similar physical and chemical properties and secondary and tertiary structures, and protein subcellular localization showed that the Mi14-3-3 family proteins were localized to the nucleus. The sequence analysis of the Mi14-3-3s showed that all Mi14-3-3 proteins contain a typical conserved PFAM00244 domain, and promoter sequence analysis showed that the Mi14-3-3 promoters contain multiple hormone-, stress-, and light-responsive cis-regulatory elements. Expression analysis showed that the 14-3-3 genes were expressed in all tissues of mango, but that their expression patterns were different. Drought, salt and low temperature stresses affected the expression levels of 14-3-3 genes, and different 14-3-3 genes had different responses to these stresses. This study provides a reference for further studies on the function and regulation of Mi14-3-3 family members.

Genome-Wide Identification and Characterization of the Abiotic-Stress-Responsive GRF Gene Family in Diploid Woodland Strawberry (Fragaria vesca)

Growth regulatory factors (GRF) are plant-specific transcription factors that play an important role in plant resistance to stress. This gene family in strawberry has not been investigated previously. In this study, 10 GRF genes were identified in the genome of the diploid woodland strawberry (Fragaria vesca). Chromosome analysis showed that the 10 FvGRF genes were unevenly distributed on five chromosomes. Phylogenetic analysis resolved the FvGRF proteins into five groups. Genes of similar structure were placed in the same group, which was indicative of functional redundance. Whole-genome duplication/segmental duplication and dispersed duplication events effectively promoted expansion of the strawberry GRF gene family. Quantitative reverse transcription-PCR analysis suggested that FvGRF genes played potential roles in the growth and development of vegetative organs. Expression profile analysis revealed that FvGRF3, FvGRF5, and FvGRF7 were up-regulated under low-temperature stress, FvGRF4 and FvGRF9 were up-regulated under high-temperature stress, FvGRF6 and FvGRF8 were up-regulated under drought stress, FvGRF3, FvGRF6, and FvGRF8 were up-regulated under salt stress, FvGRF2, FvGRF7, and FvGRF9 were up-regulated under salicylic acid treatment, and FvGRF3, FvGRF7, FvGRF9, and FvGRF10 were up-regulated under abscisic acid treatment. Promoter analysis indicated that FvGRF genes were involved in plant growth and development and stress response. These results provide a theoretical and empirical foundation for the elucidation of the mechanisms of abiotic stress responses in strawberry.

Genome-wide identification and expression analysis of the growth regulating factor (GRF) family in Jatropha curcas

GRF genes have been confirmed to have important regulatory functions in plant growth, development and response to abiotic stress. Although the genome of Jatropha curcas is sequenced, knowledge about the identification of the species' GRF genes and their expression patterns is still lacking. In this study, we characterized the 10 JcGRF genes. A detailed investigation into the physic nut GRF gene family is performed, including analysis of the exon-intron structure, conserved domains, conserved motifs, phylogeny, chromosomal locations, potential small RNA targets and expression profiles under both normal growth and abiotic stress conditions. Phylogenetic analysis indicated that the 10 JcGRF genes were classified into five groups corresponding to group I, II, III, IV and V. The analysis of conserved domains showed that the motifs of JcGRF genes were highly conserved in Jatropha curcas. Expression analysis based on RNA-seq and qRT-PCR showed that almost all JcGRF genes had the highest expression in seeds, but very low expression was detected in the non-seed tissues tested, and four JcGRF genes responded to at least one abiotic stress at at least one treatment point. Our research will provide an important scientific basis for further research on the potential functions of JcGRF genes in Jatropha curcas growth and development, and response to abiotic stress, and will eventually provide candidate genes for the breeding of Jatropha curcas.

Characterization of cotton ARF factors and the role of GhARF2b in fiber development

BACKGROUND

Cotton fiber is a model system for studying plant cell development. At present, the functions of many transcription factors in cotton fiber development have been elucidated, however, the roles of auxin response factor (ARF) genes in cotton fiber development need be further explored.

RESULTS

Here, we identify auxin response factor (ARF) genes in three cotton species: the tetraploid upland cotton G. hirsutum, which has 73 ARF genes, and its putative extent parental diploids G. arboreum and G. raimondii, which have 36 and 35 ARFs, respectively. Ka and Ks analyses revealed that in G. hirsutum ARF genes have undergone asymmetric evolution in the two subgenomes. The cotton ARFs can be classified into four phylogenetic clades and are actively expressed in young tissues. We demonstrate that GhARF2b, a homolog of the Arabidopsis AtARF2, was preferentially expressed in developing ovules and fibers. Overexpression of GhARF2b by a fiber specific promoter inhibited fiber cell elongation but promoted initiation and, conversely, its downregulation by RNAi resulted in fewer but longer fiber. We show that GhARF2b directly interacts with GhHOX3 and represses the transcriptional activity of GhHOX3 on target genes.

CONCLUSION

Our results uncover an important role of the ARF factor in modulating cotton fiber development at the early stage.