Does phylogeographic structure relate to climatic niche divergence? A test using maritime pine (Pinus pinaster Ait.)

Aim To disentangle the effects of environmental and geographical processes driving phylogenetic distances among clades of maritime pine (Pinus pinaster). To assess the implications for conservation management of combining molecular information with species distribution models (SDMs; which predict species distribution based on known occurrence records and on environmental variables). Location Western Mediterranean Basin and European Atlantic coast. Methods We undertook two cluster analyses for eight genetically defined pine clades based on climatic niche and genetic similarities. We assessed niche similarity by means of a principal component analysis and Schoener's D metric. To calculate genetic similarity, we used the unweighted pair group method with arithmetic mean based on Nei's distance using 266 single nucleotide polymorphisms. We then assessed the contribution of environmental and geographical distances to phylogenetic distance by means of Mantel regression with variance partitioning. Finally, we compared the projection obtained from SDMs fitted from the species level (SDMsp) and composed from the eight clade-level models (SDMcm). Results Genetically and environmentally defined clusters were identical. Environmental and geographical distances explained 12.6% of the phylogenetic distance variation and, overall, geographical and environmental overlap among clades was low. Large differences were detected between SDMsp and SDMcm (57.75% of disagreement in the areas predicted as suitable). Main conclusions The genetic structure within the maritime pine subspecies complex is primarily a consequence of its demographic history, as seen by the high proportion of unexplained variation in phylogenetic distances. Nevertheless, our results highlight the contribution of local environmental adaptation in shaping the lower-order, phylogeographical distribution patterns and spatial genetic structure of maritime pine: (1) genetically and environmentally defined clusters are consistent, and (2) environment, rather than geography, explained a higher proportion of variation in phylogenetic distance. SDMs, key tools in conservation management, better characterize the fundamental niche of the species when they include molecular information.

Association between parental history and genetic risk scores for coronary heart disease prediction: The population-based CoLaus study.

BACKGROUND AND AIMS: Parental history (PH) and genetic risk scores (GRSs) are separately associated with coronary heart disease (CHD), but evidence regarding their combined effects is lacking. We aimed to evaluate the joint associations and predictive ability of PH and GRSs for incident CHD. METHODS: Data for 4283 Caucasians were obtained from the population-based CoLaus Study, over median follow-up time of 5.6 years. CHD was defined as incident myocardial infarction, angina, percutaneous coronary revascularization or bypass grafting. Single nucleotide polymorphisms for CHD identified by genome-wide association studies were used to construct unweighted and weighted versions of three GRSs, comprising of 38, 53 and 153 SNPs respectively. RESULTS: PH was associated with higher values of all weighted GRSs. After adjustment for age, sex, smoking, diabetes, systolic blood pressure, low and high density lipoprotein cholesterol, PH was significantly associated with CHD [HR 2.61, 95% CI (1.47-4.66)] and further adjustment for GRSs did not change this estimate. Similarly, one standard deviation change of the weighted 153-SNPs GRS was significantly associated with CHD [HR 1.50, 95% CI (1.26-1.80)] and remained so, after further adjustment for PH. The weighted, 153-SNPs GRS, but not PH, modestly improved discrimination [(C-index improvement, 0.016), p = 0.048] and reclassification [(NRI improvement, 8.6%), p = 0.027] beyond cardiovascular risk factors. After including both the GRS and PH, model performance improved further [(C-index improvement, 0.022), p = 0.006]. CONCLUSION: After adjustment for cardiovascular risk factors, PH and a weighted, polygenic GRS were jointly associated with CHD and provided additive information for coronary events prediction.

Local effects drive heterozygosity-fitness correlations in an outcrossing long-lived tree.

Heterozygosity-fitness correlations (HFCs) have been used to understand the complex interactions between inbreeding, genetic diversity and evolution. Although frequently reported for decades, evidence for HFCs was often based on underpowered studies or inappropriate methods, and hence their underlying mechanisms are still under debate. Here, we used 6100 genome-wide single nucleotide polymorphisms (SNPs) to test for general and local effect HFCs in maritime pine (Pinus pinaster Ait.), an iconic Mediterranean forest tree. Survival was used as a fitness proxy, and HFCs were assessed at a four-site common garden under contrasting environmental conditions (total of 16 288 trees). We found no significant correlations between genome-wide heterozygosity and fitness at any location, despite variation in inbreeding explaining a substantial proportion of the total variance for survival. However, four SNPs (including two non-synonymous mutations) were involved in significant associations with survival, in particular in the common gardens with higher environmental stress, as shown by a novel heterozygosity-fitness association test at the species-wide level. Fitness effects of SNPs involved in significant HFCs were stable across maritime pine gene pools naturally growing in distinct environments. These results led us to dismiss the general effect hypothesis and suggested a significant role of heterozygosity in specific candidate genes for increasing fitness in maritime pine. Our study highlights the importance of considering the species evolutionary and demographic history and different spatial scales and testing environments when assessing and interpreting HFCs.

Impact of a cis-associated gene expression SNP on chromosome 20q11.22 on bipolar disorder susceptibility, hippocampal structure and cognitive performance.

BackgroundBipolar disorder is a highly heritable polygenic disorder. Recent enrichment analyses suggest that there may be true risk variants for bipolar disorder in the expression quantitative trait loci (eQTL) in the brain.AimsWe sought to assess the impact of eQTL variants on bipolar disorder risk by combining data from both bipolar disorder genome-wide association studies (GWAS) and brain eQTL.MethodTo detect single nucleotide polymorphisms (SNPs) that influence expression levels of genes associated with bipolar disorder, we jointly analysed data from a bipolar disorder GWAS (7481 cases and 9250 controls) and a genome-wide brain (cortical) eQTL (193 healthy controls) using a Bayesian statistical method, with independent follow-up replications. The identified risk SNP was then further tested for association with hippocampal volume (n = 5775) and cognitive performance (n = 342) among healthy individuals.ResultsIntegrative analysis revealed a significant association between a brain eQTL rs6088662 on chromosome 20q11.22 and bipolar disorder (log Bayes factor = 5.48; bipolar disorder P = 5.85×10(-5)). Follow-up studies across multiple independent samples confirmed the association of the risk SNP (rs6088662) with gene expression and bipolar disorder susceptibility (P = 3.54×10(-8)). Further exploratory analysis revealed that rs6088662 is also associated with hippocampal volume and cognitive performance in healthy individuals.ConclusionsOur findings suggest that 20q11.22 is likely a risk region for bipolar disorder; they also highlight the informative value of integrating functional annotation of genetic variants for gene expression in advancing our understanding of the biological basis underlying complex disorders, such as bipolar disorder.

Genome-wide Association Studies Identify Genetic Loci Associated With Albuminuria in Diabetes.

Elevated concentrations of albumin in the urine, albuminuria, are a hallmark of diabetic kidney disease and are associated with an increased risk for end-stage renal disease and cardiovascular events. To gain insight into the pathophysiological mechanisms underlying albuminuria, we conducted meta-analyses of genome-wide association studies and independent replication in up to 5,825 individuals of European ancestry with diabetes and up to 46,061 without diabetes, followed by functional studies. Known associations of variants in CUBN, encoding cubilin, with the urinary albumin-to-creatinine ratio (UACR) were confirmed in the overall sample (P = 2.4 × 10(-10)). Gene-by-diabetes interactions were detected and confirmed for variants in HS6ST1 and near RAB38/CTSC. Single nucleotide polymorphisms at these loci demonstrated a genetic effect on UACR in individuals with but not without diabetes. The change in the average UACR per minor allele was 21% for HS6ST1 (P = 6.3 × 10(-7)) and 13% for RAB38/CTSC (P = 5.8 × 10(-7)). Experiments using streptozotocin-induced diabetic Rab38 knockout and control rats showed higher urinary albumin concentrations and reduced amounts of megalin and cubilin at the proximal tubule cell surface in Rab38 knockout versus control rats. Relative expression of RAB38 was higher in tubuli of patients with diabetic kidney disease compared with control subjects. The loci identified here confirm known pathways and highlight novel pathways influencing albuminuria.

Prevalence and determinants of periodic limb movements in the general population.

OBJECTIVE: Periodic limb movements during sleep (PLMS) are sleep phenomena characterized by periodic episodes of repetitive stereotyped limb movements. The aim of this study was to describe the prevalence and determinants of PLMS in a middle to older aged general population. METHODS: Data from 2,162 subjects (51.2% women, mean age = 58.4 ± 11.1 years) participating in a population-based study (HypnoLaus, Lausanne, Switzerland) were collected. Assessments included laboratory tests, sociodemographic data, personal and treatment history, and full polysomnography at home. PLMS index (PLMSI) was determined, and PLMSI > 15/h was considered as significant. RESULTS: Prevalence of PLMSI > 15/h was 28.6% (31.3% in men, 26% in women). Compared to subjects with PLMSI ≤ 15/h, subjects with PLMSI > 15/h were older (p < 0.001), were predominantly males (p = 0.007), had a higher proportion of restless legs syndrome (RLS; p < 0.001), had a higher body mass index (p = 0.001), and had a lower mean glomerular filtration rate (p < 0.001). Subjects with PLMSI > 15/h also had a higher prevalence of diabetes, hypertension, and beta-blocker or hypnotic treatments. The prevalence of antidepressant use was higher, but not statistically significant (p = 0.07). Single nucleotide polymorphisms (SNPs) within BTBD9 (rs3923809), TOX3 (rs3104788), and MEIS1 (rs2300478) genes were significantly associated with PLSMI > 15/h. Conversely, mean hemoglobin and ferritin levels were similar in both groups. In the multivariate analysis, age, male gender, antidepressant intake, RLS, and rs3923809, rs3104788, and rs2300478 SNPs were independently associated with PLMSI > 15/h. INTERPRETATION: PLMS are highly prevalent in our middle-aged European population. Age, male gender, RLS, antidepressant treatment, and specific BTBD9, TOX3, and MEIS1 SNP distribution are independent predictors of PLMSI > 15/h. ANN NEUROL 2016;79:464-474.

Hybridization Capture Using RAD Probes (hyRAD), a New Tool for Performing Genomic Analyses on Collection Specimens.

In the recent years, many protocols aimed at reproducibly sequencing reduced-genome subsets in non-model organisms have been published. Among them, RAD-sequencing is one of the most widely used. It relies on digesting DNA with specific restriction enzymes and performing size selection on the resulting fragments. Despite its acknowledged utility, this method is of limited use with degraded DNA samples, such as those isolated from museum specimens, as these samples are less likely to harbor fragments long enough to comprise two restriction sites making possible ligation of the adapter sequences (in the case of double-digest RAD) or performing size selection of the resulting fragments (in the case of single-digest RAD). Here, we address these limitations by presenting a novel method called hybridization RAD (hyRAD). In this approach, biotinylated RAD fragments, covering a random fraction of the genome, are used as baits for capturing homologous fragments from genomic shotgun sequencing libraries. This simple and cost-effective approach allows sequencing of orthologous loci even from highly degraded DNA samples, opening new avenues of research in the field of museum genomics. Not relying on the restriction site presence, it improves among-sample loci coverage. In a trial study, hyRAD allowed us to obtain a large set of orthologous loci from fresh and museum samples from a non-model butterfly species, with a high proportion of single nucleotide polymorphisms present in all eight analyzed specimens, including 58-year-old museum samples. The utility of the method was further validated using 49 museum and fresh samples of a Palearctic grasshopper species for which the spatial genetic structure was previously assessed using mtDNA amplicons. The application of the method is eventually discussed in a wider context. As it does not rely on the restriction site presence, it is therefore not sensitive to among-sample loci polymorphisms in the restriction sites that usually causes loci dropout. This should enable the application of hyRAD to analyses at broader evolutionary scales.

Host Genetics of Response to Hepatitis B Vaccine: A Systematic Review

Background. Hepatitis B virus (HBV) is an important cause of chronic viral disease worldwide and can be life threatening. While a safe and effective vaccine is widely available, 5 to 10% of healthy vaccinees fail to achieve a protective anti-hepatitis B surface antigen antibody (anti-HBs) titer (>10mIU/ml). A limited number of studies investigated host genetics of the response to HBV vaccine. To our knowledge, no comprehensive overview of genetic polymorphisms both within and outside the HLA system has been done so far. Aim. The aim of this study was to perform a systematic review of the literature of human genetics influencing immune response after hepatitis B vaccination. Methods. Literature searches using keywords were conducted in the electronic databases Medline, Embase and ISI Web of Science the cut-off date being March 2014. After selection of papers according to stringent inclusion criteria, relevant information was systematically collected from the remaining articles, including demographic data, number of patients, schedule and type of vaccine, phenotypes, genes and single nucleotide polymorphisms (SNPs) genotyping results and their association with immune response to hepatitis B vaccine. Results. The literature search produced a total of 1968 articles from which 46 studies were kept for further analyses. From these studies, data was extracted for 19 alleles from the human leukocyte antigen (HLA) region that were reported as significant at least twice. Among those alleles, 9 were firmly associated with vaccine response outcome (DQ2 [DQB1*02 and DQB1*0201], DR3 [DRB1*03 and DRB1*0301], DR7 [DRB1*07 and DRB1*0701], C4AQ0, DPB1*0401, DQ3, DQB1*06, DRB1*01 and DRB1*13 [DRB1*1301]). In addition, data was extracted for 55 different genes from which 13 extra-HLA genes had polymorphisms that were studied by different group of investigators or by the same group with a replication study. Among the 13 genes allowing comparison, 4 genes (IL-1B, IL-2, IL-4R and IL- 6) revealed no significant data, 6 genes (IL-4, IL-10, IL-12B, IL-13, TNFA, IFNG and TLR2) were explored with inconsistent results and 2 genes (CD3Z and ITGAL) yielded promising results as their association with vaccine response was confirmed by a replication approach. Furthermore, this review produced a list of 46 SNPs from 26 genes that were associated with immune response to vaccine only once, providing novel candidates to be tested in datasets from existing genome-wide association studies (GWAS). Conclusion. To the best of our knowledge, this is the first systematic review of immunogenetic studies of response to hepatitis B vaccine. While this work reassesses the role of several HLA alleles on vaccine response outcome, the associations with polymorphisms in genes outside the HLA region were rather inconsistent. Moreover, this work produced a list of 46 significant SNPs that were reported by a single group of investigators, opening up some interesting possibilities for further research.

Effect of immune pressure on hepatitis C virus evolution: insights from a single-source outbreak.

The host's immune response to hepatitis C virus (HCV) can result in the selection of characteristic mutations (adaptations) that enable the virus to escape this response. The ability of the virus to mutate at these sites is dependent on the incoming virus, the fitness cost incurred by the mutation, and the benefit to the virus in escaping the response. Studies examining viral adaptation in chronic HCV infection have shown that these characteristic immune escape mutations can be observed at the population level as human leukocyte antigen (HLA)-specific viral polymorphisms. We examined 63 individuals with chronic HCV infection who were infected from a single HCV genotype 1b source. Our aim was to determine the extent to which the host's immune pressure affects HCV diversity and the ways in which the sequence of the incoming virus, including preexisting escape mutations, can influence subsequent mutations in recipients and infection outcomes. Conclusion: HCV sequences from these individuals revealed 29 significant associations between specific HLA types within the new hosts and variations within their viruses, which likely represent new viral adaptations. These associations did not overlap with previously reported adaptations for genotypes 1a and 3a and possibly reflected a combination of constraint due to the incoming virus and genetic distance between the strains. However, these sites accounted for only a portion of the sites in which viral diversity was observed in the new hosts. Furthermore, preexisting viral adaptations in the incoming (source) virus likely influenced the outcomes in the new hosts.

Toll-like receptor 2 (TLR2) polymorphisms are associated with reversal reaction in leprosy.

BACKGROUND: Leprosy is characterized by a spectrum of clinical manifestations that depend on the type of immune response against the pathogen. Patients may undergo immunological changes known as "reactional states" (reversal reaction and erythema nodosum leprosum) that result in major clinical deterioration. The goal of the present study was to assess the effect of Toll-like receptor 2 (TLR2) polymorphisms on susceptibility to and clinical presentation of leprosy. METHODS: Three polymorphisms in TLR2 (597C--&gt;T, 1350T--&gt;C, and a microsatellite marker) were analyzed in 431 Ethiopian patients with leprosy and 187 control subjects. The polymorphism-associated risk of developing leprosy, lepromatous (vs. tuberculoid) leprosy, and leprosy reactions was assessed by multivariate logistic regression models. RESULTS: The microsatellite and the 597C--&gt;T polymorphisms both influenced susceptibility to reversal reaction. Although the 597T allele had a protective effect (odds ratio [OR], 0.34 [95% confidence interval {CI}, 0.17-0.68]; P= .002 under the dominant model), homozygosity for the 280-bp allelic length of the microsatellite strongly increased the risk of reversal reaction (OR, 5.83 [95% CI, 1.98-17.15]; P= .001 under the recessive model). These associations were consistent among 3 different ethnic groups. CONCLUSIONS: These data suggest a significant role for TLR-2 in the occurrence of leprosy reversal reaction and provide new insights into the immunogenetics of the disease.

Relationships among endogenous ouabain, alpha-adducin polymorphisms and renal sodium handling in primary hypertension.

OBJECTIVE: The basolateral Na pump drives renotubular reabsorption. In cultured renal cells, mutant adducins, as well as sub-nanomolar ouabain concentrations, stimulate the Na-K pump. METHODS: To determine whether these factors interact and affect Na handling and blood pressure (BP) in vivo, we studied 155 untreated hypertensive patients subdivided on the basis of their plasma endogenous ouabain or alpha-adducin genotype (ADD1 Gly460Trp-rs4961). RESULTS: Under basal conditions, proximal tubular reabsorption and plasma Na were higher in patients with mutated Trp ADD1 or increased endogenous ouabain (P = 0.002 and 0.05, respectively). BPs were higher in the high plasma endogenous ouabain group (P = 0.001). Following volume loading, the increment in BP (7.73 vs. 4.81 mmHg) and the slopes of the relationship between BP and Na excretion were greater [0.017 +/- 0.002 vs. 0.009 +/- 0.003 mmHg/(muEq min)] in ADD1 Trp vs. ADD1 Gly carriers (P &lt; 0.05). BP changes were similar, whereas the slopes of the relationship between BP and Na excretion were lower [0.016 +/- 0.003 vs. 0.008 +/- 0.002 mmHg/(muEq min)] in patients with low vs. high endogenous ouabain (P &lt; 0.05). In patients with high endogenous ouabain, volume loading increased the BP in the ADD1 Trp group but not in the Gly group (P &lt; 0.05). Thus, patients with ADD1 Trp alleles are sensitive to salt and tubular Na reabsorption remains elevated after volume expansion. CONCLUSION: With saline loading, BP changes are similar in high and low endogenous ouabain patients, whereas tubular Na reabsorption increases in the high endogenous ouabain group. Saline loading unmasks differences in renal Na handling in patients with mutant adducin or high endogenous ouabain and exposes an interaction of endogenous ouabain and Trp alleles on BP.

Impact of genetic SLC28 transporter and ITPA variants on ribavirin 21 serum level, hemoglobin drop and therapeutic response in patients with HCV infection

Background: T reatment o f chronic hepatitis C i s evolving, a nd direct acting antivirals ( DAAs) are now a dded to p egylated interferon-α ( Peg- INF-α) and ribavirin (RBV) for the treatment o f hepatitis C v irus ( HCV) genotype 1 infection. DAAs c ause d ifferent side effects and can even worsen RBV induced hemolytic anemia. T herefore, identifying host genetic d eterminants of R BV bioavailability and therapeutic e fficacy will remain crucial for individualized treatment. Recent d ata showed associations between R BV induced h emolytic anemia and genetic polymorphisms o f concentrative nucleoside transporters s uch as C NT3 (SLC28A3) and i nosine t riphosphatase (ITPA). T o analyze t he association of genetic variants of SLC28 transporters and ITPA with RBV induced hemolytic anemia and treatment o utcome. Methods: I n our study, 173 patients f rom t he S wiss Hepatitis C C ohort Study and 2 2 patients from Swiss Association for the Study of the Liver study 24 (61% HCV g enotype 1, 3 9% genotypes 2 o r 3) were analyzed for SLC28A2 single nucleotide p olymorphism (SNP) rs11854484, SLC28A3 rs56350726 and SLC28A3 rs10868138 as well as ITPA SNPs rs1127354 and rs7270101. RBV serum levels during treatment were measured in 49 patients. Results: SLC28A2 r s11854484 genotype TT was associated with significantly higher dosage- and body weight-adjusted RBV levels as compared to genotypes TC and CC (p=0.04 and p=0.02 at weeks 4 and 8, respectively). ITPA SNPs rs1127354 and rs7270101 were associated with h emolytic a nemia both in genotype as w ell as i n allelic a nalyses. SLC28A3 rs56350726 genotype TT (vs. AT/AA, RR=2.1; 95% CI 1.1-4.1) as well as the T allele (vs. A; RR=1.8, 95% CI 1.1-3.2) were associated with increased SVR rates. The combined analysis of overall ITPA activity and SLC28 v ariants together revealed n o significant a dditive effects on either treatment-related anemia or SVR. Conclusions: T he newly identified association between RBV serum levels a nd SLC28A2 rs11854484 genotype as well as the replicated association of ITPA and SLC28A3 g enetic p olymorphisms w ith RBV induced hemolytic anemia and treatment r esponse underpin the need for further studies on host genetic d eterminants of R BV bioavailability and therapeutic e fficacy f or individualized treatment of chronic hepatitis C.

Genome-wide patterns of latitudinal differentiation among populations of Drosophila melanogaster from North America.

Understanding the genetic underpinnings of adaptive change is a fundamental but largely unresolved problem in evolutionary biology. Drosophila melanogaster, an ancestrally tropical insect that has spread to temperate regions and become cosmopolitan, offers a powerful opportunity for identifying the molecular polymorphisms underlying clinal adaptation. Here, we use genome-wide next-generation sequencing of DNA pools ('pool-seq') from three populations collected along the North American east coast to examine patterns of latitudinal differentiation. Comparing the genomes of these populations is particularly interesting since they exhibit clinal variation in a number of important life history traits. We find extensive latitudinal differentiation, with many of the most strongly differentiated genes involved in major functional pathways such as the insulin/TOR, ecdysone, torso, EGFR, TGFβ/BMP, JAK/STAT, immunity and circadian rhythm pathways. We observe particularly strong differentiation on chromosome 3R, especially within the cosmopolitan inversion In(3R)Payne, which contains a large number of clinally varying genes. While much of the differentiation might be driven by clinal differences in the frequency of In(3R)P, we also identify genes that are likely independent of this inversion. Our results provide genome-wide evidence consistent with pervasive spatially variable selection acting on numerous loci and pathways along the well-known North American cline, with many candidates implicated in life history regulation and exhibiting parallel differentiation along the previously investigated Australian cline.

Sequence conservation in Plasmodium falciparum alpha-helical coiled coil domains proposed for vaccine development.

BACKGROUND: The availability of the P. falciparum genome has led to novel ways to identify potential vaccine candidates. A new approach for antigen discovery based on the bioinformatic selection of heptad repeat motifs corresponding to alpha-helical coiled coil structures yielded promising results. To elucidate the question about the relationship between the coiled coil motifs and their sequence conservation, we have assessed the extent of polymorphism in putative alpha-helical coiled coil domains in culture strains, in natural populations and in the single nucleotide polymorphism data available at PlasmoDB. METHODOLOGY/PRINCIPAL FINDINGS: 14 alpha-helical coiled coil domains were selected based on preclinical experimental evaluation. They were tested by PCR amplification and sequencing of different P. falciparum culture strains and field isolates. We found that only 3 out of 14 alpha-helical coiled coils showed point mutations and/or length polymorphisms. Based on promising immunological results 5 of these peptides were selected for further analysis. Direct sequencing of field samples from Papua New Guinea and Tanzania showed that 3 out of these 5 peptides were completely conserved. An in silico analysis of polymorphism was performed for all 166 putative alpha-helical coiled coil domains originally identified in the P. falciparum genome. We found that 82% (137/166) of these peptides were conserved, and for one peptide only the detected SNPs decreased substantially the probability score for alpha-helical coiled coil formation. More SNPs were found in arrays of almost perfect tandem repeats. In summary, the coiled coil structure prediction was rarely modified by SNPs. The analysis revealed a number of peptides with strictly conserved alpha-helical coiled coil motifs. CONCLUSION/SIGNIFICANCE: We conclude that the selection of alpha-helical coiled coil structural motifs is a valuable approach to identify potential vaccine targets showing a high degree of conservation.