Segregation distortion and linkage analysis in eggplant (Solanum melongena L.)

A Compendium for Novel Marker-Based Breeding Strategies in Eggplant

The worldwide production of eggplant is estimated at about 58 Mt, with China, India and Egypt being the major producing countries. Breeding efforts in the species have mainly focused on increasing productivity, abiotic and biotic tolerance/resistance, shelf-life, the content of health-promoting metabolites in the fruit rather than decreasing the content of anti-nutritional compounds in the fruit. From the literature, we collected information on mapping quantitative trait loci (QTLs) affecting eggplant's traits following a biparental or multi-parent approach as well as genome-wide association (GWA) studies. The positions of QTLs were lifted according to the eggplant reference line (v4.1) and more than 700 QTLs were identified, here organized into 180 quantitative genomic regions (QGRs). Our findings thus provide a tool to: (i) determine the best donor genotypes for specific traits; (ii) narrow down QTL regions affecting a trait by combining information from different populations; (iii) pinpoint potential candidate genes.

Genome-Wide Association Study for Agro-Morphological Traits in Eggplant Core Collection

Eggplant is one of the most economically and nutritionally important vegetables worldwide. The study of the association of phenotypic traits with genetic factors is vital for the rapid and efficient identification and selection of eggplant genetic resources for breeding purposes with desired traits. The eggplant resources (587) collected from different countries, including Korea, were used for establishing the core collection. A total of 288 accessions were selected from 587 Solanum accessions based on 52 single nucleotide polymorphisms (SNPs) markers together with 17 morphological traits. This core collection was further used to analyze the genetic associations of eggplant morphological variations. A large variation was found among the evaluated eggplant accessions for some agro-morphological traits. Stem prickles and leaf prickles showed a significant positive correlation (r = 0.83***), followed by days to flowering and days to maturity (r = 0.64***). A total of 114,981 SNPs were filtered and used for phylogenetic tree analysis, population structure analysis, and genome-wide association study (GWAS). Among the agro-morphological traits, significantly associated SNPs were found for six traits. A total of 377 significantly associated SNPs with six agro-morphological traits were identified. These six traits and the number of SNPs were: days to maturity (51), flower size (121), fruit width (20), harvest fruit color (42), leaf prickles (38), and stem prickles (105). The largest fraction of significant SNPs (11.94%) was obtained on chromosome Ch01, followed by Ch07 and Ch06 with 11.67% and 10.08%, respectively. This study will help to develop markers linked to the most important agro-morphological traits of eggplant genetic resources and support the selection of desirable traits for eggplant breeding programs.

A Genomic BSAseq Approach for the Characterization of QTLs Underlying Resistance to Fusarium oxysporum in Eggplant

Eggplant (Solanum melongena L.), similar to many other crops, suffers from soil-borne diseases, including Fusarium oxysporum f. sp. melongenae (Fom), causing wilting and heavy yield loss. To date, the genetic factors underlying plant responses to Fom are not well known. We previously developed a Recombinant Inbred Lines (RILs) population using as a female parent the fully resistant line ‘305E40’ and as a male parent the partially resistant line ‘67/3’. The fully resistant trait to Fom was introgressed from the allied species S. aethiopicum. In this work, the RIL population was assessed for the responses to Fom and by using a genomic mapping approach, two major QTLs on chromosomes CH02 and CH11 were identified, associated with the full and partial resistance trait to Fom, respectively. A targeted BSAseq procedure in which Illumina reads bulks of RILs grouped according to their resistance score was aligned to the appropriate reference genomes highlighted differentially enriched regions between resistant/susceptible progeny in the genomic regions underlying both QTLs. The characterization of such regions allowed us to identify the most reliable candidate genes for the two resistance traits. With the aim of revealing exclusive species-specific contigs and scaffolds inherited from the allied species and thus associated with the full resistance trait, a draft de-novo assembly of available Illumina sequences of the ‘305E40’ parent was developed to better resolve the non-recombining genomic region on its CH02 carrying the introgressed Fom resistance locus from S. aethiopicum.

Newly Developed MAGIC Population Allows Identification of Strong Associations and Candidate Genes for Anthocyanin Pigmentation in Eggplant

Multi-parent advanced generation inter-cross (MAGIC) populations facilitate the genetic dissection of complex quantitative traits in plants and are valuable breeding materials. We report the development of the first eggplant MAGIC population (S3 Magic EGGplant InCanum, S3MEGGIC; 8-way), constituted by the 420 S3 individuals developed from the intercrossing of seven cultivated eggplant (Solanum melongena) and one wild relative (S. incanum) parents. The S3MEGGIC recombinant population was genotyped with the eggplant 5k probes SPET platform and phenotyped for anthocyanin presence in vegetative plant tissues (PA) and fruit epidermis (FA), and for the light-insensitive anthocyanic pigmentation under the calyx (PUC). The 7,724 filtered high-confidence single-nucleotide polymorphisms (SNPs) confirmed a low residual heterozygosity (6.87%), a lack of genetic structure in the S3MEGGIC population, and no differentiation among subpopulations carrying a cultivated or wild cytoplasm. Inference of haplotype blocks of the nuclear genome revealed an unbalanced representation of the founder genomes, suggesting a cryptic selection in favour or against specific parental genomes. Genome-wide association study (GWAS) analysis for PA, FA, and PUC detected strong associations with two myeloblastosis (MYB) genes similar to MYB113 involved in the anthocyanin biosynthesis pathway, and with a COP1 gene which encodes for a photo-regulatory protein and may be responsible for the PUC trait. Evidence was found of a duplication of an ancestral MYB113 gene with a translocation from chromosome 10 to chromosome 1 compared with the tomato genome. Parental genotypes for the three genes were in agreement with the identification of the candidate genes performed in the S3MEGGIC population. Our new eggplant MAGIC population is the largest recombinant population in eggplant and is a powerful tool for eggplant genetics and breeding studies.

Doubled Haploids in Eggplant

Eggplant is a solanaceous crop cultivated worldwide for its edible fruit. Eggplant breeding programs are mainly aimed to the generation of F1 hybrids by crossing two highly homozygous, pure lines, which are traditionally obtained upon several self crossing generations, which is an expensive and time consuming process. Alternatively, fully homozygous, doubled haploid (DH) individuals can be induced from haploid cells of the germ line in a single generation. Several attempts have been made to develop protocols to produce eggplant DHs principally using anther culture and isolated microspore culture. Eggplant could be considered a moderately recalcitrant species in terms of ability for DH production. Anther culture stands nowadays as the most valuable technology to obtain eggplant DHs. However, the theoretical possibility of having plants regenerated from somatic tissues of the anther walls cannot be ruled out. For this reason, the use of isolated microspores is recommended when possible. This approach still has room for improvement, but it is largely genotype-dependent. In this review, we compile the most relevant advances made in DH production in eggplant, their application to breeding programs, and the future perspectives for the development of other, less genotype-dependent, DH technologies.

An Eggplant Recombinant Inbred Population Allows the Discovery of Metabolic QTLs Controlling Fruit Nutritional Quality

Eggplant (Solanum melongena L.) represents the third most important crop of the Solanaceae family and is an important component of our daily diet. A population of 164 F6 recombinant inbred lines (RILs), derived from two eggplant lines differing with respect to several key agronomic traits, "305E40" and "67/3," was grown to the commercial maturation stage, and fruits were harvested, separated into peel and flesh, and subjected to liquid chromatography Liquid Chromatography/Mass Spectrometry (LC/MS) analysis. Through a combination of untargeted and targeted metabolomics approaches, a number of metabolites belonging to the glycoalkaloid, anthocyanin, and polyamine classes and showing a differential accumulation in the two parental lines and F1 hybrid were identified. Through metabolic profiling of the RILs, we identified several metabolomic quantitative trait loci (mQTLs) associated with the accumulation of those metabolites. Each of the metabolic traits proved to be controlled by one or more quantitative trait loci (QTLs); for most of the traits, one major mQTL (phenotypic variation explained [PVE] ≥ 10%) was identified. Data on mQTL mapping and dominance-recessivity relationships of measured compounds in the parental lines and F1 hybrid, as well as an analysis of the candidate genes underlying the QTLs and of their sequence differences in the two parental lines, suggested a series of candidate genes underlying the traits under study.

Use of Molecular Markers for Doubled Haploid Technology: From Academia to Plant Breeding Companies

Molecular markers are employed for doubled haploid (DH) technology by researchers and applied plant breeders in many crops. In the 1990s, isozymes and RFLPs were commonly used marker technologies to characterize DHs and were later replaced by PCR- based markers (e.g., RAPDs, AFLPs, ISSRs, SSRs) and today by SNPs. Markers are used for multiple purposes in DH production, that is, for the study of genes underlying haploid induction and confirming homozygous plants of gametophytic origin. Furthermore, they are tools for investigating segregation in DH populations and for mapping simple and complex traits using DHs. The deployment of DHs and markers for developing trait-linked markers are demonstrated with examples from rapeseed, wheat, and barley. Marker development for resistance to viruses derived from genetic resources and their use in, for example, pyramiding of resistance genes, are given as an example for the combination of DH-technology and marker development in research. Today, marker systems amenable to automation are frequently used in applied plant breeding. Practical examples are given from Lantmännen (LM) ( https://Lantmannen.com ) using large-scale genotyping for variety development based on SSRs and SNPs.

Identification of Quantitative Trait Loci Controlling the Development of Prickles in Eggplant by Genome Re-sequencing Analysis

Eggplant (Solanum melongena L.) is the third most important crop in the family of Solanaceae. Prickles are considered as the undesirable traits during the plantation of eggplant and the transportation of fruits. In this study, we constructed a high-quality genetic linkage Bin map derived from the re-sequencing analysis on a cross of a prickly wild landrace, 17C01, and a cultivated variety, 17C02. The major quantitative trait locus (QTL) controlling the development of prickles on the calyx (explained 30.42% of the phenotypic variation), named as qPC.12, was identified on a ~7 kb region on chromosome 12. A gene within qPC.12, which encodes a WUSCHEL-related homeobox-like protein, with higher expression levels in 17C01 calyx and 22-bp deletion in 17C02 was probably the functional gene for prickle formation. Results from this study would ultimately facilitate uncovering the molecular regulatory mechanisms underlying the development of a prickle in eggplant.

Fine mapping of a major locus representing the lack of prickles in eggplant revealed the availability of a 0.5-kb insertion/deletion for marker-assisted selection

As prickles cause labour inefficiency during cultivation and scratches on the skin of fruits during transportation, they are considered undesirable traits of eggplant (Solanum melongena L.). Because the molecular basis of prickle emergence has not been entirely revealed in plants, we mapped an eggplant semi-dominant Prickle (Pl) gene locus, which causes the absence of prickles, on chromosome 6 of a linkage map of the F₂ population derived from crossing the no-prickly cultivar 'Togenashi-senryo-nigo' and the prickly line LS1934. By performing synteny mapping with tomato, the genomic region corresponding to the eggplant Pl locus was identified. Through bacterial artificial chromosome (BAC) screening, positive BAC clones and the contig sequence that harbour the Pl locus in the prickly eggplant genome were revealed. The BAC contig length was 133 kb, and it contained 16 predicted genes. Among them, a characteristic 0.5-kb insertion/deletion was detected. As the 0.5-kb insertion was commonly identified with the prickly phenotype worldwide, a primer pair that amplifies the insertion/deletion could be used for marker-assisted selection of the no-prickly phenotype. Such findings contribute to map-based-cloning of the Pl gene and the understanding of gene function, ultimately providing new insights into the regulatory molecular mechanisms underlying prickle emergence in plants.

A New Intra-Specific and High-Resolution Genetic Map of Eggplant Based on a RIL Population, and Location of QTLs Related to Plant Anthocyanin Pigmentation and Seed Vigour

Eggplant is the second most important solanaceous berry-producing crop after tomato. Despite mapping studies based on bi-parental progenies and GWAS approaches having been performed, an eggplant intraspecific high-resolution map is still lacking. We developed a RIL population from the intraspecific cross ‘305E40’, (androgenetic introgressed line carrying the locus Rfo-Sa1 conferring Fusarium resistance) x ‘67/3’ (breeding line whose genome sequence was recently released). One hundred and sixty-three RILs were genotyped by a genotype-by-sequencing (GBS) approach, which allowed us to identify 10,361 polymorphic sites. Overall, 267 Gb of sequencing data were generated and ~773 M Illumina paired end (PE) reads were mapped against the reference sequence. A new linkage map was developed, including 7249 SNPs assigned to the 12 chromosomes and spanning 2169.23 cM, with iaci@liberoan average distance of 0.4 cM between adjacent markers. This was used to elucidate the genetic bases of seven traits related to anthocyanin content in different organs recorded in three locations as well as seed vigor. Overall, from 7 to 17 QTLs (at least one major QTL) were identified for each trait. These results demonstrate that our newly developed map supplies valuable information for QTL fine mapping, candidate gene identification, and the development of molecular markers for marker assisted selection (MAS) of favorable alleles.

Mapping QTLs for 15 morpho-metric traits in Arabidopsis thaliana using Col-0 × Don-0 population

Genome wide quantitative trait loci (QTL) mapping was conducted in Arabidopsis thaliana using F₂ mapping population (Col-0 × Don-0) and SNPs markers. A total of five linkage groups were obtained with number of SNPs varying from 45 to 59 per linkage group. The composite interval mapping detected a total of 36 QTLs for 15 traits and the number of QTLs ranged from one (root length, root dry biomass, cauline leaf width, number of internodes and internode distance) to seven (for bolting days). The range of phenotypic variance explained (PVE) and logarithm of the odds ratio of these 36 QTLs was found be 0.19-38.17% and 3.0-6.26 respectively. Further, the epistatic interaction detected one main effect QTL and four epistatic QTLs. Five major QTLs viz. Qbd.nbri.4.3, Qfd.nbri.4.2, Qrdm.nbri.5.1, Qncl.nbri.2.2, Qtd.nbri.4.1 with PVE > 15.0% might be useful for fine mapping to identify genes associated with respective traits, and also for development of specialized population through marker assisted selection. The identification of additive and dominant effect QTLs and desirable alleles of each of above mentioned traits would also be important for future research.

Construction of a SNP-Based Genetic Map Using SLAF-Seq and QTL Analysis of Morphological Traits in Eggplant

Eggplant (Solanum melongena; 2n = 24) is an economically important fruit crop of the family Solanaceae that was domesticated in India and Southeast Asia. Construction of a high-resolution genetic map and map-based gene mining in eggplant have lagged behind other crops within the family such as tomato and potato. In this study, we conducted high-throughput single nucleotide polymorphism (SNP) discovery in the eggplant genome using specific length amplified fragment (SLAF) sequencing and constructed a high-density genetic map for the quantitative trait locus (QTL) analysis of multiple traits. An interspecific F2 population of 121 individuals was developed from the cross between cultivated eggplant "1836" and the wild relative S. linnaeanum "1809." Genomic DNA extracted from parental lines and the F2 population was subjected to high-throughput SLAF sequencing. A total of 111.74 Gb of data and 487.53 million pair-end reads were generated. A high-resolution genetic map containing 2,122 SNP markers and 12 linkage groups was developed for eggplant, which spanned 1530.75 cM, with an average distance of 0.72 cM between adjacent markers. A total of 19 QTLs were detected for stem height and fruit and leaf morphology traits of eggplant, explaining 4.08-55.23% of the phenotypic variance. These QTLs were distributed on nine linkage groups (LGs), but not on LG2, 4, and 9. The number of SNPs ranged from 2 to 11 within each QTL, and the genetic interval varied from 0.15 to 10.53 cM. Overall, the results establish a foundation for the fine mapping of complex QTLs, candidate gene identification, and marker-assisted selection of favorable alleles in eggplant breeding.

Comprehensive Characterization of Simple Sequence Repeats in Eggplant (Solanum melongena L.) Genome and Construction of a Web Resource

We have characterized the simple sequence repeat (SSR) markers of the eggplant (Solanum melongena) using a recent high quality sequence of its whole genome. We found nearly 133,000 perfect SSRs, a density of 125.5 SSRs/Mbp, and also about 178,400 imperfect SSRs. Of the perfect SSRs, 15.6% were complex, with two stretches of repeats separated by an intervening block of <100 nt. Di- and trinucleotide SSRs accounted, respectively, for 43 and 37% of the total. The SSRs were classified according to their number of repeats and overall length, and were assigned to their linkage group. We found 2,449 of the perfect SSRs in 2,086 genes, with an overall density of 18.5 SSRs/Mbp across the gene space; 3,524 imperfect SSRs were present in 2,924 genes at a density of 26.7 SSRs/Mbp. Putative functions were assigned via ontology to genes containing at least one SSR. Using this data we developed an "Eggplant Microsatellite DataBase" (EgMiDB) which permits identification of SSR markers in terms of their location on the genome, type of repeat (perfect vs. imperfect), motif type, sequence, repeat number and genomic/gene context. It also suggests forward and reverse primers. We employed an in silico PCR analysis to validate these SSR markers, using as templates two CDS sets and three assembled transcriptomes obtained from diverse eggplant accessions.

Detailed mapping of a resistance locus against Fusarium wilt in cultivated eggplant (Solanum melongena)

KEY MESSAGE

This is the first report on genetic mapping of a resistance locus against Fusarium wilt caused by the plant pathogen Fusarium oxysporum f. sp. melongenae in cultivated eggplant.

ABSTRACT

Fusarium wilt, caused by the plant pathogen Fusarium oxysporum f. sp. melongenae, is a major soil-borne disease threatening stable production in eggplant (Solanum melongena). Although three eggplant germplasms, LS1934, LS174, and LS2436, are known to be highly resistant to the pathogen, their resistance loci have not been mapped. In this study, we performed quantitative trait locus analyses in F2:3 populations and detected a resistance locus, FM1, at the end of chromosome 2, with two alleles, Fm1(L) and Fm1(E), in the F2 populations LWF2 [LS1934 × WCGR112-8 (susceptible)] and EWF2 [EPL-1 (derived from LS174) × WCGR112-8], respectively. The percentage of phenotypic variance explained by Fm1(L) derived from LS1934 was 75.0% [Logarithm of the odds (LOD) = 29.3], and that explained by Fm1(E) derived from EPL-1 was 92.2% (LOD = 65.8). Using backcrossed inbred lines, we mapped FM1 between two simple sequence repeat markers located ~4.881 cM apart from each other. Comparing the location of the above locus to those of previously reported ones, the resistance locus Rfo-sa1 from an eggplant ally (Solanum aethiopicum gr. Gilo) was mapped very close to FM1, whereas another resistance locus, from LS2436, was mapped to the middle of chromosome 4. This is the first report of mapping of a Fusarium resistance locus in cultivated eggplant. The availability of resistance-linked markers will enable the application of marker-assisted selection to overcome problems posed by self-incompatibility and introduction of negative traits because of linkage drag, and will lead to clear understanding of genetic mechanism of Fusarium resistance.

Mapping Quantitative Trait Loci Affecting Biochemical and Morphological Fruit Properties in Eggplant (Solanum melongena L.)

Eggplant berries are a source of health-promoting metabolites including antioxidant and nutraceutical compounds, mainly anthocyanins and chlorogenic acid; however, they also contain some anti-nutritional compounds such as steroidal glycoalkaloids (SGA) and saponins, which are responsible for the bitter taste of the flesh and with potential toxic effects on humans. Up to now, Quantitative Trait Loci (QTL) for the metabolic content are far from being characterized in eggplant, thus hampering the application of breeding programs aimed at improving its fruit quality. Here we report on the identification of some QTL for the fruit metabolic content in an F2 intraspecific mapping population of 156 individuals, obtained by crossing the eggplant breeding lines "305E40" × "67/3." The same population was previously employed for the development of a RAD-tag based linkage map and the identification of QTL associated to morphological and physiological traits. The mapping population was biochemically characterized for both fruit basic qualitative data, like dry matter, °Brix, sugars, and organic acids, as well as for health-related compounds such chlorogenic acid, (the main flesh monomeric phenol), the two peel anthocyanins [i.e., delphinidin-3-rutinoside (D3R) and delphinidin-3-(p- coumaroylrutinoside)-5-glucoside (nasunin)] and the two main steroidal glycoalkaloids, solasonine, and solamargine. For most of the traits, one major QTL (PVE ≥10%) was spotted and putative orthologies with other Solanaceae crops are discussed. The present results supply valuable information to eggplant breeders on the inheritance of key fruit quality traits, thus providing potential tools to assist future breeding programs.

Genetic Map Construction and Quantitative Trait Locus (QTL) Detection of Six Economic Traits Using an F2 Population of the Hybrid from Saccharina longissima and Saccharina japonica

Saccharina (Laminaria) is one of the most important economic seaweeds. Previously, four genetic linkage maps of Saccharina have been constructed and five QTLs have been identified. However, they were not enough for its breeding. In this work, Saccharina longissima (♀) and Saccharina japonica (♂), which showed obvious differences in morphology and genetics, were applied in hybridization to yield the F₂ mapping population with 102 individuals. Using these 102 F₂ hybrids, the genetic linkage map of Saccharina was constructed by MapMaker software based on 37 amplified fragment length polymorphisms (AFLPs), 22 sequence-related amplified polymorphisms (SRAPs) and 139 simple sequence repeats (SSRs) markers. Meanwhile, QTL analysis was performed for six economic traits. The linkage map constructed in this research consisted of 422 marker loci (137 AFLPs, 57 SRAPs and 228 SSRs), which formed 45 linkage groups (LGs) with an average marker space of 7.92 cM; they spanned a total length of 2233.1 cM, covering the whole estimated genome size. A total of 29 QTLs were identified for six economic traits, which explained 1.06 to 64.00% of phenotypic variation, including three QTLs for frond length (FL) and raw weight (RW), five QTLs for frond width (FW), two QTLs for frond fascia width (FFW) and frond thickness (FT), and fourteen QTLs for base shape (BS). The results of this research will improve the breeding efficiency and be beneficial for marker-assisted selection (MAS) schemes in Saccharina breeding.

Location of chlorogenic acid biosynthesis pathway and polyphenol oxidase genes in a new interspecific anchored linkage map of eggplant

BACKGROUND

Eggplant is a powerful source of polyphenols which seems to play a key role in the prevention of several human diseases, such as cancer and diabetes. Chlorogenic acid is the polyphenol most present in eggplant, comprising between the 70% and 90% of the total polyphenol content. Introduction of the high chlorogenic acid content of wild relatives, such as S. incanum, into eggplant varieties will be of great interest. A potential side effect of the increased level polyphenols could be a decrease on apparent quality due to browning caused by the polyphenol oxidase enzymes mediated oxidation of polyphenols. We report the development of a new interspecific S. melongena × S. incanum linkage map based on a first backcross generation (BC1) towards the cultivated S. melongena as a tool for introgressing S. incanum alleles involved in the biosynthesis of chlorogenic acid in the genetic background of S. melongena.

RESULTS

The interspecific genetic linkage map of eggplant developed in this work anchor the most informative previously published genetic maps of eggplant using common markers. The 91 BC1 plants of the mapping population were genotyped with 42 COSII, 99 SSRs, 88 AFLPs, 9 CAPS, 4 SNPs and one morphological polymorphic markers. Segregation marker data resulted in a map encompassing 1085 cM distributed in 12 linkage groups. Based on the syntheny with tomato, the candidate genes involved in the core chlorogenic acid synthesis pathway in eggplant (PAL, C4H, 4CL, HCT, C3'H, HQT) as well as five polyphenol oxidase (PPO1, PPO2, PPO3, PPO4, PPO5) were mapped. Except for 4CL and HCT chlorogenic acid genes were not linked. On the contrary, all PPO genes clustered together. Candidate genes important in domestication such as fruit shape (OVATE, SISUN1) and prickliness were also located.

CONCLUSIONS

The achievements in location of candidate genes will allow the search of favorable alleles employing marker-assisted selection in order to develop new varieties with higher chlorogenic content alongside a lower polyphenol oxidase activity. This will result into an enhanced product showing a lower fruit flesh browning with improved human health properties.

Linkage disequilibrium and genome-wide association analysis for anthocyanin pigmentation and fruit color in eggplant

Background

The genome-wide association (GWA) approach represents an alternative to biparental linkage mapping for determining the genetic basis of trait variation. Both approaches rely on recombination to re-arrange the genome, and seek to establish correlations between phenotype and genotype. The major advantages of GWA lie in being able to sample a much wider range of the phenotypic and genotypic variation present, in being able to exploit multiple rounds of historical recombination in many different lineages and to include multiple accessions of direct relevance to crop improvement.

Results

A 191 accessions eggplant (Solanum melongena L.) association panel, comprising a mixture of breeding lines, old varieties and landrace selections originating from Asia and the Mediterranean Basin, was SNP genotyped and scored for anthocyanin pigmentation and fruit color at two locations over two years. The panel formed two major clusters, reflecting geographical provenance and fruit type. The global level of linkage disequilibrium was 3.4 cM. A mixed linear model appeared to be the most appropriate for GWA. A set of 56 SNP locus/phenotype associations was identified and the genomic regions harboring these loci were distributed over nine of the 12 eggplant chromosomes. The associations were compared with the location of known QTL for the same traits.

Conclusion

The GWA mapping approach was effective in validating a number of established QTL and, thanks to the wide diversity captured by the panel, was able to detect a series of novel marker/trait associations.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-896) contains supplementary material, which is available to authorized users.

Linkage group correction using epistatic distorted markers in F2 and backcross populations

Epistasis has been frequently observed in all types of mapping populations. However, relatively little is known about the effect of epistatic distorted markers on linkage group construction. In this study, a new approach was proposed to correct the recombination fraction between epistatic distorted markers in backcross and F2 populations under the framework of fitness and liability models. The information for three or four markers flanking with an epistatic segregation distortion locus was used to estimate the recombination fraction by the maximum likelihood method, implemented via an expectation-maximisation algorithm. A set of Monte Carlo simulation experiments along with a real data analysis in rice was performed to validate the new method. The results showed that the estimates from the new method are unbiased. In addition, five statistical properties for the new method in a backcross were summarised and confirmed by theoretical, simulated and real data analyses.

QTL mapping in eggplant reveals clusters of yield-related loci and orthology with the tomato genome

In spite of its widespread cultivation and nutritional and economic importance, the eggplant (Solanum melongena L.) genome has not been extensively explored. A lack of knowledge of the patterns of inheritance of key agronomic traits has hindered the exploitation of marker technologies to accelerate its genetic improvement. An already established F₂ intraspecific population of eggplant bred from the cross ‘305E40’ x ‘67/3’ was phenotyped for 20 agronomically relevant traits at two sites. Up to seven quantitative trait loci (QTL) per trait were identified and the percentage of the phenotypic variance (PV) explained per QTL ranged from 4 to 93%. Not all the QTL were detectable at both sites, but for each trait at least one major QTL (PV explained ≥10%) was identified. Although no detectable QTL x environment interaction was found, some QTL identified were location-specific. Many of the fruit-related QTL clustered within specific chromosomal regions, reflecting either linkage and/or pleiotropy. Evidence for putative tomato orthologous QTL/genes was obtained for several of the eggplant QTL. Information regarding the inheritance of key agronomic traits was obtained. Some of the QTL, along with their respective linked markers, may be useful in the context of marker-assisted breeding.

The population structure and diversity of eggplant from Asia and the Mediterranean Basin

A collection of 238 eggplant breeding lines, heritage varieties and selections within local landraces provenanced from Asia and the Mediterranean Basin was phenotyped with respect to key plant and fruit traits, and genotyped using 24 microsatellite loci distributed uniformly throughout the genome. STRUCTURE analysis based on the genotypic data identified two major sub-groups, which to a large extent mirrored the provenance of the entries. With the goal to identify true-breeding types, 38 of the entries were discarded on the basis of microsatellite-based residual heterozygosity, along with a further nine which were not phenotypically uniform. The remaining 191 entries were scored for a set of 19 fruit and plant traits in a replicated experimental field trial. The phenotypic data were subjected to principal component and hierarchical principal component analyses, allowing three major morphological groups to be identified. All three morphological groups were represented in both the “Occidental” and the “Oriental” germplasm, so the correlation between the phenotypic and the genotypic data sets was quite weak. The relevance of these results for evolutionary studies and the further improvement of eggplant are discussed. The population structure of the core set of germplasm shows that it can be used as a basis for an association mapping approach.

The maternal cytoplasmic environment may be involved in the viability selection of gametes and zygotes

Segregation distortion is the phenomenon whereby the observed genotypic frequencies of a locus fall outside the expected Mendelian segregation ratio, and it is increasingly recognised as a potentially powerful evolutionary force. The main reason for segregation distortion is a difference in the viability of gametes and zygotes caused by viability loci in the segregating progeny. However, the maternal cytoplasm may also be involved in the viability selection of gametes and zygotes. The objectives of this study were to map the segregation distortion loci (SDL) in maize and to test the hypothesis that the viability of gametes and zygotes may also be associated with the maternal cytoplasmic environment. In the present study, a reciprocal mating design was conducted to generate an F2-segregating population. A linkage map was constructed with 126 microsatellite markers. A whole-genome scan was performed to detect the SDL in segregating populations with different maternal cytoplasm environments. Altogether, 14 SDL with strong LOD (logarithm (base 10) of odds) supports were identified in the specifically designed F2 populations. Interestingly, we found dramatic changes in the genotypic frequencies of the SDL in the two maternal cytoplasmic backgrounds, which indicated a change in the viability of gametes and zygotes in different cytoplasmic environments. Furthermore, in the JB cytoplasmic background, most of the detected SDL and complete distortion markers exhibited similar bias patterns favouring the Y53 alleles. These results suggested that selfish cytoplasmic elements may have an important role in shaping the patterns of segregation distortion in F2 populations through selective viability of gametes and zygotes.

SNP-based analysis of genetic diversity in anther-derived rice by whole genome sequencing

BACKGROUND

Anther culture has advantage to obtain a homozygous progeny by induced doubling of haploid chromosomes and to improve selection efficiency for invaluable agronomical traits. Therefore, anther culturing is widely utilized to breed new varieties and to induce genetic variations in several crops including rice. Genome sequencing technologies allow the detection of a massive number of DNA polymorphism such as SNPs and Indels between closely related cultivars. These DNA polymorphisms permit the rapid identification of genetic diversity among cultivars and genomic locations of heritable traits. To estimate sequence diversity derived from anther culturing, we performed whole-genome resequencing of five Korean rice accessions, including three anther culture lines (BLB, HY-04 and HY-08), their progenitor cultivar (Hwayeong), and an additional japonica cultivar (Dongjin).

RESULTS

A total of 1,165 × 106 raw reads were generated with over 58× coverage that detected 1,154,063 DNA polymorphisms between the Korean rice accessions and Nipponbare. We observed that in Hwayeong and its progenies, 0.64 SNP was found per one kb of Nipponbare genome, while Dongjin, bred by a conventional breeding method, had a lower number of SNPs (0.45 SNP/kb). Among 1,154,063 DNA polymorphisms, 29,269 non-synonymous SNPs located on 30,013 genes and these genes were functionally classified based on gene ontology (GO). We also analyzed line-specific SNPs which were estimated 1 ~ 3% of the total SNPs. The frequency of non-synonymous SNPs in each accession ranged from 26 SNPs in Hwayeong to 214 SNPs in HY-04.

CONCLUSIONS

The genetic difference we detected between the progenies derived from anther culture and their mother cultivar is due to somaclonal variation during tissue culture process, such as karyotype change, chromosome rearrangement, gene amplification and deletion, transposable element, and DNA methylation. Detection of genome-wide DNA polymorphisms by high-throughput sequencer enabled to identify sequence diversity derived from anther culturing and genomic locations of heritable traits. Furthermore, it will provide an invaluable resource to identify molecular markers and genes associated with diverse traits of agronomical importance.

A RAD tag derived marker based eggplant linkage map and the location of QTLs determining anthocyanin pigmentation

Both inter- and intra-specific maps have been developed in eggplant (Solanum melongena L.). The former benefit from an enhanced frequency of marker polymorphism, but their relevance to marker-assisted crop breeding is limited. Combining the restriction-site associated DNA strategy with high throughput sequencing has facilitated the discovery of a large number of functional single nucleotide polymorphism (SNP) markers discriminating between the two eggplant mapping population parental lines '305E40' and '67/3'. A set of 347 de novo SNPs, together with 84 anchoring markers, were applied to the F(2) mapping population bred from the cross '305E40' x '67/3' to construct a linkage map. In all, 415 of the 431 markers were assembled into twelve major and one minor linkage group, spanning 1,390 cM, and the inclusion of established markers allowed each linkage group to be assigned to one of the 12 eggplant chromosomes. The map was then used to discover the genetic basis of seven traits associated with anthocyanin content. Each of the traits proved to be controlled by between one and six quantitative trait loci (QTL), of which at least one was a major QTL. Exploitation of syntenic relationships between the eggplant and tomato genomes facilitated the identification of potential candidate genes for the eggplant QTLs related to anthocyanin accumulation. The intra-specific linkage map should have utility for elucidating the genetic basis of other phenotypic traits in eggplant.

Development of gene-based markers and construction of an integrated linkage map in eggplant by using Solanum orthologous (SOL) gene sets

We constructed an integrated DNA marker linkage map of eggplant (Solanum melongena L.) using DNA marker segregation data sets obtained from two independent intraspecific F(2) populations. The linkage map consisted of 12 linkage groups and encompassed 1,285.5 cM in total. We mapped 952 DNA markers, including 313 genomic SSR markers developed by random sequencing of simple sequence repeat (SSR)-enriched genomic libraries, and 623 single-nucleotide polymorphisms (SNP) and insertion/deletion polymorphisms (InDels) found in eggplant-expressed sequence tags (ESTs) and related genomic sequences [introns and untranslated regions (UTRs)]. Because of their co-dominant inheritance and their highly polymorphic and multi-allelic nature, the SSR markers may be more versatile than the SNP and InDel markers for map-based genetic analysis of any traits of interest using segregating populations derived from any intraspecific crosses of practical breeding materials. However, we found that the distribution of microsatellites in the genome was biased to some extent, and therefore a considerable part of the eggplant genome was first detected when gene-derived SNP and InDel markers were mapped. Of the 623 SNP and InDel markers mapped onto the eggplant integrated map, 469 were derived from eggplant unigenes contained within Solanum orthologous (SOL) gene sets (i.e., sets of orthologous unigenes from eggplant, tomato, and potato). Out of the 469 markers, 326 could also be mapped onto the tomato map. These common markers will be informative landmarks for the transfer of tomato's more saturated genomic information to eggplant and will also provide comparative information on the genome organization of the two solanaceous species. The data are available from the DNA marker database of vegetables, VegMarks (http://vegmarks.nivot.affrc.go.jp).

Development of selective markers linked to a major QTL for parthenocarpy in eggplant (Solanum melongena L.)

Parthenocarpy, the ability to set fruits without pollination, is a useful trait for setting fruit under unfavorable conditions. To identify the loci controlling parthenocarpy in eggplant (Solanum melongena L.), we constructed linkage maps by using co-dominant simple sequence repeat and single nucleotide polymorphism markers in F(2) populations derived from intraspecific crosses between two non-parthenocarpic lines (LS1934 and Nakate-Shinkuro) and a parthenocarpic line (AE-P03). Total map distances were 1,414.6 cM (ALF2: LS1934 x AE-P03) and 1,153.8 cM (NAF2: Nakate-Shinkuro x AE-P03), respectively. Quantitative trait locus (QTL) analyses revealed two QTLs on chromosomes 3 and 8, which we denoted as Controlling parthenocarpy3.1 (Cop3.1) and Cop8.1, respectively. The percentage of phenotypic variance explained (PVE) of Cop3.1 was 6.3% in ALF2 (LOD = 4.2) and 10.6% in NAF2 (LOD = 3.0). The PVE of Cop8.1 was 45.7% in ALF2 (LOD = 23.8) and 29.7% in NAF2 (LOD = 7.9). Using a population of backcross inbred lines, we confirmed the effect of Cop8.1, but there was no evidence to support the contribution of Cop3.1. We need to verify the effect of Cop3.1 under various temperature conditions. In addition, we clarified the effectiveness of selective SSR markers, emf21H22 and emh11J10, mapped on each side of Cop8.1 in other F(2) populations derived from various parental combinations. This is the first report concerning QTL analysis of parthenocarpy in eggplant using molecular markers. It will be useful in marker-assisted selection and in revealing the genomic mechanism underlying parthenocarpy in eggplant.

Identification of SNP and SSR markers in eggplant using RAD tag sequencing

Background

The eggplant (Solanum melongena L.) genome is relatively unexplored, especially compared to those of the other major Solanaceae crops tomato and potato. In particular, no SNP markers are publicly available; on the other hand, over 1,000 SSR markers were developed and publicly available. We have combined the recently developed Restriction-site Associated DNA (RAD) approach with Illumina DNA sequencing for rapid and mass discovery of both SNP and SSR markers for eggplant.

Results

RAD tags were generated from the genomic DNA of a pair of eggplant mapping parents, and sequenced to produce ~17.5 Mb of sequences arrangeable into ~78,000 contigs. The resulting non-redundant genomic sequence dataset consisted of ~45,000 sequences, of which ~29% were putative coding sequences and ~70% were in common between the mapping parents. The shared sequences allowed the discovery of ~10,000 SNPs and nearly 1,000 indels, equivalent to a SNP frequency of 0.8 per Kb and an indel frequency of 0.07 per Kb. Over 2,000 of the SNPs are likely to be mappable via the Illumina GoldenGate assay. A subset of 384 SNPs was used to successfully fingerprint a panel of eggplant germplasm, producing a set of informative diversity data. The RAD sequences also included nearly 2,000 putative SSRs, and primer pairs were designed to amplify 1,155 loci.

Conclusion

The high throughput sequencing of the RAD tags allowed the discovery of a large number of DNA markers, which will prove useful for extending our current knowledge of the genome organization of eggplant, for assisting in marker-aided selection and for carrying out comparative genomic analyses within the Solanaceae family.