Identification and transcript profiles of citrus growth-regulating factor genes involved in the regulation of leaf and fruit development

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.

Association analyses reveal both anthropic and environmental selective events during eggplant domestication

Eggplant (Solanum melongena) is one of the four most important Solanaceous crops, widely cultivated and consumed in Asia, the Mediterranean basin, and Southeast Europe. We studied the genome-wide association of historical genebank phenotypic data on a genotyped worldwide collection of 3449 eggplant accessions. Overall, 334 significant associations for key agronomic traits were detected. Significant correlations were obtained between different types of phenotypic data, some of which were not obvious, such as between fruit size/yield and fruit color components, suggesting simultaneous anthropic selection for genetically unrelated traits. Anthropic selection of traits like leaf prickles, fruit color, and yield, acted on distinct genomic regions in the two domestication centers (India and Southeast Asia), further confirming the multiple domestication of eggplant. To discriminate anthropic from environmental selection in domestication centers, we conducted a genotype-environment association (GEA) on a subset of georeferenced accessions from the Indian subcontinent. The population structure in this area revealed four genetic clusters, corresponding to a latitudinal gradient, and environmental factors explained 31% of the population structure when the effect of spatial distances was removed. GEA and outlier association identified 305 candidate regions under environmental selection, containing genes for abiotic stress responses, plant development, and flowering transition. Finally, in the Indian domestication center anthropic and environmental selection acted largely independently, and on different genomic regions. These data allow a better understanding of the different effects of environmental and anthropic selection during domestication of a crop, and the different world regions where some traits were initially selected by humans.

The MADS-box transcription factor CsAGL9 plays essential roles in seed setting in Camellia sinensis

The number of seed setting (NSS) is an important biological trait that affects tea propagation and yield. In this study, the NSS of an F1 tea population (n = 324) generated via a cross between 'Longjing 43' and 'Baihaozao' was investigated at two locations in two consecutive years. Quantitative trait locus (QTL) mapping of the NSS was performed, and 10 major QTLs were identified. In total, 318 genes were found in these 10 QTLs intervals, and 11 key candidate genes were preliminarily identified. Among them, the MADS-box transcription factor AGAMOUS LIKE 9 (CsAGL9, CSS0037962) located in the most stable QTL (qNSS2) was identified as a key gene affecting the NSS. CsAGL9 overexpression in Arabidopsis promoted early flowering and significantly decreased the length and number of pods and number of seeds per pod. Transcriptome analysis demonstrated that the auxin pathway, a key hormone pathway regulating plant reproduction, was highly affected in the transgenic lines. The auxin pathway was likewise the most prominent in the gene co-expression network study of CsAGL9 in tea plants. In summary, we identified CsAGL9 is essential for seed setting using QTL mapping integrated with RNA-seq, which shed a new light on the mechanism NSS of seed setting in tea plants.

Overview of molecular mechanisms of plant leaf development: a systematic review

Leaf growth initiates in the peripheral region of the meristem at the apex of the stem, eventually forming flat structures. Leaves are pivotal organs in plants, serving as the primary sites for photosynthesis, respiration, and transpiration. Their development is intricately governed by complex regulatory networks. Leaf development encompasses five processes: the leaf primordium initiation, the leaf polarity establishment, leaf size expansion, shaping of leaf, and leaf senescence. The leaf primordia starts from the side of the growth cone at the apex of the stem. Under the precise regulation of a series of genes, the leaf primordia establishes adaxial-abaxial axes, proximal-distal axes and medio-lateral axes polarity, guides the primordia cells to divide and differentiate in a specific direction, and finally develops into leaves of a certain shape and size. Leaf senescence is a kind of programmed cell death that occurs in plants, and as it is the last stage of leaf development. Each of these processes is meticulously coordinated through the intricate interplay among transcriptional regulatory factors, microRNAs, and plant hormones. This review is dedicated to examining the regulatory influences of major regulatory factors and plant hormones on these five developmental aspects of leaves.

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.

Genome-Wide Identification and Characterization of the GRF Gene Family in Melastoma dodecandrum

Growth-regulating factor (GRF) is a kind of transcription factor unique to plants, playing an important role in the flowering regulation, growth, and development of plants. Melastoma dodecandrum is an important member of Melastomataceae, with ornamental, medicinal, and edible benefits. The identification of the GRF gene family in M. dodecandrum can help to improve their character of flavor and continuous flowering. The members of the GRF gene family were identified from the M. dodecandrum genome, and their bioinformatics, selective pressure, and expression patterns were analyzed. The results showed that there were 20 GRF genes in M. dodecandrum. Phylogenetic analysis showed that the 71 GRF genes from M. dodecandrum, Arabidopsis thaliana, Camellia sinensis, and Oryza sativa can be divided into three clades and six subclades. The 20 GRF genes of M. dodecandrum were distributed in twelve chromosomes and one contig. Furthermore, the gene structure and motif analysis showed that the intron and motif within each clade were very similar, but there were great differences among different clades. The promoter contained cis-acting elements related to hormone induction, stress, and growth and development. Different transcriptomic expression of MdGRFs indicated that MdGRFs may be involved in regulating the growth and development of M. dodecandrum. The results laid a foundation for further study on the function and molecular mechanism of the M. dodecandrum GRF gene family.

Comprehensive Analysis for GRF Transcription Factors in Sacred Lotus (Nelumbo nucifera)

Sacred lotus (Nelumbo nucifera) is an aquatic perennial plant with essential food, ornamental, and pharmacological value. Growth-regulating factor (GRF) is a transcription factor (TF) family that plays an important role in regulating the growth and development of plants. In this study, a comprehensive analysis of the GRF family in N. nucifera was performed, and its role in N. nucifera development was studied. A total of eight GRF genes were identified in the N. nucifera genome. Phylogenetic analysis divided the 38 GRF genes into six clades, while the NuGRFs only contained five clades. The analyses of gene structures, motifs, and cis-acting regulatory elements of the GRF gene family were performed. In addition, the chromosome location and collinearity were analyzed. The expression pattern based on transcriptomic data and real-time reverse transcription-quantitative PCR (qRT-PCR) revealed that the GRF genes were expressed in multiple organs and were abundant in actively growing tissues, and the expression levels decreased as the age of N. nucifera increased. Then, 3D structures of the NuGRF proteins were predicted by homology modeling. Finally, the subcellular localization of GRF1 was ascertained in the tobacco leaf through a vector. Therefore, this study provides a comprehensive overview of the GRF TF family in N. nucifera.

High efficient de novo root-to-shoot organogenesis in Citrus jambhiri Lush.: Gene expression, genetic stability and virus indexing

A protocol for high-frequency direct organogenesis from root explants of Kachai lemon (Citrus jambhiri Lush.) was developed. Full-length roots (~3 cm) were isolated from the in vitro grown seedlings and cultured on Murashige and Skoog basal medium supplemented with Nitsch vitamin (MSN) with different concentrations of cytokinin [6-benzylaminopurine, (BAP)] and gibberellic acid (GA₃). The frequency of multiple shoot proliferation was very high, with an average of 34.3 shoots per root explant when inoculated on the MSN medium supplemented with BAP (1.0 mg L^–1) and GA₃ (1.0 mg L^–1). Optimal rooting was induced in the plantlets under half strength MSN medium supplemented with indole-3-acetic acid (IAA, 0.5–1.0 mg L^–1). IAA induced better root structure than 1-naphthaleneacetic acid (NAA), which was evident from the scanning electron microscopy (SEM). The expressions of growth regulating factor genes (GRF1 and GRF5) and GA₃ signaling genes (GA2OX1 and KO1) were elevated in the regenerants obtained from MSN+BAP (1.0 mg L^-1)+GA₃ (1.0 mg L^-1). The expressions of auxin regulating genes were high in roots obtained in ½ MSN+IAA 1.0 mg L^-1. Furthermore, indexing of the regenerants confirmed that there was no amplicons detected for Huanglongbing bacterium and Citrus tristeza virus. Random amplified polymorphic DNA (RAPD) and inter simple sequence repeat (ISSR) markers detected no polymorphic bands amongst the regenerated plants. This is the first report that describes direct organogenesis from the root explant of Citrus jambhiri Lush. The high-frequency direct regeneration protocol in the present study provides an enormous significance in Citrus organogenesis, its commercial cultivation and genetic conservation.

Genome-wide analysis of growth-regulating factors (GRFs) in Triticum aestivum

The Growth-Regulating Factor (GRF) family encodes a type of plant-specific transcription factor (TF). GRF members play vital roles in plant development and stress response. Although GRF family genes have been investigated in a variety of plants, they remain largely unstudied in bread wheat (Triticum aestivum L.). The present study was conducted to comprehensively identify and characterize the T. aestivum GRF (TaGRF) gene family members. We identified 30 TaGRF genes, which were divided into four groups based on phylogenetic relationship. TaGRF members within the same subgroup shared similar motif composition and gene structure. Synteny analysis suggested that duplication was the dominant reason for family member expansion. Expression pattern profiling showed that most TaGRF genes were highly expressed in growing tissues, including shoot tip meristems, stigmas and ovaries, suggesting their key roles in wheat growth and development. Further qRT-PCR analysis revealed that all 14 tested TaGRFs were significantly differentially expressed in responding to drought or salt stresses, implying their additional involvement in stress tolerance of wheat. Our research lays a foundation for functional determination of TaGRFs, and will help to promote further scrutiny of their regulatory network in wheat development and stress response.

Characterization of poplar growth-regulating factors and analysis of their function in leaf size control

BACKGROUND

Growth-regulating factors (GRFs) are plant-specific transcription factors that control organ size. Nineteen GRF genes were identified in the Populus trichocarpa genome and one was reported to control leaf size mainly by regulating cell expansion. In this study, we further characterize the roles of the other poplar GRFs in leaf size control in a similar manner.

RESULTS

The 19 poplar GRF genes were clustered into six groups according to their phylogenetic relationship with Arabidopsis GRFs. Bioinformatic analysis, degradome, and transient transcription assays showed that 18 poplar GRFs were regulated by miR396, with GRF12b the only exception. The functions of PagGRF6b (Pag, Populus alba × P. glandulosa), PagGRF7a, PagGRF12a, and PagGRF12b, representing three different groups, were investigated. The results show that PagGRF6b may have no function on leaf size control, while PagGRF7a functions as a negative regulator of leaf size by regulating cell expansion. By contrast, PagGRF12a and PagGRF12b may function as positive regulators of leaf size control by regulating both cell proliferation and expansion, primarily cell proliferation.

CONCLUSIONS

The diversity of poplar GRFs in leaf size control may facilitate the specific, coordinated regulation of poplar leaf development through fine adjustment of cell proliferation and expansion.

Genome-wide identification of GRF transcription factors in soybean and expression analysis of GmGRF family under shade stress

BACKGROUND

The Growth-regulating factor (GRF) family encodes plant-specific transcription factors which contain two conserved domains, QLQ and WRC. Members of this family play vital roles in plant development and stress response processes. Although GRFs have been identified in various plant species, we still know little about the GRF family in soybean (Glycine max).

RESULTS

In the present study, 22 GmGRFs distributed on 14 chromosomes and one scaffold were identified by searching soybean genome database and were clustered into five subgroups according to their phylogenetic relationships. GmGRFs belonging to the same subgroup shared a similar motif composition and gene structure. Synteny analysis revealed that large-scale duplications played key roles in the expansion of the GmGRF family. Tissue-specific expression data showed that GmGRFs were strongly expressed in growing tissues, including the shoot apical meristems, developing seeds and flowers, indicating that GmGRFs play critical roles in plant growth and development. On the basis of expression analysis of GmGRFs under shade conditions, we found that all GmGRFs responded to shade stress. Most GmGRFs were down-regulated in soybean leaves after shade treatment.

CONCLUSIONS

Taken together, this research systematically analyzed the characterization of the GmGRF family and its primary roles in soybean development and shade stress response. Further studies of the function of the GmGRFs in the growth, development and stress tolerance of soybean, especially under shade stress, will be valuable.

Global Methylome and gene expression analysis during early Peanut pod development

BACKGROUND

Early peanut pod development is an important process of peanut reproductive development. Modes of DNA methylation during early peanut pod development are still unclear, possibly because its allotetraploid genome may cause difficulty for the methylome analysis.

RESULTS

To investigate the functions of the dynamic DNA methylation during the early development of the peanut pod, global methylome and gene expression analyses were carried out by Illumina high throughput sequencing. A novel mapping strategy of reads was developed and used for methylome and gene expression analysis. Differentially methylated genes, such as nodulin, cell number regulator-like protein, and senescence-associated genes, were identified during the early developmental stages of the peanut pod. The expression levels of gibberellin-related genes changed during this period of pod development. From the stage one (S1) gynophore to the stage two (S2) gynophore, the expression levels of two key methyltransferase genes, DRM2 and MET1, were up-regulated, which may lead to global DNA methylation changes between these two stages. The differentially methylated and expressed genes identified in the S1, S2, and stage 3 (S3) gynophore are involved in different biological processes such as stem cell fate determination, response to red, blue, and UV light, post-embryonic morphogenesis, and auxin biosynthesis. The expression levels of many genes were co-related by their DNA methylation levels. In addition, our results showed that the abundance of some 24-nucleotide siRNAs and miRNAs were positively associated with DNA methylation levels of their target loci in peanut pods.

CONCLUSION

A novel mapping strategy of reads was described and verified in this study. Our results suggest that the methylated modes of the S1, S2, and S3 gynophore are different. The methylation changes that were identified during early peanut pod development provide useful information for understanding the roles of epigenetic regulation in peanut pod development.

Comparative Analysis of Fruit Ripening-Related miRNAs and Their Targets in Blueberry Using Small RNA and Degradome Sequencing

MicroRNAs (miRNAs) play vital roles in the regulation of fruit development and ripening. Blueberry is an important small berry fruit crop with economical and nutritional value. However, nothing is known about the miRNAs and their targets involved in blueberry fruit ripening. In this study, using high-throughput sequencing of small RNAs, 84 known miRNAs belonging to 28 families and 16 novel miRNAs were identified in white fruit (WF) and blue fruit (BF) libraries, which represent fruit ripening onset and in progress, respectively. Among them, 41 miRNAs were shown to be differentially expressed during fruit maturation, and 16 miRNAs representing 16 families were further chosen to validate the sRNA sequencing data by stem-loop qRT-PCR. Meanwhile, 178 targets were identified for 41 known and 7 novel miRNAs in WF and BF libraries using degradome sequencing, and targets of miR160 were validated using RLM-RACE (RNA Ligase-Mediated (RLM)-Rapid Amplification of cDNA Ends) approach. Moreover, the expression patterns of 6 miRNAs and their targets were examined during fruit development and ripening. Finally, integrative analysis of miRNAs and their targets revealed a complex miRNA-mRNA regulatory network involving a wide variety of biological processes. The findings will facilitate future investigations of the miRNA-mediated mechanisms that regulate fruit development and ripening in blueberry.

Genome-wide identification of leaf abscission associated microRNAs in sugarcane (Saccharum officinarum L.)

Background

Sugarcane (Saccharum officinarum L.) is an economically important crop, mainly due to the production of sugar and biofuel (Azevedo RA, Carvalho RF, Cia MC, & Gratão PL, Trop Plant Biol 4:42-51, 2011). Grown mainly in tropical and subtropical countries, sugarcane is a highly polyploid plant with up to ten copies of each chromosome, which increases the difficulties of genome assembly and genetic, physiologic and biochemical analyses. The increasing demands of sugar and the increasing cost of sugarcane harvest require sugarcane varieties which can shed their leaves during the maturity time, so it is important to study the mechanism of leaf abscission in sugarcane.

Results

To improve the understanding of miRNA roles in sugarcane leaf abscission, we reported the genome-wide characterization of miRNAs and their putative targets in sugarcane using deep sequencing for six small RNA libraries. In total, 93 conserved miRNAs and 454 novel miRNAs were identified in sugarcane using previously reported transcriptome as reference. Among them, 25 up-regulated and 13 down-regulated miRNAs were identified in leaf abscission sugarcane plants (LASP) compared to leaf packaging sugarcane plants (LPSP). Target prediction revealed several miRNA-mRNA modules including miR156-SPL, miR319-TPR2, miR396-GRF and miR408-LAC3 might be involved in the sugarcane leaf abscission. KEGG pathway enrichment analysis showed differentially expressed miRNAs may regulate pathways like “plant hormone signal transduction” and “plant-pathogen interaction”, which is consistent with previous transcriptome study. In addition, we identified 96 variant miRNAs with 135 single nucleotide polymorphisms (SNPs). The expression of sugarcane miRNAs and variant miRNAs were confirmed by qRT-PCR. We identified a possible poaceae specific miRNA called miR5384 for the first time in sugarcane.

Conclusions

We not only reported miR5384, a possible poaceae specific miRNA, for the first time in sugarcane but also presented some miRNA-mRNA modules including miR156-SPL, miR319-TPR2, miR396-GRF and miR408-LAC in sugarcane. These modules might be involved in the regulation of sugarcane leaf abscission during the maturity time. All of these findings may lay ground work for future application of sugarcane breeding program and benefit research studies of sugarcane miRNAs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4053-3) contains supplementary material, which is available to authorized users.

Molecular Characterization and Expression Profiling of Tomato GRF Transcription Factor Family Genes in Response to Abiotic Stresses and Phytohormones

Growth regulating factors (GRFs) are plant-specific transcription factors that are involved in diverse biological and physiological processes, such as growth, development and stress and hormone responses. However, the roles of GRFs in vegetative and reproductive growth, development and stress responses in tomato (Solanum lycopersicum) have not been extensively explored. In this study, we characterized the 13 SlGRF genes. In silico analysis of protein motif organization, intron–exon distribution, and phylogenetic classification confirmed the presence of GRF proteins in tomato. The tissue-specific expression analysis revealed that most of the SlGRF genes were preferentially expressed in young and growing tissues such as flower buds and meristems, suggesting that SlGRFs are important during growth and development of these tissues. Some of the SlGRF genes were preferentially expressed in fruits at distinct developmental stages suggesting their involvement in fruit development and the ripening process. The strong and differential expression of different SlGRFs under NaCl, drought, heat, cold, abscisic acid (ABA), and jasmonic acid (JA) treatment, predict possible functions for these genes in stress responses in addition to their growth regulatory functions. Further, differential expression of SlGRF genes upon gibberellic acid (GA3) treatment indicates their probable function in flower development and stress responses through a gibberellic acid (GA)-mediated pathway. The results of this study provide a basis for further functional analysis and characterization of this important gene family in tomato.

Systems biology study of transcriptional and post-transcriptional co-regulatory network sheds light on key regulators involved in important biological processes in Citrus sinensis

Transcriptional and post-transcriptional regulators including transcription regulator, transcription factor and miRNA are the main endogenous molecular elements which control complex cellular mechanisms such as development, growth and response to biotic and abiotic stresses in a coordinated manner in plants. Utilizing the most recent information on such relationships in a plant species, obtained from high-throughput experimental technologies and advanced computational tools, we can reconstruct its co-regulatory network which consequently sheds light on key regulators involved in its important biological processes. In this article, combined systems biology approaches such as mining the literatures, various databases and network reconstruction, analysis, and visualization tools were employed to infer and visualize the coregulatory relationships between miRNAs and transcriptional regulators in Citrus sinensis. Using computationally and experimentally verified miRNA-target interactions and constructed co-expression networks on array-based data, 10 coregulatory networks and 10 corresponding subgraphs include FFL motifs were obtained for 10 distinct tissues/conditions. Then PPI subnetworks were extracted for transcripts/genes included in mentioned subgraphs in order to the functional analysis of extracted coregulatory circuits. These proposed coregulatory connections shed light on precisely identifying C. sinensis metabolic pathways key switches, which are demanded for ultimate goals such as genome editing.