Linkage disequilibrium for different scales and applications

Assessing the patterns of linkage disequilibrium (LD) has become an important issue in both evolutionary biology and medical genetics since the rapid accumulation of densely spaced DNA sequence variation data in several organisms. LD deals with the correlation of genetic variation at two or more loci or sites in the genome within a given population. There are a variety of LD measures which range from traditional pairwise LD measures such as D' or r2 to entropy-based multi-locus measures or haplotype-specific approaches. Understanding the evolutionary forces (in particular recombination) that generate the observed variation of LD patterns across genomic regions is addressed by model-based LD analysis. Marker type and its allelic composition also influence the observed LD pattern, microsatellites having a greater power to detect LD in population isolates than SNPs. This review aims to explain basic LD measures and their application properties.     

Linkage disequilibrium structures in cattle and their application to breed identification testing

We examined the extent of linkage disequilibrium (LD) block lengths in four breed populations: Japanese Black, Angus, Hereford and Holstein. Three chromosomal regions in which QTL were previously mapped in Japanese Black populations were scanned with 84 microsatellite markers. The estimated LD lengths in these four purebred populations varied from 535 to 683 kb, which is much shorter than the values reported previously. Our findings suggest that QTL can be mapped in sub-centimorgan regions in these populations using an LD-mapping method. We also developed breed identification methods to distinguish Japanese Black from Angus, Hereford, Holstein and F₁ animals (Japanese Black × Holstein) respectively using the haplotypic frequencies of a pair of markers in the breed populations. After assessing the distributions of posterior probabilities to be Japanese Black, we obtained several pairs of markers that completely distinguished Japanese Black from the other breeds. We also obtained several combinations of six markers that completely distinguished Japanese Black animals from F₁ animals.     

Linkage disequilibrium and signatures of selection on chromosomes 19 and 29 in beef and dairy cattle

The objective of this study was to quantify the extent of linkage disequilibrium (LD) on bovine chromosomes 19 and 29 and to study the pattern of selection signatures in beef and dairy breeds (Angus and Holstein) of Bos taurus. The extent of LD was estimated for 370 and 186 single nucleotide polymorphism markers on BTA19 and 29 respectively using the square of the correlation coefficient (r(2)) among alleles at pairs of loci. A comparison of the extent of LD found that the decline of LD followed a similar pattern in both breeds. We observed long-range LD and found that LD dissipates to background levels at a locus separation of about 20 Mb on both chromosomes. Along each chromosome, patterns of LD were variable in both breeds. We find that a minimum of 30 000 informative and evenly spaced markers would be required for whole-genome association studies in cattle. In addition, we have identified chromosomal regions that show some evidence of selection for economically important traits in Angus and Holstein cattle. The results of this study are of importance for the design and application of association studies.     

Linkage disequilibrium and association studies in higher plants: Present status and future prospects

During the last two decades, DNA-based molecular markers have been extensively utilized for a variety of studies in both plant and animal systems. One of the major uses of these markers is the construction of genome-wide molecular maps and the genetic analysis of simple and complex traits. However, these studies are generally based on linkage analysis in mapping populations, thus placing serious limitations in using molecular markers for genetic analysis in a variety of plant systems. Therefore, alternative approaches have been suggested, and one of these approaches makes use of linkage disequilibrium (LD)-based association analysis. Although this approach of association analysis has already been used for studies on genetics of complex traits (including different diseases) in humans, its use in plants has just started. In the present review, we first define and distinguish between LD and association mapping, and then briefly describe various measures of LD and the two methods of its depiction. We then give a list of different factors that affect LD without discussing them, and also discuss the current issues of LD research in plants. Later, we also describe the various uses of LD in plant genomics research and summarize the present status of LD research in different plant genomes. In the end, we discuss briefly the future prospects of LD research in plants, and give a list of softwares that are useful in LD research, which is available as electronic supplementary material (ESM) Electronic supplementary material Electronic supplementary material is available for this article at and accessible for authorised users.     

Extensive intragenic recombination and patterns of linkage disequilibrium at the CSN3 locus in European rabbit

Kappa-casein (CSN3) plays an important role in stabilising the Ca-sensitive caseins in the micelle. The European rabbit (Oryctolagus cuniculus) CSN3 has previously been shown to possess two alleles (A and B), which differ deeply in their intronic regions (indels of 100 and 1550 nucleotides in introns 1 and 4, respectively). Furthermore, a correlation between several reproductive performance traits and the different alleles was described. However, all these data were exclusively collected in rabbit domestic breeds, preventing a deeper understanding of the extensive polymorphism observed in the CSN3 gene. Additionally, the techniques available for the typing of both indel polymorphisms were until now not suitable for large-scale studies. In this report, we describe a simple, PCR-based typing method to distinguish rabbit CSN3 alleles. We analyse both ancient wild rabbit populations from the Iberian Peninsula and France, and the more recently derived English wild rabbits and domestic stocks. A new allele (C) showing another major indel (250 bp) in intron 1 was found, but exclusively detected in Iberian wild rabbits. In addition, our survey revealed the occurrence of new haplotypes in wild populations, suggesting that intragenic recombination is important in creating genetic diversity at this locus. This easy and low cost single-step PCR-based method results in an improvement over previous described techniques, can be easily set up in a routine molecular laboratory and would probably be a valuable tool in the management of rabbit domestic breeds.     

Linkage disequilibrium between alleles at highly polymorphic mini- and micro-satellite loci of Theileria parva isolated from cattle in three regions of Kenya

Theileria parva schizont-infected lymphocyte culture isolates from western, central and coastal Kenya were analysed for size polymorphism at 30 T. parva-specific variable number tandem repeat (VNTR) loci using a panel of mini- and micro-satellite markers. The mean number of alleles ranged from 3 to 11 at individual loci and 183 distinct alleles were observed in total, indicating high genetic diversity within the T. parva gene pool in Kenyan cattle. The frequency distribution of the length variation of specific alleles among isolates ranged from normal to markedly discontinuous. Genetic relationships between isolates were analysed using standard indices of genetic distance. Genetic distances and dendrograms derived from these using neighbour-joining algorithms did not indicate significant clustering on a geographical basis. Analysis of molecular variance demonstrated that the genetic variation between individual isolates was 72%, but only 2.3% when isolates from different regions were pooled. Both these observations suggest minimal genetic sub-structuring relative to geographical origin. Linkage disequilibrium was observed between pairs of loci within populations, as in certain Ugandan T. parva populations. A novel observation was that disequilibrium was also detected between alleles at three individual pairs of VNTR loci when isolates from the three regional meta-populations were pooled for analysis.     

The heterogeneous levels of linkage disequilibrium in white spruce genes and comparative analysis with other conifers

In plants, knowledge about linkage disequilibrium (LD) is relevant for the design of efficient single-nucleotide polymorphism arrays in relation to their use in population and association genomics studies. Previous studies of conifer genes have shown LD to decay rapidly within gene limits, but exceptions have been reported. To evaluate the extent of heterogeneity of LD among conifer genes and its potential causes, we examined LD in 105 genes of white spruce (Picea glauca) by sequencing a panel of 48 haploid megagametophytes from natural populations and further compared it with LD in other conifer species. The average pairwise r(2) value was 0.19 (s.d.=0.19), and LD dropped quickly with a half-decay being reached at a distance of 65 nucleotides between sites. However, LD was significantly heterogeneous among genes. A first group of 29 genes had stronger LD (mean r(2)=0.28), and a second group of 38 genes had weaker LD (mean r(2)=0.12). While a strong relationship was found with the recombination rate, there was no obvious relationship between LD and functional classification. The level of nucleotide diversity, which was highly heterogeneous across genes, was also not significantly correlated with LD. A search for selection signatures highlighted significant deviations from the standard neutral model, which could be mostly attributed to recent demographic changes. Little evidence was seen for hitchhiking and clear relationships with LD. When compared among conifer species, on average, levels of LD were similar in genes from white spruce, Norway spruce and Scots pine, whereas loblolly pine and Douglas fir genes exhibited a significantly higher LD.     

An extensive candidate gene approach to speciation: diversity,divergence and linkage disequilibrium in candidate pigmentation genes across the European crow hybrid zone

Colouration patterns have an important role in adaptation and speciation. The European crow system, in which all-black carrion crows and grey-coated hooded crows meet in a narrow hybrid zone, is a prominent example. The marked phenotypic difference is maintained by assortative mating in the absence of neutral genetic divergence, suggesting the presence of few pigmentation genes of major effect. We made use of the rich phenotypic and genetic resources in mammals and identified a comprehensive panel of 95 candidate pigmentation genes for birds. Based on functional annotation, we chose a subset of the most promising 37 candidates, for which we developed a marker system that demonstrably works across the avian phylogeny. In total, we sequenced 107 amplicons (∼3 loci per gene, totalling 60 kb) in population samples of crows (n=23 for each taxon). Tajima''s D, Fu''s F_S, DHEW and HKA (Hudson–Kreitman–Aguade) statistics revealed several amplicons that deviated from neutrality; however, none of these showed significantly elevated differentiation between the two taxa. Hence, colour divergence in this system may be mediated by uncharacterized pigmentation genes or regulatory regions outside genes. Alternatively, the observed high population recombination rate (4N_er∼0.03), with overall linkage disequilibrium dropping rapidly within the order of few 100 bp, may compromise the power to detect causal loci with nearby markers. Our results add to the debate as to the utility of candidate gene approaches in relation to genomic features and the genetic architecture of the phenotypic trait in question.     

The extent of linkage disequilibrium in a large cattle population of western Africa and its consequences for association studies

Several previous studies concluded that linkage disequilibrium (LD) in livestock populations from developed countries originated from the impact of strong selection. Here, we assessed the extent of LD in a cattle population from western Africa that was bred in an extensive farming system. The analyses were performed on 363 individuals in a Bos indicus x Bos taurus population using 42 microsatellite markers on BTA04, BTA07 and BTA13. A high level of expected heterozygosity (0.71), a high mean number of alleles per locus (9.7) and a mild shift in Hardy-Weinberg equilibrium were found. Linkage disequilibrium extended over shorter distances than what has been observed in cattle from developed countries. Effective population size was assessed using two methods; both methods produced large values: 1388 when considering heterozygosity (assuming a mutation rate of 10(-3)) and 2344 when considering LD on whole linkage groups (assuming a constant population size over generations). However, analysing the decay of LD as a function of marker spacing indicated a decreasing trend in effective population size over generations. This decrease could be explained by increasing selective pressure and/or by an admixture process. Finally, LD extended over small distances, which suggested that whole-genome scans will require a large number of markers. However, association studies using such populations will be effective.     

Genetic variation, population structure, and linkage disequilibrium in European elite germplasm of perennial ryegrass

Perennial ryegrass (Lolium perenne L.) is a highly valued temperate climate grass species grown as forage crop and for amenity uses. Due to its outbreeding nature and recent domestication, a high degree of genetic diversity is expected among cultivars. The aim of this study was to assess the extent of linkage disequilibrium (LD) within European elite germplasm and to evaluate the appropriate methodology for genetic association mapping in perennial ryegrass. A high level of genetic diversity was observed in a set of 380 perennial ryegrass elite genotypes when genotyped with 40 SSRs and 2 STS markers. A Bayesian structure analysis identified two subpopulations, which were confirmed by principal coordinate analysis (PCoA). One subpopulation consisted mainly of genotypes originating from the UK, while germplasm mostly from Continental Europe was grouped into the second subpopulation. LD (r²) decay was rapid and occurred within 0.4 cM across European varieties, when population structure was taken into consideration. However, an extended LD of up to 6.6 cM was detected within the variety Aberdart. High genetic diversity and rapid LD decay provide means for high resolution association mapping in elite materials of perennial ryegrass. However, different strategies need to be applied depending on the material used. Genome-wide association study (GWAS) with several hundred markers can be applied within synthetic varieties to identify large (up to 10 cM) genomic regions affecting trait variation. A combination of available and novel DNA markers is needed to achieve resolution required for GWAS in elite breeding materials. An even higher marker density of several million SNPs might be needed for GWAS in diverse ecotype collections, potentially resulting in quantitative trait polymorphism (QTP) identification.     

Stenosternus Karsch,a possible link between Neotropical and Afrotropical Orphninae (Coleoptera,Scarabaeidae)

The monotypical orphnine genus Stenosternus Karsch is known from a single specimen of Stenosternus costatus collected on the São Tomé island (Gulf of Guinea). The holotype of Stenosternus costatus Karsch is re-examined and its characters are discussed and illustrated. Although the genus was implicitly placed by Paulian (1984) in the Old World tribe Orphnini Erichson, re-examination of the holotype of Stenosternus costatus shows that it has characters similar to those of the members of the New World tribe Aegidiini Paulian. Placement of Stenosternus in the Aegidiini is supported by the metepisternum widened posteriorly (forming posterior metepisternal lock for closed elytra) and a keel separating basal and anterolateral parts of propleurae. Relationships of Stenosternus with other orphnine taxa and possible ways of origin of São Toméan orphnine fauna are discussed.     

Linkage Disequilibrium with Linkage Analysis of Multiline Crosses Reveals Different Multiallelic QTL for Hybrid Performance in the Flint and Dent Heterotic Groups of Maize

Multiparental designs combined with dense genotyping of parents have been proposed as a way to increase the diversity and resolution of quantitative trait loci (QTL) mapping studies, using methods combining linkage disequilibrium information with linkage analysis (LDLA). Two new nested association mapping designs adapted to European conditions were derived from the complementary dent and flint heterotic groups of maize (Zea mays L.). Ten biparental dent families (N = 841) and 11 biparental flint families (N = 811) were genotyped with 56,110 single nucleotide polymorphism markers and evaluated as test crosses with the central line of the reciprocal design for biomass yield, plant height, and precocity. Alleles at candidate QTL were defined as (i) parental alleles, (ii) haplotypic identity by descent, and (iii) single-marker groupings. Between five and 16 QTL were detected depending on the model, trait, and genetic group considered. In the flint design, a major QTL (R² = 27%) with pleiotropic effects was detected on chromosome 10, whereas other QTL displayed milder effects (R² < 10%). On average, the LDLA models detected more QTL but generally explained lower percentages of variance, consistent with the fact that most QTL display complex allelic series. Only 15% of the QTL were common to the two designs. A joint analysis of the two designs detected between 15 and 21 QTL for the five traits. Of these, between 27 for silking date and 41% for tasseling date were significant in both groups. Favorable allelic effects detected in both groups open perspectives for improving biomass production.     

Single nucleotide polymorphisms generated by genotyping by sequencing to characterize genome-wide diversity,linkage disequilibrium,and selective sweeps in cultivated watermelon

Background

A large single nucleotide polymorphism (SNP) dataset was used to analyze genome-wide diversity in a diverse collection of watermelon cultivars representing globally cultivated, watermelon genetic diversity. The marker density required for conducting successful association mapping depends on the extent of linkage disequilibrium (LD) within a population. Use of genotyping by sequencing reveals large numbers of SNPs that in turn generate opportunities in genome-wide association mapping and marker-assisted selection, even in crops such as watermelon for which few genomic resources are available. In this paper, we used genome-wide genetic diversity to study LD, selective sweeps, and pairwise F_ST distributions among worldwide cultivated watermelons to track signals of domestication.

Results

We examined 183 Citrullus lanatus var. lanatus accessions representing domesticated watermelon and generated a set of 11,485 SNP markers using genotyping by sequencing. With a diverse panel of worldwide cultivated watermelons, we identified a set of 5,254 SNPs with a minor allele frequency of ≥ 0.05, distributed across the genome. All ancestries were traced to Africa and an admixture of various ancestries constituted secondary gene pools across various continents. A sliding window analysis using pairwise F_ST values was used to resolve selective sweeps. We identified strong selection on chromosomes 3 and 9 that might have contributed to the domestication process. Pairwise analysis of adjacent SNPs within a chromosome as well as within a haplotype allowed us to estimate genome-wide LD decay. LD was also detected within individual genes on various chromosomes. Principal component and ancestry analyses were used to account for population structure in a genome-wide association study. We further mapped important genes for soluble solid content using a mixed linear model.

Conclusions

Information concerning the SNP resources, population structure, and LD developed in this study will help in identifying agronomically important candidate genes from the genomic regions underlying selection and for mapping quantitative trait loci using a genome-wide association study in sweet watermelon.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-767) contains supplementary material, which is available to authorized users.     

Genome wide linkage disequilibrium in Chinese asparagus bean (Vigna. unguiculata ssp. sesquipedialis) germplasm: implications for domestication history and genome wide association studies

Association mapping of important traits of crop plants relies on first understanding the extent and patterns of linkage disequilibrium (LD) in the particular germplasm being investigated. We characterize here the genetic diversity, population structure and genome wide LD patterns in a set of asparagus bean (Vigna. unguiculata ssp. sesquipedialis) germplasm from China. A diverse collection of 99 asparagus bean and normal cowpea accessions were genotyped with 1127 expressed sequence tag-derived single nucleotide polymorphism markers (SNPs). The proportion of polymorphic SNPs across the collection was relatively low (39%), with an average number of SNPs per locus of 1.33. Bayesian population structure analysis indicated two subdivisions within the collection sampled that generally represented the 'standard vegetable' type (subgroup SV) and the 'non-standard vegetable' type (subgroup NSV), respectively. Level of LD (r(2)) was higher and extent of LD persisted longer in subgroup SV than in subgroup NSV, whereas LD decayed rapidly (0-2 cM) in both subgroups. LD decay distance varied among chromosomes, with the longest (≈ 5 cM) five times longer than the shortest (≈ 1 cM). Partitioning of LD variance into within- and between-subgroup components coupled with comparative LD decay analysis suggested that linkage group 5, 7 and 10 may have undergone the most intensive epistatic selection toward traits favorable for vegetable use. This work provides a first population genetic insight into domestication history of asparagus bean and demonstrates the feasibility of mapping complex traits by genome wide association study in asparagus bean using a currently available cowpea SNPs marker platform.     

Expression,SNV identification,linkage disequilibrium,and combined genotype association analysis of the muscle-specific gene CSRP3 in Chinese cattle

The cysteine and glycine-rich protein 3 (CSRP3) plays an important role in the myofiber differentiation. Here, we identified five SNVs in all exon and intron regions of the CSRP3 gene using DNA sequencing, PCR-RFLP and forced-PCR-RFLP methods in 554 cattle. Four of the five SNVs were significantly associated with growth performance and carcass traits of the cattle. In addition, we evaluated haplotype frequency and linkage disequilibrium coefficient of five sequence variants. The result of haplotype analysis demonstrated 28 haplotypes present in Qinchuan and two haplotypes in Chinese Holstein. Only haplotypes 1 and 8 were being shared by two populations, haplotype 14 had the highest haplotype frequency in Qinchuan (17.4%) and haplotype 8 had the highest haplotype frequency in Chinese Holstein (94.4%). Statistical analyses of combined genotypes indicated that some combined genotypes were significantly or highly significantly associated with growth and carcass traits in the Qinchuan cattle population. qPCR analyses also showed that bovine CSRP3 gene was exclusively expressed in longissimus dorsi muscle and heart tissues. The data support the high potential of the CSRP3 as a marker gene for the improvement of growth performance and carcass traits in selection programs.     

Linkage exclusion analysis of two candidate regions on chromosomes 7 and 11: leptin and UCP2/UCP3 are not QTLs for obesity in US Caucasians

Leptin (LEP) and the uncoupling proteins 2 and 3 (UCP2/UCP3) are key molecules involved in the regulation of food intake and energy expenditure. However, their contribution to variation of obesity phenotypes in the general population remains controversial. The present study is to investigate whether chromosomal regions 7q and 11q, which contain LEP and UCP2/UCP3, respectively, can be excluded for linkage with obesity phenotypes. The obesity phenotypes include body mass index (BMI), fat mass, and percentage fat mass (PFM), with the latter two measured by dual-energy X-ray absorptiometry. We conducted exclusion linkage analyses using a variance component approach in a sample of 1816 individuals coming from 79 extended Caucasian pedigrees. In this study, we were able to exclude chromosomal region 7q containing LEP as having an effect on fat mass and PFM at effect sizes of 5% or greater, and on BMI at effect sizes of 10% or greater. We were able to exclude chromosomal region 11q containing UCP2/UCP3 as having an effect on fat mass and PFM at effect sizes of 10% or greater, and on BMI at effect sizes of 5% or greater. Our results suggest that the LEP and UCP2/UCP3 genes are unlikely to have a substantial effect on variation in obesity phenotypes in this particular US Caucasian population.     

The identification and quantification of three new 6alpha-hydroxylated corticosteroids in human neonatal urine

The 24 hours urines of six two days old fullterm newborn infants were investigated for polar corticosteroids. 6alpha-hydroxy-tetrahydrocortisone, 6alpha-hydroxy-20alpha-cortolone and 6alpha-hydroxy-20beta-cortolone were identified by gas chromatographic-mass spectrometric comparison of the urinary steroids to compounds synthesized previously. These 6alpha-hydroxylated corticosteroids as well as seven other polar corticosteroids were quantified by gas chromatography or mass fragmentography. It was shown that the newly identified steroids constituted a quantitatively important part of the neonatal urinary corticosteroids. The unconjugated- and glucuronic acid conjugated steroids were quantified separately. It was found that the extent of glucuronoconjugation decreased with increasing polarity of the steroid moiety.