WUSCHEL-related homeobox (WOX) transcription factors (TFs) are well known for their role in plant development but are rarely studied in citrus. In this study, we identified 11 putative genes from the sweet orange genome and divided the citrus WOX genes into three clades (modern/WUSCHEL(WUS), intermediate, and ancient). Subsequently, we performed syntenic relationship, intron-exon organization, motif composition, and cis-element analysis. Co-expression analysis based on RNA-seq and tissue-specific expression patterns revealed that CsWOX gene expression has multiple intrinsic functions. CsWUS homolog of AtWUS functions as a transcriptional activator and binds to specific DNA. Overexpression of CsWUS in tobacco revealed dramatic phenotypic changes, including malformed leaves and reduced gynoecia with no seed development. Silencing of CsWUS in lemon using the virus-induced gene silencing (VIGS) system implied the involvement of CsWUS in cells of the plant stem. In addition, CsWUS was found to interact with CsCYCD3, an ortholog in Arabidopsis (AtCYCD3,1). Yeast one-hybrid screening and dual luciferase activity revealed that two TFs (CsRAP2.12 and CsHB22) bind to the promoter of CsWUS and regulate its expression. Altogether, these results extend our knowledge of the WOX gene family along with CsWUS function and provide valuable findings for future study on development regulation and comprehensive data of WOX members in citrus.

SQUAMOSA-promoter binding protein (SBP)-box genes encode a family of plant-specific transcription factors that play vital roles in plant growth and development. In this study, 15 SBP-box genes were identified and isolated from Citrus clementina (CclSBPs), where 10 of these genes were predicted to be putative targets of Citrus clementina microRNA156 (CclmiR156). The 15 CclSBP genes could be classified into six groups based on phylogenetic analysis, diverse intron–exon structure, and motif prediction, similar to the SQUAMOSA promoter binding protein-like (SPL) gene family of Populus trichocarpa and Arabidopsis thaliana. Furthermore, CclSBPs classified into a group/subgroup have similar gene structures and conserved motifs, implying their functional redundancy. Tissue-specific expression analysis of CclSBPs demonstrated their diversified expression patterns. To further explore the potential role of CclSBPs during floral inductive water deficits, the dynamic changes of the 15 CclSBPs were investigated during floral inductive water deficits, and the results showed that some CclSBPs were associated with floral induction. Among these genes, CclSBP6 was not homologous to the Arabidopsis SBP-box gene family, and CclSBP7 was regulated by being alternatively spliced. Therefore, CclSBP6 and CclSBP7 were genetically transformed in Arabidopsis. Overexpression of the two genes changed the flowering time of Arabidopsis.

A high-density genetic linkage map is essential for genetic and genomic studies including QTL mapping, genome assembly, and comparative genomic analysis. Here, we constructed a citrus high-density linkage map using SSR and SNP markers, which are evenly distributed across the citrus genome. The integrated linkage map contains 4163 markers with an average distance of 1.12 cM. The female and male linkage maps contain 1478 and 2976 markers with genetic lengths of 1093.90 cM and 1227.03 cM, respectively. Meanwhile, a genetic map comparison demonstrates that the linear order of common markers is highly conserved between the clementine mandarin and Poncirus trifoliata. Based on this high-density integrated citrus genetic map and two years of deciduous phenotypic data, two loci conferring leaf abscission phenotypic variation were detected on scaffold 1 (including 36 genes) and scaffold 8 (including 107 genes) using association analysis. Moreover, the expression patterns of 30 candidate genes were investigated under cold stress conditions because cold temperature is closely linked with the deciduous trait. The developed high-density genetic map will facilitate QTL mapping and genomic studies, and the localization of the leaf abscission deciduous trait will be valuable for understanding the mechanism of this deciduous trait and citrus breeding.

To identify the genetic variation between Citrus trifoliata and Citrus clementina, we performed genome resequencing on the two citrus species. Compared with the citrus reference genome, a total of 9,449,204 single-nucleotide polymorphisms (SNPs) and 846,615 insertion/deletion polymorphisms (InDels) were identified in the two citrus species, while 1,868,115 (19.77%) of the SNPs and 190,199 (22.47%) of the InDels from the two citrus species were located in the genic regions. Meanwhile, a total of 8,091,407 specific SNPs and 692,654 specific InDels were identified in the two citrus genotypes, yielding an average of 27.32 SNPs/kb and 2.34 InDels/kb. We identified and characterized the patterns of gene exchanges in the grafted citrus plants by using specific genetic variation from genome resequencing. A total of 4396 transporting genes across graft junctions was identified. Some specific genetic variation and mobile genes was also confirmed by Sanger sequencing. Furthermore, these mobile genes could move directionally or bidirectionally between the scions and the rootstocks. In addition, a total of 1581 and 2577 differentially expressed genes were found in the scions and the rootstocks after grafting compared with the control, respectively. These genetic variations provide fundamental information on the genetic basis of important traits between C. trifoliata and C. clementina, as the transport of genes would be applicable to horticulture crops.