Cyclocarya paliurus is an important medical plant owing to the diverse bioactive compounds in its leaves. However, the heterodichogamy with female and male functions segregation within protandry (PA) or protogyny (PG) may greatly affect seed quality and its plantations for medicinal use. To speculate on the factor playing the dominant role in regulating heterodichogamy in C. paliurus, based on phenotypic observations, our study performed a multi comparison transcriptome analysis on female and male buds (PG and PA types) using RNA-seq. For the female and male bud comparisons, a total of 6753 differentially expressed genes (DEGs) were detected. In addition, functional analysis revealed that these DEGs were significantly enriched in floral development, hormone, and GA-related pathways. As the dominant hormones responsible for floral differentiation and development, gibberellins (GAs) in floral buds from PG and PA types were quantified using HPLC-MS. Among the tested GAs, GA₃ positively regulated the physiological differentiation (S0) and germination (S2) of floral buds. The dynamic changes of GA₃ content and floral morphological features were consistent with the expression levels of GA-related genes. Divergences of GA₃ contents at S0 triggered the asynchronism of physiological differentiation between male and female buds of intramorphs (PA-M vs. PA-F and PG-F vs. PG-M). A significant difference in GA₃ content enlarged this asynchronism at S2. Thus, we speculate that GA₃ plays the dominant role in the formation of heterodichogamy in C. paliurus. Meanwhile, the expression patterns of GA-related DEGs, including CPS, KO, GA20ox, GA2_OX, GID1, and DELLA genes, which play central roles in regulating flower development, coincided with heterodichogamous characteristics. These results support our speculations well, which should be further confirmed.

The GRAS gene family, derived from GAI, RGA, and SCR, plays a crucial role in plant growth and development. In the diploid Cyclocarya paliurus (2n = 2x = 32) with heterodichogamous characteristics, 51 CpGRAS genes were identified and phylogenetically classified into 10 subfamilies. Structural analysis revealed that CpGRAS genes possessed a canonical GRAS domain, but 70% lacked introns. WGD/segmental duplication was the major driver in the expansion of the CpGRAS family. In addition, a Ka/Ks ratio below 1 for these genes implied functional constraints and evolutionary conservation. Transcriptional profiling revealed significant differential expressions of CpGRAS genes between male and female flowers from two mating types, protogyny (PG) and protandry (PA). Notably, members of the DELLA subfamily exhibited significant upregulation in female flowers at the inflorescence elongation (S3) stage. The expression level of CpSCL6-2 in late-flowering samples (PA-F and PG-M) was higher than in early-flowering ones (PA-M and PG-F). Co-expression analysis identified that CpRGL1 and CpGAI-2 of the DELLA subfamily, along with CpSCL6-2, acted as hub genes, implying their crucial roles in floral development and potential involvement in the heterodichogamous flowering mechanism in C. paliurus. These findings broaden our understanding of CpGRAS genes and provide new insights into the molecular basis of heterodichogamy.