Linkage disequilibrium in European elite maize germplasm investigated with SSRs

Information about the extent and genomic distribution of linkage disequilibrium (LD) is of fundamental importance for association mapping. The main objectives of this study were to (1) investigate genetic diversity within germplasm groups of elite European maize (Zea mays L.) inbred lines, (2) examine the population structure of elite European maize germplasm, and (3) determine the extent and genomic distribution of LD between pairs of simple sequence repeat (SSR) markers. We examined genetic diversity and LD in a cross section of European and US elite breeding material comprising 147 inbred lines genotyped with 100 SSR markers. For gene diversity within each group, significant (P<0.05) differences existed among the groups. The LD was significant (P<0.05) for 49% of the SSR marker pairs in the 80 flint lines and for 56% of the SSR marker pairs in the 57 dent lines. The ratio of linked to unlinked loci in LD was 1.1 for both germplasm groups. The high incidence of LD suggests that the extent of LD between SSR markers should allow the detection of marker-phenotype associations in a genome scan. However, our results also indicate that a high proportion of the observed LD is generated by forces, such as relatedness, population stratification, and genetic drift, which cause a high risk of detecting false positives in association mapping.     

Recombination disequilibrium in ideal and natural populations

Following Hardy-Weinberg disequilibrium (HWD) occurring at a single locus and linkage disequilibrium (LD) between two loci in generations, we here proposed the third genetic disequilibrium in a population: recombination disequilibrium (RD). RD is a measurement of crossover interference among multiple loci in a random mating population. In natural populations besides recombination interference, RD may also be due to selection, mutation, gene conversion, drift and/or migration. Therefore, similarly to LD, RD will also reflect the history of natural selection and mutation. In breeding populations, RD purely results from recombination interference and hence can be used to build or evaluate and correct a linkage map. Practical examples from F₂, testcross and human populations indeed demonstrate that RD is useful for measuring recombination interference between two short intervals and evaluating linkage maps. As with LD, RD will be important for studying genetic mapping, association of haplotypes with disease, plant breading and population history.     

Comparison of Linkage Disequilibrium in Elite European Maize Inbred Lines using AFLP and SSR Markers

Application of association mapping to plant breeding populations has the potential to revolutionize plant genetics. The main objectives of this study were to (i) investigate the extent and genomic distribution of linkage disequilibrium (LD) between pairs of amplified fragment length polymorphism (AFLP) markers, (ii) compare these results with those obtained with simple sequence repeat (SSR) markers, and (iii) compare the usefulness of AFLP and SSR markers for genomewide association mapping in plant breeding populations. We examined LD in a cross-section of 72 European elite inbred lines genotyped with 452 AFLP and 93 SSR markers. LD was significant (p < 0.05) for about 15% of the AFLP marker pairs and for about 49% of the SSR marker pairs in each of the two germplasm groups, flint and dent. In both germplasm groups the ratio of linked to unlinked loci pairs in LD was higher for AFLPs than for SSRs. The observation of LD due to linkage for both marker types suggested that genome-wide association mapping should be possible using either AFLPs or SSRs. The results of our study indicated that SSRs should be favored over AFLPs but the opposite applies to populations with a long history of recombination.     

Linkage disequilibrium in two European F<Subscript>2</Subscript> flint maize populations under modified recurrent full-sib selection

According to quantitative genetic theory, linkage disequilibrium (LD) can hamper the short- and long-term selection response in recurrent selection (RS) programs. We analyzed LD in two European flint maize populations, KW1265 x D146 (A x B) and D145 x KW1292 (C x D), under modified recurrent full-sib selection. Our objectives were to investigate (1) the decay of initial parental LD present in F(2) populations by three generations of intermating, (2) the generation of new LD in four (A x B) and seven (C x D) selection cycles, and (3) the relationship between LD changes and estimates of the additive genetic variance. We analyzed the F(2) and the intermated populations as well as all selection cycles with 104 (A x B) and 101 (C x D) simple sequence repeat (SSR) markers with a uniform coverage of the entire maize genome. The LD coefficient D and the composite LD measure Delta were estimated and significance tests for LD were performed. LD was reduced by intermating as expected from theory. A directional generation of negative LD between favorable alleles could not be observed during the selection cycles. However, considerable undirectional changes in D were observed, which we attributed to genetic sampling due to the finite population size used for recombination. Consequently, a long-term reduction of the additive genetic variance due to negative LD was not observed. Our experimental results support the hypothesis that in practical RS programs with maize, LD generated by selection is not a limiting factor for obtaining a high selection response.     

A new test for family-based association mapping with inbred lines from plant breeding programs

Association mapping holds great promise for the detection of quantitative trait loci (QTL) in plant breeding populations. The main objectives of this study were to (1) adapt the quantitative pedigree disequilibrium test to typical pedigrees of inbred lines produced in plant breeding programs, (2) compare the newly developed quantitative inbred pedigree disequilibrium test (QIPDT) with the commonly employed logistic regression ratio test (LRRT), with respect to the power and type I error rate of QTL detection, and (3) demonstrate the use of the QIPDT by applying it to flowering data of European elite maize inbreds. QIPDT and LRRT were compared based on computer simulations modeling 55 years of hybrid maize breeding in Central Europe. Furthermore, we applied QIPDT to a cross-section of 49 European elite maize inbred lines genotyped with 722 amplified fragment length polymorphism markers and phenotyped in four environments for days to anthesis. Compared to LRRT, the power to detect QTL was higher with QIPDT when using data collected routinely in plant breeding programs. Application of QIPDT to the 49 European maize inbreds resulted in a significant (P < 0.05) association located at a position for which a consensus QTL was detected in a previous study. The results of our study suggested that QIPDT is a promising QTL detection method for data collected routinely in plant breeding programs.     

The Genomic Impacts of Drift and Selection for Hybrid Performance in Maize

Although maize is naturally an outcrossing organism, modern breeding utilizes highly inbred lines in controlled crosses to produce hybrids. The U.S. Department of Agriculture’s reciprocal recurrent selection experiment between the Iowa Stiff Stalk Synthetic (BSSS) and the Iowa Corn Borer Synthetic No. 1 (BSCB1) populations represents one of the longest running experiments to understand the response to selection for hybrid performance. To investigate the genomic impact of this selection program, we genotyped the progenitor lines and >600 individuals across multiple cycles of selection using a genome-wide panel of ∼40,000 SNPs. We confirmed previous results showing a steady temporal decrease in genetic diversity within populations and a corresponding increase in differentiation between populations. Thanks to detailed historical information on experimental design, we were able to perform extensive simulations using founder haplotypes to replicate the experiment in the absence of selection. These simulations demonstrate that while most of the observed reduction in genetic diversity can be attributed to genetic drift, heterozygosity in each population has fallen more than expected. We then took advantage of our high-density genotype data to identify extensive regions of haplotype fixation and trace haplotype ancestry to single founder inbred lines. The vast majority of regions showing such evidence of selection differ between the two populations, providing evidence for the dominance model of heterosis. We discuss how this pattern is likely to occur during selection for hybrid performance and how it poses challenges for dissecting the impacts of modern breeding and selection on the maize genome.     

Dissection of complex traits in forest trees — opportunities for marker-assisted selection

Due to their long reproductive cycles and the time to expression of mature traits, marker-assisted selection is particularly attractive for tree breeding. In this review, we discuss different approaches used for developing markers and propose a method for application of markers in low linkage disequilibrium (LD) populations. Identification of useful markers for application in tree breeding is mainly based on two approaches, quantitative trait locus (QTL) mapping and association genetic studies. While several studies have identified significant markers, effect of the individual markers is low making it difficult to utilize them in breeding programs. Recently, genomic selection (GS) was proposed for overcoming some of these difficulties. In GS, high density markers are used for predicting phenotypes from genotypes. Currently small effective populations with high LD are being tested for GS in tree breeding. For wider application, GS needs to be applied in low LD populations which are found in many tree breeding programs. Here we propose an approach in which the significant markers from association studies may be used for developing prediction models in low LD populations using the same methods as in GS. Preliminary analyses indicate that a modest numbers of markers may be sufficient for developing prediction models in low LD populations. GS based on large numbers of random markers or small numbers of associated markers is poised to make marker-assisted selection a reality in forest tree breeding.     

Randomly distributed crossovers may generate block-like patterns of linkage disequilibrium: an act of genetic drift

There is considerable interest in identifying and characterizing block-like patterns of linkage disequilibrium (LD; haplotype blocks) in the human genome as these may facilitate the identification of complex disease genes via genome-wide association studies. Although recombination hot-spots have been suggested as the primary mechanism to explain the block-like pattern of LD, other forces, such as genetic drift, may also be important. To this end, we have studied the effect of various recombination models on patterns of LD by using extensive simulations. As expected, haplotype blocks were observed under a model allowing recombination hot-spots. However, we also observed similar block-like patterns in the models where recombination crossovers are randomly and uniformly distributed, and we demonstrate that these blocks are generated by genetic drift. We caution that genetic drift may be an alternative mechanism (in addition to recombination hot-spots) that can lead to block-like patterns of LD. Our findings highlight the necessity of characterizing haplotype blocks in world-wide populations.     

Evidence for a large-scale population structure of Arabidopsis thaliana from genome-wide single nucleotide polymorphism markers

Population-based methods for the genetic mapping of adaptive traits and the analysis of natural selection require that the population structure and demographic history of a species are taken into account. We characterized geographic patterns of genetic variation in the model plant Arabidopsis thaliana by genotyping 115 genome-wide single nucleotide polymorphism (SNP) markers in 351 accessions from the whole species range using a matrix-assisted laser desorption/ionization time-of-flight assay, and by sequencing of nine unlinked short genomic regions in a subset of 64 accessions. The observed frequency distribution of SNPs is not consistent with a constant-size neutral model of sequence polymorphism due to an excess of rare polymorphisms. There is evidence for a significant population structure as indicated by differences in genetic diversity between geographic regions. Accessions from Central Asia have a low level of polymorphism and an increased level of genome-wide linkage disequilibrium (LD) relative to accessions from the Iberian Peninsula and Central Europe. Cluster analysis with the structure program grouped Eurasian accessions into K=6 clusters. Accessions from the Iberian Peninsula and from Central Asia constitute distinct populations, whereas Central and Eastern European accessions represent admixed populations in which genomes were reshuffled by historical recombination events. These patterns likely result from a rapid postglacial recolonization of Eurasia from glacial refugial populations. Our analyses suggest that mapping populations for association or LD mapping should be chosen from regional rather than a species-wide sample or identified genetically as sets of individuals with similar average genetic distances. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

玉米株高和穗位高的遗传基础与分子机制*

株高和穗位高是玉米重要育种性状,直接影响植株的养分利用效率及抗倒伏性,进而影响玉米产量。玉米株高和穗位高属于典型数量性状,目前通过数量性状位点(quantitative trait loci mapping,QTL)定位和全基因组关联分析(genome-wide association study, GWAS)等方法已挖掘到较多相关遗传位点,通过QTL精细定位及利用突变体克隆了一些调控株高和穗位高关键基因。但是由于各研究组所利用的群体类型和大小、标记类型和密度以及统计方法不同,所鉴定QTL差异较大,单个研究难以揭示玉米株高和穗位高遗传结构。早期QTL定位的结果多以遗传距离来展示,不同时期GWAS研究所使用参考基因组版本不同,这进一步增加了借鉴和利用前人研究结果的难度。首次将目前已鉴定株高和穗位高遗传定位信息统一锚定至玉米自交系B73参考基因组V4版本,构建了株高和穗位高性状定位的一致性图谱,并鉴定出可被多个独立研究定位的热点区间。进一步对已克隆玉米株高和穗位高调控基因进行总结与分类,揭示株高和穗位高性状调控机制,对深度解析株高和穗位高遗传结构、指导基因克隆和利用分子标记辅助选择优化玉米株高和穗位高性状均具有重要意义。     

Linkage disequilibrium in two European F<Subscript>2</Subscript> flint maize populations under modified recurrent full-sib selection

According to quantitative genetic theory, linkage disequilibrium (LD) can hamper the short- and long-term selection response in recurrent selection (RS) programs. We analyzed LD in two European flint maize populations, KW1265 × D146 (A × B) and D145 × KW1292 (C × D), under modified recurrent full-sib selection. Our objectives were to investigate (1) the decay of initial parental LD present in F₂ populations by three generations of intermating, (2) the generation of new LD in four (A × B) and seven (C × D) selection cycles, and (3) the relationship between LD changes and estimates of the additive genetic variance. We analyzed the F₂ and the intermated populations as well as all selection cycles with 104 (A × B) and 101 (C × D) simple sequence repeat (SSR) markers with a uniform coverage of the entire maize genome. The LD coefficient D and the composite LD measure Δ were estimated and significance tests for LD were performed. LD was reduced by intermating as expected from theory. A directional generation of negative LD between favorable alleles could not be observed during the selection cycles. However, considerable undirectional changes in D were observed, which we attributed to genetic sampling due to the finite population size used for recombination. Consequently, a long-term reduction of the additive genetic variance due to negative LD was not observed. Our experimental results support the hypothesis that in practical RS programs with maize, LD generated by selection is not a limiting factor for obtaining a high selection response.     

Association mapping in an elite maize breeding population

Association mapping (AM) is a powerful approach to dissect the genetic architecture of quantitative traits. The main goal of our study was to empirically compare several statistical methods of AM using data of an elite maize breeding program with respect to QTL detection power and possibility to correct for population stratification. These models were based on the inclusion of cofactors (Model A), cofactors and population effect (Model B), and SNP effects nested within populations (Model C). A total of 930 testcross progenies of an elite maize breeding population were field-evaluated for grain yield and grain moisture in multi-location trials and fingerprinted with 425 SNP markers. For grain yield, population stratification was effectively controlled by Model A. For grain moisture with a high ratio of variance among versus within populations, Model B should be applied in order to avoid potential false positives. Model C revealed large differences among allele substitution effects for trait-associated SNPs across multiple plant breeding populations. This heterogeneous SNP allele substitution effects have a severe impact for genomic selection studies, where SNP effects are often assumed to be independent of the genetic background.     

Genetic diversity of maize kernel starch-synthesis genes with SNAPs.

Measuring genetic diversity in populations of a crop species is very important for understanding the genetic structure of and subsequently improving the crop species by genetic manipulation. Single-nucleotide amplified polymorphisms (SNAPs) among and within maize populations of waxy, dent, and sweet corns at 25 single-nucleotide polymorphism (SNP) sites in 6 kernel starch-synthesis genes (sh2, bt2, su1, ae1, wx1, and sh1) were determined. Because of the intensive selection of some favorable alleles in starch-synthesis genes during the breeding process, and the resultant strong linkage disequilibrium (LD), the number of haplotypes in each population was far less than expected. Subsequent phenetic clustering analysis with the SNAPs indicated that the dent, waxy, and sweet corns formed distinct subclusters, except in a few incidences. LD was surveyed among SNAPs of intragenic, intergenic, and intrachromosomal SNPs in whole and subpopulations, which revealed that some SNAPs showed high LD with many other SNAPs, but some SNAPs showed low or no significant LD with others, depending on the subpopulation, indicating that these starch genes have undergone different selection in each subpopulation during the breeding process. Because the starch synthesis genes used in this study are important in maize breeding, the genetic diversity, LD, and accessions having rare SNAP alleles might be valuable in maize improvement programs.     

Genetic and Morphometric Divergence of an Invasive Bird: The Introduced House Sparrow (Passer domesticus) in Brazil

Introduced species are interesting systems for the study of contemporary evolution in new environments because of their spatial and temporal scales. For this study we had three aims: (i) to determine how genetic diversity and genetic differentiation of introduced populations of the house sparrow (Passer domesticus) in Brazil varies with range expansion, (ii) to determine how genetic diversity and differentiation in Brazil compares to ancestral European populations; and (iii) to determine whether selection or genetic drift has been more influential on phenotypic divergence. We used six microsatellite markers to genotype six populations from Brazil and four populations from Europe. We found slightly reduced levels of genetic diversity in Brazilian compared to native European populations. However, among introduced populations of Brazil, we found no association between genetic diversity and time since introduction. Moreover, overall genetic differentiation among introduced populations was low indicating that the expansion took place from large populations in which genetic drift effects would likely have been weak. We found significant phenotypic divergence among sites in Brazil. Given the absence of a spatial genetic pattern, divergent selection and not genetic drift seems to be the main force behind most of the phenotypic divergence encountered. Unravelling whether microevolution (e.g., allele frequency change), phenotypic plasticity, or both mediated phenotypic divergence is challenging and will require experimental work (e.g., common garden experiments or breeding programs).     

Genome-Wide Association Mapping for Yield and Other Agronomic Traits in an Elite Breeding Population of Tropical Rice (Oryza sativa)

Genome-wide association mapping studies (GWAS) are frequently used to detect QTL in diverse collections of crop germplasm, based on historic recombination events and linkage disequilibrium across the genome. Generally, diversity panels genotyped with high density SNP panels are utilized in order to assay a wide range of alleles and haplotypes and to monitor recombination breakpoints across the genome. By contrast, GWAS have not generally been performed in breeding populations. In this study we performed association mapping for 19 agronomic traits including yield and yield components in a breeding population of elite irrigated tropical rice breeding lines so that the results would be more directly applicable to breeding than those from a diversity panel. The population was genotyped with 71,710 SNPs using genotyping-by-sequencing (GBS), and GWAS performed with the explicit goal of expediting selection in the breeding program. Using this breeding panel we identified 52 QTL for 11 agronomic traits, including large effect QTLs for flowering time and grain length/grain width/grain-length-breadth ratio. We also identified haplotypes that can be used to select plants in our population for short stature (plant height), early flowering time, and high yield, and thus demonstrate the utility of association mapping in breeding populations for informing breeding decisions. We conclude by exploring how the newly identified significant SNPs and insights into the genetic architecture of these quantitative traits can be leveraged to build genomic-assisted selection models.     

The effect of gametic-phase disequilibrium on the prediction of response to recurrent selection in plants

Selection reduces additive genetic variation by generating gametic-phase disequilibrium, a phenomenon largely ignored when predicting response in plant breeding programs. The development of gametic-phase disequilibrium is here taken into account when predicting the response to selection for various schemes of recurrent selection applicable to plant populations. A general program permits prediction of response to selection from schemes of recurrent selection in which two or more rounds of selection occur in each cycle. An example from Sugar Beet, with alternate rounds of half-sib and S1 family selection, is illustrated. It is shown that failure to take into account the effects of gametic-phase disequilibrium can result in substantial overestimation of the response to selection as well as to changes in rank of the merits of alternative breeding schemes. For a given scheme, ignoring gametic-phase disequilibrium has only small effects on defining the optimum allocation of plots and the numbers of families tested.     

Genomic Diversity in Pig (Sus scrofa) and its Comparison with Human and other Livestock

We have reviewed the current pig (Sus scrofa) genomic diversity within and between sites and compared them with human and other livestock. The current Porcine 60K single nucleotide polymorphism (SNP) panel has an average SNP distance in a range of 30 - 40 kb. Most of genetic variation was distributed within populations, and only a small proportion of them existed between populations. The average heterozygosity was lower in pig than in human and other livestock. Genetic inbreeding coefficient (F(IS)), population differentiation (F(ST)), and Nei's genetic distance between populations were much larger in pig than in human and other livestock. Higher average genetic distance existed between European and Asian populations than between European or between Asian populations. Asian breeds harboured much larger variability and higher average heterozygosity than European breeds. The samples of wild boar that have been analyzed displayed more extensive genetic variation than domestic breeds. The average linkage disequilibrium (LD) in improved pig breeds extended to 1 - 3 cM, much larger than that in human (~ 30 kb) and cattle (~ 100 kb), but smaller than that in sheep (~ 10 cM). European breeds showed greater LD that decayed more slowly than Asian breeds. We briefly discuss some processes for maintaining genomic diversity in pig, including migration, introgression, selection, and drift. We conclude that, due to the long time of domestication, the pig possesses lower heterozygosity, higher F(IS), and larger LD compared with human and cattle. This implies that a smaller effective population size and less informative markers are needed in pig for genome wide association studies.     

Genetic drift and selection effects of modified recurrent full-sib selection programs in two F2 populations of European flint maize

Selection response of a modified recurrent full-sib (FS) selection scheme conducted in two European flint F₂ maize (Zea mays L.) populations was re-evaluated. Our objectives were to (1) determine the selection response for per se and testcross performance in both populations and (2) separate genetic effects due to selection from those due to random genetic drift. Modified recurrent FS selection was conducted at three locations using an effective population size N _e = 32 and a selection rate of 25% for a selection index, based on grain yield and grain moisture. Recombination was performed according to a pseudo-factorial mating scheme. Selection response was assessed using a population diallel including the source population and advanced selection cycles, as well as testcrosses with unrelatesd inbred line testers and the parental F₁ generation. Selection response per cycle was significant for grain yield and grain moisture in both populations. Effects of random genetic drift caused only a small reduction in the selection response. No significant selection response was observed for testcrosses, suggesting that for heterotic traits, such as grain yield, a high frequency of favorable alleles in the elite tester masked the effects of genes segregating in the populations. We conclude that our modified recurrent FS selection is an alternative to other commonly applied intrapopulation recurrent selection schemes, and some of its features may also be useful for increasing the efficiency of interpopulation recurrent selection programs.     

High-density linkage mapping of vitamin E content in maize grain

Vitamin E refers to eight distinct compounds collectively known as tocochromanols and can be further divided into two classes, tocotrienols and tocopherols. Tocochromanols are the major lipid-soluble antioxidants in maize (Zea mays L.) grain. Enhancing vitamin E content of maize through plant breeding has important implications for human and animal nutrition. Four inbred lines exhibiting unique variation for tocochromanol compounds were chosen from the Goodman maize diversity panel to construct two biparental mapping populations (N6xNC296 and E2558xCo125). The N6xNC296 population was developed to analyze segregation for α-tocopherol and α-tocotrienol content. The E2558WxCo125 population was developed to analyze segregation for the ratio of total tocotrienols to tocopherols. The tocochromanol variation in two replicates of each population was quantified using liquid chromatography-diode array detection. Using high-density linkage mapping, novel quantitative trait loci (QTL) in the N6xNC296 population were mapped using tocopherol ratio traits. These QTL contain the candidate gene homogentisate phytyltransferase (ZmVTE2) within the respective support intervals. This locus was not mapped in a previous genome-wide association study that analyzed tocochromanols in the Goodman diversity panel. Transgressive segregation was observed for γ- and α-tocochromanols in these populations, which facilitated QTL identification. These QTL and transgressive segregant families can be used in selection programs for vitamin E enhancement in maize. This work illustrates the complementary nature of biparental mapping populations and genome-wide association studies to further characterize genetic variation of tocochromanol content in maize grain.     

Allelic Association: Linkage Disequilibrium Structure and Gene Mapping

The linkage disequilibrium (LD) structure of the human genome is now well understood and characterised for a number of human populations. The LD structure underpins the design and execution of candidate gene and genome-wide association mapping studies. Successful association mapping studies completed to date provide vital new insights into the genetic influences on common diseases, such as diabetes, some cancers and heart disease. The LD structure also presents new avenues of research into the genetic history of human populations, the effects of natural selection and the impact of recombination on the genomic landscape. This review introduces this exciting and complex field by encompassing this range of topics.