Citrus species are among the most economically important fruit crops. Physiological characteristics and molecular mechanisms associated with de-etiolation have been partially revealed. However, little is known about the mechanisms controlling the expression and function of genes associated with photosynthesis and chlorophyll biosynthesis in variegated citrus seedlings. The lower biomass, chlorophyll contents, and photosynthetic parameter values recorded for the variegated seedlings suggested that chlorophyll biosynthesis was partially inhibited. Additionally, roots of the variegated seedlings were longer than the roots of green seedlings. We obtained 567.07 million clean reads and 85.05 Gb of RNA-sequencing data, with more than 94.19% of the reads having a quality score of Q30 (sequencing error rate = 0.1%). Furthermore, we detected 4,786 and 7,007 differentially expressed genes (DEGs) between variegated and green Shiranuhi and Huangguogan seedlings. Thirty common pathways were differentially regulated, including pathways related to photosynthesis (GO: 0015979) and the chloroplast (GO: 0009507). Photosynthesis (44 and 63 DEGs), photosynthesis-antenna proteins (14 and 29 DEGs), and flavonoid biosynthesis (16 and 29 DEGs) pathways were the most common KEGG pathways detected in two analyzed libraries. Differences in the expression patterns of PsbQ, PetF, PetB, PsaA, PsaN, PsbP, PsaF, Cluster-2274.8338 (ZIP1), Cluster-2274.38688 (PTC52), and Cluster-2274.78784 might be responsible for the variegation in citrus seedlings. We completed a physiological- and transcriptome-level comparison of the Shiranuhi and Huangguogan cultivars that differ in terms of seedling variegation. We performed mRNA-seq analyses of variegated and green Shiranuhi and Huangguogan seedlings to explore the genes and regulatory pathways involved in the inhibition of chlorophyll biosynthesis and decreases in Chl a and Chl b contents. The candidate genes described herein should be investigated in greater detail to further characterize variegated citrus seedlings.

Fruit cracking affects both the yield and economic efficiency of citrus; however, the underlying mechanism remains unclear. Therefore, this study focused on resistant and susceptible cultivars to identify the mechanisms underlying fruit cracking. The results showed that in ‘Mingrijian’, pectin morphological transformation and hemicellulose and lignin degradation in the pericarp were important contributing factors. During the critical fruit-cracking period (115–150 days after flowering), the water-soluble pectin, protopectin, and lignin contents in the pericarp of ‘Daya’ presented inverse changes relative to those in ‘Mingrijian’, thus enhancing the mechanical properties and resistance of pericarp. From 115 to 150 days after flowering, the soluble sugar content in the pulp of ‘Mingrijian’ increased rapidly by 97.35%, aiding in pulp water absorption and expansion. Moreover, the soluble protein content in the pericarp of ‘Mingrijian’ exhibited a declining trend and was lower than that of ‘Daya’, thus affecting the overall metabolism. The superoxide dismutase (SOD) activity in the pericarp of ‘Mingrijian’ gradually decreased from 115 to 180 days after flowering, while the peroxidase (POD) activity remained at a low level, resulting in weaker antioxidant capacity and lower environmental resistance. This study provides valuable insights into the mechanisms of citrus fruit cracking, laying the foundation for preventive and control strategies.

The interstock, a stock between the rootstock and scion, has a significant regulatory effect on the stock and scion, and its function is highly dependent on graft compatibility. To assess the graft compatibility of the interstock and scion, ‘Yuanxiaochun’ was top grafted onto ‘Ponkan’, ‘Shiranuhi’, ‘Harumi’, ‘Tarocco’, and ‘Kumquat’. The results showed significant differences in the survival ratio and preservation ratio among different combinations. Grafting compatibility in the early stages of grafting was associated with the genetic relationship. The biomass accumulation revealed that the interstock could influence both the rootstock and the scion. The physiological and biochemical traits analysis suggested that SPAD reading and soluble sugar could be employed as preselected indices to evaluate graft compatibility in the late stage of grafting. These results indicated that the evaluation of graft compatibility was a dynamic process. The findings provided a new approach for studying the stock and scion interaction mechanisms mediated by interstock, and directly provided a theoretical and practical basis for the high-grafting of ‘Yuanxiaochun’ citrus.

Xyloglucan endotransglycosylase (XET) genes are widely distributed in most plants, but the codon usage bias of XET genes has remained uncharacterized. Thus, we analyzed the codon usage bias using 4500 codons of 20 XET genes to elucidate the genetic and evolutionary patterns. Phylogenetic and hierarchical cluster analyses revealed that the 20 XET genes belonged to two groups. The closer the genetic distance, the more similar the codon usage preference. The codon usage bias of most XET genes was weak, but there was also some codon usage bias. AGA, AGG, AUC, and GUG were the top four codons (RSCU > 1.5) in the 20 XET genes. CitXET had a stronger codon usage bias, and there were eight optimal codons of CitXET (i.e., AGA, AUU, UCU, CUU, CCA, GCU, GUU, and AAA). The RSCU values underwent a correspondence analysis. The two main factors affecting codon usage bias (i.e., Axes 1 and 2) accounted for 54.8% and 17.6% of the total variation, respectively. Multiple correspondence analysis revealed that XET genes were widely distributed, with Group 1 genes being closer to Axis 1 than Group 2 genes, which were closer to Axis 2. Codons with A/U at the third codon position were distributed closer to Axis 1 than codons with G/C at the third codon position. PgXET, ZmXET, VlXET, VrXET, and PcXET were biased toward codons ending with G/C. In contrast, CitXET, DpXET, and BrpXET were strongly biased toward codons ending with A/U, indicating that these XET genes have a strong codon usage bias. Translational selection and base composition (especially A and U at the third codon position), followed by mutation pressure and natural selection, may be the most important factors affecting codon usage of 20 XET genes. These results may be useful in clarifying the codon usage bias of XET genes and the relevant evolutionary characteristics.