Insights into the molecular genetics of celiac disease : Applying family- and population-based strategies

Celiac disease, or gluten intolerance, is triggered by dietary glutens in genetically susceptible individuals and it affects approximately 1% of the Caucasian population. The best known genetic risk factors for celiac disease are HLA DQ2 and DQ8 heterodimers, which are necessary for the development of the disease. However, they alone are not sufficient for disease induction, other risk factors are required. This thesis investigated genetic factors for celiac disease, concentrating on susceptibility loci on chromosomes 5q31-q33, 19p13 and 2q12 previously reported in genome-wide linkage and association studies. In addition, a novel genotyping method for the detection of HLA DQ2 and DQ8 coding haplotypes was validated. This study was conducted using Finnish and Hungarian family materials, and Finnish, Hungarian and Italian case-control materials. Genetic linkage and association were analysed in these materials using candidate gene and fine-mapping approaches. The results confirmed linkage to celiac disease on the chromosomal regions 5q31-q33 and 19p13. Fine-mapping on chromosome 5q31-q33 revealed several modest associations in the region, and highlighted the need for further investigations to locate the causal risk variants. The MYO9B gene on chromosome 19p13 showed evidence for linkage and association particularly with dermatitis herpetiformis, the skin manifestation of celiac disease. This implies a potential difference in the genetic background of the intestinal and skin forms of the disease, although studies on larger samplesets are required. The IL18RAP locus on chromosome 2q12, shown to be associated with celiac disease in a previous genome-wide association study and a subsequent follow-up, showed association in the Hungarian population in this study. The expression of IL18RAP was further investigated in small intestinal tissue and in peripheral blood mononuclear cells. The results showed that IL18RAP is expressed in the relevant tissues. Two putative isoforms of IL18RAP were detected by Western blot analysis, and the results suggested that the ratios and total levels of these isoforms may contribute to the aetiology of celiac disease. A novel genotyping method for celiac disease-associated HLA haplotypes was also validated in this thesis. The method utilises single-nucleotide polymorphisms tagging these HLA haplotypes with high sensitivity and specificity. Our results suggest that this method is transferable between populations, and it is suitable for large-scale analysis. In conclusion, this doctorate study provides an insight into the roles of the 5q31-q33, MYO9B, IL18RAP and HLA loci in the susceptibility to celiac disease in the Finnish, Hungarian and Italian populations, highlighting the need for further studies at these genetic loci and examination of the function of the candidate genes.

Interleukin 9: A candidate gene for asthma

Asthma is a complex heritable inflammatory disorder of the airways associated with clinical signs of atopy and bronchial hyperresponsiveness. Recent studies localized a major gene for asthma to chromosome 5q31-q33 in humans. Thus, this segment of the genome represents a candidate region for genes that determine susceptibility to bronchial hyperresponsiveness and atopy in animal models. Homologs of candidate genes on human chromosome 5q31-q33 are found in four regions in the mouse genome, two on chromosome 18, and one each on chromosomes 11 and 13. We assessed bronchial responsiveness as a quantitative trait in mice and found it linked to chromosome 13. Interleukin 9 (IL-9) is located in the linked region and was analyzed as a gene candidate. The expression of IL-9 was markedly reduced in bronchial hyporesponsive mice, and the level of expression was determined by sequences within the qualitative trait locus (QTL). These data suggest a role for IL-9 in the complex pathogenesis of bronchial hyperresponsiveness as a risk factor for asthma.

Major quantitative trait locus for eosinophil count is located on chromosome 2q

Background: Eosinophils are granulocytic white blood cells implicated in asthma and atopic disease. The degree of eosinophilia in the blood of patients with asthma correlates with the severity of asthmatic symptoms. Quantitative trait loci (QTL) linkage analysis of eosinophil count may be a more powerful strategy of mapping genes involved in asthma than linkage analysis using affected relative pairs. 1 Objective: To identify QTLs responsible for variation in eosinophil count in adolescent twins. Methods: We measured eosinophil count longitudinally in 738 pairs of twins at 12, 14, and 16 years of age. We typed 757 highly polymorphic microsatellite markers at an average spacing of similar to5 centimorgans across the genome. We then used multipoint variance components linkage analysis to test for linkage between marker loci and eosinophil concentrations at each age across the genome. Results: We found highly significant linkage on chromosome 2q33 in 12-year-old twins (logarithm of the odds = 4.6; P = .000002) and suggestive evidence of linkage in the same region in 14-year-olds (logarithm of the odds = 1.0; P = .016). We also found suggestive evidence of linkage at other areas of the genome, including regions on chromosomes 2, 3, 4, 8, 9, 11, 12, 17, 20, and 22. Conclusion: A QTL for eosinophil count is present on chromosome 2q33. This QTL might represent a gene involved in asthma pathophysiology.

Familial aggregation of human infection with Schistosoma japonicum in the Poyang Lake region, China

Despite the success of extensive control measures that have been implemented in China for over 50 years, the number of individuals infected with Schistosoma japonicum remains high in the existing endemic areas. A variance components analysis was undertaken to estimate the heritable and environmental components that contribute to S. japonicum infection in the Poyang Lake region of Jiangxi Province, PR China. The total target population was 3148 from four separate administrative villages. Two thousand seven hundred and five of these comprised 400 families ranging in size from 3 to 188. After adjustments were made for gender, water contact and past history of having had schistosomiasis, the heritable component was estimated to account for as much as 58% of the phenotype variation under the polygenic model. Household was not shown to be an important environmental factor. Incorporating village effects indicated that the results were valid for the total population. We conclude that genetic heritability in this region is high and plays an important role in determining risk of infection with S. japonicum. (c) 2005 Australian Society for Parasitology Inc. Published by Elsevier Ltd. All rights reserved.

Linkage to chromosome 2q32.2-q33.3 in familial serrated neoplasia (Jass syndrome)

Abstract Causative genetic variants have to date been identified for only a small proportion of familial colorectal cancer (CRC). While conditions such as Familial Adenomatous Polyposis and Lynch syndrome have well defined genetic causes, the search for variants underlying the remainder of familial CRC is plagued by genetic heterogeneity. The recent identification of families with a heritable predisposition to malignancies arising through the serrated pathway (familial serrated neoplasia or Jass syndrome) provides an opportunity to study a subset of familial CRC in which heterogeneity may be greatly reduced. A genome-wide linkage screen was performed on a large family displaying a dominantly-inherited predisposition to serrated neoplasia genotyped using the Affymetrix GeneChip Human Mapping 10 K SNP Array. Parametric and nonparametric analyses were performed and resulting regions of interest, as well as previously reported CRC susceptibility loci at 3q22, 7q31 and 9q22, were followed up by finemapping in 10 serrated neoplasia families. Genome-wide linkage analysis revealed regions of interest at 2p25.2-p25.1, 2q24.3-q37.1 and 8p21.2-q12.1. Finemapping linkage and haplotype analyses identified 2q32.2-q33.3 as the region most likely to harbour linkage, with heterogeneity logarithm of the odds (HLOD) 2.09 and nonparametric linkage (NPL) score 2.36 (P = 0.004). Five primary candidate genes (CFLAR, CASP10, CASP8, FZD7 and BMPR2) were sequenced and no segregating variants identified. There was no evidence of linkage to previously reported loci on chromosomes 3, 7 and 9.

Partial epilepsy syndrome in a Gypsy family linked to 5q31.3-q32

PURPOSE The restricted genetic diversity and homogeneous molecular basis of Mendelian disorders in isolated founder populations have rarely been explored in epilepsy research. Our long-term goal is to explore the genetic basis of epilepsies in one such population, the Gypsies. The aim of this report is the clinical and genetic characterization of a Gypsy family with a partial epilepsy syndrome. METHODS Clinical information was collected using semistructured interviews with affected subjects and informants. At least one interictal electroencephalography (EEG) recording was performed for each patient and previous data obtained from records. Neuroimaging included structural magnetic resonance imaging (MRI). Linkage and haplotype analysis was performed using the Illumina IVb Linkage Panel, supplemented with highly informative microsatellites in linked regions and Affymetrix SNP 5.0 array data. RESULTS We observed an early-onset partial epilepsy syndrome with seizure semiology strongly suggestive of temporal lobe epilepsy (TLE), with mild intellectual deficit co-occurring in a large proportion of the patients. Psychiatric morbidity was common in the extended pedigree but did not cosegregate with epilepsy. Linkage analysis definitively excluded previously reported loci, and identified a novel locus on 5q31.3-q32 with an logarithm of the odds (LOD) score of 3 corresponding to the expected maximum in this family. DISCUSSION The syndrome can be classified as familial temporal lobe epilepsy (FTLE) or possibly a new syndrome with mild intellectual deficit. The linked 5q region does not contain any ion channel-encoding genes and is thus likely to contribute new knowledge about epilepsy pathogenesis. Identification of the mutation in this family and in additional patients will define the full phenotypic spectrum.

Association between the 5q31.1 gene neurogenin1 and schizophrenia

Multiple lines of evidence suggest that schizophrenia results from aberrant neurodevelopment. The neurogenin1 gene (neurog1) consists of a single 1,666 bp exon that encodes a basic helix-loop-helix (bHLH) transcription factor that causes neuronal differentiation and induces cortical and glutamatergic differentiation programs. Because of its function and its location in 5q31.1, which has been linked to schizophrenia in multiple samples, we tested it for association with the disorder. We sequenced neurog1 in 25 affected subjects from the Irish Study of High-Density Schizophrenia Families. We observed a 5'-UTR SNP at position -60, already present in databases as rs8192558, and tested it along with rs2344485, rs8192559, and rs2344484. Narrow, intermediate, and broad diagnostic definitions were used. The major alleles of rs8192558 and rs2344484 were over-transmitted to affected subjects using both Pedigree Disequilibrium Test (PDT) (0.01 <or = P <or = 0.06) and FBAT (0.02 <or = P <or = 0.07). A haplotype consisting of the major alleles of all four SNPs was significantly over-transmitted in FBAT to the broad definition (P = 0.049), with trend significance to the narrow and intermediate definitions, and with trend significance in PDT. In confirmatory tests using 657 cases and 411 controls, this haplotype was slightly but not significantly over-represented in cases (81% vs. 77%, P = 0.21). These results, along with a priori evidence for the involvement of neurog1 in neurodevelopment, suggest that variants in neurog1 might have a small effect on susceptibility to schizophrenia. This gene should be tested in additional and larger samples.

Bayesian latent trait modeling of migraine symptom data

Definition of disease phenotype is a necessary preliminary to research into genetic causes of a complex disease. Clinical diagnosis of migraine is currently based on diagnostic criteria developed by the International Headache Society. Previously, we examined the natural clustering of these diagnostic symptoms using latent class analysis (LCA) and found that a four-class model was preferred. However, the classes can be ordered such that all symptoms progressively intensify, suggesting that a single continuous variable representing disease severity may provide a better model. Here, we compare two models: item response theory and LCA, each constructed within a Bayesian context. A deviance information criterion is used to assess model fit. We phenotyped our population sample using these models, estimated heritability and conducted genome-wide linkage analysis using Merlin-qtl. LCA with four classes was again preferred. After transformation, phenotypic trait values derived from both models are highly correlated (correlation = 0.99) and consequently results from subsequent genetic analyses were similar. Heritability was estimated at 0.37, while multipoint linkage analysis produced genome-wide significant linkage to chromosome 7q31-q33 and suggestive linkage to chromosomes 1 and 2. We argue that such continuous measures are a powerful tool for identifying genes contributing to migraine susceptibility.

Co-movements in commodity prices : a note based on network analysis

This article analyses co-movements in a wide group of commodity prices during the time period 1992–2010. Our methodological approach is based on the correlation matrix and the networks inside. Through this approach we are able to summarize global interaction and interdependence, capturing the existing heterogeneity in the degrees of synchronization between commodity prices. Our results produce two main findings: (a) we do not observe a persistent increase in the degree of co-movement of the commodity prices in our time sample, however from mid-2008 to the end of 2009 co-movements almost doubled when compared with the average correlation; (b) we observe three groups of commodities which have exhibited similar price dynamics (metals, oil and grains, and oilseeds) and which have increased their degree of co-movement during the sampled period.

The molecular basis of Marinesco-Sjögren syndrome

Marinesco-Sjögren syndrome (MSS) is a rare autosomal recessive neurodegenerative disorder characterized by cerebellar ataxia due to cerebellar cortical atrophy, infantile- or childhood-onset bilateral cataracts, progressive myopathy, and mild to severe mental retardation. Additional features include hypergonadotropic hypogonadism, various skeletal abnormalities, short stature, and strabismus. The neuroradiologic hallmarks are hypoplasia of both the vermis and cerebellar hemispheres. The histopathologic findings include severe cerebellar atrophy and loss of Purkinje and granule cells. The common pathologic findings in muscle biopsy are variation in muscle fiber size, atrophic fibers, fatty replacement, and rimmed vacuole formation. The presence of marked cerebellar atrophy with myopathy distinguishes MSS from another rare syndrome, the congenital cataracts, facial dysmorphism, and neuropathy syndrome (CCFDN). Previously, work by others had resulted in the identification of an MSS locus on chromosome 5q31. A subtype of MSS with myoglobinuria and neuropathy had been linked to the CCFDN locus on chromosome 18qter, at which mutations in the CTDP1 gene had been identified. We confirmed linkage to the previously identified locus on chromosome 5q31 in two Finnish families with eight affected individuals, reduced the critical region by fine-mapping, and identified SIL1 as a gene underlying MSS. We found a common homozygous founder mutation in all Finnish patients. The same mutation was also present in patient samples from Norway and Sweden. Altogether, we identified eight mutations in SIL1, including nonsense, frameshift, splice site alterations, and one missense mutation. SIL1 encodes a nucleotide exchange factor for the endoplasmic reticulum (ER) resident heat-shock protein 70 chaperone GRP78. GRP78 functions in protein synthesis and quality control of the newly synthesized polypeptides. It senses and responds to stressful cellular conditions. We showed that in mice, SIL1 and GRP78 show highly similar spatial and temporal tissue expression in developing and mature brain, eye, and muscle. Studying endogenous proteins in mouse primary hippocampal neurons, we found that SIL1 and GRP78 colocalize and that SIL1 localizes to the ER. We studied the subcellular localization of two mutant proteins, a missense mutant found in two patients and an artificial mutant lacking the ER retrieval signal, and found that both mutant proteins formed aggregates within the ER. Well in line with our findings and the clinical features of MSS, recent work by Zhao et al. showed that a truncation of SIL1 causes ataxia and cerebellar Purkinje cell loss in the naturally occurring woozy mutant mouse. Prior to Purkinje cell degeneration, the unfolded protein response is initiated and abnormal protein accumulations are present. MSS thus joins the group of protein misfolding and accumulation diseases. These findings highlight the importance of SIL1 and the role of the ER in neuronal function and survival. The results presented in this thesis provide tools for the molecular genetic diagnostics of MSS and give a basis for future studies on the molecular pathogenesis of MSS. Understanding the mechanisms behind this pleiotropic syndrome may provide insights into more common forms of ataxia, myopathy, and neurodegeneration.

Is the histidine triad nucleotide-binding protein 1 (HINT1) gene a candidate for schizophrenia?

Background: The histidine triad nucleotide-binding protein 1, HINT1, hydrolyzes adenosine 5'monophosphoramidate substrates such as AMP-morpholidate. The human HINT1 gene is located on chromosome 5q31.2, a region implicated in linkage studies of schizophrenia. HINT1 had been shown to have different expression in postmortem brains between schizophrenia patients and unaffected controls. It was also found to be associated with the dysregulation of postsynaptic dopamine transmission, thus suggesting a potential role in several neuropsychiatric diseases.

Downregulation of the Wnt inhibitor CXXC5 predicts a better prognosis in acute myeloid leukemia

The gene CXXC5 on 5q31 is frequently deleted in acute myeloid leukemia (AML) with del(5q), suggesting that inactivation of CXXC5 might play a role in leukemogenesis. Here, we investigated the functional and prognostic implications of CXXC5 expression in AML. CXXC5 mRNA was downregulated in AML with MLL rearrangements, t(8;21) and GATA2 mutations. As a mechanism of CXXC5 inactivation, we found evidence for epigenetic silencing by promoter methylation. Patients with CXXC5 expression below the median level had a lower relapse rate (45% vs 59%; P = .007) and a better overall survival (OS, 46% vs 28%; P < .001) and event-free survival (EFS, 36% vs 21%; P < .001) at 5 years, independent of cytogenetic risk groups and known molecular risk factors. In gene-expression profiling, lower CXXC5 expression was associated with upregulation of cell-cycling genes and codownregulation of genes implicated in leukemogenesis (WT1, GATA2, MLL, DNMT3B, RUNX1). Functional analyses demonstrated CXXC5 to inhibit leukemic cell proliferation and Wnt signaling and to affect the p53-dependent DNA damage response. In conclusion, our data suggest a tumor suppressor function of CXXC5 in AML. Inactivation of CXXC5 is associated with different leukemic pathways and defines an AML subgroup with better outcome.

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.

Recurrent genomic alterations in sequential progressive leukoplakia and oral cancer: drivers of oral tumorigenesis?

A significant proportion (up to 62) of oral squamous cell carcinomas (OSCCs) may arise from oral potential malignant lesions (OPMLs), such as leukoplakia. Patient outcomes may thus be improved through detection of lesions at a risk for malignant transformation, by identifying and categorizing genetic changes in sequential, progressive OPMLs. We conducted array comparative genomic hybridization analysis of 25 sequential, progressive OPMLs and same-site OSCCs from five patients. Recurrent DNA copy number gains were identified on 1p in 20/25 cases (80) with minimal, high-level amplification regions on 1p35 and 1p36. Other regions of gains were frequently observed: 11q13.4 (68), 9q34.13 (64), 21q22.3 (60), 6p21 and 6q25 (56) and 10q24, 19q13.2, 22q12, 5q31.2, 7p13, 10q24 and 14q22 (48). DNA losses were observed in 20 of samples and mainly detected on 5q31.2 (35), 16p13.2 (30), 9q33.1 and 9q33.29 (25) and 17q11.2, 3p26.2, 18q21.1, 4q34.1 and 8p23.2 (20). Such copy number alterations (CNAs) were mapped in all grades of dysplasia that progressed, and their corresponding OSCCs, in 70 of patients, indicating that these CNAs may be associated with disease progression. Amplified genes mapping within recurrent CNAs (KHDRBS1, PARP1, RAB1A, HBEGF, PAIP2, BTBD7) were selected for validation, by quantitative real-time PCR, in an independent set of 32 progressive leukoplakia, 32 OSSCs and 21 non-progressive leukoplakia samples. Amplification of BTBD7, KHDRBS1, PARP1 and RAB1A was exclusively detected in progressive leukoplakia and corresponding OSCC. BTBD7, KHDRBS1, PARP1 and RAB1A may be associated with OSCC progression. Proteinprotein interaction networks were created to identify possible pathways associated with OSCC progression.