Generation of an epigenetic signature by chronic hypoxia in prostate cells

Increasing levels of tissue hypoxia have been reported as a natural feature of the aging prostate gland and may be a risk factor for the development of prostate cancer. In this study, we have used PwR-1E benign prostate epithelial cells and an equivalently aged hypoxia-adapted PwR-1E sub-line to identify phenotypic and epigenetic consequences of chronic hypoxia in prostate cells. We have identified a significantly altered cellular phenotype in response to chronic hypoxia as characterized by increased receptor-mediated apoptotic resistance, the induction of cellular senescence, increased invasion and the increased secretion of IL-1 beta, IL6, IL8 and TNFalpha cytokines. In association with these phenotypic changes and the absence of HIF-1 alpha protein expression, we have demonstrated significant increases in global levels of DNA methylation and H3K9 histone acetylation in these cells, concomitant with the increased expression of DNA methyltransferase DMNT3b and gene-specific changes in DNA methylation at key imprinting loci. In conclusion, we have demonstrated a genome-wide adjustment of DNA methylation and histone acetylation under chronic hypoxic conditions in the prostate. These epigenetic signatures may represent an additional mechanism to promote and maintain a hypoxic-adapted cellular phenotype with a potential role in tumour development.

Alleles that increase risk for type 2 diabetes mellitus are not associated with increased risk for Alzheimer's disease

Although epidemiological studies suggest that type 2 diabetes mellitus (T2DM) increases the risk of late-onset Alzheimer's disease (LOAD), the biological basis of this relationship is not well understood. The aim of this study was to examine the genetic comorbidity between the 2 disorders and to investigate whether genetic liability to T2DM, estimated by a genotype risk scores based on T2DM associated loci, is associated with increased risk of LOAD. This study was performed in 2 stages. In stage 1, we combined genotypes for the top 15 T2DM-associated polymorphisms drawn from approximately 3000 individuals (1349 cases and 1351 control subjects) with extracted and/or imputed data from 6 genome-wide studies (>10,000 individuals; 4507 cases, 2183 controls, 4989 population controls) to form a genotype risk score and examined if this was associated with increased LOAD risk in a combined meta-analysis. In stage 2, we investigated the association of LOAD with an expanded T2DM score made of 45 well-established variants drawn from the 6 genome-wide studies. Results were combined in a meta-analysis. Both stage 1 and stage 2 T2DM risk scores were not associated with LOAD risk (odds ratio = 0.988; 95% confidence interval, 0.972-1.004; p = 0.144 and odds ratio = 0.993; 95% confidence interval, 0.983-1.003; p = 0.149 per allele, respectively). Contrary to expectation, genotype risk scores based on established T2DM candidates were not associated with increased risk of LOAD. The observed epidemiological associations between T2DM and LOAD could therefore be a consequence of secondary disease processes, pleiotropic mechanisms, and/or common environmental risk factors. Future work should focus on well-characterized longitudinal cohorts with extensive phenotypic and genetic data relevant to both LOAD and T2DM.

DATA DRIVEN APPROACHES TO IDENTIFY DETERMINANTS OF HEART DISEASES AND CANCER RESISTANCE

Cancer and cardio-vascular diseases are the leading causes of death world-wide. Caused by systemic genetic and molecular disruptions in cells, these disorders are the manifestation of profound disturbance of normal cellular homeostasis. People suffering or at high risk for these disorders need early diagnosis and personalized therapeutic intervention. Successful implementation of such clinical measures can significantly improve global health. However, development of effective therapies is hindered by the challenges in identifying genetic and molecular determinants of the onset of diseases; and in cases where therapies already exist, the main challenge is to identify molecular determinants that drive resistance to the therapies. Due to the progress in sequencing technologies, the access to a large genome-wide biological data is now extended far beyond few experimental labs to the global research community. The unprecedented availability of the data has revolutionized the capabilities of computational researchers, enabling them to collaboratively address the long standing problems from many different perspectives. Likewise, this thesis tackles the two main public health related challenges using data driven approaches. Numerous association studies have been proposed to identify genomic variants that determine disease. However, their clinical utility remains limited due to their inability to distinguish causal variants from associated variants. In the presented thesis, we first propose a simple scheme that improves association studies in supervised fashion and has shown its applicability in identifying genomic regulatory variants associated with hypertension. Next, we propose a coupled Bayesian regression approach -- eQTeL, which leverages epigenetic data to estimate regulatory and gene interaction potential, and identifies combinations of regulatory genomic variants that explain the gene expression variance. On human heart data, eQTeL not only explains a significantly greater proportion of expression variance in samples, but also predicts gene expression more accurately than other methods. We demonstrate that eQTeL accurately detects causal regulatory SNPs by simulation, particularly those with small effect sizes. Using various functional data, we show that SNPs detected by eQTeL are enriched for allele-specific protein binding and histone modifications, which potentially disrupt binding of core cardiac transcription factors and are spatially proximal to their target. eQTeL SNPs capture a substantial proportion of genetic determinants of expression variance and we estimate that 58% of these SNPs are putatively causal. The challenge of identifying molecular determinants of cancer resistance so far could only be dealt with labor intensive and costly experimental studies, and in case of experimental drugs such studies are infeasible. Here we take a fundamentally different data driven approach to understand the evolving landscape of emerging resistance. We introduce a novel class of genetic interactions termed synthetic rescues (SR) in cancer, which denotes a functional interaction between two genes where a change in the activity of one vulnerable gene (which may be a target of a cancer drug) is lethal, but subsequently altered activity of its partner rescuer gene restores cell viability. Next we describe a comprehensive computational framework --termed INCISOR-- for identifying SR underlying cancer resistance. Applying INCISOR to mine The Cancer Genome Atlas (TCGA), a large collection of cancer patient data, we identified the first pan-cancer SR networks, composed of interactions common to many cancer types. We experimentally test and validate a subset of these interactions involving the master regulator gene mTOR. We find that rescuer genes become increasingly activated as breast cancer progresses, testifying to pervasive ongoing rescue processes. We show that SRs can be utilized to successfully predict patients' survival and response to the majority of current cancer drugs, and importantly, for predicting the emergence of drug resistance from the initial tumor biopsy. Our analysis suggests a potential new strategy for enhancing the effectiveness of existing cancer therapies by targeting their rescuer genes to counteract resistance. The thesis provides statistical frameworks that can harness ever increasing high throughput genomic data to address challenges in determining the molecular underpinnings of hypertension, cardiovascular disease and cancer resistance. We discover novel molecular mechanistic insights that will advance the progress in early disease prevention and personalized therapeutics. Our analyses sheds light on the fundamental biological understanding of gene regulation and interaction, and opens up exciting avenues of translational applications in risk prediction and therapeutics.

Factor analysis applied to genome prediction for high-dimensional phenotypes in pigs.

The aim of the present study was to propose and evaluate the use of factor analysis (FA) in obtaining latent variables (factors) that represent a set of pig traits simultaneously, for use in genome-wide selection (GWS) studies. We used crosses between outbred F2 populations of Brazilian Piau X commercial pigs. Data were obtained on 345 F2 pigs, genotyped for 237 SNPs, with 41 traits. FA allowed us to obtain four biologically interpretable factors: ?weight?, ?fat?, ?loin?, and ?performance?. These factors were used as dependent variables in multiple regression models of genomic selection (Bayes A, Bayes B, RR-BLUP, and Bayesian LASSO). The use of FA is presented as an interesting alternative to select individuals for multiple variables simultaneously in GWS studies; accuracy measurements of the factors were similar to those obtained when the original traits were considered individually. The similarities between the top 10% of individuals selected by the factor, and those selected by the individual traits, were also satisfactory. Moreover, the estimated markers effects for the traits were similar to those found for the relevant factor.

Genome-Wiide estimates of coancestry, in breeding and effective population size in the Spanish Holstein population

Estimates of effective population size in the Holstein cattle breed have usually been low despite the large number of animals that constitute this breed. Effective population size is inversely related to the rates at which coancestry and inbreeding increase and these rates have been high as a consequence of intense and accurate selection. Traditionally, coancestry and inbreeding coefficients have been calculated from pedigree data. However, the development of genome-wide single nucleotide polymorphisms has increased the interest of calculating these coefficients from molecular data in order to improve their accuracy. In this study, genomic estimates of coancestry, inbreeding and effective population size were obtained in the Spanish Holstein population and then compared with pedigree-based estimates. A total of 11,135 animals genotyped with the Illumina BovineSNP50 BeadChip were available for the study. After applying filtering criteria, the final genomic dataset included 36,693 autosomal SNPs and 10,569 animals. Pedigree data from those genotyped animals included 31,203 animals. These individuals represented only the last five generations in order to homogenise the amount of pedigree information across animals. Genomic estimates of coancestry and inbreeding were obtained from identity by descent segments (coancestry) or runs of homozygosity (inbreeding). The results indicate that the percentage of variance of pedigree-based coancestry estimates explained by genomic coancestry estimates was higher than that for inbreeding. Estimates of effective population size obtained from genome-wide and pedigree information were consistent and ranged from about 66 to 79. These low values emphasize the need of controlling the rate of increase of coancestry and inbreeding in Holstein selection programmes.

Xylella fastidiosa gene expression analysis by DNA microarrays

Xylella fastidiosa genome sequencing has generated valuable data by identifying genes acting either on metabolic pathways or in associated pathogenicity and virulence. Based on available information on these genes, new strategies for studying their expression patterns, such as microarray technology, were employed. A total of 2,600 primer pairs were synthesized and then used to generate fragments using the PCR technique. The arrays were hybridized against cDNAs labeled during reverse transcription reactions and which were obtained from bacteria grown under two different conditions (liquid XDM2 and liquid BCYE). All data were statistically analyzed to verify which genes were differentially expressed. In addition to exploring conditions for X. fastidiosa genome-wide transcriptome analysis, the present work observed the differential expression of several classes of genes (energy, protein, amino acid and nucleotide metabolism, transport, degradation of substances, toxins and hypothetical proteins, among others). The understanding of expressed genes in these two different media will be useful in comprehending the metabolic characteristics of X. fastidiosa, and in evaluating how important certain genes are for the functioning and survival of these bacteria in plants.

Heteroduplex formation and S1 digestion for mapping alternative splicing sites

The identification of alternatively spliced transcripts has contributed to a better comprehension of developmental mechanisms, tissue-specific physiological processes and human diseases. Polymerase chain reaction amplification of alternatively spliced variants commonly leads to the formation of heteroduplexes as a result of base pairing involving exons common between the two variants. S1 nuclease cleaves single-stranded loops of heteroduplexes and also nicks the opposite DNA strand. In order to establish a strategy for mapping alternative splice-prone sites in the whole transcriptome, we developed a method combining the formation of heteroduplexes between 2 distinct splicing variants and S1 nuclease digestion. For 20 consensuses identified here using this methodology, 5 revealed a conserved splice site after inspection of the cDNA alignment against the human genome (exact splice sites). For 8 other consensuses, conserved splice sites were mapped at 2 to 30 bp from the border, called proximal splice sites; for the other 7 consensuses, conserved splice sites were mapped at 40 to 800 bp, called distal splice sites. These latter cases showed a nonspecific activity of S1 nuclease in digesting double-strand DNA. From the 20 consensuses identified here, 5 were selected for reverse transcription-polymerase chain reaction validation, confirming the splice sites. These data showed the potential of the strategy in mapping splice sites. However, the lack of specificity of the S1 nuclease enzyme is a significant obstacle that impedes the use of this strategy in large-scale studies.

Familial partial epilepsy with variable foci: A new partial epilepsy syndrome with suggestion of linkage to chromosome 2

Familial partial epilepsy with variable foci (FPEVF) joins the recently recognized group of inherited partial epilepsies. We describe an Australian family with 10 individuals with partial seizures over four generations. Detailed electroclinical studies were performed on all affected and 17 clinically unaffected family members. The striking finding was that the clinical features of the seizures and interictal electroencephalographic foci differed among family members and included frontal, temporal, occipital, and centroparietal seizures. Mean age of seizure onset was 13 years (range, 0.75-43 years). Two individuals without seizures had epileptiform abnormalities on electroencephalographic studies. Penetrance of seizures was 62%. A genome-wide search failed to demonstrate definitive linkage, but a suggestion of linkage was found on chromosome 2q with a LOD score of 2.74 at recombination fraction of zero with the marker D2S133. FPEVF differs from the other inherited partial epilepsies where partial seizures in different family members are clinically similar. The inherited nature of this new syndrome may be overlooked because of relatively low penetrance and because of the variability in age at onset and electroclinical features between affected family members.

A novel genetic locus for low renin hypertension: familial hyperaldosteronism type II maps to chromosome 7 (7p22)

Familial hyperaldosteronism type II (FH-II) is caused by adrenocortical hyperplasia or aldosteronoma or both and is frequently transmitted in an autosomal dominant fashion. Unlike FH type I (FI-I-I), which results from fusion of the CYP11B1 and CYP11B2 genes, hyperaldosteronism in FH-II is not glucocorticoid remediable. A large family with FH-II was used for a genome wide search and its members were evaluated by measuring the aldosterone:renin ratio. In those with an increased ratio, FH-II was confirmed by fludrocortisone suppression testing. After excluding most of the genome, genetic linkage was identified with a maximum two point lod score of 3.26 at theta =0, between FH-II in this family and the polymorphic markers D7S511, D7S517, and GATA24F03 on chromosome 7,a region that corresponds to cytogenetic band 7p22. This is the first identified locus for FH-II; its molecular elucidation may provide further insight into the aetiology of primary aldosteronism.

A new clinical and molecular form of Unverricht-Lundborg disease localized by homozygosity mapping

Progressive myoclonus epilepsy (PME) has a number of causes, of which Unverricht-Lundborg disease (ULD) is the most common. ULD has previously been mapped to a locus on chromosome 21 (EPM1). Subsequently, mutations in the cystatin B gene have been found in most cases. In the present work we identified an inbred Arab family with a clinical pattern compatible with ULD, but mutations in the cystatin B gene were absent. We sought to characterize the clinical and molecular features of the disorder. The family was studied by multiple field trips to their town to clarify details of the complex consanguineous relationships and to personally examine the family. DNA was collected for subsequent molecular analyses from 21 individuals. A genome-wide screen was performed using 811 microsatellite markers. Homozygosity mapping was used to identify loci of interest. There were eight affected individuals. Clinical onset was at 7.3 +/- 1.5 years with myoclonic or tonic-clonic seizures. All had myoclonus that progressed in severity over time and seven had tonic-clonic seizures. Ataxia, in addition to myoclonus, occurred in all. Detailed cognitive assessment was not possible, but there was no significant progressive dementia. There was intrafamily variation in severity; three required wheelchairs in adult life; the others could walk unaided. MRI, muscle and skin biopsies on one individual were unremarkable. We mapped the family to a 15-megabase region at the pericentromeric region of chromosome 12 with a maximum lod score of 6.32. Although the phenotype of individual subjects was typical of ULD, the mean age of onset (7.3 years versus 11 years for ULD) was younger. The locus on chromosome 12 does not contain genes for any other form of PME, nor does it have genes known to be related to cystatin B. This represents a new form of PME and we have designated the locus as EPM1B.

Cytogenetic and Molecular Evaluation and 20-Year Follow-Up of a Patient With Ring Chromosome 14

We present a 20-year follow-up on a patient with a ring chromosome 14. The ring chromosome was studied by fluorescence in-situ hybridization (FISH), multiplex-ligation probe amplification (MLPA), and genome wide SNP array, and no deletions of chromosome 14 were detected, although the telomeric repeat sequence was absent from the ring chromosome. The patient had skeletal abnormalities, and susceptibility to infections, as well as seizures and retinal pigmentation, which are commonly found in individuals with a ring 14. Our patient corroborates the idea that even when no genes are lost during ring formation, a complete ring chromosome can produce phenotypic alterations, which presumably result from ring instability or gene silencing due to the new chromosomal architecture. (C) 2010 Wiley-Liss, Inc.

De novo copy number variants identify new genes and loci in isolated sporadic tetralogy of Fallot

Tetralogy of Fallot (TOF), the most common severe congenital heart malformation, occurs sporadically, without other anomaly, and from unknown cause in 70% of cases. Through a genome-wide survey of 114 subjects with TOF and their unaffected parents, we identified 11 de novo copy number variants (CNVs) that were absent or extremely rare (<0.1%) in 2,265 controls. We then examined a second, independent TOF cohort (n = 398) for additional CNVs at these loci. We identified CNVs at chromosome 1q21.1 in 1% (5/512, P = 0.0002, OR = 22.3) of nonsyndromic sporadic TOF cases. We also identified recurrent CNVs at 3p25.1, 7p21.3 and 22q11.2. CNVs in a single subject with TOF occurred at six loci, two that encode known (NOTCH1, JAG1) disease-associated genes. Our findings predict that at least 10% (4.5-15.5%, 95% confidence interval) of sporadic nonsyndromic TOF cases result from de novo CNVs and suggest that mutations within these loci might be etiologic in other cases of TOF.

Epstein-Barr virus detection in invasive and pre-invasive lesions of the uterine cervix

In the present study, our aim was to investigate whether EBV DNA could be found in association with invasive and pre-invasive cervical cancer lesions. We hypothesize that EBV is not merely a commensal agent when present in malignant cervical lesions. DNA was extracted from cervical scrapings followed by nested PCR-based amplification. The patients were 66 women with high grade cervical intraepithelial neoplasia and 14 women with invasive cervical cancer. The control group consisted of 89 women with a normal Pap smear and colposcopy as well as a negative HPV DNA test. Analysis of our results, in conjunction with the work of other authors, leads us to propose that EBV is not merely a commensal agent when present in malignant cervical lesions. The presence of DNA from EBV is significantly associated with cervical cancer.

A transposon toolkit for gene transfer and mutagenesis in protozoan parasites

Protozoan parasites affect millions of people around the world. Treatment and control of these diseases are complicated partly due to the intricate biology of these organisms. The interactions of species of Plasmodium, Leishmania and trypanosomes with their hosts are mediated by an unusual control of gene expression that is not fully understood. The availability of the genome sequence of these protozoa sets the stage for using more comprehensive, genome-wide strategies to study gene function. Transposons are effective tools for the systematic introduction of genetic alterations and different transposition systems have been adapted to study gene function in these human pathogens. A mariner transposon toolkit for use in vivo or in vitro in Leishmania parasites has been developed and can be used in a variety of applications. These modified mariner elements not only permit the inactivation of genes, but also mediate the rescue of translational gene fusions, bringing a major contribution to the investigation of Leishmania gene function. The piggyBac and Tn5 transposons have also been shown to mobilize across Plasmodium spp. genomes circumventing the current limitations in the genetic manipulation of these organisms.