TriMEDB: A database to integrate transcribed markers and facilitate genetic studies of the tribe Triticeae

Background  

The recent rapid accumulation of sequence resources of various crop species ensures an improvement in the genetics approach, including quantitative trait loci (QTL) analysis as well as the holistic population analysis and association mapping of natural variations. Because the tribe Triticeae includes important cereals such as wheat and barley, integration of information on the genetic markers in these crops should effectively accelerate map-based genetic studies on Triticeae species and lead to the discovery of key loci involved in plant productivity, which can contribute to sustainable food production. Therefore, informatics applications and a semantic knowledgebase of genome-wide markers are required for the integration of information on and further development of genetic markers in wheat and barley in order to advance conventional marker-assisted genetic analyses and population genomics of Triticeae species.     

Simple sequence repeat markers useful for sorghum downy mildew (Peronosclerospora sorghi) and related species

Background

Association mapping, initially developed in human disease genetics, is now being applied to plant species. The model species Arabidopsis provided some of the first examples of association mapping in plants, identifying previously cloned flowering time genes, despite high population sub-structure. More recently, association genetics has been applied to barley, where breeding activity has resulted in a high degree of population sub-structure. A major genotypic division within barley is that between winter- and spring-sown varieties, which differ in their requirement for vernalization to promote subsequent flowering. To date, all attempts to validate association genetics in barley by identifying major flowering time loci that control vernalization requirement (VRN-H1 and VRN-H2) have failed. Here, we validate the use of association genetics in barley by identifying VRN-H1 and VRN-H2, despite their prominent role in determining population sub-structure.

Results

By taking barley as a typical inbreeding crop, and seasonal growth habit as a major partitioning phenotype, we develop an association mapping approach which successfully identifies VRN-H1 and VRN-H2, the underlying loci largely responsible for this agronomic division. We find a combination of Structured Association followed by Genomic Control to correct for population structure and inflation of the test statistic, resolved significant associations only with VRN-H1 and the VRN-H2 candidate genes, as well as two genes closely linked to VRN-H1 (HvCSFs1 and HvPHYC).

Conclusion

We show that, after employing appropriate statistical methods to correct for population sub-structure, the genome-wide partitioning effect of allelic status at VRN-H1 and VRN-H2 does not result in the high levels of spurious association expected to occur in highly structured samples. Furthermore, we demonstrate that both VRN-H1 and the candidate VRN-H2 genes can be identified using association mapping. Discrimination between intragenic VRN-H1 markers was achieved, indicating that candidate causative polymorphisms may be discerned and prioritised within a larger set of positive associations. This proof of concept study demonstrates the feasibility of association mapping in barley, even within highly structured populations. A major advantage of this method is that it does not require large numbers of genome-wide markers, and is therefore suitable for fine mapping and candidate gene evaluation, especially in species for which large numbers of genetic markers are either unavailable or too costly.     

A new set of BXD recombinant inbred lines from advanced intercross populations in mice

Background  

Recombinant inbred (RI) strains are an important resource for mapping complex traits in many species. While large RI panels are available for Arabidopsis, maize, C. elegans, and Drosophila, mouse RI panels typically consist of fewer than 30 lines. This is a severe constraint on the power and precision of mapping efforts and greatly hampers analysis of epistatic interactions.     

Association mapping of resistance to Puccinia hordei in Australian barley breeding germplasm

Key message

“To find stable resistance using association mapping tools, QTL with major and minor effects on leaf rust reactions were identified in barley breeding lines by assessing seedlings and adult plants.”

Abstract

Three hundred and sixty (360) elite barley (Hordeum vulgare L.) breeding lines from the Northern Region Barley Breeding Program in Australia were genotyped with 3,244 polymorphic diversity arrays technology markers and the results used to map quantitative trait loci (QTL) conferring a reaction to leaf rust (Puccinia hordei Otth). The F_3:5 (Stage 2) lines were derived or sourced from different geographic origins or hubs of international barley breeding ventures representing two breeding cycles (2009 and 2011 trials) and were evaluated across eight environments for infection type at both seedling and adult plant stages. Association mapping was performed using mean scores for disease reaction, accounting for family effects using the eigenvalues from a matrix of genotype correlations. In this study, 15 QTL were detected; 5 QTL co-located with catalogued leaf rust resistance genes (Rph1, Rph3/19, Rph8/14/15, Rph20, Rph21), 6 QTL aligned with previously reported genomic regions and 4 QTL (3 on chromosome 1H and 1 on 7H) were novel. The adult plant resistance gene Rph20 was identified across the majority of environments and pathotypes. The QTL detected in this study offer opportunities for breeding for more durable resistance to leaf rust through pyramiding multiple genomic regions via marker-assisted selection.     

A first step toward the development of a barley NAM population and its utilization to detect QTLs conferring leaf rust seedling resistance

Key message

We suggest multi-parental nested association mapping as a valuable innovation in barley genetics, which increases the power to map quantitative trait loci and assists in extending genetic diversity of the elite barley gene pool.

Abstract

Plant genetic resources are a key asset to further improve crop species. The nested association mapping (NAM) approach was introduced to identify favorable genes in multi-parental populations. Here, we report toward the development of the first explorative barley NAM population and demonstrate its usefulness in a study on mapping quantitative trait loci (QTLs) for leaf rust resistance. The NAM population HEB-5 was developed from crossing and backcrossing five exotic barley donors with the elite barley cultivar ‘Barke,’ resulting in 295 NAM lines in generation BC₁S₁. HEB-5 was genetically characterized with 1,536 barley SNPs. Across HEB-5 and within the NAM families, no deviation from the expected genotype and allele frequencies was detected. Genetic similarity between ‘Barke’ and the NAM families ranged from 78.6 to 83.1 %, confirming the backcrossing step during population development. To explore its usefulness, a screen for leaf rust (Puccinia hordei) seedling resistance was conducted. Resistance QTLs were mapped to six barley chromosomes, applying a mixed model genome-wide association study. In total, four leaf rust QTLs were detected across HEB-5 and four QTLs within family HEB-F23. Favorable exotic QTL alleles reduced leaf rust symptoms on two chromosomes by 33.3 and 36.2 %, respectively. The located QTLs may represent new resistance loci or correspond to new alleles of known resistance genes. We conclude that the exploratory population HEB-5 can be applied to mapping and utilizing exotic QTL alleles of agronomic importance. The NAM concept will foster the evaluation of the genetic diversity, which is present in our primary barley gene pool.     

Comparative analysis of haplotype association mapping algorithms

Background  

Finding the genetic causes of quantitative traits is a complex and difficult task. Classical methods for mapping quantitative trail loci (QTL) in miceuse an F2 cross between two strains with substantially different phenotype and an interval mapping method to compute confidence intervals at each position in the genome. This process requires significant resources for breeding and genotyping, and the data generated are usually only applicable to one phenotype of interest. Recently, we reported the application of a haplotype association mapping method which utilizes dense genotyping data across a diverse panel of inbred mouse strains and a marker association algorithm that is independent of any specific phenotype. As the availability of genotyping data grows in size and density, analysis of these haplotype association mapping methods should be of increasing value to the statistical genetics community.     

Plant traits correlated with generation time directly affect inbreeding depression and mating system and indirectly genetic structure

Background  

Understanding the mechanisms that control species genetic structure has always been a major objective in evolutionary studies. The association between genetic structure and species attributes has received special attention. As species attributes are highly taxonomically constrained, phylogenetically controlled methods are necessary to infer causal relationships. In plants, a previous study controlling for phylogenetic signal has demonstrated that Wright's F _ST, a measure of genetic differentiation among populations, is best predicted by the mating system (outcrossing, mixed-mating or selfing) and that plant traits such as perenniality and growth form have only an indirect influence on F _ST via their association with the mating system. The objective of this study is to further outline the determinants of plant genetic structure by distinguishing the effects of mating system on gene flow and on genetic drift. The association of biparental inbreeding and inbreeding depression with population genetic structure, mating system and plant traits are also investigated.     

Multipoint-likelihood maximization mapping on 4 segregating populations to achieve an integrated framework map for QTL analysis in pot azalea (Rhododendron simsii hybrids)

Background  

Azalea (Rhododendron simsii hybrids) is the most important flowering pot plant produced in Belgium, being exported world-wide. In the breeding program, flower color is the main feature for selection, only in later stages cultivation related plant quality traits are evaluated. As a result, plants with attractive flowering are kept too long in the breeding cycle. The inheritance of flower color has been well studied; information on the heritability of cultivation related quality traits is lacking. For this purpose, QTL mapping in diverse genetic backgrounds appeared to be a must and therefore 4 mapping populations were made and analyzed.     

Modelling the genetic architecture of flowering time control in barley through nested association mapping

Background

Barley, globally the fourth most important cereal, provides food and beverages for humans and feed for animal husbandry. Maximizing grain yield under varying climate conditions largely depends on the optimal timing of flowering. Therefore, regulation of flowering time is of extraordinary importance to meet future food and feed demands. We developed the first barley nested association mapping (NAM) population, HEB-25, by crossing 25 wild barleys with one elite barley cultivar, and used it to dissect the genetic architecture of flowering time.

Results

Upon cultivation of 1,420 lines in multi-field trials and applying a genome-wide association study, eight major quantitative trait loci (QTL) were identified as main determinants to control flowering time in barley. These QTL accounted for 64% of the cross-validated proportion of explained genotypic variance (p_G). The strongest single QTL effect corresponded to the known photoperiod response gene Ppd-H1. After sequencing the causative part of Ppd-H1, we differentiated twelve haplotypes in HEB-25, whereof the strongest exotic haplotype accelerated flowering time by 11 days compared to the elite barley haplotype. Applying a whole genome prediction model including main effects and epistatic interactions allowed predicting flowering time with an unmatched accuracy of 77% of cross-validated p_G.

Conclusions

The elaborated causal models represent a fundamental step to explain flowering time in barley. In addition, our study confirms that the exotic biodiversity present in HEB-25 is a valuable toolbox to dissect the genetic architecture of important agronomic traits and to replenish the elite barley breeding pool with favorable, trait-improving exotic alleles.

     

Analysis of molecular diversity,population structure and linkage disequilibrium in a worldwide survey of cultivated barley germplasm (<Emphasis Type="Italic">Hordeum vulgare</Emphasis> L.)

Background

The goal of our study was a systematic survey of the molecular diversity in barley genetic resources. To this end 953 cultivated barley accessions originating from all inhabited continents except Australia were genotyped with 48 SSR markers. Molecular diversity was evaluated with routine statistics (allelic richness, gene diversity, allele frequency, heterozygosity and unique alleles), Principal Coordinate Analysis (PCoA), and analysis of genome-wide linkage disequilibrium.

Results

A genotyping database for 953 cultivated barley accessions profiled with 48 SSR markers was established. The PCoA revealed structuring of the barley population with regard to (i) geographical regions and (ii) agronomic traits. Geographic origin contributed most to the observed molecular diversity. Genome-wide linkage disequilibrium (LD) was estimated as squared correlation of allele frequencies (r²). The values of LD for barley were comparable to other plant species (conifers, poplar, maize). The pattern of intrachromosomal LD with distances between the genomic loci ranging from 1 to 150 cM revealed that in barley LD extended up to distances as long as 50 cM with r² > 0.05, or up to 10 cM with r² > 0.2. Few loci mapping to different chromosomes showed significant LD with r² > 0.05. The number of loci in significant LD as well as the pattern of LD were clearly dependent on the population structure. The LD in the homogenous group of 207 European 2-rowed spring barleys compared to the highly structured worldwide barley population was increased in the number of loci pairs with r² > 0.05 and had higher values of r², although the percentage of intrachromosomal loci pairs in significant LD based on P < 0.001 was 100% in the whole set of varieties, but only 45% in the subgroup of European 2-rowed spring barleys. The value of LD also varied depending on the polymorphism of the loci selected for genotyping. The 17 most polymorphic loci (PIC > 0.80) provided higher LD values as compared to 19 low polymorphic loci (PIC < 0.73) in both structured (all accessions) and non-structured (European 2-rowed spring varieties) barley populations.

Conclusion

A global population of cultivated barley accessions was highly structured. Clustering highlighted the accessions with the same geographic origin, as well as accessions possessing similar agronomic characters. LD in barley extended up to 50 cM, and was strongly dependent on the population structure. The data on LD were summarized as a genome-wide LD map for barley.

     

Adaptation and diversity along an altitudinal gradient in Ethiopian barley (<Emphasis Type="Italic">Hordeum vulgare</Emphasis> L.) landraces revealed by molecular analysis

Background  

Among the cereal crops, barley is the species with the greatest adaptability to a wide range of environments. To determine the level and structure of genetic diversity in barley (Hordeum vulgare L.) landraces from the central highlands of Ethiopia, we have examined the molecular variation at seven nuclear microsatellite loci.     

Statistical epistasis between candidate gene alleles for complex tuber traits in an association mapping population of tetraploid potato

Association mapping using DNA-based markers is a novel tool in plant genetics for the analysis of complex traits. Potato tuber yield, starch content, starch yield and chip color are complex traits of agronomic relevance, for which carbohydrate metabolism plays an important role. At the functional level, the genes and biochemical pathways involved in carbohydrate metabolism are among the best studied in plants. Quantitative traits such as tuber starch and sugar content are therefore models for association genetics in potato based on candidate genes. In an association mapping experiment conducted with a population of 243 tetraploid potato varieties and breeding clones, we previously identified associations between individual candidate gene alleles and tuber starch content, starch yield and chip quality. In the present paper, we tested 190 DNA markers at 36 loci scored in the same association mapping population for pairwise statistical epistatic interactions. Fifty marker pairs were associated mainly with tuber starch content and/or starch yield, at a cut-off value of q ≤ 0.20 for the experiment-wide false discovery rate (FDR). Thirteen marker pairs had an FDR of q ≤ 0.10. Alleles at loci encoding ribulose-bisphosphate carboxylase/oxygenase activase (Rca), sucrose phosphate synthase (Sps) and vacuolar invertase (Pain1) were most frequently involved in statistical epistatic interactions. The largest effect on tuber starch content and starch yield was observed for the paired alleles Pain1-8c and Rca-1a, explaining 9 and 10% of the total variance, respectively. The combination of these two alleles increased the means of tuber starch content and starch yield. Biological models to explain the observed statistical epistatic interactions are discussed.     

Linkage Disequilibrium Mapping of Morphological, Resistance, and Other Agronomically Relevant Traits in Modern Spring Barley Cultivars

A set of 148 modern spring barley cultivars was explored for the extent of linkage disequilibrium (LD) between genes governing traits and nearby marker alleles. Associations of agronomically relevant traits (days to heading, plant height), resistance traits (leaf rust, barley yellow dwarf virus (BYD)), and morphological traits (rachilla hair length, lodicule size) with AFLP markers and SSR markers were found. Known major genes and QTLs were confirmed, but also new putative QTLs were found. The LD mapping clearly indicated the common occurrence of Rph3, a gene for hypersensitivity resistance against Puccinia hordei, and also confirmed the QTL Rphq2 for prolonging latency period of P. hordei in seedlings. We also found strong indication for a hitherto not reported gene for resistance or tolerance to BYD on chromosome 2, linked to SSR marker HVM054. Our conclusion is that LD mapping is a valuable additional tool in the search for applicable marker associations with major genes and QTLs. Electronic supplementary material Electronic supplementary material is available for this article at and accessible for authorised users.     

Resistance to <Emphasis Type="Italic">Rhynchosporium commune</Emphasis> in a collection of European spring barley germplasm

Key message

Association analyses of resistance to Rhynchosporium commune in a collection of European spring barley germplasm detected 17 significant resistance quantitative trait loci. The most significant association was confirmed as Rrs1.

Abstract

Rhynchosporium commune is a fungal pathogen of barley which causes a highly destructive and economically important disease known as rhynchosporium. Genome-wide association mapping was used to investigate the genetic control of host resistance to R. commune in a collection of predominantly European spring barley accessions. Multi-year disease nursery field trials revealed 8 significant resistance quantitative trait loci (QTL), whilst a separate association mapping analysis using historical data from UK national and recommended list trials identified 9 significant associations. The most significant association identified in both current and historical data sources, collocated with the known position of the major resistance gene Rrs1. Seedling assays with R. commune single-spore isolates expressing the corresponding avirulence protein NIP1 confirmed that this locus is Rrs1. These results highlight the significant and continuing contribution of Rrs1 to host resistance in current elite spring barley germplasm. Varietal height was shown to be negatively correlated with disease severity, and a resistance QTL was identified that co-localised with the semi-dwarfing gene sdw1, previously shown to contribute to disease escape. The remaining QTL represent novel resistances that are present within European spring barley accessions. Associated markers to Rrs1 and other resistance loci, identified in this study, represent a set of tools that can be exploited by breeders for the sustainable deployment of varietal resistance in new cultivars.

     

Association mapping utilizing diverse barley lines reveals net form net blotch seedling resistance/susceptibility loci

Key message

A diverse collection of barley lines was phenotyped with three North American Pyrenophora teres f. teres isolates and association analyses detected 78 significant marker-trait associations at 16 genomic loci.

Abstract

Pyrenophora teres f. teres is a necrotrophic fungal pathogen and the causal agent of the economically important foliar disease net form net blotch (NFNB) of barley. The deployment of effective and durable resistance against P. teres f. teres has been hindered by the complexity of quantitative resistance and susceptibility. Several bi-parental mapping populations have been used to identify QTL associated with NFNB disease on all seven barley chromosomes. Here, we report the first genome-wide association study (GWAS) to detect marker-trait associations for resistance or susceptibility to P. teres f. teres. Geographically diverse barley genotypes from a world barley core collection (957) were genotyped with the Illumina barley iSelect chip and phenotyped with three P. teres f. teres isolates collected in two geographical regions of the USA (15A, 6A and LDNH04Ptt19). The best of nine regression models tested were identified for each isolate and used for association analysis resulting in the identification of 78 significant marker-trait associations (MTA; ?log10p value?>3.0). The MTA identified corresponded to 16 unique genomic loci as determined by analysis of local linkage disequilibrium between markers that did not meet a correlation threshold of R ²?≥?0.1, indicating that the markers represented distinct loci. Five loci identified represent novel QTL and were designated QRptts-3HL, QRptts-4HS, QRptts-5HL.1, QRptts-5HL.2, and QRptts-7HL.1. In addition, 55 of the barley lines examined exhibited a high level of resistance to all three isolates and the SNP markers identified will provide useful genetic resources for barley breeding programs.

     

Association mapping of drought tolerance-related traits in barley to complement a traditional biparental QTL mapping study

Key message

Association mapping of drought-related traits in barley was used to increase the density of existing QTL maps without recreating mapping populations.

Abstract

We used 109 spring barley genotypes exhibiting high or low drought tolerance to elucidate the associations between diversity array technology sequencing (DArTseq) and single nucleotide polymorphism (SNP) markers and various physiological parameters related to plant responses to drought conditions. The study was performed in controlled conditions (growth chambers), drought tolerance was phenotyped in the four-leaf seedlings. We identified 58 associations including 34 new markers (i.e., 16 DArTseq and 18 SNP markers). The results for three markers were consistent with the data obtained in an earlier traditional biparental QTL mapping study. The regions neighboring markers on linkage group 2H contained the highest number of significant marker–trait associations. Five markers related to the photosynthetic activity of photosystem II were detected on chromosome 4H. The lowest number of associations were observed for the sequences neighboring DArT markers on linkage group 6H. A chromosome 3H region related to water use efficiency and net photosynthesis rate in both biparental QTL, and association study, may be particularly valuable, as these parameters correspond to the ability of plants to remain highly productive under water deficit stress. Our findings confirm that association mapping can increase the density of existing QTL maps without recreating mapping populations.

     

Association of common polymorphisms in known susceptibility genes with rheumatoid arthritis in a Slovak population using osteoarthritis patients as controls

Introduction  

Both genetic and environmental factors contribute to rheumatoid arthritis (RA), a common and complex autoimmune disease. As well as the major susceptibility gene HLA-DRB1, recent genome-wide and candidate-gene studies reported additional evidence for association of single nucleotide polymorphism (SNP) markers in the PTPN22, STAT4, OLIG3/TNFAIP3 and TRAF1/C5 loci with RA. This study was initiated to investigate the association between defined genetic markers and RA in a Slovak population. In contrast to recent studies, we included intensively-characterized osteoarthritis (OA) patients as controls.     

Retrospective genomic analysis of sorghum adaptation to temperate-zone grain production

Background

Sorghum is a tropical C₄ cereal that recently adapted to temperate latitudes and mechanized grain harvest through selection for dwarfism and photoperiod-insensitivity. Quantitative trait loci for these traits have been introgressed from a dwarf temperate donor into hundreds of diverse sorghum landraces to yield the Sorghum Conversion lines. Here, we report the first comprehensive genomic analysis of the molecular changes underlying this adaptation.

Results

We apply genotyping-by-sequencing to 1,160 Sorghum Conversion lines and their exotic progenitors, and map donor introgressions in each Sorghum Conversion line. Many Sorghum Conversion lines carry unexpected haplotypes not found in either presumed parent. Genome-wide mapping of introgression frequencies reveals three genomic regions necessary for temperate adaptation across all Sorghum Conversion lines, containing the Dw1, Dw2, and Dw3 loci on chromosomes 9, 6, and 7 respectively. Association mapping of plant height and flowering time in Sorghum Conversion lines detects significant associations in the Dw1 but not the Dw2 or Dw3 regions. Subpopulation-specific introgression mapping suggests that chromosome 6 contains at least four loci required for temperate adaptation in different sorghum genetic backgrounds. The Dw1 region fractionates into separate quantitative trait loci for plant height and flowering time.

Conclusions

Generating Sorghum Conversion lines has been accompanied by substantial unintended gene flow. Sorghum adaptation to temperate-zone grain production involves a small number of genomic regions, each containing multiple linked loci for plant height and flowering time. Further characterization of these loci will accelerate the adaptation of sorghum and related grasses to new production systems for food and fuel.     

Mapping gene expression quantitative trait loci by singular value decomposition and independent component analysis

Background  

The combination of gene expression profiling with linkage analysis has become a powerful paradigm for mapping gene expression quantitative trait loci (eQTL). To date, most studies have searched for eQTL by analyzing gene expression traits one at a time. As thousands of expression traits are typically analyzed, this can reduce power because of the need to correct for the number of hypothesis tests performed. In addition, gene expression traits exhibit a complex correlation structure, which is ignored when analyzing traits individually.