Computational Identification of Recessive Mutations in Cancers using High Throughput SNP-arrays

This thesis presents a highly sensitive genome wide search method for recessive mutations. The method is suitable for distantly related samples that are divided into phenotype positives and negatives. High throughput genotype arrays are used to identify and compare homozygous regions between the cohorts. The method is demonstrated by comparing colorectal cancer patients against unaffected references. The objective is to find homozygous regions and alleles that are more common in cancer patients. We have designed and implemented software tools to automate the data analysis from genotypes to lists of candidate genes and to their properties. The programs have been designed in respect to a pipeline architecture that allows their integration to other programs such as biological databases and copy number analysis tools. The integration of the tools is crucial as the genome wide analysis of the cohort differences produces many candidate regions not related to the studied phenotype. CohortComparator is a genotype comparison tool that detects homozygous regions and compares their loci and allele constitutions between two sets of samples. The data is visualised in chromosome specific graphs illustrating the homozygous regions and alleles of each sample. The genomic regions that may harbour recessive mutations are emphasised with different colours and a scoring scheme is given for these regions. The detection of homozygous regions, cohort comparisons and result annotations are all subjected to presumptions many of which have been parameterized in our programs. The effect of these parameters and the suitable scope of the methods have been evaluated. Samples with different resolutions can be balanced with the genotype estimates of their haplotypes and they can be used within the same study.

Genomic Profiling of Gastric Cancer

Gastric cancer is the fourth most common cancer and the second most common cause of cancer-related death worldwide. Due to lack of early symptoms, gastric cancer is characterized by late stage diagnosis and unsatisfactory options for curative treatment. Several genomic alterations have been identified in gastric cancer, but the major factors contributing to initiation and progression of gastric cancer remain poorly known. Gene copy number alterations play a key role in the development of gastric cancer, and a change in gene copy number is one of the fundamental mechanisms for a cancer cell to control the expression of potential oncogenes and tumor suppressor genes. This thesis aims at clarifying the complex genomic alterations of gastric cancer to identify novel molecular biomarkers for diagnostic purposes as well as for targeted treatment. To highlight genes of potential biological and clinical relevance, we carried out a systematic microarray-based survey of gene expression and copy number levels in primary gastric tumors and gastric cancer cell lines. Results were validated using immunohistochemistry, real-time qRT-PCR, and affinity capture-based transcript (TRAC) assay. Altogether 192 clinical gastric tissue samples and 7 gastric cancer cell lines were included in this study. Multiple chromosomal regions with recurrent copy number alterations were detected. The most frequent chromosomal alterations included gains at 7q, 8q, 17q, 19q, and 20q and losses at 9p, 18q, and 21q. Distinctive patterns of copy number alterations were detected for different histological subtypes (intestinal and diffuse) and for cancers located in different parts of the stomach. The impact of copy number alterations on gene expression was significant, as 6-10% of genes located in the regions of gains and losses also showed concomitant alterations in their expression. By combining the information from the DNA- and RNA-level analyses many novel gastric cancer-related genes, such as ALPK2, ENAH, HHIPL2, and OSMR, were identified. Independent genome-wide gene expression analysis of Finnish and Japanese gastric tumors revealed an additional set of genes that was differentially expressed in cancerous gastric tissues compared with normal tissue. Overexpression of one of these genes, CXCL1, was associated with an improved survival of gastric cancer. Thus, using an integrative microarray analysis, several novel genes were identified that may be critically important for gastric carcinogenesis. Further studies of these genes may lead to novel biomarkers for gastric cancer diagnosis and targeted therapy.

Molecular Alterations in Asbestos-Related Lung Cancer

Background: Asbestos is a well known cancer-causing mineral fibre, which has a synergistic effect on lung cancer risk in combination with tobacco smoking. Several in vitro and in vivo experiments have demonstrated that asbestos can evoke chromosomal damage and cause alterations as well as gene expression changes. Lung tumours, in general, have very complex karyotypes with several recurrently gained and lost chromosomal regions and this has made it difficult to identify specific molecular changes related primarily to asbestos exposure. The main aim of these studies has been to characterize asbestos-related lung cancer at a molecular level. Methods: Samples from asbestos-exposed and non-exposed lung cancer patients were studied using array comparative genomic hybridization (aCGH) and fluorescent in situ hybridization (FISH) to detect copy number alterations (CNA) as well as microsatellite analysis to detect allelic imbalance (AI). In addition, asbestos-exposed cell lines were studied using gene expression microarrays. Results: Eighteen chromosomal regions showing differential copy number in the lung tumours of asbestos-exposed patients compared to those of non-exposed patients were identified. The most significant differences were detected at 2p21-p16.3, 5q35.3, 9q33.3-q34.11, 9q34.13-q34.3, 11p15.5, 14q11.2 and 19p13.1-p13.3 (p<0.005). The alterations at 2p and 9q were validated and characterized in detail using AI and FISH analysis in a larger study population. Furthermore, in vitro studies were performed to examine the early gene expression changes induced by asbestos in three different lung cell lines. The results revealed specific asbestos-associated gene expression profiles and biological processes as well as chromosomal regions enriched with genes believed to contribute to the common asbestos-related responses in the cell lines. Interestingly, the most significant region enriched with asbestos-response genes was identified at 2p22, close to the previously identified region showing asbestos-related CNA in lung tumours. Additionally, in this thesis, the dysregulated biological processes (Gene Ontology terms) detected in the cell line experiment were compared to dysregulated processes identified in patient samples in a later study (Ruosaari et al., 2008a). Commonly affected processes such as those related to protein ubiquitination, ion transport and surprisingly sensory perception of smell were identified. Conclusions: The identification of specific CNA and dysregulated biological processes shed some light on the underlying genes acting as mediators in asbestos-related lung carcinogenesis. It is postulated that the combination of several asbestos-specific molecular alterations could be used to develop a diagnostic method for the identification of asbestos-related lung cancer.

Acute Lymphoblastic Leukemia in Adolescents and Young Adults in Finland

In recent reports, adolescents and young adults (AYA) with acute lymphoblastic leukemia (ALL) have had a better outcome with pediatric treatment than with adult protocols. ALL can be classified into biologic subgroups according to immunophenotype and cytogenetics, with different clinical characteristics and outcome. The proportions of the subgroups are different in children and adults. ALL subtypes in AYA patients are less well characterized. In this study, the treatment and outcome of ALL in AYA patients aged 10-25 years in Finland on pediatric and adult protocols was retrospectively analyzed. In total, 245 patients were included. The proportions of biologic subgroups in different age groups were determined. Patients with initially normal or failed karyotype were examined with oligonucleotide microarray-based comparative genomic hybridization (aCGH). Also deletions and instability of chromosome 9p were screened in ALL patients. In addition, patients with other hematologic malignancies were screened for 9p instability. aCGH data were also used to determine a gene set that classifies AYA patients at diagnosis according to their risk of relapse. Receiver operating characteristic analysis was used to assess the value of the set of genes as prognostic classifiers. The 5-year event-free survival of AYA patients treated with pediatric or adult protocols was 67% and 60% (p=0.30), respectively. White blood cell count larger than 100x109/l was associated with poor prognosis. Patients treated with pediatric protocols and assigned to an intermediate-risk group fared significantly better than those of the pediatric high-risk or adult treatment groups. Deletions of 9p were detected in 46% of AYA ALL patients. The chromosomal region 9p21.3 was always affected, and the CDKN2A gene was always deleted. In about 15% of AYA patients, the 9p21.3 deletion was smaller than 200 kb in size, and therefore, probably undetectable with conventional methods. Deletion of 9p was the most common aberration of AYA ALL patients with initially normal karyotype. Instability of 9p, defined as multiple separate areas of copy number loss or homozygous loss within a larger heterozygous area in 9p, was detected in 19% (n=27) of ALL patients. This abnormality was restricted to ALL; none of the patients with other hematologic malignancies had the aberration. The prognostic model identification procedure resulted in a model of four genes: BAK1, CDKN2B, GSTM1, and MT1F. The copy number profile combinations of these genes differentiated between AYA ALL patients at diagnosis depending on their risk of relapse. Deletions of CDKN2B and BAK1 in combination with amplification of GSTM1 and MT1F were associated with a higher probability of relapse. Unlike all previous studies, we found that the outcome of AYA patients with ALL treated using pediatric or adult therapeutic protocols was comparable. The success of adult ALL therapy emphasizes the benefit of referral of patients to academic centers and adherence to research protocols. 9p deletions and instability are common features of ALL and may act together with oncogene-activating translocations in leukemogenesis. New and more sensitive methods of molecular cytogenetics can reveal previously cryptic genetic aberrations with an important role in leukemic development and prognosis and that may be potential targets of therapy. aCGH also provides a viable approach for model design aiming at evaluation of risk of relapse in ALL.

Data integration from two microarray platforms identifies bi-allelic genetic inactivation of RIC8A in a breast cancer cell line

Background: Using array comparative genomic hybridization (aCGH), a large number of deleted genomic regions have been identified in human cancers. However, subsequent efforts to identify target genes selected for inactivation in these regions have often been challenging. Methods: We integrated here genome-wide copy number data with gene expression data and non-sense mediated mRNA decay rates in breast cancer cell lines to prioritize gene candidates that are likely to be tumour suppressor genes inactivated by bi-allelic genetic events. The candidates were sequenced to identify potential mutations. Results: This integrated genomic approach led to the identification of RIC8A at 11p15 as a putative candidate target gene for the genomic deletion in the ZR-75-1 breast cancer cell line. We identified a truncating mutation in this cell line, leading to loss of expression and rapid decay of the transcript. We screened 127 breast cancers for RIC8A mutations, but did not find any pathogenic mutations. No promoter hypermethylation in these tumours was detected either. However, analysis of gene expression data from breast tumours identified a small group of aggressive tumours that displayed low levels of RIC8A transcripts. qRT-PCR analysis of 38 breast tumours showed a strong association between low RIC8A expression and the presence of TP53 mutations (P = 0.006). Conclusion: We demonstrate a data integration strategy leading to the identification of RIC8A as a gene undergoing a classical double-hit genetic inactivation in a breast cancer cell line, as well as in vivo evidence of loss of RIC8A expression in a subgroup of aggressive TP53 mutant breast cancers.

Array-based gene expression, CGH and tissue data defines a 12q24 gain in neuroblastic tumors with prognostic implication

Neuroblastoma has successfully served as a model system for the identification of neuroectoderm-derived oncogenes. However, in spite of various efforts, only a few clinically useful prognostic markers have been found. Here, we present a framework, which integrates DNA, RNA and tissue data to identify and prioritize genetic events that represent clinically relevant new therapeutic targets and prognostic biomarkers for neuroblastoma.

Transposable elements in a marginal plant population

Background
How new forms arise in nature has engaged evolutionary biologists since Darwin's seminal treatise on the origin of species. Transposable elements (TEs) may be among the most important internal sources for intraspecific variability. Thus, we aimed to explore the temporal dynamics of several TEs in individual genotypes from a small, marginal population of Aegilops speltoides. A diploid cross-pollinated grass species, it is a wild relative of the various wheat species known for their large genome sizes contributed by an extraordinary number of TEs, particularly long terminal repeat (LTR) retrotransposons. The population is characterized by high heteromorphy and possesses a wide spectrum of chromosomal abnormalities including supernumerary chromosomes, heterozygosity for translocations, and variability in the chromosomal position or number of 45S and 5S ribosomal DNA (rDNA) sites. We propose that variability on the morphological and chromosomal levels may be linked to variability at the molecular level and particularly in TE proliferation.

Results
Significant temporal fluctuation in the copy number of TEs was detected when processes that take place in small, marginal populations were simulated. It is known that under critical external conditions, outcrossing plants very often transit to self-pollination. Thus, three morphologically different genotypes with chromosomal aberrations were taken from a wild population of Ae. speltoides, and the dynamics of the TE complex traced through three rounds of selfing. It was discovered that: (i) various families of TEs vary tremendously in copy number between individuals from the same population and the selfed progenies; (ii) the fluctuations in copy number are TE-family specific; (iii) there is a great difference in TE copy number expansion or contraction between gametophytes and sporophytes; and (iv) a small percentage of TEs that increase in copy number can actually insert at novel locations and could serve as a bona fide mutagen.

Conclusions
We hypothesize that TE dynamics could promote or intensify morphological and karyotypical changes, some of which may be potentially important for the process of microevolution, and allow species with plastic genomes to survive as new forms or even species in times of rapid climatic change.

Molecular characterization of sediment bacterial communities affected by fish farming

Fish farming introduces nutrients, microbes and a wide variety of chemicals such as heavy metals, antifoulants and antibiotics to the surrounding environment. Introduction of antibiotics has been linked with the increased incidence of antibiotic resistant pathogenic bacteria in the farm vicinities. In this thesis molecular methods such as quantitative PCR and DNA sequencing were applied to analyze bacterial communities in sediments from fish farms and pristine locations. Altogether four farms and four pristine sites were sampled in the Baltic Sea. Two farm and two pristine locations were sampled over a surveillance period of four years. Furthermore, a new methodology was developed as a part of the study that permits amplifying single microbial genomes and capturing them according to any genetic traits, including antibiotic resistance genes. The study revealed that several resistance genes for tetracycline were found at the sediment underneath the aquaculture farms. The copy number of these genes remained elevated even at a farm that had not used any antibiotics since year 2000, six years before this study started. Similarly, an increase in the amount of mercury resistance gene merA was observed at the aquaculture sediment. The persistence of the resistance genes in absence of any selection pressure from antibiotics or heavy metals suggests that the genes may be introduced to the sediment by the farming process. This is also supported by the diversity pattern of the merA gene between farm and pristine sediments. The bacterial community-level changes in response to fish farming were very complex and no single phylogenetic groups were found that would be typical to fish farm sediments. However, the community structures had some correlation with the exposure to fish farming. Our studies suggest that the established approaches to deal with antibiotic resistance at the aquaculture, such as antibiotic cycling, are fundamentally flawed because they cannot prevent the introduction of the resistance genes and resistant bacteria to the farm area by the farming process. Further studies are required to study the entire fish farming process to identify the sources of the resistance genes and the resistant bacteria. The results also suggest that in order to prevent major microbiological changes in the surrounding aquatic environment, the farms should not be founded in shallow water where currents do not transport sedimenting matter from the farms. Finally, the technique to amplify and select microbial genomes will potentially have a considerable impact in microbial ecology and genomics.

Uniparental DNA markers and forensic genetics in Finland

Over the years, a wide range of methods to verify identity have been developed. Molecular markers have been used for identification since the 1920s, commencing with blood types and culminating with the advent of DNA techniques in the 1980s. Identification is required by authorities in many occasions, e.g. in disputed paternity cases, identification of deceased, or crime investigation. To clarify maternal and paternal lineages, uniparental DNA markers in mtDNA and Y-chromosome can be utilized. These markers have several advantages: male specific Y-chromosome can be used to identify a male from a mixture of male and female cells, e.g. in rape cases. MtDNA is durable and has a high copy number, allowing analyses even from old or degraded samples. However, both markers are lineage-specific, not individualizing, and susceptible to genetic drift. Prior to the application of any DNA marker in forensic casework, it is of utmost importance to investigate its qualities and peculiarities in the target population. Earlier studies on the Finnish population have shown reduced variation in the Y-chromosome, but in mtDNA results have been ambiguous. The obtained results confirmed the low diversity in Y-chromosome in Finland. Detailed population analysis revealed large regional differences, and extremely reduced diversity especially in East Finland. Analysis of the qualities affecting Y-chromosomal short tandem repeat (Y-STR) variation and mutation frequencies, and search of new polymorphic markers resulted a set of Y-STRs with especially high diversity in Finland. Contrary to Y-chromosome, neither reduced diversity nor regional differences were found in mtDNA within Finland. In fact, mtDNA diversity was found similar to other European populations. The revealed peculiarities in the uniparental markers are a legacy of the Finnish population history. The obtained results challenge the traditional explanation which emphasizes relatively recent founder effects creating the observed east-west patterns. Uniparentally inherited markers, both mtDNA and Y-chromosome, are applicable for identification purposes in Finland. By adjusting the analysed Y marker set to meet the characteristics of Finnish population, Y-chromosomal diversity increases and the regional differentiation decreases, resulting increase in discrimination power and thus usefulness of Y-chromosomal analysis in forensic casework.

Genomic alterations in Myeloproliferative Neoplasms and Myelodysplastic Syndromes

Myeloproliferative neoplasms (MPN) and myelodysplastic syndromes (MDS) are a heterogeneous group of clonal hematopoietic disorders whose etiology and molecular pathogenesis are poorly understood. During the past decade, enormous developments in microarray technology and bioinformatics methods have made it possible to mine novel molecular alterations in a large number of malignancies, including MPN and MDS, which has facilitated the detection of new prognostic, predictive and therapeutic biomarkers for disease stratification. By applying novel microarray techniques, we profiled copy number alterations and microRNA (miRNA) expression changes in bone marrow aspirate and blood samples. In addition, we set up and validated an miRNA expression test for bone marrow core biopsies in order to utilize the large archive material available in many laboratories. We also tested JAK2 mutation status and compare it with the in vitro growth pattern of hematologic progenitors cells. In the study focusing on 100 MPN cases, we detected a Janus kinase 2 (JAK2) mutation in 71 cases. We observed spontaneous erythroid colony growth in all mutation-positive cases in addition to nine mutation negative cases. Interestingly, seven JAK2V167F negative ET cases showed spontaneous megakaryocyte colony formation, one case of which also harbored a myeloproliferative leukemia virus oncogene (MPL) mutation. We studied copy number alterations in 35 MPN and 37 MDS cases by using oligonucleotide-based array comparative hybridization (array CGH). Only one essential thrombocythemia (ET) case presented copy number alterations in chromosomes 1q and 13q. In contrast, MDS cases were characterized by numerous novel cryptic chromosomal aberrations with the most common copy number losses at 5q21.3q33.1 and 7q22.1q33, while the most common copy number gain was trisomy 8. As for the study of the bone marrow core biopsy samples, we showed that even though these samples were embedded in paraffin and underwent decalcification, they were reliable sources of miRNA and suitable for array expression analysis. Further, when studying the miRNA expression profiles of the 19 MDS cases, we found that, compared to controls, two miRNAs (one human Epstein-Barr virus (miR-BART13) miRNA and one human (has-miR-671-5p) miRNA) were downregulated, whereas two other miRNAs (hsa-miR-720 and hsa-miR-21) were upregulated. However, we could find no correlation between copy number alterations and microRNA expression when integrating these two data. This thesis brings to light new information about genomic changes implicated in the development of MPN and MDS, and also underlines the power of applying genome-wide array screening techniques in neoplasias. Rapid advances in molecular techniques and the integration of different genomic data will enable the discovery of the biological contexts of many complex disorders, including myeloid neoplasias.

DNA damage recognition in the normal epithelium of human prostate and seminal vesicles

Prostate cancer is one of the most prevalent cancer types in men. The development of prostate tumors is known to require androgen exposure, and several pathways governing cell growth are deregulated in prostate tumorigenesis. Recent genetic studies have revealed that complex gene fusions and copy - number alterations are frequent in prostate cancer, a unique feature among solid tumors. These chromosomal aberrations are though to arise as a consequence of faulty repair of DNA double strand breaks (DSB). Most repair mechanisms have been studied in detail in cancer cell lines, but how DNA damage is detected and repaired in normal differentiated human cells has not been widely addressed. The events leading to the gene fusions in prostate cancer are under rigorous studies, as they not only shed light on the basic pathobiologic mechanisms but may also produce molecular targets for prostate cancer treatment and prevention. Prostate and seminal vesicles are part of the male reproductive system. They share similar structure and function but differ dramatically in their cancer incidence. Approximately fifty primary seminal vesicle carcinomas have been reported worldwide. Surprisingly, only little is known on why seminal vesicles are resistant to neoplastic changes. As both tissues are androgen dependent, it is a mystery that androgen signaling would only lead to tumors in prostate tissue. In this work, we set up novel ex vivo human tissue culture models of prostate and seminal vesicles, and used them to study how DNA damage is recognized in normal epithelium. One of the major DNA - damage inducible pathways, mediated by the ATM kinase, was robustly activated in all main cell types of both tissues. Interestingly, we discovered that secretory epithelial cells had less histone variant H2A.X and after DNA damage lower levels of H2AX were phosphorylated on serine 139 (γH2AX) than in basal or stromal cells. γH2AX has been considered essential for efficient DSB repair, but as there were no significant differences in the γH2AX levels between the two tissues, it seems more likely that the role of γH2AX is less important in postmitotic cells. We also gained insight into the regulation of p53, an important transcription factor that protects genomic integrity via multiple mechanisms, in human tissues. DSBs did not lead to a pronounced activation of p53, but treatments causing transcriptional stress, on the other hand, were able to launch a notable p53 response in both tissue types. In general, ex vivo culturing of human tissues provided unique means to study differentiated cells in their relevant tissue context, and is suited for testing novel therapeutic drugs before clinical trials. In order to study how prostate and seminal vesicle epithelial cells are able to activate DNA damage induced cell cycle checkpoints, we used primary cultures of prostate and seminal vesicle epithelial cells. To our knowledge, we are the first to report isolation of human primary seminal vesicle cells. Surprisingly, human prostate epithelial cells did not activate cell cycle checkpoints after DSBs in part due to low levels of Wee1A, a kinase regulating CDK activity, while primary seminal vesicle epithelial cells possessed proficient cell cycle checkpoints and expressed high levels of Wee1A. Similarly, seminal vesicle cells showed a distinct activation of the p53 - pathway after DSBs that did not occur in prostate epithelial cells. This indicates that p53 protein function is under different control mechanisms in the two cell types, which together with proficient cell cycle checkpoints may be crucial in protecting seminal vesicles from endogenous and exogenous DNA damaging factors and, as a consequence, from carcinogenesis. These data indicate that two very similar organs of male reproductive system do not respond to DNA damage similarly. The differentiated, non - replicating cells of both tissues were able to recognize DSBs, but under proliferation human prostate epithelial cells had deficient activation of the DNA damage response. This suggests that prostate epithelium is most vulnerable to accumulating genomic aberrations under conditions where it needs to proliferate, for example after inflammatory cellular damage.

Variations in paediatric upper respiratory surgery rates within and between two Nordic countries

Lasten ylähengitystiekirurgia (kita-nielurisojen poisto ja tärykalvon putkitus) on länsimaissa erittäin yleistä. Leikkausten lukumäärät vaihtelevat niin kansallisesti kuin kansainvälisestikin, mutta selvää syytä näille eroille ei tiedetä. Hoitosuositusten merkitys käytäntöihin on kyseenalaistettu ja voi olla, ettei hoitosuosituksia noudateta. Leikkaukset saattavat aiheuttaa lapsipotilaille psykologisen vamman, ja lisäksi niihin sisältyy komplikaatioiden, jopa kuoleman, vaara. Jotta haittoja voidaan välttää, on tärkeää tunnistaa ne lapset, jotka hyötyvät leikkauksesta. Ongelma on paitsi lääketieteellinen, myös taloudellinen: ylähengitystiekirurgiasta aiheutuu merkittäviä kuluja. Leikkausmäärien arvioiminen on tärkeää, jotta leikkauskäytäntöjä voidaan järkeistää. Tässä väitöskirjatyössä tutkittiin ylähengitystieleikkausten määriä Suomessa ja Norjassa sekä näiden kahden maan välillä. Aiempaa tutkimusta aiheesta ei kummassakaan maassa ole tehty. Kitarisanpoiston, välikorvan putkituksen, tärykalvopiston, nielurisanpoiston ja kita- ja nielurisanpoiston leikkausmäärät saatiin kansallisista tietokannoista. Lukuja verrattiin ko. maan lasten lukumäärään, maantieteelliseen sijoittumiseen sekä lasten ikään ja sukupuoleen. Lisäksi leikkausmääriä arvioitiin suhteessa korva-, nenä- ja kurkkulääkäreiden sekä yleislääkäreiden määrään, maantieteelliseen sijoittumiseen ja lääkäreiden ikään ja sukupuoleen. Leikkausten määrissä havaittiin suurta vaihtelua niin Suomessa kuin Norjassa. Suomessa suurimmat erot leikkausmäärissä löydettiin läntisen ja itäisen miljoonapiirin välillä. Läntisessä piirissä tehtiin lähes kaksin kertaa enemmän leikkauksia kuin itäisessä piirissä. Norjassa suurimmat erot olivat pohjoisen ja itäisen piirin välillä. Pohjoisessa piirissä tehtiin kaksinkertainen määrä leikkauksia itäiseen piirrin verrattuna. Suomessa tehtiin tutkimuksen koko aikavälillä enemmän kitarisanpoistoja kuin Norjassa, mutta ko. leikkausten määrä oli maassamme selvästi laskussa. Vuonna 2002 Suomessa tehtiin 2,5 kertaa enemmän kitarisanpoistoja kuin Norjassa. (Kita)nielurisanpoistoja tehtiin kuitenkin Suomessa vähemmän kuin Norjassa. Näiden leikkausten määrät pysyivät tutkimuksen aikavälillä Suomessa samalla tasolla, kun Norjassa leikkausmäärät hieman nousivat. Suomalaisia lapsia leikattiin keskimäärin paljon nuorempina kuin norjalaisia lapsia. Tutkimuksessa ei löydetty selitystä ylähengitystieleikkausten määrän suurelle vaihtelulle Suomessa ja Norjassa tai maiden välillä. Kuitenkin Suomessa tehtyjen kitarisanpoistojen huomattavan vähenemisen myötä maiden ylähengitystieleikkausten määrät lähenivät toisiaan.