EST-SSR genetic maps for Citrus sinensis and Poncirus trifoliata

Whole-genome sequencing and variant discovery of Citrus reticulata "Kinnow" from Pakistan

Citrus is a source of nutritional and medicinal advantages, cultivated worldwide with major groups of sweet oranges, mandarins, grapefruits, kumquats, lemons and limes. Pakistan produces all major citrus groups with mandarin (Citrus reticulata) being the prominent group that includes local commercial cultivars Feutral's Early, Dancy, Honey, and Kinnow. The present study designed to understand the genetic architecture of this unique variety of Citrus reticulata 'Kinnow.' The whole-genome resequencing and variant calling was performed to map the genomic variability that might be responsible for its particular characteristics like taste, seedlessness, juice content, thickness of peel, and shelf-life. A total of 139,436,350 raw sequence reads were generated with 20.9 Gb data in Fastq format having 98% effectiveness and 0.2% base call error rate. Overall, 3,503,033 SNPs, 176,949 MNPs, 323,287 INS, and 333,083 DEL were identified using the GATK4 variant calling pipeline against Citrus clementina. Furthermore, g:Profiler was applied for annotating the newly found variants, harbor genes/transcripts and their involved pathways. A total of 73,864 transcripts harbors 4,336,352 variants, most of the observed variants were predicted in non-coding regions and 1009 transcripts were found well annotated by different databases. Out of total aforementioned transcripts, 588 involved in biological processes, 234 in molecular functions and 167 transcripts in cellular components. In a nutshell, 18,153 high impact variants and 216 genic variants found in the current study, which may be used after its functional validation for marker-assisted breeding programs of "Kinnow" to propagate its valued traits for the improvement of contemporary citrus varieties in the region.

Meiotic Behaviors of Allotetraploid Citrus Drive the Interspecific Recombination Landscape, the Genetic Structures, and Traits Inheritance in Tetrazyg Progenies Aiming to Select New Rootstocks

Sexual breeding at the tetraploid level is a promising strategy for rootstock breeding in citrus. Due to the interspecific origin of most of the conventional diploid citrus rootstocks that produced the tetraploid germplasm, the optimization of this strategy requires better knowledge of the meiotic behavior of the tetraploid parents. This work used Genotyping By Sequencing (GBS) data from 103 tetraploid hybrids to study the meiotic behavior and generate a high-density recombination landscape for their tetraploid intergenic Swingle citrumelo and interspecific Volkamer lemon progenitors. A genetic association study was performed with root architecture traits. For citrumelo, high preferential chromosome pairing was revealed and led to an intermediate inheritance with a disomic tendency. Meiosis in Volkamer lemon was more complex than that of citrumelo, with mixed segregation patterns from disomy to tetrasomy. The preferential pairing resulted in low interspecific recombination levels and high interspecific heterozygosity transmission by the diploid gametes. This meiotic behavior affected the efficiency of Quantitative Trait Loci (QTL) detection. Nevertheless, it enabled a high transmission of disease and pest resistance candidate genes from P. trifoliata that are heterozygous in the citrumelo progenitor. The tetrazyg strategy, using doubled diploids of interspecific origin as parents, appears to be efficient in transferring the dominant traits selected at the parental level to the tetraploid progenies.

Construction of a high density genetic linkage map to define the locus conferring seedlessness from Mukaku Kishu mandarin

Mukaku Kishu ('MK'), a small sized mandarin, is an important source of seedlessness in citrus breeding. Identification and mapping the gene(s) governing 'MK' seedlessness will expedite seedless cultivar development. In this study, two 'MK'-derived mapping populations- LB8-9 Sugar Belle^® ('SB') × 'MK' (N=97) and Daisy ('D') × 'MK' (N=68) were genotyped using an Axiom_Citrus56 Array encompassing 58,433 SNP probe sets, and population specific male and female parent linkage maps were constructed. The parental maps of each population were integrated to produce sub-composite maps, which were further merged to develop a consensus linkage map. All the parental maps (except 'MK_D') had nine major linkage groups, and contained 930 ('SB'), 810 ('MK_SB'), 776 ('D') and 707 ('MK_D') SNPs. The linkage maps displayed 96.9 ('MK_D') to 98.5% ('SB') chromosomal synteny with the reference Clementine genome. The consensus map was comprised of 2588 markers including a phenotypic seedless (Fs)-locus and spanned a genetic distance of 1406.84 cM, with an average marker distance of 0.54 cM, which is substantially lower than the reference Clementine map. For the phenotypic Fs-locus, the distribution of seedy and seedless progenies in both 'SB' × 'MK' (55:42, χ² = 1.74) and 'D' × 'MK' populations (33:35, χ² = 0.06) followed a test cross pattern. The Fs-locus mapped on chromosome 5 with SNP marker 'AX-160417325' at 7.4 cM in 'MK_SB' map and between two SNP markers 'AX-160536283' and 'AX-160906995' at a distance of 2.4 and 4.9 cM, respectively in 'MK_D' map. The SNPs 'AX-160417325' and 'AX-160536283' correctly predicted seedlessness of 25-91.9% progenies in this study. Based on the alignment of flanking SNP markers to the Clementine reference genome, the candidate gene for seedlessness hovered in a ~ 6.0 Mb region between 3.97 Mb (AX-160906995) to 10.00 Mb (AX-160536283). This region has 131 genes of which 13 genes (belonging to seven gene families) reportedly express in seed coat or developing embryo. The findings of the study will prove helpful in directing future research for fine mapping this region and eventually underpinning the exact causative gene governing seedlessness in 'MK'.

Potential application of simple easy-to-use insertion-deletion (InDel) markers in citrus cultivar identification

We previously developed insertion-deletion (InDel) markers that distinguish three genotypes (two homozygous and one heterozygous) of diverse citrus cultivars. These InDel markers were codominant and could be clearly detected by using simple agarose gel electrophoresis. We sought to establish a method for cultivar identification using these 28 InDel markers to genotype 31 citrus cultivars. The results revealed that a minimum of 6 markers were required to identify individuals using the three-genotype classification method. Furthermore, we found that a simple method for distinguishing between two genotypes (homozygous and heterozygous) could be used to identify individuals using a minimum of 7 markers. Our findings provide a basis for the development of simple and rapid citrus cultivar identification methods.

Developing functional markers for vitamin E biosynthesis in oil palm

Vitamin E is essential for human health and plays positive roles in anti-oxidation. Previously, we detected large variation in vitamin E content among 161 oil palm accessions. In this study, twenty oil palm accessions with distinct variation in vitamin E contents (171.30 to 1 258.50 ppm) were selected for genetic variation analysis and developing functional markers associated with vitamin E contents. Thirty-seven homologous genes in oil palm belonging to vitamin E biosynthesis pathway were identified via BLASTP analysis, the lengths of which ranged from 426 to 25 717 bp (average 7 089 bp). Multiplex PCR sequencing for the 37 genes found 1 703 SNPs and 85 indels among the 20 oil palm accessions, with 226 SNPs locating in the coding regions. Clustering analysis for these polymorphic loci showed that the 20 oil palm accessions could be divided into five groups. Among these groups, group I included eight oil palm accessions whose vitamin E content (mean value: 893.50 ppm) was far higher than other groups (mean value 256.29 to 532.94 ppm). Correlation analysis between the markers and vitamin E traits showed that 134 SNP and 7 indel markers were significantly (p < 0.05) related with total vitamin E content. Among these functional markers, the indel EgTMT-1-24 was highly correlated with variation in vitamin E content, especially tocotrienol content. Our study identified a number of candidate function associated markers and provided clues for further research into molecular breeding for high vitamin E content oil palm.

Utilization of somatic fusion techniques for the development of HLB tolerant breeding resources employing the Australian finger lime (Citrus australasica)

The Australian finger lime is a unique citrus species that has gained importance due to its unique fruit characteristics and perceived tolerance to Huanglongbing (HLB), an often-fatal disease of citrus trees. In this study, we developed allotetraploid finger lime hybrids and cybrids by utilizing somatic cell fusion techniques to fuse diploid ‘OLL8’ sweet orange or ‘Page’ tangelo callus-derived protoplasts with finger lime (FL) mesophyll-derived protoplasts. Six somatic fusions were regenerated from the ‘OLL8’ + FL fusion, while three putative cybrids were regenerated from the ‘Page’ + FL fusion. Ploidy levels and nuclear-expressed sequence tag derived simple sequence repeat (EST-SSR) markers confirmed the somatic hybrid production, and mitochondrial DNA primer sets confirmed the cybrid nature. Several trees produced by the somatic fusion remained HLB negative even after 6 years of growth in an HLB-endemic environment. Pathogenesis related (PR) and other genes that are often upregulated in HLB-tolerant trees were also upregulated in our somatic fusions. These newly developed somatic fusions and cybrids could potentially be used as breeding parents to develop the next generation of improved HLB-tolerant rootstocks and scions.

High-Density Genetic Map Construction and Identification of QTLs Controlling Leaf Abscission Trait in Poncirus trifoliata

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.

Expression Quantitative Trait Loci (eQTL) mapping for callose synthases in intergeneric hybrids of Citrus challenged with the bacteria Candidatus Liberibacter asiaticus

Citrus plants have been extremely affected by Huanglongbing (HLB) worldwide, causing economic losses. HLB disease causes disorders in citrus plants, leading to callose deposition in the phloem vessel sieve plates. Callose is synthesized by callose synthases, which are encoded by 12 genes (calS1– calS12)in Arabidopsis thaliana. We evaluated the expression of eight callose synthase genes from Citrus in hybrids between Citrus sunki and Poncirus trifoliata infected with HLB. The objective of this work was to identify possible tolerance loci combining the expression quantitative trait loci (eQTL) of different callose synthases and genetic Single-Nucleotide Polymorphism (SNP) maps of C. sunki and P. trifoliata. The expression data from all CscalS ranged widely among the hybrids. Furthermore, the data allowed the detection of 18 eQTL in the C. sunki map and 34 eQTL in the P. trifoliata map. In both maps, some eQTL for different CscalS were overlapped; thus, a single region could be associated with the regulation of more than one CscalS. The regions identified in this work can be interesting targets for future studies of Citrus breeding programs to manipulate callose synthesis during HLB infection.

Analysis of Crossover Events and Allele Segregation Distortion in Interspecific Citrus Hybrids by Single Pollen Genotyping

In citrus, a classical method of studying crossovers and segregation distortion (SD) is the genetic analysis of progenies. A new strategy combining fluorescence-activated cell sorting and whole genome amplification of haploid pollen nuclei with a large set of molecular markers, offers the opportunity to efficiently determine the frequency of crossovers and the identification of SD without the need to generate segregating populations. Here we have analyzed meiotic crossover events in a pollen nuclei population from "Eureka" lemon and the allelic SD was evaluated in a pollen nuclei population from a clementine × sweet orange hybrid ("CSO"). Data obtained from the "CSO" pollen nuclei population were compared to those obtained from genotyping of a segregating population ("RTSO") arising from a hand-made sexual hybridization between diploid non apomictic selected tangor (mandarin × sweet orange; "RTO" tangor) as female parent pollinated with "CSO" tangor as male parent. The analysis of crossovers rates on chromosome 1 revealed the presence of up to five crossovers events on one arm and four on the corresponding other arm, with an average of 1.97 crossovers per chromosome while no crossover events were observed in five "Eureka" lemon pollen nuclei. The rate of SD observed in "CSO" pollen nuclei (13.8%) was slightly lower than that recovered in the "RTSO" population (20.7%). In the pollen nuclei population, SD was found on linkage group (LG) 2, while the "RTSO" population showed SD on LGs 2 and 7. Potential male gametic selection mechanisms were distinguished in pollen grains, while in the population, mechanisms of gametophytic selection and/or zygotic selection were observed. This methodology is a very useful tool to facilitate research focused on the reproductive biology of citrus and study the mechanisms that affect crossovers and SD.

SLAF-Based Construction of a High-Density Genetic Map and Its Application in QTL Mapping of Carotenoids Content in Citrus Fruit

Carotenoids are important antioxidant components in the human diet. To develop carotenoid-rich agricultural products by genetic intervention, understanding the genetic basis of carotenoids variation is essential. In this study, we constructed a high-density integrated genetic map with 3817 molecular markers using specific locus amplified fragment (SLAF) sequencing from a C. reticulata × P. trifoliata F₁ pseudotestcross population. A total of 17 significant quantitative trait loci (QTLs) distributed on Chromosomes (Chr) 2, 3, 5, 6, and 9 were detected to determine the carotenoid variation in the population. In particular, three QTL colocalizations for multiple carotenoid constituents were observed on Chr 2, 3, and 9, one of which was located on Chr2:34,654,608-35430715 accounted for 20.1-25.4% of the variation of luteoxanthin, auroxanthin, lutein, violaxanthin, and total carotenoid content. Overall, this study provides a genetic foundation for marker-assisted selection (MAS) breeding of nutritionally enhanced citrus fruit.

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.

Preferential Homologous Chromosome Pairing in a Tetraploid Intergeneric Somatic Hybrid (Citrus reticulata + Poncirus trifoliata) Revealed by Molecular Marker Inheritance

The creation of intergeneric somatic hybrids between Citrus and Poncirus is an efficient approach for citrus rootstock breeding, offering the possibility of combining beneficial traits from both genera into novel rootstock lineages. These somatic hybrids are also used as parents for further tetraploid sexual breeding. In order to optimize these latter breeding schemes, it is essential to develop knowledge on the mode of inheritance in the intergeneric tetraploid hybrids. We assessed the meiotic behavior of an intergeneric tetraploid somatic hybrid resulting from symmetric protoplast fusion of diploid Citrus reticulata and diploid Poncirus trifoliata. The analysis was based on the segregation patterns of 16 SSR markers and 9 newly developed centromeric/pericentromeric SNP markers, representing all nine linkage groups of the Citrus genetic map. We found strong but incomplete preferential pairing between homologues of the same ancestral genome. The proportion of gametes that can be explained by random meiotic chromosome associations (τ) varied significantly between chromosomes, from 0.09 ± 0.02 to 0.47 ± 0.09, respectively, in chromosome 2 and 1. This intermediate inheritance between strict disomy and tetrasomy, with global preferential disomic tendency, resulted in a high level of intergeneric heterozygosity of the diploid gametes. Although limited, intergeneric recombinations occurred, whose observed rates, ranging from 0.09 to 0.29, respectively, in chromosome 2 and 1, were significantly correlated with τ. Such inheritance is of particular interest for rootstock breeding because a large part of the multi-trait value selected at the teraploid parent level is transmitted to the progeny, while the potential for some intergeneric recombination offers opportunities for generating plants with novel allelic combinations that can be targeted by selection.

QTL mapping of male sterility and transmission pattern in progeny of Satsuma mandarin

Seedlessness is one of the important traits in citrus breeding. Male sterility derived from Satsuma mandarin (Citrus unshiu) has been used in Japanese citrus breeding programs to obtain seedless cultivars. The efficiency of seedless cultivar breeding would be improved by developing a selection marker linked to seedlessness. In this study, we performed QTL mapping in ‘Okitsu No. 46’ × ‘Okitsu No. 56’ (O46-O56) crosses for the number of pollen grains per anther (NPG) and apparent pollen fertility (APF), two traits used as an index of male sterility, and detected a candidate QTL for NPG (MS-P1) on linkage group 8 with a significant LOD score (7.31) and 47% of variance explained. The QTL for APF (MS-F1) was detected on linkage group 6 with a significant LOD score (5.71) and 63.6% of variance explained. The role of both MS-P1 in reducing NPG and MS-F1 in decreasing APF were confirmed with the ‘Okitsu No.46’ × ‘Kara’ (O46-K) cross. Pedigree analysis inferred that both MS-P1 and MS-F1 in ‘Okitsu No. 46’ were derived from kunenbo (Citrus nobilis) through hassaku (C. hassaku) and ‘Sweet Spring’. Cytoplasm analysis revealed that both male-sterile ‘Sweet Spring’ and ‘Okitsu No. 46’ have cytoplasm derived from Kishu (C. kinokuni hort. ex Tanaka), but the cytoplasm of male-sterile kunenbo and hassaku were derived from other varieties rather than Kishu. These results suggest that MS-P1 and MS-F1 primarily reduce the NPG and decrease APF, but their expression requires a cytoplasm derived from Kishu. These findings will improve our understanding of the molecular mechanism of male sterility in citrus and help to develop a DNA marker for seedless breeding in citrus.

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.

QTL mapping for fruit quality in Citrus using DArTseq markers

Background

Citrus breeding programs have many limitations associated with the species biology and physiology, requiring the incorporation of new biotechnological tools to provide new breeding possibilities. Diversity Arrays Technology (DArT) markers, combined with next-generation sequencing, have wide applicability in the construction of high-resolution genetic maps and in quantitative trait locus (QTL) mapping. This study aimed to construct an integrated genetic map using full-sib progeny derived from Murcott tangor and Pera sweet orange and DArTseq™ molecular markers and to perform QTL mapping of twelve fruit quality traits. A controlled Murcott x Pera crossing was conducted at the Citrus Germplasm Repository at the Sylvio Moreira Citrus Centre of the Agronomic Institute (IAC) located in Cordeirópolis, SP, in 1997. In 2012, 278 F₁ individuals out of a family of 312 confirmed hybrid individuals were analyzed for fruit traits and genotyped using the DArTseq markers. Using OneMap software to obtain the integrated genetic map, we considered only the DArT loci that showed no segregation deviation. The likelihood ratio and the genomic information from the available Citrus sinensis L. Osbeck genome were used to determine the linkage groups (LGs).

Results

The resulting integrated map contained 661 markers in 13 LGs, with a genomic coverage of 2,774 cM and a mean density of 0.23 markers/cM. The groups were assigned to the nine Citrus haploid chromosomes; however, some of the chromosomes were represented by two LGs due the lack of information for a single integration, as in cases where markers segregated in a 3:1 fashion. A total of 19 QTLs were identified through composite interval mapping (CIM) of the 12 analyzed fruit characteristics: fruit diameter (cm), height (cm), height/diameter ratio, weight (g), rind thickness (cm), segments per fruit, total soluble solids (TSS, %), total titratable acidity (TTA, %), juice content (%), number of seeds, TSS/TTA ratio and number of fruits per box. The genomic sequence (pseudochromosomes) of C. sinensis was compared to the genetic map, and synteny was clearly identified. Further analysis of the map regions with the highest LOD scores enabled the identification of putative genes that could be associated with the fruit quality characteristics.

Conclusion

An integrated linkage map of Murcott tangor and Pera sweet orange using DArTseq™ molecular markers was established and it was useful to perform QTL mapping of twelve fruit quality traits. The next generation sequences data allowed the comparison between the linkage map and the genomic sequence (pseudochromosomes) of C. sinensis and the identification of genes that may be responsible for phenotypic traits in Citrus. The obtained linkage map was used to assign sequences that had not been previously assigned to a position in the reference genome.

Electronic supplementary material

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

Hybrid Origins of Citrus Varieties Inferred from DNA Marker Analysis of Nuclear and Organelle Genomes

Most indigenous citrus varieties are assumed to be natural hybrids, but their parentage has so far been determined in only a few cases because of their wide genetic diversity and the low transferability of DNA markers. Here we infer the parentage of indigenous citrus varieties using simple sequence repeat and indel markers developed from various citrus genome sequence resources. Parentage tests with 122 known hybrids using the selected DNA markers certify their transferability among those hybrids. Identity tests confirm that most variant strains are selected mutants, but we find four types of kunenbo (Citrus nobilis) and three types of tachibana (Citrus tachibana) for which we suggest different origins. Structure analysis with DNA markers that are in Hardy-Weinberg equilibrium deduce three basic taxa coinciding with the current understanding of citrus ancestors. Genotyping analysis of 101 indigenous citrus varieties with 123 selected DNA markers infers the parentages of 22 indigenous citrus varieties including Satsuma, Temple, and iyo, and single parents of 45 indigenous citrus varieties, including kunenbo, C. ichangensis, and Ichang lemon by allele-sharing and parentage tests. Genotyping analysis of chloroplast and mitochondrial genomes using 11 DNA markers classifies their cytoplasmic genotypes into 18 categories and deduces the combination of seed and pollen parents. Likelihood ratio analysis verifies the inferred parentages with significant scores. The reconstructed genealogy identifies 12 types of varieties consisting of Kishu, kunenbo, yuzu, koji, sour orange, dancy, kobeni mikan, sweet orange, tachibana, Cleopatra, willowleaf mandarin, and pummelo, which have played pivotal roles in the occurrence of these indigenous varieties. The inferred parentage of the indigenous varieties confirms their hybrid origins, as found by recent studies.

Patterns of genetic structure and evidence of gene flow among Tunisian Citrus species based on informative nSSR markers

This study investigates the extent of genetic diversity, phylogenetic relationships and the amount of gene flow among Tunisian Citrus species based on a set of 15 informative nuclear SSR molecular markers. Genotyping data highlighted an allelic richness among Tunisian Citrus species and has allowed the detection of 168 alleles among them 104.19 were effective. The partition of the total genetic diversity (HT=0.832) showed that the highest amount of variation within the Citrus species is HS=0.550, while the relative amount of the between-species genetic diversity GST does not exceed 0.338. This pattern of genetic structure was supported by low-to-moderate FST pairwise values and the presence of a gene flow (Nm) among the eight Citrus species. The lowest genetic differentiation was revealed between the species C. sinensis and C. insitorum (FST=0.111, Nm=1.99), while the highest genetic differentiation was recorded between the species C. aurantifolia and C. paradisi (FST=0.367, Nm=0.43). The established Neighbor Joining analysis showed that all genotypes were widely discriminated and clearly pooled according to their species of origin, with minor exceptions.

Single-pollen genotyping to estimate mode of unreduced pollen formation in Citrus tamurana cv. Nishiuchi Konatsu

2n pollen formed by FDR in citrus. The Japanese local citrus cultivar, Nishiuchi Konatsu (Citrus tamurana hort. ex Tanaka; NK hereafter), has the ability to produce unreduced 2n pollen grains, allowing generation of polyploid progenies via sexual polyploidization. In this study, we developed a method of single-pollen genotyping for citrus and applied it to the analysis of transmission of heterozygosity in NK 2n pollen grains. Heterozygosity transmission was expressed as the percentage inheritance of a set of heterozygous alleles from the parent to the 2n gamete. The pathway of 2n pollen development was investigated by applying the observed heterozygosity transmission and genetic distance to two different map functions, for first division restitution (FDR) and second division restitution (SDR). The fit of the values observed for both functions was calculated, while virtually moving the centromere position. We screened for six heterozygous SSR (codominant microsatellite marker loci) in NK, all of which were expected to lie within the same linkage group. Pollen germination prior to DNA extraction was essential for this work, and 6-h incubation proved to be optimal for subsequent PCR amplification. Single-pollen genotyping unreduced NK 2n pollen grains revealed that heterozygosity transmission exceeded 50 % in all six alleles, and fitness tests indicated that the FDR map function better fitted the heterozygosity transmission observed rather than the SDR function. Our data thus strongly indicate that 2n pollen in NK is a result of first division restitution.

Comparative in silico analysis of EST-SSRs in angiosperm and gymnosperm tree genera

BACKGROUND

Simple Sequence Repeats (SSRs) derived from Expressed Sequence Tags (ESTs) belong to the expressed fraction of the genome and are important for gene regulation, recombination, DNA replication, cell cycle and mismatch repair. Here, we present a comparative analysis of the SSR motif distribution in the 5'UTR, ORF and 3'UTR fractions of ESTs across selected genera of woody trees representing gymnosperms (17 species from seven genera) and angiosperms (40 species from eight genera).

RESULTS

Our analysis supports a modest contribution of EST-SSR length to genome size in gymnosperms, while EST-SSR density was not associated with genome size in neither angiosperms nor gymnosperms. Multiple factors seem to have contributed to the lower abundance of EST-SSRs in gymnosperms that has resulted in a non-linear relationship with genome size diversity. The AG/CT motif was found to be the most abundant in SSRs of both angiosperms and gymnosperms, with a relative increase in AT/AT in the latter. Our data also reveals a higher abundance of hexamers across the gymnosperm genera.

CONCLUSIONS

Our analysis provides the foundation for future comparative studies at the species level to unravel the evolutionary processes that control the SSR genesis and divergence between angiosperm and gymnosperm tree species.

Sequence analysis reveals genomic factors affecting EST-SSR primer performance and polymorphism

This study was to explore genomic factors affecting the performance and polymorphism of 340 randomly selected EST-SSR (expressed sequence tag-simple sequence repeat) primers through BLAST of primer sequences to a reference genome. Genotyping showed 111 failed and 229 succeeded. The failed types included "no peaks" (NP, 69 primers), "weak peaks" (WP, 30), and "multiple peaks" (MP, 12). The successful types were divided into HM (homozygous between two selected parents, 78 primers) and HT (heterozygous at least in one parent, 151 primers). The BLAST revealed primer alignment status, genomic amplicon size (GAS), and genomic and expressed amplicon size difference (ASD). The alignment status was categorized as: "no hits found" (NHF); "multiple partial alignments" (MPA); "single partial alignment" (SPA); "multiple full alignments" (MFA); and "single full alignment" (SFA). NHF and partial alignment (PA) mainly resulted from discrepant nucleotides in contig-derived primers. The ASD separated 247 non-NHF primers into: "deletion", "same size", "insertion", "intron (GAS ≤500)", "intron (GAS >500)", and "error" categories. Most SFA primers were successful. About 88 % "error", 53 % NHF primers, and 47 % "intron (GAS >500)" failed. The "deletion" and "insertion" primers had the higher HT rates, and the "same size" had the highest HM rate. Optimized primer selection criteria are discussed.

Candidate gene markers for Candidatus Liberibacter asiaticus for detecting citrus greening disease

Citrus Huanglongbing (HLB) also known as citrus greening is one of the most devastating diseases of citrus worldwide. The disease is caused by Candidatus Liberibacter bacterium, vectored by the psyllid Diaphorina citri Kuwayama and Trioza erytreae Del Guercio. Citrus plants infected by the HLB bacterium may not show visible symptoms sometimes for years following infection. The aim of this study was to develop effective gene-specific primer pairs for polymerase chain reaction based method for quick screening of HLB disease. Thirty-two different gene-specific primer pairs, across the Ca. Liberibacter asiaticus genome, were successfully developed. The possibility of these primer pairs for cross-genome amplification across 'Ca. Liberibacter africanus' and 'Ca. Liberibacter americanus' were tested. The applicability of these primer pairs for detection and differentiation of Ca Liberibacter spp. is discussed.

Development and characterization of genomic and expressed SSRs in citrus by genome-wide analysis

Microsatellites or simple sequence repeats (SSRs) are one of the most popular sources of genetic markers and play a significant role in plant genetics and breeding. In this study, we identified citrus SSRs in the genome of Clementine mandarin and analyzed their frequency and distribution in different genomic regions. A total of 80,708 SSRs were detected in the genome with an overall density of 268 SSRs/Mb. While di-nucleotide repeats were the most frequent microsatellites in genomic DNA sequence, tetra-nucleotides, which had more repeat units than any other SSR types, had the highest cumulative sequence length. We identified 6,834 transcripts as containing 8,989 SSRs in 33,929 Clementine mandarin transcripts, among which, tri-nucleotide motifs (36.0%) were the most common, followed by di-nucleotide (26.9%) and hexa-nucleotide motifs (15.1%). The motif AG (16.7%) was most abundant among these SSRs, while motifs AAG (6.6%), AAT (5.0%), and TAG (2.2%) were most common among tri-nucleotides. Functional categorization of transcripts containing SSRs revealed that 5,879 (86.0%) of such transcripts had homology with known proteins, GO and KEGG annotation revealed that transcripts containing SSRs were those implicated in diverse biological processes in plants, including binding, development, transcription, and protein degradation. When 27 genomic and 78 randomly selected SSRs were tested on Clementine mandarin, 95 SSRs revealed polymorphism. These 95 SSRs were further deployed on 18 genotypes of the three generas of Rutaceae for the genetic diversity assessment, genomic SSRs generally show low transferability in comparison to SSRs developed from expressed sequences. These transcript-markers identified in our study may provide a valuable genetic and genomic tool for further genetic research and varietal development in citrus, such as diversity study, QTL mapping, molecular breeding, comparative mapping and other genetic analyses.

Cytological and molecular characterization of three gametoclones of Citrus clementina

BACKGROUND

Three gametoclonal plants of Citrus clementina Hort. ex Tan., cv. Nules, designated ESP, FRA, and ITA (derived from three labs in Spain, France, and Italy, respectively), were selected for cytological and molecular characterization in order to elucidate genomic rearrangements provoked by haploidization. The study included comparisons of their ploidy, homozygosity, genome integrity, and gene dosage, using chromosome counting, flow cytometry, SSR marker genotyping, and array-Comparative Genomic Hybridization (array-CGH).

RESULTS

Chromosome counting and flow cytometry revealed that ESP and FRA were haploid, but ITA was tri-haploid. Homozygous patterns, represented by a single peak (allele), were observed among the three plants at almost all SSR loci distributed across the entire diploid donor genome. Those few loci with extra peaks visualized as output from automated sequencing runs, generally low or ambiguous, might result from amplicons of paralogous members at the locus, non-specific sites, or unexpected recombinant alleles. No new alleles were found, suggesting the genomes remained stable and intact during gametogenesis and regeneration. The integrity of the haploid genome also was supported by array-CGH studies, in which genomic profiles were comparable to the diploid control.

CONCLUSIONS

The presence of few gene hybridization abnormalities, corroborated by gene dosage measurements, were hypothetically due to the segregation of hemizygous alleles and minor genomic rearrangements occurring during the haploidization procedure. In conclusion, these plants that are valuable genetic and breeding materials contain completely homozygous and essentially intact genomes.

EST-SSR Marker Sets for Practical Authentication of All Nine Registered Ginseng Cultivars in Korea

Panax ginseng has been cultivated for centuries, and nine commercial cultivars have been registered in Korea. However, these nine elite cultivars are grown in less than 10% of ginseng fields, and there is no clear authentication system for each cultivar even though their values are higher than those of local landraces. Here, we have developed 19 microsatellite markers using expressed gene sequences and established an authentication system for all nine cultivars. Five cultivars, ‘Chunpoong’, ‘Sunpoong’, ‘Gumpoong’, ‘Sunun’, and ‘Sunone’, can each be identified by one cultivar-unique allele, gm47n-a, gm47n-c, gm104-a, gm184-a (or gm129-a), and gm175-c, respectively. ‘Yunpoong’ can be identified by the co-appearance of gm47n-b and gm129-c. ‘Sunhyang’ can be distinguished from the other eight cultivars by the co-appearance of gm47n-b, gm129-b, and gm175-a. The two other cultivars, ‘Gopoong’ and ‘Cheongsun’, can be identified by their specific combinations of five marker alleles. This marker set was successfully utilized to identify the cultivars among 70 ginseng individuals and to select true F1 hybrid plants between two cultivars. We further analyzed the homogeneity of each cultivar and phylogenetic relationships among cultivars using these markers. This marker system will be useful to the seed industry and for breeding of ginseng.

Identification and characterization of gene-based SSR markers in date palm (Phoenix dactylifera L.)

BACKGROUND

Date palm (Phoenix dactylifera L.) is an important tree in the Middle East and North Africa due to the nutritional value of its fruit. Molecular Breeding would accelerate genetic improvement of fruit tree through marker assisted selection. However, the lack of molecular markers in date palm restricts the application of molecular breeding.

RESULTS

In this study, we analyzed 28,889 EST sequences from the date palm genome database to identify simple-sequence repeats (SSRs) and to develop gene-based markers, i.e. expressed sequence tag-SSRs (EST-SSRs). We identified 4,609 ESTs as containing SSRs, among which, trinucleotide motifs (69.7%) were the most common, followed by tetranucleotide (10.4%) and dinucleotide motifs (9.6%). The motif AG (85.7%) was most abundant in dinucleotides, while motifs AGG (26.8%), AAG (19.3%), and AGC (16.1%) were most common among trinucleotides. A total of 4,967 primer pairs were designed for EST-SSR markers from the computational data. In a follow up laboratory study, we tested a sample of 20 random selected primer pairs for amplification and polymorphism detection using genomic DNA from date palm cultivars. Nearly one-third of these primer pairs detected DNA polymorphism to differentiate the twelve date palm cultivars used. Functional categorization of EST sequences containing SSRs revealed that 3,108 (67.4%) of such ESTs had homology with known proteins.

CONCLUSION

Date palm EST sequences exhibits a good resource for developing gene-based markers. These genic markers identified in our study may provide a valuable genetic and genomic tool for further genetic research and varietal development in date palm, such as diversity study, QTL mapping, and molecular breeding.

A mixed model QTL analysis for sugarcane multiple-harvest-location trial data

Sugarcane-breeding programs take at least 12 years to develop new commercial cultivars. Molecular markers offer a possibility to study the genetic architecture of quantitative traits in sugarcane, and they may be used in marker-assisted selection to speed up artificial selection. Although the performance of sugarcane progenies in breeding programs are commonly evaluated across a range of locations and harvest years, many of the QTL detection methods ignore two- and three-way interactions between QTL, harvest, and location. In this work, a strategy for QTL detection in multi-harvest-location trial data, based on interval mapping and mixed models, is proposed and applied to map QTL effects on a segregating progeny from a biparental cross of pre-commercial Brazilian cultivars, evaluated at two locations and three consecutive harvest years for cane yield (tonnes per hectare), sugar yield (tonnes per hectare), fiber percent, and sucrose content. In the mixed model, we have included appropriate (co)variance structures for modeling heterogeneity and correlation of genetic effects and non-genetic residual effects. Forty-six QTLs were found: 13 QTLs for cane yield, 14 for sugar yield, 11 for fiber percent, and 8 for sucrose content. In addition, QTL by harvest, QTL by location, and QTL by harvest by location interaction effects were significant for all evaluated traits (30 QTLs showed some interaction, and 16 none). Our results contribute to a better understanding of the genetic architecture of complex traits related to biomass production and sucrose content in sugarcane.

Characterizing the citrus cultivar Carrizo genome through 454 shotgun sequencing

The citrus cultivar Carrizo is the single most important rootstock to the US citrus industry and has resistance or tolerance to a number of major citrus diseases, including citrus tristeza virus, foot rot, and Huanglongbing (HLB, citrus greening). A Carrizo genomic sequence database providing approximately 3.5× genome coverage (haploid genome size approximately 367 Mb) was populated through 454 GS FLX shotgun sequencing. Analysis of the repetitive DNA fraction indicated a total interspersed repeat fraction of 36.5%. Assembly and characterization of abundant citrus Ty3/gypsy elements revealed a novel type of element containing open reading frames encoding a viral RNA-silencing suppressor protein (RNA binding protein, rbp) and a plant cytokinin riboside 5′-monophosphate phosphoribohydrolase-related protein (LONELY GUY, log). Similar gypsy elements were identified in the Populus trichocarpa genome. Gene-coding region analysis indicated that 24.4% of the nonrepetitive reads contained genic regions. The depth of genome coverage was sufficient to allow accurate assembly of constituent genes, including a putative phloem-expressed gene. The development of the Carrizo database ( http://citrus.pw.usda.gov/ ) will contribute to characterization of agronomically significant loci and provide a publicly available genomic resource to the citrus research community.

Multilocus half-tetrad analysis and centromere mapping in citrus: evidence of SDR mechanism for 2n megagametophyte production and partial chiasma interference in mandarin cv 'Fortune'

The genetic structure of 2n gametes and, particularly, the parental heterozygosity restitution at each locus depends on the meiotic process by which they originated, with first-division restitution and second-division restitution (SDR) being the two major mechanisms. The origin of 2n gametes in citrus is still controversial, although sexual polyploidisation is widely used for triploid seedless cultivar development. In this study, we report the analysis of 2n gametes of mandarin cv 'Fortune' by genotyping 171 triploid hybrids with 35 simple sequence repeat markers. The microsatellite DNA allele counting-peak ratios method for allele-dosage evaluation proved highly efficient in segregating triploid progenies and allowed half-tetrad analysis (HTA) by inferring the 2n gamete allelic configuration. All 2n gametes arose from the female genitor. The observed maternal heterozygosity restitution varied between 10 and 82%, depending on the locus, thus SDR appears to be the mechanism underlying 2n gamete production in mandarin cv 'Fortune'. A new method to locate the centromere, based on the best fit between observed heterozygosity restitution within a linkage group and theoretical functions under either partial or no chiasmata interference hypotheses was successfully applied to linkage group II. The maximum value of heterozygosity restitution and the pattern of restitution along this linkage group would suggest there is partial chiasma interference. The implications of such a restitution mechanism for citrus breeding are discussed.

Global gene expression of Poncirus trifoliata, Citrus sunki and their hybrids under infection of Phytophthora parasitica

Background

Gummosis and root rot caused by Phytophthora are among the most economically important diseases in citrus. Four F₁ resistant hybrids (Pool R), and four F₁ susceptible hybrids (Pool S) to P. parasitica, were selected from a cross between susceptible Citrus sunki and resistant Poncirus trifoliata cv. Rubidoux. We investigated gene expression in pools of four resistant and four susceptible hybrids in comparison with their parents 48 hours after P. parasitica inoculation. We proposed that genes differentially expressed between resistant and susceptible parents and between their resistant and susceptible hybrids provide promising candidates for identifying transcripts involved in disease resistance. A microarray containing 62,876 UniGene transcripts selected from the CitEST database and prepared by NimbleGen Systems was used for analyzing global gene expression 48 hours after infection with P. parasitica.

Results

Three pairs of data comparisons (P. trifoliata/C. sunki, Pool R/C. sunki and Pool R/Pool S) were performed. With a filter of false-discovery rate less than 0.05 and fold change greater than 3.0, 21 UniGene transcripts common to the three pairwise comparative were found to be up-regulated, and 3 UniGene transcripts were down-regulated. Among them, our results indicated that the selected transcripts were probably involved in the whole process of plant defense responses to pathogen attack, including transcriptional regulation, signaling, activation of defense genes participating in HR, single dominant genes (R gene) such as TIR-NBS-LRR and RPS4 and switch of defense-related metabolism pathway. Differentially expressed genes were validated by RT-qPCR in susceptible and resistant plants and between inoculated and uninoculated control plants

Conclusions

Twenty four UniGene transcripts were identified as candidate genes for Citrus response to P. parasitica. UniGene transcripts were likely to be involved in disease resistance, such as genes potentially involved in secondary metabolite synthesis, intracellular osmotic adjustment, signal transduction pathways of cell death, oxidative burst and defense gene expression. Furthermore, our microarray data suggest another type of resistance in Citrus-Phytophthora interaction conferred by single dominant genes (R gene) since we encountered two previously reported R genes (TIR-NBS-LRR and RPS4) upregulated in the resistant genotypes relative to susceptible. We identified 7 transcripts with homology in other plants but yet unclear functional characterization which are an interesting pool for further analyses and 3 transcripts where no significant similarity was found. This is the first microarray study addressing an evaluation of transcriptional changes in response to P. parasitica in Citrus.

Lack of Evidence for Transmission of 'Candidatus Liberibacter asiaticus' Through Citrus Seed Taken from Affected Fruit

Citrus huanglongbing, putatively caused by the associated bacterium 'Candidatus Liberibacter asiaticus', is the greatest threat to the world citrus industry today. The bacterium is spread locally and regionally by the citrus psyllid Diaphorina citri, and also can be disseminated by propagation of contaminated scion budwood that is grafted to the appropriate rootstock. The planting of 'Ca. Liberibacter asiaticus'-free trees is a component of a comprehensive strategy to manage huanglongbing. In contrast to the scion budwood, the rootstocks used to produce these trees are grown from seed. This research was undertaken to provide evidence as to whether or not 'Ca. L. asiaticus' can be transmitted through seed. Two groups of 360 or more seedlings each of various citrus species were grown from seed removed from fruit on trees that were symptomatic for huanglongbing and confirmed to be infected with 'Ca. L. asiaticus' by polymerase chain reaction (PCR) tests. These seedlings were tested multiple times over periods of up to 3 years. No symptoms typical of huanglongbing, such as blotchy leaf mottle, chlorotic shoots, or dieback of branches, were observed in these seedlings, and none of these 723 seedlings tested positive for the presence of 'Ca. L. asiaticus' even after repeated testing by sensitive quantitative PCR assays. Some sour orange seedlings did have quite pronounced and atypical growth, including stunting and mild to severe leaf malformation. These atypical growth habits were limited to seedlings that arose from zygotic embryos as determined by expressed-sequence tag simple-sequence repeat analyses. Thus, no evidence of transmission of 'Ca. L. asiaticus' via seed was obtained, and an earlier report of transmission of the pathogen through seed was not confirmed.

Analysis of expressed sequence tags from Prunus mume flower and fruit and development of simple sequence repeat markers

Background

Expressed Sequence Tag (EST) has been a cost-effective tool in molecular biology and represents an abundant valuable resource for genome annotation, gene expression, and comparative genomics in plants.

Results

In this study, we constructed a cDNA library of Prunus mume flower and fruit, sequenced 10,123 clones of the library, and obtained 8,656 expressed sequence tag (EST) sequences with high quality. The ESTs were assembled into 4,473 unigenes composed of 1,492 contigs and 2,981 singletons and that have been deposited in NCBI (accession IDs: GW868575 - GW873047), among which 1,294 unique ESTs were with known or putative functions. Furthermore, we found 1,233 putative simple sequence repeats (SSRs) in the P. mume unigene dataset. We randomly tested 42 pairs of PCR primers flanking potential SSRs, and 14 pairs were identified as true-to-type SSR loci and could amplify polymorphic bands from 20 individual plants of P. mume. We further used the 14 EST-SSR primer pairs to test the transferability on peach and plum. The result showed that nearly 89% of the primer pairs produced target PCR bands in the two species. A high level of marker polymorphism was observed in the plum species (65%) and low in the peach (46%), and the clustering analysis of the three species indicated that these SSR markers were useful in the evaluation of genetic relationships and diversity between and within the Prunus species.

Conclusions

We have constructed the first cDNA library of P. mume flower and fruit, and our data provide sets of molecular biology resources for P. mume and other Prunus species. These resources will be useful for further study such as genome annotation, new gene discovery, gene functional analysis, molecular breeding, evolution and comparative genomics between Prunus species.

FG. Citrus genomics

Citrus is one of the most widespread fruit crops globally, with great economic and health value. It is among the most difficult plants to improve through traditional breeding approaches. Currently, there is risk of devastation by diseases threatening to limit production and future availability to the human population. As technologies rapidly advance in genomic science, they are quickly adapted to address the biological challenges of the citrus plant system and the world's industries. The historical developments of linkage mapping, markers and breeding, EST projects, physical mapping, an international citrus genome sequencing project, and critical functional analysis are described. Despite the challenges of working with citrus, there has been substantial progress. Citrus researchers engaged in international collaborations provide optimism about future productivity and contributions to the benefit of citrus industries worldwide and to the human population who can rely on future widespread availability of this health-promoting and aesthetically pleasing fruit crop.