Genome-wide association of early-onset myocardial infarction with single nucleotide polymorphisms and copy number variants

We conducted a genome-wide association study testing single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) for association with early-onset myocardial infarction in 2,967 cases and 3,075 controls. We carried out replication in an independent sample with an effective sample size of up to 19,492. SNPs at nine loci reached genome-wide significance: three are newly identified (21q22 near MRPS6-SLC5A3-KCNE2, 6p24 in PHACTR1 and 2q33 in WDR12) and six replicated prior observations1-4 (9p21, 1p13 near CELSR2-PSRC1-SORT1, 10q11 near CXCL12, 1q41 in MIA3, 19p13 near LDLR and 1p32 near PCSK9). We tested 554 common copy number polymorphisms (>1% allele frequency) and none met the pre-specified threshold for replication (P < 10-3). We identified 8,065 rare CNVs but did not detect a greater CNV burden in cases compared to controls, in genes compared to the genome as a whole, or at any individual locus. SNPs at nine loci were reproducibly associated with myocardial infarction, but tests of common and rare CNVs failed to identify additional associations with myocardial infarction risk.

Deep sequencing reveals predominant expression of miR-21 amongst the small non-coding RNAs in retinal microvascular endothelial cells

The retinal vascular endothelium is essential for angiogenesis and is involved in maintaining barrier selectivity and vascular tone. The aim of this study was to identify and quantify microRNAs and other small regulatory non-coding RNAs (ncRNAs) which may regulate these crucial functions. Primary bovine retinal microvascular endothelial cells (RMECs) provide a well-characterized in vitro system for studying angiogenesis. RNA extracted from RMECs was used to prepare a small RNA library for deep sequencing (Illumina Genome Analyzer). A total of 6.8 million reads were mapped to 250 known microRNAs in miRBase (release 16). In many cases, the most frequent isomiR differed from the sequence reported in miRBase. In addition, five novel microRNAs, 13 novel bovine orthologs of known human microRNAs and multiple new members of the miR-2284/2285 family were detected. Several similar to 30 nucleotide sno-miRNAs were identified, with the most highly expressed being derived from snoRNA U78. Highly expressed microRNAs previously associated with endothelial cells included miR-126 and miR-378, but the most highly expressed was miR-21, comprising more than one-third of all mapped reads. Inhibition of miR-21 with an LNA inhibitor significantly reduced proliferation, migration, and tube-forming capacity of RMECs. The independence from prior sequence knowledge provided by deep sequencing facilitates analysis of novel microRNAs and other small RNAs. This approach also enables quantitative evaluation of microRNA expression, which has highlighted the predominance of a small number of microRNAs in RMECs. Knockdown of miR-21 suggests a role for this microRNA in regulation of angiogenesis in the retinal microvasculature. J. Cell. Biochem. 113: 20982111, 2012. (C) 2012 Wiley Periodicals, Inc.

Genome-wide analysis of p63 binding sites identifies AP-2 factors as co-regulators of epidermal differentiation

The p63 transcription factor (TP63) is critical in development, growth and differentiation of stratifying epithelia. This is highlighted by the severity of congenital abnormalities caused by TP63 mutations in humans, the dramatic phenotypes in knockout mice and de-regulation of TP63 expression in neoplasia altering the tumour suppressive roles of the TP53 family. In order to define the normal role played by TP63 and provide the basis for better understanding how this network is perturbed in disease, we used chromatin immunoprecipitation combined with massively parallel sequencing (ChIP-seq) to identify >7500 high-confidence TP63-binding regions across the entire genome, in primary human neonatal foreskin keratinocytes (HFKs). Using integrative strategies, we demonstrate that only a subset of these sites are bound by TP53 in response to DNA damage. We identify a role for TP63 in transcriptional regulation of multiple genes genetically linked to cleft palate and identify AP-2alpha (TFAP2A) as a co-regulator of a subset of these genes. We further demonstrate that AP-2gamma (TFAP2C) can bind a subset of these regions and that acute depletion of either TFAP2A or TFAP2C alone is sufficient to reduce terminal differentiation of organotypic epidermal skin equivalents, indicating overlapping physiological functions with TP63.

Sequencing and its consequences:Path dependence and the relationships between genetics and medicalization

Both advocacy for and critiques of the Human Genome Project assume a self-sustaining relationship between genetics and. medicalization. However, this assumption ignores the ways in which the meanings of genetic research are conditional on its position in sequences of events. Based, on analyses of three conditions for which at least one putative gene or genetic marker has been identified, this article argues that critical junctures in the institutional stabilization of phenotypes and the mechanisms that sustain such classifications over time configure the practices and meanings of genetic research. Path dependence is critical to understanding the lack of consistent fit between genetics and medlcalization.

Validation of Next Generation Sequencing Technologies in Comparison to Current Diagnostic Gold Standards for BRAF, EGFR and KRAS Mutational Analysis

Next Generation Sequencing (NGS) has the potential of becoming an important tool in clinical diagnosis and therapeutic decision-making in oncology owing to its enhanced sensitivity in DNA mutation detection, fast-turnaround of samples in comparison to current gold standard methods and the potential to sequence a large number of cancer-driving genes at the one time. We aim to test the diagnostic accuracy of current NGS technology in the analysis of mutations that represent current standard-of-care, and its reliability to generate concomitant information on other key genes in human oncogenesis. Thirteen clinical samples (8 lung adenocarcinomas, 3 colon carcinomas and 2 malignant melanomas) already genotyped for EGFR, KRAS and BRAF mutations by current standard-of-care methods (Sanger Sequencing and q-PCR), were analysed for detection of mutations in the same three genes using two NGS platforms and an additional 43 genes with one of these platforms. The results were analysed using closed platform-specific proprietary bioinformatics software as well as open third party applications. Our results indicate that the existing format of the NGS technology performed well in detecting the clinically relevant mutations stated above but may not be reliable for a broader unsupervised analysis of the wider genome in its current design. Our study represents a diagnostically lead validation of the major strengths and weaknesses of this technology before consideration for diagnostic use.

Genome-wide association study of urinary albumin excretion rate in patients with type 1 diabetes

Aims/hypothesis: An abnormal urinary albumin excretion rate (AER) is often the first clinically detectable manifestation of diabetic nephropathy. Our aim was to estimate the heritability and to detect genetic variation associated with elevated AER in patients with type 1 diabetes. 
 Methods: The discovery phase genome-wide association study (GWAS) included 1,925 patients with type 1 diabetes and with data on 24 h AER. AER was analysed as a continuous trait and the analysis was stratified by the use of antihypertensive medication. Signals with a p value <10^{−4 were followed up in 3,750 additional patients with type 1 diabetes from seven studies.} 
 Results: The narrow-sense heritability, captured with our genotyping platform, was estimated to explain 27.3% of the total AER variability, and 37.6% after adjustment for covariates. In the discovery stage, five single nucleotide polymorphisms in the GLRA3 gene were strongly associated with albuminuria (p < 5 × 10^{−8). In the replication group, a nominally significant association (p = 0.035) was observed between albuminuria and rs1564939 in GLRA3, but this was in the opposite direction. Sequencing of the surrounding genetic region in 48 Finnish and 48 UK individuals supported the possibility that population-specific rare variants contribute to the synthetic association observed at the common variants in GLRA3. The strongest replication (p = 0.026) was obtained for rs2410601 between the PSD3 and SH2D4A genes. Pathway analysis highlighted natural killer cell mediated immunity processes.} 
 Conclusions/interpretation: This study suggests novel pathways and molecular mechanisms for the pathogenesis of albuminuria in type 1 diabetes.

Genome-wide characterization reveals complex interplay between TP53 and TP63 in response to genotoxic stress

In response to genotoxic stress the TP53 tumour suppressor activates target gene expression to induce cell cycle arrest or apoptosis depending on the extent of DNA damage. These canonical activities can be repressed by TP63 in normal stratifying epithelia to maintain proliferative capacity or drive proliferation of squamous cell carcinomas, where TP63 is frequently overexpressed/amplified. Here we use ChIP-sequencing, integrated with microarray analysis, to define the genome-wide interplay between TP53 and TP63 in response to genotoxic stress in normal cells. We reveal that TP53 and TP63 bind to overlapping, but distinct cistromes of sites through utilization of distinctive consensus motifs and that TP53 is constitutively bound to a number of sites. We demonstrate that cisplatin and adriamycin elicit distinct effects on TP53 and TP63 binding events, through which TP53 can induce or repress transcription of an extensive network of genes by direct binding and/or modulation of TP63 activity. Collectively, this results in a global TP53-dependent repression of cell cycle progression, mitosis and DNA damage repair concomitant with activation of anti-proliferative and pro-apoptotic canonical target genes. Further analyses reveal that in the absence of genotoxic stress TP63 plays an important role in maintaining expression of DNA repair genes, loss of which results in defective repair.

Genome Wide Association Mapping of Grain Arsenic, Copper, Molybdenum and Zinc in Rice (Oryza sativa L.) Grown at Four International Field Sites

The mineral concentrations in cereals are important for human health, especially for individuals who consume a cereal subsistence diet. A number of elements, such as zinc, are required within the diet, while some elements are toxic to humans, for example arsenic. In this study we carry out genome-wide association (GWA) mapping of grain concentrations of arsenic, copper, molybdenum and zinc in brown rice using an established rice diversity panel of,300 accessions and 36.9 k single nucleotide polymorphisms (SNPs). The study was performed across five environments: one field site in Bangladesh, one in China and two in the US, with one of the US sites repeated over two years. GWA mapping on the whole dataset and on separate subpopulations of rice revealed a large number of loci significantly associated with variation in grain arsenic, copper, molybdenum and zinc. Seventeen of these loci were detected in data obtained from grain cultivated in more than one field location, and six co-localise with previously identified quantitative trait loci. Additionally, a number of candidate genes for the uptake or transport of these elements were located near significantly associated SNPs (within 200 kb, the estimated global linkage disequilibrium previously employed in this rice panel). This analysis highlights a number of genomic regions and candidate genes for further analysis as well as the challenges faced when mapping environmentally-variable traits in a highly genetically structured diversity panel.

Characterisation of a novel panel of polymorphic microsatellite loci for the liver fluke, Fasciola hepatica, using a next generation sequencing approach

The liver fluke, Fasciola hepatica is an economically important pathogen of sheep and cattle and has been described by the WHO as a re-emerging zoonosis. Control is heavily reliant on the use of drugs, particularly triclabendazole and as a result resistance has now emerged. The population structure of F. hepatica is not well known, yet it can impact on host-parasite interactions and parasite control with drugs, particularly regarding the spread of triclabendazole resistance. We have identified 2448 potential microsatellites from 83Mb of F. hepatica genome sequence using msatfinder. Thirty-five loci were developed and optimised for microsatellite PCR, resulting in a panel of 15 polymorphic loci, with a range of three to 15 alleles. This panel was validated on genomic DNA from 46 adult F. hepatica; 38 liver flukes sourced from a Northwest abattoir, UK and 8 liver flukes from an established isolate (Shrewsbury; Ridgeway Research). Evidence for null alleles was found at four loci (Fh_1, Fh_8, Fh_13 and Fh_14), which showed markedly higher levels of homozygosity than the remaining 11 loci. Of the 38 liver flukes isolated from cattle livers (n=10) at the abattoir, 37 genotypes were identified. Using a multiplex approach all 15 loci could be amplified from several life cycle stages that typically yield low amounts of DNA, including metacercariae, the infective life cycle stage present on pasture, highlighting the utility of this multiplex microsatellite panel. This study reports the largest panel of microsatellite markers available to date for population studies of F. hepatica and the first multiplex panel of microsatellite markers that can be used for several life cycle stages.

Genome analysis of a simultaneously predatory and prey-independent, novel Bdellovibrio bacteriovorus from the River Tiber, supports in silico predictions of both ancient and recent lateral gene transfer from diverse bacteria

BACKGROUND: Evolution equipped Bdellovibrio bacteriovorus predatory bacteria to invade other bacteria, digesting and replicating, sealed within them thus preventing nutrient-sharing with organisms in the surrounding environment. Bdellovibrio were previously described as "obligate predators" because only by mutations, often in gene bd0108, are 1 in ~1x10(7) of predatory lab strains of Bdellovibrio converted to prey-independent growth. A previous genomic analysis of B. bacteriovorus strain HD100 suggested that predatory consumption of prey DNA by lytic enzymes made Bdellovibrio less likely than other bacteria to acquire DNA by lateral gene transfer (LGT). However the Doolittle and Pan groups predicted, in silico, both ancient and recent lateral gene transfer into the B. bacteriovorus HD100 genome.

RESULTS: To test these predictions, we isolated a predatory bacterium from the River Tiber- a good potential source of LGT as it is rich in diverse bacteria and organic pollutants- by enrichment culturing with E. coli prey cells. The isolate was identified as B. bacteriovorus and named as strain Tiberius. Unusually, this Tiberius strain showed simultaneous prey-independent growth on organic nutrients and predatory growth on live prey. Despite the prey-independent growth, the homolog of bd0108 did not have typical prey-independent-type mutations. The dual growth mode may reflect the high carbon content of the river, and gives B. bacteriovorus Tiberius extended non-predatory contact with the other bacteria present. The HD100 and Tiberius genomes were extensively syntenic despite their different cultured-terrestrial/freshly-isolated aquatic histories; but there were significant differences in gene content indicative of genomic flux and LGT. Gene content comparisons support previously published in silico predictions for LGT in strain HD100 with substantial conservation of genes predicted to have ancient LGT origins but little conservation of AT-rich genes predicted to be recently acquired.

CONCLUSIONS: The natural niche and dual predatory, and prey-independent growth of the B. bacteriovorus Tiberius strain afforded it extensive non-predatory contact with other marine and freshwater bacteria from which LGT is evident in its genome. Thus despite their arsenal of DNA-lytic enzymes; Bdellovibrio are not always predatory in natural niches and their genomes are shaped by acquiring whole genes from other bacteria.

Exome sequencing of prostate cancer supports the hypothesis of independent tumour origins

BACKGROUND: Prostate cancer (PCa) is a clinically and pathologically heterogeneous disease. The rapid development of sequencing technology has the potential to deliver new biomarkers with emphasis on aggressive disease and to revolutionise personalised cancer treatment. However, a prostate harbouring cancer commonly contains multiple separate tumour foci, with the potential to aggravate tumour sampling. The level of intraprostatic tumour heterogeneity remains to be determined.

OBJECTIVE: To determine the level of intraprostatic tumour heterogeneity through genome-wide, high-resolution profiling of multiple tumour samples from the same individual.

DESIGN, SETTINGS, AND PARTICIPANTS: Multiple tumour samples were obtained from four individuals following radical prostatectomy. One individual (SWE-1) contained >70% cancer cells in all tumour samples, whereas the other three (SWE-2 to SWE-4) required the use of laser capture microdissection for tumour cell enrichment. Subsequently, DNA was extracted from all tissue samples, and exome sequencing was performed. All tumour foci of SWE-1 were also profiled using a high-resolution array for the identification of copy number alterations (CNA).

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Shared somatic high-frequency single nucleotide variants (SNV) and CNAs were used to infer the level of intraprostatic tumour heterogeneity.

RESULTS AND LIMITATIONS: No high-frequency mutations, common for the three tumour samples of SWE-1, were identified. Ten randomly chosen positions were validated with Sanger sequencing in all foci, which verified the exome data. The high level of intraprostatic heterogeneity was consistent in all individuals. In total, three out of four individuals harboured tumours without an apparent common somatic denominator. Although we cannot exclude the presence of common structural rearrangements, a high-density array was used for the detection of deletions and amplifications in SWE-1, which agreed with the exome data.

CONCLUSIONS: We present evidence for the presence of somatically independent tumours within the same prostate. This finding will have implications for personalised cancer treatment and biomarker discovery.

An integrated and comparative approach to genetic selection of the aquaculture species Dicentrarchus labrax

The European sea bass, Dicentrarchus labrax, is one of the most important marine species cultivated in Southern Europe and has not benefited from selective breeding. One of the major goals in the sea bass (D. labrax) aquaculture industry is to understand and control the complexity of growth associated traits. The aim of the methodology developed for the studies reported in the thesis was not only to establish genetic and genomic resources for sea bass, but to also develop a conceptual strategy to efficiently create knowledge in a research environment that can easily be transferred to the aquaculture industry. The strategy involved; i) establishing an annotated sea bass transcriptome and then using it to, ii) identify new genetic markers for target QTL regions so that, iii) new QTL analysis could be performed and marker based resolution of the DNA regions of interest increased, and then iv) to merge the linkage map and the physical map in order to map the QTL confidence intervals to the sea bass genome and identify genes underlying the targeted traits. Finally to test if genes in the QTL regions that are candidates for divergent growth phenotypes have modified patterns of transcription that reflects the modified whole organism physiology SuperSAGE-SOLiD4 gene expression was used with sea bass with high growth heterogeneity. The SuperSAGE contributed to significantly increase the transcriptome information for sea bass muscle, brain and liver and also led to the identification of putative candidate genes lying in the genomic region of growth related QTL. Lastly all differentially expressed transcripts in brain, liver and muscle of the European sea bass with divergent specific growth rates were mapped to gene pathways and networks and the regulatory pathways most affected identified and established the tissue specific changes underlying the divergent SGR. Owing to the importance of European sea bass to Mediterranean aquaculture and the developed genomics resources from the present thesis and from other studies it should be possible to implement genetic selection programs using marker assisted selection.

Next-generation sequencing is highly sensitive for the detection of beta-catenin mutations in desmoid-type fibromatoses

Desmoid-type fibromatoses are locally aggressive and frequently recurrent tumours, and an accurate diagnosis is essential for patient management. The majority of sporadic lesions harbour beta-catenin (CTNNB1) mutations. We used next-generation sequencing to detect CTNNB1 mutations and to compare the sensitivity and specificity of next-generation sequencing with currently employed mutation detection techniques: mutation-specific restriction enzyme digestion and polymerase chain reaction amplification. DNA was extracted from formalin-fixed paraffin-embedded needle biopsy or resection tissue sections from 144 patients with sporadic desmoid-type fibromatoses, four patients with syndrome-related desmoid-type fibromatoses and 11 morphological mimics. Two primer pairs were designed for CTNNB1 mutation hotspots. Using ≥10 ng of DNA, libraries were generated by Fluidigm and sequenced on the Ion Torrent Personal Genome Machine. Next-generation sequencing had a sensitivity of 92.36 % (133/144, 95 % CIs: 86.74 to 96.12 %) and a specificity of 100 % for the detection of CTNNB1 mutations in desmoid-type fibromatoses-like spindle cell lesions. All mutations detected by mutation-specific restriction enzyme digestion were identified by next-generation sequencing. Next-generation sequencing identified additional mutations in 11 tumours that were not detected by mutation-specific restriction enzyme digestion, two of which have not been previously described. Next-generation sequencing is highly sensitive for the detection of CTNNB1 mutations. This multiplex assay has the advantage of detecting additional mutations compared to those detected by mutation-specific restriction enzyme digestion (sensitivity 82.41 %). The technology requires minimal DNA and is time- and cost-efficient.