Haplotype block and superblock structures of the alpha1-adrenergic receptor genes reveal echoes from the chromosomal past

A significant proportion of the human genome is contained within haplotype blocks across which pairwise linkage disequilibrium (LD) is very high. However, LD is also often high between markers at more remote distances, and within different haplotype blocks. Here, we evaluate the origins of haplotype block structure in the three genes for alpha1 adrenergic receptors (alpha1-AR) in the human genome ( ADRA1A, ADRA1B and ADRA1D) by genotyping dense single-nucleotide polymorphism (SNP) marker maps, and show that LD signals between distant markers are due to the presence of extended haplotype superblocks in individuals with ancient chromosomes which have escaped historic recombination. ARs mediate the physiological effects of epinephrine and norepinephrine, and are targets of many therapeutic drugs. This work has identified haplotype backgrounds of alpha1-AR missense variants, haplotype block structures in US Caucasians and African Americans, and haplotype tag SNPs for each block, and we present strong evidence for ancient haplotype block superstructure at these genes which has been partially disrupted by recombination, and evidence for reinstatement of linkage disequilibrium by subsequent recombination events. ADRA1A is comprised of four haplotype blocks in US Caucasians, while in African Americans Block 1 is split. ADRA1B has four blocks in US Caucasians, but in African Americans only the first two blocks are present. ADRA1D has two blocks in US Caucasians, and the first block is replaced by two smaller blocks in African Americans. For both ADRA1A and ADRA1B, haplotype superstructures may represent a novel, higher-level hierarchy in the human genome, which may reduce redundancy of testing by further aggregation of genotype data.Electronic Supplementary Material Supplementary material is available in the online version of this article at Communicated by W. R. McCombie     

Genetic assignment methods for the direct, real-time estimation of migration rate: a simulation-based exploration of accuracy and power

Genetic assignment methods use genotype likelihoods to draw inference about where individuals were or were not born, potentially allowing direct, real-time estimates of dispersal. We used simulated data sets to test the power and accuracy of Monte Carlo resampling methods in generating statistical thresholds for identifying F0 immigrants in populations with ongoing gene flow, and hence for providing direct, real-time estimates of migration rates. The identification of accurate critical values required that resampling methods preserved the linkage disequilibrium deriving from recent generations of immigrants and reflected the sampling variance present in the data set being analysed. A novel Monte Carlo resampling method taking into account these aspects was proposed and its efficiency was evaluated. Power and error were relatively insensitive to the frequency assumed for missing alleles. Power to identify F0 immigrants was improved by using large sample size (up to about 50 individuals) and by sampling all populations from which migrants may have originated. A combination of plotting genotype likelihoods and calculating mean genotype likelihood ratios (DLR) appeared to be an effective way to predict whether F0 immigrants could be identified for a particular pair of populations using a given set of markers.     

Major histocompatibility complex variation at three class II loci in the northern elephant seal

Northern elephant seals were hunted to near extinction in the 19th century, yet have recovered remarkably and now number around 175,000. We surveyed 110 seals for single-strand conformation polymorphism (SSCP) and sequence variation at three major histocompatibility (MHC) class II loci (DQA, DQB and DRB) to evaluate the genetic consequences of the population bottleneck at these loci vs. other well-studied genes. We found very few alleles at each MHC locus, significant variation among breeding sites for the DQA locus, and linkage disequilibrium between the DQB and DRB loci. Northern elephant seals are evidently inbred, although there is as yet no evidence of correlative reductions in fitness.     

Linkage disequilibrium and founder effect analysis of the NF1 gene in French Canadians from the Quebec population

We genotyped 19 neurofibromatosis type 1 (NF1) families from French Canadians of the Quebec population with four intragenic microsatellites (IVS26-2.3, IVS27AC28.4, IVS27AC33.1, and IVS38GT53.0). Linkage analysis of the four microsatellite markers among the 19 NF1 families indicates that the four microsatellites are strongly linked with NF1 disease (LOD = 2.76-3.64). The four markers are associated (P = 0-0.077) except marker pair IVS26-2.3/IVS27AC33.1 (P = 0.18 or 0.17). However, perhaps due to the high mutation rate of the NF1 gene, no founder effect for NF1 was detected in the Quebec French Canadians.     

Variation in recombination rate across the genome: evidence and implications

Recent data from humans and other species provide convincing evidence of variation in recombination rate in different genomic regions. Comparison of physical and genetic maps reveals variation on a scale of megabases, with substantial differences between sexes. Recombination is often suppressed near centromeres and elevated near telomeres, but neither of these observations is true for all chromosomes. In humans, patterns of linkage disequilibrium and experimental measures of recombination from sperm-typing reveal dramatic hotspots of recombination on a scale of kilobases. Genome-wide variation in the amount of crossing-over may be due to variation in the density of hotspots, the intensity of hotspots, or both. Theoretical models of selection and linkage predict that genetic variation will be reduced in regions of low recombination, and this prediction is supported by data from several species. Heterogeneity in rates of crossing-over provides both an opportunity and a challenge for identifying disease genes: as associations occur in blocks, genomic regions containing disease loci may be identified with relatively few markers, yet identifying the causal mutations is unlikely to be achieved through associations alone.     

The effect of single nucleotide polymorphism identification strategies on estimates of linkage disequilibrium

At present there is tremendous interest in characterizing the magnitude and distribution of linkage disequilibrium (LD) throughout the human genome, which will provide the necessary foundation for genome-wide LD analyses and facilitate detailed evolutionary studies. To this end, a human high-density single-nucleotide polymorphism (SNP) marker map has been constructed. Many of the SNPs on this map, however, were identified by sampling a small number of chromosomes from a single population, and inferences drawn from studies using such SNPs may be influenced by ascertainment bias (AB). Through extensive simulations, we have found that AB is a potentially significant problem in estimating and comparing LD within and between populations. Specifically, the magnitude of AB is a function of the SNP discovery strategy, number of chromosomes used for SNP discovery, population genetic characteristics of the particular genomic region considered, amount of gene flow between populations, and demographic history of the populations. We demonstrate that a balanced SNP discovery strategy (where equal numbers of chromosomes are sampled from multiple subpopulations) is the optimal study design for generating broadly applicable SNP resources. Finally, we validate our theoretical predictions by comparing our results to publicly available data from ten genes sequenced in 24 African American and 23 European American individuals.     

Patterns of polymorphism detected in the chloroplast and nuclear genomes of barley landraces sampled from Syria and Jordan

In order to examine how molecular polymorphism in barley landraces, sampled from five different ecogeographical regions of Syria and Jordan, is organised and partitioned, genetic variability at 21 nuclear and 10 chloroplast microsatellite loci were examined. Chloroplast polymorphism was detected, with most variation being ascribed to differences between the five regions (F_st 0.45) and to within sites within each region (F_st 0.44). Moreover, the distribution of chloroplast polymorphism is structured and not distributed randomly across the barley landraces sampled. From a total of 125 landrace accessions (five lines from each of five sites from each of five regions) genotyped with 21 SSRs a total of 244 alleles were detected, of which 38 were common to the five regions sampled. Most nuclear variation was detected within sites. Significant differentiation between sites (F_st 0.29) was detected with nuclear SSRs and this partially mirrored polymorphism in the chloroplast genome. Strong statistical associations/interaction was also detected between the chloroplast and nuclear SSRs, together with non-random association (linkage disequilibrium) of alleles at both linked and unlinked SSR loci. These results are discussed in the context of adaptation of landraces to the extreme environment, the concept of 'adapted gene complexes' and the exploitation of landraces in breeding programmes.Communicated by P. Langridge     

Identification and characterization of coding single-nucleotide polymorphisms within a human olfactory receptor gene cluster

Single-nucleotide polymorphisms (SNPs) were studied in 15 olfactory receptor (OR) coding regions, one control region and two noncoding sequences all residing within a 412 kb OR gene cluster on human chromosome 17p13.3, as well as in other G-protein coupled receptors (GPCRs). A total of 26 SNPs were identified in ORs, 21 of which are coding SNPs (cSNPs). The mean nucleotide diversity of OR coding regions was 0.078% (ranging from 0 to 0.16%), which is about twice higher than that of other GPCRs, and similar to the nucleotide diversity levels of noncoding regions along the human genome. The high polymorphism level in the OR coding regions might be due to a weak positive selection pressure acting on the OR genes. In two cases, OR genes have been found to share the same cSNP. This could be explained by recent gene conversion events, which might be a part of a concerted evolution mechanism acting on the OR superfamily. Using the genotype data of 85 unrelated individuals in 15 SNPs, we found linkage disequilibrium (LD) between pairs of SNPs located on the centromeric part of the cluster. On the other hand, no LD was found between SNPs located on the telomeric part of the cluster, suggesting the presence of several hot-spots for recombination within this cluster. Thus, different regions of this gene cluster may have been subject to different recombination rates.     

Quantitative genetic variation in Daphnia: temporal changes in genetic architecture

Nonadditive genetic variation and genetic disequilibrium are two important factors that influence the evolutionary trajectory of natural populations. We assayed quantitative genetic variation in a temporary-pond-dwelling population of Daphnia pulex over a full season to examine the role of nonadditive genetic variation and genetic disequilibrium in determining the short-term evolutionary trajectory of a cyclic parthenogen. Quantitative traits were influenced by three factors: (1) clonal selection significantly changed the population mean phenotype during the course of the growing season; (2) sexual reproduction and recombination led to significant changes in life-history trait means and the levels of expressed genetic variation, implying the presence of substantial nonadditive genetic variation and genetic disequilibrium; and (3) Egg-bank effects were found to be an important component of the realized year-to-year change. Additionally, we examined the impact of genetic disequilibria induced by clonal selection on the genetic (co)variance structure with a common principal components model. Clonal selection caused significant changes in the (co)variance structure that were eliminated by a single bout of random mating, suggesting that a build-up of disequilibria was the primary source of changes in the (co)variance structure. The results of this study highlight the complexity of natural selection operating on populations that undergo alternating phases of sexual and asexual reproduction.     

The D' measure of overall gametic disequilibrium between pairs of multiallelic loci

Abstract The D ' coefficient is one of the most commonly used measures of the extent of gametic disequilibrium between multiallelic loci. It has been suggested that the range of the D ' measure of overall disequilibrium between pairs of multiallelic loci depends on allele frequencies, except under some very restricted conditions. Nevertheless, the problem of dependence of the range of D ' has not been characterized under a wide set of possible polymorphisms. Evaluation of the utility of D ' as a measure of the strength of overall disequilibrium between all possible pairs of alleles at two multiallelic loci requires better knowledge of its range than is currently available. In this work, the conditions of polymorphism under which the range of D ' is frequency independent are given. It is found that the range of D ' is more often independent of allelic frequencies than is commonly thought. Furthermore, the range of D ' undergoes only small fluctuations as a function of the polymorphisms at the loci. Numerical cases and microsatellite data from humans are used for illustration. These observations indicate that the D ' coefficient is a useful tool for the estimation and comparison of the extent of overall disequilibrium across pairs of multiallelic loci.