Building Ultra-High-Density Linkage Maps Based on Efficient Filtering of Trustable Markers.

The study is focused on addressing the problem of building genetic maps in the presence of ∼103-104 of markers per chromosome. We consider a spectrum of situations with intrachromosomal heterogeneity of recombination rate, different level of genotyping errors, and missing data. In the ideal scenario of the absence of errors and missing data, the majority of markers should appear as groups of cosegregating markers ("twins") representing no challenge for map construction. The central aspect of the proposed approach is to take into account the structure of the marker space, where each twin group (TG) and singleton markers are represented as points of this space. The confounding effect of genotyping errors and missing data leads to reduction of TG size, but upon a low level of these effects surviving TGs can still be used as a source of reliable skeletal markers. Increase in the level of confounding effects results in a considerable decrease in the number or even disappearance of usable TGs and, correspondingly, of skeletal markers. Here, we show that the paucity of informative markers can be compensated by detecting kernels of markers in the marker space using a clustering procedure, and demonstrate the utility of this approach for high-density genetic map construction on simulated and experimentally obtained genotyping datasets.

Genome-Wide Differentiation of Various Melon Horticultural Groups for Use in GWAS for Fruit Firmness and Construction of a High Resolution Genetic Map.

Melon (Cucumis melo L.) is a phenotypically diverse eudicot diploid (2n = 2x = 24) has climacteric and non-climacteric morphotypes and show wide variation for fruit firmness, an important trait for transportation and shelf life. We generated 13,789 SNP markers using genotyping-by-sequencing (GBS) and anchored them to chromosomes to understand genome-wide fixation indices (Fst) between various melon morphotypes and genomewide linkage disequilibrium (LD) decay. The FST between accessions of cantalupensis and inodorus was 0.23. The FST between cantalupensis and various agrestis accessions was in a range of 0.19-0.53 and between inodorus and agrestis accessions was in a range of 0.21-0.59 indicating sporadic to wide ranging introgression. The EM (Expectation Maximization) algorithm was used for estimation of 1436 haplotypes. Average genome-wide LD decay for the melon genome was noted to be 9.27 Kb. In the current research, we focused on the genome-wide divergence underlying diverse melon horticultural groups. A high-resolution genetic map with 7153 loci was constructed. Genome-wide segregation distortion and recombination rate across various chromosomes were characterized. Melon has climacteric and non-climacteric morphotypes and wide variation for fruit firmness, a very important trait for transportation and shelf life. Various levels of QTLs were identified with high to moderate stringency and linked to fruit firmness using both genome-wide association study (GWAS) and biparental mapping. Gene annotation revealed some of the SNPs are located in ß-D-xylosidase, glyoxysomal malate synthase, chloroplastic anthranilate phosphoribosyltransferase, and histidine kinase, the genes that were previously characterized for fruit ripening and softening in other crops.

Experimental evolution of recombination and crossover interference in Drosophila caused by directional selection for stress-related traits.

BACKGROUND: Population genetics predicts that tight linkage between new and/or pre-existing beneficial and deleterious alleles should decrease the efficiency of natural selection in finite populations. By decoupling beneficial and deleterious alleles and facilitating the combination of beneficial alleles, recombination accelerates the formation of high-fitness genotypes. This may impose indirect selection for increased recombination. Despite the progress in theoretical understanding, interplay between recombination and selection remains a controversial issue in evolutionary biology. Even less satisfactory is the situation with crossover interference, which is a deviation of double-crossover frequency in a pair of adjacent intervals from the product of recombination rates in the two intervals expected on the assumption of crossover independence. Here, we report substantial changes in recombination and interference in three long-term directional selection experiments with Drosophila melanogaster: for desiccation (~50 generations), hypoxia, and hyperoxia tolerance (>200 generations each). RESULTS: For all three experiments, we found a high interval-specific increase of recombination frequencies in selection lines (up to 40-50% per interval) compared to the control lines. We also discovered a profound effect of selection on interference as expressed by an increased frequency of double crossovers in selection lines. Our results show that changes in interference are not necessarily coupled with increased recombination. CONCLUSIONS: Our results support the theoretical predictions that adaptation to a new environment can promote evolution toward higher recombination. Moreover, this is the first evidence of selection for different recombination-unrelated traits potentially leading, not only to evolution toward increased crossover rates, but also to changes in crossover interference, one of the fundamental features of recombination.

Fast and accurate construction of ultra-dense consensus genetic maps using evolution strategy optimization.

Our aim was to develop a fast and accurate algorithm for constructing consensus genetic maps for chip-based SNP genotyping data with a high proportion of shared markers between mapping populations. Chip-based genotyping of SNP markers allows producing high-density genetic maps with a relatively standardized set of marker loci for different mapping populations. The availability of a standard high-throughput mapping platform simplifies consensus analysis by ignoring unique markers at the stage of consensus mapping thereby reducing mathematical complicity of the problem and in turn analyzing bigger size mapping data using global optimization criteria instead of local ones. Our three-phase analytical scheme includes automatic selection of ~100-300 of the most informative (resolvable by recombination) markers per linkage group, building a stable skeletal marker order for each data set and its verification using jackknife re-sampling, and consensus mapping analysis based on global optimization criterion. A novel Evolution Strategy optimization algorithm with a global optimization criterion presented in this paper is able to generate high quality, ultra-dense consensus maps, with many thousands of markers per genome. This algorithm utilizes "potentially good orders" in the initial solution and in the new mutation procedures that generate trial solutions, enabling to obtain a consensus order in reasonable time. The developed algorithm, tested on a wide range of simulated data and real world data (Arabidopsis), outperformed two tested state-of-the-art algorithms by mapping accuracy and computation time.

Evaluation of marker-assisted selection for the stripe rust resistance gene Yr15, introgressed from wild emmer wheat.

Stripe rust disease is caused by the fungus Puccinia striiformis f. sp. tritici and severely threatens wheat worldwide, repeatedly breaking resistance conferred by resistance genes and evolving more aggressive strains. Wild emmer wheat, Triticum dicoccoides, is an important source for novel stripe rust resistance (Yr) genes. Yr15, a major gene located on chromosome 1BS of T. dicoccoides, was previously reported to confer resistance to a broad spectrum of stripe rust isolates, at both seedling and adult plant stages. Introgressions of Yr15 into cultivated T. aestivum bread wheat and T. durum pasta wheat that began in the 1980s are widely used. In the present study, we aimed to validate SSR markers from the Yr15 region as efficient tools for marker-assisted selection (MAS) for introgression of Yr15 into wheat and to compare the outcome of gene introgression by MAS and by conventional phenotypic selection. Our findings establish the validity of MAS for introgression of Yr15 into wheat. We show that the size of the introgressed segment, defined by flanking markers, varies for both phenotypic selection and MAS. The genetic distance of the MAS marker from Yr15 and the number of backcross steps were the main factors affecting the length of the introgressed donor segments. Markers Xbarc8 and Xgwm493, which are the nearest flanking markers studied, were consistent and polymorphic in all 34 introgressions reported here and are therefore the most recommended markers for the introgression of Yr15 into wheat cultivars. Introgression directed by markers, rather than by phenotype, will facilitate simultaneous selection for multiple stripe rust resistant genes and will help to avoid escapees during the selection process.

High-resolution genetic map for understanding the effect of genome-wide recombination rate on nucleotide diversity in watermelon.

We used genotyping by sequencing to identify a set of 10,480 single nucleotide polymorphism (SNP) markers for constructing a high-resolution genetic map of 1096 cM for watermelon. We assessed the genome-wide variation in recombination rate (GWRR) across the map and found an association between GWRR and genome-wide nucleotide diversity. Collinearity between the map and the genome-wide reference sequence for watermelon was studied to identify inconsistency and chromosome rearrangements. We assessed genome-wide nucleotide diversity, linkage disequilibrium (LD), and selective sweep for wild, semi-wild, and domesticated accessions of Citrullus lanatus var. lanatus to track signals of domestication. Principal component analysis combined with chromosome-wide phylogenetic study based on 1563 SNPs obtained after LD pruning with minor allele frequency of 0.05 resolved the differences between semi-wild and wild accessions as well as relationships among worldwide sweet watermelon. Population structure analysis revealed predominant ancestries for wild, semi-wild, and domesticated watermelons as well as admixture of various ancestries that were important for domestication. Sliding window analysis of Tajima's D across various chromosomes was used to resolve selective sweep. LD decay was estimated for various chromosomes. We identified a strong selective sweep on chromosome 3 consisting of important genes that might have had a role in sweet watermelon domestication.

Construction of a microsatellites-based linkage map for the white grouper (Epinephelus aeneus).

The white grouper (Epinephelus aeneus) is a promising candidate for domestication and aquaculture due to its fast growth, excellent taste, and high market price. A linkage map is an essential framework for mapping quantitative trait loci for economic traits and the study of genome evolution. DNA of a single individual was deep-sequenced, and microsatellite markers were identified in 177 of the largest scaffolds of the sequence assembly. The success rate of developing polymorphic homologous markers was 94.9% compared with 63.1% of heterologous markers from other grouper species. Of the 12 adult mature fish present in the broodstock tank, two males and two females were identified as parents of the assigned offspring by parenthood analysis using 34 heterologous markers. A single full-sib family of 48 individuals was established for the construction of first-generation linkage maps based on genotyping data of 222 microsatellites. The markers were assigned to 24 linkage groups in accordance to the 24 chromosomal pairs. The female and male maps consisting of 203 and 202 markers spanned 1053 and 886 cM, with an average intermarker distance of 5.8 and 5.0 cM, respectively. Mapping of markers to linkage groups ends was enriched by using markers originating from scaffolds harboring telomeric repeat-containing RNA. Comparative mapping showed high synteny relationships among the white grouper, kelp grouper (E. bruneus), orange-spotted grouper (E. coioides), and Nile tilapia (Oreochromis niloticus). Thus, it would be useful to integrate the markers that were developed for different groupers, depending on sharing of sequence data, into a comprehensive consensus map.

Some ways to improve QTL mapping accuracy.

In this paper, we review some approaches for QTL mapping developed by our research group in collaboration with, following the recommendation of, or under inspiration of Moshe Soller. Specifically, we explain at a simple intuitive level the main principles and ideas of: (a) QTL mapping by fractioned DNA pooling, (b) increasing the detection power of QTL mapping (in the case of individual genotyping) by multiple-trait analysis, and (c) the role of variance-covariance effects in QTL mapping. On each of these themes we had long and deep discussions with Soller on the statistical aspects of the proposed procedures. We hope that together we made important contributions towards making QTL mapping procedures easier and more effective.

High-density genetic map of durum wheat x wild emmer wheat based on SSR and DArT markers.

A genetic linkage map of tetraploid wheat was constructed based on a cross between durum wheat [Triticum turgidum ssp. durum (Desf.) MacKey] cultivar Langdon and wild emmer wheat [T. turgidum ssp. dicoccoides (Körn.) Thell.] accession G18-16. One hundred and fifty-two single-seed descent derived F(6) recombinant inbred lines (RILs) were analyzed with a total of 690 loci, including 197 microsatellite and 493 DArT markers. Linkage analysis defined 14 linkage groups. Most markers were mapped to the B-genome (60%), with an average of 57 markers per chromosome and the remaining 40% mapped to the A-genome, with an average of 39 markers per chromosome. To construct a stabilized (skeleton) map, markers interfering with map stability were removed. The skeleton map consisted of 307 markers with a total length of 2,317 cM and average distance of 7.5 cM between adjacent markers. The length of individual chromosomes ranged between 112 cM for chromosome 4B to 217 cM for chromosome 3B. A fraction (30.1%) of the markers deviated significantly from the expected Mendelian ratios; clusters of loci showing distorted segregation were found on chromosomes 1A, 1BL, 2BS, 3B, and 4B. DArT markers showed high proportion of clustering, which may be indicative of gene-rich regions. Three hundred and fifty-two new DArT markers were mapped for the first time on the current map. This map provides a useful groundwork for further genetic analyses of important quantitative traits, positional cloning, and marker-assisted selection, as well as for genome comparative genomics and genome organization studies in wheat and other cereals.

Genetic dissection of intermated recombinant inbred lines using a new genetic map of maize.

A new genetic map of maize, ISU-IBM Map4, that integrates 2029 existing markers with 1329 new indel polymorphism (IDP) markers has been developed using intermated recombinant inbred lines (IRILs) from the intermated B73xMo17 (IBM) population. The website http://magi.plantgenomics.iastate.edu provides access to IDP primer sequences, sequences from which IDP primers were designed, optimized marker-specific PCR conditions, and polymorphism data for all IDP markers. This new gene-based genetic map will facilitate a wide variety of genetic and genomic research projects, including map-based genome sequencing and gene cloning. The mosaic structures of the genomes of 91 IRILs, an important resource for identifying and mapping QTL and eQTL, were defined. Analyses of segregation data associated with markers genotyped in three B73/Mo17-derived mapping populations (F2, Syn5, and IBM) demonstrate that allele frequencies were significantly altered during the development of the IBM IRILs. The observations that two segregation distortion regions overlap with maize flowering-time QTL suggest that the altered allele frequencies were a consequence of inadvertent selection. Detection of two-locus gamete disequilibrium provides another means to extract functional genomic data from well-characterized plant RILs.

Domestication quantitative trait loci in Triticum dicoccoides, the progenitor of wheat.

Wild emmer wheat, Triticum dicoccoides, is the progenitor of modern tetraploid and hexaploid cultivated wheats. Our objective was to map domestication-related quantitative trait loci (QTL) in T. dicoccoides. The studied traits include brittle rachis, heading date, plant height, grain size, yield, and yield components. Our mapping population was derived from a cross between T. dicoccoides and Triticum durum. Approximately 70 domestication QTL effects were detected, nonrandomly distributed among and along chromosomes. Seven domestication syndrome factors were proposed, each affecting 5-11 traits. We showed: (i) clustering and strong effects of some QTLs; (ii) remarkable genomic association of strong domestication-related QTLs with gene-rich regions; and (iii) unexpected predominance of QTL effects in the A genome. The A genome of wheat may have played a more important role than the B genome during domestication evolution. The cryptic beneficial alleles at specific QTLs derived from T. dicoccoides may contribute to wheat and cereal improvement.