Application of Multi-SNP Approaches Bayesian LASSO and AUC-RF to Detect Main Effects of Inflammatory-Gene Variants Associated with Bladder Cancer Risk

The relationship between inflammation and cancer is well established in several tumor types, including bladder cancer. We performed an association study between 886 inflammatory-gene variants and bladder cancer risk in 1,047 cases and 988 controls from the Spanish Bladder Cancer (SBC)/EPICURO Study. A preliminary exploration with the widely used univariate logistic regression approach did not identify any significant SNP after correcting for multiple testing. We further applied two more comprehensive methods to capture the complexity of bladder cancer genetic susceptibility: Bayesian Threshold LASSO (BTL), a regularized regression method, and AUC-Random Forest, a machine-learning algorithm. Both approaches explore the joint effect of markers. BTL analysis identified a signature of 37 SNPs in 34 genes showing an association with bladder cancer. AUC-RF detected an optimal predictive subset of 56 SNPs. 13 SNPs were identified by both methods in the total population. Using resources from the Texas Bladder Cancer study we were able to replicate 30% of the SNPs assessed. The associations between inflammatory SNPs and bladder cancer were reexamined among non-smokers to eliminate the effect of tobacco, one of the strongest and most prevalent environmental risk factor for this tumor. A 9 SNP-signature was detected by BTL. Here we report, for the first time, a set of SNP in inflammatory genes jointly associated with bladder cancer risk. These results highlight the importance of the complex structure of genetic susceptibility associated with cancer risk.

Copy number variants, diseases and gene expression.

Copy number variation (CNV) has recently gained considerable interest as a source of genetic variation likely to play a role in phenotypic diversity and evolution. Much effort has been put into the identification and mapping of regions that vary in copy number among seemingly normal individuals in humans and a number of model organisms, using bioinformatics or hybridization-based methods. These have allowed uncovering associations between copy number changes and complex diseases in whole-genome association studies, as well as identify new genomic disorders. At the genome-wide scale, however, the functional impact of CNV remains poorly studied. Here we review the current catalogs of CNVs, their association with diseases and how they link genotype and phenotype. We describe initial evidence which revealed that genes in CNV regions are expressed at lower and more variable levels than genes mapping elsewhere, and also that CNV not only affects the expression of genes varying in copy number, but also have a global influence on the transcriptome. Further studies are warranted for complete cataloguing and fine mapping of CNVs, as well as to elucidate the different mechanisms by which they influence gene expression.

No interactions between previously associated 2-hour glucose gene variants and physical activity or BMI on 2-hour glucose levels.

Gene-lifestyle interactions have been suggested to contribute to the development of type 2 diabetes. Glucose levels 2 h after a standard 75-g glucose challenge are used to diagnose diabetes and are associated with both genetic and lifestyle factors. However, whether these factors interact to determine 2-h glucose levels is unknown. We meta-analyzed single nucleotide polymorphism (SNP) × BMI and SNP × physical activity (PA) interaction regression models for five SNPs previously associated with 2-h glucose levels from up to 22 studies comprising 54,884 individuals without diabetes. PA levels were dichotomized, with individuals below the first quintile classified as inactive (20%) and the remainder as active (80%). BMI was considered a continuous trait. Inactive individuals had higher 2-h glucose levels than active individuals (β = 0.22 mmol/L [95% CI 0.13-0.31], P = 1.63 × 10(-6)). All SNPs were associated with 2-h glucose (β = 0.06-0.12 mmol/allele, P ≤ 1.53 × 10(-7)), but no significant interactions were found with PA (P > 0.18) or BMI (P ≥ 0.04). In this large study of gene-lifestyle interaction, we observed no interactions between genetic and lifestyle factors, both of which were associated with 2-h glucose. It is perhaps unlikely that top loci from genome-wide association studies will exhibit strong subgroup-specific effects, and may not, therefore, make the best candidates for the study of interactions.

Natural genetic variation in Arabidopsis: tools, traits and prospects for evolutionary ecology.

BACKGROUND: The model plant Arabidopsis thaliana (Arabidopsis) shows a wide range of genetic and trait variation among wild accessions. Because of its unparalleled biological and genomic resources, the potential of Arabidopsis for molecular genetic analysis of this natural variation has increased dramatically in recent years. SCOPE: Advanced genomics has accelerated molecular phylogenetic analysis and gene identification by quantitative trait loci (QTL) mapping and/or association mapping in Arabidopsis. In particular, QTL mapping utilizing natural accessions is now becoming a major strategy of gene isolation, offering an alternative to artificial mutant lines. Furthermore, the genomic information is used by researchers to uncover the signature of natural selection acting on the genes that contribute to phenotypic variation. The evolutionary significance of such genes has been evaluated in traits such as disease resistance and flowering time. However, although molecular hallmarks of selection have been found for the genes in question, a corresponding ecological scenario of adaptive evolution has been difficult to prove. Ecological strategies, including reciprocal transplant experiments and competition experiments, and utilizing near-isogenic lines of alleles of interest will be a powerful tool to measure the relative fitness of phenotypic and/or allelic variants. CONCLUSIONS: As the plant model organism, Arabidopsis provides a wealth of molecular background information for evolutionary genetics. Because genetic diversity between and within Arabidopsis populations is much higher than anticipated, combining this background information with ecological approaches might well establish Arabidopsis as a model organism for plant evolutionary ecology.

Genetic variation in GIPR influences the glucose and insulin responses to an oral glucose challenge.

Glucose levels 2 h after an oral glucose challenge are a clinical measure of glucose tolerance used in the diagnosis of type 2 diabetes. We report a meta-analysis of nine genome-wide association studies (n = 15,234 nondiabetic individuals) and a follow-up of 29 independent loci (n = 6,958-30,620). We identify variants at the GIPR locus associated with 2-h glucose level (rs10423928, beta (s.e.m.) = 0.09 (0.01) mmol/l per A allele, P = 2.0 x 10(-15)). The GIPR A-allele carriers also showed decreased insulin secretion (n = 22,492; insulinogenic index, P = 1.0 x 10(-17); ratio of insulin to glucose area under the curve, P = 1.3 x 10(-16)) and diminished incretin effect (n = 804; P = 4.3 x 10(-4)). We also identified variants at ADCY5 (rs2877716, P = 4.2 x 10(-16)), VPS13C (rs17271305, P = 4.1 x 10(-8)), GCKR (rs1260326, P = 7.1 x 10(-11)) and TCF7L2 (rs7903146, P = 4.2 x 10(-10)) associated with 2-h glucose. Of the three newly implicated loci (GIPR, ADCY5 and VPS13C), only ADCY5 was found to be associated with type 2 diabetes in collaborating studies (n = 35,869 cases, 89,798 controls, OR = 1.12, 95% CI 1.09-1.15, P = 4.8 x 10(-18)).

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.

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.

Inheritance of late flowering in natural variants of soybean cultivars under short-day conditions

The objective of this work was to determine the inheritance of the long juvenile period trait in natural variants of the Doko, BR 9 (Savana), Davis, Embrapa 1 (IAS 5RC), and BR 16 soybean cultivars. Complete diallel crosses were made between the Doko and BR 16 cultivars and their variants. A 3:1 segregation ratio was observed in the F2 populations of the 'Doko' x Doko-18T, 'Doko' x Doko-Milionária, 'Davis' x São Carlos, and 'BR 9 (Savana)' x MABR92-836 (Savanão) crosses, indicating that the long juvenile period trait is controlled by a pair of recessive genes. The difference in late flowering between the Doko cultivar and both of its variants was caused by a recessive spontaneous mutation at the same genetic locus. However, the variants Doko-18T and Doko-Milionária are identical mutants that share a pair of genes that control the long juvenile period under short-day conditions. These mutants can be used in breeding programs to develop cultivars adapted to low-latitude tropical regions.

Genetic associations in the detection of QTLs for wheat spike-related traits

The objective of this work was to assess the genetic diversity and population structure of wheat genotypes, to detect significant and stable genetic associations, as well as to evaluate the efficiency of statistical models to identify chromosome regions responsible for the expression of spike-related traits. Eight important spike characteristics were measured during five growing seasons in Serbia. A set of 30 microsatellite markers positioned near important agronomic loci was used to evaluate genetic diversity, resulting in a total of 349 alleles. The marker-trait associations were analyzed using the general linear and mixed linear models. The results obtained for number of allelic variants per locus (11.5), average polymorphic information content value (0.68), and average gene diversity (0.722) showed that the exceptional level of polymorphism in the genotypes is the main requirement for association studies. The population structure estimated by model-based clustering distributed the genotypes into six subpopulations according to log probability of data. Significant and stable associations were detected on chromosomes 1B, 2A, 2B, 2D, and 6D, which explained from 4.7 to 40.7% of total phenotypic variations. The general linear model identified a significantly larger number of marker-trait associations (192) than the mixed linear model (76). The mixed linear model identified nine markers associated to six traits.

Modification of the monomer composition of polyhydroxyalkanoate synthesized in Saccharomyces cerevisiae expressing variants of the beta-oxidation-associated multifunctional enzyme.

Expression by Saccharomyces cerevisiae of a polyhydroxyalkanoate (PHA) synthase modified at the carboxy end by the addition of a peroxisome targeting signal derived from the last 34 amino acids of the Brassica napus isocitrate lyase (ICL) and containing the terminal tripeptide Ser-Arg-Met resulted in the synthesis of PHA. The ability of the terminal peptide Ser-Arg-Met and of the 34-amino-acid peptide from the B. napus ICL to target foreign proteins to the peroxisome of S. cerevisiae was demonstrated with green fluorescent protein fusions. PHA synthesis was found to be dependent on the presence of both the enzymes generating the beta-oxidation intermediate 3-hydroxyacyl-coenzyme A (3-hydroxyacyl-[CoA]) and the peroxin-encoding PEX5 gene, demonstrating the requirement for a functional peroxisome and a beta-oxidation cycle for PHA synthesis. Using a variant of the S. cerevisiae beta-oxidation multifunctional enzyme with a mutation inactivating the B domain of the R-3-hydroxyacyl-CoA dehydrogenase, it was possible to modify the PHA monomer composition through an increase in the proportion of the short-chain monomers of five and six carbons.

Genome-wide meta-analysis identifies 11 new loci for anthropometric traits and provides insights into genetic architecture.

Approaches exploiting trait distribution extremes may be used to identify loci associated with common traits, but it is unknown whether these loci are generalizable to the broader population. In a genome-wide search for loci associated with the upper versus the lower 5th percentiles of body mass index, height and waist-to-hip ratio, as well as clinical classes of obesity, including up to 263,407 individuals of European ancestry, we identified 4 new loci (IGFBP4, H6PD, RSRC1 and PPP2R2A) influencing height detected in the distribution tails and 7 new loci (HNF4G, RPTOR, GNAT2, MRPS33P4, ADCY9, HS6ST3 and ZZZ3) for clinical classes of obesity. Further, we find a large overlap in genetic structure and the distribution of variants between traits based on extremes and the general population and little etiological heterogeneity between obesity subgroups.

Organisation of genetic variation in multinucleate arbuscular mycorrhizal fungi

Les Champignons Endomycorhiziens Arbusculaires (CEA) forment une symbiose racinaire avec environ 80% des espèces connues de plantes vasculaires. Ils occupent une position écologique très importante liée aux bénéfices qu'ils confèrent aux plantes. Des études moléculaires effectuées sur des gènes ribosomaux ont révélé un très grand polymorphisme, tant à l'intérieur des espèces qu'entre celles-ci. Ces champignons étant coenocytiques et multinucléés, l'organisation de cette variabilité génétique intraspécifique pourrait avoir différentes origines. Ce travail se propose d'examiner l'organisation et l'évolution de cette variabilité. Sur la base de fossiles, l'existence des CEA remonte à au moins 450 millions d'années. Cette symbiose peut donc être considérée comme ancienne. Les premières données moléculaires n'indiquant pas de reproduction sexuée, une hypothèse fut élaborée stipulant que les CEA seraient des asexués ancestraux. La première partie de cette thèse (chapitre 2) met en évidence l'existence de recombinaison dans différents CEA mais montre également que celle-ci est insuffisante pour purger les mutations accumulées. La reproduction étant essentiellement asexuée, on peut prédire que les nombreux noyaux ont probablement divergé génétiquement. En collaboration avec M. Hijri nous avons pu vérifier cette hypothèse (chapitre 2). Dans le chapitre 3 j'ai cherché à comprendre si le polymorphisme était également présent dans une population naturelle du CEA Glomus intraradices au niveau intraspécifique, ce qui n'avait encore jamais été examiné. En comparant les empreintes génétiques d'individus obtenus chacun à partir d'une spore mise en culture, j'ai clairement démontré que d'importantes différences génétiques existent entre ceux-ci. Un résultat similaire, portant sur des traits quantitatifs d'individus de la même population, a été trouvé par A. Koch. Les deux études en ensemble montre que le polymorphisme génétique dans cette population est suffisamment grand pour être important au niveau écologique. Dans le chapitre 4, j'ai cherché a examiner le polymorphisme des séquences du gène BiP au sein d'un individu. C'est la première étude qui examine la diversité génétique du génome de CEA avec un autre marqueur que l'ADN ribosomique. J'ai trouvé 31 types de séquences différentes du gène BiP issu d'un isolat de G. intraradices mis en culture à partir d'une seule spore. Cette variation n'était pas restreinte à des zones sélectivement neutres du BiP. Mes résultats montrent qu'il y a un grand nombre de variants non-fonctionnels, proportionnellement au faible nombre de copies attendues par noyau. Ceci va dans le sens d'une partition de l'information génétique entre les noyaux.<br/><br/>Arbuscular mycorrhizal fungi (AMF) are root symbionts with about 80% of all known species of vascular land plants. AMF are ecologically important because of the benefits that they confer to plants. Molecular studies on AMF showed that rDNA sequences were highly variable between species and within species. Because AMF are coenocytic and multinucleate there are several possibilities how this intraspecific genetic variation could be organized. Therefore, the organization and evolution of this variation in AMF were investigated in the present work. Based on fossil records the AMF symbiosis has existed for 450 Million years and is therefore considered ancient. First molecular data indicated no evident sexual reproduction and gave rise to the hypothesis that AMF might be ancient asexuals. The first part of this thesis (Chapter 2) shows evidence for recombination in different AMF but also indicates that it has not been frequent enough to purge accumulated mutations. Given asexual reproduction, it has been predicted that the many nuclei in AMF should diverge leading to genetically different nuclei. This hypothesis has been confirmed by an experiment of M. Hijri and is also included in chapter 2 as the results were published together. In chapter 3 I then investigated whether intraspecific genetic variation also exists in a field population of the AMF Glomus intraradices. Comparing genetic fingerprints of individuals derived from single spores I could clearly show that large genetic differences exist. A similar result, based on quantitative genetic traits, was found for the same population by A. Koch. The two studies taken together show that the genetic variation observed in the population is high enough to be of ecological relevance. Lastly, in chapter 4, I investigated within individual genetic variation among BiP gene sequences. It is the first study that has analyzed genetic diversity in the AMF genome in a region of DNA other than rDNA. I found 31 sequence variants of the BiP gene in one G. intraradices isolate that originated from one spore. Genetic variation was not only restricted to selectively neutral parts of BiP. A high number of predicted non-functional variants compared to a likely low number of copies per nucleus indicated that functional genetic information might even be partitioned among nuclei. The results of this work contribute to our understanding of potential evolutionary strategies of ancient asexuals, they also suggest that genetic differences in a population might be ecologically relevant and they show that this variation even occurs in functional regions of the AMF genome.

Analysis of the genetic basis of periodic fever with aphthous stomatitis, pharyngitis, and cervical adenitis (PFAPA) syndrome.

PFAPA syndrome is the most common autoinflammatory syndrome in children from Western countries. In spite of its strong familial clustering, its genetic basis and inheritance pattern are still unknown. We performed a comprehensive genetic study on 68 individuals from 14 families. Linkage analysis suggested a susceptibility locus on chromosome 8, but direct molecular sequencing did not support this initial statistical finding. Exome sequencing revealed the absence of any gene that was mutated in all patients. Exhaustive screening of genes involved in other autoinflammatory syndromes or encoding components of the human inflammasome showed no DNA variants that could be linked to PFAPA molecular pathology. Among these, the previously-reported missense mutation V198M in the NLRP3 gene was clearly shown not to co-segregate with PFAPA. Our results on this relatively large cohort indicate that PFAPA syndrome is unlikely to be a monogenic condition. Moreover, none of the several genes known to be involved in inflammation or in autoinflammatory disorders seem to be relevant, alone, to its etiology, suggesting that PFAPA results from oligogenic or complex inheritance of variants in multiple disease genes and/or non-genetic factors.