Constructing large-scale genetic maps using an evolutionary strategy algorithm

Evolutionary genomics of socially polymorphic populations of Pogonomyrmex californicus

BACKGROUND

Social insects vary considerably in their social organization both between and within species. In the California harvester ant, Pogonomyrmex californicus (Buckley 1867), colonies are commonly founded and headed by a single queen (haplometrosis, primary monogyny). However, in some populations in California (USA), unrelated queens cooperate not only during founding (pleometrosis) but also throughout the life of the colony (primary polygyny). The genetic architecture and evolutionary dynamics of this complex social niche polymorphism (haplometrosis vs pleometrosis) have remained unknown.

RESULTS

We provide a first analysis of its genomic basis and evolutionary history using population genomics comparing individuals from a haplometrotic population to those from a pleometrotic population. We discovered a recently evolved (< 200 k years), 8-Mb non-recombining region segregating with the observed social niche polymorphism. This region shares several characteristics with supergenes underlying social polymorphisms in other socially polymorphic ant species. However, we also find remarkable differences from previously described social supergenes. Particularly, four additional genomic regions not in linkage with the supergene show signatures of a selective sweep in the pleometrotic population. Within these regions, we find for example genes crucial for epigenetic regulation via histone modification (chameau) and DNA methylation (Dnmt1).

CONCLUSIONS

Altogether, our results suggest that social morph in this species is a polygenic trait involving a potential young supergene. Further studies targeting haplo- and pleometrotic individuals from a single population are however required to conclusively resolve whether these genetic differences underlie the alternative social phenotypes or have emerged through genetic drift.

QTL and PACE analyses identify candidate genes for anthracnose resistance in tomato

Anthracnose, caused by the fungal pathogen Colletotrichum spp., is one of the most significant tomato diseases in the United States and worldwide. No commercial cultivars with anthracnose resistance are available, limiting resistant breeding. Cultivars with genetic resistance would significantly reduce crop losses, reduce the use of fungicides, and lessen the risks associated with chemical application. A recombinant inbred line (RIL) mapping population (N=243) has been made from a cross between the susceptible US28 cultivar and the resistant but semiwild and small-fruited 95L368 to identify quantitative trait loci (QTLs) associated with anthracnose resistance. The RIL population was phenotyped for resistance by inoculating ripe field-harvested tomato fruits with Colletotrichum coccodes for two seasons. In this study, we identified twenty QTLs underlying resistance, with a range of phenotypic variance of 4.5 to 17.2% using a skeletal linkage map and a GWAS. In addition, a QTLseq analysis was performed using deep sequencing of extreme bulks that validated QTL positions identified using traditional mapping and resolved candidate genes underlying various QTLs. We further validated AP2-like ethylene-responsive transcription factor, N-alpha-acetyltransferase (NatA), cytochrome P450, amidase family protein, tetratricopeptide repeat, bHLH transcription factor, and disease resistance protein RGA2-like using PCR allelic competitive extension (PACE) genotyping. PACE assays developed in this study will enable high-throughput screening for use in anthracnose resistance breeding in tomato.

Progressive Learning Hill Climbing Algorithm with Energy-Map-Based Initialization for Image Reconstruction

Image reconstruction is an interesting yet challenging optimization problem that has several potential applications. The task is to reconstruct an image using a fixed number of transparent polygons. Traditional gradient-based algorithms cannot be applied to the problem since the optimization objective has no explicit expression and cannot be represented by computational graphs. Metaheuristic search algorithms are powerful optimization techniques for solving complex optimization problems, especially in the context of incomplete information or limited computational capability. In this paper, we developed a novel metaheuristic search algorithm named progressive learning hill climbing (ProHC) for image reconstruction. Instead of placing all the polygons on a blank canvas at once, ProHC starts from one polygon and gradually adds new polygons to the canvas until reaching the number limit. Furthermore, an energy-map-based initialization operator was designed to facilitate the generation of new solutions. To assess the performance of the proposed algorithm, we constructed a benchmark problem set containing four different types of images. The experimental results demonstrated that ProHC was able to produce visually pleasing reconstructions of the benchmark images. Moreover, the time consumed by ProHC was much shorter than that of the existing approach.

The First Genetic Map for a Psoraleoid Legume (Bituminaria bituminosa) Reveals Highly Conserved Synteny with Phaseoloid Legumes

We present the first genetic map of tedera (Bituminaria bituminosa (L.) C.H. Stirton), a drought-tolerant forage legume from the Canary Islands with useful pharmaceutical properties. It is also the first genetic map for any species in the tribe Psoraleeae (Fabaceae). The map comprises 2042 genotyping-by-sequencing (GBS) markers distributed across 10 linkage groups, consistent with the haploid chromosome count for this species (n = 10). Sequence tags from the markers were used to find homologous matches in the genome sequences of the closely related species in the Phaseoleae tribe: soybean, common bean, and cowpea. No tedera linkage groups align in their entirety to chromosomes in any of these phaseoloid species, but there are long stretches of collinearity that could be used in tedera research for gene discovery purposes using the better-resourced phaseoloid species. Using Ks analysis of a tedera transcriptome against five legume genomes provides an estimated divergence time of 17.4 million years between tedera and soybean. Genomic information and resources developed here will be invaluable for breeding tedera varieties for forage and pharmaceutical purposes.

Grain protein content and thousand kernel weight QTLs identified in a durum × wild emmer wheat mapping population tested in five environments

Genetic dissection of GPC and TKW in tetraploid durum × WEW RIL population, based on high-density SNP genetic map, revealed 12 GPC QTLs and 11 TKW QTLs, with favorable alleles for 11 and 5 QTLs, respectively, derived from WEW. Wild emmer wheat (Triticum turgidum ssp. dicoccoides, WEW) was shown to exhibit high grain protein content (GPC) and therefore possess a great potential for improvement of cultivated wheat nutritional value. Genetic dissection of thousand kernel weight (TKW) and grain protein content (GPC) was performed using a high-density genetic map constructed based on a recombinant inbred line (RIL) population derived from a cross between T. durum var. Svevo and WEW acc. Y12-3. Genotyping of 208 F₆ RILs with a 15 K wheat single nucleotide polymorphism (SNP) array yielded 4166 polymorphic SNP markers, of which 1510 were designated as skeleton markers. A total map length of 2169 cM was obtained with an average distance of 1.5 cM between SNPs. A total of 12 GPC QTLs and 11 TKW QTLs were found under five different environments. No significant correlations were found between GPC and TKW across all environments. Four major GPC QTLs with favorable alleles from WEW were found on chromosomes 4BS, 5AS, 6BS and 7BL. The 6BS GPC QTL coincided with the physical position of the NAC transcription factor TtNAM-B1, underlying the cloned QTL, Gpc-B1. Comparisons of the physical intervals of the GPC QTLs described here with the results previously reported in other durum × WEW RIL population led to the discovery of seven novel GPC QTLs. Therefore, our research emphasizes the importance of GPC QTL dissection in diverse WEW accessions as a source of novel alleles for improvement of GPC in cultivated wheat.

Grain protein content and thousand kernel weight QTLs identified in a durum × wild emmer wheat mapping population tested in five environments

Genetic dissection of GPC and TKW in tetraploid durum × WEW RIL population, based on high-density SNP genetic map, revealed 12 GPC QTLs and 11 TKW QTLs, with favorable alleles for 11 and 5 QTLs, respectively, derived from WEW. Wild emmer wheat (Triticum turgidum ssp. dicoccoides, WEW) was shown to exhibit high grain protein content (GPC) and therefore possess a great potential for improvement of cultivated wheat nutritional value. Genetic dissection of thousand kernel weight (TKW) and grain protein content (GPC) was performed using a high-density genetic map constructed based on a recombinant inbred line (RIL) population derived from a cross between T. durum var. Svevo and WEW acc. Y12-3. Genotyping of 208 F₆ RILs with a 15 K wheat single nucleotide polymorphism (SNP) array yielded 4166 polymorphic SNP markers, of which 1510 were designated as skeleton markers. A total map length of 2169 cM was obtained with an average distance of 1.5 cM between SNPs. A total of 12 GPC QTLs and 11 TKW QTLs were found under five different environments. No significant correlations were found between GPC and TKW across all environments. Four major GPC QTLs with favorable alleles from WEW were found on chromosomes 4BS, 5AS, 6BS and 7BL. The 6BS GPC QTL coincided with the physical position of the NAC transcription factor TtNAM-B1, underlying the cloned QTL, Gpc-B1. Comparisons of the physical intervals of the GPC QTLs described here with the results previously reported in other durum × WEW RIL population led to the discovery of seven novel GPC QTLs. Therefore, our research emphasizes the importance of GPC QTL dissection in diverse WEW accessions as a source of novel alleles for improvement of GPC in cultivated wheat.

Characterization of a dominant mutation for the liguleless trait: Aegilops tauschii liguleless (Lgt)

BACKGROUND

Leaves of Poaceae have a unique morphological feature: they consist of a proximal sheath and a distal blade separated by a ligular region. The sheath provides structural support and protects young developing leaves, whereas the main function of the blade is photosynthesis. The auricles allow the blade to tilt back for optimal photosynthesis and determine the angle of a leaf, whereas the ligule protects the stem from the entry of water, microorganisms, and pests. Liguleless variants have an upright leaf blade that wraps around the culm. Research on liguleless mutants of maize and other cereals has led to identification of genes that are involved in leaf patterning and differentiation.

RESULTS

We characterized an induced liguleless mutant (LM) of Aegilops tauschii Coss., a donor of genome D of bread wheat Triticum aestivum L.. The liguleless phenotype of LM is under dominant monogenic control (Lg^t). To determine precise position of Lg^t on the Ae. tauschii genetic map, highly saturated genetic maps were constructed containing 887 single-nucleotide polymorphism (SNP) markers derived via diversity arrays technology (DArT)seq. The Lg^t gene was mapped to chromosome 5DS. Taking into account coordinates of the SNP markers, flanking Lg^t, on the pseudomolecule 5D, a chromosomal region that contains this gene was determined, and a list of candidate genes was identified. Morphological features of the LM phenotype suggest that Lg^t participates in the control of leaf development, mainly, in leaf proximal-distal patterning, and its dominant mutation causes abnormal ligular region but does not affect reproductive development.

CONCLUSIONS

Here we report characterization of a liguleless Ae. tauschii mutant, whose phenotype is under control of a dominant mutation of Lg^t. The dominant mode of inheritance of the liguleless trait in a Triticeae species is reported for the first time. The position of the Lg^t locus on chromosome 5DS allowed us to identify a list of candidate genes. This list does not contain Ae. tauschii orthologs of any well-characterized cereal genes whose mutations cause liguleless phenotypes. Thus, the characterized Lg^t mutant represents a new model for further investigation of plant leaf patterning and differentiation.

A transcriptome-SNP-derived linkage map of Apios americana (potato bean) provides insights about genome re-organization and synteny conservation in the phaseoloid legumes

KEY MESSAGE

We report a linkage map for Apios americana and describe synteny with selected warm-season legumes. A translocation event in common bean and soybean is confirmed against Apios and Vigna species. Apios (Apios americana; "apios"), a tuberous perennial legume in the Phaseoleae tribe, was widely used as a food by Native Americans. Work in the last 40 years has led to several improved breeding lines. Aspects of the pollination biology (complex floral structure and tripping mechanism) have made controlled crosses difficult, and the previous reports indicated that the plant is likely primarily an outcrosser. We used a pseudo-testcross strategy to construct a genetic map specific to the maternal parent. The map was built using single-nucleotide polymorphism markers identified by comparing the expressed sequences of individuals in the mapping population against a de novo maternal reference transcriptome assembly. The apios map consists of 11 linkage groups and 1121 recombinationally distinct loci, covering ~ 938.6 cM. By sequencing the transcriptomes of all potential pollen parents, we were able to identify the probable pollen donors and to discover new aspects of the pollination biology in apios. No selfing was observed, but multiple pollen parents were seen within individual pods. Comparisons with genome sequences in other species in the Phaseoleae showed extended synteny for most apios linkage groups. This synteny supports the robustness of the map, and also sheds light on the history of the Phaseoleae, as apios is relatively early diverging in this tribe. We detected a translocation event that separates apios and two Vigna species from Phaseolus vulgaris and Glycine max. This apios mapping work provides a general protocol for sequencing-based construction of high-density linkage maps in outcrossing species with heterogeneous pollen parents.

TSPmap, a tool making use of traveling salesperson problem solvers in the efficient and accurate construction of high-density genetic linkage maps

Background

Recent advances in nucleic acid sequencing technologies have led to a dramatic increase in the number of markers available to generate genetic linkage maps. This increased marker density can be used to improve genome assemblies as well as add much needed resolution for loci controlling variation in ecologically and agriculturally important traits. However, traditional genetic map construction methods from these large marker datasets can be computationally prohibitive and highly error prone.

Results

We present TSPmap, a method which implements both approximate and exact Traveling Salesperson Problem solvers to generate linkage maps. We demonstrate that for datasets with large numbers of genomic markers (e.g. 10,000) and in multiple population types generated from inbred parents, TSPmap can rapidly produce high quality linkage maps with low sensitivity to missing and erroneous genotyping data compared to two other benchmark methods, JoinMap and MSTmap. TSPmap is open source and freely available as an R package.

Conclusions

With the advancement of low cost sequencing technologies, the number of markers used in the generation of genetic maps is expected to continue to rise. TSPmap will be a useful tool to handle such large datasets into the future, quickly producing high quality maps using a large number of genomic markers.

Electronic supplementary material

The online version of this article (10.1186/s13040-017-0158-0) contains supplementary material, which is available to authorized users.

Effect of Co-segregating Markers on High-Density Genetic Maps and Prediction of Map Expansion Using Machine Learning Algorithms

Advances in sequencing and genotyping methods have enable cost-effective production of high throughput single nucleotide polymorphism (SNP) markers, making them the choice for linkage mapping. As a result, many laboratories have developed high-throughput SNP assays and built high-density genetic maps. However, the number of markers may, by orders of magnitude, exceed the resolution of recombination for a given population size so that only a minority of markers can accurately be ordered. Another issue attached to the so-called 'large p, small n' problem is that high-density genetic maps inevitably result in many markers clustering at the same position (co-segregating markers). While there are a number of related papers, none have addressed the impact of co-segregating markers on genetic maps. In the present study, we investigated the effects of co-segregating markers on high-density genetic map length and marker order using empirical data from two populations of wheat, Mohawk × Cocorit (durum wheat) and Norstar × Cappelle Desprez (bread wheat). The maps of both populations consisted of 85% co-segregating markers. Our study clearly showed that excess of co-segregating markers can lead to map expansion, but has little effect on markers order. To estimate the inflation factor (IF), we generated a total of 24,473 linkage maps (8,203 maps for Mohawk × Cocorit and 16,270 maps for Norstar × Cappelle Desprez). Using seven machine learning algorithms, we were able to predict with an accuracy of 0.7 the map expansion due to the proportion of co-segregating markers. For example in Mohawk × Cocorit, with 10 and 80% co-segregating markers the length of the map inflated by 4.5 and 16.6%, respectively. Similarly, the map of Norstar × Cappelle Desprez expanded by 3.8 and 11.7% with 10 and 80% co-segregating markers. With the increasing number of markers on SNP-chips, the proportion of co-segregating markers in high-density maps will continue to increase making map expansion unavoidable. Therefore, we suggest developers improve linkage mapping algorithms for efficient analysis of high-throughput data. This study outlines a practical strategy to estimate the IF due to the proportion of co-segregating markers and outlines a method to scale the length of the map accordingly.

Identification and evaluation of resistance to powdery mildew and yellow rust in a wheat mapping population

Deployment of cultivars with genetic resistance is an effective approach to control the diseases of powdery mildew (PM) and yellow rust (YR). Chinese wheat cultivar XK0106 exhibits high levels of resistance to both diseases, while cultivar E07901 has partial, adult plant resistance (APR). The aim of this study was to map resistance loci derived from the two cultivars and analyze their effects against PM and YR in a range of environments. A doubled haploid population (388 lines) was used to develop a framework map consisting of 117 SSR markers, while a much higher density map using the 90K Illumina iSelect SNP array was produced with a subset of 80 randomly selected lines. Seedling resistance was characterized against a range of PM and YR isolates, while field scores in multiple environments were used to characterize APR. Composite interval mapping (CIM) of seedling PM scores identified two QTLs (QPm.haas-6A and QPm.haas-2A), the former being located at the Pm21 locus. These QTLs were also significant in field scores, as were Qpm.haas-3A and QPm.haas-5A. QYr.haas-1B-1 and QYr.haas-2A were identified in field scores of YR and were located at the Yr24/26 and Yr17 chromosomal regions respectively. A second 1B QTL, QYr.haas-1B-2 was also identified. QPm.haas-2A and QYr.haas-1B-2 are likely to be new QTLs that have not been previously identified. Effects of the QTLs were further investigated in multiple environments through the testing of selected lines predicted to contain various QTL combinations. Significant additive interactions between the PM QTLs highlighted the ability to pyramid these loci to provide higher level of resistance. Interactions between the YR QTLs gave insights into the pathogen populations in the different locations as well as showing genetic interactions between these loci.

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 ∼10³-10⁴ 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.

QTL mapping identifies a major locus for resistance in wheat to Sunn pest (Eurygaster integriceps) feeding at the vegetative growth stage

This research provides the first report of a major locus controlling wheat resistance to Sunn pest. It developed and validated SNP markers that will be useful for marker-assisted selection. Sunn pest (Eurygaster integriceps Puton) is the most destructive insect pest of bread wheat and durum wheat in West and Central Asia and East Europe. Breeding for resistance at the vegetative stage of growth is vital in reducing the damage caused by overwintered adult populations that feed on shoot and leaves of seedlings, and in reducing the next generation of pest populations (nymphs and adults), which can cause damage to grain quality by feeding on spikes. In the present study, two doubled haploid (DH) populations involving resistant landraces from Afghanistan were genotyped with the 90k SNP iSelect assay and candidate gene-based KASP markers. The DH lines and parents were phenotyped for resistance to Sunn pest feeding, using artificial infestation cages at Terbol station, in Lebanon, over three years. Quantitative trait locus (QTL) analysis identified a single major locus on chromosome 4BS in the two populations, with the resistance allele derived from the landrace accessions, IG139431 and IG139883. The QTL explained a maximum of 42 % of the phenotypic variation in the Cham6 × IG139431 and 56 % in the Cham6 × IG139883 populations. SNP markers closest to the QTL showed high similarity to rice genes that putatively encode proteins for defense response to herbivory and wounding. The markers were validated in a large, unrelated population of parental wheat genotypes. All wheat lines carrying the 'C-G' haplotype at the identified SNPs were resistant, suggesting that selection based on a haplotype of favourable alleles would be effective in predicting resistance status of unknown genotypes.

The potential of probabilistic graphical models in linkage map construction

Probabilistic graphical models show great potential for robust and reliable construction of linkage maps. We show how to use probabilistic graphical models to construct high-quality linkage maps in the face of data perturbations caused by genotyping errors and reciprocal translocations. It has been shown that linkage map construction can be hampered by the presence of genotyping errors and chromosomal rearrangements such as inversions and translocations. Here, we report a novel method for linkage map construction using probabilistic graphical models. The method is proven, both theoretically and practically, to be effective in filtering out markers that contain genotyping errors. In particular, it carries out marker filtering and ordering simultaneously, and is therefore superior to the standard post hoc filtering using nearest-neighbour stress. Furthermore, we demonstrate empirically that the proposed method offers a promising solution to linkage map construction in the case of a reciprocal translocation.

A High Resolution Radiation Hybrid Map of Wheat Chromosome 4A

Bread wheat has a large and complex allohexaploid genome with low recombination level at chromosome centromeric and peri-centromeric regions. This significantly hampers ordering of markers, contigs of physical maps and sequence scaffolds and impedes obtaining of high-quality reference genome sequence. Here we report on the construction of high-density and high-resolution radiation hybrid (RH) map of chromosome 4A supported by high-density chromosome deletion map. A total of 119 endosperm-based RH lines of two RH panels and 15 chromosome deletion bin lines were genotyped with 90K iSelect single nucleotide polymorphism (SNP) array. A total of 2316 and 2695 markers were successfully mapped to the 4A RH and deletion maps, respectively. The chromosome deletion map was ordered in 19 bins and allowed precise identification of centromeric region and verification of the RH panel reliability. The 4A-specific RH map comprises 1080 mapping bins and spans 6550.9 cR with a resolution of 0.13 Mb/cR. Significantly higher mapping resolution in the centromeric region was observed as compared to recombination maps. Relatively even distribution of deletion frequency along the chromosome in the RH panel was observed and putative functional centromere was delimited within a region characterized by two SNP markers.

Identification and genetic mapping of PmAF7DS a powdery mildew resistance gene in bread wheat (Triticum aestivum L.)

Gene PmAF7DS confers resistance to wheat powdery mildew (isolate Bgt#211 ); it was mapped to a 14.6-cM interval ( Xgwm350 a- Xbarc184 ) on chromosome 7DS. The flanking markers could be applied in MAS breeding. Wheat powdery mildew (Pm) is caused by the biotrophic pathogen Blumeria graminis tritici (DC.) (Bgt). An ongoing threat of breakdown of race-specific resistance to Pm requires a continuous effort to discover new alleles in the wheat gene pool. Developing new cultivars with improved disease resistance is an economically and environmentally safe approach to reduce yield losses. To identify and characterize genes for resistance against Pm in bread wheat we used the (Arina × Forno) RILs population. Initially, the two parental lines were screened with a collection of 61 isolates of Bgt from Israel. Three Pm isolates Bgt#210 , Bgt#211 and Bgt#213 showed differential reactions in the parents: Arina was resistant (IT = 0), whereas Forno was moderately susceptible (IT = -3). Isolate Bgt#211 was then used to inoculate the RIL population. The segregation pattern of plant reactions among the RILs indicates that a single dominant gene controls the conferred resistance. A genetic map of the region containing this gene was assembled with DNA markers and assigned to the 7D physical bin map. The gene, temporarily designated PmAF7DS, was located in the distal region of chromosome arm 7DS. The RILs were also inoculated with Bgt#210 and Bgt#213. The plant reactions to these isolates showed high identity with the reaction to Bgt#211, indicating the involvement of the same gene or closely linked, but distinct single genes. The genomic location of PmAF7DS, in light of other Pm genes on 7DS is discussed.

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 F_ST between accessions of cantalupensis and inodorus was 0.23. The F_ST 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.

Genetic evidence for differential selection of grain and embryo weight during wheat evolution under domestication

Wheat is one of the Neolithic founder crops domesticated ~10 500 years ago. Following the domestication episode, its evolution under domestication has resulted in various genetic modifications. Grain weight, embryo weight, and the interaction between those factors were examined among domesticated durum wheat and its direct progenitor, wild emmer wheat. Experimental data show that grain weight has increased over the course of wheat evolution without any parallel change in embryo weight, resulting in a significantly reduced (30%) embryo weight/grain weight ratio in domesticated wheat. The genetic factors associated with these modifications were further investigated using a population of recombinant inbred substitution lines that segregated for chromosome 2A. A cluster of loci affecting grain weight and shape was identified on the long arm of chromosome 2AL. Interestingly, a novel locus controlling embryo weight was mapped on chromosome 2AS, on which the wild emmer allele promotes heavier embryos and greater seedling vigour. To the best of our knowledge, this is the first report of a QTL for embryo weight in wheat. The results suggest a differential selection of grain and embryo weight during the evolution of domesticated wheat. It is argued that conscious selection by early farmers favouring larger grains and smaller embryos appears to have resulted in a significant change in endosperm weight/embryo weight ratio in the domesticated wheat. Exposing the genetic factors associated with endosperm and embryo size improves our understanding of the evolutionary dynamics of wheat under domestication and is likely to be useful for future wheat-breeding efforts.

Characterizing uncertainty in high-density maps from multiparental populations

Multiparental populations are of considerable interest in high-density genetic mapping due to their increased levels of polymorphism and recombination relative to biparental populations. However, errors in map construction can have significant impact on QTL discovery in later stages of analysis, and few methods have been developed to quantify the uncertainty attached to the reported order of markers or intermarker distances. Current methods are computationally intensive or limited to assessing uncertainty only for order or distance, but not both simultaneously. We derive the asymptotic joint distribution of maximum composite likelihood estimators for intermarker distances. This approach allows us to construct hypothesis tests and confidence intervals for simultaneously assessing marker-order instability and distance uncertainty. We investigate the effects of marker density, population size, and founder distribution patterns on map confidence in multiparental populations through simulations. Using these data, we provide guidelines on sample sizes necessary to map markers at sub-centimorgan densities with high certainty. We apply these approaches to data from a bread wheat Multiparent Advanced Generation Inter-Cross (MAGIC) population genotyped using the Illumina 9K SNP chip to assess regions of uncertainty and validate them against the recently released pseudomolecule for the wheat chromosome 3B.

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.

Transcriptome sequencing of different narrow-leafed lupin tissue types provides a comprehensive uni-gene assembly and extensive gene-based molecular markers

Narrow-leafed lupin (NLL; Lupinus angustifolius L.) is an important grain legume crop that is valuable for sustainable farming and is becoming recognized as a human health food. NLL breeding is directed at improving grain production, disease resistance, drought tolerance and health benefits. However, genetic and genomic studies have been hindered by a lack of extensive genomic resources for the species. Here, the generation, de novo assembly and annotation of transcriptome datasets derived from five different NLL tissue types of the reference accession cv. Tanjil are described. The Tanjil transcriptome was compared to transcriptomes of an early domesticated cv. Unicrop, a wild accession P27255, as well as accession 83A:476, together being the founding parents of two recombinant inbred line (RIL) populations. In silico predictions for transcriptome-derived gene-based length and SNP polymorphic markers were conducted and corroborated using a survey assembly sequence for NLL cv. Tanjil. This yielded extensive indel and SNP polymorphic markers for the two RIL populations. A total of 335 transcriptome-derived markers and 66 BAC-end sequence-derived markers were evaluated, and 275 polymorphic markers were selected to genotype the reference NLL 83A:476 × P27255 RIL population. This significantly improved the completeness, marker density and quality of the reference NLL genetic map.

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.

QTL dissection of the loss of green colour during post-anthesis grain maturation in two-rowed barley

Ability to genetically manipulate the loss of green colour during grain maturation has potentials for increasing productivity, disease resistance, and drought and heat tolerance in crop plants. Two doubled haploid, two-rowed barley populations (Vlamingh × Buloke and VB9524 × ND11231*12) were monitored over 2 years for loss of green colour during grain filling using a portable active sensor. The aims were to determine the genomic regions that control trait heritability by quantitative trait locus (QTL) analysis, and to examine patterns of QTL-environment interactions under different conditions of water stress. In the Vlamingh × Buloke cross, broad-sense heritability estimate for loss of green colour (measured as the difference in sensor readings taken at anthesis and maturity, ∆SRI) was 0.68, and 0.78 for the VB9524 × ND11231*12 population. In the VB9524 × ND11231*12 population, rapid loss of green colour was positively associated with grain yield and percent plump grains, but in the Vlamingh × Buloke population, a slower loss of green colour (low ∆SRI) was associated with increased grain plumpness. With the aid of a dense array of single nucleotide polymorphisms (SNPs) and EST-derived SSR markers, a total of nine QTLs were detected across the two populations. Of these, a single major locus on the short arm of barley chromosome 5H was consistently linked with trait variation across the populations and multiple environments. The QTL was independent of flowering time and explained between 5.4 and 15.4 % of the variation observed in both populations, depending on the environment, and although a QTL × E interaction was detected, it was largely due to a change in the magnitude of the effect, rather than a change in direction. The results suggest that loss of green colour during grain maturation may be under the control of a simple genetic architecture, but a careful study of target populations and environments would be required for breeding purposes.

Genotype-specific SNP map based on whole chromosome 3B sequence information from wheat cultivars Arina and Forno

Agronomically important traits are frequently controlled by rare, genotype-specific alleles. Such genes can only be mapped in a population derived from the donor genotype. This requires the development of a specific genetic map, which is difficult in wheat because of the low level of polymorphism among elite cultivars. The absence of sufficient polymorphism, the complexity of the hexaploid wheat genome as well as the lack of complete sequence information make the construction of genetic maps with a high density of reproducible and polymorphic markers challenging. We developed a genotype-specific genetic map of chromosome 3B from winter wheat cultivars Arina and Forno. Chromosome 3B was isolated from the two cultivars and then sequenced to 10-fold coverage. This resulted in a single-nucleotide polymorphisms (SNP) database of the complete chromosome. Based on proposed synteny with the Brachypodium model genome and gene annotation, sequences close to coding regions were used for the development of 70 SNP-based markers. They were mapped on a Arina × Forno Recombinant Inbred Lines population and found to be spread over the complete chromosome 3B. While overall synteny was well maintained, numerous exceptions and inversions of syntenic gene order were identified. Additionally, we found that the majority of recombination events occurred in distal parts of chromosome 3B, particularly in hot-spot regions. Compared with the earlier map based on SSR and RFLP markers, the number of markers increased fourfold. The approach presented here allows fast development of genotype-specific polymorphic markers that can be used for mapping and marker-assisted selection.

Development of an ultra-dense genetic map of the sunflower genome based on single-feature polymorphisms

The development of ultra-dense genetic maps has the potential to facilitate detailed comparative genomic analyses and whole genome sequence assemblies. Here we describe the use of a custom Affymetrix GeneChip containing nearly 2.4 million features (25 bp sequences) targeting 86,023 unigenes from sunflower (Helianthus annuus L.) and related species to test for single-feature polymorphisms (SFPs) in a recombinant inbred line (RIL) mapping population derived from a cross between confectionery and oilseed sunflower lines (RHA280×RHA801). We then employed an existing genetic map derived from this same population to rigorously filter out low quality data and place 67,486 features corresponding to 22,481 unigenes on the sunflower genetic map. The resulting map contains a substantial fraction of all sunflower genes and will thus facilitate a number of downstream applications, including genome assembly and the identification of candidate genes underlying QTL or traits of interest.

Two-phase analysis in consensus genetic mapping

Numerous mapping projects conducted on different species have generated an abundance of mapping data. Consequently, many multilocus maps have been constructed using diverse mapping populations and marker sets for the same organism. The quality of maps varies broadly among populations, marker sets, and software used, necessitating efforts to integrate the mapping information and generate consensus maps. The problem of consensus genetic mapping (MCGM) is by far more challenging compared with genetic mapping based on a single dataset, which by itself is also cumbersome. The additional complications introduced by consensus analysis include inter-population differences in recombination rate and exchange distribution along chromosomes; variations in dominance of the employed markers; and use of different subsets of markers in different labs. Hence, it is necessary to handle arbitrary patterns of shared sets of markers and different level of mapping data quality. In this article, we introduce a two-phase approach for solving MCGM. In phase 1, for each dataset, multilocus ordering is performed combined with iterative jackknife resampling to evaluate the stability of marker orders. In this phase, the ordering problem is reduced to the well-known traveling salesperson problem (TSP). Namely, for each dataset, we look for order that gives minimum sum of recombination distances between adjacent markers. In phase 2, the optimal consensus order of shared markers is selected from the set of allowed orders and gives the minimal sum of total lengths of nonconflicting maps of the chromosome. This criterion may be used in different modifications to take into account the variation in quality of the original data (population size, marker quality, etc.). In the foregoing formulation, consensus mapping is considered as a specific version of TSP that can be referred to as “synchronized TSP.” The conflicts detected after phase 1 are resolved using either a heuristic algorithm over the entire chromosome or an exact/heuristic algorithm applied subsequently to the revealed small non-overlapping regions with conflicts separated by non-conflicting regions. The proposed approach was tested on a wide range of simulated data and real datasets from maize.

Identification and characterization of a novel powdery mildew resistance gene PmG3M derived from wild emmer wheat, Triticum dicoccoides

Powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt) is one of the most important wheat diseases worldwide. Wild emmer wheat, Triticum turgidum ssp. dicoccoides, the tetraploid ancestor (AABB) of domesticated bread and durum wheat, harbors many important alleles for resistance to various diseases, including powdery mildew. In the current study, two tetraploid wheat mapping populations, derived from a cross between durum wheat (cv. Langdon) and wild emmer wheat (accession G-305-3M), were used to identify and map a novel powdery mildew resistance gene. Wild emmer accession G-305-3M was resistant to all 47 Bgt isolates tested, from Israel and Switzerland. Segregation ratios of F(2) progenies and F(6) recombinant inbred line (RIL) mapping populations, in their reactions to inoculation with Bgt, revealed a Mendelian pattern (3:1 and 1:1, respectively), indicating the role of a single dominant gene derived from T. dicoccoides accession G-305-3M. This gene, temporarily designated PmG3M, was mapped on chromosome 6BL and physically assigned to chromosome deletion bin 6BL-0.70-1.00. The F(2) mapping population was used to construct a genetic map of the PmG3M gene region consisted of six simple sequence repeats (SSR), 11 resistance gene analog (RGA), and two target region amplification polymorphism (TRAP) markers. A second map, constructed based on the F(6) RIL population, using a set of skeleton SSR markers, confirmed the order of loci and distances obtained for the F(2) population. The discovery and mapping of this novel powdery mildew resistance gene emphasize the importance of the wild emmer wheat gene pool as a source for crop improvement.

Divergent functions of orthologous NAC transcription factors in wheat and rice

The wheat GPC-B1 gene located on chromosome 6B is an early regulator of senescence and affects remobilization of protein and minerals to the grain. GPC-B1 is a NAC transcription factor and has a paralogous copy on chromosome 2B in wheat, GPC-B2. The closest rice homolog to both wheat GPC genes is Os07g37920 which is located on rice chromosome 7 and is colinear with GPC-B2. Since rice is a diploid species with a sequenced genome, we initiated the study of Os07g37920 to develop a simpler model to study senescence and mineral remobilization in cereals. We developed eleven independent RNA interference transgenic rice lines (Os07g37920-RNAi) and 10 over-expressing transgenic lines (Os07g37920-OE), but none of them showed differences in senescence. Transgenic Os07g37920-RNAi rice plants had reduced proportions of viable pollen grains and were male-sterile, but were able to produce seeds by cross pollination. Analysis of the flower morphology of the transgenic rice plants showed that anthers failed to dehisce. Transgenic Os07g37920-OE lines showed no sterility or anther dehiscence problems. Os07g37920 transcript levels were higher in stamens compared to leaves and significantly reduced in the transgenic Os07g37920-RNAi plants. Wheat GPC genes showed the opposite transcription profile (higher transcript levels in leaves than in flowers) and plants carrying knock-out mutations of all GPC-1 and GPC-2 genes exhibited delayed senescence but normal anther dehiscence and fertility. These results indicate a functional divergence of the homologous wheat and rice NAC genes and suggest the need for separate studies of the function and targets of these transcription factors in wheat and rice.

Diversity array technology markers: genetic diversity analyses and linkage map construction in rapeseed (Brassica napus L.)

We developed Diversity Array Technology (DArT) markers for application in genetic studies of Brassica napus and other Brassica species with A or C genomes. Genomic representation from 107 diverse genotypes of B. napus L. var. oleifera (rapeseed, AACC genomes) and B. rapa (AA genome) was used to develop a DArT array comprising 11 520 clones generated using PstI/BanII and PstI/BstN1 complexity reduction methods. In total, 1547 polymorphic DArT markers of high technical quality were identified and used to assess molecular diversity among 89 accessions of B. napus, B. rapa, B. juncea, and B. carinata collected from different parts of the world. Hierarchical cluster and principal component analyses based on genetic distance matrices identified distinct populations clustering mainly according to their origin/pedigrees. DArT markers were also mapped in a new doubled haploid population comprising 131 lines from a cross between spring rapeseed lines ‘Lynx-037DH’ and ‘Monty-028DH’. Linkage groups were assigned on the basis of previously mapped simple sequence repeat (SSRs), intron polymorphism (IP), and gene-based markers. The map consisted of 437 DArT, 135 SSR, 6 IP, and 6 gene-based markers and spanned 2288 cM. Our results demonstrate that DArT markers are suitable for genetic diversity analysis and linkage map construction in rapeseed.

New Diversity Arrays Technology (DArT) markers for tetraploid oat (Avena magna Murphy et Terrell) provide the first complete oat linkage map and markers linked to domestication genes from hexaploid A. sativa L

Nutritional benefits of cultivated oat (Avena sativa L., 2n = 6x = 42, AACCDD) are well recognized; however, seed protein levels are modest and resources for genetic improvement are scarce. The wild tetraploid, A. magna Murphy et Terrell (syn A. maroccana Gdgr., 2n = 4x = 28, CCDD), which contains approximately 31% seed protein, was hybridized with cultivated oat to produce a domesticated A. magna. Wild and cultivated accessions were crossed to generate a recombinant inbred line (RIL) population. Although these materials could be used to develop domesticated, high-protein oat, mapping and quantitative trait loci introgression is hindered by a near absence of genetic markers. Objectives of this study were to develop high-throughput, A. magna-specific markers; generate a genetic linkage map based on the A. magna RIL population; and map genes controlling oat domestication. A Diversity Arrays Technology (DArT) array derived from 10 A. magna genotypes was used to generate 2,688 genome-specific probes. These, with 12,672 additional oat clones, produced 2,349 polymorphic markers, including 498 (21.2%) from A. magna arrays and 1,851 (78.8%) from other Avena libraries. Linkage analysis included 974 DArT markers, 26 microsatellites, 13 SNPs, and 4 phenotypic markers, and resulted in a 14-linkage-group map. Marker-to-marker correlation coefficient analysis allowed classification of shared markers as unique or redundant, and putative linkage-group-to-genome anchoring. Results of this study provide for the first time a collection of high-throughput tetraploid oat markers and a comprehensive map of the genome, providing insights to the genome ancestry of oat and affording a resource for study of oat domestication, gene transfer, and comparative genomics.

Comparisons of four approximation algorithms for large-scale linkage map construction

Efficient construction of large-scale linkage maps is highly desired in current gene mapping projects. To evaluate the performance of available approaches in the literature, four published methods, the insertion (IN), seriation (SER), neighbor mapping (NM), and unidirectional growth (UG) were compared on the basis of simulated F(2) data with various population sizes, interferences, missing genotype rates, and mis-genotyping rates. Simulation results showed that the IN method outperformed, or at least was comparable to, the other three methods. These algorithms were also applied to a real data set and results showed that the linkage order obtained by the IN algorithm was superior to the other methods. Thus, this study suggests that the IN method should be used when constructing large-scale linkage maps.

Stochastic deletion-insertion algorithm to construct dense linkage maps

In this study, we proposed a stochastic deletion-insertion (SDI) algorithm for constructing large-scale linkage maps. This SDI algorithm was compared with three published approximation approaches, the seriation (SER), neighbor mapping (NM), and unidirectional growth (UG) approaches, on the basis of simulated F(2) data with different population sizes, missing genotype rates, and numbers of markers. Simulation results showed that the SDI method had a similar or higher percentage of correct linkage orders than the other three methods. This SDI algorithm was also applied to a real dataset and compared with the other three methods. The total linkage map distance (cM) obtained by the SDI method (148.08 cM) was smaller than the distance obtained by SER (225.52 cM) and two published distances (150.11 cM and 150.38 cM). Since this SDI algorithm is stochastic, a more accurate linkage order can be quickly obtained by repeating this algorithm. Thus, this SDI method, which combines the advantages of accuracy and speed, is an important addition to the current linkage mapping toolkit for constructing improved linkage maps.

LTC: a novel algorithm to improve the efficiency of contig assembly for physical mapping in complex genomes

Background

Physical maps are the substrate of genome sequencing and map-based cloning and their construction relies on the accurate assembly of BAC clones into large contigs that are then anchored to genetic maps with molecular markers. High Information Content Fingerprinting has become the method of choice for large and repetitive genomes such as those of maize, barley, and wheat. However, the high level of repeated DNA present in these genomes requires the application of very stringent criteria to ensure a reliable assembly with the FingerPrinted Contig (FPC) software, which often results in short contig lengths (of 3-5 clones before merging) as well as an unreliable assembly in some difficult regions. Difficulties can originate from a non-linear topological structure of clone overlaps, low power of clone ordering algorithms, and the absence of tools to identify sources of gaps in Minimal Tiling Paths (MTPs).

Results

To address these problems, we propose a novel approach that: (i) reduces the rate of false connections and Q-clones by using a new cutoff calculation method; (ii) obtains reliable clusters robust to the exclusion of single clone or clone overlap; (iii) explores the topological contig structure by considering contigs as networks of clones connected by significant overlaps; (iv) performs iterative clone clustering combined with ordering and order verification using re-sampling methods; and (v) uses global optimization methods for clone ordering and Band Map construction. The elements of this new analytical framework called Linear Topological Contig (LTC) were applied on datasets used previously for the construction of the physical map of wheat chromosome 3B with FPC. The performance of LTC vs. FPC was compared also on the simulated BAC libraries based on the known genome sequences for chromosome 1 of rice and chromosome 1 of maize.

Conclusions

The results show that compared to other methods, LTC enables the construction of highly reliable and longer contigs (5-12 clones before merging), the detection of "weak" connections in contigs and their "repair", and the elongation of contigs obtained by other assembly methods.

Identification and mapping of PmG16, a powdery mildew resistance gene derived from wild emmer wheat

The gene-pool of wild emmer wheat, Triticum turgidum ssp. dicoccoides, harbors a rich allelic repertoire for disease resistance. In the current study, we made use of tetraploid wheat mapping populations derived from a cross between durum wheat (cv. Langdon) and wild emmer (accession G18-16) to identify and map a new powdery mildew resistance gene derived from wild emmer wheat. Initially, the two parental lines were screened with a collection of 42 isolates of Blumeria graminis f. sp. tritici (Bgt) from Israel and 5 isolates from Switzerland. While G18-16 was resistant to 34 isolates, Langdon was resistant only to 5 isolates and susceptible to 42 isolates. Isolate Bgt#15 was selected to differentiate between the disease reactions of the two genotypes. Segregation ratio of F(2-3) and recombinant inbreed line (F(7)) populations to inoculation with isolate Bgt#15 indicated the role of a single dominant gene in conferring resistance to Bgt#15. This gene, temporarily designated PmG16, was located on the distal region of chromosome arm 7AL. Genetic map of PmG16 region was assembled with 32 simple sequence repeat (SSR), sequence tag site (STS), Diversity array technology (DArT) and cleaved amplified polymorphic sequence (CAPS) markers and assigned to the 7AL physical bin map (7AL-16). Using four DNA markers we established colinearity between the genomic region spanning the PmG16 locus within the distal region of chromosome arm 7AL and the genomic regions on rice chromosome 6 and Brachypodium Bd1. A comparative analysis was carried out between PmG16 and other known Pm genes located on chromosome arm 7AL. The identified PmG16 may facilitate the use of wild alleles for improvement of powdery mildew resistance in elite wheat cultivars via marker-assisted selection.

Molecular mapping and candidate gene identification of the Rf2 gene for pollen fertility restoration in sorghum [Sorghum bicolor (L.) Moench]

The A1 cytoplasmic-nuclear male sterility system in sorghum is used almost exclusively for the production of commercial hybrid seed and thus, the dominant genes that restore male fertility in F(1) hybrids are of critical importance to commercial seed production. The genetics of fertility restoration in sorghum can appear complex, being controlled by at least two major genes with additional modifiers and additional gene-environment interaction. To elucidate the molecular processes controlling fertility restoration and to develop a marker screening system for this important trait, two sorghum recombinant inbred line populations were created by crossing a restorer and a non-restoring inbred line, with fertility phenotypes evaluated in hybrid combination with three unique cytoplasmic male sterile lines. In both populations, a single major gene segregated for restoration which was localized to chromosome SBI-02 at approximately 0.5 cM from microsatellite marker, Xtxp304. In the two populations we observed that approximately 85 and 87% of the phenotypic variation in seed set was associated with the major Rf gene on SBI-02. Some evidence for modifier genes was also observed since a continuum of partial restored fertility was exhibited by lines in both RIL populations. With the prior report (Klein et al. in Theor Appl Genet 111:994-1012, 2005) of the cloning of the major fertility restoration gene Rf1 in sorghum, the major fertility restorer locus identified in this study was designated Rf2. A fine-mapping population was used to resolve the Rf2 locus to a 236,219-bp region of chromosome SBI-02, which spanned ~31 predicted open reading frames including a pentatricopeptide repeat (PPR) gene family member. The PPR gene displayed high homology with rice Rf1. Progress towards the development of a marker-assisted screen for fertility restoration is discussed.

Aligning a new reference genetic map of Lupinus angustifolius with the genome sequence of the model legume, Lotus japonicus

We have developed a dense reference genetic map of Lupinus angustifolius (2n = 40) based on a set of 106 publicly available recombinant inbred lines derived from a cross between domesticated and wild parental lines. The map comprised 1090 loci in 20 linkage groups and three small clusters, drawing together data from several previous mapping publications plus almost 200 new markers, of which 63 were gene-based markers. A total of 171 mainly gene-based, sequence-tagged site loci served as bridging points for comparing the Lu. angustifolius genome with the genome sequence of the model legume, Lotus japonicus via BLASTn homology searching. Comparative analysis indicated that the genomes of Lu. angustifolius and Lo. japonicus are highly diverged structurally but with significant regions of conserved synteny including the region of the Lu. angustifolius genome containing the pod-shatter resistance gene, lentus. We discuss the potential of synteny analysis for identifying candidate genes for domestication traits in Lu. angustifolius and in improving our understanding of Fabaceae genome evolution.

A comparison of recombination frequencies in intraspecific versus interspecific mapping populations of Nasonia

We present the first intraspecific linkage map for Nasonia vitripennis based on molecular markers. The map consists of 36 new microsatellite markers, extracted from the Nasonia genome sequence, and spans 515 cM. The five inferred linkage groups correspond to the five chromosomes of Nasonia. Comparison of recombination frequencies of the marker intervals spread over the whole genome (N=33 marker intervals) between the intraspecific N. vitripennis map and an interspecific N. vitripennis x N. giraulti map revealed a slightly higher (1.8%) recombination frequency in the intraspecific cross. We further considered an N. vitripennis x N. longicornis map with 29 microsatellite markers spanning 430 cM. Recombination frequencies in the two interspecific crosses differed neither between reciprocal crosses nor between mapping populations of embryos and adults. No major chromosomal rearrangements were found for the analyzed genomic segments. The observed differential F(2) hybrid male mortality has no significant effect on the genome-wide recombination frequency in Nasonia. We conclude that interspecific crosses between the different Nasonia species, a hallmark of Nasonia genetics, are generally suitable for mapping quantitative and qualitative trait loci for species differences.

Recombination and its impact on the genome of the haplodiploid parasitoid wasp Nasonia

Homologous meiotic recombination occurs in most sexually reproducing organisms, yet its evolutionary advantages are elusive. Previous research explored recombination in the honeybee, a eusocial hymenopteran with an exceptionally high genome-wide recombination rate. A comparable study in a non-social member of the Hymenoptera that would disentangle the impact of sociality from Hymenoptera-specific features such as haplodiploidy on the evolution of the high genome-wide recombination rate in social Hymenoptera is missing. Utilizing single-nucleotide polymorphisms (SNPs) between two Nasonia parasitoid wasp genomes, we developed a SNP genotyping microarray to infer a high-density linkage map for Nasonia. The map comprises 1,255 markers with an average distance of 0.3 cM. The mapped markers enabled us to arrange 265 scaffolds of the Nasonia genome assembly 1.0 on the linkage map, representing 63.6% of the assembled N. vitripennis genome. We estimated a genome-wide recombination rate of 1.4-1.5 cM/Mb for Nasonia, which is less than one tenth of the rate reported for the honeybee. The local recombination rate in Nasonia is positively correlated with the distance to the center of the linkage groups, GC content, and the proportion of simple repeats. In contrast to the honeybee genome, gene density in the parasitoid wasp genome is positively associated with the recombination rate; regions of low recombination are characterized by fewer genes with larger introns and by a greater distance between genes. Finally, we found that genes in regions of the genome with a low recombination frequency tend to have a higher ratio of non-synonymous to synonymous substitutions, likely due to the accumulation of slightly deleterious non-synonymous substitutions. These findings are consistent with the hypothesis that recombination reduces interference between linked sites and thereby facilitates adaptive evolution and the purging of deleterious mutations. Our results imply that the genomes of haplodiploid and of diploid higher eukaryotes do not differ systematically in their recombination rates and associated parameters.

EST-SSR markers derived from an elite barley cultivar (Hordeum vulgare L. 'Morex'): polymorphism and genetic marker potential

Microsatellites or simple sequence repeats have become the markers of choice for marker-assisted selection because of their low template DNA requirement, high reproducibility, and high level of polymorphism. This study investigated a new set of barley (Hordeum vulgare L.) EST-derived SSR markers designed to target gene sequences expressed during grain development, as they are more likely to be important in determining grain quality. The EST sequences (HVSMEh and HVSMEi) were derived from cDNA libraries of the elite six-rowed cultivar Morex, made from spikes harvested at 5 to 45 days after pollination. Approximately half of the 110 SSR markers derived from the ESTs were polymorphic in a panel of 8 diverse barley genotypes, with PIC values between 0.19 and 0.79. Twenty of the new markers were mapped to chromosomal locations using 2 doubled haploid populations. To demonstrate marker potential, quantitative trait locus (QTL) analyses were carried out with phenotypic data on wort beta-glucan content and beta-glucanase activity, two traits with a long history of genetic studies. Most of the EST-SSR markers mapped to within 10 cM of the cellulose synthase (HvCesA) and cellulose synthase-like (HvCslF) genes, which provides highly informative functional markers for tracking these genes in breeding programs. It was also observed that on any given chromosome, the QTL for beta-glucan content and beta-glucanase activity were rarely coincident but tended to occur in adjacent intervals along chromosomal regions, which agreed with their independent genetic basis; the adjacent localization may be important for coordination of cell wall degradation during germination and malting.

Mu transposon insertion sites and meiotic recombination events co-localize with epigenetic marks for open chromatin across the maize genome

The Mu transposon system of maize is highly active, with each of the ∼50–100 copies transposing on average once each generation. The approximately one dozen distinct Mu transposons contain highly similar ∼215 bp terminal inverted repeats (TIRs) and generate 9-bp target site duplications (TSDs) upon insertion. Using a novel genome walking strategy that uses these conserved TIRs as primer binding sites, Mu insertion sites were amplified from Mu stocks and sequenced via 454 technology. 94% of ∼965,000 reads carried Mu TIRs, demonstrating the specificity of this strategy. Among these TIRs, 21 novel Mu TIRs were discovered, revealing additional complexity of the Mu transposon system. The distribution of >40,000 non-redundant Mu insertion sites was strikingly non-uniform, such that rates increased in proportion to distance from the centromere. An identified putative Mu transposase binding consensus site does not explain this non-uniformity. An integrated genetic map containing more than 10,000 genetic markers was constructed and aligned to the sequence of the maize reference genome. Recombination rates (cM/Mb) are also strikingly non-uniform, with rates increasing in proportion to distance from the centromere. Mu insertion site frequencies are strongly correlated with recombination rates. Gene density does not fully explain the chromosomal distribution of Mu insertion and recombination sites, because pronounced preferences for the distal portion of chromosome are still observed even after accounting for gene density. The similarity of the distributions of Mu insertions and meiotic recombination sites suggests that common features, such as chromatin structure, are involved in site selection for both Mu insertion and meiotic recombination. The finding that Mu insertions and meiotic recombination sites both concentrate in genomic regions marked with epigenetic marks of open chromatin provides support for the hypothesis that open chromatin enhances rates of both Mu insertion and meiotic recombination.

Genomic insertion sites of Mu transposons were amplified and sequenced via next generation technology, revealing more than 40,000 non-redundant Mu insertion sites that are non-uniformly distributed across the maize genome and within genes. Along chromosomes, frequencies of Mu transposon insertions are strongly correlated with recombination rates. Although both Mu and recombination occur preferentially in genes, gene density does not fully explain these patterns. Instead, the finding that Mu insertions and meiotic recombination sites both concentrate in genomic regions marked with epigenetic marks of open chromatin provides support for the hypothesis that open chromatin enhances rates of both Mu insertion and meiotic recombination.

High-throughput genetic mapping of mutants via quantitative single nucleotide polymorphism typing

Advances in next-generation sequencing technology have facilitated the discovery of single nucleotide polymorphisms (SNPs). Sequenom-based SNP-typing assays were developed for 1359 maize SNPs identified via comparative next-generation transcriptomic sequencing. Approximately 75% of these SNPs were successfully converted into genetic markers that can be scored reliably and used to generate a SNP-based genetic map by genotyping recombinant inbred lines from the intermated B73 x Mo17 population. The quantitative nature of Sequenom-based SNP assays led to the development of a time- and cost-efficient strategy to genetically map mutants via quantitative bulked segregant analysis. This strategy was used to rapidly map the loci associated with several dozen recessive mutants. Because a mutant can be mapped using as few as eight multiplexed sets of SNP assays on a bulk of as few as 20 mutant F(2) individuals, this strategy is expected to be widely adopted for mapping in many species.

Evaluation of algorithms used to order markers on genetic maps

When building genetic maps, it is necessary to choose from several marker ordering algorithms and criteria, and the choice is not always simple. In this study, we evaluate the efficiency of algorithms try (TRY), seriation (SER), rapid chain delineation (RCD), recombination counting and ordering (RECORD) and unidirectional growth (UG), as well as the criteria PARF (product of adjacent recombination fractions), SARF (sum of adjacent recombination fractions), SALOD (sum of adjacent LOD scores) and LHMC (likelihood through hidden Markov chains), used with the RIPPLE algorithm for error verification, in the construction of genetic linkage maps. A linkage map of a hypothetical diploid and monoecious plant species was simulated containing one linkage group and 21 markers with fixed distance of 3 cM between them. In all, 700 F(2) populations were randomly simulated with 100 and 400 individuals with different combinations of dominant and co-dominant markers, as well as 10 and 20% of missing data. The simulations showed that, in the presence of co-dominant markers only, any combination of algorithm and criteria may be used, even for a reduced population size. In the case of a smaller proportion of dominant markers, any of the algorithms and criteria (except SALOD) investigated may be used. In the presence of high proportions of dominant markers and smaller samples (around 100), the probability of repulsion linkage increases between them and, in this case, use of the algorithms TRY and SER associated to RIPPLE with criterion LHMC would provide better results.

A consensus genetic map of sorghum that integrates multiple component maps and high-throughput Diversity Array Technology (DArT) markers

BACKGROUND

Sorghum genome mapping based on DNA markers began in the early 1990s and numerous genetic linkage maps of sorghum have been published in the last decade, based initially on RFLP markers with more recent maps including AFLPs and SSRs and very recently, Diversity Array Technology (DArT) markers. It is essential to integrate the rapidly growing body of genetic linkage data produced through DArT with the multiple genetic linkage maps for sorghum generated through other marker technologies. Here, we report on the colinearity of six independent sorghum component maps and on the integration of these component maps into a single reference resource that contains commonly utilized SSRs, AFLPs, and high-throughput DArT markers.

RESULTS

The six component maps were constructed using the MultiPoint software. The lengths of the resulting maps varied between 910 and 1528 cM. The order of the 498 markers that segregated in more than one population was highly consistent between the six individual mapping data sets. The framework consensus map was constructed using a "Neighbours" approach and contained 251 integrated bridge markers on the 10 sorghum chromosomes spanning 1355.4 cM with an average density of one marker every 5.4 cM, and were used for the projection of the remaining markers. In total, the sorghum consensus map consisted of a total of 1997 markers mapped to 2029 unique loci (1190 DArT loci and 839 other loci) spanning 1603.5 cM and with an average marker density of 1 marker/0.79 cM. In addition, 35 multicopy markers were identified. On average, each chromosome on the consensus map contained 203 markers of which 58.6% were DArT markers. Non-random patterns of DNA marker distribution were observed, with some clear marker-dense regions and some marker-rare regions.

CONCLUSION

The final consensus map has allowed us to map a larger number of markers than possible in any individual map, to obtain a more complete coverage of the sorghum genome and to fill a number of gaps on individual maps. In addition to overall general consistency of marker order across individual component maps, good agreement in overall distances between common marker pairs across the component maps used in this study was determined, using a difference ratio calculation. The obtained consensus map can be used as a reference resource for genetic studies in different genetic backgrounds, in addition to providing a framework for transferring genetic information between different marker technologies and for integrating DArT markers with other genomic resources. DArT markers represent an affordable, high throughput marker system with great utility in molecular breeding programs, especially in crops such as sorghum where SNP arrays are not publicly available.

Genetics of phenotypic plasticity: QTL analysis in barley, Hordeum vulgare

Phenotypic plasticity is the variation in phenotypic traits produced by a genotype in different environments. In contrast, environmental canalization is defined as the insensitivity of a genotype's phenotype to variation in environments. Despite the extensive literature on the evolutionary significance and potential genetic mechanisms driving plasticity and canalization, few studies tried to unravel the genetic basis of this phenomenon. Using both simulations and real data from barley (Hordeum vulgare), we used QTL mapping to obtain insights into the genetics of phenotypic plasticity. We explored two ways of quantifying phenotypic plasticity, namely the phenotypic variance across environments and the Finlay-Wilkinson's regression slope. Each relates to a different concept of stability. Through QTL detection with real and simulated data, we show that each measure of plasticity detects specific types of plasticity QTL. Most of the plasticity QTLs were detected in the data set with the lowest number of environments. All plasticity QTL co-located with loci showing QTL x E interaction and there were no QTL that only affected plasticity. The number of environments that are considered and their homogeneity is a key to interpret the genetic control of phenotypic plasticity. Regulatory pathways of plasticity may vary from one set of environments to another due to unique features of each environment. Therefore, with an increasing number of environments, it may become impossible to detect a single 'consistent' regulatory pathway for all environments.

QTL analysis of ergot resistance in sorghum

Sorghum ergot, caused predominantly by Claviceps africana Frederickson, Mantle, de Milliano, is a significant threat to the sorghum industry worldwide. The objectives of this study were firstly, to identify molecular markers linked to ergot resistance and to two pollen traits, pollen quantity (PQ) and pollen viability (PV), and secondly, to assess the relationship between the two pollen traits and ergot resistance in sorghum. A genetic linkage map of sorghum RIL population R931945-2-2 x IS 8525 (resistance source) was constructed using 303 markers including 36 SSR, 117 AFLP , 148 DArT and two morphological trait loci. Composite interval mapping identified nine, five, and four QTL linked to molecular markers for percentage ergot infection (PCERGOT), PQ and PV, respectively, at a LOD >2.0. Co-location/linkage of QTL were identified on four chromosomes while other QTL for the three traits mapped independently, indicating that both pollen and non pollen-based mechanisms of ergot resistance were operating in this sorghum population. Of the nine QTL identified for PCERGOT, five were identified using the overall data set while four were specific to the group data sets defined by temperature and humidity. QTL identified on SBI-02 and SBI-06 were further validated in additional populations. This is the first report of QTL associated with ergot resistance in sorghum. The markers reported herein could be used for marker-assisted selection for this important disease of sorghum.

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 and QTL analyses of seed dormancy and preharvest sprouting resistance in the wheat germplasm CN10955

The inheritance and genetic linkage analysis for seed dormancy and preharvest sprouting (PHS) resistance were carried out in an F8 recombinant inbred lines (RILs) derived from the cross between "CN19055" (white-grained, PHS-resistant) with locally adapted Australian cultivar "Annuello" (white-grained, PHS-susceptible). Seed dormancy was assessed as germination index (GI7) while assessment for preharvest sprouting resistance was based on whole head assay (sprouting index, SI) and visibly sprouted seeds (VI). Segregation analysis of the F2, F3 data from the glasshouse and the RIL population in 2004 and 2005 field data sets indicated that seed dormancy and PHS resistance in CN19055 is controlled by at least two genes. Heritabilities for GI7 and VI were high and moderate for SI. The most accurate method for assessing PHS resistance was achieved using VI and GI7 while SI exhibited large genotype by environment interaction. Two quantitative trait loci (QTLs) QPhs.dpivic.4A.1 and QPhs.dpivic.4A.2 were identified. On pooled data across four environments, the major QTL, QPhs.dpivic.4A.2, explained 45% of phenotypic variation for GI7, 43% for VI and 20% for SI, respectively. On the other hand, QPhs.dpivic.4A.1 which accounted for 31% of the phenotypic variation in GI7 in 2004 Horsham field trial, was not stable across environments. Physical mapping of two SSR markers, Xgwm937 and Xgwm894 linked to the major QTL for PHS resistance, using Chinese Spring deletions lines for chromosome 4AS and 4AL revealed that the markers were located in the deletion bins 4AL-12 and 4AL-13. The newly identified SSR markers (Xgwm937/Xgwm894) showed strong association with seed dormancy and PHS resistance in a range of wheat lines reputed to possess PHS resistance. The results suggest that Xgwm937/Xgwm894 could be used in marker-assisted selection (MAS) for incorporating preharvest sprouting resistance into elite wheat cultivars susceptible to PHS.