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

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.

Genome-Wide Survey and Development of the First Microsatellite Markers Database (AnCorDB) in Anemone coronaria L

Anemone coronaria L. (2n = 2x = 16) is a perennial, allogamous, highly heterozygous plant marketed as a cut flower or in gardens. Due to its large genome size, limited efforts have been made in order to develop species-specific molecular markers. We obtained the first draft genome of the species by Illumina sequencing an androgenetic haploid plant of the commercial line “MISTRAL^® Magenta”. The genome assembly was obtained by applying the MEGAHIT pipeline and consisted of 2 × 10⁶ scaffolds. The SciRoKo SSR (Simple Sequence Repeats)-search module identified 401.822 perfect and 188.987 imperfect microsatellites motifs. Following, we developed a user-friendly “Anemone coronaria Microsatellite DataBase” (AnCorDB), which incorporates the Primer3 script, making it possible to design couples of primers for downstream application of the identified SSR markers. Eight genotypes belonging to eight cultivars were used to validate 62 SSRs and a subset of markers was applied for fingerprinting each cultivar, as well as to assess their intra-cultivar variability. The newly developed microsatellite markers will find application in Breeding Rights disputes, developing genetic maps, marker assisted breeding (MAS) strategies, as well as phylogenetic studies.

Genetic Diversity Revealed by Microsatellites in Genus Carya

The genus Carya consists of 17 species divided into 3 sections: Carya or the true hickories, Apocarya or the pecan hickories, and Sinocarya or the Asian hickories. Interspecific hybrids exist and have been used in pecan cultivar development. Nuclear and plastid microsatellite or SSR markers have been useful in distinguishing species, sections, and populations. They provide evidence for hybridity between species and can confirm heredity within crosses. As more sophisticated methods of genomic evaluation are cooperatively developed for use in pecan breeding and selection, the use of these methods will be supplemented and informed by the lessons provided by microsatellite markers, as interpreted across broad germplasm collections. In this study, over 400 Carya accessions from diverse diploid and tetraploid taxa and their interspecific hybrids, maintained at the USDA National Collection of Genetic Resources for Carya (NCGR-Carya), were analyzed using 14 nuclear and 3 plastid microsatellite markers. Principal coordinate analysis showed clear taxonomic classifications at multiple taxonomic levels along with patterns of interspecific hybridity. Evidence was also found for genetic differences associated with geographic distribution. The results indicate that this group of markers is useful in examining and characterizing populations and hybrids in the genus Carya and may help delineate the composition of a core collection to help characterize the NCGR-Carya repository collection for use in its pecan breeding program. The SSR fingerprints of the inventories of the USDA NCGR-Carya repository can also be used as a reference for identifying unknown pecan trees for growers.

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.

Identification of a large dataset of SNPs in hair-fin anchovy (Setipinna tenuifilis) based on RAD-seq

Hair-fin anchovy (Setipinna tenuifilis) is an economically important fish distributed in the West Indian Ocean and the Northwest Pacific Ocean. In this study, 154 individuals in eight populations of S. tenuifilis were sequenced and 850 million raw reads were obtained using restriction site-associated DNA sequencing (RAD-seq). First, we identified 14 012 044 hypothetical SNP markers. A dataset of 199 903 high-quality SNPs was collected after further screening. These SNPs have a strong ability to test the genetic diversity between the eight populations. The differentiation and genetic law between samples were explored based on SNPs in populations of S. tenuifilis. The results of this study will provide data for protecting the genetic resources of the species.

High-Density Genetic Linkage Map Construction and White Rot Resistance Quantitative Trait Loci Mapping for Genus Vitis Based on Restriction Site-Associated DNA Sequencing

Grape white rot (Coniothyrium diplodiella) is a major fungal disease affecting grape yield and quality. Quantitative trait locus (QTL) analysis is an important method for studying important horticultural traits of grapevine. This study was conducted to construct a high-density map and conduct QTL mapping for grapevine white rot resistance. A mapping population with 177 genotypes was developed from interspecific hybridization of a white rot-resistant cultivar (Vitis vinifera × V. labrusca 'Zhuosexiang') and white rot-susceptible cultivar (V. vinifera 'Victoria'). Single-nucleotide polymorphism (SNP) markers were developed by restriction site-associated DNA sequencing. The female, male, and integrated maps contained 2,501, 4,110, and 6,249 SNP markers with average genetic distances of adjacent markers of 1.25, 0.77, and 0.50 cM, respectively. QTL mapping was conducted based on white rot resistance identification of 177 individuals in July and August of 2017 and 2018. Notably, one stable QTL related to white rot resistance was detected and located on linkage group LG14. The phenotypic variance ranged from 12.93 to 13.43%. An SNP marker (chr14_3929380), which cosegregated with white rot resistance, was discovered and shows potential for use in marker-assisted selection to generate new grapevine cultivars with resistance to white rot.

Resolving phylogeny and polyploid parentage using genus-wide genome-wide sequence data from birch trees

Numerous plant genera have a history including frequent hybridisation and polyploidisation (allopolyploidisation), which means that their phylogeny is a network of reticulate evolution that cannot be accurately depicted as a bifurcating tree with a single tip per species. The genus Betula, which contains many ecologically important tree species, is a case in point. We generated genome-wide sequence reads for 27 diploid and 36 polyploid Betula species or subspecies using restriction site associated DNA (RAD) sequences. These reads were assembled into contigs with a mean length of 675 bp. We reconstructed the evolutionary relationships among diploid Betula species using both supermatrix (concatenation) and species tree methods. We identified the closest diploid relatives of the polyploids according to the relative rates at which reads from polyploids mapped to contigs from different diploid species within a concatenated reference sequence. By mapping reads from allopolyploids to their different putative diploid relatives we assembled contigs from the putative sub-genomes of allopolyploid taxa. We used these to build new phylogenies that included allopolyploid sub-genomes as separate tips. This approach yielded a highly evidenced phylogenetic hypothesis for the genus Betula, including the complex reticulate origins of the majority of its polyploid taxa. Our phylogeny divides the genus into two well supported clades, which, interestingly, differ in their seed-wing morphology. We therefore propose to split Betula into two subgenera.

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.

High-density genetic linkage map construction and cane cold hardiness QTL mapping for Vitis based on restriction site-associated DNA sequencing

Background

Cold hardiness is an important agronomic trait and can significantly affect grape production and quality. Until now, there are no reports focusing on cold hardiness quantitative trait loci (QTL) mapping. In this study, grapevine interspecific hybridisation was carried out with the maternal parent ‘Cabernet sauvignon’ and paternal parent ‘Zuoyouhong’. A total of 181 hybrid offspring and their parents were used as samples for restriction-site associated DNA sequencing (RAD). Grapevine cane phloem and xylem cold hardiness of the experimental material was detected using the low-temperature exotherm method in 2016, 2017 and 2018. QTL mapping was then conducted based on the integrated map.

Results

We constructed a high-density genetic linkage map with 16,076, 11,643, and 25,917 single-nucleotide polymorphism (SNP) markers anchored in the maternal, paternal, and integrated maps, respectively. The average genetic distances of adjacent markers in the maps were 0.65 cM, 0.77 cM, and 0.41 cM, respectively. Colinearity analysis was conducted by comparison with the grape reference genome and showed good performance. Six QTLs were identified based on the phenotypic data of 3 years and they were mapped on linkage group (LG) 2, LG3, and LG15. Based on QTL results, candidate genes which may be involved in grapevine cold hardiness were selected.

Conclusions

High-density linkage maps can facilitate grapevine fine QTL mapping, genome comparison, and sequence assembly. The cold hardiness QTL mapping and candidate gene discovery performed in this study provide an important reference for molecular-assisted selection in grapevine cold hardiness breeding.

Transcriptional dynamics of Zn-accumulation in developing kernels of maize reveals important Zn-uptake mechanisms

In the present study, transcriptomic analysis of 10-days old baby kernels of two contrasting maize genotypes, namely VQL-2 (high kernel Zn accumulator) and CM-145 (low kernel Zn accumulator), under low- and optimum- soil Zn conditions generated 1948 differentially expressed transcripts. Among these, 666 and 437 transcripts were up-regulated and down-regulated respectively in VQL-2; whereas, 437 and 408 transcripts were up-regulated and down-regulated respectively in CM-145. Remarkably, 135 transcription factors and 77 known Zn transporters expressed differentially. By comparing the transcripts differentially expressed between the optimum-Zn and low-Zn libraries of the contrasting genotypes, we identified 21,986 and 26,871 SNPs, respectively. Similarly, 6810 and 8192 InDels were found between optimum- and low-Zn growing conditions, respectively. Further, 21 differentially expressed genes were co-localized with already known QTLs associated with Zn uptake, such as qZn10, CQZnK9-1 and YNZnK6. These findings will be useful to develop high Zn-accumulator maize through marker-assisted breeding in future.

Role of NGS and SNP genotyping methods in sugarcane improvement programs

Sugarcane (Saccharum spp.) is one of the most economically significant crops because of its high sucrose content and it is a promising biomass feedstock for biofuel production. Sugarcane genome sequencing and analysis is a difficult task due to its heterozygosity and polyploidy. Long sequence read technologies, PacBio Single-Molecule Real-Time (SMRT) sequencing, the Illumina TruSeq, and the Oxford Nanopore sequencing could solve the problem of genome assembly. On the applications side, next generation sequencing (NGS) technologies played a major role in the discovery of single nucleotide polymorphism (SNP) and the development of low to high throughput genotyping platforms. The two mainstream high throughput genotyping platforms are the SNP microarray and genotyping by sequencing (GBS). This paper reviews the NGS in sugarcane genomics, genotyping methodologies, and the choice of these methods. Array-based SNP genotyping is robust, provides consistent SNPs, and relatively easier downstream data analysis. The GBS method identifies large scale SNPs across the germplasm. A combination of targeted GBS and array-based genotyping methods should be used to increase the accuracy of genomic selection and marker-assisted breeding.

The taxonomic significance of ddRADseq based microsatellite markers in the closely related species of Heracleum (Apiaceae)

Many studies on Heracleum have shown poor correspondence between observed molecular clusters and established taxonomic classification amongst closely related species. This might reflect both unresolved taxonomy but perhaps also a lack of good genetic markers. This lack of appropriate and cost effective species-specific genetic markers hinders a resolved relationship for the species complex, and this in turn causes profound management challenges for a genus that contains both endemic species, with important ecological roles, and species with an invasive potential. Microsatellites are traditionally considered markers of choice for comprehensive, yet inexpensive, analyses of genetic variation, including examination of population structure, species identity, linkage map construction and cryptic speciation. In this study, we have used double digest restriction site associated DNA sequencing (ddRADseq) to develop microsatellite markers in Heracleum rechingeri. Genomic DNA from three individuals were digested with Sbf1 and Nde1 and size selected for library construction. The size-selected fragments were sequenced on an Ion Torrent sequencer and a total of 54 microsatellite sequences were bioinformatically confirmed. Twenty five loci were then tested for amplification, resulting in 19 of these being successfully amplified across eight species, comprising both the so-called thick-stemmed species (H. persicum, H. rechingeri, H. gorganicum and H. lasiopetalum), and thin-stemmed species (H. anisactis, H. pastinasifolium and H. transcaucasicum). Both Bayesian and distance-based clustering, and principal coordinate analyses clearly separated these into two groups. Surprisingly, three H. pastinacifolium populations were not separated from populations of the morphologically similar endemic species, H. anisactis, suggesting lack of genetic differentiation. Likewise, high genetic similarity was found between H. persicum and H. rechingeri populations, questioning taxonomic separation at the species level between these taxa. Further analyses are needed to re-evaluate the taxonomic significance of observed morphological variability currently applied to distinguish these sister taxa. Nevertheless, our results represent progress in the effort to develop cost-efficient molecular tools for species discrimination in this genus.

Genetic variation in eggplant for Nitrogen Use Efficiency under contrasting NO3 - supply

Eggplant (Solanum melongena L.) yield is highly sensitive to N fertilization, the excessive use of which is responsible for environmental and human health damage. Lowering N input together with the selection of improved Nitrogen-Use-Efficiency (NUE) genotypes, more able to uptake, utilize, and remobilize N available in soils, can be challenging to maintain high crop yields in a sustainable agriculture. The aim of this study was to explore the natural variation among eggplant accessions from different origins, in response to Low (LN) and High (HN) Nitrate (NO₃ ^- ) supply, to identify NUE-contrasting genotypes and their NUE-related traits, in hydroponic and greenhouse pot experiments. Two eggplants, AM222 and AM22, were identified as N-use efficient and inefficient, respectively, in hydroponic, and these results were confirmed in a pot experiment, when crop yield was also evaluated. Overall, our results indicated the key role of N-utilization component (NUtE) to confer high NUE. The remobilization of N from leaves to fruits may be a strategy to enhance NUtE, suggesting glutamate synthase as a key enzyme. Further, omics technologies will be used for focusing on C-N metabolism interacting networks. The availability of RILs from two other selected NUE-contrasting genotypes will allow us to detect major genes/quantitative trait loci related to NUE.

A First Insight into a Draft Genome of Silver Sillago (Sillago sihama) via Genome Survey Sequencing

Sillago sihama has high economic value and is one of the most attractive aquaculture species in China. Despite its economic importance, studies of its genome have barely been performed. In this study, we conducted a first genomic survey of S. sihama using next-generation sequencing (NGS). In total, 45.063 Gb of high-quality sequence data were obtained. For the 17-mer frequency distribution, the genome size was estimated to be 508.50 Mb. The sequence repeat ratio was calculated to be 21.25%, and the heterozygosity ratio was 0.92%. Reads were assembled into 1,009,363 contigs, with a N50 length of 1362 bp, and then into 814,219 scaffolds, with a N50 length of 2173 bp. The average Guanine and Cytosine (GC) content was 45.04%. Dinucleotide repeats (56.55%) were the dominant form of simple sequence repeats (SSR).

Optimizing ddRADseq in Non-Model Species: A Case Study in Eucalyptus dunnii Maiden

Restriction site-associated DNA sequencing (RADseq) and its derived protocols, such as double digest RADseq (ddRADseq), offer a flexible and highly cost-effective strategy for efficient plant genome sampling. This has become one of the most popular genotyping approaches for breeding, conservation, and evolution studies in model and non-model plant species. However, universal protocols do not always adapt well to non-model species. Herein, this study reports the development of an optimized and detailed ddRADseq protocol in Eucalyptus dunnii, a non-model species, which combines different aspects of published methodologies. The initial protocol was established using only two samples by selecting the best combination of enzymes and through optimal size selection and simplifying lab procedures. Both single nucleotide polymorphisms (SNPs) and simple sequence repeats (SSRs) were determined with high accuracy after applying stringent bioinformatics settings and quality filters, with and without a reference genome. To scale it up to 24 samples, we added barcoded adapters. We also applied automatic size selection, and therefore obtained an optimal number of loci, the expected SNP locus density, and genome-wide distribution. Reliability and cross-sequencing platform compatibility were verified through dissimilarity coefficients of 0.05 between replicates. To our knowledge, this optimized ddRADseq protocol will allow users to go from the DNA sample to genotyping data in a highly accessible and reproducible way.

Genetic Analysis for Fruit Phenolics Content, Flesh Color, and Browning Related Traits in Eggplant (Solanum melongena)

Eggplant varieties rich in bioactive chlorogenic acid along with less browning are preferred by consumers. Therefore, genetics of fruit phenolics, fruit flesh colour, and browning related traits were studied in the genotypes of eggplant, comprising of nine cultivated varieties and one accession of eggplant‘s primary genepool wild relative Solanum insanum (INS2). These accessions were genotyped based on the 7335 polymorphic single-nucleotide polymorphisms (SNP) markers. After that, genotypes were crossed in half diallel fashion to produce 45 hybrids. The INS2 displayed the highest values for the total phenolics and chlorogenic acid content (CGA). For all of the biochemical traits studied, significant values of general and specific combining ability (GCA and SCA) effects were determined. The baker ratio estimates were high (>0.75) for all of the traits. Highly significant and positive heterosis (%) was determined for the dry matter, total phenolics, CGA, and area (%) of CGA content. The phenolics content of the fruit (total phenolics and CGA) was not significantly correlated with flesh colour and browning related traits. However, when the path coefficient analysis was performed considering the CGA as a dependent variable, it was determined that the flesh colour related traits most considerably affected the CGA. The genetic distance showed a diminutive correlation with the hybrid means, heterosis, and SCA values. Overall, this study provides important information regarding the underlying genetics of important biochemical traits of eggplant fruit.

Whole-Genome Resequencing of Seven Eggplant (Solanum melongena) and One Wild Relative (S. incanum) Accessions Provides New Insights and Breeding Tools for Eggplant Enhancement

Whole-genome resequencing provides information of great relevance for crop genetics, evolution, and breeding. Here, we present the first whole-genome resequencing study using seven eggplant (Solanum melongena) and one wild relative (Solanum incanum) accessions. These eight accessions were selected for displaying a high phenotypic and genetic diversity and for being the founder parents of an eggplant multiparent advanced generation intercrosses population. By resequencing at an average depth of 19.8× and comparing to the high-quality reference genome "67/3" over 10 million high-reliable polymorphisms were discovered, of which over 9 million (84.5%) were single nucleotide polymorphisms and more than 700,000 (6.5%) InDels. However, while for the S. melongena accessions, the variants identified ranged from 0.8 to 1.3 million, over 9 million were detected for the wild S. incanum. This confirms the narrow genetic diversity of the domesticated eggplant and suggests that introgression breeding using wild relatives can efficiently contribute to broadening the genetic basis of this crop. Differences were observed among accessions for the enrichment in genes regulating important biological processes. By analyzing the distribution of the variants, we identified the potential footprints of old introgressions from wild relatives that can help to unravel the unclear domestication and breeding history. The comprehensive annotation of these eight genomes and the information provided in this study represents a landmark in eggplant genomics and allows the development of tools for eggplant genetics and breeding.

Analysis of evolution and genetic diversity of sweetpotato and its related different polyploidy wild species I. trifida using RAD-seq

BACKGROUND

Sweetpotato (Ipomoea batatas (L.) Lam.) is one of the most important crops from the family of Convolvulaceae. It is widely reported that cultivated sweetpotato was originated from Ipomoea trifida. However, diploid, tetraploid and hexaploid I. trifida were found in nature. The relationship, between them, and among them and sweetpotato, is remaining unclear.

RESULTS

In the present study, we detected the genome diversity and relationship of sweetpotato and different polyploidy types I. trifida using Restriction-site Associated DNA Sequencing (RAD-seq). A total of 38,605 RAD-tags containing 832,204 SNPs had been identified. These tags were annotated using five public databases, about 11,519 tags were aligned to functional genes in various pathways. Based on SNP genotype, phylogenetic relation analysis results confirmed that cultivated sweetpotato has a closer relationship with I. trifida 6× than with I. trifida 4X and I. trifida 2×. Besides, 5042 SSRs were detected in I. trifida 6×, and 3202 pairs of high-quality SSR primers were developed. A total of 68 primers were randomly selected and synthesized, of which 61 were successfully amplified.

CONCLUSION

These results provided new evidence that cultivated sweetpotato originated from I. trifida 6×, and that I. trifida 6× evolved from I. trifida 4X and I. trifida 2×. Therefore, using I. trifida 6× as the model plant of sweetpotato research should be more practical than using I. trifida 2× in the future. Meanwhile, sequence information and markers from the present study will be helpful for sweetpotato and I. trifida studies in the future.

Mapping QTLs for water-use efficiency reveals the potential candidate genes involved in regulating the trait in apple under drought stress

BACKGROUND

Improvement of water-use efficiency (WUE) can effectively reduce production losses caused by drought stress. A better understanding of the genetic determination of WUE in crops under drought stress has great potential value for developing cultivars adapted to arid regions. To identify the genetic loci associated with WUE and reveal genes responsible for the trait in apple, we aim to map the quantitative trait loci (QTLs) for carbon isotope composition, the proxy for WUE, applying two contrasting irrigating regimes over the two-year experiment and search for the candidate genes encompassed in the mapped QTLs.

RESULTS

We constructed a high-density genetic linkage map with 10,172 markers of apple, using single nucleotide polymorphism (SNP) markers obtained through restriction site-associated DNA sequencing (RADseq) and a final segregating population of 350 seedlings from the cross of Honeycrisp and Qinguan. In total, 33 QTLs were identified for carbon isotope composition in apple under both well-watered and drought-stressed conditions. Three QTLs were stable over 2 years under drought stress on linkage groups LG8, LG15 and LG16, as validated by Kompetitive Allele-Specific PCR (KASP) assays. In those validated QTLs, 258 genes were screened according to their Gene Ontology functional annotations. Among them, 28 genes were identified, which exhibited significant responses to drought stress in 'Honeycrisp' and/or 'Qinguan'. These genes are involved in signaling, photosynthesis, response to stresses, carbohydrate metabolism, protein metabolism and modification, hormone metabolism and transport, transport, respiration, transcriptional regulation, and development regulation. They, especially those for photoprotection and relevant signal transduction, are potential candidate genes connected with WUE regulation in drought-stressed apple.

CONCLUSIONS

We detected three stable QTLs for carbon isotope composition in apple under drought stress over 2 years, and validated them by KASP assay. Twenty-eight candidate genes encompassed in these QTLs were identified. These stable genetic loci and series of genes provided here serve as a foundation for further studies on marker-assisted selection of high WUE and regulatory mechanism of WUE in apple exposed to drought conditions, respectively.

A nuclear DNA barcode for eastern North American oaks and application to a study of hybridization in an Arboretum setting

DNA barcoding has proved difficult in a number of woody plant genera, including the ecologically important oak genus Quercus. In this study, we utilized restrictionsite-associated DNA sequencing (RAD-seq) to develop an economical single nucleotide polymorphism (SNP) DNA barcoding system that suffices to distinguish eight common, sympatric eastern North American white oak species. Two de novo clustering pipelines, PyRAD and Stacks, were used in combination with postclustering bioinformatic tools to generate a list of 291 potential SNPs, 80 of which were included in a barcoding toolkit that is easily implemented using MassARRAY mass spectrometry technology. As a proof-of-concept, we used the genotyping toolkit to infer potential hybridization between North American white oaks transplanted outside of their native range (Q. michauxii, Q. montana, Q muehlenbergii/Q. prinoides, and Q. stellata) into a horticultural collection surrounded by natural forests of locally native trees (Q. alba and Q. macrocarpa) in the living collection at The Morton Arboretum (Lisle, IL, USA). Phylogenetic and clustering analyses suggested low rates of hybridization between cultivated and native species, with the exception of one Q. michauxii mother tree, the acorns of which exhibited high admixture from either Q. alba or Q. stellata and Q. macrocarpa, and a hybrid between Q. stellata that appears to have backcrossed almost exclusively to Q. alba. Together, RAD-seq and MassARRAY technologies allow for efficient development and implementation of a multispecies barcode for one of the more challenging forest tree genera.

Genome re-sequencing, SNP analysis, and genetic mapping of the parental lines of a commercial F1 hybrid cultivar of Chinese cabbage

The genome-wide characterization of single nucleotide polymorphism (SNP) between cultivars or between inbred lines contributes to the creation of genetic markers that are important for plant breeding. Functional markers derived from polymorphisms within genes that affect phenotypic variation are especially valuable in plant breeding. Here, we report on the genome re-sequencing and analysis of the two parental inbred lines of the commercial F₁ hybrid Chinese cabbage cultivar "W77". Through the genome-wide identification and classification of the SNPs and indels present in each parental line, we identified about 1,500 putative non-functional genes in each parent. We designed cleaved amplified polymorphic sequence (CAPS) markers using specific mutations found at Eco RI restriction sites in the parental lines and confirmed their Mendelian segregation by constructing a linkage map using 96 F₂ plants derived from the F₁ hybrid cultivar, "W77". Our results and data will be a useful genomic resource for future studies of gene function and metagenomic studies in Chinese cabbage.

Genetic Evaluation of Natural Populations of the Endangered Conifer Thuja koraiensis Using Microsatellite Markers by Restriction-Associated DNA Sequencing

Thuja koraiensis Nakai is an endangered conifer of high economic and ecological value in Jilin Province, China. However, studies on its population structure and conservation genetics have been limited by the lack of genomic data. Here, 37,761 microsatellites (simple sequence repeat, SSR) were detected based on 875,792 de novo-assembled contigs using a restriction-associated DNA (RAD) approach. Among these SSRs, 300 were randomly selected to test for polymorphisms and 96 obtained loci were able to amplify a fragment of expected size. Twelve polymorphic SSR markers were developed to analyze the genetic diversity and population structure of three natural populations. High genetic diversity (mean NA = 5.481, HE = 0.548) and moderate population differentiation (pairwise Fst = 0.048–0.078, Nm = 2.940–4.958) were found in this species. Molecular variance analysis suggested that most of the variation (83%) existed within populations. Combining the results of STRUCTURE, principal coordinate, and neighbor-joining analysis, the 232 individuals were divided into three genetic clusters that generally correlated with their geographical distributions. Finally, appropriate conservation strategies were proposed to protect this species. This study provides genetic information for the natural resource conservation and utilization of T. koraiensis and will facilitate further studies of the evolution and phylogeography of the species.

Resolving relationships and phylogeographic history of the Nyssa sylvatica complex using data from RAD-seq and species distribution modeling

Nyssa sylvatica complex consists of several woody taxa occurring in eastern North America. These taxa were recognized as two or three species including three or four varieties by different authors. Due to high morphological similarities and complexity of morphological variation, classification and delineation of taxa in the group have been difficult and controversial. Here we employ data from RAD-seq to elucidate the genetic structure and phylogenetic relationships within the group. Using the genetic evidence, we evaluate previous classifications and delineate species. We also employ Species Distribution Modeling (SDM) to evaluate impacts of climatic changes on the ranges of the taxa and to gain insights into the relevant refugia in eastern North America. Results from Molecular Variance Analysis (AMOVA), STRUCTURE, phylogenetic analyses using Maximum likelihood, Bayesian Inference, and Splittree methods of RAD-seq data strongly support a two-clade pattern, largely separating samples of N. sylvatica from those of N. biflora-N. ursina mix. Divergence time analysis with BEAST suggests the two clades diverged in the mid Miocene. The ancestor of the present trees of N. sylvatica was suggested to be in the Pliocene and that of N. biflora-N. ursina mix in the end of the Miocene. Results from SDM predicted a smaller range in the southern part of the species present range of each clade during the Last Glacial Maximum (LGM). A northward expansion of the ranges during interglacial period and a northward shift of the ranges in the future under a model of global warming were also predicted. Our results support the recognition of two species in the complex, N. sylvatica and N. biflora, following the phylogenetic species concept. We found no genetic evidence supporting recognitions of intraspecific taxa. However, we propose subsp. ursina and subsp. biflora within N. biflora due to their distinction in habits, distributions, and habitats. Our results further support movements of trees in eastern North America in response to climatic changes. Finally, our study demonstrates that RAD-seq data and a combination of population genomics and SDM are valuable in resolving relationship and biogeographic history of closely related species that are taxonomically difficult.

Genome-wide genetic variation and comparison of fruit-associated traits between kumquat (Citrus japonica) and Clementine mandarin (Citrus clementina)

The comprehensive genetic variation of two citrus species were analyzed at genome and transcriptome level. A total of 1090 differentially expressed genes were found during fruit development by RNA-sequencing. Fruit size (fruit equatorial diameter) and weight (fresh weight) are the two most important components determining yield and consumer acceptability for many horticultural crops. However, little is known about the genetic control of these traits. Here, we performed whole-genome resequencing to reveal the comprehensive genetic variation of the fruit development between kumquat (Citrus japonica) and Clementine mandarin (Citrus clementina). In total, 5,865,235 single-nucleotide polymorphisms (SNPs) and 414,447 insertions/deletions (InDels) were identified in the two citrus species. Based on integrative analysis of genome and transcriptome of fruit, 640,801 SNPs and 20,733 InDels were identified. The features, genomic distribution, functional effect, and other characteristics of these genetic variations were explored. RNA-sequencing identified 1090 differentially expressed genes (DEGs) during fruit development of kumquat and Clementine mandarin. Gene Ontology revealed that these genes were involved in various molecular functional and biological processes. In addition, the genetic variation of 939 DEGs and 74 multiple fruit development pathway genes from previous reports were also identified. A global survey identified 24,237 specific alternative splicing events in the two citrus species and showed that intron retention is the most prevalent pattern of alternative splicing. These genome variation data provide a foundation for further exploration of citrus diversity and gene-phenotype relationships and for future research on molecular breeding to improve kumquat, Clementine mandarin and related species.

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.

MiSNPDb: a web-based genomic resources of tropical ecology fruit mango (Mangifera indica L.) for phylogeography and varietal differentiation

Mango is one of the most important fruits of tropical ecological region of the world, well known for its nutritive value, aroma and taste. Its world production is >45MT worth >200 billion US dollars. Genomic resources are required for improvement in productivity and management of mango germplasm. There is no web-based genomic resources available for mango. Hence rapid and cost-effective high throughput putative marker discovery is required to develop such resources. RAD-based marker discovery can cater this urgent need till whole genome sequence of mango becomes available. Using a panel of 84 mango varieties, a total of 28.6 Gb data was generated by ddRAD-Seq approach on Illumina HiSeq 2000 platform. A total of 1.25 million SNPs were discovered. Phylogenetic tree using 749 common SNPs across these varieties revealed three major lineages which was compared with geographical locations. A web genomic resources MiSNPDb, available at http://webtom.cabgrid.res.in/mangosnps/ is based on 3-tier architecture, developed using PHP, MySQL and Javascript. This web genomic resources can be of immense use in the development of high density linkage map, QTL discovery, varietal differentiation, traceability, genome finishing and SNP chip development for future GWAS in genomic selection program. We report here world’s first web-based genomic resources for genetic improvement and germplasm management of mango.

Identification of massive molecular markers in Echinochloa phyllopogon using a restriction-site associated DNA approach

Echinochloa phyllopogon proliferation seriously threatens rice production worldwide. We combined a restriction-site associated DNA (RAD) approach with Illumina DNA sequencing for rapid and mass discovery of simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers for E. phyllopogon. RAD tags were generated from the genomic DNA of two E. phyllopogon plants, and sequenced to produce 5197.7 Mb and 5242.9 Mb high quality sequences, respectively. The GC content of E. phyllopogon was 45.8%, which is high for monocots. In total, 4710 putative SSRs were identified in 4132 contigs, which permitted the design of PCR primers for E. phyllopogon. Most repeat motifs among the SSRs identified were dinucleotide (>82%), and most of these SSRs were four motif-repeats (>75%). The most frequent motif was AT, accounting for 36.3%-37.2%, followed by AG and AC. In total, 78 putative polymorphic SSR loci were found. A total of 49,179 SNPs were discovered between the two samples of E. phyllopogon, 67.1% of which were transversions and 32.9% were transitions. We used eight SSRs to study the genetic diversity of four E. phyllopogon populations collected from rice fields in China and all eight loci tested were polymorphic.

Toward Identification of Black Lemma and Pericarp Gene Blp1 in Barley Combining Bulked Segregant Analysis and Specific-Locus Amplified Fragment Sequencing

Black barley is caused by phytomelanin synthesized in lemma and/or pericarp and the trait is controlled by one dominant gene Blp1. The gene is mapped on chromosome 1H by molecular markers, but it is yet to be isolated. Specific-locus amplified fragment sequencing (SLAF-seq) is an effective method for large-scale de novo single nucleotide polymorphism (SNP) discovery and genotyping. In the present study, SLAF-seq with bulked segregant analysis (BSA) was employed to obtain sufficient markers to fine mapping Blp1 gene in an F2 population derived from Hatiexi No.1 × Zhe5819. Based on SNP screening criteria, a total of 77,542 polymorphic SNPs met the requirements for association analysis. Combining two association analysis methods, the overlapped region with a size of 32.41 Mb on chromosome 1H was obtained as the candidate region of Blp1 gene. According to SLAF-seq data, markers were developed in the target region and were used for mapping the Blp1 gene. Linkage analysis showed that Blp1 co-segregated with HZSNP34 and HZSNP36, and was delimited by two markers (HZSNP35 and HZSNP39) spanning 8.1 cM in 172 homozygous yellow grain F2 plants of Hatiexi No.1 × Zhe5819. More polymorphic markers were screened in the reduced target region and were used to genotype the population. As a result, Blp1 was delimited within a 1.66 Mb on chromosome 1H by the upstream marker HZSNP63 and the downstream marker HZSNP59. Our results demonstrated the utility of SLAF-seq-BSA approach to identify the candidate region and discover polymorphic markers at the specific targeted genomic region.

Mining sequence variations in representative polyploid sugarcane germplasm accessions

Background

Sugarcane (Saccharum spp.) is one of the most important economic crops because of its high sugar production and biofuel potential. Due to the high polyploid level and complex genome of sugarcane, it has been a huge challenge to investigate genomic sequence variations, which are critical for identifying alleles contributing to important agronomic traits. In order to mine the genetic variations in sugarcane, genotyping by sequencing (GBS), was used to genotype 14 representative Saccharum complex accessions. GBS is a method to generate a large number of markers, enabled by next generation sequencing (NGS) and the genome complexity reduction using restriction enzymes.

Results

To use GBS for high throughput genotyping highly polyploid sugarcane, the GBS analysis pipelines in 14 Saccharum complex accessions were established by evaluating different alignment methods, sequence variants callers, and sequence depth for single nucleotide polymorphism (SNP) filtering. By using the established pipeline, a total of 76,251 non-redundant SNPs, 5642 InDels, 6380 presence/absence variants (PAVs), and 826 copy number variations (CNVs) were detected among the 14 accessions. In addition, non-reference based universal network enabled analysis kit and Stacks de novo called 34,353 and 109,043 SNPs, respectively. In the 14 accessions, the percentages of single dose SNPs ranged from 38.3% to 62.3% with an average of 49.6%, much more than the portions of multiple dosage SNPs. Concordantly called SNPs were used to evaluate the phylogenetic relationship among the 14 accessions. The results showed that the divergence time between the Erianthus genus and the Saccharum genus was more than 10 million years ago (MYA). The Saccharum species separated from their common ancestors ranging from 0.19 to 1.65 MYA.

Conclusions

The GBS pipelines including the reference sequences, alignment methods, sequence variant callers, and sequence depth were recommended and discussed for the Saccharum complex and other related species. A large number of sequence variations were discovered in the Saccharum complex, including SNPs, InDels, PAVs, and CNVs. Genome-wide SNPs were further used to illustrate sequence features of polyploid species and demonstrated the divergence of different species in the Saccharum complex. The results of this study showed that GBS was an effective NGS-based method to discover genomic sequence variations in highly polyploid and heterozygous species.

Quantitative trait loci from identification to exploitation for crop improvement

Advancement in the field of genetics and genomics after the discovery of Mendel's laws of inheritance has led to map the genes controlling qualitative and quantitative traits in crop plant species. Mapping of genomic regions controlling the variation of quantitatively inherited traits has become routine after the advent of different types of molecular markers. Recently, the next generation sequencing methods have accelerated the research on QTL analysis. These efforts have led to the identification of more closely linked molecular markers with gene/QTLs and also identified markers even within gene/QTL controlling the trait of interest. Efforts have also been made towards cloning gene/QTLs or identification of potential candidate genes responsible for a trait. Further new concepts like crop QTLome and QTL prioritization have accelerated precise application of QTLs for genetic improvement of complex traits. In the past years, efforts have also been made in exploitation of a number of QTL for improving grain yield or other agronomic traits in various crops through markers assisted selection leading to cultivation of these improved varieties at farmers' field. In present article, we reviewed QTLs from their identification to exploitation in plant breeding programs and also reviewed that how improved cultivars developed through introgression of QTLs have improved the yield productivity in many crops.

Construction of a high-density genetic map for grape using specific length amplified fragment (SLAF) sequencing

Genetic maps are important tools in plant genomics and breeding. We report a large-scale discovery of single nucleotide polymorphisms (SNPs) using the specific length amplified fragment sequencing (SLAF-seq) technique for the construction of high-density genetic maps for two elite wine grape cultivars, ‘Chardonnay’ and ‘Beibinghong’, and their 130 F₁ plants. A total of 372.53 M paired-end reads were obtained after preprocessing. The average sequencing depth was 33.81 for ‘Chardonnay’ (the female parent), 48.20 for ‘Beibinghong’ (the male parent), and 12.66 for the F₁ offspring. We detected 202,349 high-quality SLAFs of which 144,972 were polymorphic; 10,042 SNPs were used to construct a genetic map that spanned 1,969.95 cM, with an average genetic distance of 0.23 cM between adjacent markers. This genetic map contains the largest molecular marker number of the grape maps so far reported. We thus demonstrate that SLAF-seq is a promising strategy for the construction of high-density genetic maps; the map that we report here is a good potential resource for QTL mapping of genes linked to major economic and agronomic traits, map-based cloning, and marker-assisted selection of grape.

Comparative proteomic analysis of eggplant (Solanum melongena L.) heterostylous pistil development

Heterostyly is a common floral polymorphism, but the proteomic basis of this trait is still largely unexplored. In this study, self- and cross-pollination of L-morph and S-morph flowers and comparison of embryo sac development in eggplant (Solanum melongena L.) suggested that lower fruit set from S-morph flowers results from stigma-pollen incompatibility. To explore the molecular mechanism underlying heterostyly development, we conducted isobaric tags for relative and absolute quantification (iTRAQ) proteomic analysis of eggplant pistils for L- and S-morph flowers. A total of 5,259 distinct proteins were identified during heterostyly development. Compared S-morph flowers with L-morph, we discovered 57 and 184 differentially expressed proteins (DEPs) during flower development and maturity, respectively. Quantitative real time polymerase chain reactions were used for nine genes to verify DEPs from the iTRAQ approach. During flower development, DEPs were mainly involved in morphogenesis, biosynthetic processes, and metabolic pathways. At flower maturity, DEPs primarily participated in biosynthetic processes, metabolic pathways, and the formation of ribosomes and proteasomes. Additionally, some proteins associated with senescence and programmed cell death were found to be upregulated in S-morph pistils, which may lead to the lower fruit set in S-morph flowers. Although the exact roles of these related proteins are not yet known, this was the first attempt to use an iTRAQ approach to analyze proteomes of heterostylous eggplant flowers, and these results will provide insights into biochemical events taking place during the development of heterostyly.

Footprints of domestication revealed by RAD-tag resequencing in loquat: SNP data reveals a non-significant domestication bottleneck and a single domestication event

Background

The process of crop domestication has long been a major area of research to gain insights into the history of human civilization and to understand the process of evolution. Loquat (Eriobotrya japonica Lindl.) is one of the typical subtropical fruit trees, which was domesticated in China at least 2000 years ago. In the present study, we re-sequenced the genome of nine wild loquat accessions collected from wide geographical range and 10 representative cultivated loquat cultivars by using RAD-tag tacit to exploit the molecular footprints of domestication.

Results

We obtained 26.4 Gb clean sequencing data from 19 loquat accessions, with an average of 32.64 M reads per genotype. We identified more than 80,000 SNPs distributed throughout the loquat genome. The SNP density and numbers were slightly higher in the wild loquat populations than that in the cultivated populations. All cultivars were clustered together by structure, phylogenetic and PCA analyses.

Conclusion

The modern loquat cultivars have experienced a non-significant genetic bottleneck during domestication, and originated from a single domesticated event. Moreover, our study revealed that Hubei province of China is probably the origin center of cultivated loquat.

Electronic supplementary material

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

A user guide to the Brassica 60K Illumina Infinium™ SNP genotyping array

The Brassica napus 60K Illumina Infinium™ SNP array has had huge international uptake in the rapeseed community due to the revolutionary speed of acquisition and ease of analysis of this high-throughput genotyping data, particularly when coupled with the newly available reference genome sequence. However, further utilization of this valuable resource can be optimized by better understanding the promises and pitfalls of SNP arrays. We outline how best to analyze Brassica SNP marker array data for diverse applications, including linkage and association mapping, genetic diversity and genomic introgression studies. We present data on which SNPs are locus-specific in winter, semi-winter and spring B. napus germplasm pools, rather than amplifying both an A-genome and a C-genome locus or multiple loci. Common issues that arise when analyzing array data will be discussed, particularly those unique to SNP markers and how to deal with these for practical applications in Brassica breeding applications.

Reduced representation genome sequencing reveals patterns of genetic diversity and selection in apple

Identifying DNA sequence variations is a fundamental step towards deciphering the genetic basis of traits of interest. Here, a total of 20 cultivated and 10 wild apples were genotyped using specific-locus amplified fragment sequencing, and 39,635 single nucleotide polymorphisms with no missing genotypes and evenly distributed along the genome were selected to investigate patterns of genome-wide genetic variations between cultivated and wild apples. Overall, wild apples displayed higher levels of genetic diversity than cultivated apples. Linkage disequilibrium (LD) decays were observed quite rapidly in cultivated and wild apples, with an r² -value below 0.2 at 440 and 280 bp, respectively. Moreover, bidirectional gene flow and different distribution patterns of LD blocks were detected between domesticated and wild apples. Most LD blocks unique to cultivated apples were located within QTL regions controlling fruit quality, thus suggesting that fruit quality had probably undergone selection during apple domestication. The genome of the earliest cultivated apple in China, Nai, was highly similar to that of Malus sieversii, and contained a small portion of genetic material from other wild apple species. This suggested that introgression could have been an important driving force during initial domestication of apple. These findings will facilitate future breeding and genetic dissection of complex traits in apple.

Genotyping-by-sequencing provides the first well-resolved phylogeny for coffee (Coffea) and insights into the evolution of caffeine content in its species: GBS coffee phylogeny and the evolution of caffeine content

A comprehensive and meaningful phylogenetic hypothesis for the commercially important coffee genus (Coffea) has long been a key objective for coffee researchers. For molecular studies, progress has been limited by low levels of sequence divergence, leading to insufficient topological resolution and statistical support in phylogenetic trees, particularly for the major lineages and for the numerous species occurring in Madagascar. We report here the first almost fully resolved, broadly sampled phylogenetic hypothesis for coffee, the result of combining genotyping-by-sequencing (GBS) technology with a newly developed, lab-based workflow to integrate short read next-generation sequencing for low numbers of additional samples. Biogeographic patterns indicate either Africa or Asia (or possibly the Arabian Peninsula) as the most likely ancestral locality for the origin of the coffee genus, with independent radiations across Africa, Asia, and the Western Indian Ocean Islands (including Madagascar and Mauritius). The evolution of caffeine, an important trait for commerce and society, was evaluated in light of our phylogeny. High and consistent caffeine content is found only in species from the equatorial, fully humid environments of West and Central Africa, possibly as an adaptive response to increased levels of pest predation. Moderate caffeine production, however, evolved at least one additional time recently (between 2 and 4Mya) in a Madagascan lineage, which suggests that either the biosynthetic pathway was already in place during the early evolutionary history of coffee, or that caffeine synthesis within the genus is subject to convergent evolution, as is also the case for caffeine synthesis in coffee versus tea and chocolate.

Characterization of Leaf Transcriptome in Banksia hookeriana

Banksia is a significant element in vegetation of southwestern Australia, a biodiversity hotspot with global significance. In particular, Banksia hookeriana represents a species with significant economic and ecological importance in the region. For better conservation and management, we reported an overview of transcriptome of B. hookeriana using RNA-seq and de novo assembly. We have generated a total of 202.7 million reads (18.91 billion of nucleotides) from four leaf samples in four plants of B. hookeriana, and assembled 59,063 unigenes (average size = 1098 bp) through de novotranscriptome assembly. Among them, 39,686 unigenes were annotated against the Swiss-Prot, Clusters of Orthologous Groups (COG), and NCBI non-redundant (NR) protein databases. We showed that there was approximately one single nucleotide polymorphism (SNP) per 5.6–7.1 kb in the transcriptome, and the ratio of transitional to transversional polymorphisms was approximately 1.82. We compared unigenes of B. hookeriana to those of Arabidopsis thaliana and Nelumbo nucifera through sequence homology, Gene Ontology (GO) annotation, and KEGG pathway analyses. The comparative analysis revealed that unigenes of B. hookeriana were closely related to those of N. nucifera. B. hookeriana, N. nucifera, and A. thaliana shared similar GO annotations but different distributions in KEGG pathways, indicating that B. hookeriana has adapted to dry-Mediterranean type shrublands via regulating expression of specific genes. In total 1927 potential simple sequence repeat (SSR) markers were discovered, which could be used in the genotype and genetic diversity studies of the Banksia genus. Our results provide valuable sequence resource for further study in Banksia.

Genotyping-by-sequencing approaches to characterize crop genomes: choosing the right tool for the right application

In the last decade, the revolution in sequencing technologies has deeply impacted crop genotyping practice. New methods allowing rapid, high-throughput genotyping of entire crop populations have proliferated and opened the door to wider use of molecular tools in plant breeding. These new genotyping-by-sequencing (GBS) methods include over a dozen reduced-representation sequencing (RRS) approaches and at least four whole-genome resequencing (WGR) approaches. The diversity of methods available, each often producing different types of data at different cost, can make selection of the best-suited method seem a daunting task. We review the most common genotyping methods used today and compare their suitability for linkage mapping, genomewide association studies (GWAS), marker-assisted and genomic selection and genome assembly and improvement in crops with various genome sizes and complexity. Furthermore, we give an outline of bioinformatics tools for analysis of genotyping data. WGR is well suited to genotyping biparental cross populations with complex, small- to moderate-sized genomes and provides the lowest cost per marker data point. RRS approaches differ in their suitability for various tasks, but demonstrate similar costs per marker data point. These approaches are generally better suited for de novo applications and more cost-effective when genotyping populations with large genomes or high heterozygosity. We expect that although RRS approaches will remain the most cost-effective for some time, WGR will become more widespread for crop genotyping as sequencing costs continue to decrease.

SNP-Discovery by RAD-Sequencing in a Germplasm Collection of Wild and Cultivated Grapevines (V. vinifera L.)

Whole-genome comparisons of Vitis vinifera subsp. sativa and V. vinifera subsp. sylvestris are expected to provide a better estimate of the valuable genetic diversity still present in grapevine, and help to reconstruct the evolutionary history of a major crop worldwide. To this aim, the increase of molecular marker density across the grapevine genome is fundamental. Here we describe the SNP discovery in a grapevine germplasm collection of 51 cultivars and 44 wild accessions through a novel protocol of restriction-site associated DNA (RAD) sequencing. By resequencing 1.1% of the grapevine genome at a high coverage, we recovered 34K BamHI unique restriction sites, of which 6.8% were absent in the ‘PN40024’ reference genome. Moreover, we identified 37,748 single nucleotide polymorphisms (SNPs), 93% of which belonged to the 19 assembled chromosomes with an average of 1.8K SNPs per chromosome. Nearly half of the SNPs fell in genic regions mostly assigned to the functional categories of metabolism and regulation, whereas some nonsynonymous variants were identified in genes related with the detection and response to environmental stimuli. SNP validation was carried-out, showing the ability of RAD-seq to accurately determine genotypes in a highly heterozygous species. To test the usefulness of our SNP panel, the main diversity statistics were evaluated, highlighting how the wild grapevine retained less genetic variability than the cultivated form. Furthermore, the analysis of Linkage Disequilibrium (LD) in the two subspecies separately revealed how the LD decays faster within the domesticated grapevine compared to its wild relative. Being the first application of RAD-seq in a diverse grapevine germplasm collection, our approach holds great promise for exploiting the genetic resources available in one of the most economically important fruit crops.

RNA-Sequencing Reveals Biological Networks during Table Grapevine ('Fujiminori') Fruit Development

Grapevine berry development is a complex and genetically controlled process, with many morphological, biochemical and physiological changes occurring during the maturation process. Research carried out on grapevine berry development has been mainly concerned with wine grape, while barely focusing on table grape. 'Fujiminori' is an important table grapevine cultivar, which is cultivated in most provinces of China. In order to uncover the dynamic networks involved in anthocyanin biosynthesis, cell wall development, lipid metabolism and starch-sugar metabolism in 'Fujiminori' fruit, we employed RNA-sequencing (RNA-seq) and analyzed the whole transcriptome of grape berry during development at the expanding period (40 days after full bloom, 40DAF), véraison period (65DAF), and mature period (90DAF). The sequencing depth in each sample was greater than 12×, and the expression level of nearly half of the expressed genes were greater than 1. Moreover, greater than 64% of the clean reads were aligned to the Vitis vinifera reference genome, and 5,620, 3,381, and 5,196 differentially expressed genes (DEGs) were identified between different fruit stages, respectively. Results of the analysis of DEGs showed that the most significant changes in various processes occurred from the expanding stage to the véraison stage. The expression patterns of F3'H and F3'5'H were crucial in determining red or blue color of the fruit skin. The dynamic networks of cell wall development, lipid metabolism and starch-sugar metabolism were also constructed. A total of 4,934 SSR loci were also identified from 4,337 grapevine genes, which may be helpful for the development of phylogenetic analysis in grapevine and other fruit trees. Our work provides the foundation for developmental research of grapevine fruit as well as other non-climacteric fruits.

Genome-wide SNP discovery and population structure analysis in pepper (Capsicum annuum) using genotyping by sequencing

Background

Knowledge on population structure and genetic diversity in vegetable crops is essential for association mapping studies and genomic selection. Genotyping by sequencing (GBS) represents an innovative method for large scale SNP detection and genotyping of genetic resources. Herein we used the GBS approach for the genome-wide identification of SNPs in a collection of Capsicum spp. accessions and for the assessment of the level of genetic diversity in a subset of 222 cultivated pepper (Capsicum annum) genotypes.

Results

GBS analysis generated a total of 7,568,894 master tags, of which 43.4% uniquely aligned to the reference genome CM334. A total of 108,591 SNP markers were identified, of which 105,184 were in C. annuum accessions. In order to explore the genetic diversity of C. annuum and to select a minimal core set representing most of the total genetic variation with minimum redundancy, a subset of 222 C. annuum accessions were analysed using 32,950 high quality SNPs. Based on Bayesian and Hierarchical clustering it was possible to divide the collection into three clusters. Cluster I had the majority of varieties and landraces mainly from Southern and Northern Italy, and from Eastern Europe, whereas clusters II and III comprised accessions of different geographical origins. Considering the genome-wide genetic variation among the accessions included in cluster I, a second round of Bayesian (K = 3) and Hierarchical (K = 2) clustering was performed. These analysis showed that genotypes were grouped not only based on geographical origin, but also on fruit-related features.

Conclusions

GBS data has proven useful to assess the genetic diversity in a collection of C. annuum accessions. The high number of SNP markers, uniformly distributed on the 12 chromosomes, allowed the accessions to be distinguished according to geographical origin and fruit-related features. SNP markers and information on population structure developed in this study will undoubtedly support genome-wide association mapping studies and marker-assisted selection programs.

Electronic supplementary material

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

Development of SSR markers for a Tibetan medicinal plant, Lancea tibetica (Phrymaceae), based on RAD sequencing

PREMISE OF THE STUDY

Lancea tibetica (Phrymaceae), a Tibetan medicinal plant, is endemic to the Qinghai-Tibet Plateau. The over-exploitation of wild L. tibetica has led to the destruction of many populations. To enhance protection and management, biological research, especially population genetic studies, should be carried out on L. tibetica. Simple sequence repeat (SSR) markers of L. tibetica were developed to analyze population diversity.

METHODS AND RESULTS

Four thousand four hundred and forty-one SSR loci were identified for L. tibetica based on restriction-site associated DNA (RAD) sequencing on the Illumina HiSeq platform. One hundred SSR loci were arbitrarily selected for primer design, and 38 of them were successfully amplified. These markers were tested on 56 individuals from three populations of L. tibetica, and 10 markers displayed polymorphisms. The total number of alleles per locus ranged from three to eight, and observed and expected heterozygosities ranged from 0.200 to 1.000 and 0.683 to 0.879, respectively. We tested for cross-amplification of these 10 markers in the related species L. hirsuta and found that nine could be successfully amplified.

CONCLUSIONS

The SSR markers characterized here are the first to be developed and tested in L. tibetica. They will be useful for future population genetic studies on L. tibetica and closely related species.

Construction of high-resolution genetic maps of Zoysia matrella (L.) Merrill and applications to comparative genomic analysis and QTL mapping of resistance to fall armyworm

BACKGROUND

Zoysia matrella, widely used in lawns and sports fields, is of great economic and ecological value. Z. matrella is an allotetraploid species (2n = 4x = 40) in the genus zoysia under the subfamily Chloridoideae. Despite its ecological impacts and economic importance, the subfamily Chloridoideae has received little attention in genomics studies. As a result, limited genetic and genomic information are available for this subfamily, which have impeded progress in understanding evolutionary history of grasses in this important lineage. The lack of a high-resolution genetic map has hampered efforts to improve zoysiagrass using molecular genetic tools.

RESULTS

We used restriction site-associated DNA sequencing (RADSeq) approach and a segregating population developed from the cross between Z. matrella cultivars 'Diamond' and 'Cavalier' to construct high-resolution genetic maps of Z. matrella. The genetic map of Diamond consists of 2,375 Single Nucleotide Polymorphism (SNP) markers mapped on 20 linkage groups (LGs) with a total length of 1754.48 cM and an average distance between adjacent markers at 0.74 cM. The genetic map of Cavalier contains 3,563 SNP markers on 20 LGs, covering 1824.92 cM, with an average distance between adjacent markers at 0.51 cM. A higher level of genome collinearity between Z. matrella and rice than that between Z. matrella and sorghum was revealed by comparative genomic analysis. Pairwise comparison revealed that two independent nested chromosome fusion events occurred after Z. matrella and sorghum split from a common ancestor. The high-resolution linkage maps were applied into mapping QTLs associated with fall armyworm (FAW) resistance and six loci located on LGs 8 and 20 were detected to be significantly associated with FAW resistance.

CONCLUSION

The high-resolution linkage maps provide anchor points for comparative genomics analysis between Z. matrella and other grass species. Our comparative genomic analysis suggested that the chromosome number reduction from 12 to 10 had occurred independently via a single-step in the subfamilies Chloridoideae and Panicoideae. The high-resolution genetic maps provide an essential framework for mapping QTLs associated with economically and agronomically important traits. The major QTLs mapped on LG8 of the Cavalier map provide a starting point for cloning FAW resistance genes and further studies for a better understanding of FAW resistance in zoysiagrass.

Development of a universal and simplified ddRAD library preparation approach for SNP discovery and genotyping in angiosperm plants

BACKGROUND

The double digest restriction-site associated DNA sequencing technology (ddRAD-seq) is a reduced representation sequencing technology by sampling genome-wide enzyme loci developed on the basis of next-generation sequencing. ddRAD-seq has been widely applied to SNP marker development and genotyping on animals, especially on marine animals as the original ddRAD protocol is mainly built and trained based on animal data. However, wide application of ddRAD-seq technology in plant species has not been achieved so far. Here, we aim to develop an optimized ddRAD library preparation protocol be accessible to most angiosperm plant species without much startup pre-experiment and costs.

RESULTS

We first tested several combinations of enzymes by in silico analysis of 23 plant species covering 17 families of angiosperm and 1 family of bryophyta and found AvaII + MspI enzyme pair produced consistently higher number of fragments in a broad range of plant species. Then we removed two purifying and one quantifying steps of the original protocol, replaced expensive consumables and apparatuses by conventional experimental apparatuses. Besides, we shortened P1 adapter from 37 to 25 bp and designed a new barcode-adapter system containing 20 pairs of barcodes of varying length. This is an optimized ddRAD strategy for angiosperm plants that is economical, time-saving and requires little technical expertise or investment in laboratory equipment. We refer to this simplified protocol as MiddRAD and we demonstrated the utility and flexibility of our approach by resolving phylogenetic relationships of two genera of woody bamboos (Dendrocalamus and Phyllostachys). Overall our results provide empirical evidence for using this method on different model and non-model plants to produce consistent data.

CONCLUSIONS

As MiddRAD adopts an enzyme pair that works for a broad range of angiosperm plants, simplifies library constructing procedure and requires less DNA input, it will greatly facilitate designing a ddRAD project. Our optimization of this method may make ddRAD be widely used in fields of plant population genetics, phylogenetics, phylogeography and molecular breeding.