gemBS: high throughput processing for DNA methylation data from bisulfite sequencing.

MOTIVATION: DNA methylation is essential for normal embryogenesis and development in mammals and can be captured at single base pair resolution by whole genome bisulfite sequencing (WGBS). Current available analysis tools are becoming rapidly outdated as they lack sensible functionality and efficiency to handle large amounts of data now commonly created. RESULTS: We developed gemBS, a fast high-throughput bioinformatics pipeline specifically designed for large scale BS-Seq analysis that combines a high performance BS-mapper (GEM3) and a variant caller specifically for BS-Seq data (BScall). gemBS provides genotype information and methylation estimates for all genomic cytosines in different contexts (CpG and non-CpG) and a set of quality reports for comprehensive and reproducible analysis. gemBS is highly modular and can be easily automated, while producing robust and accurate results. AVAILABILITY AND IMPLEMENTATION: gemBS is released under the GNU GPLv3+ license. Source code and documentation are freely available from www.statgen.cat/gemBS. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Distinct DNA methylomes of newborns and centenarians.

Human aging cannot be fully understood in terms of the constrained genetic setting. Epigenetic drift is an alternative means of explaining age-associated alterations. To address this issue, we performed whole-genome bisulfite sequencing (WGBS) of newborn and centenarian genomes. The centenarian DNA had a lower DNA methylation content and a reduced correlation in the methylation status of neighboring cytosine--phosphate--guanine (CpGs) throughout the genome in comparison with the more homogeneously methylated newborn DNA. The more hypomethylated CpGs observed in the centenarian DNA compared with the neonate covered all genomic compartments, such as promoters, exonic, intronic, and intergenic regions. For regulatory regions, the most hypomethylated sequences in the centenarian DNA were present mainly at CpG-poor promoters and in tissue-specific genes, whereas a greater level of DNA methylation was observed in CpG island promoters. We extended the study to a larger cohort of newborn and nonagenarian samples using a 450,000 CpG-site DNA methylation microarray that reinforced the observation of more hypomethylated DNA sequences in the advanced age group. WGBS and 450,000 analyses of middle-age individuals demonstrated DNA methylomes in the crossroad between the newborn and the nonagenarian/centenarian groups. Our study constitutes a unique DNA methylation analysis of the extreme points of human life at a single-nucleotide resolution level.

Genetic variants associated with Lp(a) lipoprotein level and coronary disease.

BACKGROUND: An increased level of Lp(a) lipoprotein has been identified as a risk factor for coronary artery disease that is highly heritable. The genetic determinants of the Lp(a) lipoprotein level and their relevance for the risk of coronary disease are incompletely understood. METHODS: We used a novel gene chip containing 48,742 single-nucleotide polymorphisms (SNPs) in 2100 candidate genes to test for associations in 3145 case subjects with coronary disease and 3352 control subjects. Replication was tested in three independent populations involving 4846 additional case subjects with coronary disease and 4594 control subjects. RESULTS: Three chromosomal regions (6q26-27, 9p21, and 1p13) were strongly associated with the risk of coronary disease. The LPA locus on 6q26-27 encoding Lp(a) lipoprotein had the strongest association. We identified a common variant (rs10455872) at the LPA locus with an odds ratio for coronary disease of 1.70 (95% confidence interval [CI], 1.49 to 1.95) and another independent variant (rs3798220) with an odds ratio of 1.92 (95% CI, 1.48 to 2.49). Both variants were strongly associated with an increased level of Lp(a) lipoprotein, a reduced copy number in LPA (which determines the number of kringle IV-type 2 repeats), and a small Lp(a) lipoprotein size. Replication studies confirmed the effects of both variants on the Lp(a) lipoprotein level and the risk of coronary disease. A meta-analysis showed that with a genotype score involving both LPA SNPs, the odds ratios for coronary disease were 1.51 (95% CI, 1.38 to 1.66) for one variant and 2.57 (95% CI, 1.80 to 3.67) for two or more variants. After adjustment for the Lp(a) lipoprotein level, the association between the LPA genotype score and the risk of coronary disease was abolished. CONCLUSIONS: We identified two LPA variants that were strongly associated with both an increased level of Lp(a) lipoprotein and an increased risk of coronary disease. Our findings provide support for a causal role of Lp(a) lipoprotein in coronary disease.

Impact of DNA methylation on 3D genome structure.

Determining the effect of DNA methylation on chromatin structure and function in higher organisms is challenging due to the extreme complexity of epigenetic regulation. We studied a simpler model system, budding yeast, that lacks DNA methylation machinery making it a perfect model system to study the intrinsic role of DNA methylation in chromatin structure and function. We expressed the murine DNA methyltransferases in Saccharomyces cerevisiae and analyzed the correlation between DNA methylation, nucleosome positioning, gene expression and 3D genome organization. Despite lacking the machinery for positioning and reading methylation marks, induced DNA methylation follows a conserved pattern with low methylation levels at the 5' end of the gene increasing gradually toward the 3' end, with concentration of methylated DNA in linkers and nucleosome free regions, and with actively expressed genes showing low and high levels of methylation at transcription start and terminating sites respectively, mimicking the patterns seen in mammals. We also see that DNA methylation increases chromatin condensation in peri-centromeric regions, decreases overall DNA flexibility, and favors the heterochromatin state. Taken together, these results demonstrate that methylation intrinsically modulates chromatin structure and function even in the absence of cellular machinery evolved to recognize and process the methylation signal.

Investigation of the fine structure of European populations with applications to disease association studies.

An investigation into fine-scale European population structure was carried out using high-density genetic variation on nearly 6000 individuals originating from across Europe. The individuals were collected as control samples and were genotyped with more than 300 000 SNPs in genome-wide association studies using the Illumina Infinium platform. A major East-West gradient from Russian (Moscow) samples to Spanish samples was identified as the first principal component (PC) of the genetic diversity. The second PC identified a North-South gradient from Norway and Sweden to Romania and Spain. Variation of frequencies at markers in three separate genomic regions, surrounding LCT, HLA and HERC2, were strongly associated with this gradient. The next 18 PCs also accounted for a significant proportion of genetic diversity observed in the sample. We present a method to predict the ethnic origin of samples by comparing the sample genotypes with those from a reference set of samples of known origin. These predictions can be performed using just summary information on the known samples, and individual genotype data are not required. We discuss issues raised by these data and analyses for association studies including the matching of case-only cohorts to appropriate pre-collected control samples for genome-wide association studies.

Information recovery from low coverage whole-genome bisulfite sequencing.

The cost of whole-genome bisulfite sequencing (WGBS) remains a bottleneck for many studies and it is therefore imperative to extract as much information as possible from a given dataset. This is particularly important because even at the recommend 30X coverage for reference methylomes, up to 50% of high-resolution features such as differentially methylated positions (DMPs) cannot be called with current methods as determined by saturation analysis. To address this limitation, we have developed a tool that dynamically segments WGBS methylomes into blocks of comethylation (COMETs) from which lost information can be recovered in the form of differentially methylated COMETs (DMCs). Using this tool, we demonstrate recovery of ∼30% of the lost DMP information content as DMCs even at very low (5X) coverage. This constitutes twice the amount that can be recovered using an existing method based on differentially methylated regions (DMRs). In addition, we explored the relationship between COMETs and haplotypes in lymphoblastoid cell lines of African and European origin. Using best fit analysis, we show COMETs to be correlated in a population-specific manner, suggesting that this type of dynamic segmentation may be useful for integrated (epi)genome-wide association studies in the future.

Single-nucleotide polymorphism versus microsatellite markers in a combined linkage and segregation analysis of a quantitative trait.

Increasingly, single-nucleotide polymorphism (SNP) markers are being used in preference to microsatellite markers. However, methods developed for microsatellites may be problematic when applied to SNP markers. We evaluated the results of using SNPs vs. microsatellites in Monte Carlo Markov chain (MCMC) oligogenic combined segregation and linkage analysis methods. These methods were developed with microsatellite markers in mind. We selected chromosome 7 from the Collaborative Study on the Genetics of Alcoholism dataset for analysis because linkage to an electrophysiological trait had been reported there. We found linkage in the same region of chromosome 7 with the Affymetrix SNP data, the Illumina SNP data, and the microsatellite marker data. The MCMC sampler appears to mix with both types of data. The sampler implemented in this MCMC oligogenic combined segregation and linkage analysis appears to handle SNP data as well as microsatellite data and it is possible that the localizations with the SNP data are better.

A comprehensive assessment of somatic mutation detection in cancer using whole-genome sequencing.

As whole-genome sequencing for cancer genome analysis becomes a clinical tool, a full understanding of the variables affecting sequencing analysis output is required. Here using tumour-normal sample pairs from two different types of cancer, chronic lymphocytic leukaemia and medulloblastoma, we conduct a benchmarking exercise within the context of the International Cancer Genome Consortium. We compare sequencing methods, analysis pipelines and validation methods. We show that using PCR-free methods and increasing sequencing depth to ∼ 100 × shows benefits, as long as the tumour:control coverage ratio remains balanced. We observe widely varying mutation call rates and low concordance among analysis pipelines, reflecting the artefact-prone nature of the raw data and lack of standards for dealing with the artefacts. However, we show that, using the benchmark mutation set we have created, many issues are in fact easy to remedy and have an immediate positive impact on mutation detection accuracy.

Candidate genes involved in cardiovascular risk factors by a family-based association study on the island of Kosrae, Federated States of Micronesia.

Altered plasma levels of lipids and lipoproteins, obesity, hypertension, and diabetes are major risk factors for atherosclerotic cardiovascular disease. To identify genes that affect these traits and disorders, we looked for association between markers in candidate genes (apolipoprotein AII (apo AII), apolipoprotein AI-CIII-AIV gene cluster (apo AI-CIII-AIV), apolipoprotein E (apo E), cholesteryl ester transfer protein (CETP), cholesterol 7alpha-hydroxylase (CYP7a), hepatic lipase (HL), and microsomal triglyceride transfer protein (MTP)) and known risk factors (triglycerides (Tg), total cholesterol (TC), apolipoprotein AI (apo AI), apolipoprotein AII (apo AII), apolipoprotein B (apo B), body mass index (BMI), blood pressure (BP), leptin, and fasting blood sugar (FBS) levels.) A total of 1,102 individuals from the Pacific island of Kosrae were genotyped for the following markers: Apo AII/MspI, Apo CIII/SstI, Apo AI/XmnI, Apo E/HhaI, CETP/TaqIB, CYP7a/BsaI, HL/DraI, and MTP/HhpI. After testing for population stratification, family-based association analysis was carried out. Novel associations found were: 1) the apo AII/MspI with apo AI and BP levels, 2) the CYP7a/BsaI with apo AI and BMI levels. We also confirmed the following associations: 1) the apo AII/MspI with Tg level; 2) the apo CIII/SstI with Tg, TC, and apo B levels; 3) the Apo E/HhaI E2, E3, and E4 alleles with TC, apo AI, and apo B levels; and 4) the CETP/TaqIB with apo AI level. We further confirmed the connection between the apo AII gene and Tg level by a nonparametric linkage analysis. We therefore conclude that many of these candidate genes may play a significant role in susceptibility to heart disease.

Tuning of natural killer cell reactivity by NKp46 and Helios calibrates T cell responses.

Natural killer (NK) cells are lymphocytes involved in antimicrobial and antitumoral immune responses. Using N-ethyl-N-nitrosourea mutagenesis in mice, we identified a mutant with increased resistance to viral infections because of the presence of hyperresponsive NK cells. Whole-genome sequencing and functional analysis revealed a loss-of-function mutation in the Ncr1 gene encoding the activating receptor NKp46. The down-regulation of NK cell activity by NKp46 was associated with the silencing of the Helios transcription factor in NK cells. NKp46 was critical for the subsequent development of antiviral and antibacterial T cell responses, which suggests that the regulation of NK cell function by NKp46 allows for the optimal development of adaptive immune responses. NKp46 blockade enhanced NK cell reactivity in vivo, which could enable the design of immunostimulation strategies in humans.

Genome-wide association study identifies eight loci associated with blood pressure.

Elevated blood pressure is a common, heritable cause of cardiovascular disease worldwide. To date, identification of common genetic variants influencing blood pressure has proven challenging. We tested 2.5 million genotyped and imputed SNPs for association with systolic and diastolic blood pressure in 34,433 subjects of European ancestry from the Global BPgen consortium and followed up findings with direct genotyping (N ≤ 71,225 European ancestry, N ≤ 12,889 Indian Asian ancestry) and in silico comparison (CHARGE consortium, N = 29,136). We identified association between systolic or diastolic blood pressure and common variants in eight regions near the CYP17A1 (P = 7 × 10(-24)), CYP1A2 (P = 1 × 10(-23)), FGF5 (P = 1 × 10(-21)), SH2B3 (P = 3 × 10(-18)), MTHFR (P = 2 × 10(-13)), c10orf107 (P = 1 × 10(-9)), ZNF652 (P = 5 × 10(-9)) and PLCD3 (P = 1 × 10(-8)) genes. All variants associated with continuous blood pressure were associated with dichotomous hypertension. These associations between common variants and blood pressure and hypertension offer mechanistic insights into the regulation of blood pressure and may point to novel targets for interventions to prevent cardiovascular disease.

Genetic factors of susceptibility and of severity in primary biliary cirrhosis.

BACKGROUND/AIMS: In primary biliary cirrhosis (PBC), pathogenesis is influenced by genetic factors that remain poorly elucidated up to now. We investigated the impact of sequence diversity in candidate genes involved in immunity (CTLA-4 and TNFalpha), in bile formation (10 hepatobiliary transporter genes) and in the adaptative response to cholestasis (three nuclear receptor genes) on the susceptibility and severity of PBC. METHODS: A total of 42 Ht SNPs were identified and compared in 258 PBC patients and two independent groups of 286 and 269 healthy controls. All participants were white continental individuals with French ancestry. RESULTS: Ht SNPs of CTLA-4 and TNFalpha genes were significantly associated with susceptibility to PBC. The progression rate of liver disease under ursodeoxycholic acid (UDCA) therapy was significantly linked to SNPs of TNFalpha and SLC4A2/anion exchanger 2 (AE2) genes. A multivariate Cox regression analysis including clinical and biochemical parameters showed that SLC4A2/AE2 variant was an independent prognostic factor. CONCLUSIONS: These data point to a primary role of genes encoding regulators of the immune system in the susceptibility to PBC. They also demonstrate that allelic variations in TNFalpha and SLC4A2/AE2 have a significant impact on the evolutive profile of PBC under UDCA therapy.

Newly identified loci that influence lipid concentrations and risk of coronary artery disease.

To identify genetic variants influencing plasma lipid concentrations, we first used genotype imputation and meta-analysis to combine three genome-wide scans totaling 8,816 individuals and comprising 6,068 individuals specific to our study (1,874 individuals from the FUSION study of type 2 diabetes and 4,184 individuals from the SardiNIA study of aging-associated variables) and 2,758 individuals from the Diabetes Genetics Initiative, reported in a companion study in this issue. We subsequently examined promising signals in 11,569 additional individuals. Overall, we identify strongly associated variants in eleven loci previously implicated in lipid metabolism (ABCA1, the APOA5-APOA4-APOC3-APOA1 and APOE-APOC clusters, APOB, CETP, GCKR, LDLR, LPL, LIPC, LIPG and PCSK9) and also in several newly identified loci (near MVK-MMAB and GALNT2, with variants primarily associated with high-density lipoprotein (HDL) cholesterol; near SORT1, with variants primarily associated with low-density lipoprotein (LDL) cholesterol; near TRIB1, MLXIPL and ANGPTL3, with variants primarily associated with triglycerides; and a locus encompassing several genes near NCAN, with variants strongly associated with both triglycerides and LDL cholesterol). Notably, the 11 independent variants associated with increased LDL cholesterol concentrations in our study also showed increased frequency in a sample of coronary artery disease cases versus controls.

Associations of the IL2Ralpha, IL4Ralpha, IL10Ralpha, and IFN (gamma) R1 cytokine receptor genes with AIDS progression in a French AIDS cohort.

We have performed an extensive analysis of Th1/Th2 cytokine receptors IL2Ralpha, IL4Ralpha, IL10Ralpha, and IFNgammaR1 gene polymorphisms to evaluate their impact on AIDS progression. The coding regions and promoters of these genes were sequenced in the genetics of resistance to immunodeficiency virus cohort, composed of 327 HIV-1-positive patients with extreme progression phenotypes, slow and rapid progressors, and of 446 healthy control subjects, all of them of Caucasian descent. Overall, 104 single nucleotide polymorphisms and four insertions/deletions with a minor allelic frequency higher than 1% were identified, 21 of them being newly characterized. We observed weak associations for 13 polymorphisms of IL2Ralpha, IL4Ralpha, IL10Ralpha, and IFNgammaR1, and 11 haplotypes of IL2Ralpha, IL4Ralpha, and IFNgammaR1. However, we could not relate these positive signals to any relevant biological information on the gene function. To affirm these putative associations in AIDS, further confirmation on other AIDS cohorts will be needed. This complete catalog of polymorphisms in IL2Ralpha, IL4Ralpha, IL10Ralpha, and IFNgammaR1 cytokine receptor genes should also be useful for investigating associations in other immune-related diseases.

Linkage analysis of quantitative traits for obesity, diabetes, hypertension, and dyslipidemia on the island of Kosrae, Federated States of Micronesia.

Obesity, diabetes, hypertension, and heart disease are highly heritable conditions that in aggregate are the major causes of morbidity and mortality in the developed world and are growing problems in developing countries. To map the causal genes, we conducted a population screen for these conditions on the Pacific Island of Kosrae. Family history and genetic data were used to construct a pedigree for the island. Analysis of the pedigree showed highly significant heritability for the metabolic traits under study. DNA samples from 2,188 participants were genotyped with 405 microsatellite markers with an average intermarker distance of 11 cM. A protocol using loki, a Markov chain Monte Carlo sampling method, was developed to analyze the Kosraen pedigree for height, a model quantitative trait. Robust quantitative trait loci for height were found on 10q21 and 1p31. This protocol was used to map a set of metabolic traits, including plasma leptin to chromosome region 5q35; systolic blood pressure to 20p12; total cholesterol to 19p13, 12q24, and 16qter; hip circumference to 10q25 and 4q23; body mass index to 18p11 and 20q13; apolipoprotein B to 2p24-25; weight to 18q21; and fasting blood sugar to 1q31-1q43. Several of these same chromosomal regions have been identified in previous studies validating the use of loki. These studies add information about the genetics of the metabolic syndrome and establish an analytical approach for linkage analysis of complex pedigrees. These results also lay the foundation for whole genome scans with dense sets of SNPs aimed to identifying causal genes.

Phenotypical characterisation of the isolated norfolk island population focusing on epidemiological indicators of cardiovascular disease.

OBJECTIVES: Only 193 people from Pitcairn Island, all descended from 9 'Bounty' mutineers and 12 Tahitian women, moved to the uninhabited Norfolk Island in 1856. Our objective was to assess the population of Norfolk Island, several thousand km off the eastern coast of Australia, as a genetic isolate of potential use for cardiovascular disease (CVD) gene mapping. METHODS: A total of 602 participants, approximately two thirds of the island's present adult population, were characterized for a panel of CVD risk factors. Statistical power and heritability were calculated. RESULTS: Norfolk Islander's possess an increased prevalence of hypertension, obesity and multiple CVD risk factors when compared to outbred Caucasian populations. 64% of the study participants were descendents of the island's original founder population. Triglycerides, cholesterol, and blood pressures all had heritabilities above 0.2. CONCLUSIONS: The Norfolk Island population is a potentially useful genetic isolate for gene mapping studies aimed at identifying CVD risk factor quantitative trait loci (QTL).

Exhaustive genotyping of the CEM15 (APOBEC3G) gene and absence of association with AIDS progression in a French cohort.

CEM15 (or APOBEC3G) has recently been identified as an inhibitor of human immunodeficiency virus type 1 (HIV-1) replication in vitro. To evaluate the impact of its genetic variations on the progression of acquired immunodeficiency syndrome (AIDS), we have performed an extensive genetic analysis of CEM15. We have sequenced CEM15 in a cohort of 327 HIV-1-seropositive patients with extreme disease progression phenotypes--either slow progression or rapid progression--and in 446 healthy control subjects, all of white descent. We have identified 29 polymorphisms with allele frequencies >1%, 14 of which were newly characterized. There were no significant associations between the polymorphisms or haplotypes of CEM15 and a disease progression phenotype in our cohort.

Expression profiling of lineage differentiation in pluripotential human embryonal carcinoma cells.

Pluripotential human embryonal carcinoma (EC) cell linesundergo differentiation programs resembling those occurring in embryonal stem cells during development. Expression profiling was performed during the terminal differentiation of the EC cell line, NTera2/Clone D1 by all-trans-retinoic acid. Time-response analysis via clustering of >12,000 human transcripts revealed distinct stages in the transition from an EC cell to neuronal progenitor cells expressing patterning markers compatible with posterior hindbrain fates followed by the appearance of immature postmitotic neurons with an evolving synaptic apparatus. Global analysis of gene expression allows monitoring cell fate and differentiation of EC cells in vitro and may provide insight into human embryonal stem cell development.

Genetic variation at the 22q11 PRODH2/DGCR6 locus presents an unusual pattern and increases susceptibility to schizophrenia.

The location of a schizophrenia susceptibility locus at chromosome 22q11 has been suggested by genome-wide linkage studies. Additional support was provided by the observation of a higher-than-expected frequency of 22q11 microdeletions in patients with schizophrenia and the demonstration that approximately 20-30% of individuals with 22q11 microdeletions develop schizophrenia or schizoaffective disorder in adolescence and adulthood. Analysis of the extent of these microdeletions by using polymorphic markers afforded further refinement of this locus to a region of approximately 1.5 Mb. Recently, a high rate of 22q11 microdeletions was also reported for a cohort of 47 patients with Childhood Onset Schizophrenia, a rare and severe form of schizophrenia with onset by age 13. It is therefore likely that this 1.5-Mb region contains one or more genes that predispose to schizophrenia. In three independent samples, we provide evidence for a contribution of the PRODH2/DGCR6 locus in 22q11-associated schizophrenia. We also uncover an unusual pattern of PRODH2 gene variation that mimics the sequence of a linked pseudogene. Several of the pseudogene-like variants we identified result in missense changes at conserved residues and may prevent synthesis of a fully functional enzyme. Our results have implications for understanding the genetic basis of the 22q11-associated psychiatric phenotypes and provide further insights into the genomic instability of this region.