Genetic Variation and Population Structure in Native Americans

We examined genetic diversity and population structure in the American landmass using 678 autosomal microsatellite markers genotyped in 422 individuals representing 24 Native American populations sampled from North, Central, and South America. These data were analyzed jointly with similar data available in 54 other indigenous populations worldwide, including an additional five Native American groups. The Native American populations have lower genetic diversity and greater differentiation than populations from other continental regions. We observe gradients both of decreasing genetic diversity as a function of geographic distance from the Bering Strait and of decreasing genetic similarity to Siberians-signals of the southward dispersal of human populations from the northwestern tip of the Americas. We also observe evidence of: (1) a higher level of diversity and lower level of population structure in western South America compared to eastern South America, (2) a relative lack of differentiation between Mesoamerican and Andean populations, (3) a scenario in which coastal routes were easier for migrating peoples to traverse in comparison with inland routes, and (4) a partial agreement on a local scale between genetic similarity and the linguistic classification of populations. These findings offer new insights into the process of population dispersal and differentiation during the peopling of the Americas.

Genetic variation and demographic history of the Haplochromis laparogramma group of Lake Victoria-An analysis based on SINEs and mitochondrial DNA

More than 500 endemic haplochromine cichlid species inhabit Lake Victoria. This striking species diversity is a classical example of recent explosive adaptive radiation thought to have happened within the last similar to 15,000 years. In this study, we examined the population structure and historical demography of 3 pelagic haplochromine cichlid species that resemble in morphology and have similar niche, Haplochromis (Yssichromis) laparogramma, Haplochromis (Y.) pyrrhocephalus, and Haplochromis (Y.) sp. "glaucocephalus". We investigated the sequences of the mitochondrial DNA control region and the insertion patterns of short interspersed elements (SINEs) of 759 individuals. We show that sympatric forms are genetically differentiated in 4 of 6 cases, but we also found apparent weakening of the genetic differentiation in areas with turbid water. We estimated the timings of population expansion and species divergence to coincide with the refilling of the lake at the Pleistocene/Holocene boundary. We also found that estimates can be altered significantly by the choice of the shape of the molecular clock. If we employ the nonlinear clock model of evolutionary rates in which the rates are higher towards the recent, the population expansion was dated at around the event of desiccation of the lake ca. 17,000 YBP. Thus, we succeeded in clarifying the species and population structure of closely related Lake Victoria cichlids and in showing the importance of applying appropriate clock calibrations in elucidating recent evolutionary events. (C) 2009 Elsevier B.V. All rights reserved.

A genome-to-genome analysis of associations between human genetic variation, HIV-1 sequence diversity, and viral control

HIV-1 sequence diversity is affected by selection pressures arising from host genomic factors. Using paired human and viral data from 1071 individuals, we ran >3000 genome-wide scans, testing for associations between host DNA polymorphisms, HIV-1 sequence variation and plasma viral load (VL), while considering human and viral population structure. We observed significant human SNP associations to a total of 48 HIV-1 amino acid variants (p<2.4 × 10−12). All associated SNPs mapped to the HLA class I region. Clinical relevance of host and pathogen variation was assessed using VL results. We identified two critical advantages to the use of viral variation for identifying host factors: (1) association signals are much stronger for HIV-1 sequence variants than VL, reflecting the ‘intermediate phenotype’ nature of viral variation; (2) association testing can be run without any clinical data. The proposed genome-to-genome approach highlights sites of genomic conflict and is a strategy generally applicable to studies of host–pathogen interaction.

Origins of Shared Genetic Variation in African Cichlids

Cichlid fishes have evolved tremendous morphological and behavioral diversity in the waters of East Africa. Within each of the Great Lakes Tanganyika, Malawi, and Victoria, the phenomena of hybridization and retention of ancestral polymorphism explain allele sharing across species. Here, we explore the sharing of single nucleotide polymorphisms (SNPs) between the major East African cichlid assemblages. A set of approximately 200 genic and nongenic SNPs was ascertained in five Lake Malawi species and genotyped in a diverse collection of 160 species from across Africa. We observed segregating polymorphism outside of the Malawi lineage for more than 50% of these loci; this holds similarly for genic versus nongenic SNPs, as well as for SNPs at putative CpG versus non-CpG sites. Bayesian and principal component analyses of genetic structure in the data demonstrate that the Lake Malawi endemic flock is not monophyletic and that river species have likely contributed significantly to Malawi genomes. Coalescent simulations support the hypothesis that river cichlids have transported polymorphism between lake assemblages. We observed strong genetic differentiation between Malawi lineages for approximately 8% of loci, with contributions from both genic and nongenic SNPs. Notably, more than half of these outlier loci between Malawi groups are polymorphic outside of the lake. Cichlid fishes have evolved diversity in Lake Malawi as new mutations combined with standing genetic variation shared across East Africa.

Population amalgamation and genetic variation: observations on artificially agglomerated tribal populations of Central and South America.

The interpretation of data on genetic variation with regard to the relative roles of different evolutionary factors that produce and maintain genetic variation depends critically on our assumptions concerning effective population size and the level of migration between neighboring populations. In humans, recent population growth and movements of specific ethnic groups across wide geographic areas mean that any theory based on assumptions of constant population size and absence of substructure is generally untenable. We examine the effects of population subdivision on the pattern of protein genetic variation in a total sample drawn from an artificial agglomerate of 12 tribal populations of Central and South America, analyzing the pooled sample as though it were a single population. Several striking findings emerge. (1) Mean heterozygosity is not sensitive to agglomeration, but the number of different alleles (allele count) is inflated, relative to neutral mutation/drift/equilibrium expectation. (2) The inflation is most serious for rare alleles, especially those which originally occurred as tribally restricted "private" polymorphisms. (3) The degree of inflation is an increasing function of both the number of populations encompassed by the sample and of the genetic divergence among them. (4) Treating an agglomerated population as though it were a panmictic unit of long standing can lead to serious biases in estimates of mutation rates, selection pressures, and effective population sizes. Current DNA studies indicate the presence of numerous genetic variants in human populations. The findings and conclusions of this paper are all fully applicable to the study of genetic variation at the DNA level as well.

Genetic variation in genes for the xenobiotic-metabolizing enzymes CYP1A1, EPHX1, GSTM1, GSTT1, and GSTP1 and susceptibility to colorectal cancer in Lynch syndrome.

Individuals with Lynch syndrome are predisposed to cancer due to an inherited DNA mismatch repair gene mutation. However, there is significant variability observed in disease expression likely due to the influence of other environmental, lifestyle, or genetic factors. Polymorphisms in genes encoding xenobiotic-metabolizing enzymes may modify cancer risk by influencing the metabolism and clearance of potential carcinogens from the body. In this retrospective analysis, we examined key candidate gene polymorphisms in CYP1A1, EPHX1, GSTT1, GSTM1, and GSTP1 as modifiers of age at onset of colorectal cancer among 257 individuals with Lynch syndrome. We found that subjects heterozygous for CYP1A1 I462V (c.1384A>G) developed colorectal cancer 4 years earlier than those with the homozygous wild-type genotype (median ages, 39 and 43 years, respectively; log-rank test P = 0.018). Furthermore, being heterozygous for the CYP1A1 polymorphisms, I462V and Msp1 (g.6235T>C), was associated with an increased risk for developing colorectal cancer [adjusted hazard ratio for AG relative to AA, 1.78; 95% confidence interval, 1.16-2.74; P = 0.008; hazard ratio for TC relative to TT, 1.53; 95% confidence interval, 1.06-2.22; P = 0.02]. Because homozygous variants for both CYP1A1 polymorphisms were rare, risk estimates were imprecise. None of the other gene polymorphisms examined were associated with an earlier onset age for colorectal cancer. Our results suggest that the I462V and Msp1 polymorphisms in CYP1A1 may be an additional susceptibility factor for disease expression in Lynch syndrome because they modify the age of colorectal cancer onset by up to 4 years.

Nuclear markers reveal unexpected genetic variation and a Congolese-Nilotic origin of the Lake Victoria cichlid species flock

Phylogenetic analyses based on mitochondrial (mt) DNA have indicated that the cichlid species flock of the Lake Victoria region is derived from a single ancestral species found in East African rivers, closely related to the ancestor of the Lake Malawi cichlid species flock. The Lake Victoria flock contains ten times less mtDNA variation than the Lake Malawi radiation, consistent with current estimates of the ages of the lakes. We present results of a phylogenetic investigation using nuclear (amplified fragment length polymorphism) markers and a wider coverage of riverine haplochromines. We demonstrate that the Lake Victoria–Edward flock is derived from the morphologically and ecologically diverse cichlid genus Thoracochromis from the Congo and Nile, rather than from the phenotypically conservative East African Astatotilapia. This implies that the ability to express much of the morphological diversity found in the species flock may by far pre–date the origin of the flock. Our data indicate that the nuclear diversity of the Lake Victoria–Edward species flock is similar to that of the Lake Malawi flock, indicating that the genetic diversity is considerably older than the 15 000 years that have passed since the lake began to refill. Most of this variation is manifested in trans–species polymorphisms, indicating very recent cladogenesis from a genetically very diverse founder stock. Our data do not confirm strict monophyly of either of the species flocks, but raise the possibility that these flocks have arisen from hybrid swarms.

The interplay between host genetic variation, viral replication and microbial translocation in untreated HIV-infected individuals.

Systemic immune activation, a major determinant of HIV disease progression, is the result of a complex interplay between viral replication, dysregulation of the immune system, and microbial translocation due to gut mucosal damage. While human genetic variants influencing HIV viral load have been identified, it is unknown to what extent the host genetic background contributes to inter-individual differences in other determinants of HIV pathogenesis like gut damage and microbial translocation. Using samples and data from 717 untreated participants in the Swiss HIV Cohort Study and a genome-wide association study design, we searched for human genetic determinants of plasma levels of intestinal fatty-acid binding protein (I-FABP/FABP2), a marker of gut damage, and of soluble sCD14 (sCD14), a marker of LPS bioactivity and microbial translocation. We also assessed the correlations between HIV viral load, sCD14 and I-FABP. While we found no genome-wide significant determinant of the tested plasma markers, we observed strong associations between sCD14 and both HIV viral load and I-FABP, shedding new light on the relationships between processes that drive progression of untreated HIV infection.

Identification, quantification, spatiotemporal distribution and genetic variation of major latex secondary metabolites in the common dandelion (Taraxacum officinale agg.)

The secondary metabolites in the roots, leaves and flowers of the common dandelion (Taraxacum officinale agg.) have been studied in detail. However, little is known about the specific constituents of the plant’s highly specialized laticifer cells. Using a combination of liquid and gas chromatography, mass spectrometry and nuclear magnetic resonance spectrometry, we identified and quantified the major secondary metabolites in the latex of different organs across different growth stages in three genotypes, and tested the activity of the metabolites against the generalist root herbivore Diabrotica balteata. We found that common dandelion latex is dominated by three classes of secondary metabolites: phenolic inositol esters (PIEs), triterpene acetates (TritAc) and the sesquiterpene lactone taraxinic acid β-d-glucopyranosyl ester (TA-G). Purification and absolute quantification revealed concentrations in the upper mg g−1 range for all compound classes with up to 6% PIEs, 5% TritAc and 7% TA-G per gram latex fresh weight. Contrary to typical secondary metabolite patterns, concentrations of all three classes increased with plant age. The highest concentrations were measured in the main root. PIE profiles differed both quantitatively and qualitatively between plant genotypes, whereas TritAc and TA-G differed only quantitatively. Metabolite concentrations were positively correlated within and between the different compound classes, indicating tight biosynthetic co-regulation. Latex metabolite extracts strongly repelled D. balteata larvae, suggesting that the latex constituents are biologically active.

Influence of CYP3A5 genetic variation on everolimus maintenance dosing after cardiac transplantation

BACKGROUND Everolimus (ERL) has become an alternative to calcineurin inhibitors (CNIs) due to its renal-sparing properties, especially in heart transplant (HTx) recipients with kidney dysfunction. However, ERL dosing is challenging due to its narrow therapeutic window combined with high inter-individual pharmacokinetic variability. Our aim was to evaluate the effect of clinical and genetic factors on ERL dosing in a pilot cohort of 37 HTx recipients. METHODS Variants in CYP3A5, CYP3A4, CYP2C8, POR, NR1I2, and ABCB1 were genotyped and clinical data were retrieved from patient charts. RESULTS While ERL trough concentration (C0 ) was within the targeted range for most patients, over 30-fold variability in the dose-adjusted ERL C0 was observed. Regression analysis revealed a significant effect of the non-functional CYP3A5*3 variant on the dose-adjusted ERL C0 (P = 0.031). ERL dose requirement was 0.02 mg/kg/day higher in patients with CYP3A5*1/*3 genotype compared to patients with CYP3A5*3/*3 to reach the targeted C0 (P = 0.041). ERL therapy substantially improved estimated glomerular filtration rate (28.6 ± 6.6 ml/min/1.73m(2) ) in patients with baseline kidney dysfunction. CONCLUSION ERL pharmacokinetics in HTx recipients is highly variable. Our preliminary data on patients on a CNI-free therapy regimen suggest that CYP3A5 genetic variation may contribute to this variability. This article is protected by copyright. All rights reserved.

Expansion load: recessive mutations and the role of standing genetic variation

Expanding populations incur a mutation burden – the so-called expansion load. Previous studies of expansion load have focused on codominant mutations. An important consequence of this assumption is that expansion load stems exclusively from the accumulation of new mutations occurring in individuals living at the wave front. Using individual-based simulations, we study here the dynamics of standing genetic variation at the front of expansions, and its consequences on mean fitness if mutations are recessive. We find that deleterious genetic diversity is quickly lost at the front of the expansion, but the loss of deleterious mutations at some loci is compensated by an increase of their frequencies at other loci. The frequency of deleterious homozygotes therefore increases along the expansion axis, whereas the average number of deleterious mutations per individual remains nearly constant across the species range. This reveals two important differences to codominant models: (i) mean fitness at the front of the expansion drops much faster if mutations are recessive, and (ii) mutation load can increase during the expansion even if the total number of deleterious mutations per individual remains constant. We use our model to make predictions about the shape of the site frequency spectrum at the front of range expansion, and about correlations between heterozygosity and fitness in different parts of the species range. Importantly, these predictions provide opportunities to empirically validate our theoretical results. We discuss our findings in the light of recent results on the distribution of deleterious genetic variation across human populations and link them to empirical results on the correlation of heterozygosity and fitness found in many natural range expansions.

Identification of genetic variation influencing apolipoprotein E levels: Follow-up of genome-wide linkage scans in the GENOA study

Apolipoprotein E (ApoE) plays a major role in the metabolism of high density and low density lipoproteins (HDL and LDL). Its common protein isoforms (E2, E3, E4) are risk factors for coronary artery disease (CAD) and explain between 16 to 23% of the inter-individual variation in plasma apoE levels. Linkage analysis has been completed for plasma apoE levels in the GENOA study (Genetic Epidemiology Network of Atherosclerosis). After stratification of the population by lipoprotein levels and body mass index (BMI) to create more homogeneity with regard to biological context for apoE levels, Hispanic families showed significant linkage on chromosome 17q for two strata (LOD=2.93 at 104 cM for a low cholesterol group, LOD=3.04 at 111 cM for a low cholesterol, high HDLC group). Replication of 17q linkage was observed for apoB and apoE levels in the unstratified Hispanic and African-American populations, and for apoE levels in African-American families. Replication of this 17q linkage in different populations and strata provides strong support for the presence of gene(s) in this region with significant roles in the determination of inter-individual variation in plasma apoE levels. Through a positional and functional candidate gene approach, ten genes were identified in the 17q linked region, and 62 polymorphisms in these genes were genotyped in the GENOA families. Association analysis was performed with FBAT, GEE, and variance-component based tests followed by conditional linkage analysis. Association studies with partial coverage of TagSNPs in the gene coding for apolipoprotein H (APOH) were performed, and significant results were found for 2 SNPs (APOH_20951 and APOH_05407) in the Hispanic low cholesterol strata accounting for 3.49% of the inter-individual variation in plasma apoE levels. Among the other candidate genes, we identified a haplotype block in the ACE1 gene that contains two major haplotypes associated with apoE levels as well as total cholesterol, apoB and LDLC levels in the unstratified Hispanic population. Identifying genes responsible for the remaining 60% of inter-individual variation in plasma apoE level, will yield new insights into the understanding of genetic interactions involved in the lipid metabolism, and a more precise understanding of the risk factors leading to CAD. ^

Permanence or change? The meaning of genetic variation

Selected aspects of the evolutionary process and more specifically of the genetic variation are considered, with an emphasis in studies performed by my group. One key aspect of evolution seems to be the concomitant occurrence of dichotomic, contradictory (dialect) processes. Genetic variation is structured, and the dynamics of change at one level is not necessarily paralleled by that in another. The pathogenesis-related protein superfamily can be cited as an example in which permanence (the maintenance of certain key genetic features) coexists with change (modifications that led to different functions in different classes of organisms). Relationships between structure and function are exemplified by studies with hemoglobin Porto Alegre. The genetic structure of tribal populations may differ in important aspects from that of industrialized societies. Evolutionary histories also may differ when considered through the investigation of patrilineal or matrilineal lineages. Global evaluations taking into consideration all of these aspects are needed if we really want to understand the meaning of genetic variation.

Functional and physiological consequences of genetic variation at phosphoglucose isomerase: Heat shock protein expression is related to enzyme genotype in a montane beetle

Allele frequency variation at the phosphoglucose isomerase (PGI) locus in Californian populations of the beetle Chrysomela aeneicollis suggests that PGI may be undergoing natural selection. We quantified (i) apparent Michaelis-Menten constant (Km) of fructose 6-phosphate at different temperatures and (ii) thermal stability for three common PGI genotypes (1–1, 1–4, and 4–4). We also measured air temperature (Ta) and beetle body temperature (Tb) in three montane drainages in the Sierra Nevada, California. Finally, we measured 70-kDa heat shock protein (Hsp70) expression in field-collected and laboratory-acclimated beetles. We found that PGI allele 1 predominated in the northernmost drainage, Rock Creek (RC), which was also significantly cooler than the southernmost drainage, Big Pine Creek (BPC), where PGI allele 4 predominated. Allele frequencies and air temperatures were intermediate in the middle drainage, Bishop Creek (BC). Differences among genotypes in Km (1–1 > 1–4 > 4–4) and thermal stability (4–4 > 1–4 > 1–1) followed a pattern consistent with temperature adaptation. In nature, Tb was closely related to Ta. Hsp70 expression in adult beetles decreased with elevation and differed among drainages (BPC > BC > RC). After laboratory acclimation (8 days, 20°C day, 4°C night) and heat shock (4 h, 28–36°C), Hsp70 expression was greater for RC than BPC beetles. In RC, field-collected beetles homozygous for PGI 1–1 had higher Hsp70 levels than heterozygotes or a 4–4 homozygote. These results reveal functional and physiological differences among PGI genotypes, which suggest that montane populations of this beetle are locally adapted to temperature.