Genotyping for genetic association studies: Methods and applications

In this thesis, two separate single nucleotide polymorphism (SNP) genotyping techniques were set up at the Finnish Genome Center, pooled genotyping was evaluated as a screening method for large-scale association studies, and finally, the former approaches were used to identify genetic factors predisposing to two distinct complex diseases by utilizing large epidemiological cohorts and also taking environmental factors into account. The first genotyping platform was based on traditional but improved restriction-fragment-length-polymorphism (RFLP) utilizing 384-microtiter well plates, multiplexing, small reaction volumes (5 µl), and automated genotype calling. We participated in the development of the second genotyping method, based on single nucleotide primer extension (SNuPeTM by Amersham Biosciences), by carrying out the alpha- and beta tests for the chemistry and the allele-calling software. Both techniques proved to be accurate, reliable, and suitable for projects with thousands of samples and tens of markers. Pooled genotyping (genotyping of pooled instead of individual DNA samples) was evaluated with Sequenom s MassArray MALDI-TOF, in addition to SNuPeTM and PCR-RFLP techniques. We used MassArray mainly as a point of comparison, because it is known to be well suited for pooled genotyping. All three methods were shown to be accurate, the standard deviations between measurements being 0.017 for the MassArray, 0.022 for the PCR-RFLP, and 0.026 for the SNuPeTM. The largest source of error in the process of pooled genotyping was shown to be the volumetric error, i.e., the preparation of pools. We also demonstrated that it would have been possible to narrow down the genetic locus underlying congenital chloride diarrhea (CLD), an autosomal recessive disorder, by using the pooling technique instead of genotyping individual samples. Although the approach seems to be well suited for traditional case-control studies, it is difficult to apply if any kind of stratification based on environmental factors is needed. Therefore we chose to continue with individual genotyping in the following association studies. Samples in the two separate large epidemiological cohorts were genotyped with the PCR-RFLP and SNuPeTM techniques. The first of these association studies concerned various pregnancy complications among 100,000 consecutive pregnancies in Finland, of which we genotyped 2292 patients and controls, in addition to a population sample of 644 blood donors, with 7 polymorphisms in the potentially thrombotic genes. In this thesis, the analysis of a sub-study of pregnancy-related venous thromboses was included. We showed that the impact of factor V Leiden polymorphism on pregnancy-related venous thrombosis, but not the other tested polymorphisms, was fairly large (odds ratio 11.6; 95% CI 3.6-33.6), and increased multiplicatively when combined with other risk factors such as obesity or advanced age. Owing to our study design, we were also able to estimate the risks at the population level. The second epidemiological cohort was the Helsinki Birth Cohort of men and women who were born during 1924-1933 in Helsinki. The aim was to identify genetic factors that might modify the well known link between small birth size and adult metabolic diseases, such as type 2 diabetes and impaired glucose tolerance. Among ~500 individuals with detailed birth measurements and current metabolic profile, we found that an insertion/deletion polymorphism of the angiotensin converting enzyme (ACE) gene was associated with the duration of gestation, and weight and length at birth. Interestingly, the ACE insertion allele was also associated with higher indices of insulin secretion (p=0.0004) in adult life, but only among individuals who were born small (those among the lowest third of birth weight). Likewise, low birth weight was associated with higher indices of insulin secretion (p=0.003), but only among carriers of the ACE insertion allele. The association with birth measurements was also found with a common haplotype of the glucocorticoid receptor (GR) gene. Furthermore, the association between short length at birth and adult impaired glucose tolerance was confined to carriers of this haplotype (p=0.007). These associations exemplify the interaction between environmental factors and genotype, which, possibly due to altered gene expression, predisposes to complex metabolic diseases. Indeed, we showed that the common GR gene haplotype associated with reduced mRNA expression in thymus of three individuals (p=0.0002).

Molecular basis of the kidney filtration barrier : Role of the nephrin protein complex

The kidney filtration barrier consists of fenestrated endothelial cell layer, glomerular basement membrane and slit diaphragm (SD), the specialized junction between glomerular viscelar epithelial cells (podocytes). Podocyte injury is associated with the development of proteinuria, and if not reversed the injury will lead to permanent deterioration of the glomerular filter. The early events are characterized by disruption of the integrity of the SD, but the molecular pathways involved are not fully understood. Congenital nephrotic syndrome of the Finnish type (CNF) is caused by mutations in NPHS1, the gene encoding the SD protein nephrin. Lack of nephrin results in loss of the SD and massive proteinuria beginning before birth. Furthermore, nephrin expression is decreased in acquired human kidney diseases including diabetic nephropathy. This highlights the importance of nephrin and consequently SD in regulating the kidney filtration function. However, the precise molecular mechanism of how nephrin is involved in the formation of the SD is unknown. This thesis work aimed at clarifying the role of nephrin and its interaction partners in the formation of the SD. The purpose was to identify novel proteins that associate with nephrin in order to define the essential molecular complex required for the establishment of the SD. The aim was also to decipher the role of novel nephrin interacting proteins in podocytes. Nephrin binds to nephrin-like proteins Neph1 and Neph2, and to adherens junction protein P-cadherin. These interactions have been suggested to play a role in the formation of the SD. In this thesis work, we identified densin as a novel interaction partner for nephrin. Densin was localized to the SD and it was shown to bind to adherens junction protein beta-catenin. Furthermore, densin was shown to behave in a similar fashion as adherens junction proteins in cell-cell contacts. These results indicate that densin may play a role in cell adhesion and, therefore, may contribute to the formation of the SD together with nephrin and adherens junction proteins. Nephrin was also shown to bind to Neph3, which has been previously localized to the SD. Neph3 and Neph1 were shown to induce cell adhesion alone, whereas nephrin needed to trans-interact with Neph1 or Neph3 from the opposite cell surface in order to make cell-cell contacts. This was associated with the decreased tyrosine phosphorylation of nephrin. These data extend the current knowledge of the molecular composition of the nephrin protein complex at the SD and also provide novel insights of how the SD may be formed. This thesis work also showed that densin was up-regulated in the podocytes of CNF patients. Neph3 was up-regulated in nephrin deficient mouse kidneys, which share similar podocyte alterations and lack of the SD as observed in CNF patients podocytes. These data suggest that densin and Neph3 may have a role in the formation of morphological alterations in podocytes detected in CNF patients. Furthermore, this thesis work showed that deletion of beta-catenin specifically from adult mouse podocytes protected the mice from the development of adriamycin-induced podocyte injury and proteinuria compared to wild-type mice. These results show that beta-catenin play a role in the adriamycin induced podocyte injury. Podocyte injury is a hallmark in many kidney diseases and the changes observed in the podocytes of CNF patient share characteristics with injured podocytes observed in chronic kidney diseases. Therefore, the results obtained in this thesis work suggest that densin, Neph3 and beta-catenin participate in the molecular pathways which result in morphological alterations commonly detected in injured podocytes in kidney diseases.

Molecular Genetic Background of Tumours in Hereditary Leiomyomatosis and Renal Cell Cancer Syndrome

Hereditary Leiomyomatosis and Renal Cell Cancer (HLRCC) is a hereditary tumour predisposition syndrome. Its phenotype includes benign cutaneous and uterine leiomyomas (CLM, ULM) with high penetrance and rarer renal cell cancer (RCC), most commonly of papillary type 2 subtype. Over 130 HLRCC families have been identified world-wide but the RCC phenotype seems to concentrate in families from Finland and North America for unknown reasons. HLRCC is caused by heterozygous germline mutations in the fumarate hydratase (FH) gene. FH encodes the enzyme fumarase from mitochondrial citric acid cycle. Fumarase enzyme activity or type or site of the FH mutation are unassociated with disease phenotype. The strongest evidence for tumourigenesis mechanism in HLRCC supports a hypoxia inducible factor driven process called pseudohypoxia resulting from accumulation of the fumarase substrate fumarate. In this study, to assess the importance of gene- or exon-level deletions or amplifications of FH in patients with HLRCC-associated phenotypes, multiplex ligation-dependent probe amplification (MLPA) method was used. One novel FH mutation, deletion of exon 1, was found in a Swedish male patient with an evident HLRCC phenotype with CLM, RCC, and a family history of ULM and RCC. Six other patients with CLM and 12 patients with only RCC or uterine leiomyosarcoma (ULMS) remained FH mutation-negative. These results suggest that copy number aberrations of FH or its exons are an infrequent cause of HLRCC and that only co-occurrence of benign tumour types justifies FH-mutation screening in RCC or ULMS patients. Determination of the genomic profile of 11 HLRCC-associated RCCs from Finnish patients was performed by array comparative genomic hybridization. The most common copy number aberrations were gains of 2, 7, and 17 and losses of 13q12.3-q21.1, 14, 18, and X. When compared to aberrations of sporadic papillary RCCs, HLRCC-associated RCCs harboured a distinct DNA copy number profile and lacked many of the changes characterizing the sporadic RCCs. The findings suggest a divergent molecular pathway for tumourigenesis of papillary RCCs in HLRCC. In order to find a genetic modifier of RCC risk in HLRCC, genome-wide linkage and identical by descent (IBD) analysis studies were performed in Finnish HLRCC families with microsatellite marker mapping and SNP-array platforms. The linkage analysis identified only one locus of interest, the FH gene locus in 1q43, but no mutations were found in the genes of the region. IBD analysis yielded no convincing haplotypes shared by RCC patients. Although these results do not exclude the existence of a genetic modifier for RCC risk in HLRCC, they emphasize the role of FH mutations in the malignant tumourigenesis of HLRCC. To study the benign tumours in HLRCC, genome-wide DNA copy number and gene expression profiles of sporadic and HLRCC ULMs were defined with modern SNP- and gene-expression array platforms. The gene expression array suggests novel genes involved in FH-deficient ULM tumourigenesis and novel genes with putative roles in propagation of sporadic ULM. Both the gene expression and copy number profiles of HLRCC ULMs differed from those of sporadic ULMs indicating distinct molecular basis of the FH-deficient HLRCC tumours.

Genomic microarrays in chromosomal analysis of leukemia

Chromosomal alterations in leukemia have been shown to have prognostic and predictive significance and are also important minimal residual disease (MRD) markers in the follow-up of leukemia patients. Although specific oncogenes and tumor suppressors have been discovered in some of the chromosomal alterations, the role and target genes of many alterations in leukemia remain unknown. In addition, a number of leukemia patients have a normal karyotype by standard cytogenetics, but have variability in clinical course and are often molecularly heterogeneous. Cytogenetic methods traditionally used in leukemia analysis and diagnostics; G-banding, various fluorescence in situ hybridization (FISH) techniques, and chromosomal comparative genomic hybridization (cCGH), have enormously increased knowledge about the leukemia genome, but have limitations in resolution or in genomic coverage. In the last decade, the development of microarray comparative genomic hybridization (array-CGH, aCGH) for DNA copy number analysis and the SNP microarray (SNP-array) method for simultaneous copy number and loss of heterozygosity (LOH) analysis has enabled investigation of chromosomal and gene alterations genome-wide with high resolution and high throughput. In these studies, genetic alterations were analyzed in acute myeloid leukemia (AML) and chronic lymphocytic leukemia (CLL). The aim was to screen and characterize genomic alterations that could play role in leukemia pathogenesis by using aCGH and SNP-arrays. One of the most important goals was to screen cryptic alterations in karyotypically normal leukemia patients. In addition, chromosomal changes were evaluated to narrow the target regions, to find new markers, and to obtain tumor suppressor and oncogene candidates. The work presented here shows the capability of aCGH to detect submicroscopic copy number alterations in leukemia, with information about breakpoints and genes involved in the alterations, and that genome-wide microarray analyses with aCGH and SNP-array are advantageous methods in the research and diagnosis of leukemia. The most important findings were the cryptic changes detected with aCGH in karyotypically normal AML and CLL, characterization of amplified genes in 11q marker chromosomes, detection of deletion-based mechanisms of MLL-ARHGEF12 fusion gene formation, and detection of LOH without copy number alteration in karyotypically normal AML. These alterations harbor candidate oncogenes and tumor suppressors for further studies.

Fumarate Hydratase and Succinate Dehydrogenase in Neoplasia

Germline mutations in fumarate hydratase (FH) cause hereditary leiomyomatosis and renal cell cancer (HLRCC). FH is a nuclear encoded enzyme which functions in the Krebs tricarboxylic acid cycle, and homozygous mutation in FH lead to severe developmental defects. Both uterine and cutaneous leiomyomas are components of the HLRCC phenotype. Most of these tumours show loss of the wild-type allele and, also, the mutations reduce FH enzyme activity, which indicate that FH is a tumour suppressor gene. The renal cell cancers associated with HLRCC are of rare papillary type 2 histology. Other genes involved in the Krebs cycle, which are also implicated in neoplasia are 3 of the 4 subunits encoding succinate dehydrogenase (SDH); mutations in SHDB, SDHC, and SDHD predispose to paraganglioma and phaeochromocytoma. Although uterine leiomyomas (or fibroids) are very common, the estimations of affected women ranging from 25% to 77%, not much is known about their genetic background. Cytogenetic studies have revealed that rearrangements involving chromosomes 6, 7, 12 and 14 are most commonly seen in fibroids. Deletions on the long arm of chromosome 7 have been reported to be involved in about 17 to 34 % of leiomyomas and the small commonly deleted region on 7q22 suggests that there might be an underlying tumour suppressor gene in that region. The purpose of this study was to investigate the genetic mechanisms behind the development of tumours associated with HLRCC, both renal cell cancer and uterine fibroids. Firstly, a database search at the Finnish cancer registry was conducted in order to identify new families with early-onset RCC and to test if the family history was compatible with HLRCC. Secondly, sporadic uterine fibroids were tested for deletions on 7q in order to define the minimal deleted 7q-region, followed by mutation analysis of the candidate genes. Thirdly, oligonucleotide chips were utilised to study the global gene expression profiles of uterine fibroids in order to test whether 7q-deletions and FH mutations significantly affected fibroid biology. In the screen for early-onset RCC, 214 families were identified. Subsequently, the pedigrees were constructed and clinical data obtained. One of the index cases (RCC at the age of 28) had a mother who had been diagnosed with a heart tumour, which in further investigation turned out to be a paraganglioma. This lead to an alternative hypothesis that SDH, instead of FH, could be involved. SDHA, SDHB, SDHC and SDHD were sequenced from these individuals; a germline SDHB R27X mutation was detected with loss of the wild-type allele in both tumours. These results suggest that germline mutations in the SDHB gene predispose to early-onset RCC establishing a novel form of hereditary RCC. This has immediate clinical implications in the surveillance of patients suffering from early-onset RCC and phaeochromocytoma/paraganglioma. For the studies on sporadic uterine fibroids, a set of 166 fibroids from 51 individuals were collected. The 7q LOH mapping defined a commonly deleted region of about 3.2 mega bases in 11 of the 166 tumours. The deletion was consistent with previously reported allelotyping studies of leiomyomas and it therefore suggested the presence of a tumour suppressor gene in the deleted region. Furthermore, the high-resolution aCGH-chip analysis refined the deleted region to only 2.79Mb. When combined with previous data, the commonly deleted region was only 2.3Mb. The mutation screening of the known genes within the commonly deleted region did not reveal pathogenic mutations, however. The expression microarray analysis revealed that FH-deficient fibroids, both sporadic and familial, had their distinct gene expression profile as they formed their own group in the unsupervised clustering. On the other hand, the presence or absence of 7q-deletions did not significantly alter the global gene expression pattern of fibroids, suggesting that these two groups do not have different biological backgrounds. Multiple differentially expressed genes were identified between FH wild-type and FH-mutant fibroids, and the most significant increase was seen in the expression of carbohydrate metabolism-related and hypoxia inducible factor (HIF) target genes.

The Myopathic Protein Myotilin in Developing Mouse and in Muscle Function

Skeletal muscle cells are highly specialised in order to accomplish their function. During development, the fusion of hundreds of immature myoblasts creates large syncytial myofibres with a highly ordered cytoplasm filled with packed myofibrils. The assembly and organisation of contractile myofibrils must be tightly controlled. Indeed, the number of proteins involved in sarcomere building is impressive, and the role of many of them has only recently begun to be elucidated. Myotilin was originally identified as a high affinity a-actinin binding protein in yeast twohybrid screen. It was then found to interact also with filamin C, actin, ZASP and FATZ-1. Human myotilin is mainly expressed in striated muscle and induces efficient actin bundling in vitro and in cells. Moreover, mutations in myotilin cause different forms of muscle disease, now collectively known as myotilinopathies. In this thesis, consisting of three publications, the work on the mouse orthologue is presented. First, the cloning and molecular characterisation of the mouse myotilin gene showed that human and mouse myotilin share high sequence homology and a similar expression pattern and gene regulation. Functional analysis of the mouse promoter revealed the myogenic factor-binding elements that are required for myotilin gene transcription. Secondly, expression of myotilin was studied during mouse embryogenesis. Surprisingly, myotilin was expressed in a wide array of tissues at some stages of development; its expression pattern became more restricted at perinatal stages and in adult life. Immunostaining of human embryos confirmed broader myotilin expression compared to the sarcomeric marker titin. Finally, in the third article, targeted deletion of myotilin gene in mice revealed that it is not essential for muscle development and function. These data altogether indicate that the mouse can be used as a model for human myotilinopathy and that loss of myotilin does not alter significantly muscle structure and function. Therefore, disease-associated mutant myotilin may act as a dominant myopathic factor.

Genetic basis of hereditary colorectal cancers: Hereditary nonpolyposis colorectal cancer and Familial adenomatous polyposis

Hereditary nonpolyposis colorectal cancer (HNPCC) and familial adenomatous polyposis (FAP) are characterized by a high risk and early onset of colorectal cancer (CRC). HNPCC is due to a germline mutation in one of the following MMR genes: MLH1, MSH2, MSH6 and PMS2. A majority of FAP and attenuated FAP (AFAP) cases are due to germline mutations of APC, causing the development of multiple colorectal polyps. To date, over 450 MMR gene mutations and over 800 APC mutations have been identified. Most of these mutations lead to a truncated protein, easily detected by conventional mutation detection methods. However, in about 30% of HNPCC and FAP, and about 90% of AFAP families, mutations remain unknown. We aimed to clarify the genetic basis and genotype-phenotype correlation of mutation negative HNPCC and FAP/AFAP families by advanced mutation detection methods designed to detect large genomic rearrangements, mRNA and protein expression alterations, promoter mutations, phenotype linked haplotypes, and tumoral loss of heterozygosity. We also aimed to estimate the frequency of HNPCC in Uruguayan CRC patients. Our expression based analysis of mutation negative HNPCC divided these families into two categories: 1) 42% of families linked to the MMR genes with a phenotype resembling that of mutation positive, and 2) 58% of families likely to be associated with other susceptibility genes. Unbalanced mRNA expression of MLH1 was observed in two families. Further studies revealed that a MLH1 nonsense mutation, R100X was associated with aberrant splicing of exons not related to the mutation and an MLH1 deletion (AGAA) at nucleotide 210 was associated with multiple exon skipping, without an overall increase in the frequency of splice events. APC mutation negative FAP/AFAP families were divided into four groups according to the genetic basis of their predisposition. Four (14%) families displayed a constitutional deletion of APC with profuse polyposis, early age of onset and frequent extracolonic manifestations. Aberrant mRNA expression of one allele was observed in seven (24%) families with later onset and less frequent extracolonic manifestations. In 15 (52%) families the involvement of APC could neither be confirmed nor excluded. In three (10%) of the families a germline mutation was detected in genes other than APC: AXIN2 in one family, and MYH in two families. The families with undefined genetic basis and especially those with AXIN2 or MYH mutations frequently displayed AFAP or atypical polyposis. Of the Uruguayan CRC patients, 2.6% (12/461) fulfilled the diagnostic criteria for HNPCC and 5.6% (26/461) were associated with increased risk of cancer. Unexpectedly low frequency of molecularly defined HNPCC cases may suggest a different genetic profile in the Uruguayan population and the involvement of novel susceptibility genes. Accurate genetic and clinical characterization of families with hereditary colorectal cancers, and the definition of the genetic basis of "mutation negative" families in particular, facilitate proper clinical management of such families.

Microarrays in Molecular Profiling of Ewing Sarcoma Family of Tumors and CDKN2A Aberrations

Background: The Ewing sarcoma family of tumors (ESFT) are rare but highly malignant neoplasms that occur mainly in bone or but also in soft tissue. ESFT affects patients typically in their second decade of life, whereby children and adolescents bear the heaviest incidence burden. Despite recent advances in the clinical management of ESFT patients, their prognosis and survival are still disappointingly poor, especially in cases with metastasis. No targeted therapy for ESFT patients is currently available. Moreover, based merely on current clinical and biological characteristics, accurate classification of ESFT patients often fails at the time of diagnosis. Therefore, there is a constant need for novel molecular biomarkers to be applied in tandem with conventional parameters to further intensify ESFT risk-stratification and treatment selection, and ultimately to develop novel targeted therapies. In this context, a greater understanding of the genetics and immune characteristics of ESFT is needed. Aims: This study sought to open novel insights into gene copy number changes and gene expression in ESFT and, further, to enlighten the role of inflammation in ESFT. For this purpose, microarrays were used to provide gene-level information on a genomewide scale. In addition, this study focused on screening of 9p21.3 deletion sizes and frequencies in ESFT and, in another pediatric cancer, acute lymphocytic leukemia (ALL), in order to define more exact criteria for highrisk patient selection and to provide data for developing a more reliable diagnostic method to detect CDKN2A deletions. Results: In study I, 20 novel ESFT-associated suppressor genes and oncogenes were pinpointed using combined array CGH and expression analysis. In addition, interesting chromosomal rearrangements were identified: (1) Duplication of derivative chromosome der(22)(11;22) was detected in three ESFT patients. This duplication included the EWSR1-FLI1 fusion gene leading to increase in its copy number; (2) Cryptic amplifications on chromosomes 20 and 22 were detected, suggesting a novel translocation between chromosomes 20 and 22, which most probably produces a fusion between EWSR1 and NFATC2. In study II, bioinformatic analysis of ESFT expression profiles showed that inflammatory gene activation is detectable in ESFT patient samples and that the activation is characterized by macrophage gene expression. Most interestingly, ESFT patient samples were shown to express certain inflammatory genes that were prognostically significant. High local expression of C5 and JAK1 at the tumor site was shown to associate with favorable clinical outcome, whereas high local expression of IL8 was shown to be detrimental. Studies III and IV showed that the smallest overlapping region of deletion in 9p21.3 includes CDKN2A in all cases and that the length of this region is 12.2 kb in both Ewing sarcoma and ALL. Furthermore, our results showed that the most widely used commercial CDKN2A FISH probe creates false negative results in the narrowest microdeletion cases (<190 kb). Therefore, more accurate methods should be developed for the detection of deletions in the CDKN2A locus. Conclusions: This study provides novel insights into the genetic changes involved in the biology of ESFT, in the interaction between ESFT cells and immune system, and in the inactivation of CDKN2A. Novel ESFT biomarker genes identified in this study serve as a useful resource for future studies and in developing novel therapeutic strategies to improve the survival of patients with ESFT.