Insights into recombination from patterns of linkage disequilibrium in humans

An ability to predict levels of linkage disequilibrium (LD) between linked markers would facilitate the design of association studies and help to distinguish between evolutionary models. Unfortunately, levels of LD depend crucially on the rate of recombination, a parameter that is difficult to measure. In humans, rates of genetic exchange between markers megabases apart can be estimated from a comparison of genetic and physical maps; these large-scale estimates can then be interpolated to predict LD at smaller ("local") scales. However, if there is extensive small-scale heterogeneity, as has been recently proposed, local rates of recombination could differ substantially from those averaged over much larger distances. We test this hypothesis by estimating local recombination rates indirectly from patterns of LD in 84 genomic regions surveyed by the SeattleSNPs project in a sample of individuals of European descent and of African-Americans. We find that LD-based estimates are significantly positively correlated with map-based estimates. This implies that large-scale, average rates are informative about local rates of recombination. Conversely, although LD-based estimates are based on a number of simplifying assumptions, it appears that they capture considerable information about the underlying recombination rate or at least about the ordering of regions by recombination rate. Using LD-based estimators, we also find evidence for homologous gene conversion in patterns of polymorphism. However, as we demonstrate by simulation, inferences about gene conversion are unreliable, even with extensive data from homogeneous regions of the genome, and are confounded by genotyping error.     

A Recombination Hotspot in a Schizophrenia-Associated Region of GABRB2

Background

Schizophrenia is a major disorder with complex genetic mechanisms. Earlier, population genetic studies revealed the occurrence of strong positive selection in the GABRB2 gene encoding the β₂ subunit of GABA_A receptors, within a segment of 3,551 bp harboring twenty-nine single nucleotide polymorphisms (SNPs) and containing schizophrenia-associated SNPs and haplotypes.

Methodology/Principal Findings

In the present study, the possible occurrence of recombination in this ‘S1–S29’ segment was assessed. The occurrence of hotspot recombination was indicated by high resolution recombination rate estimation, haplotype diversity, abundance of rare haplotypes, recurrent mutations and torsos in haplotype networks, and experimental haplotyping of somatic and sperm DNA. The sub-segment distribution of relative recombination strength, measured by the ratio of haplotype diversity (H_d) over mutation rate (θ), was indicative of a human specific Alu-Yi6 insertion serving as a central recombining sequence facilitating homologous recombination. Local anomalous DNA conformation attributable to the Alu-Yi6 element, as suggested by enhanced DNase I sensitivity and obstruction to DNA sequencing, could be a contributing factor of the increased sequence diversity. Linkage disequilibrium (LD) analysis yielded prominent low LD points that supported ongoing recombination. LD contrast revealed significant dissimilarity between control and schizophrenic cohorts. Among the large array of inferred haplotypes, H26 and H73 were identified to be protective, and H19 and H81 risk-conferring, toward the development of schizophrenia.

Conclusions/Significance

The co-occurrence of hotspot recombination and positive selection in the S1–S29 segment of GABRB2 has provided a plausible contribution to the molecular genetics mechanisms for schizophrenia. The present findings therefore suggest that genome regions characterized by the co-occurrence of positive selection and hotspot recombination, two interacting factors both affecting genetic diversity, merit close scrutiny with respect to the etiology of common complex disorders.     

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.     

产SEB金黄色葡萄球菌α-溶血毒素基因缺失菌株的构建

构建产肠毒素B(Staphylococcal enterotoxin B ,SEB)的金黄色葡萄球菌α-溶血毒素（α-hemolysin, α-HL)缺失菌株。首先构建用于α-HL基因敲除的同源重组质粒pMHL-α,经金黄色葡萄球菌RN4220修饰后再通过原生质体转入金黄色葡萄球菌SM-01。含重组质粒pMHL-α的金黄色葡萄球菌SM-01在42℃诱导条件下培养多代,最终筛选出α-溶血毒素基因缺失菌株。经序列分析和血平板溶血实验结果证明最终获得产SEB金黄色葡萄球菌α-HL缺失菌株。为野生型金黄色葡萄球菌的体内遗传操作及构建产超抗原药物金黄色葡萄球菌基因工程菌株提供了一定的理论基础和方法。     

Comparison of a high-density genetic linkage map to genome features in the model grass <Emphasis Type="Italic">Brachypodium distachyon</Emphasis>

The small annual grass Brachypodium distachyon (Brachypodium) is rapidly emerging as a powerful model system to study questions unique to the grasses. Many Brachypodium resources have been developed including a whole genome sequence, highly efficient transformation and a large germplasm collection. We developed a genetic linkage map of Brachypodium using single nucleotide polymorphism (SNP) markers and an F₂ mapping population of 476 individuals. SNPs were identified by targeted resequencing of single copy genomic sequences. Using the Illumina GoldenGate Genotyping platform we placed 558 markers into five linkage groups corresponding to the five chromosomes of Brachypodium. The unusually long total genetic map length, 1,598 centiMorgans (cM), indicates that the Brachypodium mapping population has a high recombination rate. By comparing the genetic map to genome features we found that the recombination rate was positively correlated with gene density and negatively correlated with repetitive regions and sites of ancestral chromosome fusions that retained centromeric repeat sequences. A comparison of adjacent genome regions with high versus low recombination rates revealed a positive correlation between interspecific synteny and recombination rate.     

Linkage disequilibrium and inference of ancestral recombination in 538 single-nucleotide polymorphism clusters across the human genome

The prospect of using linkage disequilibrium (LD) for fine-scale mapping in humans has attracted considerable attention, and, during the validation of a set of single-nucleotide polymorphisms (SNPs) for linkage analysis, a set of data for 4,833 SNPs in 538 clusters was produced that provides a rich picture of local attributes of LD across the genome. LD estimates may be biased depending on the means by which SNPs are first identified, and a particular problem of ascertainment bias arises when SNPs identified in small heterogeneous panels are subsequently typed in larger population samples. Understanding and correcting ascertainment bias is essential for a useful quantitative assessment of the landscape of LD across the human genome. Heterogeneity in the population recombination rate, rho=4Nr, along the genome reflects how variable the density of markers will have to be for optimal coverage. We find that ascertainment-corrected rho varies along the genome by more than two orders of magnitude, implying great differences in the recombinational history of different portions of our genome. The distribution of rho is unimodal, and we show that this is compatible with a wide range of mixtures of hotspots in a background of variable recombination rate. Although rho is significantly correlated across the three population samples, some regions of the genome exhibit population-specific spikes or troughs in rho that are too large to be explained by sampling. This result is consistent with differences in the genealogical depth of local genomic regions, a finding that has direct bearing on the design and utility of LD mapping and on the National Institutes of Health HapMap project.     

Modeling linkage disequilibrium and identifying recombination hotspots using single-nucleotide polymorphism data

We introduce a new statistical model for patterns of linkage disequilibrium (LD) among multiple SNPs in a population sample. The model overcomes limitations of existing approaches to understanding, summarizing, and interpreting LD by (i) relating patterns of LD directly to the underlying recombination process; (ii) considering all loci simultaneously, rather than pairwise; (iii) avoiding the assumption that LD necessarily has a "block-like" structure; and (iv) being computationally tractable for huge genomic regions (up to complete chromosomes). We examine in detail one natural application of the model: estimation of underlying recombination rates from population data. Using simulation, we show that in the case where recombination is assumed constant across the region of interest, recombination rate estimates based on our model are competitive with the very best of current available methods. More importantly, we demonstrate, on real and simulated data, the potential of the model to help identify and quantify fine-scale variation in recombination rate from population data. We also outline how the model could be useful in other contexts, such as in the development of more efficient haplotype-based methods for LD mapping.     

SNP identification, linkage disequilibrium, and haplotype analysis for a 200-kb genomic region in a Korean population

Understanding patterns of linkage disequilibrium (LD) across genomes may facilitate association mapping studies to localize genetic variants influencing complex diseases, a recognition that led to the International Haplotype Mapping Project (HapMap). Divergent patterns of haplotype frequency and LD across global populations require that the HapMap database be supplemented with haplotype and LD data from additional populations. We conducted a pilot study of the LD and haplotype structure of a genomic region in a Korean population. A total of 165 SNPs were identified in a 200-kb region of 22q13.2 by direct sequencing. Unphased genotype data were generated for 76 SNPs in 90 unrelated Korean individuals. LD, haplotype diversity, and recombination rates were assessed in this region and compared with the HapMap database. The pattern of LD and haplotype frequencies of Korean samples showed a high degree of similarity with Japanese data. There was a strong correlation between high LD and low recombination frequency in this region. We found considerable similarities in local LD patterns between three Asian populations (Han Chinese, Japanese, and Korean) and the CEPH population. Haplotype frequencies were, however, significantly different between them. Our results should further the understanding of distinctive Korean genomic features and assist in designing appropriate association studies.     

A Novel Computational Method Identifies Intra- and Inter-Species Recombination Events in Staphylococcus aureus and Streptococcus pneumoniae

Advances in high-throughput DNA sequencing technologies have determined an explosion in the number of sequenced bacterial genomes. Comparative sequence analysis frequently reveals evidences of homologous recombination occurring with different mechanisms and rates in different species, but the large-scale use of computational methods to identify recombination events is hampered by their high computational costs. Here, we propose a new method to identify recombination events in large datasets of whole genome sequences. Using a filtering procedure of the gene conservation profiles of a test genome against a panel of strains, this algorithm identifies sets of contiguous genes acquired by homologous recombination. The locations of the recombination breakpoints are determined using a statistical test that is able to account for the differences in the natural rate of evolution between different genes. The algorithm was tested on a dataset of 75 genomes of Staphylococcus aureus and 50 genomes comprising different streptococcal species, and was able to detect intra-species recombination events in S. aureus and in Streptococcus pneumoniae. Furthermore, we found evidences of an inter-species exchange of genetic material between S. pneumoniae and Streptococcus mitis, a closely related commensal species that colonizes the same ecological niche. The method has been implemented in an R package, Reco, which is freely available from supplementary material, and provides a rapid screening tool to investigate recombination on a genome-wide scale from sequence data.     

Population Genetic Analysis of Streptomyces albidoflavus Reveals Habitat Barriers to Homologous Recombination in the Diversification of Streptomycetes

Examining the population structure and the influence of recombination and ecology on microbial populations makes great sense for understanding microbial evolution and speciation. Streptomycetes are a diverse group of bacteria that are widely distributed in nature and a rich source of useful bioactive compounds; however, they are rarely subjected to population genetic investigations. In this study, we applied a five-gene-based multilocus sequence analysis (MLSA) scheme to 41 strains of Streptomyces albidoflavus derived from diverse sources, mainly insects, sea, and soil. Frequent recombination was detected in S. albidoflavus, supported by multiple lines of evidence from the pairwise homoplasy index (Φ_w) test, phylogenetic discordance, the Shimodaira-Hasegawa (SH) test, and network analysis, underpinning the predominance of homologous recombination within Streptomyces species. A strong habitat signal was also observed in both phylogenetic and Structure 2.3.3 analyses, indicating the importance of ecological difference in shaping the population structure. Moreover, all three habitat-associated groups, particularly the entomic group, demonstrated significantly reduced levels of gene flow with one another, generally revealing habitat barriers to recombination. Therefore, a combined effect of homologous recombination and ecology is inferred for S. albidoflavus, where dynamic evolution is at least partly balanced by the extent that differential distributions of strains among habitats limit genetic exchange. Our study stresses the significance of ecology in microbial speciation and reveals the coexistence of homologous recombination and ecological divergence in the evolution of streptomycetes.     

Demography,recombination hotspot intensity,and the block structure of linkage disequilibrium

BACKGROUND: Effective gene mapping based on genetic association data will require detailed knowledge of patterns of linkage disequilibrium (LD) in human populations. It has been recently suggested that linkage disequilibrium in humans may be organized in a block-like structure, with islands of high LD separated by regions of rapid breakdown of LD due to recombination hotspots. The experimental data to date, however, are limited, and fundamental questions remain about the implications of recombination rate heterogeneity. Here, we use computer simulations to evaluate how such heterogeneity influences patterns of LD, and we develop formal criteria to assess whether the patterns are functionally block like in the context of association mapping.RESULTS: Our analyses suggest that, even in models of extreme recombination rate heterogeneity, some human populations will have a functionally block-like structure to the pattern of LD, but others will not, depending on their precise demographic histories. In fact, for many models, we find that, following an LD-generating event, populations may move through discrete phases that can be functionally described as pre-block, block, and post-block. An analysis of observed and expected patterns of LD surrounding hotspots within the MHC Class II region confirms these theoretical expectations.CONCLUSIONS: Even if highly punctuated patterns of recombination are the rule, patterns of LD are still likely to show differences among populations and among genomic regions that are of practical importance in the design of genetic association studies. The notion that the average extent of LD is a useful concept for the design of association studies must be abandoned in light of the experimental and theoretical evidence.     

Population demographic history can cause the appearance of recombination hotspots

Although the prevailing view among geneticists suggests that recombination hotspots exist ubiquitously across the human genome, there is only limited experimental evidence from a few genomic regions to support the generality of this claim. A small number of true recombination hotspots are well supported experimentally, but the vast majority of hotspots have been identified on the basis of population genetic inferences from the patterns of linkage disequilibrium (LD) seen in the human population. These inferences are made assuming a particular model of human history, and one of the assumptions of that model is that the effective population size of humans has remained constant throughout our history. Our results show that relaxation of the constant population size assumption can create LD and variation patterns that are qualitatively and quantitatively similar to human populations without any need to invoke localized hotspots of recombination. In other words, apparent recombination hotspots could be an artifact of variable population size over time. Several lines of evidence suggest that the vast majority of hotspots identified on the basis of LD information are unlikely to have elevated recombination rates.     

Equilibrium processes cannot explain high levels of short- and medium-range linkage disequilibrium in the domesticated grass Sorghum bicolor

Patterns of linkage disequilibrium (LD) are of interest because they provide evidence of both equilibrium (e.g., mating system or long-term population structure) and nonequilibrium (e.g., demographic or selective) processes, as well as because of their importance in strategies for identifying the genetic basis of complex phenotypes. We report patterns of short and medium range (up to 100 kb) LD in six unlinked genomic regions in the partially selfing domesticated grass, Sorghum bicolor. The extent of allelic associations in S. bicolor, as assessed by pairwise measures of LD, is higher than in maize but lower than in Arabidopsis, in qualitative agreement with expectations based on mating system. Quantitative analyses of the population recombination parameter, rho, however, based on empirical estimates of rates of recombination, mutation, and self-pollination, show that LD is more extensive than expected under a neutral equilibrium model. The disparity between rho and the population mutation parameter, , is similar to that observed in other species whose population history appears to be complex. From a practical standpoint, these results suggest that S. bicolor is well suited for association studies using reasonable numbers of markers, since LD typically extends at least several kilobases but has largely decayed by 15 kb.     

Lower-than-expected linkage disequilibrium between tightly linked markers in humans suggests a role for gene conversion

Understanding the pattern of linkage disequilibrium (LD) in the human genome is important both for successful implementation of disease-gene mapping approaches and for inferences about human demographic histories. Previous studies have examined LD between loci within single genes or confined genomic regions, which may not be representative of the genome; between loci separated by large distances, where little LD is seen; or in population groups that differ from one study to the next. We measured LD in a large set of locus pairs distributed throughout the genome, with loci within each pair separated by short distances (average 124 bp). Given current models of the history of the human population, nearly all pairs of loci at such short distances would be expected to show complete LD as a consequence of lack of recombination in the short interval. Contrary to this expectation, a significant fraction of pairs showed incomplete LD. A standard model of recombination applied to these data leads to an estimate of effective human population size of 110,000. This estimate is an order of magnitude higher than most estimates based on nucleotide diversity. The most likely explanation of this discrepancy is that gene conversion increases the apparent rate of recombination between nearby loci.     

Fine-scale map of encyclopedia of DNA elements regions in the Korean population

The International HapMap Project aims to generate detailed human genome variation maps by densely genotyping single-nucleotide polymorphisms (SNPs) in CEPH, Chinese, Japanese, and Yoruba samples. This will undoubtedly become an important facility for genetic studies of diseases and complex traits in the four populations. To address how the genetic information contained in such variation maps is transferable to other populations, the Korean government, industries, and academics have launched the Korean HapMap project to genotype high-density Encyclopedia of DNA Elements (ENCODE) regions in 90 Korean individuals. Here we show that the LD pattern, block structure, haplotype diversity, and recombination rate are highly concordant between Korean and the two HapMap Asian samples, particularly Japanese. The availability of information from both Chinese and Japanese samples helps to predict more accurately the possible performance of HapMap markers in Korean disease-gene studies. Tagging SNPs selected from the two HapMap Asian maps, especially the Japanese map, were shown to be very effective for Korean samples. These results demonstrate that the HapMap variation maps are robust in related populations and will serve as an important resource for the studies of the Korean population in particular.     

Evidence for clonal evolution among highly polymorphic genes in methicillin-resistant Staphylococcus aureus

The evolution of Staphylococcus aureus has been described as predominantly clonal, based on evidence from seven housekeeping genes. We aimed to test if this was also true for more polymorphic genes. In a collection of 60 isolates including major European epidemic methicillin-resistant S. aureus (MRSA) and sporadic MRSA strains, we compared the partial gene sequences of seven housekeeping genes (arcC, aroE, glpF, gmk, pta, tpi, and yqiL), six core adhesion genes (present in all strains) (clfA, clfB, fnbA, map, sdrC, and spa), and four accessory adhesion genes (not present in all strains) (ebpS, fnbB, sdrD, and sdrE). Nucleotide diversity of adhesion genes was 2- to 10-fold higher than genes used for multilocus sequence typing. All genes showed evidence for purifying selection with a weakly reduced level among accessory adhesion genes. Among these highly variable genes, there was no evidence for a difference in molecular evolution between epidemic and sporadic strains. Gene trees constructed from concatenated sequences of housekeeping, core adhesion, and accessory adhesion genes were highly congruent, indicating clonality, despite some evidence for homologous exchange. Further evidence for clonality was found with an overall positive correlation of allelic and nucleotidic divergence for both seven housekeeping genes and six core adhesion genes. However, for small allelic differences that fit the demarcations of clonal complexes (CCs) there was no such correlation, suggesting that recombination occurred. Therefore, despite an overall clonal population structure, recombination between related isolates within CCs might have contributed to S. aureus evolution.     

Auxotrophic mutant of Staphylococcus aureus interferes with nasal colonization by the wild type

Staphylococcus aureus nasal carriage is a risk factor for infection in humans, particularly in the hospital setting. Bacterial interference was used as an alternative strategy for the prevention of upper respiratory, urogenital and gastrointestinal tract infections. This study was designed to assess if the administration of a live-attenuated aroA mutant of S. aureus is useful as a potential approach to prevent transient staphylococcal nasal carriage by virulent strains. We constructed an aroA mutant of S. aureus Newman strain by homologous recombination. The auxotrophic NK41 mutant was attenuated as determined by the increase of the LD(50) after intraperitoneal challenge. In mice, previous nasal colonization with the NK41 mutant significantly reduced the number of CFU of S. aureus (HU-71 and Hde288) clinical isolates and the parental Newman strain. The NK41 mutant was unable to induce a pro-inflammatory response and to damage the invaded human respiratory epithelial cells. Moreover, the cells previously or simultaneously infected with the NK41 mutant were invaded by virulent strains in a significantly lower degree than those of the control group. In conclusion, the attenuated NK41 mutant interfered with the colonization and establishment of pathogenic strains of S. aureus, which produce severe infections.     

Nonhomologous end-joining deficiency allows large chromosomal deletions to be produced by replacement-type recombination in Aspergillus oryzae

Gene targeting is a technique of introducing a genetic trait at a predetermined site within a genome; it is also used to eliminate undesirable chromosomal regions from the relevant genome. Thus far, replacement-type recombination between two homologous regions separated by a large nonhomologous sequence has been hardly achieved probably due to the low frequency of homologous recombination in filamentous fungi. In this study, we report the successful and highly efficient deletion by replacement-type recombination of up to 470-kb regions of chromosome 8 and 200-kb region in chromosome 3, which includes a homologue of aflatoxin gene cluster, by nonhomologous end-joining deficient strains of Aspergillus oryzae. Our study results indicate that the deficiency of nonhomologous end-joining increases the distance of nonhomologous regions in replacement-type recombination, i.e., the possible deletion range in generation of large chromosomal deletion by one cycle of replacement-type recombination is increased in nonhomologous end-joining deficient strains.     

Spectrum: joint Bayesian inference of population structure and recombination events

MOTIVATION: While genetic properties such as linkage disequilibrium (LD) and population structure are closely related under a common inheritance process, the statistical methodologies developed so far mostly deal with LD analysis and structural inference separately, using specialized models that do not capture their statistical and genetic relationships. Also, most of these approaches ignore the inherent uncertainty in the genetic complexity of the data and rely on inflexible models built on a closed genetic space. These limitations may make it difficult to infer detailed and consistent structural information from rich genomic data such as populational single nucleotide polymorphisms (SNP) profiles. RESULTS: We propose a new model-based approach to address these issues through joint inference of population structure and recombination events under a non-parametric Bayesian framework; we present Spectrum, an efficient implementation based on our new model. We validated Spectrum on simulated data and applied it to two real SNP datasets, including single-population Daly data and the four-population HapMap data. Our method performs well relative to LDhat 2.0 in estimating the recombination rates and hotspots on these datasets. More interestingly, it generates an ancestral spectrum for representing population structures which not only displays sub-structure based on population founders but also reveals details of the genetic diversity of each individual. It offers an alternative view of the population structures to that offered by Structure 2.1, which ignores chromosome-level mutation and recombination with respect to founders.     

Haplotypes at ATM identify coding-sequence variation and indicate a region of extensive linkage disequilibrium

Genetic variation in the human population may lead to functional variants of genes that contribute to risk for common chronic diseases such as cancer. In an effort to detect such possible predisposing variants, we constructed haplotypes for a candidate gene and tested their efficacy in association studies. We developed haplotypes consisting of 14 biallelic neutral-sequence variants that span 142 kb of the ATM locus. ATM is the gene responsible for the autosomal recessive disease ataxia-telangiectasia (AT). These ATM noncoding single-nucleotide polymorphisms (SNPs) were genotyped in nine CEPH families (89 individuals) and in 260 DNA samples from four different ethnic origins. Analysis of these data with an expectation-maximization algorithm revealed 22 haplotypes at this locus, with three major haplotypes having frequencies > or = .10. Tests for recombination and linkage disequilibrium (LD) show reduced recombination and extensive LD at the ATM locus, in all four ethnic groups studied. The most striking example was found in the study population of European ancestry, in which no evidence for recombination could be discerned. The potential of ATM haplotypes for detection of genetic variants through association studies was tested by analysis of 84 individuals carrying one of three ATM coding SNPs. Each coding SNP was detected by association with an ATM haplotype. We demonstrate that association studies with haplotypes for candidate genes have significant potential for the detection of genetic backgrounds that contribute to disease.