Evolution and divergence of SBP-box genes in land plants

Comparative analysis of SPL transcription factors from streptophyte algae and embryophytes reveals evolutionary trajectories of SPL family in streptophytes

SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) genes encode plant-specific transcription factors which are important regulators of diverse plant developmental processes. We took advantage of available genome sequences of streptophyte algae representatives to investigate the relationships of SPL genes between freshwater green algae and land plants. Our analysis showed that streptophyte algae, hornwort and liverwort genomes encode from one to four SPL genes which is the smallest set, in comparison to other land plants studied to date. Based on the phylogenetic analysis, four major SPL phylogenetic groups were distinguished with Group 3 and 4 being sister to Group 1 and 2. Comparative motif analysis revealed conserved protein motifs within each phylogenetic group and unique bryophyte-specific motifs within Group 1 which suggests lineage-specific protein speciation processes. Moreover, the gene structure analysis also indicated the specificity of each by identifying differences in exon-intron structures between the phylogenetic groups, suggesting their evolutionary divergence. Since current understanding of SPL genes mostly arises from seed plants, the presented comparative and phylogenetic analyzes from freshwater green algae and land plants provide new insights on the evolutionary trajectories of the SPL gene family in different classes of streptophytes.

Temporal regulation of vegetative phase change in plants

During their vegetative growth, plants reiteratively produce leaves, buds, and internodes at the apical end of the shoot. The identity of these organs changes as the shoot develops. Some traits change gradually, but others change in a coordinated fashion, allowing shoot development to be divided into discrete juvenile and adult phases. The transition between these phases is called vegetative phase change. Historically, vegetative phase change has been studied because it is thought to be associated with an increase in reproductive competence. However, this is not true for all species; indeed, heterochronic variation in the timing of vegetative phase change and flowering has made important contributions to plant evolution. In this review, we describe the molecular mechanism of vegetative phase change, how the timing of this process is controlled by endogenous and environmental factors, and its ecological and evolutionary significance.

Genome-wide identification, expression analysis, and functional study of the bZIP transcription factor family and its response to hormone treatments in pea (Pisum sativum L.)

BACKGROUND

Basic leucine zipper (bZIP) protein is a plant-specific transcription factor involved in various biological processes, including light signaling, seed maturation, flower development, cell elongation, seed accumulation protein, and abiotic and biological stress responses. However, little is known about the pea bZIP family.

RESULTS

In this study, we identified 87 bZIP genes in pea, named PsbZIP1 ~ PsbZIP87, via homology analysis using Arabidopsis. The genes were divided into 12 subfamilies and distributed unevenly in 7 pea chromosomes. PsbZIPs in the same subfamily contained similar intron/exon organization and motif composition. 1 tandem repeat event and 12 segmental duplication events regulated the expansion of the PsbZIP gene family. To better understand the evolution of the PsbZIP gene family, we conducted collinearity analysis using Arabidopsis thaliana, Oryza sativa Japonica, Fagopyrum tataricum, Solanum lycopersicum, Vitis vinifera, and Brachypodium distachyon as the related species of pea. In addition, interactions between PsbZIP proteins and promoters containing hormone- and stress-responsive cis-acting elements suggest that the regulation of PsbZIP expression was complex. We also evaluated the expression patterns of bZIP genes in different tissues and at different fruit development stages, all while subjecting them to five hormonal treatments.

CONCLUSION

These results provide a deeper understanding of PsbZIP gene family evolution and resources for the molecular breeding of pea. The findings suggested that PsbZIP genes, specifically PSbZIP49, play key roles in the development of peas and their response to various hormones.

Genome-wide identification, phylogenetic and expression pattern analysis of HSF family genes in the Rye (Secale cereale L.)

BACKGROUND

Heat shock factor (HSF), a typical class of transcription factors in plants, has played an essential role in plant growth and developmental stages, signal transduction, and response to biotic and abiotic stresses. The HSF genes families has been identified and characterized in many species through leveraging whole genome sequencing (WGS). However, the identification and systematic analysis of HSF family genes in Rye is limited.

RESULTS

In this study, 31 HSF genes were identified in Rye, which were unevenly distributed on seven chromosomes. Based on the homology of A. thaliana, we analyzed the number of conserved domains and gene structures of ScHSF genes that were classified into seven subfamilies. To better understand the developmental mechanisms of ScHSF family during evolution, we selected one monocotyledon (Arabidopsis thaliana) and five (Triticum aestivum L., Hordeum vulgare L., Oryza sativa L., Zea mays L., and Aegilops tauschii Coss.) specific representative dicotyledons associated with Rye for comparative homology mapping. The results showed that fragment replication events modulated the expansion of the ScHSF genes family. In addition, interactions between ScHSF proteins and promoters containing hormone- and stress-responsive cis-acting elements suggest that the regulation of ScHSF expression was complex. A total of 15 representative genes were targeted from seven subfamilies to characterize their gene expression responses in different tissues, fruit developmental stages, three hormones, and six different abiotic stresses.

CONCLUSIONS

This study demonstrated that ScHSF genes, especially ScHSF1 and ScHSF3, played a key role in Rye development and its response to various hormones and abiotic stresses. These results provided new insights into the evolution of HSF genes in Rye, which could help the success of molecular breeding in Rye.

The SPL transcription factor genes are potential targets for epigenetic regulation in response to drought stress in chickpea (C. arietinum L.)

BACKGROUND

Crop improvement for tolerance to various biotic and abiotic stress factors necessitates understanding the key gene regulatory mechanisms. One such mechanism of gene regulation involves changes in cytosine methylation at the gene body and flanking regulatory sequences. The present study was undertaken to identify genes which might be potential targets of drought-induced DNA methylation in chickpea.

METHODS AND RESULTS

Two chickpea genotypes, which contrast for drought tolerance, were subjected to drought stress conditions and their differential response was studied by analysing different morpho-physiological traits. Utilizing the in-house, high throughput sequencing data, the SQUAMOSA promoter-binding (SBP) protein-like (SPL) transcription factor genes were identified to be differentially methylated and expressed amongst the two genotypes, in response to drought stress. The methylation status of one of these genes was examined and validated through bisulfite PCR (BS-PCR). The identified genes could be possible homologs to known epialleles and can therefore serve as potential epialleles which can be utilized for crop improvement in chickpea.

CONCLUSION

The SPL TF genes are potential targets of epigenetic regulation in response to drought stress in chickpea. Since these are TFs, they might play important roles in controlling the expression of other genes, thus contributing to differential drought response of the two genotypes.

Genome-wide identification and expression analysis of the SPL transcription factor family and its response to abiotic stress in Quinoa (Chenopodium quinoa)

BACKGROUND

Squamous promoter binding protein-like (SPL) proteins are a class of transcription factors that play essential roles in plant growth and development, signal transduction, and responses to biotic and abiotic stresses. The rapid development of whole genome sequencing has enabled the identification and characterization of SPL gene families in many plant species, but to date this has not been performed in quinoa (Chenopodium quinoa).

RESULTS

This study identified 23 SPL genes in quinoa, which were unevenly distributed on 18 quinoa chromosomes. Quinoa SPL genes were then classified into eight subfamilies based on homology to Arabidopsis thaliana SPL genes. We selected three dicotyledonous and monocotyledonous representative species, each associated with C. quinoa, for comparative sympatric mapping to better understand the evolution of the developmental mechanisms of the CqSPL family. Furthermore, we also used 15 representative genes from eight subfamilies to characterize CqSPLs gene expression in different tissues and at different fruit developmental stages under six different abiotic stress conditions.

CONCLUSIONS

This study, the first to identify and characterize SPL genes in quinoa, reported that CqSPL genes, especially CqSPL1, play a critical role in quinoa development and in its response to various abiotic stresses.

Genome-Wide Identification and Characterization of SPL Family Genes in Chenopodium quinoa

SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) genes encode a large family of plant-specific transcription factors that play important roles in plant growth, development, and stress responses. However, there is little information available on SPL genes in Chenopodiaceae. Here, 23 SPL genes were identified and characterized in the highly nutritious crop Chenopodium quinoa. Chromosome localization analysis indicated that the 23 CqSPL genes were unevenly distributed on 12 of 18 chromosomes. Two zinc finger-like structures and a nuclear location signal were present in the SBP domains of all CqSPLs, with the exception of CqSPL21/22. Phylogenetic analysis revealed that these genes were classified into eight groups (group I–VIII). The exon–intron structure and motif composition of the genes in each group were similar. Of the 23 CqSPLs, 13 were potential targets of miR156/7. In addition, 5 putative miR156-encoding loci and 13 putative miR157-encoding loci were predicted in the quinoa genome, and they were unevenly distributed on chromosome 1–4. The expression of several Cqu-MIR156/7 loci was confirmed by reverse transcription polymerase chain reaction in seedlings. Many putative cis-elements associated with light, stress, and phytohormone responses were identified in the promoter regions of CqSPLs, suggesting that CqSPL genes are likely involved in the regulation of key developmental processes and stress responses. Expression analysis revealed highly diverse expression patterns of CqSPLs among tissues. Many CqSPLs were highly expressed in leaves, flowers, and seeds, and their expression levels were low in the roots, suggesting that CqSPLs play distinct roles in the development and growth of quinoa. The expression of 13 of 23 CqSPL genes responded to salt treatment (11 up-regulated and 2 down-regulated). A total of 22 of 23 CqSPL genes responded to drought stress (21 up-regulated and 1 down-regulated). Moreover, the expression of 14 CqSPL genes was significantly altered following cadmium treatment (3 up-regulated and 11 down-regulated). CqSPL genes are thus involved in quinoa responses to salt/drought and cadmium stresses. These findings provide new insights that will aid future studies of the biological functions of CqSPLs in C. quinoa.

A genome-wide comparative evolutionary analysis of zinc finger-BED transcription factor genes in land plants

Zinc finger (Zf)-BED proteins are a novel superfamily of transcription factors that controls numerous activities in plants including growth, development, and cellular responses to biotic and abiotic stresses. Despite their important roles in gene regulation, little is known about the specific functions of Zf-BEDs in land plants. The current study identified a total of 750 Zf-BED-encoding genes in 35 land plant species including mosses, bryophytes, lycophytes, gymnosperms, and angiosperms. The gene family size was somewhat proportional to genome size. All identified genes were categorized into 22 classes based on their specific domain architectures. Of these, class I (Zf-BED_DUF-domain_Dimer_Tnp_hAT) was the most common in the majority of the land plants. However, some classes were family-specific, while the others were species-specific, demonstrating diversity at different classification levels. In addition, several novel functional domains were also predicated including WRKY and nucleotide-binding site (NBS). Comparative genomics, transcriptomics, and proteomics provided insights into the evolutionary history, duplication, divergence, gene gain and loss, species relationship, expression profiling, and structural diversity of Zf-BEDs in land plants. The comprehensive study of Zf-BEDs in Gossypium sp., (cotton) also demonstrated a clear footprint of polyploidization. Overall, this comprehensive evolutionary study of Zf-BEDs in land plants highlighted significant diversity among plant species.

The miR157-SPL-CNR module acts upstream of bHLH101 to negatively regulate iron deficiency responses in tomato

Iron (Fe) homeostasis is critical for plant growth, development, and stress responses. Fe levels are tightly controlled by intricate regulatory networks in which transcription factors (TFs) play a central role. A series of basic helix-loop-helix (bHLH) TFs have been shown to contribute to Fe homeostasis, but the regulatory layers beyond bHLH TFs remain largely unclear. Here, we demonstrate that the SQUAMOSA PROMOTER-BINDING PROTEIN-LIKE (SPL) TF SlSPL-CNR negatively regulates Fe-deficiency responses in tomato (Solanum lycopersicum) roots. Fe deficiency rapidly repressed the expression of SlSPL-CNR, and Fe deficiency responses were intensified in two clustered regularly interspaced palindromic repeats (CRISPR)/CRISPR-associated protein 9-generated SlSPL-CNR knock-out lines compared to the wild-type. Comparative transcriptome analysis identified 47 Fe deficiency-responsive genes the expression of which is negatively regulated by SlSPL-CNR, one of which, SlbHLH101, helps regulate Fe uptake genes. SlSPL-CNR localizes the nucleus and interacts with the GTAC and BOX 4 (ATTAAT) motifs in the SlbHLH101 promoter to repress its expression. Inhibition of SlSPL-CNR expression in response to Fe deficiency was well correlated with the expression of the microRNA SlymiR157. SlymiR157-overexpressing tomato lines displayed enhanced Fe deficiency responses, as did SlSPL-CNR loss-of-function mutants. We propose that the SlymiR157-SlSPL-CNR module represents a novel pathway that acts upstream of SlbHLH101 to regulate Fe homeostasis in tomato roots.

The Mulberry SPL Gene Family and the Response of MnSPL7 to Silkworm Herbivory through Activating the Transcription of MnTT2L2 in the Catechin Biosynthesis Pathway

SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) genes, as unique plant transcription factors, play important roles in plant developmental regulation and stress response adaptation. Although mulberry is a commercially valuable tree species, there have been few systematic studies on SPL genes. In this work, we identified 15 full-length SPL genes in the mulberry genome, which were distributed on 4 Morus notabilis chromosomes. Phylogenetic analysis clustered the SPL genes from five plants (Malus × domestica Borkh, Populus trichocarpa, M. notabilis, Arabidopsis thaliana, and Oryza sativa) into five groups. Two zinc fingers (Zn1 and Zn2) were found in the conserved SBP domain in all of the MnSPLs. Comparative analyses of gene structures and conserved motifs revealed the conservation of MnSPLs within a group, whereas there were significant structure differences among groups. Gene quantitative analysis showed that the expression of MnSPLs had tissue specificity, and MnSPLs had much higher expression levels in older mulberry leaves. Furthermore, transcriptome data showed that the expression levels of MnSPL7 and MnSPL14 were significantly increased under silkworm herbivory. Molecular experiments revealed that MnSPL7 responded to herbivory treatment through promoting the transcription of MnTT2L2 and further upregulating the expression levels of catechin synthesis genes (F3′H, DFR, and LAR).

Genome-Wide Identification and Functional Exploration of SBP-Box Gene Family in Black Pepper (Piper nigrum L.)

Black pepper (Piper nigrum L.), is dubbed "the King of Spices". However, the lack of genic knowledge has limited the understanding of its physiological processes and hindered the development of its molecular breeding. The SBP-box gene family is an important family in plant development and integrates multiple physiological processes. Here, we made a genome-wide identification of the pepper SBP-box gene family to provide evolutionary and functional information about this conserved transcription factor. In total, 34 SBP genes were identified in pepper. All these pepper SBP genes were clustered into eight groups, and one pepper group was not found in Arabidopsis thaliana. Segment duplications played the most important role in the expansion process of pepper SBP genes, and all these duplications were subjected to purifying selection. Half of pepper SBP genes were found miR156 target sites, and 17 miR156s were predicted. The tissue expression analysis revealed the differential expression of pepper SBP genes. Eleven SBP genes were found in four co-expression networks, and the GO enrichment further provides a functional prediction for pepper SBP genes. This study lays a foundation for further studies of pepper and provides a valuable reference for functional mining of pepper SBP genes.

Genomic analysis of SBP gene family in Saccharum spontaneum reveals their association with vegetative and reproductive development

Background

SQUAMOSA promoter binding proteins (SBPs) genes encode a family of plant-specific transcription factors involved in various growth and development processes, including flower and fruit development, leaf initiation, phase transition, and embryonic development. The SBP gene family has been identified and characterized in many species, but no systematic analysis of the SBP gene family has been carried out in sugarcane.

Results

In the present study, a total of 50 sequences for 30 SBP genes were identified by the genome-wide analysis and designated SsSBP1 to SsSBP30 based on their chromosomal distribution. According to the phylogenetic tree, gene structure and motif features, the SsSBP genes were classified into eight groups (I to VIII). By synteny analysis, 27 homologous gene pairs existed in SsSBP genes, and 37 orthologous gene pairs between sugarcane and sorghum were found. Expression analysis in different tissues, including vegetative and reproductive organs, showed differential expression patterns of SsSBP genes, indicating their functional diversity in the various developmental processes. Additionally, 22 SsSBP genes were predicted as the potential targets of miR156. The differential expression pattern of miR156 exhibited a negative correlation of transcription levels between miR156 and the SsSBP gene in different tissues.

Conclusions

The sugarcane genome possesses 30 SsSBP genes, and they shared similar gene structures and motif features in their subfamily. Based on the transcriptional and qRT-PCR analysis, most SsSBP genes were found to regulate the leaf initial and female reproductive development. The present study comprehensively and systematically analyzed SBP genes in sugarcane and provided a foundation for further studies on the functional characteristics of SsSBP genes during different development processes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08090-3.

The expansion and diversity of the CYP75 gene family in Vitaceae

The CYP75 gene family plays an important role in flavonoid biosynthesis in plants. Little is known about the evolution of the gene family within the grape family. Here, we extracted the CYP75 genes from transcriptome data of 15 grape species and 36 representative genomes from other plants to explore the evolutionary history of the CYP75 gene family in Vitaceae. The structure of the CYP75 protein sequences is highly conserved with the variation mainly occurring in the N terminal and the middle region. The evolutionary analyses suggested classifying the CYP75 gene family into three groups in Vitaceae, namely Vitaceae A1, Vitaceae A2 and Vitaceae B. The Vitaceae A1 and A2 belong to the CYP75A subfamily and the Vitaceae B belongs to the CYP75B subfamily. Within the Vitaceae A1, most Vitaceae taxa present only one copy of the CYP75A protein sequence except for Vitis vinifera with a high number of sequences, which might have originated through recent gene duplications after its split from the other species. Vitaceae A2 contain only CYP75A sequences from Vitaceae sister to one from Camellia sinensis, probably representing a relict lineage. The CYP75B proteins were found to be dominated in Vitaceae and other angiosperms. Our results provide important insights into understanding the evolutionary history of the CYP75 gene family in Vitaceae and other angiosperms.

Comprehensive Analysis of the SBP Family in Blueberry and Their Regulatory Mechanism Controlling Chlorophyll Accumulation

SQUAMOSA Promoter Binding Protein (SBP) family genes act as central players to regulate plant growth and development with functional redundancy and specificity. Addressing the diversity of the SBP family in crops is of great significance to precisely utilize them to improve agronomic traits. Blueberry is an important economic berry crop. However, the SBP family has not been described in blueberry. In the present study, twenty VcSBP genes were identified through data mining against blueberry transcriptome databases. These VcSBPs could be clustered into eight groups, and the gene structures and motif compositions are divergent among the groups and similar within each group. The VcSBPs were differentially expressed in various tissues. Intriguingly, 10 VcSBPs were highly expressed at green fruit stages and dramatically decreased at the onset of fruit ripening, implying that they are important regulators during early fruit development. Computational analysis showed that 10 VcSBPs were targeted by miR156, and four of them were further verified by degradome sequencing. Moreover, their functional diversity was studied in Arabidopsis. Noticeably, three VcSBPs significantly increased chlorophyll accumulation, and qRT-PCR analysis indicated that VcSBP13a in Arabidopsis enhanced the expression of chlorophyll biosynthetic genes such as AtDVR, AtPORA, AtPORB, AtPORC, and AtCAO. Finally, the targets of VcSBPs were computationally identified in blueberry, and the Y1H assay showed that VcSBP13a could physically bind to the promoter region of the chlorophyll-associated gene VcLHCB1. Our findings provided an overall framework for individually understanding the characteristics and functions of the SBP family in blueberry.

Genome-wide identification, phylogenetic analysis, and expression analysis of the SPL gene family in orchardgrass (Dactylis glomerata L.)

SPL (SQUAMOSA promoter binding protein-like) is a plant-specific transcription factor family that contains the conserved SBP domain, which plays a vital role in the vegetative-to-reproductive phase transition, flowering development and regulation, tillering/branching, and stress responses. Although the SPL family has been identified and characterized in various plant species, limited information about it has been obtained in orchardgrass, which is a critical forage crop worldwide. In this study, 17 putative DgSPL genes were identified among seven chromosomes, and seven groups that share similar gene structures and conserved motifs were determined by phylogenetic analysis. Of these, eight genes have potential target sites for miR156. cis-Element and gene ontology annotation analysis indicated DgSPLs may be involved in regulating development and abiotic stress responses. The expression patterns of eight DgSPL genes at five developmental stages, in five tissues, and under three stress conditions were determined by RNA-seq and qRT-PCR. These assays indicated DgSPLs are involved in vegetative-to-reproductive phase transition, floral development, and stress responses. The transient expression analysis in tobacco and heterologous expression assays in yeast indicated that miR156-targeted DG1G01828.1 and DG0G01071.1 are nucleus-localized proteins, that may respond to drought, salt, and heat stress. Our study represents the first systematic analysis of the SPL family in orchardgrass. This research provides a comprehensive assessment of the DgSPL family, which lays the foundation for further examination of the role of miR156/DgSPL in regulating development and stress responses in forages grasses.

Comparative analysis of the SPL gene family in five Rosaceae species: Fragaria vesca, Malus domestica, Prunus persica, Rubus occidentalis, and Pyrus pyrifolia

SQUAMOSA promoter-binding protein-like (SPL) transcription factors are very important for the plant growth and development. Here 15 RoSPLs were identified in Rubus occidentalis. The conserved domains and motifs, phylogenetic relationships, posttranscriptional regulation, and physiological function of the 92 SPL family genes in Fragaria vesca, Malus domestica, Prunus persica, R. occidentalis, and Pyrus pyrifolia were analyzed. Sequence alignment and phylogenetic analysis showed the SPL proteins had sequence conservation, some FvSPLs could be lost or developed, and there was a closer relationship between M. domestica and P. pyrifolia, F. vesca and R. occidentalis, respectively. Genes with similar motifs clustering together in the same group had their functional redundancy. Based on the function of SPLs in Arabidopsis thaliana, these SPLs could be involved in vegetative transition from juvenile to adult, morphological change in the reproductive phase, anthocyanin biosynthesis, and defense stress. Forty-eight SPLs had complementary sequences of miR156, of which nine PrpSPLs in P. persica and eight RoSPLs in R. occidentalis as the potential targets of miR156 were reported for the first time, suggesting the conservative regulatory effects of miR156 and indicating the roles of miR156-SPL modules in plant growth, development, and defense response. It provides a basic understanding of SPLs in Rosaceae plants.

CRISPR-based assessment of genomic structure in the conserved SQUAMOSA promoter-binding-like gene clusters in rice

Although SQUAMOSA promoter-binding-like (SPL) transcription factors are important regulators of development in rice (Oryza sativa), prior assessments of the SPL family have been limited to single genes. A functional comparison across the full gene family in standardized genetic backgrounds has not been reported previously. Here, we demonstrate that the SPL gene family in rice is enriched due to the most recent whole genome duplication (WGD). Notably, 10 of 19 rice SPL genes (52%) cluster in four units that have persisted for at least 50 million years. We show that SPL gene grouping and retention following WGD is widespread in angiosperms, suggesting the conservatism and importance of this gene arrangement. We used Cas9 editing to generate transformation lines for all 19 SPL genes in a common set of backgrounds, and found that knockouts of 14 SPL genes exhibited defects in plant height, 10 exhibited defects in panicle size, and nine had altered grain lengths. We observed subfunctionalization of genes in the paleoduplicated pairs, but little evidence of neofunctionalization. Expression of OsSPL3 was negatively correlated with that of its closest neighbor in its synteny group, OsSPL4, and its sister paired gene, OsSPL12, in the opposing group. Nucleotide diversity was lower in eight of the nine singleton genes in domesticated rice, relative to wild rice, whereas the reverse was true for the paired genes. Together, these results provide functional information on eight previously unexamined OsSPL family members and suggest that paleoduplicate pair redundancy benefits plant survival and innovation.

Global analysis of SBP gene family in Brachypodium distachyon reveals its association with spike development

SQUAMOSA-promoter binding like proteins (SBPs/SPLs) are plant specific transcription factors targeted by miR156 and involved in various biological pathways, playing multi-faceted developmental roles. This gene family is not well characterized in Brachypodium. We identified a total of 18 SBP genes in B. distachyon genome. Phylogenetic analysis revealed that SBP gene family in Brachypodium expanded through large scale duplication. A total of 10 BdSBP genes were identified as targets of miR156. Transcript cleavage analysis of selected BdSBPs by miR156 confirmed their antagonistic connection. Alternative splicing was observed playing an important role in BdSBPs and miR156 interaction. Characterization of T-DNA Bdsbp9 mutant showed reduced plant growth and spike length, reflecting its involvement in the spike development. Expression of a majority of BdSBPs elevated during spikelet initiation. Specifically, BdSBP1 and BdSBP3 differentially expressed in response to vernalization. Differential transcript abundance of BdSBP1, BdSBP3, BdSBP8, BdSBP9, BdSBP14, BdSBP18 and BdSBP23 genes was observed during the spike development under high temperature. Co-expression network, protein-protein interaction and biological pathway analysis indicate that BdSBP genes mainly regulate transcription, hormone, RNA and transport pathways. Our work reveals the multi-layered control of SBP genes and demonstrates their association with spike development and temperature sensitivity in Brachypodium.

Genome-wide identification, phylogeny and expression analysis of the SPL gene family in wheat

BACKGROUND

Members of the plant-specific SPL gene family (squamosa promoter-binding protein -like) contain the SBP conserved domain and are involved in the regulation of plant growth and development, including the development of plant flowers and plant epidermal hair, the plant stress response, and the synthesis of secondary metabolites. This family has been identified in various plants. However, there is no systematic analysis of the SPL gene family at the genome-wide level of wheat.

RESULTS

In this study, 56 putative TaSPL genes were identified using the comparative genomics method; we renamed them TaSPL001 - TaSPL056 on their chromosomal distribution. According to the un-rooted neighbor joining phylogenetic tree, gene structure and motif analyses, the 56 TaSPL genes were divided into 8 subgroups. A total of 81 TaSPL gene pairs were designated as arising from duplication events and 64 interacting protein branches were identified as involve in the protein interaction network. The expression patterns of 21 randomly selected TaSPL genes in different tissues (roots, stems, leaves and inflorescence) and under 4 treatments (abscisic acid, gibberellin, drought and salt) were detected by quantitative real-time polymerase chain reaction (qRT-PCR).

CONCLUSIONS

The wheat genome contains 56 TaSPL genes and those in same subfamily share similar gene structure and motifs. TaSPL gene expansion occurred through segmental duplication events. Combining the results of transcriptional and qRT-PCR analyses, most of these TaSPL genes were found to regulate inflorescence and spike development. Additionally, we found that 13 TaSPLs were upregulated by abscisic acid, indicating that TaSPL genes play a positive role in the abscisic acid-mediated pathway of the seedling stage. This study provides comprehensive information on the SPL gene family of wheat and lays a solid foundation for elucidating the biological functions of TaSPLs and improvement of wheat yield.

Multi-omics sequencing provides insight into floral transition in Catalpa bungei. C.A. Mey

BACKGROUND

Floral transition plays an important role in development, and proper time is necessary to improve the value of valuable ornamental trees. The molecular mechanisms of floral transition remain unknown in perennial woody plants. "Bairihua" is a type of C. bungei that can undergo floral transition in the first planting year.

RESULTS

Here, we combined short-read next-generation sequencing (NGS) and single-molecule real-time (SMRT) sequencing to provide a more complete view of transcriptome regulation during floral transition in C. bungei. The circadian rhythm-plant pathway may be the critical pathway during floral transition in early flowering (EF) C. bungei, according to horizontal and vertical analysis in EF and normal flowering (NF) C. bungei. SBP and MIKC-MADS-box were seemingly involved in EF during floral transition. A total of 61 hub genes were associated with floral transition in the MEturquoise model with Weighted Gene Co-expression Network Analysis (WGCNA). The results reveal that ten hub genes had a close connection with the GASA homologue gene (Cbu.gene.18280), and the ten co-expressed genes belong to five flowering-related pathways. Furthermore, our study provides new insights into the complexity and regulation of alternative splicing (AS). The ratio or number of isoforms of some floral transition-related genes is different in different periods or in different sub-genomes.

CONCLUSIONS

Our results will be a useful reference for the study of floral transition in other perennial woody plants. Further molecular investigations are needed to verify our sequencing data.

Molecular and functional characterization of the SBP-box transcription factor SPL-CNR in tomato fruit ripening and cell death

SlSPL-CNR, an SBP-box transcription factor (TF) gene residing at the epimutant Colourless non-ripening (Cnr) locus, is involved in tomato ripening. This epimutant provides a unique model to investigate the (epi)genetic basis of fruit ripening. Here we report that SlSPL-CNR is a nucleus-localized protein with a distinct monopartite nuclear localization signal (NLS). It consists of four consecutive residues ' 30KRKR33' at the N-terminus of the protein. Mutation of the NLS abolishes SlSPL-CNR's ability to localize in the nucleus. SlSPL-CNR comprises two zinc-finger motifs (ZFMs) within the C-terminal SBP-box domain. Both ZFMs contribute to zinc-binding activity. SlSPL-CNR can induce cell death in tomato and tobacco, dependent on its nuclear localization. However, the two ZFMs have differential impacts on SlSPL-CNR's induction of severe necrosis or mild necrotic ringspot. NLS and ZFM mutants cannot complement Cnr fruits to ripen. SlSPL-CNR interacts with SlSnRK1. Virus-induced SlSnRK1 silencing leads to reduction in expression of ripening-related genes and inhibits ripening in tomato. We conclude that SlSPL-CNR is a multifunctional protein that consists of a distinct monopartite NLS, binds to zinc, and interacts with SlSnRK1 to affect cell death and tomato fruit ripening.

Roles of transcription factor SQUAMOSA promoter binding protein-like gene family in papaya (Carica papaya) development and ripening

SQUAMOSA promoter binding protein-like (SPL) family plays vital regulatory roles in plant growth and development. The SPL family in climacteric fruit Carica papaya has not been reported. This study identified 14 papaya SPLs (CpSPL) from papaya genome and analyzed their sequence features, phylogeny, intron/exon structure, conserved motif, miR156-mediated posttranscriptional regulation, and expression patterns. 14 CpSPLs were clustered into 8 groups, and two distinct expression patterns were revealed for miR156-targeted and nontargeted CpSPLs in different tissues and fruit development stages. The expression changes of CpSPLs in ethephon and 1-MCP treated fruit during ripening suggested that the CpSPLs guided by CpmiR156 play crucial roles in ethylene signaling pathway. This study sheds light on the new function of SPL family in fruit development and ripening, providing insights on understanding evolutionary divergence of the members of SPL family among plant species.

Potential function of CbuSPL and gene encoding its interacting protein during flowering in Catalpa bungei

BACKGROUND

"Bairihua", a variety of the Catalpa bungei, has a large amount of flowers and a long flowering period which make it an excellent material for flowering researches in trees. SPL is one of the hub genes that regulate both flowering transition and development.

RESULTS

SPL homologues CbuSPL9 was cloned using degenerate primers with RACE. Expression studies during flowering transition in "Bairihua" and ectopic expression in Arabidopsis showed that CbuSPL9 was functional similarly with its Arabidopsis homologues. In the next step, we used Y2H to identify the proteins that could interact with CbuSPL9. HMGA, an architectural transcriptional factor, was identified and cloned for further research. BiFC and BLI showed that CbuSPL9 could form a heterodimer with CbuHMGA in the nucleus. The expression analysis showed that CbuHMGA had a similar expression trend to that of CbuSPL9 during flowering in "Bairihua". Intriguingly, ectopic expression of CbuHMGA in Arabidopsis would lead to aberrant flowers, but did not effect flowering time.

CONCLUSIONS

Our results implied a novel pathway that CbuSPL9 regulated flowering development, but not flowering transition, with the participation of CbuHMGA. Further investments need to be done to verify the details of this pathway.

Evolutionary and Comparative Expression Analyses of TCP Transcription Factor Gene Family in Land Plants

The plant-specific Teosinte-branched 1/Cycloidea/Proliferating (TCP) transcription factor genes are involved in plants' development, hormonal pathways, and stress response but their evolutionary history is uncertain. The genome-wide analysis performed here for 47 plant species revealed 535 TCP candidates in terrestrial plants and none in aquatic plants, and that TCP family genes originated early in the history of land plants. Phylogenetic analysis divided the candidate genes into Classes I and II, and Class II was further divided into CYCLOIDEA (CYC) and CINCINNATA (CIN) clades; CYC is more recent and originated from CIN in angiosperms. Protein architecture, intron pattern, and sequence characteristics were conserved in each class or clade supporting this classification. The two classes significantly expanded through whole-genome duplication during evolution. Expression analysis revealed the conserved expression of TCP genes from lower to higher plants. The expression patterns of Class I and CIN genes in different stages of the same tissue revealed their function in plant development and their opposite effects in the same biological process. Interaction network analysis showed that TCP proteins tend to form protein complexes, and their interaction networks were conserved during evolution. These results contribute to further functional studies on TCP family genes.

Genome-wide characterization of SPL family in Medicago truncatula reveals the novel roles of miR156/SPL module in spiky pod development

BACKGROUND

SQUAMOSA Promoter Binding Protein-Likes (SPLs) proteins are plant-specific transcription factors that play many crucial roles in plant growth and development. However, there is little information about SPL family in the model legume Medicago truncatula.

RESULTS

In this study, a total of 23 MtSPL genes were identified in M. truncatula genome, in which 17 of the MtSPLs contained the putative MtmiR156 binding site at the coding or 3' UTR regions. Tissue-specific expression pattern analysis showed that most MtmiR156-targeted MtSPLs were highly expressed in seed and pod. The observation of MtmiR156B-overexpressing plants reveals that MtmiR156/MtSPL modules are not only involved in the development of leaves and branches, but also in the seed pod development, especially the formation of spine on pod.

CONCLUSION

The spines on pods are developed in many plant species, which allow pods to adhere to the animals, and then be transported on the outside. This study sheds light on the new function of SPL family in seed dispersal by controlling the formation of spiky pod, and provides insights on understanding evolutionary divergence of the members of SPL gene family among plant species.

The architecture of the GhD7 promoter reveals the roles of GhD7 in growth, development and the abiotic stress response in rice

Grain number, plant height and heading date 7 (GhD7) is considered to be one of the key yield-related genes in the production of high-yielding and climate-ready super rice varieties. GhD7 delays the plant's flowering under long-day conditions, which ultimately results in increased yield. Recent findings indicate that GhD7 also plays a major role in the abiotic stress response; however, the fine regulatory mechanisms controlling Ghd7 expression have yet to be uncovered. This study was carried out to explore the transcription factor binding site (TFBS) architecture of the GhD7 promoter to identify the regulatory dynamics of GhD7 transcription. The promoter sequence (-2000 to +200 base pairs from the transcription start site) was retrieved from the PlantPAN 2.0 database. Ab initio promoter analysis, DNase I hypersensitive site (DHS) analysis, and methylation analysis were carried out to identify TFBSs. The TFBS diversity among rice varieties was also assessed. In addition to the previously identified 8 cis-elements, 448 novel cis-elements were identified in the GhD7 promoter that provide binding sites for 25 transcription factor families. Furthermore, a DNase I hypersensitive site and a CpG island were also identified. The identified transcription factor families include key transcription factors involved in both development and abiotic stress responses, revealing the regulatory dynamics of GhD7. Comparative analysis of multiple GhD7 promoters identified 31 single-nucleotide polymorphisms that result in TFBS variations among rice accessions. These variations are mostly found in relation to flowering and abiotic stress responsive TFBSs on the promoter. This study supports the model that GhD7 acts as a central regulator of rice growth, development, and the abiotic stress response.

Evolutionary Analyses Reveal Diverged Patterns of SQUAMOSA Promoter Binding Protein-Like (SPL) Gene Family in Oryza Genus

The SPL (SQUAMOSA promoter binding protein-like) gene family is one of the plant-specific transcription factor families and controls a considerable number of biological functions, including floral development, phytohormone signaling, and toxin resistance. However, the evolutionary patterns and driving forces of SPL genes in the Oryza genus are still not well-characterized. In this study, we investigated a total of 105 SPL genes from six AA genome Oryza representative species (O. barthii, O. glumipatula, O. nivara, O. rufipogon, O. glaberrima, and O. sativa). Phylogenetic and motif analyses indicated that SPL proteins could be divided into two distinct lineages (I and II), and further studies showed lineage II consisted of three clades (IIA, IIB, and IIC). We found that clade I had comparable structural features with clade IIA, whereas genes in clade IIC displayed intrinsic differences, such as lower exon numbers and the presence of miR156 regulation elements. Nineteen orthologous groups of OsSPLs in Oryza were also identified, and most exons within those genes maintained constant length, whereas length of intron changed relatively. All groups were constrained by stronger purifying selection and diversified continually including alterative gene number, intron length, and miR156 regulation. Subsequently, cis-acting element analyses revealed the potential role of SPLs in wild rice, which might participate in light-responsive, phytohormone response, and plant growth and development. Our results shed light on that different evolutionary rates and duplication events might result in divergent evolutionary patterns in each lineage of SPL genes, providing a guide in exploring diverse function in the rice gene family among six closely related Oryza species.

Identification and Expression Analyses of SBP-Box Genes Reveal Their Involvement in Abiotic Stress and Hormone Response in Tea Plant (Camellia sinensis)

The SQUAMOSA promoter binding protein (SBP)-box gene family is a plant-specific transcription factor family. This family plays a crucial role in plant growth and development. In this study, 20 SBP-box genes were identified in the tea plant genome and classified into six groups. The genes in each group shared similar exon-intron structures and motif positions. Expression pattern analyses in five different tissues demonstrated that expression in the buds and leaves was higher than that in other tissues. The cis-elements and expression patterns of the CsSBP genes suggested that the CsSBP genes play active roles in abiotic stress responses; these responses may depend on the abscisic acid (ABA), gibberellic acid (GA), and methyl jasmonate (MeJA) signaling pathways. Our work provides a comprehensive understanding of the CsSBP family and will aid in genetically improving tea plants.

Genome-wide identification and characterization of the SBP-box gene family in Petunia

Background

SQUAMOSA PROMOTER BINDING PROTEIN (SBP)-box genes encode a family of plant-specific transcription factors (TFs) that play important roles in many growth and development processes including phase transition, leaf initiation, shoot and inflorescence branching, fruit development and ripening etc. The SBP-box gene family has been identified and characterized in many species, but has not been well studied in Petunia, an important ornamental genus.

Results

We identified 21 putative SPL genes of Petunia axillaris and P. inflata from the reference genome of P. axillaris N and P. inflata S6, respectively, which were supported by the transcriptome data. For further confirmation, all the 21 genes were also cloned from P. hybrida line W115 (Mitchel diploid). Phylogenetic analysis based on the highly conserved SBP domains arranged PhSPLs in eight groups, analogous to those from Arabidopsis and tomato. Furthermore, the Petunia SPL genes had similar exon-intron structure and the deduced proteins contained very similar conserved motifs within the same subgroup. Out of 21 PhSPL genes, fourteen were predicted to be potential targets of PhmiR156/157, and the putative miR156/157 response elements (MREs) were located in the coding region of group IV, V, VII and VIII genes, but in the 3’-UTR regions of group VI genes. SPL genes were also identified from another two wild Petunia species, P. integrifolia and P. exserta, based on their transcriptome databases to investigate the origin of PhSPLs. Phylogenetic analysis and multiple alignments of the coding sequences of PhSPLs and their orthologs from wild species indicated that PhSPLs were originated mainly from P. axillaris. qRT-PCR analysis demonstrated differential spatiotemperal expression patterns of PhSPL genes in petunia and many were expressed predominantly in the axillary buds and/or inflorescences. In addition, overexpression of PhSPL9a and PhSPL9b in Arabidopsis suggested that these genes play a conserved role in promoting the vegetative-to-reproductive phase transition.

Conclusion

Petunia genome contains at least 21 SPL genes, and most of the genes are expressed in different tissues. The PhSPL genes may play conserved and diverse roles in plant growth and development, including flowering regulation, leaf initiation, axillary bud and inflorescence development. This work provides a comprehensive understanding of the SBP-box gene family in Petunia and lays a significant foundation for future studies on the function and evolution of SPL genes in petunia.

Electronic supplementary material

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

Annotation of the Corymbia terpene synthase gene family shows broad conservation but dynamic evolution of physical clusters relative to Eucalyptus

Terpenes are economically and ecologically important phytochemicals. Their synthesis is controlled by the terpene synthase (TPS) gene family, which is highly diversified throughout the plant kingdom. The plant family Myrtaceae are characterised by especially high terpene concentrations, and considerable variation in terpene profiles. Many Myrtaceae are grown commercially for terpene products including the eucalypts Corymbia and Eucalyptus. Eucalyptus grandis has the largest TPS gene family of plants currently sequenced, which is largely conserved in the closely related E. globulus. However, the TPS gene family has been well studied only in these two eucalypt species. The recent assembly of two Corymbia citriodora subsp. variegata genomes presents an opportunity to examine the conservation of this important gene family across more divergent eucalypt lineages. Manual annotation of the TPS gene family in C. citriodora subsp. variegata revealed a similar overall number, and relative subfamily representation, to that previously reported in E. grandis and E. globulus. Many of the TPS genes were in physical clusters that varied considerably between Eucalyptus and Corymbia, with several instances of translocation, expansion/contraction and loss. Notably, there was greater conservation in the subfamilies involved in primary metabolism than those involved in secondary metabolism, likely reflecting different selective constraints. The variation in cluster size within subfamilies and the broad conservation between the eucalypts in the face of this variation are discussed, highlighting the potential contribution of selection, concerted evolution and stochastic processes. These findings provide the foundation to better understand terpene evolution within the ecologically and economically important Myrtaceae.

Functional Characterization of the Steroid Reductase Genes GmDET2a and GmDET2b form Glycine max

Brassinosteroids are important phytohormones for plant growth and development. In soybean (Glycine max), BR receptors have been identified, but the genes encoding BR biosynthesis-related enzymes remain poorly understood. Here, we found that the soybean genome encodes eight steroid reductases (GmDET2a to GmDET2h). Phylogenetic analysis grouped 105 steroid reductases from moss, fern and higher plants into five subgroups and indicated that the steroid reductase family has experienced purifying selection. GmDET2a and GmDET2b, homologs of the Arabidopsis thaliana steroid 5 α -reductase AtDET2, are proteins of 263 amino acids. Ectopic expression of GmDET2a and GmDET2b rescued the defects of the Atdet2-1 mutant in both darkness and light. Compared to the mutant, the hypocotyl length and plant height of the transgenic lines GmDET2a and GmDET2b increased significantly, in both darkness and light, and the transcript levels of the BR biosynthesis-related genes CPD, DWF4, BR6ox-1 and BR6ox-2 were downregulated in GmDET2aOX-23 and GmDET2bOX-16 lines compared to that in Atdet2-1. Quantitative real-time PCR revealed that GmDET2a and GmDET2b are ubiquitously expressed in all tested soybean organs, including roots, leaves and hypocotyls. Moreover, epibrassinosteroid negatively regulated GmDET2a and GmDET2b expression. Sulfate deficiency downregulated GmDET2a in leaves and GmDET2b in leaves and roots; by contrast, phosphate deficiency upregulated GmDET2b in roots and leaves. Taken together, our results revealed that GmDET2a and GmDET2b function as steroid reductases.

Genome-Wide Analysis Characterization and Evolution of SBP Genes in Fragaria vesca, Pyrus bretschneideri, Prunus persica and Prunus mume

The SQUAMOSA promoter binding protein (SBP)-box proteins are plant-specific transcriptional factors in plants. SBP TFs are known to play important functions in a diverse development process and also related in the process of evolutionary novelties. SBP gene family has been characterized in several plant species, but little is known about molecular evolution, functional divergence and comprehensive study of SBP gene family in Rosacea. We carried out genome-wide investigations and identified 14, 32, 17, and 17 SBP genes from four Rosacea species (Fragaria vesca, Pyrus bretschneideri, Prunus persica and Prunus mume, respectively). According to phylogenetic analysis arranged the SBP protein sequences in seven groups. Localization of SBP genes presented an uneven distribution on corresponding chromosomes of Rosacea species. Our analyses designated that the SBP genes duplication events (segmental and tandem) and divergence. In addition, due to highly conserved structure pattern of SBP genes, recommended that highly conserved region of microsyneteny in the Rosacea species. Type I and II functional divergence was detected among various amino acids in SBP proteins, while there was no positive selection according to substitutional model analysis using PMAL software. These results recommended that the purifying selection might be leading force during the evolution process and dominate conservation of SBP genes in Rosacea species according to environmental selection pressure analysis. Our results will provide basic understanding and foundation for future research insights on the evolution of the SBP genes in Rosacea.

QTL mapping and candidate gene analysis of ferrous iron and zinc toxicity tolerance at seedling stage in rice by genome-wide association study

BACKGROUND

Ferrous iron (Fe) and zinc (Zn) at high concentration in the soil cause heavy metal toxicity and greatly affect rice yield and quality. To improve rice production, understanding the genetic and molecular resistance mechanisms to excess Fe and Zn in rice is essential. Genome-wide association study (GWAS) is an effective way to identify loci and favorable alleles governing Fe and Zn toxicty as well as dissect the genetic relationship between them in a genetically diverse population.

RESULTS

A total of 29 and 31 putative QTL affecting shoot height (SH), root length (RL), shoot fresh weight (SFW), shoot dry weight (SDW), root dry weight (RDW), shoot water content (SWC) and shoot ion concentrations (SFe or SZn) were identified at seedling stage in Fe and Zn experiments, respectively. Five toxicity tolerance QTL (qSdw3a, qSdw3b, qSdw12 and qSFe5 / qSZn5) were detected in the same genomic regions under the two stress conditions and 22 candidate genes for 10 important QTL regions were also determined by haplotype analyses.

CONCLUSION

Rice plants share partial genetic overlaps of Fe and Zn toxicity tolerance at seedling stage. Candidate genes putatively affecting Fe and Zn toxicity tolerance identified in this study provide valuable information for future functional characterization and improvement of rice tolerance to Fe and Zn toxicity by marker-assisted selection or designed QTL pyramiding.

Functional Conservation and Divergence among Homoeologs of TaSPL20 and TaSPL21, Two SBP-Box Genes Governing Yield-Related Traits in Hexaploid Wheat

Maintaining high and stable yields has become an increasing challenge in wheat breeding due to climate change. Although Squamosa-promoter binding protein (SBP)-box genes have important roles in plant development, very little is known about the actual biological functions of wheat SBP-box family members. Here, we dissect the functional conservation, divergence, and exploitation of homoeologs of two paralogous TaSPL wheat loci during domestication and breeding. TaSPL20 and TaSPL21 were highly expressed in the lemma and palea. Ectopic expressions of TaSPL20/21 in rice exhibited similar functions in terms of promoting panicle branching but had different functions during seed development. We characterized all six TaSPL20/21 genes located across the three homoeologous (A, B, and D) genomes. According to the functional analysis of naturally occurring variants in 20 environments, four favorable haplotypes were identified. Together, they reduced plant height by up to 27.5%, and TaSPL21-6D-HapII increased 1000-grain weight by 9.73%. Our study suggests that TaSPL20 and TaSPL21 homoeologs underwent diversification in function with each evolving its own distinctive characteristics. During domestication and breeding of wheat in China, favorable haplotypes of each set were selected and exploited to varying degrees due to their large effects on plant height and 1000-grain weight.

Transcriptome-Wide Identification and Prediction of miRNAs and Their Targets in Paris polyphylla var. yunnanensis by High-Throughput Sequencing Analysis

Long dormancy period of seeds limits the large-scale artificial cultivation of the scarce Paris polyphylla var. yunnanensis, an important traditional Chinese medicine. Characterizing miRNAs and their targets is crucial to understanding the role of miRNAs during seed dormancy in this species. Considering the limited genome information of this species, we first sequenced and assembled the transcriptome data of dormant seeds and their seed coats as the reference genome. A total of 146,671 unigenes with an average length of 923 bp were identified and showed functional diversity based on different annotation methods. Two small RNA libraries from respective seeds and seed coats were sequenced and the combining data indicates that 263 conserved miRNAs belonging to at least 83 families and 768 novel miRNAs in 1174 transcripts were found. The annotations of the predicted putative targets of miRNAs suggest that these miRNAs were mainly involved in the cell, metabolism and genetic information processing by direct and indirect regulation patterns in dormant seeds of P. polyphylla var. yunnanensis. Therefore, we provide the first known miRNA profiles and their targets, which will assist with further study of the molecular mechanism of seed dormancy in P. polyphylla var. yunnanensis.

Comparative genome-wide phylogenetic and expression analysis of SBP genes from potato (Solanum tuberosum)

Flowering time is a very important phase in transition to reproductive stage of life in higher plants. SQUAMOSA promoter-binding protein (SBP) gene family encodes plant-specific transcription factors that are involved in regulation of several developmental processes, especially flowering. Although SBP-box genes have been identified in different plants, there have been no study indicating the regulatory effect of SBP box in potato flowering. Here, we report for the first time the identification and characterization of SBP-box transcription factors as well as determination of expression level of SBP-box genes in Solanum tuberosum L. an important crop worldwide. Fifteen different SBP-box transcription factor genes were identified in potato genome. They were found to be distributed in nine chromosomes and eight of them included miR156 and miRNA157 target sites. Characterization of amino acid sequences were performed and protein interactions were predicted. In addition, expression levels of five S. tuberosum SBP-box genes were analysed by both in silico and qRT-PCR. All these results provide a better understanding of functional role of SBP-box gene family members in flowering time in potato.

Genomic identification, characterization and differential expression analysis of SBP-box gene family in Brassica napus

BACKGROUND

SBP-box genes belong to one of the largest families of transcription factors. Though members of this family have been characterized to be important regulators of diverse biological processes, information of SBP-box genes in the third most important oilseed crop Brassica napus is largely undefined.

RESULTS

In the present study, by whole genome bioinformatics analysis and transcriptional profiling, 58 putative members of SBP-box gene family in oilseed rape (Brassica napus L.) were identified and their expression pattern in different tissues as well as possible interaction with miRNAs were analyzed. In addition, B. napus lines with contrasting branch angle were used for investigating the involvement of SBP-box genes in plant architecture regulation. Detailed gene information, including genomic organization, structural feature, conserved domain and phylogenetic relationship of the genes were systematically characterized. By phylogenetic analysis, BnaSBP proteins were classified into eight distinct groups representing the clear orthologous relationships to their family members in Arabidopsis and rice. Expression analysis in twelve tissues including vegetative and reproductive organs showed different expression patterns among the SBP-box genes and a number of the genes exhibit tissue specific expression, indicating their diverse functions involved in the developmental process. Forty-four SBP-box genes were ascertained to contain the putative miR156 binding site, with 30 and 14 of the genes targeted by miR156 at the coding and 3'UTR region, respectively. Relative expression level of miR156 is varied across tissues. Different expression pattern of some BnaSBP genes and the negative correlation of transcription levels between miR156 and its target BnaSBP gene were observed in lines with different branch angle.

CONCLUSIONS

Taken together, this study represents the first systematic analysis of the SBP-box gene family in Brassica napus. The data presented here provides base foundation for understanding the crucial roles of BnaSBP genes in plant development and other biological processes.

Functional Characterization of Soybean Glyma04g39610 as a Brassinosteroid Receptor Gene and Evolutionary Analysis of Soybean Brassinosteroid Receptors

Brassinosteroids (BR) play important roles in plant growth and development. Although BR receptors have been intensively studied in Arabidopsis, the BR receptors in soybean remain largely unknown. Here, in addition to the known receptor gene Glyma06g15270 (GmBRI1a), we identified five putative BR receptor genes in the soybean genome: GmBRI1b, GmBRL1a, GmBRL1b, GmBRL2a, and GmBRL2b. Analysis of their expression patterns by quantitative real-time PCR showed that they are ubiquitously expressed in primary roots, lateral roots, stems, leaves, and hypocotyls. We used rapid amplification of cDNA ends (RACE) to clone GmBRI1b (Glyma04g39160), and found that the predicted amino acid sequence of GmBRI1b showed high similarity to those of AtBRI1 and pea PsBRI1. Structural modeling of the ectodomain also demonstrated similarities between the BR receptors of soybean and Arabidopsis. GFP-fusion experiments verified that GmBRI1b localizes to the cell membrane. We also explored GmBRI1b function in Arabidopsis through complementation experiments. Ectopic over-expression of GmBRI1b in Arabidopsis BR receptor loss-of-function mutant (bri1-5 bak1-1D) restored hypocotyl growth in etiolated seedlings; increased the growth of stems, leaves, and siliques in light; and rescued the developmental defects in leaves of the bri1-6 mutant, and complemented the responses of BR biosynthesis-related genes in the bri1-5 bak1-D mutant grown in light. Bioinformatics analysis demonstrated that the six BR receptor genes in soybean resulted from three gene duplication events during evolution. Phylogenetic analysis classified the BR receptors in dicots and monocots into three subclades. Estimation of the synonymous (Ks) and the nonsynonymous substitution rate (Ka) and selection pressure (Ka/Ks) revealed that the Ka/Ks of BR receptor genes from dicots and monocots were less than 1.0, indicating that BR receptor genes in plants experienced purifying selection during evolution.

EcoTILLING-Based Association Mapping Efficiently Delineates Functionally Relevant Natural Allelic Variants of Candidate Genes Governing Agronomic Traits in Chickpea

The large-scale mining and high-throughput genotyping of novel gene-based allelic variants in natural mapping population are essential for association mapping to identify functionally relevant molecular tags governing useful agronomic traits in chickpea. The present study employs an alternative time-saving, non-laborious and economical pool-based EcoTILLING approach coupled with agarose gel detection assay to discover 1133 novel SNP allelic variants from diverse coding and regulatory sequence components of 1133 transcription factor (TF) genes by genotyping in 192 diverse desi and kabuli chickpea accessions constituting a seed weight association panel. Integrating these SNP genotyping data with seed weight field phenotypic information of 192 structured association panel identified eight SNP alleles in the eight TF genes regulating seed weight of chickpea. The associated individual and combination of all SNPs explained 10-15 and 31% phenotypic variation for seed weight, respectively. The EcoTILLING-based large-scale allele mining and genotyping strategy implemented for association mapping is found much effective for a diploid genome crop species like chickpea with narrow genetic base and low genetic polymorphism. This optimized approach thus can be deployed for various genomics-assisted breeding applications with optimal expense of resources in domesticated chickpea. The seed weight-associated natural allelic variants and candidate TF genes delineated have potential to accelerate marker-assisted genetic improvement of chickpea.

Switchgrass SBP-box transcription factors PvSPL1 and 2 function redundantly to initiate side tillers and affect biomass yield of energy crop

BACKGROUND

Switchgrass (Panicum virgatum L.) is a dedicated lignocellulosic feedstock for bioenergy production. The SQUAMOSA PROMOTER-BINDING PROTEIN (SBP-box)-LIKE transcription factors (SPLs) change plant architecture and vegetative-to-reproductive phase transition significantly, and as such, they are promising candidates for genetic improvement of switchgrass biomass yield. However, the genome-wide identification and functional characterization of SPL genes have yet to be investigated in herbaceous energy crops.

RESULTS

We identified 35 full-length SPL genes in the switchgrass genome. The phylogenetic relationship and expression pattern of PvSPLs provided baseline information for their function characterization. Based on the global overview of PvSPLs, we explored the biological function of miR156-targeted PvSPL1 and PvSPL2, which are closely related members of SPL family in switchgrass. Our results showed that PvSPL1 and PvSPL2 acted redundantly to modulate side tiller initiation, whereas they did not affect phase transition and internode initiation. Consistently, overexpression of the miR156-resistant rPvSPL2 in the miR156-overexpressing transgenic plants greatly reduced tiller initiation, but did not rescue the delayed flowering and increased internode numbers. Furthermore, suppression of PvSPL2 activity in switchgrass increased biomass yield and reduced lignin accumulation, which thereby elevated the total amount of solubilized sugars.

CONCLUSIONS

Our results indicate that different miR156-targeted PvSPL subfamily genes function predominantly in certain biological processes in switchgrass. We suggest that PvSPL2 and its paralogs can be utilized as the valuable targets in molecular breeding of energy crops for developing novel germplasms with high biofuel production.