Phylogenetic analysis of IDD gene family and characterization of its expression in response to flower induction in Malus

Insights into flowering mechanisms in apple (Malus × domestica Borkh.) amidst climate change: An exploration of genetic and epigenetic factors

Apple (Malus × domestica Borkh.) holds a prominent position among global temperate fruit crops, with flowering playing a crucial role in both production and breeding. This review delves into the intricate mechanisms governing apple flowering amidst the backdrop of climate change, acknowledging the profound influence of external and internal factors on biennial bearing, flower bud quality, and ultimately, fruit quality. Notably, the challenge faced in major apple production regions is not an inadequacy of flowers but an excess, leading to compromised fruit quality necessitating thinning practices. Climate change exacerbates these challenges, rendering apple trees more susceptible to crop failure due to unusual weather events, such as reduced winter snowfall, early spring cold weather, and hailstorms during flowering and fruit setting. Altered climatic conditions, exemplified by increased spring warming coupled with sub-freezing temperatures, negatively impact developing flower buds and decrease overall crop production. Furthermore, changing winter conditions affect chilling accumulation, disrupting flower development and synchronicity. Although the physiological perception of apple flowering has been reviewed in the past, the genetic, epigenetic, and multi-omics regulatory mechanisms governing floral induction and flowering are still rarely discussed in the case of apple flowering. This article comprehensively reviews the latest literature encompassing all aspects of apple flowering, aiming to broaden our understanding and address flowering challenges while also laying a solid foundation for future research in developing cultivars that are ideally adapted to climate change.

Comprehensive Analysis of the INDETERMINATE DOMAIN (IDD) Gene Family and Their Response to Abiotic Stress in Zea mays

Transcription factors (TFs) are important regulators of numerous gene expressions due to their ability to recognize and combine cis-elements in the promoters of target genes. The INDETERMINATE DOMAIN (IDD) gene family belongs to a subfamily of C2H2 zinc finger proteins and has been identified only in terrestrial plants. Nevertheless, little study has been reported concerning the genome-wide analysis of the IDD gene family in maize. In total, 22 ZmIDD genes were identified, which can be distributed on 8 chromosomes in maize. On the basis of evolutionary relationships and conserved motif analysis, ZmIDDs were categorized into three clades (1, 2, and 3), each owning 4, 6, and 12 genes, respectively. We analyzed the characteristics of gene structure and found that 3 of the 22 ZmIDD genes do not contain an intron. Cis-element analysis of the ZmIDD promoter showed that most ZmIDD genes possessed at least one ABRE or MBS cis-element, and some ZmIDD genes owned the AuxRR-core, TCA-element, TC-rich repeats, and LTR cis-element. The Ka:Ks ratio of eight segmentally duplicated gene pairs demonstrated that the ZmIDD gene families had undergone a purifying selection. Then, the transcription levels of ZmIDDs were analyzed, and they showed great differences in diverse tissues as well as abiotic stresses. Furthermore, regulatory networks were constructed through the prediction of ZmIDD-targeted genes and miRNAs, which can inhibit the transcription of ZmIDDs. In total, 6 ZmIDDs and 22 miRNAs were discovered, which can target 180 genes and depress the expression of 9 ZmIDDs, respectively. Taken together, the results give us valuable information for studying the function of ZmIDDs involved in plant development and climate resilience in maize.

Genome-Wide Identification of the Highly Conserved INDETERMINATE DOMAIN (IDD) Zinc Finger Gene Family in Moso Bamboo (Phyllostachys edulis)

INDETERMINATE DOMAIN (IDD) proteins, a family of transcription factors unique to plants, function in multiple developmental processes. Although the IDD gene family has been identified in many plants, little is known about it in moso bamboo. In this present study, we identified 32 PheIDD family genes in moso bamboo and randomly sequenced the full-length open reading frames (ORFs) of ten PheIDDs. All PheIDDs shared a highly conserved IDD domain that contained two canonical C2H2-ZFs, two C2HC-ZFs, and a nuclear localization signal. Collinearity analysis showed that segmental duplication events played an important role in expansion of the PheIDD gene family. Synteny analysis indicated that 30 PheIDD genes were orthologous to those of rice (Oryza sativa). Thirty PheIDDs were expressed at low levels, and most PheIDDs exhibited characteristic organ-specific expression patterns. Despite their diverse expression patterns in response to exogenous plant hormones, 8 and 22 PheIDDs responded rapidly to IAA and 6-BA treatments, respectively. The expression levels of 23 PheIDDs were closely related to the outgrowth of aboveground branches and 20 PheIDDs were closely related to the awakening of underground dormant buds. In addition, we found that the PheIDD21 gene generated two products by alternative splicing. Both isoforms interacted with PheDELLA and PheSCL3. Furthermore, both isoforms could bind to the cis-elements of three genes (PH02Gene17121, PH02Gene35441, PH02Gene11386). Taken together, our work provides valuable information for studying the molecular breeding mechanism of lateral organ development in moso bamboo.

Genome-Wide Analysis of Genes Involved in the GA Signal Transduction Pathway in 'duli' Pear (Pyrus betulifolia Bunge)

Gibberellic acid (GA) is an important phytohormone that regulates every aspect of plant growth and development. While elements involved in GA signaling have been identified and, hence, their functions have been well studied in model plants, such as Arabidopsis and rice, very little is known in pear. We, therefore, analyzed the genes related to GA signaling from the recently sequenced genome of the wildtype ‘duli’ pear (Pyrus betulifolia Bunge), a widely used rootstock for grafting in pear cultivation in China due to its vigorous growth and resistance to abiotic and biotic stress. In total, 15 genes were identified, including five GA receptors PbGID1s (GA-INSENSTIVE DWARF 1), six GA negative regulators, PbDELLAs, and four GA positive regulators, PbSLYs. Exogenous application of GA could promote the expression of PbGID1s but inhibit that of PbDELLAs and PbSLYs in tissue culture ‘duli’ pear seedlings. The expression profiles of these genes in field-grown trees under normal growth conditions, as well as in tissue-cultured seedlings treated with auxin (IAA), GA, paclobutrazol (PAC), abscisic acid (ABA), and sodium chloride (NaCl), were also studied, providing further evidence of the involvement of these genes in GA signaling in ‘duli’ pear plants. The preliminary results obtained in this report lay a good foundation for future research into GA signaling pathways in pear. Importantly, the identification and preliminary functional verification of these genes could guide molecular breeding in order to obtain the highly desired dwarf pear rootstocks for high-density plantation to aid easy orchard management and high yielding of pear fruits.

Multiple indeterminate domain (IDD)-DELLA1 complexes participate in gibberellin feedback regulation in peach

Peach encodes 14 INDETERMINATE DOMAIN (IDD) transcription factors. PpIDD4, -12 and -13 mediated PpDELLA1 binding to the PpGA20ox1 promoter. Each of these three PpIDD-DELLA1 complexes activated transcription of PpGA20ox1. PpTPR1 and -4 interrupted the interaction of PpIDDs with PpDELLA1. The plant growth regulator gibberellin (GA) plays an important role in the rapid growth of annual shoots in peach. Our previous study showed that the peach cultivar 'FenHuaShouXingTao' (FHSXT), a gibberellic acid receptor (gid1) mutant, accumulates active GAs in annual shoot tips. This mutant enhances GA feedback regulation in peach. The results of this study suggested that the PpIDD-DELLA1 complex is the underlying mechanism of GA feedback regulation in peach. Fourteen IDD genes were identified in peach, and three PpIDDs (PpIDD4, -12 and -13, all from group IV) interacted with PpDELLA1, an important component in GA signaling pathway. Truncation, segmentation and site mutation of the promoter of PpGA20ox1 (a GA biosynthesis gene) showed that all three PpIDD proteins recognized the core motif TTGTC. PpIDD4 and -13 mainly bind to site 3, while PpIDD12 binds to site 5 of the PpGA20ox1 promoter. All three PpIDD-DELLA1 complexes activated the PpGA20ox1 promoter-LUC fusion. These data suggested that PpIDDs bridge PpDELLA1 and the promoter of PpGA20ox1, which then activated the transcription of PpGA20ox1. In addition, PpTPR1 and -4 disrupted the interaction of PpIDDs with PpDELLA1. Our research will be helpful for understanding and possibly modifying the regulation of annual shoot growth and GA biosynthesis.

Identification of Characteristic Genes in Whole Blood of Intervertebral Disc Degeneration Patients by Weighted Gene Coexpression Network Analysis (WGCNA)

Intervertebral disc degeneration (IDD) is a major cause of lower back pain. However, to date, the molecular mechanism of the IDD remains unclear. Gene expression profiles and clinical traits were downloaded from the Gene Expression Omnibus (GEO) database. Firstly, weighted gene coexpression network analysis (WGCNA) was used to screen IDD-related genes. Moreover, least absolute shrinkage and selection operator (LASSO) logistic regression and support vector machine (SVM) algorithms were used to identify characteristic genes. Furthermore, we further investigated the immune landscape by the Cell-type Identification By Estimating Relative Subsets Of RNA Transcripts (CIBERSORT) algorithm and the correlations between key characteristic genes and infiltrating immune cells. Finally, a competing endogenous RNA (ceRNA) network was established to show the regulatory mechanisms of characteristic genes. A total of 2458 genes were identified by WGCNA, and 48 of them were disordered. After overlapping the genes obtained by LASSO and SVM-RFE algorithms, genes including LINC01347, ASAP1-IT1, lnc-SEPT7L-1, B3GNT8, CHRNB3, CLEC4F, LOC102724000, SERINC2, and LOC102723649 were identified as characteristic genes of IDD. Moreover, differential analysis further identified ASAP1-IT1 and SERINC2 as key characteristic genes. Furthermore, we found that the expression of both ASAP1-IT1 and SERINC2 was related to the proportions of T cells gamma delta and Neutrophils. Finally, a ceRNA network was established to show the regulatory mechanisms of ASAP1-IT1 and SERINC2. In conclusion, the present study identified ASAP1-IT1 and SERINC2 as the key characteristic genes of IDD through integrative bioinformatic analyses, which may contribute to the diagnosis and treatment of IDD.

Mdm-MIR393b-mediated adventitious root formation by targeted regulation of MdTIR1A expression and weakened sensitivity to auxin in apple rootstock

Adventitious root (AR) formation is of great significance for apple rootstock breeding. It is widely accepted that miR393 influences AR formation in many plant species; however, the molecular mechanism by which factors regulate AR formation remains insufficient. In this study, the evolutionary relationship of mdm-miR393 and candidate target genes MdTIR1/AFB was systematically identified, and the expression patterns were analysed. Multisequence alignment analysis of miR393 family members suggests that miR393 conservatively evolved between different species. The evolutionary relationship of the TIR1/AFBs can be divided into G1, G2 and G3 subgroups. During AR formation, the expression level of mdm-miR393a/b/c was significantly upregulated at 1 d and 7 d by exogenous auxin treatment. Furthermore, the expression levels of MdTIR1A, MdTIR1D, MdAFB1, MdAFB2, MdAFB3, MdAFB4 and MdAFB8 also appeared to be significantly changed by exogenous auxin induction. Subsequently, tissue-specific expression analysis showed that the expression levels of mdm-miR393 and MdTIR1/AFBs in different tissues exhibited significant differences. The promoter of mdm-miR393 contains multiple elements that respond to ABA, adversity and light signals; auxin treatment can activate the mdm-MIR393b promoter but is obviously inhibited by NPA treatment. The targeting relationship between mdm-MIR393b and MdTIR1A was verified by expression patterns, degradation group data, transient tobacco conversion results, and genes functions experiments. Heterologous overexpression of mdm-MIR393b (35S::mdm-MIR393b) decreased the number of ARs in the phenotype and reduced the expression level of the target gene NtTIR1 in tobacco. Compared to the wild type, the 35S::mdm-MIR393b transgenic plants demonstrated insensitivity to auxin. Furthermore, tir1 mutant exhibited reduced root system structure relative to the control. The above results illustrated that mdm-MIR393b is involved in mediating AR formation by targeted regulation of MdTIR1A expression in apple rootstock.

Genomic identification and expression analysis of nuclear pore proteins in Malus domestica

The nuclear pore complex (NPC), comprised of individual nucleoporin (Nup) proteins, controls nucleo-cytoplasmic transport of RNA and protein, and is important for regulating plant growth and development. However, there are no reports on this complex in fruit tree species. In this study, we identified 38 apple Nups and named them based on the known Arabidopsis thaliana homologs. We also completed bioinformatics analyses of the intron and exon structural data for apple Nups. The proteins encoded by the apple Nups lacked a universally conserved domain. Moreover, a phylogenetic analysis separated the apple and A. thaliana Nups into three groups. The phylogenetic tree indicated that MdNup54 and MdNup62 are most closely related to genes in other Rosaceae species. To characterize the 38 candidate Malus domestica Nups, we measured their stage-specific expression levels. Our tests revealed these proteins were differentially expressed among diverse tissues. We analyzed the expression levels of seven apple Nups in response to an indole-3-acetic acid (IAA) treatment. The phytohormone treatment significantly inhibited apple flowering. A qRT-PCR analysis proved that an IAA treatment significantly inhibited the expression of these seven genes. A preliminary study regarding two members of the Nup62 subcomplex, MdNup54 and MdNup62, confirmed these two proteins can interact with each other. A yeast two-hybrid assay verified that MdNup54 can interact with MdKNAT4 and MdKNAT6. On the basis of the study results, we identified apple NPC and predicted its structure and function. The data generated in this investigation provide important reference material for follow-up research.

Step by step evolution of Indeterminate Domain (IDD) transcriptional regulators: from algae to angiosperms

INTRODUCTION

The Indeterminate Domain (IDD) proteins are a plant-specific subclass of C2H2 Zinc Finger transcription factors. Some of these transcription factors play roles in diverse aspects of plant metabolism and development, but the function of most of IDD genes is unknown and the molecular evolution of the subfamily has not been explored in detail.

METHODS

In this study, we mined available genome sequences of green plants (Viridiplantae) to reconstruct the phylogeny and then described the motifs/expression patterns of IDD genes.

KEY RESULTS

We identified the complete set of IDD genes of 16 Streptophyta genomes. We found that IDD and its sister clade STOP arose by a duplication at the base of Streptophyta. Once on land, the IDD genes duplicated extensively, giving rise to at least ten lineages. Some of these lineages were lost in extant non-vascular plants and gymnosperms, but all of them were retained in angiosperms, duplicating profoundly in dicots and monocots and acquiring, at the same time, surprising heterogeneity in their C-terminal regions and expression patterns.

CONCLUSIONS

IDDs were present in the last common ancestor of Streptophyta. On land, IDDs duplicated extensively, leading to ten lineages. Later, IDDs were recruited by angiosperms where they diversified greatly in number, C-terminal and expression patterns. Interestingly, such diversification occurred during the evolution of novel traits of the plant body. This study provides a solid framework of the orthology relationships of green land plant IDD transcription factors, thus increasing the accuracy of orthologue identification in model and non-model species and facilitating the identification of agronomically important genes related to plant metabolism and development.

Genome-Wide Identification of the MdKNOX Gene Family and Characterization of Its Transcriptional Regulation in Malus domestica

Knotted1-like Homeobox (KNOX) proteins play important roles in regulating plant growth, development, and other biological processes. However, little information is available on the KNOX gene family in apple (Malus domestica Borkh.). In this study, 22 KNOX genes were identified in the apple genome. The gene structure, protein characteristics, and promoter region were characterized. The MdKNOX family members were divided into three classes based on their phylogenetic relationships. Quantitative real-time PCR analysis revealed that the majority of MdKNOX genes exhibited strongly preferential expression in buds and were significantly up-regulated during the flower induction period. The transcript levels of MdKNOX genes were responsive to treatments with flowering- and stress-related hormones. The putative upstream regulation factor MdGRF could directly bind to the promoter of MdKNOX15 and MdKNOX19, and inhibit their transcriptional activities, which were confirmed by yeast one-hybrid and dual-luciferase assays. The results provide an important foundation for future analysis of the regulation and functions of the MdKNOX gene family.

Comparative genomic analysis of the IDD genes in five Rosaceae species and expression analysis in Chinese white pear (Pyrus bretschneideri)

The INDETERMINATE DOMAIN (IDD) gene family encodes hybrid transcription factors with distinct zinc finger motifs and appears to be found in all higher plant genomes. IDD genes have been identified throughout the genomes of the model plants Arabidopsis thaliana and Oryza sativa, and the functions of many members of this gene family have been studied. However, few studies have investigated the IDD gene family in Rosaceae species (among these species, a genome-wide identification of the IDD gene family has only been completed in Malus domestica). This study focuses on a comparative genomic analysis of the IDD gene family in five Rosaceae species (Pyrus bretschneideri, Fragaria vesca, Prunus mume, Rubus occidentalis and Prunus avium). We identified a total of 68 IDD genes: 16 genes in Chinese white pear, 14 genes in F. vesca, 13 genes in Prunus mume, 14 genes in R. occidentalis and 11 genes in Prunus avium. The evolution of the IDD genes in these five Rosaceae species was revealed by constructing a phylogenetic tree, tracking gene duplication events, and performing a sliding window analysis and a conserved microsynteny analysis. The expression analysis of different organs showed that most of the pear IDD genes are found at a very high transcription level in fruits, flowers and buds. Based on our results with those obtained in previous research, we speculated that PbIDD2 and PbIDD8 might participate in flowering induction in pear. A temporal expression analysis showed that the expression patterns of PbIDD3 and PbIDD5 were completely opposite to the accumulation pattern of fruit lignin and the stone cell content. The results of the composite phylogenetic tree and expression pattern analysis indicated that PbIDD3 and PbIDD5 might be involved in the metabolism of lignin and secondary cell wall (SCW) formation. In summary, we provide basic information about the IDD genes in five Rosaceae species and thereby provide a theoretical basis for studying the function of these IDD genes.

Mediation of Flower Induction by Gibberellin and its Inhibitor Paclobutrazol: mRNA and miRNA Integration Comprises Complex Regulatory Cross-Talk in Apple

Guaranteeing successful flowering is very important in economic plant species, especially apple (Malus domestica Borkh.), which is difficult to induce to flower. However, the gene expression and networks involved in flowering have not been totally characterized. Here, we employed mRNA and microRNA (miRNA) sequencing to understand the different responses to gibberellin- and its inhibitor paclobutrazol- (PAC) mediated flower induction. Significant opposite cytological and morphological changes were observed in treated terminal buds, which led to a reduced flowering rate under gibberellin and an increased flowering rate under PAC. We also found that the differentially expressed mRNAs, miRNAs and miRNA target genes participated in different biological networks including hormones, photosynthesis, redox state and other metabolic processes, which provided important clues to understand the complex networks involved in apple flower induction. Additionally, we subsequently focused on one important candidate, MdSPL3, which is one of 31 apple SPL gene family members and whose transcription was inhibited by gibberellin but promoted by PAC. Functional investigation showed that MdSPL3 was located in the nucleus, and ectopic MdSPL3 activated floral meristem identity genes, promoted the formation of floral primordia and led to an earlier flowering phenotype in Arabidopsis. Our research identified critical mRNA and miRNA responsive to gibberellin or PAC, and provided a candidate framework for flower induction. This carefully orchestrated regulatory cross-talk highlighted potential targets for developing regulatory techniques and genetic improvement of flower induction in apple.

Identification and characterization of histone modification gene family reveal their critical responses to flower induction in apple

BACKGROUND

Histone methylation and acetylation regulate biological processes in plants through various histone modifications (HMs) gene families. However, knowledge of HMs genes is limited in horticultural deciduous trees, including apple (Malus domestica).

RESULTS

Here, a comprehensive study of identifying and investigating HMs genes was performed using the recently published apple genome. In total, 198 MdHMs were identified, including 71 histone methyltransferases, 44 histone demethylases, 57 histone acetylases, and 26 histone deacetylases. Detailed analysis of the MdHMs, including chromosomes locations, gene structures, protein motif and protein-protein interactions were performed, and their orthologous genes were also predicted against nine plant species. Meanwhile, a syntenic analysis revealed that tandem, segmental, and whole genome duplications were involved in the evolution and expansion of the MdHMs gene family. Most MdHMs underwent purifying selection. The expression profiles of 198 MdHMs were investigated in response to 6-BA treatment and different flowering varieties (easy-flowering 'Yanfu No.6' and difficult-flowering 'Nagafu No.2') using transcriptome sequencing data, and most MdHMs were involved in flower induction processes. Subsequent quantitative real-time PCR was then performed to confirm the expression levels of candidate MdHMs under different flowering-related circumstances.

CONCLUSION

MdHMs were involved in, and responsive to, flower induction in apple. This study established an MdHMs platform that provided valuable information and presented enriched biological theories on flower induction in apple. The data could also be used to study the evolutionary history and functional prospects of MdHMs genes, as well as other trees.

Identification and characterization of histone modification gene family reveal their critical responses to flower induction in apple

BACKGROUND

Histone methylation and acetylation regulate biological processes in plants through various histone modifications (HMs) gene families. However, knowledge of HMs genes is limited in horticultural deciduous trees, including apple (Malus domestica).

RESULTS

Here, a comprehensive study of identifying and investigating HMs genes was performed using the recently published apple genome. In total, 198 MdHMs were identified, including 71 histone methyltransferases, 44 histone demethylases, 57 histone acetylases, and 26 histone deacetylases. Detailed analysis of the MdHMs, including chromosomes locations, gene structures, protein motif and protein-protein interactions were performed, and their orthologous genes were also predicted against nine plant species. Meanwhile, a syntenic analysis revealed that tandem, segmental, and whole genome duplications were involved in the evolution and expansion of the MdHMs gene family. Most MdHMs underwent purifying selection. The expression profiles of 198 MdHMs were investigated in response to 6-BA treatment and different flowering varieties (easy-flowering 'Yanfu No.6' and difficult-flowering 'Nagafu No.2') using transcriptome sequencing data, and most MdHMs were involved in flower induction processes. Subsequent quantitative real-time PCR was then performed to confirm the expression levels of candidate MdHMs under different flowering-related circumstances.

CONCLUSION

MdHMs were involved in, and responsive to, flower induction in apple. This study established an MdHMs platform that provided valuable information and presented enriched biological theories on flower induction in apple. The data could also be used to study the evolutionary history and functional prospects of MdHMs genes, as well as other trees.

Genome Identification of B-BOX Gene Family Members in Seven Rosaceae Species and Their Expression Analysis in Response to Flower Induction in Malus domestica

BBX proteins play important roles in regulating plant growth and development including photomorphogenesis, photoperiodic regulation of flowering, and responses to biotic and abiotic stresses. At present, the genomes of seven Rosaceae fruit species have been fully sequenced. However, little is known about the BBX gene family and their evolutionary history in these Rosaceae species. Therefore, in this study total, 212 BBX genes were investigated from seven Rosaceae species (67 from Malus × domestica, 40 from Pyruscommunis, 22 from Rosa Chinesis, 20 from Prunuspersica, 21 from Fragariavesca, 22 from Prunusavium, and 20 from Rubusoccidentalis). The chemical properties, gene structures, and evolutionary relationships of the BBX genes were also studied. All the BBX genes were grouped into six subfamilies on the basis of their phylogenetic relationships and structural features. Analysis of gene structure, segmental and tandem duplication, gene phylogeny, and tissue-specific expression with the ArrayExpress database showed their diversification in function, quantity, and structure. The expression profiles of 19 MdBBX genes in different tissues were evaluated through qRT-PCR. These genes showed distinct transcription level among the tested tissues (bud, flower, fruit, stem, and leaf). Moreover, expression patterns of 19 MdBBX genes were examined during flowering induction time under flowering-related hormones and treatments (GA3, 6-BA, and sucrose). The expressions of the candidates BBX genes were affected and showed diverse expression profile. Furthermore, changes in response to these flowering-related hormones and treatment specifying their potential involvement in flowering induction. Based on these findings, BBX genes could be used as potential genetic markers for the growth and development of plants particularly in the area of functional analysis, and their involvement in flower induction in fruit plants.

Comprehensive analysis of GASA family members in the Malus domestica genome: identification, characterization, and their expressions in response to apple flower induction

BACKGROUND

The plant-specific gibberellic acid stimulated Arabidopsis (GASA) gene family is critical for plant development. However, little is known about these genes, particularly in fruit tree species.

RESULTS

We identified 15 putative Arabidopsis thaliana GASA (AtGASA) and 26 apple GASA (MdGASA) genes. The identified genes were then characterized (e.g., chromosomal location, structure, and evolutionary relationships). All of the identified A. thaliana and apple GASA proteins included a conserved GASA domain and exhibited similar characteristics. Specifically, the MdGASA expression levels in various tissues and organs were analyzed based on an online gene expression profile and by qRT-PCR. These genes were more highly expressed in the leaves, buds, and fruits compared with the seeds, roots, and seedlings. MdGASA genes were also responsive to gibberellic acid (GA₃) and abscisic acid treatments. Additionally, transcriptome sequencing results revealed seven potential flowering-related MdGASA genes. We analyzed the expression levels of these genes in response to flowering-related treatments (GA₃, 6-benzylaminopurine, and sugar) and in apple varieties that differed in terms of flowering ('Nagafu No. 2' and 'Yanfu No. 6') during the flower induction period. These candidate MdGASA genes exhibited diverse expression patterns. The expression levels of six MdGASA genes were inhibited by GA₃, while the expression of one gene was up-regulated. Additionally, there were expression-level differences induced by the 6-benzylaminopurine and sugar treatments during the flower induction stage, as well as in the different flowering varieties.

CONCLUSION

This study represents the first comprehensive investigation of the A. thaliana and apple GASA gene families. Our data may provide useful clues for future studies and may support the hypotheses regarding the role of GASA proteins during the flower induction stage in fruit tree species.

Identification, Classification, and Expression Analysis of GRAS Gene Family in Malus domestica

GRAS genes encode plant-specific transcription factors that play important roles in plant growth and development. However, little is known about the GRAS gene family in apple. In this study, 127 GRAS genes were identified in the apple (Malus domestica Borkh.) genome and named MdGRAS1 to MdGRAS127 according to their chromosomal locations. The chemical characteristics, gene structures and evolutionary relationships of the MdGRAS genes were investigated. The 127 MdGRAS genes could be grouped into eight subfamilies based on their structural features and phylogenetic relationships. Further analysis of gene structures, segmental and tandem duplication, gene phylogeny and tissue-specific expression with ArrayExpress database indicated their diversification in quantity, structure and function. We further examined the expression pattern of MdGRAS genes during apple flower induction with transcriptome sequencing. Eight higher MdGRAS (MdGRAS6, 26, 28, 44, 53, 64, 107, and 122) genes were surfaced. Further quantitative reverse transcription PCR indicated that the candidate eight genes showed distinct expression patterns among different tissues (leaves, stems, flowers, buds, and fruits). The transcription levels of eight genes were also investigated with various flowering related treatments (GA₃, 6-BA, and sucrose) and different flowering varieties (Yanfu No. 6 and Nagafu No. 2). They all were affected by flowering-related circumstance and showed different expression level. Changes in response to these hormone or sugar related treatments indicated their potential involvement during apple flower induction. Taken together, our results provide rich resources for studying GRAS genes and their potential clues in genetic improvement of apple flowering, which enriches biological theories of GRAS genes in apple and their involvement in flower induction of fruit trees.