Pharmacogenetic Variation at CYP2D6, CYP2C9, and CYP2C19: Population Genetic and Forensic Aspects

Pharmacogenetics deals with genetically determined variation in drug response. In this context, three phase I drug-metabolizing enzymes, CYP2D6, CYP2C9, and CYP2C19, have a central role, affecting the metabolism of about 20-30% of clinically used drugs. Since genes coding for these enzymes in human populations exhibit high genetic polymorphism, they are of major pharmacogenetic importance. The aims of this study were to develop new genotyping methods for CYP2D6, CYP2C9, and CYP2C19 that would cover the most important genetic variants altering the enzyme activity, and, for the first time, to describe the distribution of genetic variation at these loci on global and microgeographic scales. In addition, pharmacogenetics was applied to a postmortem forensic setting to elucidate the role of genetic variation in drug intoxications, focusing mainly on cases related to tricyclic antidepressants, which are commonly involved in fatal drug poisonings in Finland. Genetic variability data were obtained by genotyping new population samples by the methods developed based on PCR and multiplex single-nucleotide primer extension reaction, as well as by collecting data from the literature. Data consisted of 138, 129, and 146 population samples for CYP2D6, CYP2C9, and CYP2C19, respectively. In addition, over 200 postmortem forensic cases were examined with respect to drug and metabolite concentrations and genotypic variation at CYP2D6 and CYP2C19. The distribution of genetic variation within and among human populations was analyzed by descriptive statistics and variance analysis and by correlating the genetic and geographic distances using Mantel tests and spatial autocorrelation. The correlation between phenotypic and genotypic variation in drug metabolism observed in postmortem cases was also analyzed statistically. The genotyping methods developed proved to be informative, technically feasible, and cost-effective. Detailed molecular analysis of CYP2D6 genetic variation in a global survey of human populations revealed that the pattern of variation was similar to those of neutral genomic markers. Most of the CYP2D6 diversity was observed within populations, and the spatial pattern of variation was best described as clinal. On the other hand, genetic variants of CYP2D6, CYP2C9, and CYP2C19 associated with altered enzymatic activity could reach extremely high frequencies in certain geographic regions. Pharmacogenetic variation may also be significantly affected by population-specific demographic histories, as seen within the Finnish population. When pharmacogenetics was applied to a postmortem forensic setting, a correlation between amitriptyline metabolic ratios and genetic variation at CYP2D6 and CYP2C19 was observed in the sample material, even in the presence of confounding factors typical for these cases. In addition, a case of doxepin-related fatal poisoning was shown to be associated with a genetic defect at CYP2D6. Each of the genes studied showed a distinct variation pattern in human populations and high frequencies of altered activity variants, which may reflect the neutral evolution and/or selective pressures caused by dietary or environmental exposure. The results are relevant also from the clinical point of view since the genetic variation at CYP2D6, CYP2C9, and CYP2C19 already has a range of clinical applications, e.g. in cancer treatment and oral anticoagulation therapy. This study revealed that pharmacogenetics may also contribute valuable information to the medicolegal investigation of sudden, unexpected deaths.

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).

Membrane Insertion of C-tail Anchored Proteins

The correct localization of proteins is essential for cell viability. In order to achieve correct protein localization to cellular membranes, conserved membrane targeting and translocation mechanisms have evolved. The focus of this work was membrane targeting and translocation of a group of proteins that circumvent the known targeting and translocation mechanisms, the C-tail anchored protein family. Members of this protein family carry out a wide range of functions, from protein translocation and recognition events preceding membrane fusion, to the regulation of programmed cell death. In this work, the mechanisms of membrane insertion and targeting of two C-tail anchored proteins were studied utilizing in vivo and in vitro methods, in yeast and mammalian cell systems. The proteins studied were cytochrome b(5), a well characterized C-tail anchored model protein, and N-Bak, a novel member of the Bcl-2 family of regulators of programmed cell death. Membrane insertion of cytochrome b(5) into the endoplasmic reticulum membrane was found to occur independently of the known protein conducting channels, through which signal peptide-containing polypeptides are translocated. In fact, the membrane insertion process was independent of any protein components and did not require energy. Instead membrane insertion was observed to be dependent on the lipid composition of the membrane. The targeting of N-Bak was found to depend on the cellular context. Either the mitochondrial or endoplasmic reticulum membranes were targeted, which resulted in morphological changes of the target membranes. These findings indicate the existence of a novel membrane insertion mechanism for C-tail anchored proteins, in which membrane integration of the transmembrane domain, and the translocation of C-terminal fragments, appears to be spontaneous. This mode of membrane insertion is regulated by the target membrane fluidity, which depends on the lipid composition of the bilayer, and the hydrophobicity of the transmembrane domain of the C-tail anchored protein, as well as by the availability of the C-tail for membrane integration. Together these mechanisms enable the cell to achieve spatial and temporal regulation of sub-cellular localization of C-tail anchored proteins.

Archaea, Bacteria, and methane production along environmental gradients in fens and bogs

Metanogeenit ovat hapettomissa oloissa eläviä arkkien pääryhmään kuuluvia mikrobeja, joiden ainutlaatuisen aineenvaihdunnan seurauksena syntyy metaania. Ilmakehässä metaani on voimakas kasvihuonekaasu. Yksi suurimmista luonnon metaanilähteistä ovat kosteikot. Pohjoisten soiden metaanipäästöt vaihtelevat voimakkaasti eri soiden välillä ja yhden suon sisälläkin, riippuen muun muassa vuodenajasta, suotyypistä ja kasvillisuudesta. Väitöskirjatyössä tutkittiin metaanipäästöjen vaihtelun mikrobiologista taustaa. Tutkimuksessa selvitettiin suotyypin, vuodenajan, tuhkalannoituksen ja turvesyvyyden vaikutusta metanogeeniyhteisöihin sekä metaanintuottoon kolmella suomalaisella suolla. Lisäksi tutkittiin ei-metanogeenisia arkkeja ja bakteereita, koska ne muodostavat metaanin tuoton lähtöaineet osana hapetonta hajotusta. Mikrobiyhteisöt analysoitiin DNA- ja RNA-lähtöisillä, polymeraasiketjureaktioon (PCR) perustuvilla menetelmillä. Merkkigeeneinä käytettiin metaanin tuottoon liittyvää mcrA-geeniä sekä arkkien ja bakteerien ribosomaalista 16S RNA-geeniä. Metanogeeniyhteisöt ja metaanintuotto erosivat huomattavasti happaman ja vähäravinteisen rahkasuon sekä ravinteikkaampien sarasoiden välillä. Rahkasuolta löytyi lähes yksinomaan Methanomicrobiales-lahkon metanogeeneja, jotka tuottavat metaania vedystä ja hiilidioksidista. Sarasoiden metanogeeniyhteisöt olivat monimuotoisempia, ja niillä esiintyi myös asetaattia käyttäviä metanogeeneja. Vuodenaika vaikutti merkittävästi metaanintuottoon. Talvella havaittiin odottamattoman suuri metaanintuottopotentiaali sekä viitteitä aktiivisista metanogeeneista. Arkkiyhteisön koostumus sen sijaan vaihteli vain vähän. Tuhkalannoitus, jonka tarkoituksena on edistää puiden kasvua ojitetuilla soilla, ei merkittävästi vaikuttanut metaanintuottoon tai -tuottajiin. Ojitetun suon yhteisöt kuitenkin muuttuivat turvesyvyyden mukaan. Vertailtaessa erilaisia PCR-menetelmiä todettiin, että kolmella mcrA-geeniin kohdistuvalla alukeparilla havaittiin pääosin samat ojitetun suon metanogeenit, mutta lajien runsaussuhteet riippuvat käytetyistä alukkeista. Soilla havaitut bakteerit kuuluivat pääjaksoihin Deltaproteobacteria, Acidobacteria ja Verrucomicrobia. Lisäksi löydettiin Crenarchaeota-pääjakson ryhmiin 1.1c ja 1.3 kuuluvia ei-metanogeenisia arkkeja. Tulokset ryhmien esiintymisestä hapettomassa turpeessa antavat lähtökohdan selvittää niiden mahdollisia vuorovaikutuksia metanogeenien kanssa. Tutkimuksen tulokset osoittivat, että metanogeeniyhteisön koostumus heijastaa metaanintuottoon vaikuttavia kemiallisia tai kasvillisuuden vaihteluita kuten suotyyppiä. Soiden metanogeenien ja niiden fysiologian parempi tuntemus voi auttaa ennustamaan ympäristömuutosten vaikutusta soiden metaanipäästöihin.

Real-Time PCR : A Molecular Approach to Investigate the Role of Intestinal Microbiota in the Pathophysiology of Irritable Bowel Syndrome

Irritable bowel syndrome (IBS) is a common multifactorial functional intestinal disorder, the pathogenesis of which is not completely understood. Increasing scientific evidence suggests that microbes are involved in the onset and maintenance of IBS symptoms. The microbiota of the human gastrointestinal (GI) tract constitutes a massive and complex ecosystem consisting mainly of obligate anaerobic microorganisms making the use of culture-based methods demanding and prone to misinterpretation. To overcome these drawbacks, an extensive panel of species- and group-specific assays for an accurate quantification of bacteria from fecal samples with real-time PCR was developed, optimized, and validated. As a result, the target bacteria were detectable at a minimum concentration range of approximately 10 000 bacterial genomes per gram of fecal sample, which corresponds to the sensitivity to detect 0.000001% subpopulations of the total fecal microbiota. The real-time PCR panel covering both commensal and pathogenic microorganisms was assessed to compare the intestinal microbiota of patients suffering from IBS with a healthy control group devoid of GI symptoms. Both the IBS and control groups showed considerable individual variation in gut microbiota composition. Sorting of the IBS patients according to the symptom subtypes (diarrhea, constipation, and alternating predominant type) revealed that lower amounts of Lactobacillus spp. were present in the samples of diarrhea predominant IBS patients, whereas constipation predominant IBS patients carried increased amounts of Veillonella spp. In the screening of intestinal pathogens, 17% of IBS samples tested positive for Staphylococcus aureus, whereas no positive cases were discovered among healthy controls. Furthermore, the methodology was applied to monitor the effects of a multispecies probiotic supplementation on GI microbiota of IBS sufferers. In the placebo-controlled double-blind probiotic intervention trial of IBS patients, each supplemented probiotic strain was detected in fecal samples. Intestinal microbiota remained stable during the trial, except for Bifidobacterium spp., which increased in the placebo group and decreased in the probiotic group. The combination of assays developed and applied in this thesis has an overall coverage of 300-400 known bacterial species, along with the number of yet unknown phylotypes. Hence, it provides good means for studying the intestinal microbiota, irrespective of the intestinal condition and health status. In particular, it allows screening and identification of microbes putatively associated with IBS. The alterations in the gut microbiota discovered here support the hypothesis that microbes are likely to contribute to the pathophysiology of IBS. The central question is whether the microbiota changes described represent the cause for, rather than the effect of, disturbed gut physiology. Therefore, more studies are needed to determine the role and importance of individual microbial species or groups in IBS. In addition, it is essential that the microbial alterations observed in this study will be confirmed using a larger set of IBS samples of different subtypes, preferably from various geographical locations.

Molecular Characterization of Viruses Causing the Cassava Brown Streak Disease Epidemic in Eastern Africa

Cassava brown streak disease (CBSD) was described for the first time in Tanganyika (now Tanzania) about seven decades ago. Tanganyika (now Tanzania) about seven decades ago. It was endemic in the lowland areas of East Africa and inland parts of Malawi and caused by Cassava brown streak virus (CBSV; genus Ipomovirus; Potyviridae). However, in 1990s CBSD was observed at high altitude areas in Uganda. The causes for spread to new locations were not known.The present work was thus initiated to generate information on genetic variability, clarify the taxonomy of the virus or viruses associated with CBSD in Eastern Africa as well as to understand the evolutionary forces acting on their genes. It also sought to develop a molecular based diagnostic tool for detection of CBSD-associated virus isolates. Comparison of the CP-encoding sequences of CBSD-associated virus isolates collected from Uganda and north-western Tanzania in 2007 and the partial sequences available in Genbank revealed occurrence of two genetically distinct groups of isolates. Two isolates were selected to represent the two groups. The complete genomes of isolates MLB3 (TZ:Mlb3:07) and Kor6 (TZ:Kor6:08) obtained from North-Western (Kagera) and North-Eastern (Tanga) Tanzania, respectively, were sequenced. The genomes were 9069 and 8995 nucleotides (nt), respectively. They translated into polyproteins that were predicted to yield ten mature proteins after cleavage. Nine proteins were typical in the family Potyviridae, namely P1, P3, 6K1, CI, 6K2, VPg, NIa-Pro, NIb and CP, but the viruses did not contain HC-Pro. Interestingly, genomes of both isolates contained a Maf/HAM1-like sequence (HAM1h; 678 nucleotides, 25 kDa) recombined between the NIb and CP domains in the 3’-proximal part of the genomes. HAM1h was also identified in Euphorbia ringspot virus (EuRSV) whose sequence was in GenBank. The HAM1 gene is widely spread in both prokaryotes and eukaryotes. In yeast (Saccharomyces cerevisiae) it is known to be a nucleoside triphosphate (NTP) pyrophosphatase. Novel information was obtained on the structural variation at the N-termini of polyproteins of viruses in the genus Ipomovirus. Cucumber vein yellowing virus (CVYV) and Squash vein yellowing virus (SqVYV) contain a duplicated P1 (P1a and P1b) but lack the HC-Pro. On the other hand, Sweet potato mild mottle virus (SPMMV), has a single but large P1 and has HC-Pro. Both virus isolates (TZ:Mlb3:07 & TZ:Kor6:08) characterized in this study contained a single P1 and lacked the HC-Pro which indicates unique evolution in the family Potyviridae. Comparison of 12 complete genomes of CBSD-associated viruses which included two genomes characterized in this study, revealed genetic identity of 69.0–70.3% (nt) and amino acid (aa) identities of 73.6–74.4% at polyprotein level. Comparison was also made among 68 complete CP sequences, which indicated 69.0-70.3 and 73.6-74.4 % identity at nt and aa levels, respectively. The genetic variation was large enough for dermacation of CBSD-associated virus isolates into two distinct species. The name CBSV was retained for isolates that were related to CBSV isolates available in database whereas the new virus described for the first time in this study was named Ugandan cassava brown streak virus (UCBSV) by the International Committee on Virus Taxonomy (ICTV). The isolates TZ:Mlb3:07 and TZ:Kor6:08 belong to UCBSV and CBSV, respectively. The isolates of CBSV and UCBSV were 79.3-95.5% and 86.3-99.3 % identitical at nt level, respectively, suggesting more variation amongst CBSV isolates. The main sources of variation in plant viruses are mutations and recombination. Signals for recombination events were detected in 50% of isolates of each virus. Recombination events were detected in coding and non-coding (3’-UTR) sequences except in the 5’UTR and P3. There was no evidence for recombination between isolates of CBSV and UCBSV. The non-synonomous (dN) to synonomous (dS) nucleotide substitution ratio (ω) for the HAM1h and CP domains of both viruses were ≤ 0.184 suggesting that most sites of these proteins were evolving under strong purifying selection. However, there were individual amino acid sites that were submitted to adaptive evolution. For instance, adaptive evolution was detected in the HAM1h of UCBSV (n=15) where 12 aa sites were under positive selection (P< 0.05) but not in CBSV (n=12). The CP of CBSV (n=23) contained 12 aa sites (p<0.01) while only 5 aa sites in the CP gene of UCBSV were predicted to be submitted to positive selection pressure (p<0.01). The advantages offered by the aa sites under positive selection could not be established but occurrence of such sites in the terminal ends of UCBSV-HAMIh, for example, was interpreted as a requirement for proteolysis during polyprotein processing. Two different primer pairs that simultaneously detect UCBSV and CBSV isolates were developed in this study. They were used successfully to study distribution of CBSV, UCBSV and their mixed infections in Tanzania and Uganda. It was established that the two viruses co-infect cassava and that incidences of co-infection could be as high as 50% around Lake Victoria on the Tanzanian side. Furthermore, it was revealed for the first time that both UCBSV and CBSV were widely distributed in Eastern Africa. The primer pair was also used to confirm infection in a close relative of cassava, Manihot glaziovii (Müller Arg.) with CBSV. DNA barcoding of M. glaziovii was done by sequencing the matK gene. Two out of seven M. glaziovii from the coastal areas of Korogwe and Kibaha in north eastern Tanzania were shown to be infected by CBSV but not UCBSV isolates. Detection in M. glaziovii has an implication in control and management of CBSD as it is likely to serve as virus reservoir. This study has contributed to the understanding of evolution of CBSV and UCBSV, which cause CBSD epidemic in Eastern Africa. The detection tools developed in this work will be useful in plant breeding, verification of the phytosanitary status of materials in regional and international movement of germplasm, and in all diagnostic activities related to management of CBSD. Whereas there are still many issues to be resolved such as the function and biological significance of HAM1h and its origin, this work has laid a foundation upon which the studies on these aspects can be based.