Comparative use of InDel and SSR markers in deciphering the interspecific structure of cultivated citrus genetic diversity: a perspective for genetic association studies

Comparing Apples and Oranges: Advances in Disease Resistance Breeding of Woody Perennial Fruit Crops

Apple and citrus are perennial tree fruit crops that are vital for nutritional security and agricultural economy and to achieve the Sustainable Development Goals of the United Nations. Apple scab and fire blight, along with Huanglongbing, canker, and tristeza virus, stand out as their most notorious diseases and annually destabilize fruit supply. An environmentally sound approach to managing these diseases is improving tree resistance through breeding and biotechnology. Perennial fruit tree germplasm collections are distributed globally and offer untapped potential as sources of resistance. However, long juvenility, specific pollination and flowering habits, and extensive outcrossing hinder apple and citrus breeding. Advances in breeding approaches include trans- and cis-genesis, genome editing, and rapid-cycle breeding, which, in addition to conventional crossbreeding, can all facilitate accelerated integration of resistance into elite germplasm. In addition, the global pool of available sources of resistance can be characterized by the existing genetic mapping and gene expression studies for accurate discovery of associated loci, genes, and markers to efficiently include these sources in breeding efforts. We discuss and propose a multitude of approaches to overcome the challenges of breeding for resistance in woody perennials and outline a technical path to reduce the time required for the ultimate deployment of disease-resistant cultivars.

A Genomic Evaluation of Six Selected Inbred Lines of the Naturalized Plants of Milk Thistle (Silybum marianum L. Gaertn.) in Korea

Milk thistle (Silybum marianum) belongs to the Asteraceae family and is a medicinal plant native to the Mediterranean Basin. Silymarin in achene is a widely used herbal product for chronic liver disease. There is growing interest in natural medicine using milk thistle in Korea, but the raw material completely relies on imports. Despite its economic importance, phenotypic evaluations of native resources of milk thistle in Korea have not been carried out. In addition, genomic research and molecular marker development are very limited in milk thistle. In this study, we evaluated 220 milk thistle resources consisting of 172 accessions collected from the domestic market, and 48 accessions isolated from 6 accessions distributed by the National Agrobiodiversity Center in Korea. Six plant characteristics (height, seed weight, number of flowers, seed weight per flower, spine length, and color at harvest) were measured, and six samples (M01-M06) were selected to represent the genetic diversity of the population for genomic research. To develop PCR-based and co-dominant insertion/deletion (InDel) markers, we performed genome-wide InDel detection by comparing the whole-genome resequencing data of the six selected accessions with the reference genome sequence (GCA_001541825). As a result, 177 InDel markers with high distinguishability and reproducibility were selected from the 30,845 InDel variants. Unknowingly imported alien plant resources could easily be genetically mixed, and jeopardized seed purity can cause continuous difficulties in the development of high value-added agricultural platforms utilizing natural products. The selected plant materials and 177 validated InDel markers developed via whole-genome resequencing analysis could be valuable resources for breeding, conservation, and ecological studies of natives to Korea, along with acceleration of Silybum marianum industrialization.

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.

InDel and SCoT Markers for Genetic Diversity Analysis in a Citrus Collection from the Western Caucasus

Citrus collections from extreme growing regions can be an important source of tolerant germplasms for the breeding of cold-tolerant varieties. However, the efficient utilization of these germplasms requires their genetic background information. Thus, efficient marker systems are necessary for the characterization and identification of valuable accessions. In this study, the efficiency of 36 SCoT markers and 60 InDel markers were evaluated as part of the broad citrus collection of the Western Caucasus. The interspecific and intraspecific genetic diversity and genetic structures were analyzed for 172 accessions, including 31 species and sets of the locally derived cultivars. Single markers, such as SCoT18 (0.84), SCoT20 (0.93), SCoT23 (0.87), SCoT31 (0.88), SCoT36 (0.87) и LG 1-4 (0.94), LG 4-3 (0.86), LG 7-11 (0.98), and LG 8-10 (0.83), showed a high discriminating power, indicating the good applicability of these markers to assess intraspecific diversity of the genus Citrus. Overall, SCoT markers showed a higher level of polymorphism than InDel markers. According to analysis of population structure, SCoT and InDel markers showed K = 9 and K = 5 genetic clusters, respectively. The lowest levels of genetic admixtures and diversity were observed among the locally derived satsumas and lemons. The highest level of genetic admixtures was observed in the lime group. Phylogenetic relationships indicated a high level of interspecific genetic diversity but a low level of intraspecific diversity in locally derived satsumas and lemons. The results provide new insight into the origin of citrus germplasms and their distribution in colder regions. Furthermore, they are important for implementing conservation measures, controlling genetic erosion, developing breeding strategies, and improving breeding efficiency.

Correspondence between the Compositional and Aromatic Diversity of Leaf and Fruit Essential Oils and the Pomological Diversity of 43 Sweet Oranges (Citrus x aurantium var sinensis L.)

Orange (Citrus x aurantium var sinensis) is the most widely consumed citrus fruit, and its essential oil, which is made from the peel, is the most widely used in the food, perfume, and cosmetics industries. This citrus fruit is an interspecific hybrid that would have appeared long before our era and would result from two natural crosses between mandarin and pummelo hybrids. This single initial genotype was multiplied by apomictic reproduction and diversified by mutations to produce hundreds of cultivars selected by men essentially based on phenotypic characteristics of appearance, spread of maturity, and taste. Our study aimed to assess the diversity of essential oil composition and variability in the aroma profile of 43 orange cultivars representing all morphotypes. In agreement with the mutation-based evolution of orange trees, the genetic variability tested with 10 SSR genetic markers was null. The oils from peels and leaves extracted by hydrodistillation were analyzed for composition by GC (FID) and GC/MS and for aroma profile by the CATA (Check All That Apply) method by panelists. Oil yield varied between varieties by a factor of 3 for PEO and a factor of 14 for LEO between maximum and minimum. The composition of the oils was very similar between cultivars and was mainly dominated by limonene (>90%). However, small variations were observed as well as in the aromatic profile, with some varieties clearly distinguishing themselves from the others. This low chemical diversity contrasts with the pomological diversity, suggesting that aromatic variability has never been a selection criterion in orange trees.

Characterization of Odor-Active Compounds of Ichang Lemon (Citrus wilsonii Tan.) and Identification of Its Genetic Interspecific Origin by DNA Genotyping

Ichang lemon is a citrus fruit whose rind gives off a delicious and much appreciated fragrance and flavor. The volatile components of the fruit peel of Ichang lemon were investigated by GC-MS and GC-O (AEDA method). Simultaneously, its genetic origin was identified by using diagnostic SNP markers specific to ancestral species and multiallelic SSR and InDel markers. Ichang lemon combines three ancestral genomes (Citrus maxima, Citrus ichangensis, and Citrus reticulata) and may be a pummelo × Yuzu hybrid. Although the major compounds of the Ichang lemon aromatic profile were present in Citrus junos, a few pummelo-specific compounds were also detected, such as indole and nootkatone, in agreement with its maternal lineage. 3-Methyl-3-sulfanylbutyl acetate, reported to occur in passion fruit and brewed coffee, was identified by GC-MS, GC-QTOF-MS, and GC-FTIR for the first time in citrus. This odor-active compound has a sulfurous, tropical fruity, green note.

Strategies to Produce Grapefruit-Like Citrus Varieties With a Low Furanocoumarin Content and Distinctive Flavonoid Profiles

Pummelos and hybrids, such as grapefruits, have high furanocoumarin and low flavonoid contents. Furanocoumarins interact negatively with certain drugs, while flavonoids are antioxidant compounds with health benefits. To obtain new grapefruit-like varieties with low furanocoumarin and high flavonoid contents, diploid and triploid hybrid populations from crosses between diploid and tetraploid "Clemenules" clementine and diploid "Pink" pummelo were recovered and analyzed. With regard to furanocoumarins, triploids produce less bergapten, bergamottin and 6,7-DHB than diploids. Regarding flavonoids, triploids yielded more eriocitrin, narirutin, hesperidin and neohesperidin than diploids, whereas no differences were observed in neoeriocitrin and naringin. These results indicate that, the strategy to recover triploid hybrids by 4x × 2x crosses is more appropriate than the recovery of diploid hybrids by 2x × 2x crosses for obtaining grapefruit-like varieties of citrus with lower furanocoumarin and higher flavonoid contents.

Allelic composition of carotenoid metabolic genes in 13 founders influences carotenoid composition in juice sac tissues of fruits among Japanese citrus breeding population

To enrich carotenoids, especially β-cryptoxanthin, in juice sac tissues of fruits via molecular breeding in citrus, allele mining was utilized to dissect allelic variation of carotenoid metabolic genes and identify an optimum allele on the target loci characterized by expression quantitative trait (eQTL) analysis. SNPs of target carotenoid metabolic genes in 13 founders of the Japanese citrus breeding population were explored using the SureSelect target enrichment method. An independent allele was determined based on the presence or absence of reliable SNPs, using trio analysis to confirm inheritability between parent and offspring. Among the 13 founders, there were 7 PSY alleles, 7 HYb alleles, 11 ZEP alleles, 5 NCED alleles, and 4 alleles for the eQTL that control the transcription levels of PDS and ZDS among the ancestral species, indicating that some founders acquired those alleles from them. The carotenoid composition data of 263 breeding pedigrees in juice sac tissues revealed that the phenotypic variance of carotenoid composition was similar to that in the 13 founders, whereas the mean of total carotenoid content increased. This increase in total carotenoid content correlated with the increase in either or both β-cryptoxanthin and violaxanthin in juice sac tissues. Bayesian statistical analysis between allelic composition of target genes and carotenoid composition in 263 breeding pedigrees indicated that PSY-a and ZEP-e alleles at PSY and ZEP loci had strong positive effects on increasing the total carotenoid content, including β-cryptoxanthin and violaxanthin, in juice sac tissues. Moreover, the pyramiding of these alleles also increased the β-cryptoxanthin content. Interestingly, the offset interaction between the alleles with increasing and decreasing effects on carotenoid content and the epistatic interaction among carotenoid metabolic genes were observed and these interactions complexed carotenoid profiles in breeding population. These results revealed that allele composition would highly influence the carotenoid composition in citrus fruits. The allelic genotype information for the examined carotenoid metabolic genes in major citrus varieties and the trio-tagged SNPs to discriminate the optimum alleles (PSY-a and ZEP-e) from the rest would promise citrus breeders carotenoid enrichment in fruit via molecular breeding.

Fast-track breeding system to introduce CTV resistance of trifoliate orange into citrus germplasm, by integrating early flowering transgenic plants with marker-assisted selection

BACKGROUND

Global warming will expand the range of new and invasive pathogens in orchards, and subsequently increase the risk of disease epidemics and economic losses. The development of new resistant plant varieties can help to reduce the impact of pathogens, however, the breeding speed can be extremely slow, due to the growth rates of the plants, and the availability of resistance genes. Citrus trees are suffering immense damage from serious diseases such as citrus canker (XCC), huanglongbing (HLB), and citrus tristeza virus (CTV). A fast-track breeding system, that aimed at shortening the duration for disease resistance breeding by incorporating the resistance genes from related species to commercial varieties, has been developed using the integration of precocious transgenic trifoliate orange with the overexpression of CiFT and MAS. It was applied here to incorporate CTV resistance of trifoliate orange into citrus germplasm.

RESULTS

One generation of backcrossed breeding, that would normally take at least 5 years, was achieved in a single year by fast-track breeding system. Linkage analysis using the corresponding DNA markers revealed that CTV resistance and T-DNA integrated regions were found in different linkage groups, and they were independently segregated in the BC progenies. The CTV resistant null segregants, in which the T-DNA integrated region was removed from their genome, were feasibly obtained by MAS in each generation of the BC progenies, and their CTV resistance was confirmed by immunological analysis. Several BC₃ null segregants, whose genetic backgrounds had been substituted into citrus germplasm, except for the haplotype block of CTV resistance, were successfully obtained. CGH and NGS analyses revealed that the T-DNA integrated region was safely segregated out in null segregants.

CONCLUSION

Fast-track breeding systems are expected to shorten the required breeding time by more than one-fifth in comparison with conventional cross breeding techniques. Using this system, we obtained BC₃-8, whose genetic background was successfully substituted except for the CTV resistance locus, and could be a novel mandarin breeding material. The fast-track breeding system will be useful to introduce important traits from related species to citrus germplasm while also drastically reducing the time required for breeding.

Male and female inheritance patterns in tetraploid 'Moncada' mandarin

KEY MESSAGE

Tetraploid `Moncada´ mandarin, used as male and female in interploidy hybridizations, displays mainly tetrasomic inheritance for most LGs, with slight variations according to the direction of the crossing. Triploid-breeding programs in citrus are key tool to develop seedless cultivars. Obtaining triploid citrus hybrids may be achieved through different strategies, such as the exploitation of female unreduced gamete in crosses between diploid parents and diploid by tetraploid sexual hybridizations, in which tetraploid genotypes can be used as male or female parents. Genetic configuration of triploid populations from interploid crosses greatly depends on the chromosomic segregation mode of the tetraploid parent used. Here, we have analyzed the inheritance of the tetraploid 'Moncada' mandarin and compared the genetic structures of the resulting gametes when used as male and as female parent. The preferential chromosome pairing rate is calculated from the parental heterozygosity restitution (PHR) of codominant molecular markers, indicating the proportion between disomic and tetrasomic segregation. Tetraploid 'Moncada' both as female and male parent largely exhibited tetrasomic segregation. However, as female parent, one linkage group (LG8) showed intermediate segregation with tendency towards tetrasomic inheritance, while another linkage group (LG4) evidenced a clear intermediate segregation. On the other hand, when used as male parent two linkage groups (LG5 and LG6) showed values that fit an intermediate inheritance model with tetrasomic tendency. Significant doubled reduction (DR) rates were observed in five linkage groups as female parent, and in six linkage groups as male parent. The new knowledge generated here will serve to define crossing strategies in citrus improvement programs to efficiently obtain new varieties of interest in the global fresh consumption market.

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.

Characterization of genome-wide genetic variations between two varieties of tea plant (Camellia sinensis) and development of InDel markers for genetic research

Background

Single nucleotide polymorphisms (SNPs) and insertions/deletions (InDels) are the major genetic variations and are distributed extensively across the whole plant genome. However, few studies of these variations have been conducted in the long-lived perennial tea plant.

Results

In this study, we investigated the genome-wide genetic variations between Camellia sinensis var. sinensis ‘Shuchazao’ and Camellia sinensis var. assamica ‘Yunkang 10’, identified 7,511,731 SNPs and 255,218 InDels based on their whole genome sequences, and we subsequently analyzed their distinct types and distribution patterns. A total of 48 InDel markers that yielded polymorphic and unambiguous fragments were developed when screening six tea cultivars. These markers were further deployed on 46 tea cultivars for transferability and genetic diversity analysis, exhibiting information with an average 4.02 of the number of alleles (Na) and 0.457 of polymorphism information content (PIC). The dendrogram showed that the phylogenetic relationships among these tea cultivars are highly consistent with their genetic backgrounds or original places. Interestingly, we observed that the catechin/caffeine contents between ‘Shuchazao’ and ‘Yunkang 10’ were significantly different, and a large number of SNPs/InDels were identified within catechin/caffeine biosynthesis-related genes.

Conclusion

The identified genome-wide genetic variations and newly-developed InDel markers will provide a valuable resource for tea plant genetic and genomic studies, especially the SNPs/InDels within catechin/caffeine biosynthesis-related genes, which may serve as pivotal candidates for elucidating the molecular mechanism governing catechin/caffeine biosynthesis.

Genetic, morphological and chemical investigations reveal the genetic origin of Pompia (C. medica tuberosa Risso & Poiteau) - An old endemic Sardinian citrus fruit

Citrus fruits have been introduced to the Mediterranean area from Asia for centuries and spontaneous crosses have generated several hybrid forms, some of which have had agricultural or industrial success while others have remained niche food or ornamental products, or have disappeared. Pompia (C. medica tuberosa Risso & Poiteau) is an old endemic citrus fruit from Sardinia of unknown genetic origin. Initial phenotypic and molecular characterizations revealed a high degree of similarity with lemon (C. limon (L.) Burm.) and citron (C. medica L.). To identify the ancestors of Pompia, 70 citrus species of the Citrus genus were genotyped with 36 codominant molecular markers (SSR and InDel) of nuclear and cytoplasmic genomes. Diversity analysis and allelic comparisons between each citrus species at each locus indicated that Pompia resembles lemon and limonette of Marrakech, i.e. the result of a cross between sour orange (C. aurantium L.) and citron, where citron was the pollinator. Two Italian citron varieties were identified as potential male parents, i.e. Diamante and Common Poncire. However, we were unable to differentiate varieties of sour oranges because varietal diversification in this horticultural group resulted from DNA sequence variations that SSR or InDel markers could not reveal. Rhob el Arsa and Poncire de Collioure were found to be two synonyms of Pompia. Pompia appeared to be equally distinct from citron, lemon and sour orange based on the overall analysis of the fruit, leaf and seed phenotype, and juice chemical composition. At the leaf level, the Pompia essential oil (EO) composition is close to that of citron whereas the zest is much closer to that of sour orange.

Phylogenetic Origin of Primary and Secondary Metabolic Pathway Genes Revealed by C. maxima and C. reticulata Diagnostic SNPs

Modern cultivated Citrus species and varieties result from interspecific hybridization between four ancestral taxa. Among them, Citrus maxima and Citrus reticulata, closely associated with the pummelo and mandarin horticultural groups, respectively, were particularly important as the progenitors of sour and sweet oranges (Citrus aurantium and Citrus sinensis), grapefruits (Citrus paradisi), and hybrid types resulting from modern breeding programs (tangors, tangelos, and orangelos). The differentiation between the four ancestral taxa and the phylogenomic structure of modern varieties widely drive the phenotypic diversity's organization. In particular, strong phenotypic differences exist in the coloration and sweetness and represent important criteria for breeders. In this context, focusing on the genes of the sugar, carotenoid, and chlorophyll biosynthesis pathways, the aim of this work was to develop a set of diagnostic single-nucleotide polymorphism (SNP) markers to distinguish the ancestral haplotypes of C. maxima and C. reticulata and to provide information at the intraspecific diversity level (within C. reticulata or C. maxima). In silico analysis allowed the identification of 3,347 SNPs from selected genes. Among them, 1,024 were detected as potential differentiation markers between C. reticulata and C. maxima. A total of 115 SNPs were successfully developed using a competitive PCR technology. Their transferability among all Citrus species and the true citrus genera was very good, with only 0.87% of missing data. The ancestral alleles of the SNPs were identified, and we validated the usefulness of the developed markers for tracing the ancestral haplotype in large germplasm collections and sexually recombined progeny issued from the C. reticulata/C. maxima admixture gene pool. These markers will pave the way for targeted association studies based on ancestral haplotypes.

Genotyping by sequencing can reveal the complex mosaic genomes in gene pools resulting from reticulate evolution: a case study in diploid and polyploid citrus

BACKGROUND AND AIMS

Reticulate evolution, coupled with reproductive features limiting further interspecific recombinations, results in admixed mosaics of large genomic fragments from the ancestral taxa. Whole-genome sequencing (WGS) data are powerful tools to decipher such complex genomes but still too costly to be used for large populations. The aim of this work was to develop an approach to infer phylogenomic structures in diploid, triploid and tetraploid individuals from sequencing data in reduced genome complexity libraries. The approach was applied to the cultivated Citrus gene pool resulting from reticulate evolution involving four ancestral taxa, C. maxima, C. medica, C. micrantha and C. reticulata.

METHODS

A genotyping by sequencing library was established with the restriction enzyme ApeKI applying one base (A) selection. Diagnostic single nucleotide polymorphisms (DSNPs) for the four ancestral taxa were mined in 29 representative varieties. A generic pipeline based on a maximum likelihood analysis of the number of read data was established to infer ancestral contributions along the genome of diploid, triploid and tetraploid individuals. The pipeline was applied to 48 diploid, four triploid and one tetraploid citrus accessions.

KEY RESULTS

Among 43 598 mined SNPs, we identified a set of 15 946 DSNPs covering the whole genome with a distribution similar to that of gene sequences. The set efficiently inferred the phylogenomic karyotype of the 53 analysed accessions, providing patterns for common accessions very close to that previously established using WGS data. The complex phylogenomic karyotypes of 21 cultivated citrus, including bergamot, triploid and tetraploid limes, were revealed for the first time.

CONCLUSIONS

The pipeline, available online, efficiently inferred the phylogenomic structures of diploid, triploid and tetraploid citrus. It will be useful for any species whose reproductive behaviour resulted in an interspecific mosaic of large genomic fragments. It can also be used for the first generations of interspecific breeding schemes.

Comparative transcriptomics analysis uncovers alternative splicing events and molecular markers in cabbage (Brassica oleracea L.)

Alternative splicing (AS) events were identified and verified in cabbage by comparative transcriptome analysis. The corresponding markers were developed and the germplasm resources were identified. Alternative splicing (AS) is a central regulatory mechanism that greatly contributes to plant gene expression and transcriptome diversity. A large body of evidence has shown that AS complexity is relevant for plant development, evolution, complexity, and adaptation. Both insertion/deletion (InDel) and single nucleotide polymorphism (SNP) are typically co-dominant inheritance markers and have abundant polymorphisms. These have been widely used for marker-assisted selection, genetic mapping, and germplasm identification in plants. However, little is known about the molecular mechanisms underlying AS events and the development of markers including SNP and InDel from the cabbage transcriptome. In this study, three cabbage transcriptome datasets were collected and aligned to the cabbage reference genome to analyze AS events and marker development. 31,524 AS events were identified from three cabbage genotypes, accounting for 20.8% of the total cabbage genes. Alternative 3' splice site donor (A3SS) was the most frequent type of the four main AS events in cabbage. 70,475 InDels and 706,269 SNPs were identified with average frequencies of 1 InDel/6.9 kb and 1 SNP/0.7 kb, respectively. 71,942 potential SSRs were identified in 53,129 assembled unigenes with a density of 1 SSR/6.8 kb. The ratio of SNPs with synonymous/non-synonymous mutations was 1:0.65. 142 InDels and 36 SNPs were randomly selected and validated via Sanger sequencing and polymorphism was found among 66.2% of the InDels and 78.6% of the SNPs. Furthermore, 35 informative InDel markers were successfully used for genetic diversity analysis on 36 cabbage accessions. These results facilitate understanding of the molecular regulation mechanism underlying AS events in cabbage. They also provide molecular marker resource data for genetic mapping construction and germplasm identification, and facilitate the genetic improvement of cabbage via breeding.

Facing Climate Change: Biotechnology of Iconic Mediterranean Woody Crops

The Mediterranean basin is especially sensitive to the adverse outcomes of climate change and especially to variations in rainfall patterns and the incidence of extremely high temperatures. These two concurring adverse environmental conditions will surely have a detrimental effect on crop performance and productivity that will be particularly severe on woody crops such as citrus, olive and grapevine that define the backbone of traditional Mediterranean agriculture. These woody species have been traditionally selected for traits such as improved fruit yield and quality or alteration in harvesting periods, leaving out traits related to plant field performance. This is currently a crucial aspect due to the progressive and imminent effects of global climate change. Although complete genome sequence exists for sweet orange (Citrus sinensis) and clementine (Citrus clementina), olive tree (Olea europaea) and grapevine (Vitis vinifera), the development of biotechnological tools to improve stress tolerance still relies on the study of the available genetic resources including interspecific hybrids, naturally occurring (or induced) polyploids and wild relatives under field conditions. To this respect, post-genomic era studies including transcriptomics, metabolomics and proteomics provide a wide and unbiased view of plant physiology and biochemistry under adverse environmental conditions that, along with high-throughput phenotyping, could contribute to the characterization of plant genotypes exhibiting physiological and/or genetic traits that are correlated to abiotic stress tolerance. The ultimate goal of precision agriculture is to improve crop productivity, in terms of yield and quality, making a sustainable use of land and water resources under adverse environmental conditions using all available biotechnological tools and high-throughput phenotyping. This review focuses on the current state-of-the-art of biotechnological tools such as high throughput -omics and phenotyping on grapevine, citrus and olive and their contribution to plant breeding programs.

In the name of the rose: a roadmap for rose research in the genome era

The recent completion of the rose genome sequence is not the end of a process, but rather a starting point that opens up a whole set of new and exciting activities. Next to a high-quality genome sequence other genomic tools have also become available for rose, including transcriptomics data, a high-density single-nucleotide polymorphism array and software to perform linkage and quantitative trait locus mapping in polyploids. Rose cultivars are highly heterogeneous and diverse. This vast diversity in cultivated roses can be explained through the genetic potential of the genus, introgressions from wild species into commercial tetraploid germplasm and the inimitable efforts of historical breeders. We can now investigate how this diversity can best be exploited and refined in future breeding work, given the rich molecular toolbox now available to the rose breeding community. This paper presents possible lines of research now that rose has entered the genomics era, and attempts to partially answer the question that arises after the completion of any draft genome sequence: 'Now that we have "the" genome, what's next?'. Having access to a genome sequence will allow both (fundamental) scientific and (applied) breeding-orientated questions to be addressed. We outline possible approaches for a number of these questions.

Assessing Ploidy Level Analysis and Single Pollen Genotyping of Diploid and Euploid Citrus Genotypes by Fluorescence-Activated Cell Sorting and Whole-Genome Amplification

Flow cytometry is widely used to determine genome size and ploidy level in plants. This technique, when coupled with fluorescence-activated cell sorting (FACS), whole genome amplification and genotyping (WGA), opens up new opportunities for genetic studies of individualized nuclei. This strategy was used to analyze the genetic composition of single pollen nuclei of different citrus species. The flow cytometry and microscope observations allowed us to differentiate the populations of pollen nuclei present in the diploid and euploid genotypes analyzed, showing that citrus has binuclear pollen. We have identified in the "CSO" tangor an additional nuclei population composed by the vegetative plus generative nuclei. Genotyping of this nuclei population revealed that vegetative and generative nuclei show the same genetic configuration. In addition, we have demonstrated the presence of unreduced gametes in the diploid genotype "Mexican lime." Genomic amplification is a robust method for haploid nuclei genotyping with several molecular markers, whereas in diploid nuclei using heterozygous markers showed a bias towards one of the two alleles, limiting the use of this tool in this type of nuclei. We further discuss the importance and applications of single pollen genotyping in citrus genetic studies.

Finding invisible quantitative trait loci with missing data

Evolutionary processes during plant polyploidization and speciation have led to extensive presence-absence variation (PAV) in crop genomes, and there is increasing evidence that PAV associates with important traits. Today, high-resolution genetic analysis in major crops frequently implements simple, cost-effective, high-throughput genotyping from single nucleotide polymorphism (SNP) hybridization arrays; however, these are normally not designed to distinguish PAV from failed SNP calls caused by hybridization artefacts. Here, we describe a strategy to recover valuable information from single nucleotide absence polymorphisms (SNaPs) by population-based quality filtering of SNP hybridization data to distinguish patterns associated with genuine deletions from those caused by technical failures. We reveal that including SNaPs in genetic analyses elucidate segregation of small to large-scale structural variants in nested association mapping populations of oilseed rape (Brassica napus), a recent polyploid crop with widespread structural variation. Including SNaP markers in genomewide association studies identified numerous quantitative trait loci, invisible using SNP markers alone, for resistance to two major fungal diseases of oilseed rape, Sclerotinia stem rot and blackleg disease. Our results indicate that PAV has a strong influence on quantitative disease resistance in B. napus and that SNaP analysis using cost-effective SNP array data can provide extensive added value from 'missing data'. This strategy might also be applicable for improving the precision of genetic mapping in many important crop species.

Origin of C. latifolia and C. aurantiifolia triploid limes: the preferential disomic inheritance of doubled-diploid 'Mexican' lime is consistent with an interploid hybridization hypothesis

Background and Aims

Two main types of triploid limes are produced worldwide. The 'Tahiti' lime type (Citrus latifolia) is predominant, while the 'Tanepao' type (C. aurantiifolia) is produced to a lesser extent. Both types result from natural interspecific hybridization involving a diploid gamete of C. aurantiifolia 'Mexican' lime type (itself a direct interspecific C. micrantha × C. medica hybrid). The meiotic behaviour of a doubled-diploid 'Mexican' lime, the interspecific micrantha/medica recombination and the resulting diploid gamete structures were analysed to investigate the possibility that 'Tahiti' and 'Tanepao' varieties are derived from natural interploid hybridization.

Methods

A population of 85 tetraploid hybrids was established between a doubled-diploid clementine and a doubled-diploid 'Mexican' lime and used to infer the genotypes of 'Mexican' lime diploid gametes. Meiotic behaviour was studied through combined segregation analysis of 35 simple sequenbce repeat (SSR) and single nucleotide polymorphismn (SNP) markers covering the nine citrus chromosomes and cytogenetic studies. It was supplemented by pollen viability assessment.

Key Results

Pollen viability of the doubled-diploid Mexican lime (64 %) was much higher than that of the diploid. On average, 65 % of the chromosomes paired as bivalents and 31.4 % as tetravalents. Parental heterozygosity restitution ranged from 83 to 99 %. Disomic inheritance with high preferential pairing values was deduced for three chromosomes. Intermediate inheritances, with disomic trend, were found for five chromosomes, and an intermediate inheritance was observed for one chromosome. The average effective interspecific recombination rate was low (1.2 cM Mb-1).

Conclusion

The doubled-diploid 'Mexican' lime had predominantly disomic segregation, producing interspecific diploid gamete structures with high C. medica/C. micrantha heterozygosity, compatible with the phylogenomic structures of triploid C. latifolia and C. aurantiifolia varieties. This disomic trend limits effective interspecific recombination and diversity of the diploid gamete population. Interploid reconstruction breeding using doubled-diploid lime as one parent is a promising approach for triploid lime diversification.

Genotyping by sequencing reveals the interspecific C. maxima / C. reticulata admixture along the genomes of modern citrus varieties of mandarins, tangors, tangelos, orangelos and grapefruits

The mandarin horticultural group is an important component of world citrus production for the fresh fruit market. This group formerly classified as C. reticulata is highly polymorphic and recent molecular studies have suggested that numerous cultivated mandarins were introgressed by C. maxima (the pummelos). C. maxima and C. reticulata are also the ancestors of sweet and sour oranges, grapefruit, and therefore of all the "small citrus" modern varieties (mandarins, tangors, tangelos) derived from sexual hybridization between these horticultural groups. Recently, NGS technologies have greatly modified how plant evolution and genomic structure are analyzed, moving from phylogenetics to phylogenomics. The objective of this work was to develop a workflow for phylogenomic inference from Genotyping By Sequencing (GBS) data and to analyze the interspecific admixture along the nine citrus chromosomes for horticultural groups and recent varieties resulting from the combination of the C. reticulata and C. maxima gene pools. A GBS library was established from 55 citrus varieties, using the ApekI restriction enzyme and selective PCR to improve the read depth. Diagnostic polymorphisms (DPs) of C. reticulata/C. maxima differentiation were identified and used to decipher the phylogenomic structure of the 55 varieties. The GBS approach was powerful and revealed 30,289 SNPs and 8,794 Indels with 12.6% of missing data. 11,133 DPs were selected covering the nine chromosomes with a higher density in genic regions. GBS combined with the detection of DPs was powerful for deciphering the "phylogenomic karyotypes" of cultivars derived from admixture of the two ancestral species after a limited number of interspecific recombinations. All the mandarins, mandarin hybrids, tangelos and tangors analyzed displayed introgression of C. maxima in different parts of the genome. C. reticulata/C. maxima admixture should be a major component of the high phenotypic variability of this germplasm opening up the way for association studies based on phylogenomics.

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.

Comparative transcriptomics uncovers alternative splicing and molecular marker development in radish (Raphanus sativus L.)

Background

Alternative splicing (AS) plays important roles in gene expression and proteome diversity. Single nucleotide polymorphism (SNP) and insertion/deletion (InDel) are abundant polymorphisms and co-dominant inheritance markers, which have been widely used in germplasm identification, genetic mapping and marker-assisted selection in plants. So far, however, little information is available on utilization of AS events and development of SNP and InDel markers from transcriptome in radish.

Results

In this study, three radish transcriptome datasets were collected and aligned to the reference radish genome. A total of 56,530 AS events were identified from three radish genotypes with intron retention (IR) being the most frequent AS type, which accounted for 59.4% of the total expressed genes in radish. In all, 22,412 SNPs and 9436 InDels were identified with an average frequency of 1 SNP/17.9 kb and 1 InDel/42.5 kb, respectively. A total of 43,680 potential SSRs were identified in 31,604 assembled unigenes with a density of 1 SSR/2.5 kb. The ratio of SNPs with nonsynonymous/synonymous mutations was 1.05:1. Moreover, 35 SNPs and 200 InDels were randomly selected and validated by Sanger sequencing, 83.9% of the SNPs and 70% of the InDels exhibited polymorphism among these three genotypes. In addition, the 15 SNPs and 125 InDels were found to be unevenly distributed on 9 linkage groups. Furthermore, 40 informative InDel markers were successfully used for the genetic diversity analysis on 32 radish accessions.

Conclusions

These results would not only provide new insights into transcriptome complexity and AS regulation, but also furnish large amount of molecular marker resources for germplasm identification, genetic mapping and further genetic improvement of radish in breeding programs.

Electronic supplementary material

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

Chemistry and health effects of furanocoumarins in grapefruit

Furanocoumarins are a specific group of secondary metabolites that commonly present in higher plants, such as citrus plants. The major furanocoumarins found in grapefruits (Citrus paradisi) include bergamottin, epoxybergamottin, and 6′,7′-dihydroxybergamottin. During biosynthesis of these furanocoumarins, coumarins undergo biochemical modifications corresponding to a prenylation reaction catalyzed by the cytochrome P450 enzymes with the subsequent formation of furan rings. Because of undesirable interactions with several medications, many studies have developed methods for grapefruit furanocoumarin quantification that include high-performance liquid chromatography coupled with UV detector or mass spectrometry. The distribution of furanocoumarins in grapefruits is affected by several environmental conditions, such as processing techniques, storage temperature, and packing materials. In the past few years, grapefruit furanocoumarins have been demonstrated to exhibit several biological activities including antioxidative, -inflammatory, and -cancer activities as well as bone health promotion both in vitro and in vivo. Notably, furanocoumarins potently exerted antiproliferative activities against cancer cell growth through modulation of several molecular pathways, such as regulation of the signal transducer and activator of transcription 3, nuclear factor-κB, phosphatidylinositol-3-kinase/AKT, and mitogen-activated protein kinase expression. Therefore, based on this review, we suggest furanocoumarins may serve as bioactive components that contribute, at least in part, to the health benefits of grapefruit.

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.

Towards a molecular taxonomic key of the Aurantioideae subfamily using chloroplastic SNP diagnostic markers of the main clades genotyped by competitive allele-specific PCR

Background

Chloroplast DNA is a primary source of molecular variations for phylogenetic analysis of photosynthetic eukaryotes. However, the sequencing and analysis of multiple chloroplastic regions is difficult to apply to large collections or large samples of natural populations. The objective of our work was to demonstrate that a molecular taxonomic key based on easy, scalable and low-cost genotyping method should be developed from a set of Single Nucleotide Polymorphisms (SNPs) diagnostic of well-established clades. It was applied to the Aurantioideae subfamily, the largest group of the Rutaceae family that includes the cultivated citrus species.

Results

The publicly available nucleotide sequences of eight plastid genomic regions were compared for 79 accessions of the Aurantioideae subfamily to search for SNPs revealing taxonomic differentiation at the inter-tribe, inter-subtribe, inter-genus and interspecific levels. Diagnostic SNPs (DSNPs) were found for 46 of the 54 clade levels analysed. Forty DSNPs were selected to develop KASPar markers and their taxonomic value was tested by genotyping 108 accessions of the Aurantioideae subfamily. Twenty-seven markers diagnostic of 24 clades were validated and they displayed a very high rate of transferability in the Aurantioideae subfamily (only 1.2 % of missing data on average). The UPGMA from the validated markers produced a cladistic organisation that was highly coherent with the previous phylogenetic analysis based on the sequence data of the eight plasmid regions. In particular, the monophyletic origin of the “true citrus” genera plus Oxanthera was validated. However, some clarification remains necessary regarding the organisation of the other wild species of the Citreae tribe.

Conclusions

We validated the concept that with well-established clades, DSNPs can be selected and efficiently transformed into competitive allele-specific PCR markers (KASPar method) allowing cost-effective highly efficient cladistic analysis in large collections at subfamily level. The robustness of this genotyping method is an additional decisive advantage for network collaborative research. The availability of WGS data for the main “true citrus” species should soon make it possible to develop a set of DSNP markers allowing very fine resolution of this very important horticultural group.

Electronic supplementary material

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

Inheritance in doubled-diploid clementine and comparative study with SDR unreduced gametes of diploid clementine

Tetraploid clementine displays mainly tetrasomic inheritance. Genetic structures of 2n SDR and 2 × gametes from DD clementine are complementary and will guides triploids citrus breeding strategies. Triploid breeding is developed worldwide to create new seedless cultivars. Citrus triploid hybrids can be recovered from 2x × 2x sexual hybridizations as a consequence of the formation of unreduced gametes (2n), or from 4x × 2x interploid hybridizations in which tetraploid parents used are most often doubled-diploid (DD). Here we have analyzed the inheritance in doubled-diploid clementine and compared the genetic structures of gametes of DD clementine with SDR unreduced gametes of diploid clementine. Parental heterozygosity restitution (PHR) with DD parents depends on the rate of preferential chromosome pairing and thus the proportion of disomic versus tetrasomic segregations. Doubled-diploid clementine largely exhibited tetrasomic segregation. However, three linkage groups had intermediate segregation and one had a tendency for disomy. Significant doubled reduction rates (DR) rates were observed in six of the nine LGs. Differences of PHR between 2n SDR and 2x DD gametes were highest in the centromeric region and progressively decreased toward the distal regions where they were not significant. Over all markers, PHR was lower (two-thirds) in SDR 2n gametes than in DD-derived diploid gametes. The two strategies appear complementary in terms of genotypic variability. Interploid 4x × 2x hybridization is potentially more efficient for developing new cultivars that are phenotypically closer to the diploid parent of the DD than sexual hybridization through SDR 2n gametes. Conversely, 2x × 2x triploidisation has the potential to produce novel products with characteristics for market segmentation strategies.

Phylogenetic origin of limes and lemons revealed by cytoplasmic and nuclear markers

BACKGROUND AND AIMS

The origin of limes and lemons has been a source of conflicting taxonomic opinions. Biochemical studies, numerical taxonomy and recent molecular studies suggested that cultivated Citrus species result from interspecific hybridization between four basic taxa (C. reticulata,C. maxima,C. medica and C. micrantha). However, the origin of most lemons and limes remains controversial or unknown. The aim of this study was to perform extended analyses of the diversity, genetic structure and origin of limes and lemons.

METHODS

The study was based on 133 Citrus accessions. It combined maternal phylogeny studies based on mitochondrial and chloroplastic markers, and nuclear structure analysis based on the evaluation of ploidy level and the use of 123 markers, including 73 basic taxa diagnostic single nucleotide polymorphism (SNP) and indel markers.

KEY RESULTS

The lime and lemon horticultural group appears to be highly polymorphic, with diploid, triploid and tetraploid varieties, and to result from many independent reticulation events which defined the sub-groups. Maternal phylogeny involves four cytoplasmic types out of the six encountered in the Citrus genus. All lime and lemon accessions were highly heterozygous, with interspecific admixture of two, three and even the four ancestral taxa genomes. Molecular polymorphism between varieties of the same sub-group was very low.

CONCLUSIONS

Citrus medica contributed to all limes and lemons and was the direct male parent for the main sub-groups in combination with C. micrantha or close papeda species (for C. aurata, C. excelsa, C. macrophylla and C. aurantifolia--'Mexican' lime types of Tanaka's taxa), C. reticulata(for C. limonia, C. karna and C. jambhiri varieties of Tanaka's taxa, including popular citrus rootstocks such as 'Rangpur' lime, 'Volkamer' and 'Rough' lemons), C. aurantium (for C. limetta and C. limon--yellow lemon types--varieties of Tanaka's taxa) or the C. maxima × C. reticulate hybrid (for C. limettioides--'Palestine sweet' lime types--and C. meyeri). Among triploid limes, C. latifolia accessions ('Tahiti' and 'Persian' lime types) result from the fertilization of a haploid ovule of C. limon by a diploid gamete of C. aurantifolia, while C. aurantifolia triploid accessions ('Tanepao' lime types and 'Madagascar' lemon) probably result from an interspecific backcross (a diploid ovule of C. aurantifolia fertilized by C. medica). As limes and lemons were vegetatively propagated (apomixis, horticultural practices) the intra-sub-group phenotypic diversity results from asexual variations.

Development of Genome-Wide Insertion and Deletion Polymorphism Markers from Next-Generation Sequencing Data in Rice

BACKGROUND

Next-generation sequencing technologies enable the re-sequencing of a large number of genomes and provide an unprecedented opportunity to discover numerous DNA polymorphisms throughout the genome of a species. As the second most abundant form of genetic variation, InDels, with characteristics of co-dominance, multiple alleles and high stability and density and that are easy to genotype, have received an increasing amount attention.

RESULTS

In this work, a total of 2,329,544 InDels were identified in 1767 rice genomes; these InDels were dispersed across all 12 rice chromosomes, with one InDel marker found, on average, every 160.22 bp. There were 162,380 highly polymorphic InDels with a polymorphism information content (PIC) ≥ 0.5, contributing 1.81 % to the unique primer set. Of these highly polymorphic InDels, we also selected InDels with major allele differences (the size difference between the most and second most frequent alleles) ≥ 3 bp or 8 bp for primer design, which provided a more flexible choice for researchers. Finally, we experimentally validated 100 highly polymorphic InDels for accuracy and polymorphism. The PCR results showed that the accuracy of the InDel markers was 95.70 %, while the average PIC value was 0.56, with a range of 0.19 to 0.78; the average allele number was 3.02, with a range of 2 to 5.

CONCLUSIONS

Our genome-wide and easily used InDel markers with high polymorphism and density in both cultivated and wild rice will undoubtedly have practical implications in rice marker-assisted breeding and will also meet the need of fine-scale genetic mapping in map-based rice gene cloning.

Rapid screening for citrus canker resistance employing pathogen-associated molecular pattern-triggered immunity responses

Citrus canker, caused by the bacterial pathogen Xanthomonas citri ssp. citri (Xcc), has been attributed to millions of dollars in loss or damage to commercial citrus crops in subtropical production areas of the world. Since identification of resistant plants is one of the most effective methods of disease management, the ability to screen for resistant seedlings plays a key role in the production of a long-term solution to canker. Here, an inverse correlation between reactive oxygen species (ROS) production by the plant and the ability of Xcc to grow and form lesions on infected plants is reported. Based on this information, a novel screening method that can rapidly identify citrus seedlings that are less susceptible to early infection by Xcc was devised by measuring ROS accumulation triggered by a 22-amino acid sequence of the conserved N-terminal part of flagellin (flg22) from X. citri ssp. citri (Xcc-flg22). In addition to limiting disease symptoms, ROS production was also correlated with the expression of basal defense-related genes such as the pattern recognition receptors LRR8 and FLS2, the leucine-rich repeat receptor-like protein RLP12, and the defense-related gene PR1, indicating an important role for pathogen-associated molecular pattern-triggered immunity (PTI) in determining resistance to citrus canker. Moreover, the differential expression patterns observed amongst the citrus seedlings demonstrated the existence of genetic variations in the PTI response among citrus species/varieties.

Nuclear species-diagnostic SNP markers mined from 454 amplicon sequencing reveal admixture genomic structure of modern citrus varieties

Most cultivated Citrus species originated from interspecific hybridisation between four ancestral taxa (C. reticulata, C. maxima, C. medica, and C. micrantha) with limited further interspecific recombination due to vegetative propagation. This evolution resulted in admixture genomes with frequent interspecific heterozygosity. Moreover, a major part of the phenotypic diversity of edible citrus results from the initial differentiation between these taxa. Deciphering the phylogenomic structure of citrus germplasm is therefore essential for an efficient utilization of citrus biodiversity in breeding schemes. The objective of this work was to develop a set of species-diagnostic single nucleotide polymorphism (SNP) markers for the four Citrus ancestral taxa covering the nine chromosomes, and to use these markers to infer the phylogenomic structure of secondary species and modern cultivars. Species-diagnostic SNPs were mined from 454 amplicon sequencing of 57 gene fragments from 26 genotypes of the four basic taxa. Of the 1,053 SNPs mined from 28,507 kb sequence, 273 were found to be highly diagnostic for a single basic taxon. Species-diagnostic SNP markers (105) were used to analyse the admixture structure of varieties and rootstocks. This revealed C. maxima introgressions in most of the old and in all recent selections of mandarins, and suggested that C. reticulata × C. maxima reticulation and introgression processes were important in edible mandarin domestication. The large range of phylogenomic constitutions between C. reticulata and C. maxima revealed in mandarins, tangelos, tangors, sweet oranges, sour oranges, grapefruits, and orangelos is favourable for genetic association studies based on phylogenomic structures of the germplasm. Inferred admixture structures were in agreement with previous hypotheses regarding the origin of several secondary species and also revealed the probable origin of several acid citrus varieties. The developed species-diagnostic SNP marker set will be useful for systematic estimation of admixture structure of citrus germplasm and for diverse genetic studies.

A Phylogenetic Analysis of 34 Chloroplast Genomes Elucidates the Relationships between Wild and Domestic Species within the Genus Citrus

Citrus genus includes some of the most important cultivated fruit trees worldwide. Despite being extensively studied because of its commercial relevance, the origin of cultivated citrus species and the history of its domestication still remain an open question. Here, we present a phylogenetic analysis of the chloroplast genomes of 34 citrus genotypes which constitutes the most comprehensive and detailed study to date on the evolution and variability of the genus Citrus. A statistical model was used to estimate divergence times between the major citrus groups. Additionally, a complete map of the variability across the genome of different citrus species was produced, including single nucleotide variants, heteroplasmic positions, indels (insertions and deletions), and large structural variants. The distribution of all these variants provided further independent support to the phylogeny obtained. An unexpected finding was the high level of heteroplasmy found in several of the analyzed genomes. The use of the complete chloroplast DNA not only paves the way for a better understanding of the phylogenetic relationships within the Citrus genus but also provides original insights into other elusive evolutionary processes, such as chloroplast inheritance, heteroplasmy, and gene selection.

De novo transcriptome assembly of pummelo and molecular marker development

Pummelo (Citrus grandis) is an important fruit crop worldwide because of its nutritional value. To accelerate the pummelo breeding program, it is essential to obtain extensive genetic information and develop relative molecular markers. Here, we obtained a 12-Gb transcriptome dataset of pummelo through a mixture of RNA from seven tissues using Illumina pair-end sequencing, assembled into 57,212 unigenes with an average length of 1010 bp. The annotation and classification results showed that a total of 39,584 unigenes had similar hits to the known proteins of four public databases, and 31,501 were classified into 55 Gene Ontology (GO) functional sub-categories. The search for putative molecular markers among 57,212 unigenes identified 10,276 simple sequence repeats (SSRs) and 64,720 single nucleotide polymorphisms (SNPs). High-quality primers of 1174 SSR loci were designed, of which 88.16% were localized to nine chromosomes of sweet orange. Of 100 SSR primers that were randomly selected for testing, 87 successfully amplified clear banding patterns. Of these primers, 29 with a mean PIC (polymorphic information content) value of 0.52 were effectively applied for phylogenetic analysis. Of the 20 SNP primers, 14 primers, including 54 potential SNPs, yielded target amplifications, and 46 loci were verified via Sanger sequencing. This new dataset will be a valuable resource for molecular biology studies of pummelo and provides reliable information regarding SNP and SSR marker development, thus expediting the breeding program of pummelo.

Next generation haplotyping to decipher nuclear genomic interspecific admixture in Citrus species: analysis of chromosome 2

BACKGROUND

The most economically important Citrus species originated by natural interspecific hybridization between four ancestral taxa (Citrus reticulata, Citrus maxima, Citrus medica, and Citrus micrantha) and from limited subsequent interspecific recombination as a result of apomixis and vegetative propagation. Such reticulate evolution coupled with vegetative propagation results in mosaic genomes with large chromosome fragments from the basic taxa in frequent interspecific heterozygosity. Modern breeding of these species is hampered by their complex heterozygous genomic structures that determine species phenotype and are broken by sexual hybridisation. Nevertheless, a large amount of diversity is present in the citrus gene pool, and breeding to allow inclusion of desirable traits is of paramount importance. However, the efficient mobilization of citrus biodiversity in innovative breeding schemes requires previous understanding of Citrus origins and genomic structures. Haplotyping of multiple gene fragments along the whole genome is a powerful approach to reveal the admixture genomic structure of current species and to resolve the evolutionary history of the gene pools. In this study, the efficiency of parallel sequencing with 454 methodology to decipher the hybrid structure of modern citrus species was assessed by analysis of 16 gene fragments on chromosome 2.

RESULTS

454 amplicon libraries were established using the Fluidigm array system for 48 genotypes and 16 gene fragments from chromosome 2. Haplotypes were established from the reads of each accession and phylogenetic analyses were performed using the haplotypic data for each gene fragment. The length of 454 reads and the level of differentiation between the ancestral taxa of modern citrus allowed efficient haplotype phylogenetic assignations for 12 of the 16 gene fragments. The analysis of the mixed genomic structure of modern species and cultivars (i) revealed C. maxima introgressions in modern mandarins, (ii) was consistent with previous hypotheses regarding the origin of secondary species, and (iii) provided a new picture of the evolution of chromosome 2.

CONCLUSIONS

454 sequencing was an efficient strategy to establish haplotypes with significant phylogenetic assignations in Citrus, providing a new picture of the mixed structure on chromosome 2 in 48 citrus genotypes.

Exploration of presence/absence variation and corresponding polymorphic markers in soybean genome

This study was designed to reveal the genome-wide distribution of presence/absence variation (PAV) and to establish a database of polymorphic PAV markers in soybean. The 33 soybean whole-genome sequences were compared to each other with that of Williams 82 as a reference genome. A total of 33,127 PAVs were detected and 28,912 PAV markers with their primer sequences were designed as the database NJAUSoyPAV_1.0. The PAVs scattered on whole genome while only 518 (1.8%) overlapped with simple sequence repeats (SSRs) in BARCSOYSSR_1.0 database. In a random sample of 800 PAVs, 713 (89.13%) showed polymorphism among the 12 differential genotypes. Using 126 PAVs and 108 SSRs to test a Chinese soybean germplasm collection composed of 828 Glycine soja Sieb. et Zucc. and Glycine max (L.) Merr. accessions, the per locus allele number and its variation appeared less in PAVs than in SSRs. The distinctness among alleles/bands of PCR (polymerase chain reaction) products showed better in PAVs than in SSRs, potential in accurate marker-assisted allele selection. The association mapping results showed SSR + PAV was more powerful than any single marker systems. The NJAUSoyPAV_1.0 database has enriched the source of PCR markers, and may fit the materials with a range of per locus allele numbers, if jointly used with SSR markers.

Genome-wide characterization of insertion and deletion variation in chicken using next generation sequencing

Insertion and deletion (INDEL) is one of the main events contributing to genetic and phenotypic diversity, which receives less attention than SNP and large structural variation. To gain a better knowledge of INDEL variation in chicken genome, we applied next generation sequencing on 12 diverse chicken breeds at an average effective depth of 8.6. Over 1.3 million non-redundant short INDELs (1-49 bp) were obtained, the vast majority (92.48%) of which were novel. Follow-up validation assays confirmed that most (88.00%) of the randomly selected INDELs represent true variations. The majority (95.76%) of INDELs were less than 10 bp. Both the detected number and affected bases were larger for deletions than insertions. In total, INDELs covered 3.8 Mbp, corresponding to 0.36% of the chicken genome. The average genomic INDEL density was estimated as 0.49 per kb. INDELs were ubiquitous and distributed in a non-uniform fashion across chromosomes, with lower INDEL density in micro-chromosomes than in others, and some functional regions like exons and UTRs were prone to less INDELs than introns and intergenic regions. Nearly 620,253 INDELs fell in genic regions, 1,765 (0.28%) of which located in exons, spanning 1,358 (7.56%) unique Ensembl genes. Many of them are associated with economically important traits and some are the homologues of human disease-related genes. We demonstrate that sequencing multiple individuals at a medium depth offers a promising way for reliable identification of INDELs. The coding INDELs are valuable candidates for further elucidation of the association between genotypes and phenotypes. The chicken INDELs revealed by our study can be useful for future studies, including development of INDEL markers, construction of high density linkage map, INDEL arrays design, and hopefully, molecular breeding programs in chicken.

Utility of indels for species-level identification of a biologically complex plant group: a study with intergenic spacer in Citrus

The Consortium of Barcode of Life plant working group proposed to use the defined portion of plastid genes rbcL and matK either singly or in combination as the standard DNA barcode for plants. But DNA barcode based identification of biologically complex plant groups are always a challenging task due to the occurrence of natural hybridization. Here, we examined the use of indels polymorphism in trnH-psbA and trnL-trnF sequences for rapid species identification of citrus. DNA from young leaves of selected citrus species were isolated and matK gene (~800 bp) and trnH-psbA spacer (~450 bp) of Chloroplast DNA was amplified for species level identification. The sequences within the group taxa of Citrus were aligned using the ClustalX program. With few obvious misalignments were corrected manually using the similarity criterion. We identified a 54 bp inverted repeat or palindrome sequence (27-80 regions) and 6 multi residues indel coding regions. Large inverted repeats in cpDNA provided authentication at the higher taxonomic levels. These diagnostics indel marker from trnH-psbA were successful in identifying different species (5 out of 7) within the studied Citrus except Citrus limon and Citrus medica. These two closely related species are distinguished through the 6 bp deletion in trnL-trnF. This study demonstrated that the indel polymorphism based approach easily characterizes the Citrus species and the same may be applied in other complex groups. Likewise other indels occurring intergenic spacer of chloroplast regions may be tested for rapid identification of other secondary citrus species.

Genetic analysis and molecular characterization of Chinese sesame (Sesamum indicum L.) cultivars using insertion-deletion (InDel) and simple sequence repeat (SSR) markers

Background

Sesame is an important and ancient oil crop in tropical and subtropical areas. China is one of the most important sesame producing countries with many germplasm accessions and excellent cultivars. Domestication and modern plant breeding have presumably narrowed the genetic basis of cultivated sesame. Several modern sesame cultivars were bred with a limited number of landrace cultivars in their pedigree. The genetic variation was subsequently reduced by genetic drift and selection. Characterization of genetic diversity of these cultivars by molecular markers is of great value to assist parental line selection and breeding strategy design.

Results

Three hundred and forty nine simple sequence repeat (SSR) and 79 insertion-deletion (InDel) markers were developed from cDNA library and reduced-representation sequencing of a sesame cultivar Zhongzhi 14, respectively. Combined with previously published SSR markers, 88 polymorphic markers were used to assess the genetic diversity, phylogenetic relationships, population structure, and allele distribution among 130 Chinese sesame accessions including 82 cultivars, 44 landraces and 4 wild germplasm accessions. A total of 325 alleles were detected, with the average gene diversity of 0.432. Model-based structure analysis revealed the presence of five subgroups belonging to two main groups, which were consistent with the results from principal coordinate analysis (PCA), phylogenetic clustering and analysis of molecular variance (AMOVA). Several missing or unique alleles were identified from particular types, subgroups or families, even though they share one or both parental/progenitor lines.

Conclusions

This report presented a by far most comprehensive characterization of the molecular and genetic diversity of sesame cultivars in China. InDels are more polymorphic than SSRs, but their ability for deciphering genetic diversity compared to the later. Improved sesame cultivars have narrower genetic basis than landraces, reflecting the effect of genetic drift or selection during breeding processes. Comparative analysis of allele distribution revealed genetic divergence between improved cultivars and landraces, as well as between cultivars released in different years. These results will be useful for assessing cultivars and for marker-assisted breeding in sesame.

Citrus (Rutaceae) SNP markers based on Competitive Allele-Specific PCR; transferability across the Aurantioideae subfamily

PREMISE OF THE STUDY

Single nucleotide polymorphism (SNP) markers based on Competitive Allele-Specific PCR (KASPar) were developed from sequences of three Citrus species. Their transferability was tested in 63 Citrus genotypes and 19 relative genera of the subfamily Aurantioideae to estimate the potential of SNP markers, selected from a limited intrageneric discovery panel, for ongoing broader diversity analysis at the intra- and intergeneric levels and systematic germplasm bank characterization. •

METHODS AND RESULTS

Forty-two SNP markers were developed using KASPar technology. Forty-one were successfully genotyped in all of the Citrus germplasm, where intra- and interspecific polymorphisms were observed. The transferability and diversity decreased with increasing taxonomic distance. •

CONCLUSIONS

SNP markers based on the KASPar method developed from sequence data of a limited intrageneric discovery panel provide a valuable molecular resource for genetic diversity analysis of germplasm within a genus and should be useful for germplasm fingerprinting at a much broader diversity level.

A nuclear phylogenetic analysis: SNPs, indels and SSRs deliver new insights into the relationships in the 'true citrus fruit trees' group (Citrinae, Rutaceae) and the origin of cultivated species

BACKGROUND AND AIMS

Despite differences in morphology, the genera representing 'true citrus fruit trees' are sexually compatible, and their phylogenetic relationships remain unclear. Most of the important commercial 'species' of Citrus are believed to be of interspecific origin. By studying polymorphisms of 27 nuclear genes, the average molecular differentiation between species was estimated and some phylogenetic relationships between 'true citrus fruit trees' were clarified.

METHODS

Sanger sequencing of PCR-amplified fragments from 18 genes involved in metabolite biosynthesis pathways and nine putative genes for salt tolerance was performed for 45 genotypes of Citrus and relatives of Citrus to mine single nucleotide polymorphisms (SNPs) and indel polymorphisms. Fifty nuclear simple sequence repeats (SSRs) were also analysed.

KEY RESULTS

A total of 16 238 kb of DNA was sequenced for each genotype, and 1097 single nucleotide polymorphisms (SNPs) and 50 indels were identified. These polymorphisms were more valuable than SSRs for inter-taxon differentiation. Nuclear phylogenetic analysis revealed that Citrus reticulata and Fortunella form a cluster that is differentiated from the clade that includes three other basic taxa of cultivated citrus (C. maxima, C. medica and C. micrantha). These results confirm the taxonomic subdivision between the subgenera Metacitrus and Archicitrus. A few genes displayed positive selection patterns within or between species, but most of them displayed neutral patterns. The phylogenetic inheritance patterns of the analysed genes were inferred for commercial Citrus spp.

CONCLUSIONS

Numerous molecular polymorphisms (SNPs and indels), which are potentially useful for the analysis of interspecific genetic structures, have been identified. The nuclear phylogenetic network for Citrus and its sexually compatible relatives was consistent with the geographical origins of these genera. The positive selection observed for a few genes will help further works to analyse the molecular basis of the variability of the associated traits. This study presents new insights into the origin of C. sinensis.

A reference genetic map of C. clementina hort. ex Tan.; citrus evolution inferences from comparative mapping

BACKGROUND

Most modern citrus cultivars have an interspecific origin. As a foundational step towards deciphering the interspecific genome structures, a reference whole genome sequence was produced by the International Citrus Genome Consortium from a haploid derived from Clementine mandarin. The availability of a saturated genetic map of Clementine was identified as an essential prerequisite to assist the whole genome sequence assembly. Clementine is believed to be a 'Mediterranean' mandarin × sweet orange hybrid, and sweet orange likely arose from interspecific hybridizations between mandarin and pummelo gene pools. The primary goals of the present study were to establish a Clementine reference map using codominant markers, and to perform comparative mapping of pummelo, sweet orange, and Clementine.

RESULTS

Five parental genetic maps were established from three segregating populations, which were genotyped with Single Nucleotide Polymorphism (SNP), Simple Sequence Repeats (SSR) and Insertion-Deletion (Indel) markers. An initial medium density reference map (961 markers for 1084.1 cM) of the Clementine was established by combining male and female Clementine segregation data. This Clementine map was compared with two pummelo maps and a sweet orange map. The linear order of markers was highly conserved in the different species. However, significant differences in map size were observed, which suggests a variation in the recombination rates. Skewed segregations were much higher in the male than female Clementine mapping data. The mapping data confirmed that Clementine arose from hybridization between 'Mediterranean' mandarin and sweet orange. The results identified nine recombination break points for the sweet orange gamete that contributed to the Clementine genome.

CONCLUSIONS

A reference genetic map of citrus, used to facilitate the chromosome assembly of the first citrus reference genome sequence, was established. The high conservation of marker order observed at the interspecific level should allow reasonable inferences of most citrus genome sequences by mapping next-generation sequencing (NGS) data in the reference genome sequence. The genome of the haploid Clementine used to establish the citrus reference genome sequence appears to have been inherited primarily from the 'Mediterranean' mandarin. The high frequency of skewed allelic segregations in the male Clementine data underline the probable extent of deviation from Mendelian segregation for characters controlled by heterozygous loci in male parents.