Epigenome-Wide Association Study for Parkinson's Disease

A methylation-based EWAS on carefully phenotyped individuals with Parkinson’s disease (PD) was conducted to reveal prioritised genes and pathways with statistically significant and sizable changes in PD and in the anxiety that often accompanies it. This was followed by subsequent replication of top-ranked CpG sites. Using the Infinium® HumanMethylation 450K beadchip (Illumina Inc., USA), twenty unique genes with a sizable difference in methylation (P adjusted < 0.05, Δβ ≥ 0.2), after correction for multiple testing, were identified between PD and controls, while seventeen were identified between PD with anxiety and PD without anxiety. Twelve top ranked, significantly associated loci in PD were evaluated in an independent replicate population using Sequenom EpiTYPER for 219 individuals with similar phenotypes to the cross-sectional case–control discovery design. FANCC cg14115740 and TNKS2 cg11963436 show significant differential methylation between PD cases and controls using both techniques and their Δβ values, which have the same direction of effect, are reasonable to warrant further investigation

S21 Culture Independent Identification Of Bacterial Communities In The Respiratory Tract Of Patients With Copd, Healthy Non-smokers And Healthy Smokers

Introduction and aims: The role bacteria play in the development and progression of Chronic Obstructive Pulmonary Disease (COPD) is unclear. We used culture-independent methods to describe differences and/or similarities in microbial communities in the lower airways of patients with COPD, healthy non-smokers and smokers.

Methods: Bronchial wash samples were collected from patients with COPD (GOLD 1–3; n = 18), healthy non-smokers (HV; n = 11) and healthy smokers (HS; n = 8). Samples were processed using the Illumina MiSeq platform. The Shannon-Wiener Index (SW) of diversity, lung obstruction (FEV1/FVC ratio) and ordination by Non-Metric Multidimensional Scaling (NMDS) on Bray-Curtis dissimilarity indices were analysed to evaluate how samples were related. Principal component analysis (PCA) was performed to assess the effect specific taxa had within each cohort. Characteristics of each cohort are shown in Table 1.

Results: There was no difference in taxa richness between cohorts (range: 69–71; p = 0.954). Diversity (SW Index) was significantly lower in COPD samples compared to samples from HV and HS (p = 0.009 and p = 0.033, respectively). There was no significant difference between HV and HS (p = 0.186). The FEV1/FVC ratio was significantly lower for COPD compared to HV (p = 9*10–8) and HS (p = 2*10–6), respectively. NMDS analysis showed that communities belonging to either of the healthy groups were more similar to each other than they were to samples belonging to the COPD group. PCA analysis showed that members of Streptococcus sp. and Haemophilus sp. had the largest effect on the variance explained in COPD. In HS, Haemophilus sp., Fusobaterium sp., Actinomyces sp., Prevotella sp. and Veillonella sp. had the largest effect on the variance explained, while in HV Neisseria sp., Porphyromonas sp., Actinomyces sp., Atopobium sp., Prevotella and Veillonella sp. had the largest effect on the variance explained.

Conclusions: The study demonstrates that microbial communities in the lower airways of patients with COPD are significantly different from that seen in healthy comparison groups. Patients with COPD had lower microbial diversity than either of the healthy comparison groups, higher relative abundance of members of Streptococcus sp. and lower relative abundance of a number of key anaerobes.Characteristics

DNA methylation subgroups and the CpG island methylator phenotype in gastric cancer: a comprehensive profiling approach

BACKGROUND: Methylation-induced silencing of promoter CpG islands in tumor suppressor genes plays an important role in human carcinogenesis. In colorectal cancer, the CpG island methylator phenotype (CIMP) is defined as widespread and elevated levels of DNA methylation and CIMP+ tumors have distinctive clinicopathological and molecular features. In contrast, the existence of a comparable CIMP subtype in gastric cancer (GC) has not been clearly established. To further investigate this issue, in the present study we performed comprehensive DNA methylation profiling of a well-characterised series of primary GC.

METHODS: The methylation status of 1,421 autosomal CpG sites located within 768 cancer-related genes was investigated using the Illumina GoldenGate Methylation Panel I assay on DNA extracted from 60 gastric tumors and matched tumor-adjacent gastric tissue pairs. Methylation data was analysed using a recursively partitioned mixture model and investigated for associations with clinicopathological and molecular features including age, Helicobacter pylori status, tumor site, patient survival, microsatellite instability and BRAF and KRAS mutations.

RESULTS: A total of 147 genes were differentially methylated between tumor and matched tumor-adjacent gastric tissue, with HOXA5 and hedgehog signalling being the top-ranked gene and signalling pathway, respectively. Unsupervised clustering of methylation data revealed the existence of 6 subgroups under two main clusters, referred to as L (low methylation; 28% of cases) and H (high methylation; 72%). Female patients were over-represented in the H tumor group compared to L group (36% vs 6%; P = 0.024), however no other significant differences in clinicopathological or molecular features were apparent. CpG sites that were hypermethylated in group H were more frequently located in CpG islands and marked for polycomb occupancy.

CONCLUSIONS: High-throughput methylation analysis implicates genes involved in embryonic development and hedgehog signaling in gastric tumorigenesis. GC is comprised of two major methylation subtypes, with the highly methylated group showing some features consistent with a CpG island methylator phenotype.

Distinct methylation patterns in genes that affect mitochondrial function are associated with kidney disease in blood-derived DNA from individuals with Type 1 diabetes

AIMS: Epigenetic modifications, such as DNA methylation, can influence the risk of developing kidney disease. We studied methylation profiles in genes related to mitochondrial function to assess whether differences in these epigenetic features were associated with diabetic kidney disease in people with Type 1 diabetes.

METHODS: A case-control association study was undertaken (n = 196 individuals with diabetic kidney disease vs. n = 246 individuals without renal disease). Participants were White and diagnosed with Type 1 diabetes before 31 years of age. Genes that encode mitochondrial proteins (n = 780) were downloaded from mitoproteome. org. DNA methylation profiles from blood-derived DNA were generated using the Illumina Infinium HumanMethylation450 (262 samples) and Illumina Infinium HumanMethylation27 (192 samples) arrays. Beta values (β) were calculated and quality control was conducted, including evaluating blind duplicate DNA samples.

RESULTS: Fifty-four Cytosine-phosphate-Guanine probes across 51 unique genes were significantly associated (P ≤ 10(-8) ) with diabetic kidney disease across both the 450K and the 27K methylation arrays. A subanalysis, employing the 450K array, identified 755 Cytosine-phosphate-Guanine probes in 374 genes that were significantly associated (P ≤ 10(-8) ) with end-stage renal disease. Forty-six of the top-ranked variants for diabetic kidney disease were also identified as being differentially methylated in individuals with end-stage renal disease. The largest change in methylation (Δβ = 0.2) was observed for cg03169527 in the TAMM41 gene, chromosome 3p25.2. Three genes, PMPCB, TSFM and AUH, were observed with differential methylation at multiple Cytosine-phosphate-Guanine sites each (P < 10(-12) ).

CONCLUSIONS: Differential methylation in genes that influence mitochondrial function are associated with kidney disease in individuals with Type 1 diabetes. 

Exposure to aflatoxin B1 in utero is associated with DNA methylation in white blood cells of infants in The Gambia

BACKGROUND: Exposure to environmental toxins during embryonic development may lead to epigenetic changes that influence disease risk in later life. Aflatoxin is a contaminant of staple foods in sub-Saharan Africa, is a known human liver carcinogen and has been associated with stunting in infants.

METHODS: We have measured aflatoxin exposure in 115 pregnant women in The Gambia and examined the DNA methylation status of white blood cells from their infants at 2-8 months old (mean 3.6 ± 0.9). Aflatoxin exposure in women was assessed using an ELISA method to measure aflatoxin albumin (AF-alb) adducts in plasma taken at 1-16 weeks of pregnancy. Genome-wide DNA methylation of infant white blood cells was measured using the Illumina Infinium HumanMethylation450beadchip.

RESULTS: AF-alb levels ranged from 3.9 to 458.4 pg/mg albumin. We found that aflatoxin exposure in the mothers was associated to DNA methylation in their infants for 71 CpG sites (false discovery rate < 0.05), with an average effect size of 1.7% change in methylation. Aflatoxin-associated differential methylation was observed in growth factor genes such as FGF12 and IGF1, and immune-related genes such as CCL28, TLR2 and TGFBI. Moreover, one aflatoxin-associated methylation region (corresponding to the miR-4520b locus) was identified.

CONCLUSIONS: This study shows that maternal exposure to aflatoxin during the early stages of pregnancy is associated with differential DNA methylation patterns of infants, including in genes related to growth and immune function. This reinforces the need for interventions to reduce aflatoxin exposure, especially during critical periods of fetal and infant development.

Urothelial cancer gene regulatory networks inferred from large-scale RNAseq, Bead and Oligo gene expression data

BACKGROUND: Urothelial pathogenesis is a complex process driven by an underlying network of interconnected genes. The identification of novel genomic target regions and gene targets that drive urothelial carcinogenesis is crucial in order to improve our current limited understanding of urothelial cancer (UC) on the molecular level. The inference of genome-wide gene regulatory networks (GRN) from large-scale gene expression data provides a promising approach for a detailed investigation of the underlying network structure associated to urothelial carcinogenesis.

METHODS: In our study we inferred and compared three GRNs by the application of the BC3Net inference algorithm to large-scale transitional cell carcinoma gene expression data sets from Illumina RNAseq (179 samples), Illumina Bead arrays (165 samples) and Affymetrix Oligo microarrays (188 samples). We investigated the structural and functional properties of GRNs for the identification of molecular targets associated to urothelial cancer.

RESULTS: We found that the urothelial cancer (UC) GRNs show a significant enrichment of subnetworks that are associated with known cancer hallmarks including cell cycle, immune response, signaling, differentiation and translation. Interestingly, the most prominent subnetworks of co-located genes were found on chromosome regions 5q31.3 (RNAseq), 8q24.3 (Oligo) and 1q23.3 (Bead), which all represent known genomic regions frequently deregulated or aberated in urothelial cancer and other cancer types. Furthermore, the identified hub genes of the individual GRNs, e.g., HID1/DMC1 (tumor development), RNF17/TDRD4 (cancer antigen) and CYP4A11 (angiogenesis/ metastasis) are known cancer associated markers. The GRNs were highly dataset specific on the interaction level between individual genes, but showed large similarities on the biological function level represented by subnetworks. Remarkably, the RNAseq UC GRN showed twice the proportion of significant functional subnetworks. Based on our analysis of inferential and experimental networks the Bead UC GRN showed the lowest performance compared to the RNAseq and Oligo UC GRNs.

CONCLUSION: To our knowledge, this is the first study investigating genome-scale UC GRNs. RNAseq based gene expression data is the data platform of choice for a GRN inference. Our study offers new avenues for the identification of novel putative diagnostic targets for subsequent studies in bladder tumors.

DNA microarray profiling of genes differentially regulated by the histone deacetylase inhibitors vorinostat and LBH589 in colon cancer cell lines

BACKGROUND: Despite the significant progress made in colon cancer chemotherapy, advanced disease remains largely incurable and novel efficacious chemotherapies are urgently needed. Histone deacetylase inhibitors (HDACi) represent a novel class of agents which have demonstrated promising preclinical activity and are undergoing clinical evaluation in colon cancer. The goal of this study was to identify genes in colon cancer cells that are differentially regulated by two clinically advanced hydroxamic acid HDACi, vorinostat and LBH589 to provide rationale for novel drug combination partners and identify a core set of HDACi-regulated genes.

METHODS: HCT116 and HT29 colon cancer cells were treated with LBH589 or vorinostat and growth inhibition, acetylation status and apoptosis were analyzed in response to treatment using MTS, Western blotting and flow cytometric analyses. In addition, gene expression was analyzed using the Illumina Human-6 V2 BeadChip array and Ingenuity Pathway Analysis.

RESULTS: Treatment with either vorinostat or LBH589 rapidly induced histone acetylation, cell cycle arrest and inhibited the growth of both HCT116 and HT29 cells. Bioinformatic analysis of the microarray profiling revealed significant similarity in the genes altered in expression following treatment with the two HDACi tested within each cell line. However, analysis of genes that were altered in expression in the HCT116 and HT29 cells revealed cell-line-specific responses to HDACi treatment. In addition a core cassette of 11 genes modulated by both vorinostat and LBH589 were identified in both colon cancer cell lines analyzed.

CONCLUSION: This study identified HDACi-induced alterations in critical genes involved in nucleotide metabolism, angiogenesis, mitosis and cell survival which may represent potential intervention points for novel therapeutic combinations in colon cancer. This information will assist in the identification of novel pathways and targets that are modulated by HDACi, providing much-needed information on HDACi mechanism of action and providing rationale for novel drug combination partners. We identified a core signature of 11 genes which were modulated by both vorinostat and LBH589 in a similar manner in both cell lines. These core genes will assist in the development and validation of a common gene set which may represent a molecular signature of HDAC inhibition in colon cancer.

SNP discovery using Next Generation Transcriptomic Sequencing in Atlantic herring (Clupea harengus)

The introduction of Next Generation Sequencing (NGS) has revolutionised population genetics, providing studies of non-model species with unprecedented genomic coverage, allowing evolutionary biologists to address questions previously far beyond the reach of available resources. Furthermore, the simple mutation model of Single Nucleotide Polymorphisms (SNPs) permits cost-effective high-throughput genotyping in thousands of individuals simultaneously. Genomic resources are scarce for the Atlantic herring (Clupea harengus), a small pelagic species that sustains high revenue fisheries. This paper details the development of 578 SNPs using a combined NGS and high-throughput genotyping approach. Eight individuals covering the species distribution in the eastern Atlantic were bar-coded and multiplexed into a single cDNA library and sequenced using the 454 GS FLX platform. SNP discovery was performed by de novo sequence clustering and contig assembly, followed by the mapping of reads against consensus contig sequences. Selection of candidate SNPs for genotyping was conducted using an in silico approach. SNP validation and genotyping were performed simultaneously using an Illumina 1,536 GoldenGate assay. Although the conversion rate of candidate SNPs in the genotyping assay cannot be predicted in advance, this approach has the potential to maximise cost and time efficiencies by avoiding expensive and time-consuming laboratory stages of SNP validation. Additionally, the in silico approach leads to lower ascertainment bias in the resulting SNP panel as marker selection is based only on the ability to design primers and the predicted presence of intron-exon boundaries. Consequently SNPs with a wider spectrum of minor allele frequencies (MAFs) will be genotyped in the final panel. The genomic resources presented here represent a valuable multi-purpose resource for developing informative marker panels for population discrimination, microarray development and for population genomic studies in the wild.

Novel tools for conservation genomics: comparing two high-throughput approaches for SNP discovery in the transcriptome of the European hake

The growing accessibility to genomic resources using next-generation sequencing (NGS) technologies has revolutionized the application of molecular genetic tools to ecology and evolutionary studies in non-model organisms. Here we present the case study of the European hake (Merluccius merluccius), one of the most important demersal resources of European fisheries. Two sequencing platforms, the Roche 454 FLX (454) and the Illumina Genome Analyzer (GAII), were used for Single Nucleotide Polymorphisms (SNPs) discovery in the hake muscle transcriptome. De novo transcriptome assembly into unique contigs, annotation, and in silico SNP detection were carried out in parallel for 454 and GAII sequence data. High-throughput genotyping using the Illumina GoldenGate assay was performed for validating 1,536 putative SNPs. Validation results were analysed to compare the performances of 454 and GAII methods and to evaluate the role of several variables (e.g. sequencing depth, intron-exon structure, sequence quality and annotation). Despite well-known differences in sequence length and throughput, the two approaches showed similar assay conversion rates (approximately 43%) and percentages of polymorphic loci (67.5% and 63.3% for GAII and 454, respectively). Both NGS platforms therefore demonstrated to be suitable for large scale identification of SNPs in transcribed regions of non-model species, although the lack of a reference genome profoundly affects the genotyping success rate. The overall efficiency, however, can be improved using strict quality and filtering criteria for SNP selection (sequence quality, intron-exon structure, target region score).

Natural Killer Cell-like signature observed in locally advanced rectal cancer after neoadjuvant chemoradiotherapy in a tumour regression grade dependent manner

Background: Around 10-15% of patients with locally advanced rectal cancer (LARC) undergo a pathologically complete response (TRG4) to neoadjuvant chemoradiotherapy; the rest of patients exhibit a spectrum of tumour regression (TRG1-3). Understanding therapy-related genomic alterations may help us to identify underlying biology or novel targets associated with response that could increase the efficacy of therapy in patients that do not benefit from the current standard of care.
Methods: 48 FFPE rectal cancer biopsies and matched resections were analysed using the WG-DASL HumanHT-12_v4 Beadchip array on the illumina iScan. Bioinformatic analysis was conducted in Partek genomics suite and R studio. Limma and glmnet packages were used to identify genes differentially expressed between tumour regression grades. Validation of microarray results will be carried out using IHC, RNAscope and RT-PCR.
Results: Immune response genes were observed from supervised analysis of the biopsies which may have predictive value. Differential gene expression from the resections as well as pre and post therapy analysis revealed induction of genes in a tumour regression dependent manner. Pathway mapping and Gene Ontology analysis of these genes suggested antigen processing and natural killer mediated cytotoxicity respectively. The natural killer-like gene signature was switched off in non-responders and on in the responders. IHC has confirmed the presence of Natural killer cells through CD56+ staining.
Conclusion: Identification of NK cell genes and CD56+ cells in patients responding to neoadjuvant chemoradiotherapy warrants further investigation into their association with tumour regression grade in LARC. NK cells are known to lyse malignant cells and determining whether their presence is a cause or consequence of response is crucial. Interrogation of the cytokines upregulated in our NK-like signature will help guide future in vitro models.

Biological phosphorus removal during high-rate, low-temperature, anaerobic digestion of wastewater.

We report, for the first time, extensive biologically-mediated phosphate removal from wastewater during high-rate anaerobic digestion (AD). A hybrid sludge bed/fixed-film (packed pumice stone) reactor was employed for low-temperature (12°C) anaerobic treatment of synthetic sewage wastewater. Successful phosphate removal from the wastewater (up to 78% of influent phosphate) was observed, mediated by biofilms in the reactor. Scanning electron microscopy and energy dispersive X-ray analysis revealed the accumulation of elemental phosphorus (~2%) within the sludge bed and fixed-film biofilms. 4’, 6-diamidino-2-phenylindole (DAPI) staining indicated phosphorus accumulation was biological in nature and mediated through the formation of intracellular inorganic polyphosphate (polyP) granules within these biofilms. DAPI staining further indicated that polyP accumulation was rarely associated with free cells. Efficient and consistent chemical oxygen demand (COD) removal was recorded, throughout the 732-day trial, at applied organic loading rates between 0.4-1.5 kg COD m-3 d-1 and hydraulic retention times of 8-24 hours, while phosphate removal efficiency ranged from 28-78% on average per phase. Analysis of protein hydrolysis kinetics and the methanogenic activity profiles of the biomass revealed the development, at 12˚C, of active hydrolytic and methanogenic populations. Temporal microbial changes were monitored using Illumina Miseq analysis of bacterial and archaeal 16S rRNA gene sequences. The dominant bacterial phyla present in the biomass at the conclusion of the trial were the Proteobacteria and Firmicutes and the dominant archaeal genus was Methanosaeta. Trichococcus and Flavobacterium populations, previously associated with low temperature protein degradation, developed in the reactor biomass. The presence of previously characterised polyphosphate accumulating organisms (PAOs) such as Rhodocyclus, Chromatiales, Actinobacter and Acinetobacter was recorded at low numbers. However, it is unknown as yet if these were responsible for the luxury polyP uptake observed in this system. The possibility of efficient phosphate removal and recovery from wastewater during AD would represent a major advance in the scope for widespread application of anaerobic wastewater treatment technologies.

Community dynamics and the lower airway microbiota in stable chronic obstructive pulmonary disease, smokers and healthy non-smokers

RATIONALE: The role bacteria play in the progression of COPD has increasingly been highlighted in recent years. However, the microbial community complexity in the lower airways of patients with COPD is poorly characterised.

OBJECTIVES: To compare the lower airway microbiota in patients with COPD, smokers and non-smokers.

METHODS: Bronchial wash samples from adults with COPD (n=18), smokers with no airways disease (n=8) and healthy individuals (n=11) were analysed by extended-culture and culture-independent Illumina MiSeq sequencing. We determined aerobic and anaerobic microbiota load and evaluated differences in bacteria associated with the three cohorts. Culture-independent analysis was used to determine differences in microbiota between comparison groups including taxonomic richness, diversity, relative abundance, 'core' microbiota and co-occurrence.

MEASUREMENT AND MAIN RESULTS: Extended-culture showed no difference in total load of aerobic and anaerobic bacteria between the three cohorts. Culture-independent analysis revealed that the prevalence of members of Pseudomonas spp. was greater in the lower airways of patients with COPD; however, the majority of the sequence reads for this taxa were attributed to three patients. Furthermore, members of Bacteroidetes, such as Prevotella spp., were observed to be greater in the 'healthy' comparison groups. Community diversity (α and β) was significantly less in COPD compared with healthy groups. Co-occurrence of bacterial taxa and the observation of a putative 'core' community within the lower airways were also observed.

CONCLUSIONS: Microbial community composition in the lower airways of patients with COPD is significantly different to that found in smokers and non-smokers, indicating that a component of the disease is associated with changes in microbiological status.

Metagenomic characterisation of the viral community of lough neagh, the largest freshwater lake in Ireland

Lough Neagh is the largest and the most economically important lake in Ireland. It is also one of the most nutrient rich amongst the world's major lakes. In this study, 16S rRNA analysis of total metagenomic DNA from the water column of Lough Neagh has revealed a high proportion of Cyanobacteria and low levels of Actinobacteria, Acidobacteria, Chloroflexi, and Firmicutes. The planktonic virome of Lough Neagh has been sequenced and 2,298,791 2×300 bp Illumina reads analysed. Comparison with previously characterised lakes demonstrates that the Lough Neagh viral community has the highest level of sequence diversity. Only about 15% of reads had homologs in the RefSeq database and tailed bacteriophages (Caudovirales) were identified as a major grouping. Within the Caudovirales, the Podoviridae and Siphoviridae were the two most dominant families (34.3% and 32.8% of the reads with sequence homology to the RefSeq database), while ssDNA bacteriophages constituted less than 1% of the virome. Putative cyanophages were found to be abundant. 66,450 viral contigs were assembled with the largest one being 58,805 bp; its existence, and that of another 34,467 bp contig, in the water column was confirmed. Analysis of the contigs confirmed the high abundance of cyanophages in the water column.

Analysis of RNA-seq data from the interaction of Coffea spp. - Colletotrichum kahawae

Tese de mestrado em Bioinformática e Biologia Computacional (Bioinformática), apresentada à Universidade de Lisboa, através da Faculdade de Ciências, 2014

Detection of MegaTAL-induced HIV pol Mutation By Droplet Digital PCR

Thesis (Master's)--University of Washington, 2015