Estudio Farmacogenético en Pacientes con Cáncer de Pulmón No Microcítico Avanzado tratados con Quimioterapia basada en Platino

a partir de ADN genómico obtenido de las células nucleadas de sangre periférica de 103 pacientes con Cáncer de Pulmón No Microcítico (CPNM) avanzado tratados con quimioterapia basada en platino, hemos analizado la asociación entre supervivencia y cinco SNPs (Single Nucleotide Polymorphism) pertenecientes a dos grupos de genes: i) de la via metabólica del ácido fólico (Timidilato Sintetasa (TS), Metil-tetrahidrofolato Reductasa (MTHFR) y, ii) de la vía de reparación del ADN (Excision repair cross-complemeting group 1 (ERCC1) y Xeroderma pigmentosum group D (XPD).

Global Genomic Diversity of Human Papillomavirus 6 Based on 724 Isolates and 190 Complete Genome Sequences.

Human papillomavirus type 6 (HPV6) is the major etiological agent of anogenital warts and laryngeal papillomas and has been included in both the quadrivalent and nonavalent prophylactic HPV vaccines. This study investigated the global genomic diversity of HPV6, using 724 isolates and 190 complete genomes from six continents, and the association of HPV6 genomic variants with geographical location, anatomical site of infection/disease, and gender. Initially, a 2,800-bp E5a-E5b-L1-LCR fragment was sequenced from 492/530 (92.8%) HPV6-positive samples collected for this study. Among them, 130 exhibited at least one single nucleotide polymorphism (SNP), indel, or amino acid change in the E5a-E5b-L1-LCR fragment and were sequenced in full. A global alignment and maximum likelihood tree of 190 complete HPV6 genomes (130 fully sequenced in this study and 60 obtained from sequence repositories) revealed two variant lineages, A and B, and five B sublineages: B1, B2, B3, B4, and B5. HPV6 (sub)lineage-specific SNPs and a 960-bp representative region for whole-genome-based phylogenetic clustering within the L2 open reading frame were identified. Multivariate logistic regression analysis revealed that lineage B predominated globally. Sublineage B3 was more common in Africa and North and South America, and lineage A was more common in Asia. Sublineages B1 and B3 were associated with anogenital infections, indicating a potential lesion-specific predilection of some HPV6 sublineages. Females had higher odds for infection with sublineage B3 than males. In conclusion, a global HPV6 phylogenetic analysis revealed the existence of two variant lineages and five sublineages, showing some degree of ethnogeographic, gender, and/or disease predilection in their distribution. IMPORTANCE: This study established the largest database of globally circulating HPV6 genomic variants and contributed a total of 130 new, complete HPV6 genome sequences to available sequence repositories. Two HPV6 variant lineages and five sublineages were identified and showed some degree of association with geographical location, anatomical site of infection/disease, and/or gender. We additionally identified several HPV6 lineage- and sublineage-specific SNPs to facilitate the identification of HPV6 variants and determined a representative region within the L2 gene that is suitable for HPV6 whole-genome-based phylogenetic analysis. This study complements and significantly expands the current knowledge of HPV6 genetic diversity and forms a comprehensive basis for future epidemiological, evolutionary, functional, pathogenicity, vaccination, and molecular assay development studies.

Integration of genome-wide association studies with biological knowledge identifies six novel genes related to kidney function.

In conducting genome-wide association studies (GWAS), analytical approaches leveraging biological information may further understanding of the pathophysiology of clinical traits. To discover novel associations with estimated glomerular filtration rate (eGFR), a measure of kidney function, we developed a strategy for integrating prior biological knowledge into the existing GWAS data for eGFR from the CKDGen Consortium. Our strategy focuses on single nucleotide polymorphism (SNPs) in genes that are connected by functional evidence, determined by literature mining and gene ontology (GO) hierarchies, to genes near previously validated eGFR associations. It then requires association thresholds consistent with multiple testing, and finally evaluates novel candidates by independent replication. Among the samples of European ancestry, we identified a genome-wide significant SNP in FBXL20 (P = 5.6 × 10(-9)) in meta-analysis of all available data, and additional SNPs at the INHBC, LRP2, PLEKHA1, SLC3A2 and SLC7A6 genes meeting multiple-testing corrected significance for replication and overall P-values of 4.5 × 10(-4)-2.2 × 10(-7). Neither the novel PLEKHA1 nor FBXL20 associations, both further supported by association with eGFR among African Americans and with transcript abundance, would have been implicated by eGFR candidate gene approaches. LRP2, encoding the megalin receptor, was identified through connection with the previously known eGFR gene DAB2 and extends understanding of the megalin system in kidney function. These findings highlight integration of existing genome-wide association data with independent biological knowledge to uncover novel candidate eGFR associations, including candidates lacking known connections to kidney-specific pathways. The strategy may also be applicable to other clinical phenotypes, although more testing will be needed to assess its potential for discovery in general.

Genome-wide arrays in routine diagnostics of hematological malignancies.

Over the last three decades, cytogenetic analysis of malignancies has become an integral part of disease evaluation and prediction of prognosis or responsiveness to therapy. In most diagnostic laboratories, conventional karyotyping, in conjunction with targeted fluorescence in situ hybridization analysis, is routinely performed to detect recurrent aberrations with prognostic implications. However, the genetic complexity of cancer cells requires a sensitive genome-wide analysis, enabling the detection of small genomic changes in a mixed cell population, as well as of regions of homozygosity. The advent of comprehensive high-resolution genomic tools, such as molecular karyotyping using comparative genomic hybridization or single-nucleotide polymorphism microarrays, has overcome many of the limitations of traditional cytogenetic techniques and has been used to study complex genomic lesions in, for example, leukemia. The clinical impact of the genomic copy-number and copy-neutral alterations identified by microarray technologies is growing rapidly and genome-wide array analysis is evolving into a diagnostic tool, to better identify high-risk patients and predict patients' outcomes from their genomic profiles. Here, we review the added clinical value of an array-based genome-wide screen in leukemia, and discuss the technical challenges and an interpretation workflow in applying arrays in the acquired cytogenetic diagnostic setting.

In vitro whole-genome analysis identifies a susceptibility locus for HIV-1.

Advances in large-scale analysis of human genomic variability provide unprecedented opportunities to study the genetic basis of susceptibility to infectious agents. We report here the use of an in vitro system for the identification of a locus on HSA8q24.3 associated with cellular susceptibility to HIV-1. This locus was mapped through quantitative linkage analysis using cell lines from multigeneration families, validated in vitro, and followed up by two independent association studies in HIV-positive individuals. Single nucleotide polymorphism rs2572886, which is associated with cellular susceptibility to HIV-1 in lymphoblastoid B cells and in primary T cells, was also associated with accelerated disease progression in one of two cohorts of HIV-1-infected patients. Biological analysis suggests a role of the rs2572886 region in the regulation of the LY6 family of glycosyl-phosphatidyl-inositol (GPI)-anchored proteins. Genetic analysis of in vitro cellular phenotypes provides an attractive approach for the discovery of susceptibility loci to infectious agents.

Polymorphisms in tumor necrosis factor-α increase susceptibility to intra-abdominal Candida infection in high-risk surgical ICU patients*.

OBJECTIVES: To evaluate the influence of genetic polymorphisms on the susceptibility to Candida colonization and intra-abdominal candidiasis, a blood culture-negative life-threatening infection in high-risk surgical ICU patients. DESIGN: Prospective observational cohort study. SETTING: Surgical ICUs from two University hospitals of the Fungal Infection Network of Switzerland. PATIENTS: Eighty-nine patients at high risk for intra-abdominal candidiasis (68 with recurrent gastrointestinal perforation and 21 with acute necrotizing pancreatitis). MEASUREMENTS AND MAIN RESULTS: Eighteen single-nucleotide polymorphisms in 16 genes previously associated with development of fungal infections were analyzed from patient's DNA by using an Illumina Veracode genotyping platform. Candida colonization was defined by recovery of Candida species from at least one nonsterile site by twice weekly monitoring of cultures from oropharynx, stools, urine, skin, and/or respiratory tract. A corrected colonization index greater than or equal to 0.4 defined "heavy" colonization. Intra-abdominal candidiasis was defined by the presence of clinical symptoms and signs of peritonitis or intra-abdominal abscess and isolation of Candida species either in pure or mixed culture from intraoperatively collected abdominal samples. Single-nucleotide polymorphisms in three innate immune genes were associated with development of a Candida corrected colonization index greater than or equal to 0.4 (Toll-like receptor rs4986790, hazard ratio = 3.39; 95% CI, 1.45-7.93; p = 0.005) or occurrence of intra-abdominal candidiasis (tumor necrosis factor-α rs1800629, hazard ratio = 4.31; 95% CI, 1.85-10.1; p= 0.0007; β-defensin 1 rs1800972, hazard ratio = 3.21; 95% CI, 1.36-7.59; p = 0.008). CONCLUSION: We report a strong association between the promoter rs1800629 single-nucleotide polymorphism in tumor necrosis factor-α and an increased susceptibility to intra-abdominal candidiasis in a homogenous prospective cohort of high-risk surgical ICU patients. This finding highlights the relevance of the tumor necrosis factor-α functional polymorphism in immune response to fungal pathogens. Immunogenetic profiling in patients at clinical high risk followed by targeted antifungal interventions may improve the prevention or preemptive management of this life-threatening infection.

Malaria haplotype frequency estimation.

We present a Bayesian approach for estimating the relative frequencies of multi-single nucleotide polymorphism (SNP) haplotypes in populations of the malaria parasite Plasmodium falciparum by using microarray SNP data from human blood samples. Each sample comes from a malaria patient and contains one or several parasite clones that may genetically differ. Samples containing multiple parasite clones with different genetic markers pose a special challenge. The situation is comparable with a polyploid organism. The data from each blood sample indicates whether the parasites in the blood carry a mutant or a wildtype allele at various selected genomic positions. If both mutant and wildtype alleles are detected at a given position in a multiply infected sample, the data indicates the presence of both alleles, but the ratio is unknown. Thus, the data only partially reveals which specific combinations of genetic markers (i.e. haplotypes across the examined SNPs) occur in distinct parasite clones. In addition, SNP data may contain errors at non-negligible rates. We use a multinomial mixture model with partially missing observations to represent this data and a Markov chain Monte Carlo method to estimate the haplotype frequencies in a population. Our approach addresses both challenges, multiple infections and data errors.

Association of genetic variation with systolic and diastolic blood pressure among African Americans: the Candidate Gene Association Resource study.

The prevalence of hypertension in African Americans (AAs) is higher than in other US groups; yet, few have performed genome-wide association studies (GWASs) in AA. Among people of European descent, GWASs have identified genetic variants at 13 loci that are associated with blood pressure. It is unknown if these variants confer susceptibility in people of African ancestry. Here, we examined genome-wide and candidate gene associations with systolic blood pressure (SBP) and diastolic blood pressure (DBP) using the Candidate Gene Association Resource (CARe) consortium consisting of 8591 AAs. Genotypes included genome-wide single-nucleotide polymorphism (SNP) data utilizing the Affymetrix 6.0 array with imputation to 2.5 million HapMap SNPs and candidate gene SNP data utilizing a 50K cardiovascular gene-centric array (ITMAT-Broad-CARe [IBC] array). For Affymetrix data, the strongest signal for DBP was rs10474346 (P= 3.6 × 10(-8)) located near GPR98 and ARRDC3. For SBP, the strongest signal was rs2258119 in C21orf91 (P= 4.7 × 10(-8)). The top IBC association for SBP was rs2012318 (P= 6.4 × 10(-6)) near SLC25A42 and for DBP was rs2523586 (P= 1.3 × 10(-6)) near HLA-B. None of the top variants replicated in additional AA (n = 11 882) or European-American (n = 69 899) cohorts. We replicated previously reported European-American blood pressure SNPs in our AA samples (SH2B3, P= 0.009; TBX3-TBX5, P= 0.03; and CSK-ULK3, P= 0.0004). These genetic loci represent the best evidence of genetic influences on SBP and DBP in AAs to date. More broadly, this work supports that notion that blood pressure among AAs is a trait with genetic underpinnings but also with significant complexity.

Polymorphisms in toll-like receptor 4 and toll-like receptor 9 influence viral load in a seroincident cohort of HIV-1-infected individuals.

OBJECTIVES: Toll-like receptors (TLRs) are innate immune sensors that are integral to resisting chronic and opportunistic infections. Mounting evidence implicates TLR polymorphisms in susceptibilities to various infectious diseases, including HIV-1. We investigated the impact of TLR single nucleotide polymorphisms (SNPs) on clinical outcome in a seroincident cohort of HIV-1-infected volunteers. DESIGN: We analyzed TLR SNPs in 201 antiretroviral treatment-naive HIV-1-infected volunteers from a longitudinal seroincident cohort with regular follow-up intervals (median follow-up 4.2 years, interquartile range 4.4). Participants were stratified into two groups according to either disease progression, defined as peripheral blood CD4(+) T-cell decline over time, or peak and setpoint viral load. METHODS: Haplotype tagging SNPs from TLR2, TLR3, TLR4, and TLR9 were detected by mass array genotyping, and CD4(+) T-cell counts and viral load measurements were determined prior to antiretroviral therapy initiation. The association of TLR haplotypes with viral load and rapid progression was assessed by multivariate regression models using age and sex as covariates. RESULTS: Two TLR4 SNPs in strong linkage disequilibrium [1063 A/G (D299G) and 1363 C/T (T399I)] were more frequent among individuals with high peak viral load compared with low/moderate peak viral load (odds ratio 6.65, 95% confidence interval 2.19-20.46, P < 0.001; adjusted P = 0.002 for 1063 A/G). In addition, a TLR9 SNP previously associated with slow progression was found less frequently among individuals with high viral setpoint compared with low/moderate setpoint (odds ratio 0.29, 95% confidence interval 0.13-0.65, P = 0.003, adjusted P = 0.04). CONCLUSION: This study suggests a potentially new role for TLR4 polymorphisms in HIV-1 peak viral load and confirms a role for TLR9 polymorphisms in disease progression.

Alterations of the 5'untranslated region of SLC16A12 lead to age-related cataract.

PURPOSE. Knowledge of genetic factors predisposing to age-related cataract is very limited. The aim of this study was to identify DNA sequences that either lead to or predispose for this disease. METHODS. The candidate gene SLC16A12, which encodes a solute carrier of the monocarboxylate transporter family, was sequenced in 484 patients with cataract (134 with juvenile cataract, 350 with age-related cataract) and 190 control subjects. Expression studies included luciferase reporter assay and RT-PCR experiments. RESULTS. One patient with age-related cataract showed a novel heterozygous mutation (c.-17A>G) in the 5'untranslated region (5'UTR). This mutation is in cis with the minor G-allele of the single nucleotide polymorphism (SNP) rs3740030 (c.-42T/G), also within the 5'UTR. Using a luciferase reporter assay system, a construct with the patient's haplotype caused a significant upregulation of luciferase activity. In comparison, the SNP G-allele alone promoted less activity, but that amount was still significantly higher than the amount of the common T-allele. Analysis of SLC16A12 transcripts in surrogate tissue demonstrated striking allele-specific differences causing 5'UTR heterogeneity with respect to sequence and quantity. These differences in gene expression were mirrored in an allele-specific predisposition to age-related cataract, as determined in a Swiss population (odds ratio approximately 2.2; confidence intervals, 1.23-4.3). CONCLUSIONS. The monocarboxylate transporter SLC16A12 may contribute to age-related cataract. Sequences within the 5'UTR modulate translational efficiency with pathogenic consequences.

FTO genotype is associated with phenotypic variability of body mass index.

There is evidence across several species for genetic control of phenotypic variation of complex traits, such that the variance among phenotypes is genotype dependent. Understanding genetic control of variability is important in evolutionary biology, agricultural selection programmes and human medicine, yet for complex traits, no individual genetic variants associated with variance, as opposed to the mean, have been identified. Here we perform a meta-analysis of genome-wide association studies of phenotypic variation using ∼170,000 samples on height and body mass index (BMI) in human populations. We report evidence that the single nucleotide polymorphism (SNP) rs7202116 at the FTO gene locus, which is known to be associated with obesity (as measured by mean BMI for each rs7202116 genotype), is also associated with phenotypic variability. We show that the results are not due to scale effects or other artefacts, and find no other experiment-wise significant evidence for effects on variability, either at loci other than FTO for BMI or at any locus for height. The difference in variance for BMI among individuals with opposite homozygous genotypes at the FTO locus is approximately 7%, corresponding to a difference of ∼0.5 kilograms in the standard deviation of weight. Our results indicate that genetic variants can be discovered that are associated with variability, and that between-person variability in obesity can partly be explained by the genotype at the FTO locus. The results are consistent with reported FTO by environment interactions for BMI, possibly mediated by DNA methylation. Our BMI results for other SNPs and our height results for all SNPs suggest that most genetic variants, including those that influence mean height or mean BMI, are not associated with phenotypic variance, or that their effects on variability are too small to detect even with samples sizes greater than 100,000.

Role of β-TrCP1, variant and homologue in HIV-1 life cycle

Summary The CD4 molecule plays a key role in AIDS pathogenesis, it is required for entry of the virus into permissive cells and its subsequent down-modulation of the cell surface is a hallmark of HN-1 infected cells. The virus encodes no less than three proteins that participate in this process: Nef, Vpu and Env. Vpu protein interacts with CD4 within the endoplasmic reticulum of infected cells, where it targets CD4 for degradation through the interaction with a cellular protein named ß-TrCP1. This F-box protein functions as the substrate recognition subunit of the SCF ß-Trcr E3 ubiquitin ligase, which normally induce the ubiquitination and subsequent degradation of various proteins such as ß-catenin and IxBa. Mammals possess a homologue of ß-TrCP1, HOS, also named ß-TrCP2 which has a cytoplasmic subcellular distribution. Structural analysis of the ligand-binding domain of both homologues shows striking surface similarities. Both F-box proteins have a redundant role in a number of cellular processes; however the potential role of ß-TrCP2 in HIV-1 infected cells has not been evaluated. In the present study, we assessed the existence of génetic variants of BRTC, encoding ß-TrCP1, and evaluated whether these variants would affect CD4 down-modulation. Additionally, we determined whether ß-TrCP2 shares with its homologue structural and functional properties that would allow it to bind Vpu, modulate CD4 expression, and thus participate in HN-1 pathogenesis. We identified a single nucleotide polymorphism present in the human population with an allelic frequency of 0.03 that leads to the substitution of alanine 507 by a serine. However, we showed by transient transfection in HeLa CD4+ cells that this variant behaves as ß-TrCP1 with respect to CD4 down-modulation. We established transient expression systems in HeLa CD4+ cells to test whether ß-TrCP2 is implicated in Vpu-mediated CD4 down-modulation. We show by coimmunoprecipitation experiments that ß-TrCP2 binds Vpu and is able to induce CD4 down-modulation as efficiently as ß-TrCP1. In two different cell lines, HeLa CD4+ and Jurkat, Vpu-mediated CD4 down-modulation could not be completely reversed through the silencing of endogenous ß-TrCP 1 or ß-TrCP2 individually, but required both genes to be silenced simultaneously. We evaluated the role of ß-TrCP1 and ß-TrCP2 in HIV-1 life cycle using silencing prior to actual viral infection. Both ß-TrCP1 and ß-TrCP2 contributed to CD4 down-modulation during aone-cycle viral infection iri Ghost cells. In addition, the combined silencing of both homologues in the absence of env and nef reversed CD4 down-modulation, showing that ß-TrCP 1 and ß-TrCP2 represent the main and additive effectors of HIV-1 encoded Vpu. In addition, we showed that silencing of ß-TrCPI but not ß-TrCP2 induced a decrease of HIV-1 LTR-driven expression. In a transient transfection system with Tat and a LTR luciferase reporter, both homologues modulated LTR-driven expression. The present study revealed that ß-TrCP2 represents a novel protein participating in HIV-1 cycle and complete comprehension of the complex interplay occurring between the two F-Box will improve our understanding of HIV-1 infection. Résumé La molécule CD4 joue un rôle clef dans la pathogenèse du SIDA ; elle est requise pour l'entrée du virus dans les cellules permissives et la diminution de sa concentration au niveau de la surface cellulaire est une importante caractéristique des cellules infectées par le VIH-1. Le virus encode pas moins de trois protéines qui participent à ce processus Nef, Vpu et Env. La protéine Vpu lie CD4 au niveau du réticulum endoplasmique et induit sa dégradation en interagissant avec une protéine cellulaire nommée ß-TrCP 1. Cette protéine de type F-Box est une sous unité du complexe ubiquitine-ligase E3 SCFß-TrCP. Elle permet la reconnaissance du substrat par le complexe qui induit l'ubiquitination et la subséquente dégradation de diverses protéines cellulaires comme la ß-catenin ou IκBα. Les mammifères possèdent un homologue à ß-TrCP1appelé ß-TrCP2 (ou HOS). L'analyse comparative du domaine permettant la reconnaissance des substrats des deux homologues montre de frappantes similarités. Le rôle de ß-TrCP2 dans le cycle viral du VIH-1 n'a pas encore été évalué. Lors de cette étude, nous avons recherché l'existence de variants génétique de BTRC (codant pour ß-TrCP1) et nous avons évalué si ces variants pourraient affecter la dégradation des molécules CD4 induite par le virus. Nous avons ainsi identifié un polymorphisme présent dans la population humaine avec une fréquence allélique de 0.03 qui consiste en une substitution de l'alanine 507 par une sérine. Nous avons cependant montré par transfection dans des cellules HeLa CD4+ que ce variant se comporte comme ß-TrCP 1 en ce qui concerne la modulation de CD4. De plus, nous avons déterminé si ß-TrCP2 partageait avec son homologue des propriétés structurelles et fonctionnelles qui lui permettraient de lier Vpu, moduler la concentration de CD4 et ainsi prendre part à la pathogenèse du SIDA. Pour ce faire, nous avons établi un système d'expression temporaire dans des cellules HeLa CD4+. Par co-immunoprécipitation, nous avons montré que ß-TrCP2 lie Vpu et est capable d'induire la dégradation de CD4 aussi efficacement que ß-TrCP1. Dans deux différentes lignées cellulaires, HeLa CD4+ et Jurkat, la dégradation de CD4 n'a pu être complètement inhibée par le silencing individuel de ß-TrCP 1 ou ß-TrCP2, mais nécessitait le silencing simultané des 2 gènes. Nous avons évalué le rôle des deux homologues dans le cycle viral du VIH-1 en infectant des cellules Ghost avec le virus après avoir effectué un silencing des deux protéines. Nous avons ainsi montré que ß-TrCP 1 et ß-TrCP2 contribuent de manière additive à la dégradation de CD4 induite par une infection du VIH-1. Le silencing combiné des deux homologues inhiba complètement cette dégradation en l'absence de env et nef, prouvant qu'aucune autre voie ne participe à ce processus: En outre, nous avons montré que le silencing de ß-TrCP 1 mais pas celui de ß-TrCP2 induisait une diminution de l'expression virale sous contrôle du LTR. Nous n'avons cependant pas été en mesure de reconstituer cet effet en exprimant Tat et un gène reporteur sous contrôle du LTR dans des cellules HeLa CD4+. Le présent travail révèle que ß-TrCP2 représente une nouvelle protéine participant dans le cycle viral du VIH-1. Une complète compréhension de l'effet de chacun des deux homologues sur le cycle viral permettra d'améliorer notre compréhension de l'infection par le VIH-1.

Low Y chromosome variation in Saudi-Arabian hamadryas baboons (Papio hamadryas hamadryas).

It is important to characterise the amount of variation on the mammalian Y chromosome in order to assess its potential for use in evolutionary studies. We report very low levels of polymorphism on the Y chromosome of Saudi-Arabian hamadryas baboons, Papio hamadryas hamadryas. We found no segregating sites on the Y, despite sequence analysis of 3 kb noncontiguous intron sequence in 16 males with divergent autosomal microsatellite genotypes, and a further analysis of 1.1 kb intron sequence in 97 males from four populations by SSCP. In addition, we tested seven human-derived Y-linked microsatellites in baboons. Only four of these loci were male-specific and only one was polymorphic in our 97 male sample set. Polymorphism on the Y chromosome of Arabian hamadryas appears to be low compared to other primate species for which data are available (eg humans, chimpanzees and bonobos). Low effective population size (Ne) of paternal genes due to polygyny and female-biased adult sex ratio is a potential reason for low Y chromosome variation in this species. However, low Ne for the Y should be counterbalanced to some extent by the species' atypical pattern of male philopatry and female-biased dispersal. Allelic richness averaged over seven loci was not significantly different between an African and an Arabian population, suggesting that loss of variation during the colonisation of Arabia does not explain low Y variation. Finally, in the absence of nucleotide polymorphism, it is unclear to what extent selection could be responsible for low Y variation in this species.

Genome-wide profiling of blood pressure in adults and children.

Hypertension is an important determinant of cardiovascular morbidity and mortality and has a substantial heritability, which is likely of polygenic origin. The aim of this study was to assess to what extent multiple common genetic variants contribute to blood pressure regulation in both adults and children and to assess overlap in variants between different age groups, using genome-wide profiling. Single nucleotide polymorphism sets were defined based on a meta-analysis of genome-wide association studies on systolic blood pressure and diastolic blood pressure performed by the Cohort for Heart and Aging Research in Genome Epidemiology (n=29 136), using different P value thresholds for selecting single nucleotide polymorphisms. Subsequently, genetic risk scores for systolic blood pressure and diastolic blood pressure were calculated in an independent adult population (n=2072) and a child population (n=1034). The explained variance of the genetic risk scores was evaluated using linear regression models, including sex, age, and body mass index. Genetic risk scores, including also many nongenome-wide significant single nucleotide polymorphisms, explained more of the variance than scores based only on very significant single nucleotide polymorphisms in adults and children. Genetic risk scores significantly explained ≤1.2% (P=9.6*10(-8)) of the variance in adult systolic blood pressure and 0.8% (P=0.004) in children. For diastolic blood pressure, the variance explained was similar in adults and children (1.7% [P=8.9*10(-10)] and 1.4% [P=3.3*10(-5)], respectively). These findings suggest the presence of many genetic loci with small effects on blood pressure regulation both in adults and children, indicating also a (partly) common polygenic regulation of blood pressure throughout different periods of life.

FastEpistasis: a high performance computing solution for quantitative trait epistasis.

Motivation: Genome-wide association studies have become widely used tools to study effects of genetic variants on complex diseases. While it is of great interest to extend existing analysis methods by considering interaction effects between pairs of loci, the large number of possible tests presents a significant computational challenge. The number of computations is further multiplied in the study of gene expression quantitative trait mapping, in which tests are performed for thousands of gene phenotypes simultaneously. Results: We present FastEpistasis, an efficient parallel solution extending the PLINK epistasis module, designed to test for epistasis effects when analyzing continuous phenotypes. Our results show that the algorithm scales with the number of processors and offers a reduction in computation time when several phenotypes are analyzed simultaneously. FastEpistasis is capable of testing the association of a continuous trait with all single nucleotide polymorphism ( SNP) pairs from 500 000 SNPs, totaling 125 billion tests, in a population of 5000 individuals in 29, 4 or 0.5 days using 8, 64 or 512 processors.