How accurate is the current picture of human genetic variation?

Our understanding of the distribution of worldwide human genomic diversity has greatly increased over recent years thanks to the availability of large data sets derived from short tandem repeats (STRs), insertion deletion polymorphisms (indels) and single nucleotide polymorphisms (SNPs). A concern, however, is that the current picture of worldwide human genomic diversity may be inaccurate because of biases in the selection process of genetic markers (so-called 'ascertainment bias'). To evaluate this problem, we first compared the distribution of genomic diversity between these three types of genetic markers in the populations from the HGDP-CEPH panel for evidence of bias or incongruities. In a second step, using a very relaxed set of criteria to prevent the intrusion of bias, we developed a new set of unbiased STR markers and compared the results against those from available panels. Contrarily to recent claims, our results show that the STR markers suffer from no discernible bias, and can thus be used as a baseline reference for human genetic diversity and population differentiation. The bias on SNPs is moderate compared to that on the set of indels analysed, which we recommend should be avoided for work describing the distribution of human genetic diversity or making inference on human settlement history.

A genome-wide approach accounting for body mass index identifies genetic variants influencing fasting glycemic traits and insulin resistance.

Recent genome-wide association studies have described many loci implicated in type 2 diabetes (T2D) pathophysiology and β-cell dysfunction but have contributed little to the understanding of the genetic basis of insulin resistance. We hypothesized that genes implicated in insulin resistance pathways might be uncovered by accounting for differences in body mass index (BMI) and potential interactions between BMI and genetic variants. We applied a joint meta-analysis approach to test associations with fasting insulin and glucose on a genome-wide scale. We present six previously unknown loci associated with fasting insulin at P < 5 × 10(-8) in combined discovery and follow-up analyses of 52 studies comprising up to 96,496 non-diabetic individuals. Risk variants were associated with higher triglyceride and lower high-density lipoprotein (HDL) cholesterol levels, suggesting a role for these loci in insulin resistance pathways. The discovery of these loci will aid further characterization of the role of insulin resistance in T2D pathophysiology.

Infantile spasms is associated with deletion of the MAGI2 gene on chromosome 7q11.23-q21.11.

Infantile spasms (IS) is the most severe and common form of epilepsy occurring in the first year of life. At least half of IS cases are idiopathic in origin, with others presumed to arise because of brain insult or malformation. Here, we identify a locus for IS by high-resolution mapping of 7q11.23-q21.1 interstitial deletions in patients. The breakpoints delineate a 500 kb interval within the MAGI2 gene (1.4 Mb in size) that is hemizygously disrupted in 15 of 16 participants with IS or childhood epilepsy, but remains intact in 11 of 12 participants with no seizure history. MAGI2 encodes the synaptic scaffolding protein membrane-associated guanylate kinase inverted-2 that interacts with Stargazin, a protein also associated with epilepsy in the stargazer mouse.

A missense mutation in myelin oligodendrocyte glycoprotein as a cause of familial narcolepsy with cataplexy.

Narcolepsy is a rare sleep disorder characterized by excessive daytime sleepiness and cataplexy. Familial narcolepsy accounts for less than 10% of all narcolepsy cases. However, documented multiplex families are very rare and causative mutations have not been identified to date. To identify a causative mutation in familial narcolepsy, we performed linkage analysis in the largest ever reported family, which has 12 affected members, and sequenced coding regions of the genome (exome sequencing) of three affected members with narcolepsy and cataplexy. We successfully mapped a candidate locus on chromosomal region 6p22.1 (LOD score ¼ 3.85) by linkage analysis. Exome sequencing identified a missense mutation in the second exon of MOG within the linkage region. A c.398C>G mutation was present in all affected family members but absent in unaffected members and 775 unrelated control subjects. Transient expression of mutant myelin oligodendrocyte glycoprotein (MOG) in mouse oligodendrocytes showed abnormal subcellular localization, suggesting an altered function of the mutant MOG. MOG has recently been linked to various neuropsychiatric disorders and is considered as a key autoantigen in multiple sclerosis and in its animal model, experimental autoimmune encephalitis. Our finding of a pathogenic MOG mutation highlights a major role for myelin and oligodendrocytes in narcolepsy and further emphasizes glial involvement in neurodegeneration and neurobehavioral disorders. [corrected].

Importance of IL28B gene polymorphisms in hepatitis C virus genotype 2 and 3 infected patients.

BACKGROUND & AIMS: Genetic variation in the interleukin 28B (IL28B) gene has been associated with the response to interferon-alfa/ribavirin therapy in hepatitis C virus (HCV) genotype 1-infected patients. The importance of three IL28B single nucleotide polymorphisms (rs8099917, rs12980275 and rs12979860) for HCV genotype 2/3-infected patients is unknown. METHODS: In patients with chronic hepatitis C genotype 2/3 (n=267), IL28B host genotypes (rs8099917, rs12980275 and rs12979860) were analyzed for associations with sustained virologic response (SVR) to antiviral therapy with (pegylated) interferon-alfa and ribavirin and with respect to epidemiological, biochemical, and virological parameters. For comparison, hepatitis C genotype 1 patients (n=378) and healthy controls (n=200) were included. RESULTS: The rs12979860 CC genotype, lower age, and genotype 2 were significantly associated with SVR in HCV genotype 2/3-infected patients (p=0.01, p=0.03 and p=0.03, respectively). No association was observed for rs8099917 and rs12980275. In addition, an SVR in patients with rapid virologic response (RVR) was associated with the rs12979860 CC genotype (p=0.05), while for non-RVR no association was found. Furthermore, a significant association with a higher baseline viral load was observed for all three IL28B genotypes in genotype 1/2/3-infected patients. Finally, increasing frequencies of the rs12979860 CC genotypes were observed in genotype 1- (33.9%), genotype 3- (38.9%), and genotype 2-infected (51.9%) patients in comparison with healthy controls (49.0%) (p<0.01). CONCLUSIONS: In genotype 2/3-infected patients, rs12979860 was significantly associated with SVR. The frequency of the rs12979860 CC genotype is lower in HCV genotype 1 vs. genotype 2/3 patients. All major IL28B genotypes are associated with HCV-RNA concentration.

Impact of recombination on genetic variability within Staphylococcus aureus clonal complexes.

The population structure of Staphylococcus aureus is generally described as highly clonal and is consequently subdivided into several clonal complexes (CCs). Recent data suggested that recombination might occur more frequently within than among CCs. To test this hypothesis as well as to understand how genetic diversity is created in S. aureus, we analyzed a collection of 182 isolates with MLST and five highly variable core adhesion (ADH) genes. As expected the polymorphism of ADH genes was higher than MLST genes. However both categories of genes showed low within CCs diversity with a dominant haplotype and its single nucleotide variants. Several recombination events were detected but none involved intra-CC recombination. This did not confirm the hypothesis of higher recombination within CCs. Nevertheless, molecular analyses of variance indicated that these few recombination events have a significant impact on the genetic diversity within CCs. In addition, although most ADH genes were under purifying selection, signs of positive selection associated with a recombinant group were detected. These data highlight the importance of recombination on the evolution of the highly clonal S. aureus and suggest that recombination when combined with demographic mechanisms as well as selection might favor the rapid creation of new clonal complexes.

Genetic analyses reveal a role for vitamin D insufficiency in HCV-associated hepatocellular carcinoma development

BACKGROUND: Vitamin D insufficiency has been associated with the occurrence of various types of cancer, but causal relationships remain elusive. We therefore aimed to determine the relationship between genetic determinants of vitamin D serum levels and the risk of developing hepatitis C virus (HCV)-related hepatocellular carcinoma (HCC). METHODOLOGY/PRINCIPAL FINDINGS: Associations between CYP2R1, GC, and DHCR7 genotypes that are determinants of reduced 25-hydroxyvitamin D (25[OH]D3) serum levels and the risk of HCV-related HCC development were investigated for 1279 chronic hepatitis C patients with HCC and 4325 without HCC, respectively. The well-known associations between CYP2R1 (rs1993116, rs10741657), GC (rs2282679), and DHCR7 (rs7944926, rs12785878) genotypes and 25(OH)D3 serum levels were also apparent in patients with chronic hepatitis C. The same genotypes of these single nucleotide polymorphisms (SNPs) that are associated with reduced 25(OH)D3 serum levels were found to be associated with HCV-related HCC (P = 0.07 [OR = 1.13, 95% CI = 0.99-1.28] for CYP2R1, P = 0.007 [OR = 1.56, 95% CI = 1.12-2.15] for GC, P = 0.003 [OR = 1.42, 95% CI = 1.13-1.78] for DHCR7; ORs for risk genotypes). In contrast, no association between these genetic variations and liver fibrosis progression rate (P>0.2 for each SNP) or outcome of standard therapy with pegylated interferon-α and ribavirin (P>0.2 for each SNP) was observed, suggesting a specific influence of the genetic determinants of 25(OH)D3 serum levels on hepatocarcinogenesis. CONCLUSIONS/SIGNIFICANCE: Our data suggest a relatively weak but functionally relevant role for vitamin D in the prevention of HCV-related hepatocarcinogenesis.

RANTES/CCL5 and risk for coronary events: results from the MONICA/KORA Augsburg case-cohort, Athero-Express and CARDIoGRAM studies.

BACKGROUND: The chemokine RANTES (regulated on activation, normal T-cell expressed and secreted)/CCL5 is involved in the pathogenesis of cardiovascular disease in mice, whereas less is known in humans. We hypothesised that its relevance for atherosclerosis should be reflected by associations between CCL5 gene variants, RANTES serum concentrations and protein levels in atherosclerotic plaques and risk for coronary events. METHODS AND FINDINGS: We conducted a case-cohort study within the population-based MONICA/KORA Augsburg studies. Baseline RANTES serum levels were measured in 363 individuals with incident coronary events and 1,908 non-cases (mean follow-up: 10.2±4.8 years). Cox proportional hazard models adjusting for age, sex, body mass index, metabolic factors and lifestyle factors revealed no significant association between RANTES and incident coronary events (HR [95% CI] for increasing RANTES tertiles 1.0, 1.03 [0.75-1.42] and 1.11 [0.81-1.54]). None of six CCL5 single nucleotide polymorphisms and no common haplotype showed significant associations with coronary events. Also in the CARDIoGRAM study (>22,000 cases, >60,000 controls), none of these CCL5 SNPs was significantly associated with coronary artery disease. In the prospective Athero-Express biobank study, RANTES plaque levels were measured in 606 atherosclerotic lesions from patients who underwent carotid endarterectomy. RANTES content in atherosclerotic plaques was positively associated with macrophage infiltration and inversely associated with plaque calcification. However, there was no significant association between RANTES content in plaques and risk for coronary events (mean follow-up 2.8±0.8 years). CONCLUSIONS: High RANTES plaque levels were associated with an unstable plaque phenotype. However, the absence of associations between (i) RANTES serum levels, (ii) CCL5 genotypes and (iii) RANTES content in carotid plaques and either coronary artery disease or incident coronary events in our cohorts suggests that RANTES may not be a novel coronary risk biomarker. However, the potential relevance of RANTES levels in platelet-poor plasma needs to be investigated in further studies.

Common genetic variation and the control of HIV-1 in humans.

To extend the understanding of host genetic determinants of HIV-1 control, we performed a genome-wide association study in a cohort of 2,554 infected Caucasian subjects. The study was powered to detect common genetic variants explaining down to 1.3% of the variability in viral load at set point. We provide overwhelming confirmation of three associations previously reported in a genome-wide study and show further independent effects of both common and rare variants in the Major Histocompatibility Complex region (MHC). We also examined the polymorphisms reported in previous candidate gene studies and fail to support a role for any variant outside of the MHC or the chemokine receptor cluster on chromosome 3. In addition, we evaluated functional variants, copy-number polymorphisms, epistatic interactions, and biological pathways. This study thus represents a comprehensive assessment of common human genetic variation in HIV-1 control in Caucasians.

Genetic differentiation of continental and island populations of Bombus terrestris (Hymenoptera: Apidae) in Europe.

Ten microsatellite loci and a partial sequence of the COII mitochondrial gene were used to investigate genetic differentiation in B. terrestris, a bumble bee of interest for its high-value crop pollination. The analysis included eight populations from the European continent, five from Mediterranean islands (six subspecies altogether) and one from Tenerife (initially described as a colour form of B. terrestris but recently considered as a separate species, B. canariensis). Eight of the 10 microsatellite loci displayed high levels of polymorphism in most populations. In B. terrestris populations, the total number of alleles detected per polymorphic locus ranged from 3 to 16, with observed allelic diversity from 3.8 +/- 0.5 to 6.5 +/- 1.4 and average calculated heterozygosities from 0.41 +/- 0.09 to 0.65 +/- 0.07. B. canariensis showed a significantly lower average calculated heterozygosity (0.12 +/- 0.08) and observed allelic diversity (1.5 +/- 0.04) as compared to both continental and island populations of B. terrestris. No significant differentiation was found among populations of B. terrestris from the European continent. In contrast, island populations were all significantly and most of them strongly differentiated from continental populations. B. terrestris mitochondrial DNA is characterized by a low nucleotide diversity: 0.18% +/- 0.07%, 0.20% +/- 0.04% and 0.27% +/- 0.04% for the continental populations, the island populations and all populations together, respectively. The only haplotype found in the Tenerife population differs by a single nucleotide substitution from the most common continental haplotype of B. terrestris. This situation, identical to that of Tyrrhenian islands populations and quite different from that of B. lucorum (15 substitutions between terrestris and lucorum mtDNA) casts doubts on the species status of B. canariensis. The large genetic distance between the Tenerife and B. terrestris populations estimated from microsatellite data result, most probably, from a severe bottleneck in the Canary island population. Microsatellite and mitochondrial DNA data call for the protection of the island populations of B. terrestris against importation of bumble bees of foreign origin which are used as crop pollinators.

Predicting stroke through genetic risk functions: the CHARGE Risk Score Project.

BACKGROUND AND PURPOSE: Beyond the Framingham Stroke Risk Score, prediction of future stroke may improve with a genetic risk score (GRS) based on single-nucleotide polymorphisms associated with stroke and its risk factors. METHODS: The study includes 4 population-based cohorts with 2047 first incident strokes from 22,720 initially stroke-free European origin participants aged ≥55 years, who were followed for up to 20 years. GRSs were constructed with 324 single-nucleotide polymorphisms implicated in stroke and 9 risk factors. The association of the GRS to first incident stroke was tested using Cox regression; the GRS predictive properties were assessed with area under the curve statistics comparing the GRS with age and sex, Framingham Stroke Risk Score models, and reclassification statistics. These analyses were performed per cohort and in a meta-analysis of pooled data. Replication was sought in a case-control study of ischemic stroke. RESULTS: In the meta-analysis, adding the GRS to the Framingham Stroke Risk Score, age and sex model resulted in a significant improvement in discrimination (all stroke: Δjoint area under the curve=0.016, P=2.3×10(-6); ischemic stroke: Δjoint area under the curve=0.021, P=3.7×10(-7)), although the overall area under the curve remained low. In all the studies, there was a highly significantly improved net reclassification index (P<10(-4)). CONCLUSIONS: The single-nucleotide polymorphisms associated with stroke and its risk factors result only in a small improvement in prediction of future stroke compared with the classical epidemiological risk factors for stroke.

Association of the hemochromatosis gene with pazopanib-induced transaminase elevation in renal cell carcinoma.

BACKGROUND & AIMS: Pazopanib has demonstrated clinical benefit in patients with advanced renal cell carcinoma (RCC) and is generally well tolerated. However, transaminase elevations have commonly been observed. This 2-stage study sought to identify genetic determinants of alanine transaminase (ALT) elevations in pazopanib-treated white patients with RCC.¦METHODS: Data from two separate clinical studies were used to examine the association of genetic polymorphisms with maximum on-treatment ALT levels.¦RESULTS: Of 6852 polymorphisms in 282 candidate genes examined in an exploratory dataset of 115 patients, 92 polymorphisms in 40 genes were significantly associated with ALT elevation (p<0.01). Two markers (rs2858996 and rs707889) in the HFE gene, which are not yet known to be associated with hemochromatosis, showed evidence for replication. Because of multiple comparisons, there was a 12% likelihood the replication occurred by chance. These two markers demonstrated strong linkage disequilibrium (r(2)=0.99). In the combined dataset, median (25-75th percentile) maximum ALT values were 1.2 (0.7-1.9), 1.1 (0.8-2.5), and 5.4 (1.9-7.6)×ULN for rs2858996 GG (n=148), GT (n=82), and TT (n=1 2) genotypes, respectively. All 12 TT patients had a maximum ALT>ULN, and 8 (67%) had ALT≥3×ULN. The odds ratio (95% CI) for ALT≥3×ULN for TT genotype was 39.7 (2.2-703.7) compared with other genotypes. As a predictor of ALT≥3×ULN, the TT genotype had a negative predictive value of 0.83 and positive predictive value of 0.67. No TT patients developed liver failure.¦CONCLUSIONS: The rs2858996/rs707889 polymorphisms in the HFE gene may be associated with reversible ALT elevation in pazo-panib-treated patients with RCC.

Genome-wide association analysis of blood-pressure traits in African-ancestry individuals reveals common associated genes in African and non-African populations.

High blood pressure (BP) is more prevalent and contributes to more severe manifestations of cardiovascular disease (CVD) in African Americans than in any other United States ethnic group. Several small African-ancestry (AA) BP genome-wide association studies (GWASs) have been published, but their findings have failed to replicate to date. We report on a large AA BP GWAS meta-analysis that includes 29,378 individuals from 19 discovery cohorts and subsequent replication in additional samples of AA (n = 10,386), European ancestry (EA) (n = 69,395), and East Asian ancestry (n = 19,601). Five loci (EVX1-HOXA, ULK4, RSPO3, PLEKHG1, and SOX6) reached genome-wide significance (p < 1.0 × 10(-8)) for either systolic or diastolic BP in a transethnic meta-analysis after correction for multiple testing. Three of these BP loci (EVX1-HOXA, RSPO3, and PLEKHG1) lack previous associations with BP. We also identified one independent signal in a known BP locus (SOX6) and provide evidence for fine mapping in four additional validated BP loci. We also demonstrate that validated EA BP GWAS loci, considered jointly, show significant effects in AA samples. Consequently, these findings suggest that BP loci might have universal effects across studied populations, demonstrating that multiethnic samples are an essential component in identifying, fine mapping, and understanding their trait variability.

Molecular Proxies for Climate Maladaptation in a Long-Lived Tree (Pinus pinaster Aiton, Pinaceae).

Understanding adaptive genetic responses to climate change is a main challenge for preserving biological diversity. Successful predictive models for climate-driven range shifts of species depend on the integration of information on adaptation, including that derived from genomic studies. Long-lived forest trees can experience substantial environmental change across generations, which results in a much more prominent adaptation lag than in annual species. Here, we show that candidate-gene SNPs (single nucleotide polymorphisms) can be used as predictors of maladaptation to climate in maritime pine (Pinus pinaster Aiton), an outcrossing long-lived keystone tree. A set of 18 SNPs potentially associated with climate, 5 of them involving amino acid-changing variants, were retained after performing logistic regression, latent factor mixed models, and Bayesian analyses of SNP-climate correlations. These relationships identified temperature as an important adaptive driver in maritime pine and highlighted that selective forces are operating differentially in geographically discrete gene pools. The frequency of the locally advantageous alleles at these selected loci was strongly correlated with survival in a common garden under extreme (hot and dry) climate conditions, which suggests that candidate-gene SNPs can be used to forecast the likely destiny of natural forest ecosystems under climate change scenarios. Differential levels of forest decline are anticipated for distinct maritime pine gene pools. Geographically defined molecular proxies for climate adaptation will thus critically enhance the predictive power of range-shift models and help establish mitigation measures for long-lived keystone forest trees in the face of impending climate change.