Diogo Meyer

How selection shapes variation of the human major histocompatibility complex: a review

The nature of polymorphism and molecular sequence variation in the genes of the human major histocompatibility complex (MHC) provides strong support for the idea that these genes are under selection. With the understanding that selection shapes MHC variation new questions have become the focus of study. What is the mode of selection that accounts for MHC polymorphism? Is variation maintained by pathogen pressure or by reproductive mechanisms? Discerning between these requires drawing on information from studies on association between HLA genes and infectious diseases, reproductive success and mating preferences relative to HLA genotypes, and theoretical studies that compare the outcomes of different selection regimes. The pattern that has emerged suggests that several types of selection are plausible for the maintenance of HLA polymorphism.

Global diversity and evidence for coevolution of KIR and HLA.

The killer immunoglobulin-like receptor (KIR) gene cluster shows extensive genetic diversity, as do the HLA class I loci, which encode ligands for KIR molecules. We genotyped 1,642 individuals from 30 geographically distinct populations to examine population-level evidence for coevolution of these two functionally related but unlinked gene clusters. We observed strong negative correlations between the presence of activating KIR genes and their corresponding HLA ligand groups across populations, especially KIR3DS1 and its putative HLA-B Bw4-80I ligands (r = −0.66, P = 0.038). In contrast, we observed weak positive relationships between the various inhibitory KIR genes and their ligands. We observed a negative correlation between distance from East Africa and frequency of activating KIR genes and their corresponding ligands, suggesting a balance between selection on HLA and KIR loci. Most KIR-HLA genetic association studies indicate a primary influence of activating KIR-HLA genotypes in disease risk; concomitantly, activating receptor-ligand pairs in this study show the strongest signature of coevolution of these two complex genetic systems as compared with inhibitory receptor-ligand pairs.

PyPop: a software framework for population genomics: analyzing large-scale multi-locus genotype data.

Software to analyze multi-locus genotype data for entire populations is useful for estimating haplotype frequencies, deviation from Hardy-Weinberg equilibrium and patterns of linkage disequilibrium. These statistical results are important to both those interested in human genome variation and disease predisposition as well as evolutionary genetics. As part of the 13th International Histocompatibility and Immunogenetics Working Group (IHWG), we have developed a software framework (PyPop). The primary novelty of this package is that it allows integration of statistics across large numbers of data-sets by heavily utilizing the XML file format and the R statistical package to view graphical output, while retaining the ability to inter-operate with existing software. Largely developed to address human population data, it can, however, be used for population based data for any organism. We tested our software on the data from the 13th IHWG which involved data sets from at least 50 laboratories each of up to 1000 individuals with 9 MHC loci (both class I and class II) and found that it scales to large numbers of data sets well.

Signatures of demographic history and natural selection in the human major histocompatibility complex Loci.

Many lines of evidence show that several HLA loci have experienced balancing selection. However, distinguishing among demographic and selective explanations for patterns of variation observed with HLA genes remains a challenge. In this study we address this issue using data from a diverse set of human populations at six classical HLA loci and, employing a comparative genomics approach, contrast results for HLA loci to those for non-HLA markers. Using a variety of analytic methods, we confirm and extend evidence for selection acting on several HLA loci. We find that allele frequency distributions for four of the six HLA loci deviate from neutral expectations and show that this is unlikely to be explained solely by demographic factors. Other features of HLA variation are explained in part by demographic history, including decreased heterozygosity and increased LD for populations at greater distances from Africa and a similar apportionment of genetic variation for HLA loci compared to putatively neutral non-HLA loci. On the basis of contrasts among different HLA loci and between HLA and non-HLA loci, we conclude that HLA loci bear detectable signatures of both natural selection and demographic history.

Hybridization and introgression across different ploidy levels in the Neotropical orchids Epidendrum fulgens and E. puniceoluteum (Orchidaceae)

The hypothesis of gene flow between species with large differences in chromosome numbers has rarely been tested in the wild, mainly because species of different ploidy are commonly assumed to be reproductively isolated from each other because of instantaneous and strong postzygotic barriers. In this study, a broad‐scale survey of molecular variation was carried out between two orchid species with different ploidy levels: Epidendrum fulgens (2n = 2x = 24 chromosomes) and Epidendrum puniceoluteum (2n = 4x = 52 chromosomes). To test the strength of their reproductive barriers, we investigated the distribution of genetic variation in sympatric and allopatric populations of these two species and conducted crossing experiments. Nuclear and plastid microsatellite loci were used to genotype 463 individuals from eight populations across the geographical range of both species along the Brazilian coastal plain. All six sympatric populations analysed presented hybrid zones, indicating that hybridization between E. fulgens and E. puniceoluteum is a common phenomenon. Bayesian assignment analysis detected the presence of F1 and F2 individuals and also signs of introgression, demonstrating a high potential for interspecific gene flow. Introgression occurs preferentially from E. fulgens to E. puniceoluteum. Pure parental individuals of both species display strong genotype–habitat associations, indicating that environment‐dependent selection could be acting in all hybrid zones. This study suggests that hybridization and introgression are evolutionary processes playing a role in the diversification of Epidendrum and indicates the importance of investigations of hybrid zones in understanding reproductive barriers and speciation processes in Neotropical orchid species.

The genetic effects of Late Quaternary climatic changes over a tropical latitudinal gradient: diversification of an Atlantic Forest passerine

The increase in biodiversity from high to low latitudes is a widely recognized biogeographical pattern. According to the latitudinal gradient hypothesis (LGH), this pattern was shaped by differential effects of Late Quaternary climatic changes across a latitudinal gradient. Here, we evaluate the effects of climatic changes across a tropical latitudinal gradient and its implications to diversification of an Atlantic Forest (AF) endemic passerine. We studied the intraspecific diversification and historical demography of Sclerurus scansor, based on mitochondrial (ND2, ND3 and cytb) and nuclear (FIB7) gene sequences. Phylogenetic analyses recovered three well‐supported clades associated with distinct latitudinal zones. Coalescent‐based methods were applied to estimate divergence times and changes in effective population sizes. Estimates of divergence times indicate that intraspecific diversification took place during Middle–Late Pleistocene. Distinct demographic scenarios were identified, with the southern lineage exhibiting a clear signature of demographic expansion, while the central one remained more stable. The northern lineage, contrasting with LGH predictions, exhibited a clear sign of a recent bottleneck. Our results suggest that different AF regions reacted distinctly, even in opposite ways, under the same climatic period, producing simultaneously favourable scenarios for isolation and contact among populations.

Methods for assessing gene content diversity of KIR with examples from a global set of populations

A number of statistical methods are widely used to describe allelic variation at specific genetic loci and its implication on the evolutionary history of these loci. Although the methods were developed primarily to study allelic variation at loci that are virtually always present in the genome, they are often applied to data of gene content variation (i.e., presence/absence of multiple homologous genes) at the killer cell immunoglobulin-like receptor (KIR) gene cluster. In this paper, we discuss methodological issues involved in the analysis of gene content variation data in the KIR region and also its covariation with polymorphism at the human leukocyte antigen class I loci, which encode ligands for KIR. A comparison of several statistical methods and measures (gene frequency, haplotype frequency, and linkage disequilibrium estimation) using the Centre d’Etude du Polymorphisme Humain data will be provided using KIR haplotypes that have been determined by segregation analysis, noting the strengths and weaknesses of the methods when only the presence/absence data is considered. Finally, application of these methods to a set of globally distributed populations is described (see Single et al., Nat Genet 39:1114–1119, 2007) in order to illustrate the challenges faced when inferring the joint effects of natural selection and demographic history on these immune-related genes.

Genetic diversity of two African and sixteen South American populations determined on the basis of six hypervariable loci.

A total of 582 individuals (1,164 chromosomes) from two African, eight African-derived South American, five South American Amerindian, and three Brazilian urban populations were studied at four variable number of tandem repeat (VNTR) and two short tandem repeat (STR) hypervariable loci. These two sets of loci did not show distinct allele profiles, which might be expected if different processes promoted their molecular differentiation. The two African groups showed little difference between them, and their intrapopulational variation was similar to those obtained in the African-derived South American communities. The latter showed different degrees of interpopulation variability, despite the fact that they presented almost identical average degrees of non-African admixture. The F(ST) single locus estimates differed in the five sets of populations, probably due to genetic drift, indicating the need to consider population structure in the evaluation of their total variability. A high interpopulational diversity was found among Amerindian populations in relation to Brazilian African-derived isolated communities. This is probably a consequence of the differences in the patterns of gene flow and genetic drift that each of these semi-isolated groups experienced.

Region 8q24 is a susceptibility locus for nonsyndromic oral clefting in Brazil

BACKGROUND Nonsyndromic cleft lip with or without cleft palate is a relatively common craniofacial defect with multifactorial inheritance. The association of the rs987525 single nucleotide variant, located in a gene desert at 8q24.21 region, has been consistently replicated in European populations. We performed a structured association approach combined with transcriptional analysis of the MYC gene to dissect the role of rs987525 in oral clefting susceptibility in the ethnically admixed Brazilian population. METHODS We performed the association study conditioned on the individual ancestry proportions in a sample of 563 patients and 336 controls, and in an independent sample of 221 patients and 261 controls. The correlation between rs987525 genotypes and MYC transcriptional levels in orbicularis oris muscle mesenchymal stem cells was also investigated in 42 patients and 4 controls. RESULTS We found a significant association in the larger sample (p = 0.0016; OR = 1.80 [95% confidence interval {CI}, 1.21-2.69], for heterozygous genotype, and 2.71 [95% CI, 1.47-4.96] for homozygous genotype). We did not find a significant correlation between rs987525 genotypes and MYC transcriptional levels (p = 0.14; r = -0.22, Spearman Correlation). CONCLUSIONS We present a positive association of rs987525 in the Brazilian population for the first time, and it is likely that the European contribution to our population is driving this association. We also cannot discard a role of rs987515 in MYC regulation, because this locus behaves as an expression quantitative locus of MYC in another tissue.

Phylogenetic nomenclature and evolution of mannose-binding lectin (MBL2) haplotypes

BackgroundPolymorphisms of the mannose-binding lectin gene (MBL2) affect the concentration and functional efficiency of the protein. We recently used haplotype-specific sequencing to identify 23 MBL2 haplotypes, associated with enhanced susceptibility to several diseases.ResultsIn this work, we applied the same method in 288 and 470 chromosomes from Gabonese and European adults, respectively, and found three new haplotypes in the last group. We propose a phylogenetic nomenclature to standardize MBL2 studies and found two major phylogenetic branches due to six strongly linked polymorphisms associated with high MBL production. They presented high Fst values and were imbedded in regions with high nucleotide diversity and significant Tajimas D values. Compared to others using small sample sizes and unphased genotypic data, we found differences in haplotyping, frequency estimation, Fu and Lis D* and Fst results.ConclusionUsing extensive testing for selective neutrality, we confirmed that stochastic evolutionary factors have had a major role in shaping this polymorphic gene worldwide.