Browsing by study line "Genetik och genomik"

In recent decades, ancient DNA recovered from old and degraded samples, such as bones and fossils, has presented novel prospects in the fields of genetics, archaeology and anthropology. In Finland, ancient DNA research is constrained by the poor preservation of bones: they are quickly degraded by acidic soils, limiting the age of DNA that can be recovered from physical remains. However, some soil components can bind DNA and thus protect the molecules from degradation. Ancient DNA from soils and sediments has previously been used to reconstruct paleoenvironments, to study ancient parasites and diet and to demonstrate the presence of a species at a given site, even when there are no visible fossils present. In this pilot study, I explored the potential of archaeological sediments as an alternative source of ancient human DNA. I collected sediment samples from five Finnish Neolithic Stone Age (6,000–4,000 years ago) settlement sites, located in woodland. In addition, I analysed a lakebed sample from a submerged Mesolithic (10,000–7,000 years ago) settlement site, and a soil sample from an Iron Age burial with bones present to compare DNA yields between the two materials. Soil samples were converted into Illumina sequencing libraries and enriched for human mtDNA. I analysed the sequencing data with a customised metagenomics-based bioinformatic analysis workflow. I also tested program performance with simulated data. The results suggested that human DNA preservation in Finnish archaeological sediments may be poor or very localised. I detected small amounts of human mtDNA in three Stone Age woodland settlement sites and a submerged Mesolithic settlement site. One Stone Age sample exhibited terminal damage patterns suggestive of DNA decay, but the time of deposition is difficult to estimate. Interestingly, no human DNA was recovered from the Iron Age burial soil, suggesting that body decomposition may not serve as a significant source of sedimentary ancient DNA. Additional complications may arise from the high inhibitor content of the soil and the abundance of microbial and other non-human DNA present in environmental samples. In the future, a more refined sampling approach, such as targeting microscopic bone fragments, could be a strategy worth trialling.

The Y chromosome has an essential role in the genetic sex determination in humans and other mammals. It contains a male-specific region (MSY) which escapes recombination and is inherited exclusively through the male line. The genetic variations inherited together on the MSY can be used in classifying Y chromosomes into haplogroups. Y-chromosomal haplogroups are highly informative of genetic ancestry, thus Y chromosomes have been widely used in tracing human population history. However, given the peculiar biology and analytical challenges specific to the Y chromosome, the chromosome is routinely excluded from genetic association studies. Consequently, potential impacts of Y-chromosomal variation on complex disease remain largely uncharacterized. Lately the access to large-scale biobank data has enabled to extend the Y-chromosomal genetic association studies. A recent UK Biobank study suggested links between Y-chromosomal haplogroup I1 and coronary artery disease (CAD) in the British population, but this result has not been validated in other datasets. Since Finland harbours a notable frequency of Y-chromosomal haplogroup I1, the relationship between haplogroup I1 and CAD can further be inferred in the Finnish population using data from the FinnGen project. The first aim of this thesis was to determine the prevalence of Y-chromosomal haplogroups in Finland and characterize their geographical distributions using genotyping array data from the FinnGen project. The second aim was to assess the role between Finnish Y-chromosomal haplogroups and coronary artery disease (CAD) by logistic regression. This thesis characterized the Y-chromosomal haplogroups in Finland for 24 160 males and evaluated the association between Y-chromosomal haplogroups and CAD in Finland. The dataset used in this study was extensive, providing an opportunity to study the Y-chromosomal variation geographically in Finland and its role in complex disease more accurately compared to previous studies. The geographical distribution of the Y-chromosomal haplogroups was characterized on 20 birth regions, and between eastern and western areas of Finland. Consistent with previous studies, the results demonstrated that two major Finnish Y-chromosomal haplogroup lineages, N1c1 and I1, displayed differing distributions within regions, especially between eastern and western Finland. Results from logistic regression analysis between CAD and Y-chromosomal haplogroups suggested no significant association between haplogroup I1 and CAD. Instead, the major Finnish Y-chromosomal haplogroup N1c1 displayed a decreased risk for CAD in the association analysis when compared against other haplogroups. Moreover, this thesis also demonstrated that the association results were not straightforwardly comparable between populations. For instance, haplogroup I1 displayed a decreased risk for CAD in the FinnGen dataset when compared against haplogroup R1b, whereas the same association was reported as risk increasing for CAD in the UK Biobank. Overall, this thesis demonstrates the possibility to study the genetics of Y chromosome using data from the FinnGen project, and highlights the value of including this part of the genome in the future complex disease studies.

Uterine leiomyomas are benign tumors originating in the smooth muscle cells of the uterine wall. Leiomyomas represent one of the most common tumor types in women affecting up to 80% of pre-menopausal women. Besides having extensive implications on women´s health through the numerous symptoms they cause, leiomyomas are a cause of remarkable financial burden worldwide. Bivalent promoters are defined by the co-occurrence of two histone modifications with opposite functions: trimethylation of lysine 4 on histone 3 (H3K4me3) and trimethylation of lysine 27 on histone 3 (H3K27me3). H3K4me3 is associated with promoters of actively expressed genes, whereas H3K27me3 is frequently found at promoters of silenced genes. The genes controlled by the bivalent promoters are reversibly silenced or expressed at low levels and remain poised for fast activation or full repression as a response to external cues. Bivalent chromatin is gaining more and more importance as new roles are identified in tumorigenesis and cell differentiation. Despite this, the vast majority of data available was obtained from cell lines, and not from human tissue. The aim of this thesis work was to map the genomic location of bivalent promoters in uterine leiomyoma and myometrium tissue, and to characterize the functions of bivalently-controlled genes in differentiated tissue. This would provide novel information about bivalent promoters’ distribution in human tissues and also their potential role in myomagenesis. Chromatin immunoprecipitation followed by sequencing (ChIP-seq) against H3K4me3 and H3K27me3 was performed on fresh frozen tissue samples of uterine leiomyomas and corresponding myometrium. A promoter was defined as bivalent, if it showed overlap between H3K4me3 and H3K27me3 peaks within a 2 kb region of a gene’s transcription start site in all samples. Altogether 951 bivalent promoters were found in myometrium and leiomyoma. Strikingly, only 231 (24.3%) promoters were present in both tissue types, most bivalent promoters being tissue-specific. These findings indicated bivalent promoters regulating a substantial number of genes also in differentiated tissue and the presence of extensive alterations in bivalent promoter distribution during myomagenesis. Gene ontology analyses of the bivalently-controlled genes in myometrium revealed the highest score for developmental processes. Instead, for leiomyomas, the highest enrichment was detected in stem cell fate specification-related processes. The data presented in this thesis suggests that bivalent chromatin plays an important role during myomagenesis, as it undergoes a significant reorganization during the process. Future experiments will provide novel insights about the role for these changes, i.e.: if they underlie the process.

Endometrial polyps are one of the most common benign uterine lesions, affecting approximately 10% of all adult women. While endometrial polyps have a high prevalence, their molecular pathogenesis and genetic background are largely undefined. Accordingly, the aim of this thesis was to characterize the somatic mutational landscape of endometrial polyps – to identify mutations in cancer-associated genes, and to identify mutational signatures contributing towards the somatic mutational spectrum. The present study was conducted using whole exome sequencing of 23 endometrial polyps and 18 matching normal blood samples. Mutational signature analysis was conducted using MutationalPatterns and SigProfiler. Endometrial polyps were found to carry varying number of somatic mutations in their exomes, most of them present at a low allelic fraction. Moreover, 43% (10/23) of the polyps were identified to carry one to four cancer-associated mutations, including mutations in genes such as PIK3CA 17% (4/23), KRAS 13% (3/23) and ERBB1 9% (2/23), which are well-established cancer driver genes. Cancer-associated mutational signatures do not have a notable contribution towards the somatic mutational spectrum of endometrial polyps. However, a novel signature, ‘signature B’, characterized by T>G mutations, was found to affect a subset of polyp samples. To conclude, the whole exome sequencing of endometrial polyps revealed several mutations in cancer-associated genes and a novel mutational signature, which may contribute to the development of these benign tumours. However, further research is required to confirm and validate the novel signature, and to define the genetic alterations leading to the polyp pathogenesis.

Enhancers are important regulatory elements of DNA, that are bound by transcription factors (TFs) to regulate gene expression. Enhancers control cell type specific gene expression and they can form structures called super-enhancers, that consist of multiple normal enhancers and are bound by high numbers and variety of transcription factors. These super-enhancers are important for defining cell identity and changes in the super-enhancer landscape have been linked to different cancers. In this project, characterization of super-enhancers and their transcription factors composition between primary and cancer cells were studied using genome-wide next-generation sequencing data from multiple assays, such as ChIP-seq, RNA-seq and ATAC-seq. The focus of the project was on the data processing and analysis to identify and characterize the super-enhancers. Analyses included GSEA, heatmap binding analysis, peak and super-enhancer calling and IGV analysis. This project used pancreatic HPDE cell line for primary cells and different cancers with endodermal origin as cancer cell lines. The goal of the thesis was to try show characteristic features of super-enhancers and their features in normal and cancer cells. Data analysis showed that distinct super-enhancers can be identified in cancer cells and defined super-enhancers had typical strong binding for specific transcription factor and histone modification such as histone 3 lysine 27 acetylation (H3K27ac) mark of active enhancers. Super-enhancer regions were located in highly accessible chromatin regions of the genome, and genes that were associated with HPDE super-enhancers could be shown to have association with cell identity. Peak and super-enhancer calling counts varied between cell lines for transcription factors, histone modifications and super-enhancers. Visualization of super-enhancers was successful and could show transcription factor binding and active enhancers that establish the super-enhancer structure. Comprehensive analyses allowed us to characterize typical features of super-enhancers and show differences in the numbers of super-enhancers between primary and cancer cell lines and cancer cell lines of different organ types. Analysis of the transcription factor binding showed unique peaks on some of the super-enhancers, and these peaks might have a role in inducing the super-enhancer structure.

Cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL) is an inherited autosomal dominant disease that leads to cognitive impairment, vascular dementia and ischemic strokes. In CADASIL, vascular smooth muscle cells (VSMCs) degrade gradually and are replaced by connective tissue in the small and mid-sized arteries in the brain. Extracellular granular osmiophilic material (GOM) that surround the VSMCs are a unique feature in CADASIL. The causal gene behind CADASIL is Notch3, which encodes a transmembrane protein with a signaling function. There are over 200 cysteine-altering mutations that cause CADASIL in Notch3. The potential pathology causing mechanism is still unclear, but most likely the mechanism is linked to the aggregation of GOM deposits that are potentially toxic to VSMCs. This thesis project aimed to correct CADASIL causing c.475C>T mutation in Notch3 in different CADASIL cell lines with different CRISPR base editor systems. Another aim was to create induced pluripotent stem cell (iPSC) lines from a CADASIL patient-derived skin biopsy sample to be used in the creation of an in vitro disease model for CADASIL. RNA-based ABEmax base editor system was used to correct immortalized- and primary- CADASIL cell lines. DNA-based ABEmax base editor system was used as a positive control. Simultaneous pluripotent reprogramming and pathogenic CADASIL mutation correction were done in the same transfection during this project. The editing efficiencies were evaluated by Sanger sequencing the genomic target region before and after the transfection. The editing efficiencies were good in general compared to literature. They ranged from 27 % to 73 % target base editing efficiency depending on the editing system-, guide-RNAs - and electroporation parameters used. Confirmed proximal off-target effects were not detected, and distal off-target effects were not evaluated.

Drought events will increase in the future due to climate change and thereby threaten agricultural production. The growing world population needs stable yields nonetheless thus, solutions need to be found. One way to do so is to explore mechanisms that allow plants to withstand harsh conditions and recover from them. Autophagy is such a mechanism. It is a bulk degradation pathwaythat allows plants to recycle cell components. During biotic and abiotic stress autophagy pathways are enhanced. Even though this mechanism is very useful and is proven to enhance abiotic stress tolerance, very little is known about it in plants. This thesis aimed at defining the informative time points for drought and autophagy response during drought and recovery in barley. To do so barley plants were subjected to gradual drought by withholding irrigation. They were kept at a level of severe drought for 4 days and were then rewatered. The tip of the third leaf was sampled for RNA extraction followed by qPCR. The genes this thesis focused on were GST and HSP17 as indicators of drought stress and ATG6 as indicator for autophagy. Additionally, phenotypic data was collected via RGB imaging for monitoring wilting and thermal imaging for stomata closure. GST was upregulated during the onset of drought, which seems to make it an early response gene, while HSP17 was upregulated during severe drought. ATG6 had a high transcription rate during all stages of drought and only decreased during recovery. However, because the primer efficiencies calculated in this experiment appear to be very low, and data points were missing due to pipetting practice, the results are not reliable and need to be repeated. Thus, no prediction about the informative time points can be made. The thermal data showed good results and the expected pattern of stomatal closing during drought could be observed. But because the greenhouse is a variable environment a reference would make the data even more precise.

Multiple myeloma (MM) is a heterogeneous plasma cell cancer that results from the excessive proliferation of mutated B cells in the bone marrow and the accumulation of ineffective antibodies, monoclonal proteins, in the blood. Despite recent advances in research and novel therapeutics, MM remains incurable, mainly due to the mechanisms underlying disease progression and drug resistance. Therefore, novel biomarkers and therapeutics for the treatment of relapsed and refractory MM are urgently needed. MicroRNAs (miRNAs), short non-coding RNA molecules that play a key role in post-transcriptional gene regulation, have been found to be associated with different hallmarks of MM. Previous studies have indicated that abnormally functioning miRNA-mediated gene regulation followed by oncogene activation and tumor suppressor gene silencing results in drastic alterations in cell proliferation, apoptosis, growth, and metabolism. These changes in cellular functions have been indicated to be associated with the pathogenesis, progression, and formation of drug resistance in MM. Therefore, the role and potential of miRNAs to act as biomarkers to predict MM progression and drug sensitivity should be further investigated to ultimately improve the survival rates of patients. The aim of this master’s thesis was to investigate the relationships between drug sensitivity, disease progression and miRNA regulation in MM patients. Bioinformatically predicted miRNAs identified to be associated with sensitivity to panobinostat, a novel histone deacetylase inhibitor, and MM progression were validated in MM patient samples by using real-time quantitative reverse transcription PCR (RT-qPCR). In addition, the specific gene targets of miRNAs involved in the regulation of drug responses and MM progression were predicted by identifying statistically significant, negatively correlated interactions between the miRNA and RNA sequencing data of 45 MM patients in pairwise comparative correlation analysis. Finally, the predicted miRNA targets genes were validated in MM patient samples using RT-qPCR. Based on the bioinformatic analyses and RT-qPCR validation, mir-424 expression was significantly increased in relapsed MM patients as compared to respective patient samples taken at diagnosis, suggesting a potential role of mir-424 in MM progression. Similarly, mir-4433b expression was significantly elevated in panobinostat-resistant patients compared to sensitive patients, suggesting a potential effect of mir-4433b on the regulation of panobinostat drug response in MM patients. In addition, the RT-qPCR validation demonstrated that the disease progression and drug sensitivity associated mir-92b, mir-363 and mir-221, would potentially regulate the expression of FGF2, MFF, and TMEM248, respectively, providing novel insights into the functional roles of miRNAs in MM pathways.

The European rabbit (Oryctolagus cuniculus) is a small mammal native to the Iberian Peninsula, but introduced by humans to all continents except Antarctica. The rabbit has been a remarkably successful invasive species due to its generalist nature and fast reproduction. Its spreading has mostly been destructive to the local nature, and humans have used fatal rabbit diseases such as rabbit haemorrhagic disease (RHD) to control harmful populations. The rabbit population in Helsinki is one of the most northern annually surviving rabbit populations in the world. It is believed to have originated from escaped pet rabbits in the late 1980s, and in the early 2000s, the rabbits spread rapidly around the Helsinki area. RHD spread unintentionally to Finland in 2016, and the disease caused a significant reduction in the Helsinki rabbit population. Rabbit population genetics has previously been studied in several countries, but never before in Finland. The aim of the thesis was to examine the genetic diversity and population structure of the Helsinki rabbit population before and after the RHD epidemic, and to compare the results to similar preceding rabbit population genetic studies. Rabbit populations have previously been found to recover from major population crashes without a notable loss in genetic diversity using DNA microsatellite markers. The recent RHD epidemic in Helsinki provided an opportunity to study, whether a rabbit population can recover from a population crash even in a harsher environment without losing genetic diversity. To conduct genetic analysis, fourteen DNA microsatellite loci were genotyped from individuals caught during two distinct time periods, in 2008-2009 (n=130) and in 2019-2020 (n=59). Population structure was observed in both temporal rabbit populations with small but significant FST values. The 2019-2020 population was more diverse than the 2008-2009 population in terms of allele numbers and expected heterozygosity. This result was unexpected considering the recent RHD-epidemic but could be explained by gene flow from new escaped rabbits. Compared to other wild rabbit populations around the world, the Helsinki area rabbits exhibit significantly lower genetic diversity. Bottleneck tests showed a significant signal separately in both temporal populations, but the RHD bottleneck cannot be distinguished based on the tests. The results could be biased by new gene flow, or the initial bottleneck caused by the founder effect of only a few pet rabbits. The rabbits have demonstrated their adaptation and survival skills in the cold climate of Helsinki. The population has significantly lower genetic diversity compared to other wild populations, yet recovered from a major RHD epidemic without reduction in genetic diversity under these more extreme environmental conditions. It has been proven again; the rabbit is a thriving invasive species.

Silver birch (Betula pendula) is both the third most commercially important, as well as abundant, tree in Finland. Faster growth is of interest from an industry, as well as climate change, perspective. In addition to higher biomass yield from a commercial point of view, faster growth in trees leads to absorption of more carbon dioxide, thus increasing their effect as a carbon sink to lessen the impact of climate change. By conducting a genome-wide association study, the aim of this thesis was to investigate the genetical background of B. pendula on three first-year growth traits: maximum height, maximum growth rate and time to cessation. Whole-genome sequencing data from 405 mother trees was used to provide a high-resolution genome-wide map for assessing genes associated with the three growth traits. A high-quality genotypic dataset of 5,727,473 single-nucleotide polymorphisms (SNP), which span across the whole genome, was established through rigorous filtering. Phenotypic data was estimated with a mixed linear model by best linear unbiased prediction from mother trees’ 7,266 seedlings, which were grown in a greenhouse under controlled conditions to monitor their height trajectories. Association testing was done as both a uni- and multivariate mixed linear model with the GEMMA algorithm. Overall, 16 unique suggestively significant SNPs were discovered among the three traits. Genes putatively associated with the SNPs were related to metabolic regulation and cell transportation processes, as well as biotic and abiotic stress. On each of the three traits, SNPs with the largest effects had around 5% deviation against the average phenotypic value, which could be considered a major effect for a polygenic trait. Only 4 out of the 16 suggestively associated SNPs were within gene regions, even though the putative genes themselves had also variants present abundantly in the dataset. This could suggest that the associated SNPs could be related to gene expression regulation, or they could be in linkage disequilibrium with multiple different genes, and thus be markers for combined effects of multiple genes on the focal trait. To validate the results, further studies such as gene editing on the associated genes, or a duplicated study on a different population, would be needed.