Genome-wide expression profiling identifies deregulated miRNAs in malignant astrocytoma

Malignant astrocytoma includes anaplastic astrocytoma (grade III) and glioblastoma (grade IV). Among them, glioblastoma is the most common primary brain tumor with dismal responses to all therapeutic modalities. We performed a large-scale, genome-wide microRNA (miRNA) (n=756) expression profiling of 26 glioblastoma, 13 anaplastic astrocytoma and 7 normal brain samples with an aim to find deregulated miRNA in malignant astrocytoma. We identified several differentially regulated miRNAs between these groups, which could differentiate glioma grades and normal brain as recognized by PCA. More importantly, we identified a most discriminatory 23-miRNA expression signature, by using PAM, which precisely distinguished glioblastoma from anaplastic astrocytoma with an accuracy of 95%. The differential expression pattern of nine miRNAs was further validated by real-time RT-PCR on an independent set of malignant astrocytomas (n-72) and normal samples (n=7). Inhibition of two glioblastoma-upregulated miRNAs (miR-21 and miR-23a) and exogenous overexpression of two glioblastoma-downregulated miRNAs (miR-218 and miR-219-5p) resulted in reduced soft agar colony formation but showed varying effects on cell proliferation and chemosensitivity. Thus we have identified the miRNA expression signature for malignant astrocytoma, in particular glioblastoma, and showed the miRNA involvement and their importance in astrocytoma development. Modern Pathology (2010) 23, 1404-1417; doi:10.1038/modpathol.2010.135; published online 13 August 2010

Studies on Causes of Multiple Sclerosis : From Genes to Transcriptome

Multiple sclerosis (MS) is the most common cause of neurological disability in young adults, affecting more than two million people worldwide. It manifests as a chronic inflammation in the central nervous system (CNS) and causes demyelination and neurodegeneration. Depending on the location of the demyelinated plaques and axonal loss, a variety of symptoms can be observed including deficits in vision, coordination, balance and movement. With a typical age of onset at 20-40 years, the social and economic impacts of MS on lives of the patients and their families are considerable. Unfortunately the current treatments are relatively inefficient and the development of more effective treatments has been impeded by our limited understanding of the causes and pathogenesis of MS. Risk of MS is higher in biological relatives of MS patients than in the general population. Twin and adoption studies have shown that familial clustering of MS is explained by shared genetic factors rather than by shared familial environment. While the involvement of the human leukocyte antigen (HLA) genes was first discovered four decades ago, additional genetic risk factors have only recently been identified through genome-wide association studies (GWAS). Current evidence suggests that MS is a highly polygenic disease with perhaps hundreds of common variants with relatively modest effects contributing to susceptibility. Despite extensive research, the majority of these risk factors still remain to be identified. In this thesis the aim was to identify novel genes and pathways involved in MS. Using genome-wide microarray technology, gene expression levels in peripheral blood mononuclear cells (PBMC) from 12 MS patients and 15 controls were profiled and more than 600 genes with altered expression in MS were identified. Three of five selected findings, DEFA1A3, LILRA4 and TNFRSF25, were successfully replicated in an independent sample. Increased expression of DEFA1A3 in MS is a particularly interesting observation, because its elevated levels have previously been reported also in several other autoimmune diseases. A systematic review of seven microarray studies was then performed leading to identification of 229 genes, in which either decreased or increased expression in MS had been reported in at least two studies. In general there was relatively little overlap across the experiments: 11 of the 229 genes had been reported in three studies and only HSPA1A in four studies. Nevertheless, these 229 genes were associated with several immunological pathways including interleukin pathways related to type 2 and type 17 helper T cells and regulatory T cells. However, whether these pathways are involved in causing MS or related to secondary processes activated after disease onset remains to be investigated. The 229 genes were also compared with loci identified in published MS GWASs. Single nucleotide polymorphisms (SNP) in 17 of the 229 loci had been reported to be associated with MS with P-value less than 0.0001 including variants in CXCR4 and SAPS2, which were the only loci where evidence for correlation between the associated variant and gene expression was found. The CXCR4 variant was further tested for association with MS in a large case-control sample and the previously reported suggestive association was replicated (P-value is 0.0004). Finally, common genetic variants in candidate genes, which had been selected on the basis of showing association with other autoimmune diseases (MYO9B) or showing differential expression in MS in our study (DEFA1A3, LILRA4 and TNFRSF25), were tested for association with MS, but no evidence of association was found. In conclusion, through a systematic review of genome-wide expression studies in MS we have identified several promising candidate genes and pathways for future studies. In addition, we have replicated a previously suggested association of a SNP variant upstream of CXCR4 with MS. Keywords: autoimmune disease, common variant, CXCR4, DEFA1A3, HSPA1A,gene expression, genetic association, GWAS, MS, multiple sclerosis, systematic review

LIN28B, LIN28A, KISS1, and KISS1R in idiopathic central precocious puberty

Abstract Background Pubertal timing is a strongly heritable trait, but no single puberty gene has been identified. Thus, the genetic background of idiopathic central precocious puberty (ICPP) is poorly understood. Overall, the genetic modulation of pubertal onset most likely arises from the additive effect of multiple genes, but also monogenic causes of ICPP probably exist, as cases of familial ICPP have been reported. Mutations in KISS1 and KISSR, coding for kisspeptin and its receptor, involved in GnRH secretion and puberty onset, have been suggested causative for monogenic ICPP. Variation in LIN28B was associated with timing of puberty in genome-wide association (GWA) studies. LIN28B is a human ortholog of the gene that controls, through microRNAs, developmental timing in C. elegans. In addition, Lin28a transgenic mice manifest the puberty phenotypes identified in the human GWAS. Thus, both LIN28B and LIN28A may have a role in pubertal development and are good candidate genes for monogenic ICPP. Methods Thirty girls with ICPP were included in the study. ICPP was defined by pubertal onset before 8 yrs of age, and a pubertal LH response to GnRH testing. The coding regions of LIN28B, LIN28A, KISS1, and KISS1R were sequenced. The missense change in LIN28B was also screened in 132 control subjects. Results No rare variants were detected in KISS1 or KISS1R in the 30 subjects with ICPP. In LIN28B, one missense change, His199Arg, was found in one subject with ICPP. However, this variant was also detected in one of the 132 controls. No variation in LIN28A was found. Conclusions We did not find any evidence that mutations in LIN28B or LIN28A would underlie ICPP. In addition, we confirmed that mutations in KISS1 and KISS1R are not a common cause for ICPP.

Involvement of ANXA5 and ILKAP in Susceptibility to Malignant Melanoma

Single nucleotide-polymorphisms (SNPs) are a source of diversity among human population, which may be responsible for the different individual susceptibility to diseases and/or response to drugs, among other phenotypic traits. Several low penetrance susceptibility genes associated with malignant melanoma (MM) have been described, including genes related to pigmentation, DNA damage repair and oxidative stress pathways. In the present work, we conducted a candidate gene association study based on proteins and genes whose expression we had detected altered in melanoma cell lines as compared to normal melanocytes. The result was the selection of 88 loci and 384 SNPs, of which 314 fulfilled our quality criteria for a case-control association study. The SNP rs6854854 in ANXA5 was statistically significant after conservative Bonferroni correction when 464 melanoma patients and 400 controls were analyzed in a discovery Phase I. However, this finding could not be replicated in the validation phase, perhaps because the minor allele frequency of SNP rs6854854 varies depending on the geographical region considered. Additionally, a second SNP (rs6431588) located on ILKAP was found to be associated with melanoma after considering a combined set of 1,883 MM cases and 1,358 disease-free controls. The OR was 1.29 (95% CI 1.12-1.48; p-value= 4x10(-4)). Both SNPs, rs6854854 in ANXA5 and rs6431588 in ILKAP, show population structure, which, assuming that the Spanish population is not significantly structured, suggests a role of these loci on a specific genetic adaptation to different environmental conditions. Furthermore, the biological relevance of these genes in MM is supported by in vitro experiments, which show a decrease in the transcription levels of ANXA5 and ILKAP in melanoma cells compared to normal melanocytes.

Genetic Variants in MiRNA Processing Genes and Pre-MiRNAs Are Associated with the Risk of Chronic Lymphocytic Leukemia

Genome wide association studies (GWAS) have identified several low-penetrance susceptibility alleles in chronic lymphocytic leukemia (CLL). Nevertheless, these studies scarcely study regions that are implicated in non-coding molecules such as microRNAs (miRNAs). Abnormalities in miRNAs, as altered expression patterns and mutations, have been described in CLL, suggesting their implication in the development of the disease. Genetic variations in miRNAs can affect levels of miRNA expression if present in pre-miRNAs and in miRNA biogenesis genes or alter miRNA function if present in both target mRNA and miRNA sequences. Therefore, the present study aimed to evaluate whether polymorphisms in pre-miRNAs, and/or miRNA processing genes contribute to predisposition for CLL. A total of 91 SNPs in 107 CLL patients and 350 cancer-free controls were successfully analyzed using TaqMan Open Array technology. We found nine statistically significant associations with CLL risk after FDR correction, seven in miRNA processing genes (rs3805500 and rs6877842 in DROSHA, rs1057035 in DICER1, rs17676986 in SND1, rs9611280 in TNRC6B, rs784567 in TRBP and rs11866002 in CNOT1) and two in pre-miRNAs (rs11614913 in miR196a2 and rs2114358 in miR1206). These findings suggest that polymorphisms in genes involved in miRNAs biogenesis pathway as well as in pre-miRNAs contribute to the risk of CLL. Large-scale studies are needed to validate the current findings.

Refined genome-wide comparative map of the domestic horse, donkey and human based on cross-species chromosome painting: insight into the occasional fertility of mules

We have made a complete set of painting probes for the domestic horse by degenerate oligonucleotide-primed PCR amplification of flow-sorted horse chromosomes. The horse probes, together with a full set of those available for human, were hybridized onto metaphase chromosomes of human, horse and mule. Based on the hybridization results, we have generated genome-wide comparative chromosome maps involving the domestic horse, donkey and human. These maps define the overall distribution and boundaries of evolutionarily conserved chromosomal segments in the three genomes. Our results shed further light on the karyotypic relationships among these species and, in particular, the chromosomal rearrangements that underlie hybrid sterility and the occasional fertility of mules.

Genome-wide high-resolution mapping and functional analysis of DNA methylation in arabidopsis.

Cytosine methylation is important for transposon silencing and epigenetic regulation of endogenous genes, although the extent to which this DNA modification functions to regulate the genome is still unknown. Here we report the first comprehensive DNA methylation map of an entire genome, at 35 base pair resolution, using the flowering plant Arabidopsis thaliana as a model. We find that pericentromeric heterochromatin, repetitive sequences, and regions producing small interfering RNAs are heavily methylated. Unexpectedly, over one-third of expressed genes contain methylation within transcribed regions, whereas only approximately 5% of genes show methylation within promoter regions. Interestingly, genes methylated in transcribed regions are highly expressed and constitutively active, whereas promoter-methylated genes show a greater degree of tissue-specific expression. Whole-genome tiling-array transcriptional profiling of DNA methyltransferase null mutants identified hundreds of genes and intergenic noncoding RNAs with altered expression levels, many of which may be epigenetically controlled by DNA methylation.

Genome-wide SNP and haplotype analyses reveal a rich history underlying dog domestication

Advances in genome technology have facilitated a new understanding of the historical and genetic processes crucial to rapid phenotypic evolution under domestication(1,2). To understand the process of dog diversification better, we conducted an extensive genome-wide survey of more than 48,000 single nucleotide polymorphisms in dogs and their wild progenitor, the grey wolf. Here we show that dog breeds share a higher proportion of multi-locus haplotypes unique to grey wolves from the Middle East, indicating that they are a dominant source of genetic diversity for dogs rather than wolves from east Asia, as suggested by mitochondrial DNA sequence data(3). Furthermore, we find a surprising correspondence between genetic and phenotypic/functional breed groupings but there are exceptions that suggest phenotypic diversification depended in part on the repeated crossing of individuals with novel phenotypes. Our results show that Middle Eastern wolves were a critical source of genome diversity, although interbreeding with local wolf populations clearly occurred elsewhere in the early history of specific lineages. More recently, the evolution of modern dog breeds seems to have been an iterative process that drew on a limited genetic toolkit to create remarkable phenotypic diversity.

Genome-wide computational identification of microRNAs and their targets in the deep-branching eukaryote Giardia lamblia

Using a combined computational program. we identified 50 potential microRNAs (miRNAs) in Giardia lamblia. one of the most primitive unicellular eukaryotes. These miRNAs are unique to G. lamblia and no homologues have been found in other organisms; miRNAs.

Genome-wide analysis of restriction-modification system in unicellular and filamentous cyanobacteria

Cyanobacteria are an ancient group of gram-negative bacteria with strong genome size variation ranging from 1.6 to 9.1 Mb. Here, we first retrieved all the putative restriction-modification (RM) genes in the draft genome of Spirulina and then performed a range of comparative and bioinformatic analyses on RM genes from unicellular and filamentous cyanobacterial genomes. We have identified 6 gene clusters containing putative Type I RMs and 11 putative Type II RMs or the solitary methyltransferases (MTases). RT-PCR analysis reveals that 6 of 18 MTases are not expressed in Spirulina, whereas one hsdM gene, with a mutated cognate hsdS, was detected to be expressed. Our results indicate that the number of RM genes in filamentous cyanobacteria is significantly higher than in unicellular species, and this expansion of RM systems in filamentous cyanobacteria may be related to their wide range of ecological tolerance. Furthermore, a coevolutionary pattern is found between hsdM and hsdR, with a large number of site pairs positively or negatively correlated, indicating the functional importance of these pairing interactions between their tertiary structures. No evidence for positive selection is found for the majority of RMs, e. g., hsdM, hsdS, hsdR, and Type II restriction endonuclease gene families, while a group of MTases exhibit a remarkable signature of adaptive evolution. Sites and genes identified here to have been under positive selection would provide targets for further research on their structural and functional evaluations.

Genome-wide analysis of restriction-modification system in unicellular and filamentous cyanobacteria

Cyanobacteria are an ancient group of gram-negative bacteria with strong genome size variation ranging from 1.6 to 9.1 Mb. Here, we first retrieved all the putative restriction-modification (RM) genes in the draft genome of Spirulina and then performed a range of comparative and bioinformatic analyses on RM genes from unicellular and filamentous cyanobacterial genomes. We have identified 6 gene clusters containing putative Type I RMs and 11 putative Type II RMs or the solitary methyltransferases (MTases). RT-PCR analysis reveals that 6 of 18 MTases are not expressed in Spirulina, whereas one hsdM gene, with a mutated cognate hsdS, was detected to be expressed. Our results indicate that the number of RM genes in filamentous cyanobacteria is significantly higher than in unicellular species, and this expansion of RM systems in filamentous cyanobacteria may be related to their wide range of ecological tolerance. Furthermore, a coevolutionary pattern is found between hsdM and hsdR, with a large number of site pairs positively or negatively correlated, indicating the functional importance of these pairing interactions between their tertiary structures. No evidence for positive selection is found for the majority of RMs, e. g., hsdM, hsdS, hsdR, and Type II restriction endonuclease gene families, while a group of MTases exhibit a remarkable signature of adaptive evolution. Sites and genes identified here to have been under positive selection would provide targets for further research on their structural and functional evaluations.

Molecular analysis of inflammatory mechanisms associated with Escherichia coli and inflammatory bowel diseases

Inflammatory bowel diseases (IBD), encompasses a range of chronic, immune-mediated inflammatory disorders that are usually classified under two major relapsing conditions, Crohn’s Disease (CD) and ulcerative colitis (UC). Extensive studies in the last decades have suggested that the etiology of IBD involves environmental and genetic factors that lead to dysfunction of epithelial barrier with consequent deregulation of the mucosal immune system and inadequate responses to gut microbiota.Over the last decade, the microbial species that has attracted the most attention, with respect to CD etiology, is Eschericia coli. In CD tissue, E. coli antigens have also been identified in macrophages within the lamina propria, granulomas, and in the germinal centres of mesenteric lymph nodes of patients. They have been shown to adhere to and invade intestinal epithelial cells whilst also being able to extensively replicate within macrophages. Through the work of genome-wide association studies (GWAS), there is growing evidence to suggest that the microbial imbalance between commensal and pathogenic bacteria in the gut is aided by a defect in the innate immune system. Autophagy represents a recently investigated pathway that is believed to contribute to the pathogenesis of CD, with studies identified a variant of the autophagy gene, ATG16L1, as a susceptibility gene. The aim of my thesis was to study the cellular and molecular mechanism promoted by E.coli strains in epithelial cells and to assess their contribution to IBD pathology. To achieve this we focused on developing both an in vitro and in vivo model of AIEC infection. This allowed us to further our knowledge on possible mechanisms utilised by AIEC that promoted their survival, as well as developing a better understanding of host reactions. We demonstrate a new survival mechanism promoted by E.coli HM605, whereby it induces the expression of the anti-apoptotic proteins Bcl-XL and BCL2, all of which is exacerbated in an autophagy deficient system. We have also demonstrated the presence of AIEC-induced inflammasome responses in epithelial cells which are exacerbated in an autophagy deficient system and expression of NOD-like receptors (NLRs) which might mediate inflammasome responses in vivo. Finally, we used the Citrobacter rodentium model of infectious colitis to identify Pellino3 as an important mediator in the NOD2 pathway and regulator of intestinal inflammation. In summary, we have developed robust and versatile models of AIEC infection as well as provide new insights into AIEC mediated survival pathways. The collected data provides a new perception into why AIEC bacteria are able to prosper in conditions associated with Crohn’s disease patients with a defect in autophagy.

Genome-wide hydroxymethylcytosine pattern changes in response to oxidative stress

The TET enzymes convert methylcytosine to the newly discovered base hydroxymethylcytosine. While recent reports suggest that TETs may play a role in response to oxidative stress, this role remains uncertain, and results lack in vivo models. Here we show a global decrease of hydroxymethylcytosine in cells treated with buthionine sulfoximine, and in mice depleted for the major antioxidant enzymes GPx1 and 2. Furthermore, genome-wide profiling revealed differentially hydroxymethylated regions in coding genes, and intriguingly in microRNA genes, both involved in response to oxidative stress. These results thus suggest a profound effect of in vivo oxidative stress on the global hydroxymethylome.

A genome-wide linkage scan for genes controlling variation in urinary albumin excretion in type ll diabetes

The main hallmark of diabetic nephropathy is elevation in urinary albumin excretion. We performed a genome-wide linkage scan in 63 extended families with multiple members with type II diabetes. Urinary albumin excretion, measured as the albumin-to-creatinine ratio (ACR), was determined in 426 diabetic and 431 nondiabetic relatives who were genotyped for 383 markers. The data were analyzed using variance components linkage analysis. Heritability (h2) of ACR was significant in diabetic (h2=0.23, P=0.0007), and nondiabetic (h2=0.39, P=0.0001) relatives. There was no significant difference in genetic variance of ACR between diabetic and nondiabetic relatives (P=0.16), and the genetic correlation (rG=0.64) for ACR between these two groups was not different from 1 (P=0.12). These results suggested that similar genes contribute to variation in ACR in diabetic and nondiabetic relatives. This hypothesis was supported further by the linkage results.