Genetically based location from triploid populations and gene ontology of a 3.3-mb genome region linked to Alternaria brown spot resistance in citrus reveal clusters of resistance genes

Genomic insights into citrus domestication and its important agronomic traits

Citrus originated in Southeast Asia, and it has become one of the most important fruit crops worldwide. Citrus has a long and obscure domestication history due to its clonal propagation, long life cycle, wide sexual compatibility, and complex genetic background. As the genomic information of both wild and cultivated citrus becomes available, their domestication history and underlying traits or genes are becoming clear. This review outlines the genomic features of wild and cultivated species. We propose that the reduction of citric acid is a critical trait for citrus domestication. The genetic model representing the change during domestication may be associated with a regulatory complex known as WD-repeat-MYB-bHLH-WRKY (WMBW), which is involved in acidification and anthocyanin accumulation. The reduction in or loss of anthocyanins may be due to a hitchhiking effect of fruit acidity selection, in which mutation occurs in the common regulator of these two pathways in some domesticated types. Moreover, we have summarized the domestication traits and candidate genes for breeding purposes. This review represents a comprehensive summary of the genes controlling key traits of interest, such as acidity, metabolism, and disease resistance. It also sheds light on recent advances in early flowering from transgenic studies and provides a new perspective for fast breeding of citrus. Our review lays a foundation for future research on fruit acidity, flavor, and disease resistance in citrus.

Comparative proteomics of two citrus varieties in response to infection by the fungus Alternaria alternata

Alternaria brown spot (ABS) is a disease caused by the necrotrophic fungus Alternaria alternata, which induces necrotic lesions on fruits and young leaves due to the production of the host-specific ACT toxin by the fungus. To better understand the citrus-A. alternata interaction and to identify putative resistance proteins, as well as the receptor of the ACT toxin, citrus plants susceptible ('Minneola' mandarin) and resistant ('Clemenules' tangor) to A. alternata, infected or not (control) with the pathogen were analyzed by proteomics. Protein changes were observed between citrus genotypes after infection, and 150 candidate proteins were obtained. A general scheme of the metabolic processes involved in susceptible and resistant citrus-A. alternata interactions was designed. Susceptible plants presented a high level of proteins involved in stress response at the final stages of the infection, whereas resistant plants presented high level of ROS proteins, metabolic proteins, and proteins involved in the immune system process. Proteins like ferredoxin and cyclophilin are specific to the susceptible variety and may be good candidates as fungal effector-interacting proteins. This is the first citrus-A. alternata proteomics analysis, which has allowed a better understanding of the molecular bases of the citrus response to ABS disease.

UAV-Based High Throughput Phenotyping in Citrus Utilizing Multispectral Imaging and Artificial Intelligence

Traditional plant breeding evaluation methods are time-consuming, labor-intensive, and costly. Accurate and rapid phenotypic trait data acquisition and analysis can improve genomic selection and accelerate cultivar development. In this work, a technique for data acquisition and image processing was developed utilizing small unmanned aerial vehicles (UAVs), multispectral imaging, and deep learning convolutional neural networks to evaluate phenotypic characteristics on citrus crops. This low-cost and automated high-throughput phenotyping technique utilizes artificial intelligence (AI) and machine learning (ML) to: (i) detect, count, and geolocate trees and tree gaps; (ii) categorize trees based on their canopy size; (iii) develop individual tree health indices; and (iv) evaluate citrus varieties and rootstocks. The proposed remote sensing technique was able to detect and count citrus trees in a grove of 4,931 trees, with precision and recall of 99.9% and 99.7%, respectively, estimate their canopy size with overall accuracy of 85.5%, and detect, count, and geolocate tree gaps with a precision and recall of 100% and 94.6%, respectively. This UAV-based technique provides a consistent, more direct, cost-effective, and rapid method to evaluate phenotypic characteristics of citrus varieties and rootstocks.

Construction of High-Density Genetic Maps and Detection of QTLs Associated With Huanglongbing Tolerance in Citrus

Huanglongbing (HLB), or citrus greening, is the most devastating disease in citrus worldwide. Commercial citrus varieties including sweet orange (Citrus sinensis) are highly susceptible to HLB, and trifoliate orange (Poncirus trifoliata, a close Citrus relative) is widely considered resistant or highly tolerant to HLB. In this study, an intergeneric F1 population of sweet orange and trifoliate orange was genotyped by Genotyping-by-Sequencing, and high-density SNP-based genetic maps were constructed separately for trifoliate orange and sweet orange. The two genetic maps exhibited high synteny and high coverage of the citrus genome. Progenies of the F1 population and their parents were planted in a replicated field trial, exposed to intense HLB pressure for 3 years, and then evaluated for susceptibility to HLB over 2 years. The F1 population exhibited a wide range in severity of HLB foliar symptom and canopy damage. Genome-wide QTL analysis based on the phenotypic data of foliar symptom and canopy damage in 2 years identified three clusters of repeatable QTLs in trifoliate orange linkage groups LG-t6, LG-t8 and LG-t9. Co-localization of QTLs for two traits was observed within all three regions. Additionally, one cluster of QTLs in sweet orange (linkage group LG-s7) was also detected. The majority of the identified QTLs each explained 18-30% of the phenotypic variation, indicating their major role in determining HLB responses. These results show, for the first time, a quantitative genetic nature yet the presence of major loci for the HLB tolerance in trifoliate orange. The results suggest that sweet orange also contains useful genetic factor(s) for improving HLB tolerance in commercial citrus varieties. Findings from this study should be very valuable and timely to researchers worldwide as they are hastily searching for genetic solutions to the devastating HLB crisis through breeding, genetic engineering, or genome editing.

Construction of High-Density Genetic Maps and Detection of QTLs Associated With Huanglongbing Tolerance in Citrus

Huanglongbing (HLB), or citrus greening, is the most devastating disease in citrus worldwide. Commercial citrus varieties including sweet orange (Citrus sinensis) are highly susceptible to HLB, and trifoliate orange (Poncirus trifoliata, a close Citrus relative) is widely considered resistant or highly tolerant to HLB. In this study, an intergeneric F₁ population of sweet orange and trifoliate orange was genotyped by Genotyping-by-Sequencing, and high-density SNP-based genetic maps were constructed separately for trifoliate orange and sweet orange. The two genetic maps exhibited high synteny and high coverage of the citrus genome. Progenies of the F₁ population and their parents were planted in a replicated field trial, exposed to intense HLB pressure for 3 years, and then evaluated for susceptibility to HLB over 2 years. The F₁ population exhibited a wide range in severity of HLB foliar symptom and canopy damage. Genome-wide QTL analysis based on the phenotypic data of foliar symptom and canopy damage in 2 years identified three clusters of repeatable QTLs in trifoliate orange linkage groups LG-t6, LG-t8 and LG-t9. Co-localization of QTLs for two traits was observed within all three regions. Additionally, one cluster of QTLs in sweet orange (linkage group LG-s7) was also detected. The majority of the identified QTLs each explained 18-30% of the phenotypic variation, indicating their major role in determining HLB responses. These results show, for the first time, a quantitative genetic nature yet the presence of major loci for the HLB tolerance in trifoliate orange. The results suggest that sweet orange also contains useful genetic factor(s) for improving HLB tolerance in commercial citrus varieties. Findings from this study should be very valuable and timely to researchers worldwide as they are hastily searching for genetic solutions to the devastating HLB crisis through breeding, genetic engineering, or genome editing.

Unreduced Megagametophyte Production in Lemon Occurs via Three Meiotic Mechanisms, Predominantly Second-Division Restitution

Unreduced (2n) gametes have played a pivotal role in polyploid plant evolution and are useful for sexual polyploid breeding in various species, particularly for developing new seedless citrus varieties. The underlying mechanisms of 2n gamete formation were recently revealed for Citrus reticulata but remain poorly understood for other citrus species, including lemon (C. limon [L.] Burm. f.). Here, we investigated the frequency and causal meiotic mechanisms of 2n megagametophyte production in lemon. We genotyped 48progeny plants of two lemon genotypes, "Eureka Frost" and "Fino", using 16 Simple Sequence Repeat (SSR) and 18 Single Nucleotide Polymorphism (SNP) markers to determine the genetic origin of the progenies and the underlying mechanisms for 2n gamete formation. We utilized a maximum-likelihood method based on parental heterozygosity restitution (PHR) of centromeric markers and analysis of PHR patterns along the chromosome. The frequency of 2n gamete production was 4.9% for "Eureka Frost" and 8.3% for "Fino", with three meiotic mechanisms leading to 2n gamete formation. We performed the maximum-likelihood method at the individual level via centromeric marker analysis, finding that 88% of the hybrids arose from second-division restitution (SDR), 7% from first-division restitution (FDR) or pre-meiotic doubling (PRD), and 5% from post-meiotic genome doubling (PMD). The pattern of PHR along LG1 confirmed that SDR is the main mechanism for 2n gamete production. Recombination analysis between markers in this LG revealed partial chiasma interference on both arms. We discuss the implications of these restitution mechanisms for citrus breeding and lemon genetics.

Fine Mapping for Identification of Citrus Alternaria Brown Spot Candidate Resistance Genes and Development of New SNP Markers for Marker-Assisted Selection

Alternaria brown spot (ABS) is a serious disease affecting susceptible citrus genotypes, which is a strong concern regarding citrus breeding programs. Resistance is conferred by a recessive locus (ABSr) previously located by our group within a 3.3 Mb genome region near the centromere in chromosome III. This work addresses fine-linkage mapping of this region for identifying candidate resistance genes and develops new molecular markers for ABS-resistance effective marker-assisted selection (MAS). Markers closely linked to ABSr locus were used for fine mapping using a 268-segregating diploid progeny derived from a heterozygous susceptible × resistant cross. Fine mapping limited the genomic region containing the ABSr resistance gene to 366 kb, flanked by markers at 0.4 and 0.7 cM. This region contains nine genes related to pathogen resistance. Among them, eight are resistance (R) gene homologs, with two of them harboring a serine/threonine protein kinase domain. These two genes along with a gene encoding a S-adenosyl-L-methionine-dependent-methyltransferase protein, should be considered as strong candidates for ABS-resistance. Moreover, the closest SNP was genotyped in 40 citrus varieties, revealing very high association with the resistant/susceptible phenotype. This new marker is currently used in our citrus breeding program for ABS-resistant parent and cultivar selection, at diploid, triploid and tetraploid level.

Precocious flowering of juvenile citrus induced by a viral vector based on Citrus leaf blotch virus: a new tool for genetics and breeding

The long juvenile period of citrus trees (often more than 6 years) has hindered genetic improvement by traditional breeding methods and genetic studies. In this work, we have developed a biotechnology tool to promote transition from the vegetative to the reproductive phase in juvenile citrus plants by expression of the Arabidopsis thaliana or citrus FLOWERING LOCUS T (FT) genes using a Citrus leaf blotch virus-based vector (clbvINpr-AtFT and clbvINpr-CiFT, respectively). Citrus plants of different genotypes graft inoculated with either of these vectors started flowering within 4-6 months, with no alteration of the plant architecture, leaf, flower or fruit morphology in comparison with noninoculated adult plants. The vector did not integrate in or recombine with the plant genome nor was it pollen or vector transmissible, albeit seed transmission at low rate was detected. The clbvINpr-AtFT is very stable, and flowering was observed over a period of at least 5 years. Precocious flowering of juvenile citrus plants after vector infection provides a helpful and safe tool to dramatically speed up genetic studies and breeding programmes.

LRR-RLK family from two Citrus species: genome-wide identification and evolutionary aspects

Background

Leucine-rich repeat receptor-like kinases (LRR-RLKs) represent the largest subfamily of plant RLKs. The functions of most LRR-RLKs have remained undiscovered, and a few that have been experimentally characterized have been shown to have important roles in growth and development as well as in defense responses. Although RLK subfamilies have been previously studied in many plants, no comprehensive study has been performed on this gene family in Citrus species, which have high economic importance and are frequent targets for emerging pathogens. In this study, we performed in silico analysis to identify and classify LRR-RLK homologues in the predicted proteomes of Citrus clementina (clementine) and Citrus sinensis (sweet orange). In addition, we used large-scale phylogenetic approaches to elucidate the evolutionary relationships of the LRR-RLKs and further narrowed the analysis to the LRR-XII group, which contains several previously described cell surface immune receptors.

Results

We built integrative protein signature databases for Citrus clementina and Citrus sinensis using all predicted protein sequences obtained from whole genomes. A total of 300 and 297 proteins were identified as LRR-RLKs in C. clementina and C. sinensis, respectively. Maximum-likelihood phylogenetic trees were estimated using Arabidopsis LRR-RLK as a template and they allowed us to classify Citrus LRR-RLKs into 16 groups. The LRR-XII group showed a remarkable expansion, containing approximately 150 paralogs encoded in each Citrus genome. Phylogenetic analysis also demonstrated the existence of two distinct LRR-XII clades, each one constituted mainly by RD and non-RD kinases. We identified 68 orthologous pairs from the C. clementina and C. sinensis LRR-XII genes. In addition, among the paralogs, we identified a subset of 78 and 62 clustered genes probably derived from tandem duplication events in the genomes of C. clementina and C. sinensis, respectively.

Conclusions

This work provided the first comprehensive evolutionary analysis of the LRR-RLKs in Citrus. A large expansion of LRR-XII in Citrus genomes suggests that it might play a key role in adaptive responses in host-pathogen co-evolution, related to the perennial life cycle and domestication of the citrus crop species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2930-9) contains supplementary material, which is available to authorized users.

Maximum-likelihood method identifies meiotic restitution mechanism from heterozygosity transmission of centromeric loci: application in citrus

Polyploidisation is a key source of diversification and speciation in plants. Most researchers consider sexual polyploidisation leading to unreduced gamete as its main origin. Unreduced gametes are useful in several crop breeding schemes. Their formation mechanism, i.e., First-Division Restitution (FDR) or Second-Division Restitution (SDR), greatly impacts the gametic and population structures and, therefore, the breeding efficiency. Previous methods to identify the underlying mechanism required the analysis of a large set of markers over large progeny. This work develops a new maximum-likelihood method to identify the unreduced gamete formation mechanism both at the population and individual levels using independent centromeric markers. Knowledge of marker-centromere distances greatly improves the statistical power of the comparison between the SDR and FDR hypotheses. Simulating data demonstrated the importance of selecting markers very close to the centromere to obtain significant conclusions at individual level. This new method was used to identify the meiotic restitution mechanism in nineteen mandarin genotypes used as female parents in triploid citrus breeding. SDR was identified for 85.3% of 543 triploid hybrids and FDR for 0.6%. No significant conclusions were obtained for 14.1% of the hybrids. At population level SDR was the predominant mechanisms for the 19 parental mandarins.

Genetic mapping of centromeres in the nine Citrus clementina chromosomes using half-tetrad analysis and recombination patterns in unreduced and haploid gametes

BACKGROUND

Mapping centromere locations in plant species provides essential information for the analysis of genetic structures and population dynamics. The centromere's position affects the distribution of crossovers along a chromosome and the parental heterozygosity restitution by 2n gametes is a direct function of the genetic distance to the centromere. Sexual polyploidisation is relatively frequent in Citrus species and is widely used to develop new seedless triploid cultivars. The study's objectives were to (i) map the positions of the centromeres of the nine Citrus clementina chromosomes; (ii) analyse the crossover interference in unreduced gametes; and (iii) establish the pattern of genetic recombination in haploid clementine gametes along each chromosome and its relationship with the centromere location and distribution of genic sequences.

RESULTS

Triploid progenies were derived from unreduced megagametophytes produced by second-division restitution. Centromere positions were mapped genetically for all linkage groups using half-tetrad analysis. Inference of the physical locations of centromeres revealed one acrocentric, four metacentric and four submetacentric chromosomes. Crossover interference was observed in unreduced gametes, with variation seen between chromosome arms. For haploid gametes, a strong decrease in the recombination rate occurred in centromeric and pericentromeric regions, which contained a low density of genic sequences. In chromosomes VIII and IX, these low recombination rates extended beyond the pericentromeric regions. The genomic region corresponding to a genetic distance < 5cM from a centromere represented 47% of the genome and 23% of the genic sequences.

CONCLUSIONS

The centromere positions of the nine citrus chromosomes were genetically mapped. Their physical locations, inferred from the genetic ones, were consistent with the sequence constitution and recombination pattern along each chromosome. However, regions with low recombination rates extended beyond the pericentromeric regions of some chromosomes into areas richer in genic sequences. The persistence of strong linkage disequilibrium between large numbers of genes promotes the stability of epistatic interactions and multilocus-controlled traits over successive generations but also maintains multi-trait associations. Identification of the centromere positions will allow the development of simple methods to analyse unreduced gamete formation mechanisms in a large range of genotypes and further modelling of genetic inheritance in sexual polyploidisation breeding schemes.

QoI-Resistance Stability in Relation to Pathogenic and Saprophytic Fitness Components of Alternaria alternata from Citrus

The phenotypic stability, fitness components, and ability to cause disease of isolates of the Alternaria alternata tangerine pathotype resistant to quinone-outside inhibitors (QoIs) were studied. Stability of resistance to azoxystrobin (AZ) and pyraclostrobin (PYR) was determined after consecutive transfers on potato dextrose agar (PDA). The sensitivity to QoIs did not change significantly after 10 transfers on PDA compared with the initial sensitivity of all isolates tested. Fitness components evaluated in vitro were mycelial growth, conidial germination, and conidial production. Incubation period, number of lesions per leaf area, and virulence were determined with detached leaf assays using four cultivars: Dancy, Minneola, Murcott, and Sunburst. Variability in fitness components was observed among isolates within the same sensitivity group. As a group, no significant differences in the mean values of these fitness components were observed between resistant and sensitive phenotypes, except for virulence. Resistant isolates were significantly (P < 0.05) more virulent than the sensitive isolates on Dancy, Minneola, and Sunburst but not on Murcott (P = 0.3506). There was no significant correlation between individual fitness components and the level of sensitivity to AZ and PYR. Preventive applications of Abound (commercial formulation of AZ) at full field rates failed to control disease caused by QoI-resistant isolates under greenhouse conditions. Our results suggest that QoI resistance in A. alternata tangerine pathotype is stable in the absence of QoI selection pressure and that resistance development did not affect the fitness of resistant isolates.

Genomic variants of genes associated with three horticultural traits in apple revealed by genome re-sequencing

The apple (Malus × domestica Borkh.) cultivar 'Su Shuai' exhibits greater disease resistance, shorter internodes and lighter fruit flavor compared with its parents 'Golden Delicious' and 'Indo'. To obtain a comprehensive overview of the sequence variation in these three horticultural traits, the genomes of 'Su Shuai' and 'Indo' were resequenced using next-generation sequencing and compared to the genome of 'Golden Delicious'. A wide range of genetic variations were detected, including 2 454 406 and 18 749 349 single nucleotide polymorphism (SNP) and 59 547 and 50 143 structural variants (SVs) in the 'Indo' and 'Su Shuai' genomes, respectively. Among the SVs in 'Su Shuai', 17 genes related to disease resistance, 10 genes related to Gibberellin (GA) and 19 genes associated with fruit flavor were identified. The expression patterns of eight of the SV genes were examined using reverse transcription-quantitative polymerase chain reaction (RT-qPCR). The results of this study illustrate the genomic variation in these cultivars and provide evidence for a genetic basis for the horticultural traits of disease resistance, short internodes and lighter flavor exhibited in these cultivars. These results provide a genetic basis for the phenotypic characteristics of 'Su Shuai' and, as such, these SVs could serve as gene-specific molecular markers in maker-assisted breeding of apples.