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.

APOBEC3G genetic variants and their influence on the progression to AIDS.

The cytosine deaminase APOBEC3G, in the absence of the human immunodeficiency virus type 1 (HIV-1) accessory gene HIV-1 viral infectivity factor (vif), inhibits viral replication by introducing G-->A hypermutation in the newly synthesized HIV-1 DNA negative strand. We tested the hypothesis that genetic variants of APOBEC3G may modify HIV-1 transmission and disease progression. Single nucleotide polymorphisms were identified in the promoter region (three), introns (two), and exons (two). Genotypes were determined for 3,073 study participants enrolled in six HIV-AIDS prospective cohorts. One codon-changing variant, H186R in exon 4, was polymorphic in African Americans (AA) (f = 37%) and rare in European Americans (f < 3%) or Europeans (f = 5%). For AA, the variant allele 186R was strongly associated with decline in CD4 T cells (CD4 slope on square root scale: -1.86, P = 0.009), The 186R allele was also associated with accelerated progression to AIDS-defining conditions in AA. The in vitro antiviral activity of the 186R enzyme was not inferior to that of the common H186 variant. These studies suggest that there may be a modifying role of variants of APOBEC3G on HIV-1 disease progression that warrants further investigation.

Prospective study of insulin-like growth factor-I, insulin-like growth factor-binding protein 3, genetic variants in the IGF1 and IGFBP3 genes and risk of coronary artery disease.

Although experimental studies have suggested that insulin-like growth factor I (IGF-I) and its binding protein IGFBP-3 might have a role in the aetiology of coronary artery disease (CAD), the relevance of circulating IGFs and their binding proteins in the development of CAD in human populations is unclear. We conducted a nested case-control study, with a mean follow-up of six years, within the EPIC-Norfolk cohort to assess the association between circulating levels of IGF-I and IGFBP-3 and risk of CAD in up to 1,013 cases and 2,055 controls matched for age, sex and study enrolment date. After adjustment for cardiovascular risk factors, we found no association between circulating levels of IGF-I or IGFBP-3 and risk of CAD (odds ratio: 0.98 (95% Cl 0.90-1.06) per 1 SD increase in circulating IGF-I; odds ratio: 1.02 (95% Cl 0.94-1.12) for IGFBP-3). We examined associations between tagging single nucleotide polymorphisms (tSNPs) at the IGF1 and IGFBP3 loci and circulating IGF-I and IGFBP-3 levels in up to 1,133 cases and 2,223 controls and identified three tSNPs (rs1520220, rs3730204, rs2132571) that showed independent association with either circulating IGF-I or IGFBP-3 levels. In an assessment of 31 SNPs spanning the IGF1 or IGFBP3 loci, none were associated with risk of CAD in a meta-analysis that included EPIC-Norfolk and eight additional studies comprising up to 9,319 cases and 19,964 controls. Our results indicate that IGF-I and IGFBP-3 are unlikely to be importantly involved in the aetiology of CAD in human populations.

Effect of neighborhood and plot size on experiments with multiple-harvest oleraceous crops

The objective of this work was to determine the efficiency of the Papadakis method on the quality evaluation of experiments with multiple-harvest oleraceous crops, and on the estimate of the covariate and the ideal plot size. Data from nine uniformity trials (five with bean pod, two with zucchini, and two with sweet pepper) and from one experiment with treatments (with sweet pepper) were used. Through the uniformity trials, the best way to calculate the covariate was defined and the optimal plot size was calculated. In the experiment with treatments, analyses of variance and covariance were performed, in which the covariate was calculated by the Papadakis method, and experimental precision was evaluated based on four statistics. The use of analysis of covariance with the covariate obtained by the Papadakis method increases the quality of experiments with multiple-harvest oleraceous crops and allows the use of smaller plot sizes. The best covariate is the one that considers a neighboring plot of each side of the reference plot.

Population genomics of eusocial insects: the costs of a vertebrate-like effective population size

The evolution of reproductive division of labour and social life in social insects has lead to the emergence of several life-history traits and adaptations typical of larger organisms: social insect colonies can reach masses of several kilograms, they start reproducing only when they are several years old, and can live for decades. These features and the monopolization of reproduction by only one or few individuals in a colony should affect molecular evolution by reducing the effective population size. We tested this prediction by analysing genome-wide patterns of coding sequence polymorphism and divergence in eusocial vs. noneusocial insects based on newly generated RNA-seq data. We report very low amounts of genetic polymorphism and an elevated ratio of nonsynonymous to synonymous changes - a marker of the effective population size - in four distinct species of eusocial insects, which were more similar to vertebrates than to solitary insects regarding molecular evolutionary processes. Moreover, the ratio of nonsynonymous to synonymous substitutions was positively correlated with the level of social complexity across ant species. These results are fully consistent with the hypothesis of a reduced effective population size and an increased genetic load in eusocial insects, indicating that the evolution of social life has important consequences at both the genomic and population levels.

Mapping genetic variants associated with beta-adrenergic responses in inbred mice.

β-blockers and β-agonists are primarily used to treat cardiovascular diseases. Inter-individual variability in response to both drug classes is well recognized, yet the identity and relative contribution of the genetic players involved are poorly understood. This work is the first genome-wide association study (GWAS) addressing the values and susceptibility of cardiovascular-related traits to a selective β(1)-blocker, Atenolol (ate), and a β-agonist, Isoproterenol (iso). The phenotypic dataset consisted of 27 highly heritable traits, each measured across 22 inbred mouse strains and four pharmacological conditions. The genotypic panel comprised 79922 informative SNPs of the mouse HapMap resource. Associations were mapped by Efficient Mixed Model Association (EMMA), a method that corrects for the population structure and genetic relatedness of the various strains. A total of 205 separate genome-wide scans were analyzed. The most significant hits include three candidate loci related to cardiac and body weight, three loci for electrocardiographic (ECG) values, two loci for the susceptibility of atrial weight index to iso, four loci for the susceptibility of systolic blood pressure (SBP) to perturbations of the β-adrenergic system, and one locus for the responsiveness of QTc (p<10(-8)). An additional 60 loci were suggestive for one or the other of the 27 traits, while 46 others were suggestive for one or the other drug effects (p<10(-6)). Most hits tagged unexpected regions, yet at least two loci for the susceptibility of SBP to β-adrenergic drugs pointed at members of the hypothalamic-pituitary-thyroid axis. Loci for cardiac-related traits were preferentially enriched in genes expressed in the heart, while 23% of the testable loci were replicated with datasets of the Mouse Phenome Database (MPD). Altogether these data and validation tests indicate that the mapped loci are relevant to the traits and responses studied.