The bulbil is an important vegetative reproductive organ in triploid Lilium lancifolium whose development is promoted by cytokinins. Type-B response regulators (RRs) are critical regulators that mediate primary cytokinin responses and promote cytokinin-induced gene expression. However, the function of cytokinin type-B Arabidopsis RRs (ARRs) in regulating bulbil formation is unclear. In this study, we identified five type-B LlRRs, LlRR1, LlRR2, LlRR10, LlRR11 and LlRR12, in L. lancifolium for the first time. The five LlRRs encode proteins of 715, 675, 573, 582 and 647 amino acids. All of the regulators belong to the B-I subfamily, whose members typically contain a conserved CheY-homologous receiver (REC) domain and an Myb DNA-binding (MYB) domain at the N-terminus. As transcription factors, all five type-B LlRRs localize at the nucleus and are widely expressed in plant tissues, especially during axillary meristem (AM) formation. Functional analysis showed that type-B LlRRs are involved in bulbil formation in a functionally redundant manner and can activate LlRR9 expression. In summary, our study elucidates the process by which cytokinins regulate bulbil initiation in L. lancifolium through type-B LlRRs and lays a foundation for research on the molecular mechanism of bulbil formation in the lily.

The biosynthesis and accumulation of plant secondary metabolites are tightly regulated by environmental fluctuations, serving as a crucial interface mediating plant–environment interactions. Nevertheless, the phenotypic instability of secondary metabolism-related traits induced by environmental variability has hampered the precise breeding of stress-resistant cultivars. Pinus taeda is an key timber tree species in southern China, and its foliar catechins exhibit substantial stress-resistant potential. However, phenotypic variation driven by seasonal changes has limited the germplasm innovation and genetic selection of this species. In this study, 54 half-sib families of P. taeda were used as experimental materials. Combined with near-infrared spectroscopy (NIRS) and the BLUP model, we systematically analyzed the seasonal variation characteristics, genetic parameters of catechin content (CC), and genetic gains under different breeding strategies across four seasons. Our results demonstrated that family and season had extremely significant effects on CC (p < 0.001), whereas the season × family interaction effect was not significant, indicating that the genetic expression of CC is stable across seasons. CC was higher in spring and winter but lower in summer and autumn; specifically, the mean CC in summer was 47% lower than the peak value in spring (26.95 ± 0.46 μg·g⁻¹), reflecting a resource trade-off between growth and defense metabolism. Genetic parameter analysis revealed that family-mean heritability (0.373–0.714) was higher than individual heritability and within-family heritability, with August identified as the optimal selection season. The maximum genetic gain across the three breeding strategies (individual selection, family selection, and combined selection) reached 7.86%, among which individual selection exhibited the smallest fluctuation in genetic gain. Finally, three superior families and 14 superior individuals were screened out. This study elucidates the seasonal genetic pattern of foliar CC in P. taeda, clarifies the optimal selection stage and efficient breeding strategies, and provides theoretical guidance and material support for the genetic improvement, germplasm innovation, and resource utilization of secondary metabolic traits in this ecologically and economically important tree species.