An overview of L-2-hydroxyglutarate dehydrogenase gene (L2HGDH) variants: a genotype-phenotype study.

L-2-Hydroxyglutaric aciduria (L2HGA) is a rare, neurometabolic disorder with an autosomal recessive mode of inheritance. Affected individuals only have neurological manifestations, including psychomotor retardation, cerebellar ataxia, and more variably macrocephaly, or epilepsy. The diagnosis of L2HGA can be made based on magnetic resonance imaging (MRI), biochemical analysis, and mutational analysis of L2HGDH. About 200 patients with elevated concentrations of 2-hydroxyglutarate (2HG) in the urine were referred for chiral determination of 2HG and L2HGDH mutational analysis. All patients with increased L2HG (n=106; 83 families) were included. Clinical information on 61 patients was obtained via questionnaires. In 82 families the mutations were detected by direct sequence analysis and/or multiplex ligation dependent probe amplification (MLPA), including one case where MLPA was essential to detect the second allele. In another case RT-PCR followed by deep intronic sequencing was needed to detect the mutation. Thirty-five novel mutations as well as 35 reported mutations and 14 nondisease-related variants are reviewed and included in a novel Leiden Open source Variation Database (LOVD) for L2HGDH variants (http://www.LOVD.nl/L2HGDH). Every user can access the database and submit variants/patients. Furthermore, we report on the phenotype, including neurological manifestations and urinary levels of L2HG, and we evaluate the phenotype-genotype relationship.

The NLRP3 inflammasome is differentially activated by pneumolysin variants and contributes to host defense in pneumococcal pneumonia.

Streptococcus pneumoniae is a leading cause of pneumonia, meningitis, and sepsis. Pneumococci can be divided into >90 serotypes that show differences in the pathogenicity and invasiveness. We tested the hypotheses that the innate immune inflammasome pathway is involved in fighting pneumococcal pneumonia and that some invasive pneumococcal types are not recognized by this pathway. We show that human and murine mononuclear cells responded to S. pneumoniae expressing hemolytic pneumolysin by producing IL-1β. This IL-1β production depended on the NOD-like receptor family, pyrin domain containing 3 (NLRP3) inflammasome. Some serotype 1, serotype 8, and serotype 7F bacteria, which have previously been associated with increased invasiveness and with production of toxins with reduced hemolytic activity, or bacterial mutants lacking pneumolysin did not stimulate notable IL-1β production. We further found that NLRP3 was beneficial for mice during pneumonia caused by pneumococci expressing hemolytic pneumolysin and was involved in cytokine production and maintenance of the pulmonary microvascular barrier. Overall, the inflammasome pathway is protective in pneumonia caused by pneumococci expressing hemolytic toxin but is not activated by clinically important pneumococcal sequence types causing invasive disease. The study indicates that a virulence factor polymorphism may substantially affect the recognition of bacteria by the innate immune system.

Crystal size distribution of periclase in contact metamorphic dolomite marbles from the southern Adamello Massif, Italy

Crystal size distributions (CSD) of periclase in contact metamorphic dolomite marbles are presented for two profiles near the Cima Uzza summit in the southern Adamello Massif (Italy). The database was combined with geochemical and petrological information to deduce the controls on the periclase-forming reaction. The contact metamorphic dolomite marbles are exposed at the contact of mafic intrusive rocks and are partially surrounded by them. Brucite is retrograde and pseudomorphs spherical periclase crystals. Prograde periclase growth is the consequence of limited infiltration of water-rich fluid at T near 605C. Stable isotope data show depletion in (13)C and (18)O over a narrow region (40 cm) near the magmatic contact, whereas the periclase-forming reaction front extends up to 4 m from the contact. CSD analyses along the two profiles show that the median grain size of the periclase crystals does not change, but that there is a progressively greater distribution of grain sizes, including a greater proportion of larger grains, with increasing distance from the contact. A qualitative model, based on the textural and geochemical data, attributes these variations in grain size to changing reaction affinities along a kinetically dispersed infiltration front. This study highlights the need to invoke disequilibrium processes for metamorphic mineral growth and expands the use of CSDs to systems of mineral formation driven by fluid infiltration.

Overlap Between Common Genetic Polymorphisms Underpinning Kidney Traits and Cardiovascular Disease Phenotypes: The CKDGen Consortium.

BACKGROUND: Chronic kidney disease is associated with cardiovascular disease. We tested for evidence of a shared genetic basis to these traits. STUDY DESIGN: We conducted 2 targeted analyses. First, we examined whether known single-nucleotide polymorphisms (SNPs) underpinning kidney traits were associated with a series of vascular phenotypes. Additionally, we tested whether vascular SNPs were associated with markers of kidney damage. Significance was set to 1.5×10(-4) (0.05/325 tests). SETTING & PARTICIPANTS: Vascular outcomes were analyzed in participants from the AortaGen (20,634), CARDIoGRAM (86,995), CHARGE Eye (15,358), CHARGE IMT (31,181), ICBP (69,395), and NeuroCHARGE (12,385) consortia. Tests for kidney outcomes were conducted in up to 67,093 participants from the CKDGen consortium. PREDICTOR: We used 19 kidney SNPs and 64 vascular SNPs. OUTCOMES & MEASUREMENTS: Vascular outcomes tested were blood pressure, coronary artery disease, carotid intima-media thickness, pulse wave velocity, retinal venular caliber, and brain white matter lesions. Kidney outcomes were estimated glomerular filtration rate and albuminuria. RESULTS: In general, we found that kidney disease variants were not associated with vascular phenotypes (127 of 133 tests were nonsignificant). The one exception was rs653178 near SH2B3 (SH2B adaptor protein 3), which showed direction-consistent association with systolic (P = 9.3 ×10(-10)) and diastolic (P = 1.6 ×10(-14)) blood pressure and coronary artery disease (P = 2.2 ×10(-6)), all previously reported. Similarly, the 64 SNPs associated with vascular phenotypes were not associated with kidney phenotypes (187 of 192 tests were nonsignificant), with the exception of 2 high-correlated SNPs at the SH2B3 locus (P = 1.06 ×10(-07) and P = 7.05 ×10(-08)). LIMITATIONS: The combined effect size of the SNPs for kidney and vascular outcomes may be too low to detect shared genetic associations. CONCLUSIONS: Overall, although we confirmed one locus (SH2B3) as associated with both kidney and cardiovascular disease, our primary findings suggest that there is little overlap between kidney and cardiovascular disease risk variants in the overall population. The reciprocal risks of kidney and cardiovascular disease may not be genetically mediated, but rather a function of the disease milieu itself.

Meta-analysis of 28,141 individuals identifies common variants within five new loci that influence uric acid concentrations.

Elevated serum uric acid levels cause gout and are a risk factor for cardiovascular disease and diabetes. To investigate the polygenetic basis of serum uric acid levels, we conducted a meta-analysis of genome-wide association scans from 14 studies totalling 28,141 participants of European descent, resulting in identification of 954 SNPs distributed across nine loci that exceeded the threshold of genome-wide significance, five of which are novel. Overall, the common variants associated with serum uric acid levels fall in the following nine regions: SLC2A9 (p = 5.2x10(-201)), ABCG2 (p = 3.1x10(-26)), SLC17A1 (p = 3.0x10(-14)), SLC22A11 (p = 6.7x10(-14)), SLC22A12 (p = 2.0x10(-9)), SLC16A9 (p = 1.1x10(-8)), GCKR (p = 1.4x10(-9)), LRRC16A (p = 8.5x10(-9)), and near PDZK1 (p = 2.7x10(-9)). Identified variants were analyzed for gender differences. We found that the minor allele for rs734553 in SLC2A9 has greater influence in lowering uric acid levels in women and the minor allele of rs2231142 in ABCG2 elevates uric acid levels more strongly in men compared to women. To further characterize the identified variants, we analyzed their association with a panel of metabolites. rs12356193 within SLC16A9 was associated with DL-carnitine (p = 4.0x10(-26)) and propionyl-L-carnitine (p = 5.0x10(-8)) concentrations, which in turn were associated with serum UA levels (p = 1.4x10(-57) and p = 8.1x10(-54), respectively), forming a triangle between SNP, metabolites, and UA levels. Taken together, these associations highlight additional pathways that are important in the regulation of serum uric acid levels and point toward novel potential targets for pharmacological intervention to prevent or treat hyperuricemia. In addition, these findings strongly support the hypothesis that transport proteins are key in regulating serum uric acid levels.

Effects of sample size on the performance of species distribution models

A wide range of modelling algorithms is used by ecologists, conservation practitioners, and others to predict species ranges from point locality data. Unfortunately, the amount of data available is limited for many taxa and regions, making it essential to quantify the sensitivity of these algorithms to sample size. This is the first study to address this need by rigorously evaluating a broad suite of algorithms with independent presence-absence data from multiple species and regions. We evaluated predictions from 12 algorithms for 46 species (from six different regions of the world) at three sample sizes (100, 30, and 10 records). We used data from natural history collections to run the models, and evaluated the quality of model predictions with area under the receiver operating characteristic curve (AUC). With decreasing sample size, model accuracy decreased and variability increased across species and between models. Novel modelling methods that incorporate both interactions between predictor variables and complex response shapes (i.e. GBM, MARS-INT, BRUTO) performed better than most methods at large sample sizes but not at the smallest sample sizes. Other algorithms were much less sensitive to sample size, including an algorithm based on maximum entropy (MAXENT) that had among the best predictive power across all sample sizes. Relative to other algorithms, a distance metric algorithm (DOMAIN) and a genetic algorithm (OM-GARP) had intermediate performance at the largest sample size and among the best performance at the lowest sample size. No algorithm predicted consistently well with small sample size (n < 30) and this should encourage highly conservative use of predictions based on small sample size and restrict their use to exploratory modelling.

Effects of bisoprolol fumarate on left ventricular size, function, and exercise capacity in patients with heart failure: analysis with magnetic resonance myocardial tagging.

BACKGROUND: Recent data suggest that beta-blockers can be beneficial in subgroups of patients with chronic heart failure (CHF). For metoprolol and carvedilol, an increase in ejection fraction has been shown and favorable effects on the myocardial remodeling process have been reported in some studies. We examined the effects of bisoprolol fumarate on exercise capacity and left ventricular volume with magnetic resonance imaging (MRI) and applied a novel high-resolution MRI tagging technique to determine myocardial rotation and relaxation velocity. METHODS: Twenty-eight patients (mean age, 57 +/- 11 years; mean ejection fraction, 26 +/- 6%) were randomized to bisoprolol fumarate (n = 13) or to placebo therapy (n = 15). The dosage of the drugs was titrated to match that of the the Cardiac Insufficiency Bisoprolol Study protocol. Hemodynamic and gas exchange responses to exercise, MRI measurements of left ventricular end-systolic and end-diastolic volumes and ejection fraction, and left ventricular rotation and relaxation velocities were measured before the administration of the drug and 6 and 12 months later. RESULTS: After 1 year, heart rate was reduced in the bisoprolol fumarate group both at rest (81 +/- 12 before therapy versus 61 +/- 11 after therapy; P <.01) and peak exercise (144 +/- 20 before therapy versus 127 +/- 17 after therapy; P <.01), which indicated a reduction in sympathetic drive. No differences were observed in heart rate responses in the placebo group. No differences were observed within or between groups in peak oxygen uptake, although work rate achieved was higher (117.9 +/- 36 watts versus 146.1 +/- 33 watts; P <.05) and exercise time tended to be higher (9.1 +/- 1.7 minutes versus 11.4 +/- 2.8 minutes; P =.06) in the bisoprolol fumarate group. A trend for a reduction in left ventricular end-diastolic volume (-54 mL) and left ventricular end-systolic volume (-62 mL) in the bisoprolol fumarate group occurred after 1 year. Ejection fraction was higher in the bisoprolol fumarate group (25.0 +/- 7 versus 36.2 +/- 9%; P <.05), and the placebo group remained unchanged. Most changes in volume and ejection fraction occurred during the latter 6 months of treatment. With myocardial tagging, insignificant reductions in left ventricular rotation velocity were observed in both groups, whereas relaxation velocity was reduced only after bisoprolol fumarate therapy (by 39%; P <.05). CONCLUSION: One year of bisoprolol fumarate therapy resulted in an improvement in exercise capacity, showed trends for reductions in end-diastolic and end-systolic volumes, increased ejection fraction, and significantly reduced relaxation velocity. Although these results generally confirm the beneficial effects of beta-blockade in patients with chronic heart failure, they show differential effects on systolic and diastolic function.

Novel method to estimate the phenotypic variation explained by genome-wide association studies reveals large fraction of the missing heritability.

Genome-wide association studies (GWAS) are conducted with the promise to discover novel genetic variants associated with diverse traits. For most traits, associated markers individually explain just a modest fraction of the phenotypic variation, but their number can well be in the hundreds. We developed a maximum likelihood method that allows us to infer the distribution of associated variants even when many of them were missed by chance. Compared to previous approaches, the novelty of our method is that it (a) does not require having an independent (unbiased) estimate of the effect sizes; (b) makes use of the complete distribution of P-values while allowing for the false discovery rate; (c) takes into account allelic heterogeneity and the SNP pruning strategy. We applied our method to the latest GWAS meta-analysis results of the GIANT consortium. It revealed that while the explained variance of genome-wide (GW) significant SNPs is around 1% for waist-hip ratio (WHR), the observed P-values provide evidence for the existence of variants explaining 10% (CI=[8.5-11.5%]) of the phenotypic variance in total. Similarly, the total explained variance likely to exist for height is estimated to be 29% (CI=[28-30%]), three times higher than what the observed GW significant SNPs give rise to. This methodology also enables us to predict the benefit of future GWA studies that aim to reveal more associated genetic markers via increased sample size.

Am I shrinking? On clothing size and body size.

Comment on: Hughes LA, Schouten LJ, Goldbohm RA, van den Brandt PA, Weijenberg MP. Self-reported clothing size as a proxy measure for body size. Epidemiology. 2009 Sep;20(5):673-6. PMID: 19451821

Impact of phenotype definition on genome-wide association signals: empirical evaluation in human immunodeficiency virus type 1 infection.

Discussion on improving the power of genome-wide association studies to identify candidate variants and genes is generally centered on issues of maximizing sample size; less attention is given to the role of phenotype definition and ascertainment. The authors used genome-wide data from patients infected with human immunodeficiency virus type 1 (HIV-1) to assess whether differences in type of population (622 seroconverters vs. 636 seroprevalent subjects) or the number of measurements available for defining the phenotype resulted in differences in the effect sizes of associations between single nucleotide polymorphisms and the phenotype, HIV-1 viral load at set point. The effect estimate for the top 100 single nucleotide polymorphisms was 0.092 (95% confidence interval: 0.074, 0.110) log(10) viral load (log(10) copies of HIV-1 per mL of blood) greater in seroconverters than in seroprevalent subjects. The difference was even larger when the authors focused on chromosome 6 variants (0.153 log(10) viral load) or on variants that achieved genome-wide significance (0.232 log(10) viral load). The estimates of the genetic effects tended to be slightly larger when more viral load measurements were available, particularly among seroconverters and for variants that achieved genome-wide significance. Differences in phenotype definition and ascertainment may affect the estimated magnitude of genetic effects and should be considered in optimizing power for discovering new associations.

Long-term antiretroviral treatment initiated at primary HIV-1 infection affects the size, composition, and decay kinetics of the reservoir of HIV-1-infected CD4 T cells.

Initiation of antiretroviral therapy during the earliest stages of HIV-1 infection may limit the seeding of a long-lasting viral reservoir, but long-term effects of early antiretroviral treatment initiation remain unknown. Here, we analyzed immunological and virological characteristics of nine patients who started antiretroviral therapy at primary HIV-1 infection and remained on suppressive treatment for >10 years; patients with similar treatment duration but initiation of suppressive therapy during chronic HIV-1 infection served as controls. We observed that independently of the timing of treatment initiation, HIV-1 DNA in CD4 T cells decayed primarily during the initial 3 to 4 years of treatment. However, in patients who started antiretroviral therapy in early infection, this decay occurred faster and was more pronounced, leading to substantially lower levels of cell-associated HIV-1 DNA after long-term treatment. Despite this smaller size, the viral CD4 T cell reservoir in persons with early treatment initiation consisted more dominantly of the long-lasting central-memory and T memory stem cells. HIV-1-specific T cell responses remained continuously detectable during antiretroviral therapy, independently of the timing of treatment initiation. Together, these data suggest that early HIV-1 treatment initiation, even when continued for >10 years, is unlikely to lead to viral eradication, but the presence of low viral reservoirs and durable HIV-1 T cell responses may make such patients good candidates for future interventional studies aiming at HIV-1 eradication and cure. IMPORTANCE: Antiretroviral therapy can effectively suppress HIV-1 replication to undetectable levels; however, HIV-1 can persist despite treatment, and viral replication rapidly rebounds when treatment is discontinued. This is mainly due to the presence of latently infected CD4 T cells, which are not susceptible to antiretroviral drugs. Starting treatment in the earliest stages of HIV-1 infection can limit the number of these latently infected cells, raising the possibility that these viral reservoirs are naturally eliminated if suppressive antiretroviral treatment is continued for extremely long periods of time. Here, we analyzed nine patients who started on antiretroviral therapy within the earliest weeks of the disease and continued treatment for more than 10 years. Our data show that early treatment accelerated the decay of infected CD4 T cells and led to very low residual levels of detectable HIV-1 after long-term therapy, levels that were otherwise detectable in patients who are able to maintain a spontaneous, drug-free control of HIV-1 replication. Thus, long-term antiretroviral treatment started during early infection cannot eliminate HIV-1, but the reduced reservoirs of HIV-1 infected cells in such patients may increase their chances to respond to clinical interventions aiming at inducing a drug-free remission of HIV-1 infection.

Computational analysis of the genetic and environmental contributions to disease-related human phenotypes: From predicting adverse lipid response of HIV patients receiving antiretroviral therapy to genome-wide association studies of cardiovascular and lipid related disorders

Genetic variants influence the risk to develop certain diseases or give rise to differences in drug response. Recent progresses in cost-effective, high-throughput genome-wide techniques, such as microarrays measuring Single Nucleotide Polymorphisms (SNPs), have facilitated genotyping of large clinical and population cohorts. Combining the massive genotypic data with measurements of phenotypic traits allows for the determination of genetic differences that explain, at least in part, the phenotypic variations within a population. So far, models combining the most significant variants can only explain a small fraction of the variance, indicating the limitations of current models. In particular, researchers have only begun to address the possibility of interactions between genotypes and the environment. Elucidating the contributions of such interactions is a difficult task because of the large number of genetic as well as possible environmental factors.In this thesis, I worked on several projects within this context. My first and main project was the identification of possible SNP-environment interactions, where the phenotypes were serum lipid levels of patients from the Swiss HIV Cohort Study (SHCS) treated with antiretroviral therapy. Here the genotypes consisted of a limited set of SNPs in candidate genes relevant for lipid transport and metabolism. The environmental variables were the specific combinations of drugs given to each patient over the treatment period. My work explored bioinformatic and statistical approaches to relate patients' lipid responses to these SNPs, drugs and, importantly, their interactions. The goal of this project was to improve our understanding and to explore the possibility of predicting dyslipidemia, a well-known adverse drug reaction of antiretroviral therapy. Specifically, I quantified how much of the variance in lipid profiles could be explained by the host genetic variants, the administered drugs and SNP-drug interactions and assessed the predictive power of these features on lipid responses. Using cross-validation stratified by patients, we could not validate our hypothesis that models that select a subset of SNP-drug interactions in a principled way have better predictive power than the control models using "random" subsets. Nevertheless, all models tested containing SNP and/or drug terms, exhibited significant predictive power (as compared to a random predictor) and explained a sizable proportion of variance, in the patient stratified cross-validation context. Importantly, the model containing stepwise selected SNP terms showed higher capacity to predict triglyceride levels than a model containing randomly selected SNPs. Dyslipidemia is a complex trait for which many factors remain to be discovered, thus missing from the data, and possibly explaining the limitations of our analysis. In particular, the interactions of drugs with SNPs selected from the set of candidate genes likely have small effect sizes which we were unable to detect in a sample of the present size (<800 patients).In the second part of my thesis, I performed genome-wide association studies within the Cohorte Lausannoise (CoLaus). I have been involved in several international projects to identify SNPs that are associated with various traits, such as serum calcium, body mass index, two-hour glucose levels, as well as metabolic syndrome and its components. These phenotypes are all related to major human health issues, such as cardiovascular disease. I applied statistical methods to detect new variants associated with these phenotypes, contributing to the identification of new genetic loci that may lead to new insights into the genetic basis of these traits. This kind of research will lead to a better understanding of the mechanisms underlying these pathologies, a better evaluation of disease risk, the identification of new therapeutic leads and may ultimately lead to the realization of "personalized" medicine.

Investigating the relationship between pollination strategies and the size-advantage model in zoophilous plants using the reproductive biology of Arum cylindraceum and other European Arum species as case studies

The size-advantage model (SAM) explains the temporal variation of energetic investment on reproductive structures (i.e. male and female gametes and reproductive organs) in long-lived hermaphroditic plants and animals. It proposes that an increase in the resources available to an organism induces a higher relative investment on the most energetically costly sexual structures. In plants, pollination interactions are known to play an important role in the evolution of floral features. Because the SAM directly concerns flower characters, pollinators are expected to have a strong influence on the application of the model. This hypothesis, however, has never been tested. Here, we investigate whether the identity and diversity of pollinators can be used as a proxy to predict the application of the SAM in exclusive zoophilous plants. We present a new approach to unravel the dynamics of the model and test it on several widespread Arum (Araceae) species. By identifying the species composition, abundance and spatial variation of arthropods trapped in inflorescences, we show that some species (i.e. A. cylindraceum and A. italicum) display a generalist reproductive strategy, relying on the exploitation of a low number of dipterans, in contrast to the pattern seen in the specialist A. maculatum (pollinated specifically by two fly species only). Based on the model presented here, the application of the SAM is predicted for the first two and not expected in the latter species, those predictions being further confirmed by allometric measures. We here demonstrate that while an increase in the female zone occurs in larger inflorescences of generalist species, this does not happen in species demonstrating specific pollinators. This is the first time that this theory is both proposed and empirically tested in zoophilous plants. Its overall biological importance is discussed through its application in other non-Arum systems.