Genome-wide identification and characterization of the SPL gene family in Ziziphus jujuba

Infection with Jujube Witches' Broom Phytoplasma Alters the Expression Pattern of the Argonaute Gene Family in Ziziphus jujuba

The cultivation of jujube (Ziziphus jujuba) in China is threatened by jujube witches' broom (JWB) disease, a devastating infectious disease associated with JWB phytoplasma ('Candidatus Phytoplasma ziziphi'). In many plants, proteins in the Argonaute (AGO) family, as main components of the RNA-induced silencing complex (RISC), play important roles in RNA silencing and pathogen resistance. The jujube telomere-to-telomere genome was searched by BLAST using Arabidopsis AGOs as probes. A total of nine jujube AGO gene members were identified, with each containing the conserved N-terminal, PZA, and PIWI domains. Phylogenetic analysis revealed that the nine jujube AGOs scattered into all three Arabidopsis AGO clades. Expression patterns of the ZjAGO genes were analyzed in response to phytoplasma in transcriptome data and by qRT-PCR. The jujube-phytoplasma interaction altered the expression of jujube AGO genes. ZjAGO1 and ZjAGO8 were up-regulated in the majority of the eight sampling periods subjected to qRT-PCR analysis. In the transcriptome data, ZjAGO1 and ZjAGO8 were also up-regulated during the key stages 37 and 39 weeks after grafting (WAG) with phytoplasma-infected material. These two jujube Argonaute genes may play important roles in response to JWB phytoplasma infection.

Mechanism of zju-miR156c-mediated network in regulating witches' broom symptom of Chinese jujube

Jujube witches' broom, caused by phytoplasma, is a destructive disease of Chinese jujube. Studies have shown that zju-miR156s play an important role in phytoplasma infection in jujube, but the regulatory mechanism between zju-miR156c and witches' broom remains unexplored. In the current study, miRNA-seq and gene expression analysis showed that zju-miR156c was more highly induced in infected jujube plants than the other miRNAs and its target gene was ZjSPL3. In addition, the expression levels of thymidylate kinase gene (TMKJWB) and secreted jujube protein (SJP1JWB) in diseased materials were higher than those in healthy controls. The expression level of zju-miR156c was significantly upregulated, while ZjSPL3 was sharply downregulated and the content of cytokinin (CTK) significantly increased. Overexpression of zju-miR156c in Arabidopsis significantly reduced the expression of AtSPL10 (homologous gene of ZjSPL3) but increased the content of CTK, and the transgenic plants exhibited witches' broom symptoms. In addition, yeast two-hybrid and co-immunoprecipitation assays confirmed that SJP1JWB interacted with ZjERF18. Yeast one-hybrid analysis showed that ZjERF18 could interact with the promoter of zju-MIR156c. In conclusion, our results demonstrated a novel pathogenic module of ZjERF18-zju-miR156c-ZjSPL3-CTK has an important function in the formation of witches' broom caused by SJP1JWB.

Overexpression of MiSPL3a and MiSPL3b confers early flowering and stress tolerance in Arabidopsis thaliana

SQUAMOSA promoter-binding protein-like (SPL) family genes play an important role in regulating plant flowering and resistance to stress. However, understanding the function of the SPL family in mango is still limited. In a previous study, two MiSPL3 genes, MiSPL3a and MiSPL3b (MiSPL3a/b), were identified in 'SiJiMi' mango and exhibited the highest expression in flowers at the initial flowering stage [24]. Therefore, in this study, we further investigated the expression pattern and gene function of MiSPL3a/b. The results showed that the expression of MiSPL3a was greatest at the end of floral bud differentiation, and MiSPL3b was expressed mainly during the flowering induction and vegetative growth stages. Subcellular localization showed that MiSPL3a/b localized to the nucleus. In addition, ectopic expression of MiSPL3a/b promoted earlier flowering in Arabidopsis thaliana by 3 d-6 d than in wild-type (WT) plants, which increased the expression of SUPPRESSOR OF CONSTANS1 (AtSOC1), FRUITFULL (AtFUL), and APETALA1 (AtAP1). MiSPL3a/b transgenic lines exhibited increased tolerance to drought, GA3, and abscisic acid (ABA) treatments but were sensitive to Pro-Ca treatment. Furthermore, protein interaction analysis revealed that MiSPL3a/b could interact with several stress-related proteins, flowering-related proteins, and the bridge protein 14-3-3. Taken together, MiSPL3a and MiSPL3b acted as positive regulators of flowering time and stress tolerance in transgenic Arabidopsis.

Genome-wide identification, stress- and hormone-responsive expression characteristics, and regulatory pattern analysis of Scutellaria baicalensis SbSPLs

SQUAMOSA PROMOTER BINDING PROTEIN-LIKEs (SPLs) encode plant-specific transcription factors that regulate plant growth and development, stress response, and metabolite accumulation. However, there is limited information on Scutellaria baicalensis SPLs. In this study, 14 SbSPLs were identified and divided into 8 groups based on phylogenetic relationships. SbSPLs in the same group had similar structures. Abscisic acid-responsive (ABRE) and MYB binding site (MBS) cis-acting elements were found in the promoters of 8 and 6 SbSPLs. Segmental duplications and transposable duplications were the main causes of SbSPL expansion. Expression analysis based on transcriptional profiling showed that SbSPL1, SbSPL10, and SbSPL13 were highly expressed in roots, stems, and flowers, respectively. Expression analysis based on quantitative real-time polymerase chain reaction (RT‒qPCR) showed that most SbSPLs responded to low temperature, drought, abscisic acid (ABA) and salicylic acid (SA), among which the expression levels of SbSPL7/9/10/12 were significantly upregulated in response to abiotic stress. These results indicate that SbSPLs are involved in the growth, development and stress response of S. baicalensis. In addition, 8 Sba-miR156/157 s were identified, and SbSPL1-5 was a potential target of Sba-miR156/157 s. The results of target gene prediction and coexpression analysis together indicated that SbSPLs may be involved in the regulation of L-phenylalanine (L-Phe), lignin and jasmonic acid (JA) biosynthesis. In summary, the identification and characterization of the SbSPL gene family lays the foundation for functional research and provides a reference for improved breeding of S. baicalensis stress resistance and quality traits.

CsSPL13A directly binds and positively regulates CsFT and CsBAM to accelerate flowering in cucumber

Flowering is an important developmental transition that greatly affects the yield of many vegetable crops. In cucumber (Cucumis sativus), flowering is regulated by various factors including squamosa promoter-binding-like (SPL) family proteins. However, the role of CsSPL genes in cucumber flowering remains largely unknown. In this study, we cloned the squamosa promoter-binding-like protein 13A (CsSPL13A) gene, which encodes a highly conserved SBP-domain protein that acts as a transcription factor and localizes to the nucleus. Quantitative real-time PCR (qRT-PCR) analysis showed that CsSPL13A was mainly expressed in flowers, and its expression level increased significantly nearing the flowering stage. Additionally, compared with the wild type(WT), CsSPL13A-overexpressing transgenic cucumber plants (CsSPL13A-OE) showed considerable differences in flowering phenotypes, such as early flowering, increased number of male flowers, and longer flower stalks. CsSPL13A upregulated the expression of the flowering integrator gene Flowering Locus T (CsFT) and the sugar-mediated flowering gene β-amylase (CsBAM) in cucumber. Yeast one-hybrid and firefly enzyme reporter assays confirmed that the CsSPL13A protein could directly bind to the promoters of CsFT and CsBAM, suggesting that CsSPL13A works together with CsFT and CsBAM to mediate flowering in cucumber. Overall, our results provide novel insights into the regulatory network of flowering in cucumber as well as new ideas for the genetic improvement of cucumber varieties.

Genome-Wide Identification and Expression Analysis of the SQUAMOSA Promoter-Binding Protein-like (SPL) Transcription Factor Family in Catalpabungei

As a plant-specific transcription factor, the SPL gene family plays a critical role in plant growth and development. Although the SPL gene family has been identified in diverse plant species, there have been no genome-wide identification or systematic study reports on the SPL gene family in Catalpa bungei. In this study, we identified 19 putative SPL gene family members in the C. bungei genome. According to the phylogenetic relationship, they can be divided into eight groups, and the genes in the same group have a similar gene structure and conserved motifs. Synteny analysis showed that fragment duplication played an important role in the expansion of the CbuSPL gene family. At the same time, CbuSPL genes have cis-acting elements and functions related to light response, hormone response, growth and development, and stress response. Tissue-specific expression and developmental period-specific expression analysis showed that CbuSPL may be involved in flowering initiation and development, flowering transition, and leaf development. In addition, the ectopic expression of CbuSPL4 in Arabidopsis confirmed that it can promote early flowering and induce the expression of related flowering genes. These systematic research results will lay a foundation for further study on the functional analysis of SPL genes in C. bungei.

Differential Expression of miRNAs Involved in Response to Candidatus Liberibacter asiaticus Infection in Mexican Lime at Early and Late Stages of Huanglongbing Disease

Huanglongbing (HLB) is one of the most destructive diseases threatening citriculture worldwide. This disease has been associated with α-proteobacteria species, namely Candidatus Liberibacter. Due to the unculturable nature of the causal agent, it has been difficult to mitigate the disease, and nowadays a cure is not available. MicroRNAs (miRNAs) are key regulators of gene expression, playing an essential role in abiotic and biotic stress in plants including antibacterial responses. However, knowledge derived from non-model systems including Candidatus Liberibacter asiaticus (CLas)-citrus pathosystem remains largely unknown. In this study, small RNA profiles from Mexican lime (Citrus aurantifolia) plants infected with CLas at asymptomatic and symptomatic stages were generated by sRNA-Seq, and miRNAs were obtained with ShortStack software. A total of 46 miRNAs, including 29 known miRNAs and 17 novel miRNAs, were identified in Mexican lime. Among them, six miRNAs were deregulated in the asymptomatic stage, highlighting the up regulation of two new miRNAs. Meanwhile, eight miRNAs were differentially expressed in the symptomatic stage of the disease. The target genes of miRNAs were related to protein modification, transcription factors, and enzyme-coding genes. Our results provide new insights into miRNA-mediated regulation in C. aurantifolia in response to CLas infection. This information will be useful to understand molecular mechanisms behind the defense and pathogenesis of HLB.

The NtSPL Gene Family in Nicotiana tabacum: Genome-Wide Investigation and Expression Analysis in Response to Cadmium Stress

The SQUAMOSA promoter binding protein-like (SPL)SPL family genes play an important role in regulating plant growth and development, synthesis of secondary metabolites, and resistance to stress. Understanding of the role of the SPL family in tobacco is still limited. In this study, 42 NtSPL genes were identified from the genome of the tobacco variety TN90. According to the results of the conserved motif and phylogenetic tree, the NtSPL genes were divided into eight subgroups, and the genes in the same subgroup showed similar gene structures and conserved domains. The cis-acting element analysis of the NtSPL promoters showed that the NtSPL genes were regulated by plant hormones and stresses. Twenty-eight of the 42 NtSPL genes can be targeted by miR156. Transcriptome data and qPCR results indicated that the expression pattern of miR156-targeted NtSPL genes was usually tissue specific. The expression level of miR156 in tobacco was induced by Cd stress, and the expression pattern of NtSPL4a showed a significant negative correlation with that of miR156. These results suggest that miR156-NtSPL4a may mediate the tobacco response to Cd stress. This study lays a foundation for further research on the function of the NtSPL gene and provides new insights into the involvement of NtSPL genes in the plant response to heavy metal stress.

Genome-Wide Analysis of SQUAMOSA-Promoter-Binding Protein-like Family in Flowering Pleioblastus pygmaeus

SQUAMOSA Promoter-Binding Protein-Like (SPL) family is well-known for playing an important role in plant growth and development, specifically in the reproductive process. Bamboo plants have special reproductive characteristics with a prolonged vegetative phase and uncertain flowering time. However, the underlying functions of SPL genes in reproductive growth are undisclosed in bamboo plants. In the study, a total of 28 SPLs were screened from an ornamental dwarf bamboo species, Pleioblastus pygmaeus. Phylogenetic analysis indicates that 183 SPLs from eight plant species can be classified into nine subfamilies, and the 28 PpSPLs are distributed among eight subfamilies. Homologous analysis shows that as many as 32 pairs of homologous genes were found between P. pygmaeus and rice, and 83 pairs were found between P. pygmaeus and Moso bamboo, whose Ka/Ks values are all <1. MiRNA target prediction reveals that 13 out of the 28 PpSPLs have recognition sites complementary to miRNA156. To screen the SPLs involved in the reproductive growth of bamboo plants, the mRNA abundance of the 28 PpSPLs was profiled in the different tissues of flowering P. pygmaeus and non-flowering plants by RNA-Seq. Moreover, the relative expression level of eight PpSPLs is significantly higher in flowering P. pygmaeus than that in non-flowering plants, which was also validated by RT-qPCR. Combined with phylogenetic analysis and homologous analysis, the eight significant, differentially expressed PpSPLs were identified to be associated with the reproductive process and flower organ development. Among them, there are four potential miRNA156-targeting PpSPLs involved in the flowering process. Of significant interest in the study is the identification of 28 SPLs and the exploration of four key flowering-related SPLs from P. pygmaeus, which provides a theoretic basis for revealing the underlying functions of SPLs in the reproductive growth of bamboo plants.

Genomic constitution and relationships of Ziziphus species collected from Malakand Division, Pakistan

Polyploidization has played a major role in plant evolution and can alter plant morphology, phenology, and ecology within only one or a few generations. Ziziphus species are economically as well as nutritionally important fruit-yielding trees. Identification of genotypes with unique traits or those with higher ploidy levels or a broader genetic base could lead to further improvements within the species. The current study has assessed the ploidy levels in the Ziziphus species (Ziziphus jujuba Mill. and Ziziphus nummularia (Burm. f. Wight & Arn) with phenotypic traits, flow cytometry, and chromosomal count as well as with SSRs markers. Morphological traits were inferred to be the most important drivers of trait variations among the investigated genotypes. The total sugar, total cAMPs, titratable acid, and chlorophyll (a, b, and total) were also significantly different in contrast with diploid plants, which showed that tetraploid Ziziphus had the potential to increase nutritional contents. Out of twenty (20), five (5) Z. jujuba genotypes (ZJL-9, ZJL-12, ZJL-17, ZJL-18, and ZJL-19) were found tetraploid 2n = 4x = 48, with genome size ranging from 965.9 to1238.8 Mb that was significantly higher than the tetraploid Z. jujuba Mill. variety Dongzao. Similarly, Z. nummularia ZNL-07 to ZNL-15 have found tetraploid 2n = 4x = 72 with genomic sizes ranging from 1152.2 to 1746.8 Mb respectively. Simple sequence repeats (SSRs) marker was applied to assess the genetic relationship within Ziziphus genotypes. To the best of our understanding, this is the first report on the identification of naturalized random tetraploids within the Pakistani Ziziphus species. This study provides important insights into the genomic architecture of Ziziphus species with implications for classification, conservation, and improvements of Ziziphus germplasm resources.

Genome-Wide Analysis of the SPL Gene Family and Expression Analysis during Flowering Induction in Prunus × yedoensis 'Somei-yoshino'

SQUAMOSA Promoter-Binding Protein-Like (SPL) genes encode plant-specific transcription factors which bind to the SQUAMOSA promoter of the MADS-box genes to regulate its expression. It plays important regulatory roles in floral induction and development, fertility, light signals and hormonal transduction, and stress response in plants. In this study, 32 PySPL genes with complete SBP (squamosa promoter binding protein) conserved domain were identified from the genome of Prunus × yedoensis ‘Somei-yoshino’ and analyzed by bioinformatics. 32 PySPLs were distributed on 13 chromosomes, encoding 32 PySPL proteins with different physical and chemical properties. The phylogenetic tree constructed with Arabidopsis thaliana and Oryza sativa can be divided into 10 subtribes, indicating PySPLs of different clusters have different biological functions. The conserved motif prediction showed that the number and distribution of motifs on each PySPL is varied. The gene structure analysis revealed that PySPLs harbored exons ranging from 2 to 10. The predictive analysis of acting elements showed that the promoter of PySPLs contain a large number of light-responsive elements, as well as response elements related to hormone response, growth and development and stress response. The analysis of the PySPLs expressions in flower induction and flower organs based on qRT-PCR showed that PySPL06/22 may be the key genes of flower development, PySPL01/06 and PySPL22 may play a role in the development of sepal and pistil, respectively. The results provide a foundation for the study of SPL transcription factors of Prunus × yedoensis ‘Somei-yoshino’ and provide more reference information of the function of SPL gene in flowering.

Transcriptome Profiling and Identification of the Candidate Genes Involved in Early Ripening in Ziziphus Jujuba

The early ripening jujube is an immensely popular fresh fruit due to its high commercial value as well as rich nutrition. However, little is known about the mechanism of jujube fruit’s ripening. In this study, the transcriptome profiles were comprehensively analyzed between the ‘Lingwu Changzao’ jujube and its early-ripening mutant during the fruit development and maturity. A total of 5,376 and 762 differentially expressed genes (DEGs) were presented at 80 and 90 days after the flowering of the jujube fruit, respectively. Furthermore, 521 common DEGs were identified as candidate genes that might be associated with the fruit’s early ripening. Our findings demonstrated that in a non-climacteric jujube fruit, abscisic acid (ABA) was more greatly involved in fruit ripening than ethylene. Meanwhile, the fruit ripening of the early-ripening mutant was regulated by eight promotors of DEGs related to glucose and fructose, seven repressors of DEGs related to brassinosteroid signal transduction, and a series of transcription factor genes (MYB, Bhlh, and ERF). Additionally, the expression of 20 candidate DEGs was further validated by real-time PCR during the late fruit maturation stage. Collectively, the present study sheds light on the metabolic mechanism of the fruit’s early ripening and provides valuable candidate genes for the early-ripening mutant’s breeding.

Genome-Wide Analysis and Expression Profile of Nuclear Factor Y (NF-Y) Gene Family in Z. jujuba

Nuclear factor-Y (NF-Y) is an important transcription factor in the plant species, which potentially provides a higher level of functional diversity including for abiotic stress tolerance. The genome-wide study and expression analysis of NF-Y gene family in Ziziphus, an elite abiotic stress-tolerant species, assist bioprospecting of genes. Here, a total of 32 NF-Y (8 NF-YA, 15 NF-YB, and 9 NF-YC) genes were identified in genome-wide search of Z. jujuba genome. Physicochemical properties, cellular localization, gene structure, chromosomal location, and protein motifs were analyzed for structural and functional understanding. Identified 12 NF-Ys were responsible for the expansion of NF-Y gene family by tandem duplication in Z. jujuba. Phylogenetic and comparative physical mapping of Z. jujuba NF-Ys with its orthologs illustrated evolutionary and functional insights into NF-Y gene family. A total of 45 perfect microsatellites (20bp to 40bp) were extracted across the ZjNF-Y genes. The promoter and gene ontology study suggested that Z. jujuba NF-Y gene family is functionally diverse and could play a wide-ranging role in plant abiotic stress, development, and cellular processes. An expression study revealed that large numbers of the NF-Ys are differentially expressed in response to drought and salinity. The total 15 and 18 ZjNF-Y genes that are upregulated under drought and salinity stress, respectively, are the potential candidates for further functional analysis for development of climate-resilient crops. The present study established a base for understanding the role of NF-Ys in Z. jujuba under abiotic stress conditions and paved a way for further research.

Phytoplasma effector Zaofeng6 induces shoot proliferation by decreasing the expression of ZjTCP7 in Ziziphus jujuba

The jujube witches' broom (JWB) phytoplasma is associated with witches' broom, dwarfism, and smaller leaves in jujube, resulting in yield losses. In this study, eight putative JWB effector proteins were identified from potential mobile units of the JWB genome. Among them, Zaofeng6 induced witches' broom symptoms in Arabidopsis and jujube. Zaofeng6-overexpressing Arabidopsis and unrooted jujube transformants displayed witches' broom-like shoot proliferation. Transient expression of Zaofeng6 induced hypersensitive response like cell death and expression of hypersensitive response marker genes, like harpin-induced gene 1 (H1N1), and the pathogenesis-related genes PR1, PR2, and PR3 in transformed Nicotiana benthamiana leaves, suggesting that Zaofeng6 could be a virulence effector. Yeast two-hybrid library screening and bimolecular fluorescence complementation confirmed that Zaofeng6 interacts with ZjTCP7 through its first two α-helix domains in the cell nuclei. ZjTCP7 mRNA and protein abundance decreased in Zaofeng6 transgenic jujube seedlings. The expression of some genes in the strigolactone signaling pathway (ZjCCD7, ZjCCD8, and CYP711A1) were down-regulated in jujube shoots overexpressing Zaofeng6 and in zjtcp7 CRISPR/Cas9 mutants. Zaofeng6 induces shoot proliferation through decreased expression of ZjTCP7 at the transcriptional and translational levels.

Characterization and expression analysis of the SPL gene family during floral development and abiotic stress in pecan (Carya illinoinensis)

SQUAMOSA promoter binding protein-like (SPL) genes are a type of plant-specific transcription factors that play crucial roles in the regulation of phase transition, floral transformation, fruit development, and various stresses. Although SPLs have been characterized in several model species, no systematic analysis has been studied in pecans, an important woody oil tree species. In this study, a total of 32 SPL genes (CiSPLs) were identified in the pecan genome. After conducting phylogenetic analysis of the conserved SBP proteins from Arabidopsis, rice, and poplar, the CiSPLs were separated into eight subgroups. The CiSPL genes within the same subgroup contained very similar exon-intron structures and conserved motifs. Nine segmentally duplicated gene pairs in the pecan genome and 16 collinear gene pairs between the CiSPL and AtSPL genes were identified. Cis-element analysis showed that CiSPL genes may regulate plant meristem differentiation and seed development, participate in various biological processes, and respond to plant hormones and environmental stresses. Therefore, we focused our study on the expression profiles of CiSPL genes during flower and fruit development. Most of the CiSPL genes were predominantly expressed in buds and/or female flowers. Additionally, quantitative real time PCR (qRT-PCR) analyses confirmed that CiSPL genes showed distinct spatiotemporal expression patterns in response to drought and salt treatments. The study provides foundation for the further exploration of the function and evolution of SPL genes in pecan.

Transcription factor SmSPL7 promotes anthocyanin accumulation and negatively regulates phenolic acid biosynthesis in Salvia miltiorrhiza

Plant-specific SQUAMOSA promoter-binding protein-like (SPL) transcription factors play critical regulatory roles during plant growth and development. However, the functions of SPLs in Salvia miltiorrhiza (SmSPLs; a model medicinal plant) have not been reported. Here, the expression patterns and functions of SmSPL7 were characterized in S. miltiorrhiza. SmSPL7 was expressed in all parts of S. miltiorrhiza, with the highest expression level in the leaves, and could be inhibited by multiple hormones, including methyl jasmonate, auxin, abscisic acid, and gibberellin. SmSPL7 is localized within the nucleus and exhibits robust transcriptional activation activity. Transgenic lines overexpressing SmSPL7 demonstrated pronounced growth inhibition, accompanied by increased anthocyanin accumulation via the genetic activation of the anthocyanin biosynthesis pathway. However, SmSPL7 overexpression significantly decreased salvianolic acid B (SalB) production by inhibiting the transcripts of genes implicated in its biosynthesis pathway. Further analysis indicated that SmSPL7 directly binds to SmTAT1 and Sm4CL9 promoters and blocks their expression to inhibit the biosynthesis of SalB. Taken together, these results indicate that SmSPL7 is a negative regulator of SalB biosynthesis but positively regulates anthocyanin accumulation in S. miltiorrhiza. These findings provide new insights into the functionality of the SPL family while establishing an important foundation for further uncovering the crucial roles of SmSPL7 in the growth of S. miltiorrhiza.

Hierarchical Action of Mulberry miR156 in the Vegetative Phase Transition

The vegetative phase transition is a prerequisite for flowering in angiosperm plants. Mulberry miR156 has been confirmed to be a crucial factor in the vegetative phase transition in Arabidopsis thaliana. The over-expression of miR156 in transgenic Populus × canadensis dramatically prolongs the juvenile phase. Here, we find that the expression of mno-miR156 decreases with age in all tissues in mulberry, which led us to study the hierarchical action of miR156 in mulberry. Utilizing degradome sequencing and dual-luciferase reporter assays, nine MnSPLs were shown to be directly regulated by miR156. The results of yeast one-hybrid and dual-luciferase reporter assays also revealed that six MnSPLs could recognize the promoter sequences of mno-miR172 and activate its expression. Our results demonstrate that mno-miR156 performs its role by repressing MnSPL/mno-miR172 pathway expression in mulberry. This work uncovered a miR156/SPLs/miR172 regulation pathway in the development of mulberry and fills a gap in our knowledge about the molecular mechanism of vegetative phase transition in perennial woody plants.

SQUAMOSA Promoter Binding Protein-Like (SPL) Gene Family: TRANSCRIPTOME-Wide Identification, Phylogenetic Relationship, Expression Patterns and Network Interaction Analysis in Panax ginseng C. A. Meyer

: SPL (SQUAMOSA promoter binding protein-like) gene family is specific transcription factor in the plant that have an important function for plant growth and development. Although the SPL gene family has been widely studied and reported in many various plant species from gymnosperm to angiosperm, there are no systematic studies and reports about the SPL gene family in Panax ginseng C. A. Meyer. In this study, we conducted transcriptome-wide identification, evolutionary analysis, structure analysis, and expression characteristics analysis of SPL gene family in Panax ginseng by bioinformatics. We annotated the PgSPL gene family and found that they might involve in multiple functions including encoding structural proteins, but the main function were still focused on the binding function. The result showed that 106 PgSPL transcripts were classified into two clades - A and B, both of which respectively consisted of three groups. Besides, we profiled PgSPL transcripts' genotypic, temporal, and spatial expression characteristics. Furthermore, we calculated the correlation of PgSPL transcripts in the 14 tissues of a 4 years old ginseng and 42 farmers' cultivars farmers' cultivars of 4 years old ginsengs' roots with both results showing that SPL transcripts formed a single network, which indicated that PgSPLs inter-coordinated when performing their functions. What's more, we found that most PgSPL transcripts tended to express in older ginseng instead of younger ginseng, which was not only reflected in the expression of more types of SPL transcripts in older ginseng, but also in the higher expression of SPL transcripts in older ginseng. Additionally, we found that four PgSPL transcripts were only massively expressed in roots. According to PgSPL transcripts' expression characteristics, we found that PgSPL23-35 and PgSPL24-09 were most proper two transcripts to further study as ginseng age's molecular marker. These results provide the basis for further elucidation of the PgSPL transcripts' biological function in ginseng and ginseng genetics improvement and gene breeding in the future.

Genome-wide identification and characterization of novel microRNAs in seed development of soybean

MicroRNAs (miRNAs) are important and ubiquitous regulators of gene expression in eukaryotes. However, the information about miRNAs population and their regulatory functions involving in soybean seed development remains incomplete. Base on the Dicer-like1-mediated cleavage signals during miRNA processing could be employed for novel miRNA discovery, a genome-wide search for miRNA candidates involved in seed development was carried out. As a result, 17 novel miRNAs, 14 isoforms of miRNA (isomiRs) and 31 previously validated miRNAs were discovered. These novel miRNAs and isomiRs represented tissue-specific expression and the isomiRs showed significantly higher abundance than that of their miRNA counterparts in different tissues. After target prediction and degradome sequencing data-based validation, 13 novel miRNA-target pairs were further identified. Besides, five targets of 22-nt iso-gma-miR393h were found to be triggered to produce secondary trans-acting siRNA (ta-siRNAs). Summarily, our results could expand the repertoire of miRNAs with potentially important functions in soybean.

Genome-wide characterization and expression analyses of the MYB superfamily genes during developmental stages in Chinese jujube

The MYB transcription factor (TF) superfamily, one of the largest gene superfamilies, regulates a variety of physiological processes in plants. Although many MYB superfamily genes have been identified in plants, the MYB TFs in Chinese jujube (Ziziphus jujuba Mill.) have not been fully identified and characterized. Additionally, the functions of these genes remain unclear. In total, we identified 171 MYB superfamily genes in jujube and divided them into five subfamilies containing 99 genes of the R2R3-MYB subfamily, 58 genes of the MYB-related subfamily, four genes of the R1R2R3-MYB subfamily, one gene of the 4R-MYB subfamily, and nine genes of the atypical MYB subfamily. The 99 R2R3-MYB genes of jujube were divided into 35 groups, C1–C35, and the 58 MYB-related genes were divided into the following groups: the R-R-type, CCA1-like, I-box-binding-like, TBP-like, CPC-like, and Chinese jujube-specific groups. ZjMYB genes in jujube were well supported by additional highly conserved motifs and exon/intron structures. Most R1 repeats of MYB-related proteins comprised the R2 repeat and had highly conserved EED and EEE residue groups in jujube. Three tandem duplicated gene pairs were found on 12 chromosomes in jujube. According to an expression analysis of 126 ZjMYB genes, MYB-related genes played important roles in jujube development and fruit-related biological processes. The total flavonoid content of jujube fruit decreased as ripening progressed. A total of 93 expressed genes were identified in the RNA-sequencing data from jujube fruit, and 56 ZjMYB members presented significant correlations with total flavonoid contents by correlation analysis. Five pairs of paralogous MYB genes within jujube were composed of nine jujube MYB genes. A total of 14 ZjMYB genes had the same homology to the MYB genes of Arabidopsis and peach, indicating that these 14 MYB genes and their orthologs probably existed before the ancestral divergence of the MYB superfamily. We used a synteny analysis of MYB genes between jujube and Arabidopsis to predict that the functions of the ZjMYBs involve flavonoid/phenylpropanoid metabolism, the light signaling pathway, auxin signal transduction, and responses to various abiotic stresses (cold, drought, and salt stresses). Additionally, we speculate that ZjMYB108 is an important TF involved in the flavonoid metabolic pathway. This comprehensive analysis of MYB superfamily genes in jujube lay a solid foundation for future comprehensive analyses of ZjMYB gene functions.

Transcriptomic Analysis of Betula halophila in Response to Salt Stress

Soil salinization is a matter of concern worldwide. It can eventually lead to the desertification of land and severely damage local agricultural production and the ecological environment. Betula halophila is a tree with high salt tolerance, so it is of importance to understand and discover the salt responsive genes of B. halophila for breeding salinity resistant varieties of trees. However, there is no report on the transcriptome in response to salt stress in B. halophila. Using Illumina sequencing platform, approximately 460 M raw reads were generated and assembled into 117,091 unigenes. Among these unigenes, 64,551 unigenes (55.12%) were annotated with gene descriptions, while the other 44.88% were unknown. 168 up-regulated genes and 351 down-regulated genes were identified, respectively. These Differentially Expressed Genes (DEGs) involved in multiple pathways including the Salt Overly Sensitive (SOS) pathway, ion transport and uptake, antioxidant enzyme, ABA signal pathway and so on. The gene ontology (GO) enrichments suggested that the DEGs were mainly involved in a plant-type cell wall organization biological process, cell wall cellular component, and structural constituent of cell wall molecular function. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment showed that the top-four enriched pathways were 'Fatty acid elongation', 'Ribosome', 'Sphingolipid metabolism' and 'Flavonoid biosynthesis'. The expression patterns of sixteen DEGs were analyzed by qRT-PCR to verify the RNA-seq data. Among them, the transcription factor AT-Hook Motif Nuclear Localized gene and dehydrins might play an important role in response to salt stress in B. halophila. Our results provide an important gene resource to breed salt tolerant plants and useful information for further elucidation of the molecular mechanism of salt tolerance in B. halophila.

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.

Molecular Characterization of SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) Gene Family in Betula luminifera

As a major family of plant-specific transcription factors, SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) genes play vital regulatory roles in plant growth, development and stress responses. In this study, 18 SPL genes were identified and cloned from Betula luminifera. Two zinc finger-like structures and a nuclear location signal (NLS) segments were existed in the SBP domains of all BlSPLs. Phylogenetic analysis showed that these genes were clustered into nine groups (group I-IX). The intron/exon structure and motif composition were highly conserved within the same group. 12 of the 18 BlSPLs were experimentally verified as the targets of miR156, and two cleavage sites were detected in these miR156-targeted BlSPL genes. Many putative cis-elements, associated with light, stresses and phytohormones response, were identified in the promoter regions of BlSPLs, suggesting that BlSPL genes are probably involved in important physiological processes and developmental events. Tissue-specific expression analysis showed that miR156-targeted BlSPLs exhibited a more differential expression pattern, while most miR156-nontargeted BlSPLs tended to be constitutively expressed, suggesting the distinct roles of miR156-targeted and nontargeted BlSPLs in development and growth of B. luminifera. Further expression analysis revealed that miR156-targeted BlSPLs were dramatically up-regulated with age, whereas mature BlmiR156 level was apparently declined with age, indicating that miR156/SPL module plays important roles in vegetative phase change of B. luminifera. Moreover, yeast two-hybrid assay indicated that several miR156-targeted and nontargeted BlSPLs could interact with two DELLA proteins (BlRGA and BlRGL), which suggests that certain BlSPLs take part in the GA regulated processes through protein interaction with DELLA proteins. All these results provide an important basis for further exploring the biological functions of BlSPLs in B. luminifera.