Genetic aspects of outcome in kidney transplantation: cytokine and thrombosis associated candidate genes and gene expression biomarkers

Kidney transplantation (Tx) is the treatment of choice for end stage renal disease. Immunosuppressive medications are given to prevent an immunological rejection of the transplant. However, immunosuppressive drugs increase e.g. the risk of infection, cancer or nephrotoxicity. A major genetic contributors to immunological acceptance of the graft are human leukocyte antigen (HLA) genes. Also other non-HLA gene polymorphisms may predict the future risk of complications before Tx, possibly enabling individualised immunotherapy. Graft function after Tx is monitored using non-specific clinical symptoms and laboratory markers. The definitive diagnosis of graft rejection however relies on a biopsy of the graft. In the acute rejection (AR) diagnostics there is a need for an alternative to biopsy that would be an easily repeatable and simple method for regular use. Frequent surveillance of acute or subclinical rejection (SCR) may improve long-term function. In this thesis, associations between cytokine and thrombosis associated candidate genes and the outcome of kidney Tx were studied. Cytotoxic and co-stimulatory T lymphocyte molecule gene expression biomarkers for the diagnosis of the AR and the SCR were also investigated. We found that polymorphisms in the cytokine genes tumor necrosis factor and interleukin 10 (IL10) of the recipients were associated with AR. In addition, certain IL10 gene polymorphisms of the donors were associated with the incidence of cytomegalovirus infection and occurrence of later infection in a subpopulation of recipients. Further, polymorphisms in genes related to the risk of thrombosis and those of certain cytokines were not associated with the occurrence of thrombosis, infarction, AR or graft survival. In the study of biomarkers for AR, whole blood samples were prospectively collected from adult kidney Tx patients. With real-time quantitative PCR (RT-QPCR) gene expression quantities of CD154 and ICOS differentiated the patients with AR from those without, but not from the patients with other causes of graft dysfunction. Biomarkers for SCR were studied in paediatric kidney Tx patients. We used RT-QPCR to quantify the gene expression of immunological candidate genes in a low-density array format. In addition, we used RT-QPCR to validate the results of the microarray analysis. No gene marker differentiated patients with SCR from those without SCR. This research demonstrates the lack of robust markers among polymorphisms or biomarkers in investigated genes that could be included in routine analysis in a clinical laboratory. In genetic studies, kidney Tx can be regarded as a complex trait, i.e. several environmental and genetic factors may determine its outcome. A number of currently unknown genetic factors probably influence the results of Tx.

Applications of gene sequence polymorphisms in evolutionary genetic studies of Atlantic salmon (Salmo salar) and other teleost fish species

Evolutionary genetics incorporates traditional population genetics and studies of the origins of genetic variation by mutation and recombination, and the molecular evolution of genomes. Among the primary forces that have potential to affect the genetic variation within and among populations, including those that may lead to adaptation and speciation, are genetic drift, gene flow, mutations and natural selection. The main challenges in knowing the genetic basis of evolutionary changes is to distinguish the adaptive selection forces that cause existent DNA sequence variants and also to identify the nucleotide differences responsible for the observed phenotypic variation. To understand the effects of various forces, interpretation of gene sequence variation has been the principal basis of many evolutionary genetic studies. The main aim of this thesis was to assess different forms of teleost gene sequence polymorphisms in evolutionary genetic studies of Atlantic salmon (Salmo salar) and other species. Firstly, the level of Darwinian adaptive evolution affected coding regions of the growth hormone (GH) gene during the teleost evolution was investigated based on the sequence data existing in public databases. Secondly, a target gene approach was used to identify within population variation in the growth hormone 1 (GH1) gene in salmon. Then, a new strategy for single nucleotide polymorphisms (SNPs) discovery in salmonid fishes was introduced, and, finally, the usefulness of a limited number of SNP markers as molecular tools in several applications of population genetics in Atlantic salmon was assessed. This thesis showed that the gene sequences in databases can be utilized to perform comparative studies of molecular evolution, and some putative evidence of the existence of Darwinian selection during the teleost GH evolution was presented. In addition, existent sequence data was exploited to investigate GH1 gene variation within Atlantic salmon populations throughout its range. Purifying selection is suggested to be the predominant evolutionary force controlling the genetic variation of this gene in salmon, and some support for gene flow between continents was also observed. The novel approach to SNP discovery in species with duplicated genome fragments introduced here proved to be an effective method, and this may have several applications in evolutionary genetics with different species - e.g. when developing gene-targeted markers to investigate quantitative genetic variation. The thesis also demonstrated that only a few SNPs performed highly similar signals in some of the population genetic analyses when compared with the microsatellite markers. This may have useful applications when estimating genetic diversity in genes having a potential role in ecological and conservation issues, or when using hard biological samples in genetic studies as SNPs can be applied with relatively highly degraded DNA.

Comparative and functional genome analysis of fungi for development of the protein production host Trichoderma reesei

Filamentous fungi of the subphylum Pezizomycotina are well known as protein and secondary metabolite producers. Various industries take advantage of these capabilities. However, the molecular biology of yeasts, i.e. Saccharomycotina and especially that of Saccharomyces cerevisiae, the baker's yeast, is much better known. In an effort to explain fungal phenotypes through their genotypes we have compared protein coding gene contents of Pezizomycotina and Saccharomycotina. Only biomass degradation and secondary metabolism related protein families seem to have expanded recently in Pezizomycotina. Of the protein families clearly diverged between Pezizomycotina and Saccharomycotina, those related to mitochondrial functions emerge as the most prominent. However, the primary metabolism as described in S. cerevisiae is largely conserved in all fungi. Apart from the known secondary metabolism, Pezizomycotina have pathways that could link secondary metabolism to primary metabolism and a wealth of undescribed enzymes. Previous studies of individual Pezizomycotina genomes have shown that regardless of the difference in production efficiency and diversity of secreted proteins, the content of the known secretion machinery genes in Pezizomycotina and Saccharomycotina appears very similar. Genome wide analysis of gene products is therefore needed to better understand the efficient secretion of Pezizomycotina. We have developed methods applicable to transcriptome analysis of non-sequenced organisms. TRAC (Transcriptional profiling with the aid of affinity capture) has been previously developed at VTT for fast, focused transcription analysis. We introduce a version of TRAC that allows more powerful signal amplification and multiplexing. We also present computational optimisations of transcriptome analysis of non-sequenced organism and TRAC analysis in general. Trichoderma reesei is one of the most commonly used Pezizomycotina in the protein production industry. In order to understand its secretion system better and find clues for improvement of its industrial performance, we have analysed its transcriptomic response to protein secretion stress conditions. In comparison to S. cerevisiae, the response of T. reesei appears different, but still impacts on the same cellular functions. We also discovered in T. reesei interesting similarities to mammalian protein secretion stress response. Together these findings highlight targets for more detailed studies.

Regulation of Osteoblast Differentiation and Gene Expression by Phosphodiesterases and Bioactive Peptides

Bone is a mineralized tissue that enables multiple mechanical and metabolic functions to be carried out in the skeleton. Bone contains distinct cell types: osteoblasts (bone-forming cells), osteocytes (mature osteoblast that embedded in mineralized bone matrix) and the osteoclasts (bone-resorbing cells). Remodelling of bone begins early in foetal life, and once the skeleton is fully formed in young adults, almost all of the metabolic activity is in this form. Bone is constantly destroyed or resorbed by osteoclasts and then replaced by osteoblasts. Many bone diseases, i.e. osteoporosis, also known as bone loss, typically reflect an imbalance in skeletal turnover. The cyclic adenosine monophosphate (cAMP) and the cyclic guanosine monophosphate (cGMP) are second messengers involved in a variety of cellular responses to such extracellular agents as hormones and neurotransmitters. In the hormonal regulation of bone metabolism, i.e. via parathyroid hormone (PTH), parathyroid hormone-related peptide (PTHrp) and prostaglandin E2 signal via cAMP. cAMP and cGMP are formed by adenylate and guanylate cyclases and are degraded by phosphodiesterases (PDEs). PDEs determine the amplitudes of cyclic nucleotide-mediated hormonal responses and modulate the duration of the signal. The activities of the PDEs are regulated by multiple inputs from other signalling systems and are crucial points of cross-talk between the pathways. Food-derived bioactive peptides are reported to express a variety of functions in vivo. The angiotensin-converting enzymes (ACEs) are involved in the regulation of the specific maturation or degradation of a number of mammalian bioactive peptides. The bioactive peptides offer also a nutriceutical and a nutrigenomic aspect to bone cell biology. The aim of this study was to investigate the influence of PDEs and bioactive peptides on the activation and the differentiation of human osteoblast cells. The profile of PDEs in human osteoblast-like cells and the effect of glucocorticoids on the function of cAMP PDEs, were investigated at the mRNA and enzyme levels. The effects of PDEs on bone formation and osteoblast gene expression were determined with chemical inhibitors and siRNAs (short interfering RNAs). The influence of bioactive peptides on osteoblast gene expression and proliferation was studied at the mRNA and cellular levels. This work provides information on how PDEs are involved in the function and the differentiation of osteoblasts. The findings illustrate that gene-specific silencing with an RNA interference (RNAi) method is useful in inhibiting, the gene expression of specific PDEs and further, PDE7 inhibition upregulates several osteogenic genes and increases bALP activity and mineralization in human mesenchymal stem cells-derived osteoblasts. PDEs appear to be involved in a mechanism by which glucocorticoids affect cAMP signaling. This may provide a potential route in the formation of glucocorticoid-induced bone loss, involving the down-regulation of cAMP-PDE. PDEs may play an important role in the regulation of osteoblastic differentiation. Isoleucine-proline-proline (IPP), a bioactive peptide, possesses the potential to increase osteoblast proliferation, differentiation and signalling.

Neuronal K-Cl Cotransporter : Transcriptional Mechanisms of KCC2 Gene Regulation

K-Cl cotransporter 2 (KCC2) maintains a low intracellular Cl concentration required for fast hyperpolarizing responses of neurons to classical inhibitory neurotransmitters γ-aminobutyric acid (GABA) and glycine. Decreased Cl extrusion observed in genetically modified KCC2-deficient mice leads to depolarizing GABA responses, impaired brain inhibition, and as a consequence to epileptic seizures. Identification of mechanisms regulating activity of the SLC12A5 gene, which encodes the KCC2 cotransporter, in normal and pathological conditions is, thus, of extreme importance. Multiple reports have previously elucidated in details a spatio-temporal pattern of KCC2 expression. Among the characteristic features are an exclusive neuronal specificity, a dramatic upregulation during embryonic and early postnatal development, and a significant downregulation by neuronal trauma. Numerous studies confirmed these expressional features, however transcriptional mechanisms predetermining the SLC12A5 gene behaviour are still unknown. The aim of the presented thesis is to recognize such transcriptional mechanisms and, on their basis, to create a transcriptional model that would explain the established SLC12A5 gene behaviour. Up to recently, only one KCC2 transcript has been thought to exist. A particular novelty of the presented work is the identification of two SLC12A5 gene promoters (SLC12A5-1a and SLC12A5-1b) that produce at least two KCC2 isoforms (KCC2a and KCC2b) differing by their N-terminal parts. Even though a functional 86Rb+ assay reveals no significant difference between transport activities of the isoforms, consensus sites for several protein kinases, found in KCC2a but not in KCC2b, imply a distinct kinetic regulation. As a logical continuation, the current work presents a detailed analysis of the KCC2a and KCC2b expression patterns. This analysis shows an exclusively neuron-specific pattern and similar expression levels for both isoforms during embryonic and neonatal development in rodents. During subsequent postnatal development, the KCC2b expression dramatically increases, while KCC2a expression, depending on central nervous system (CNS) area, either remains at the same level or moderately decreases. In an attempt to explain both the neuronal specificity and the distinct expressional kinetics of the KCC2a and KCC2b isoforms during postnatal development, the corresponding SLC12A5-1a and SLC12A5-1b promoters have been subjected to a comprehensive bioinformatical analysis. Binding sites of several transcription factors (TFs), conserved in the mammalian SLC12A5 gene orthologs, have been identified that might shed light on the observed behaviour of the SLC12A5 gene. Possible roles of these TFs in the regulating of the SLC12A5 gene expression have been elucidated in subsequent experiments and are discussed in the current thesis.

Positional cloning and pathway analysis of the asthma susceptibility gene, NPSR1

In the present study, we identified a novel asthma susceptibility gene, NPSR1 (neuropeptide S receptor 1) on chromosome 7p14.3 by the positional cloning strategy. An earlier significant linkage mapping result among Finnish Kainuu asthma families was confirmed in two independent cohorts: in asthma families from Quebec, Canada and in allergy families from North Karelia, Finland. The linkage region was narrowed down to a 133-kb segment by a hierarchial genotyping method. The observed 77-kb haplotype block showed 7 haplotypes and a similar risk and nonrisk pattern in all three populations studied. All seven haplotypes occur in all three populations at frequences > 2%. Significant elevated relative risks were detected for elevated total IgE (immunoglobulin E) or asthma. Risk effects of the gene variants varied from 1.4 to 2.5. NPSR1 belongs to the G protein-coupled receptor (GPCR) family with a topology of seven transmembrane domains. NPSR1 has 9 exons, with the two main transcripts, A and B, encoding proteins of 371 and 377 amino acids, respectively. We detected a low but ubiquitous expression level of NPSR1-B in various tissues and endogenous cell lines while NPSR1-A has a more restricted expression pattern. Both isoforms were expressed in the lung epithelium. We observed aberrant expression levels of NPSR1-B in smooth muscle in asthmatic bronchi as compared to healthy. In an experimental mouse model, the induced lung inflammation resulted in elevated Npsr1 levels. Furthermore, we demonstrated that the activation of NPSR1 with its endogenous agonist, neuropeptide S (NPS), resulted in a significant inhibition of the growth of NPSR1-A overexpressing stable cell lines (NPSR1-A cells). To determine which target genes were regulated by the NPS-NPSR1 pathway, NPSR1-A cells were stimulated with NPS, and differentially expressed genes were identified using the Affymetrix HGU133Plus2 GeneChip. A total of 104 genes were found significantly up-regulated and 42 down-regulated 6 h after NPS administration. The up-regulated genes included many neuronal genes and some putative susceptibility genes for respiratory disorders. By Gene Ontology enrichment analysis, the biological process terms, cell proliferation, morphogenesis and immune response were among the most altered. The expression of four up-regulated genes, matrix metallopeptidase 10 (MMP10), INHBA (activin A), interleukin 8 (IL8) and EPH receptor A2 (EPHA2), were verified and confirmed by quantitative reverse-transcriptase-PCR. In conclusion, we identified a novel asthma susceptibility gene, NPSR1, on chromosome 7p14.3. NPS-NPSR1 represents a novel pathway that regulates cell proliferation and immune responses, and thus may have functional relevance in the pathogenesis of asthma.

Neurological manifestations and molecular genetics of acute intermittent porphyria in North Western Russia

Acute intermittent porphyria (AIP, MIM #176000) is an inherited metabolic disease due to a partial deficiency of the third enzyme, hydroxymethylbilane synthase (HMBS, EC: 4.3.1.8), in the haem biosynthesis. Neurological symptoms during an acute attack, which is the major manifestation of AIP, are variable and relatively rare, but may endanger a patient's life. In the present study, 12 Russian and two Finnish AIP patients with severe neurological manifestations during an acute attack were studied prospectively from 1995 to 2006. Autonomic neuropathy manifested as abdominal pain (88%), tachycardia (94%), hypertension (75%) and constipation (88%). The most common neurological sign was acute motor peripheral neuropathy (PNP, 81%) often associated with neuropathic sensory loss (54%) and CNS involvement (85%). Despite heterogeneity of the neurological manifestations in our patients with acute porphyria, the major pattern of PNP associated with abdominal pain, dysautonomia, CNS involvement and mild hepatopathy could be demonstrated. If more strict inclusion criteria for biochemical abnormalities (>10-fold increase in excretion of urinary PBG) are applied, neurological manifestations in an acute attack are probably more homogeneous than described previously, which suggests that some of the neurological patients described previously may not have acute porphyria but rather secondary porphyrinuria. Screening for acute porphyria using urinary PBG is useful in a selected group of neurological patients with acute PNP or encephalopathy and seizures associated with pain and dysautonomia. Clinical manifestations and the outcome of acute attacks were used as a basis for developing a 30-score scale of the severity of an acute attack. This scale can easily be used in clinical practice and to standardise the outcome of an attack. Degree of muscle weakness scored by MRC, prolonged mechanical ventilation, bulbar paralysis, impairment of consciousness and hyponatraemia were important signs of a poor prognosis. Arrhythmia was less important and autonomic dysfunction, severity of pain and mental symptoms did not affect the outcome. The delay in the diagnosis and repeated administrations of precipitating factors were the main cause of proceeding of an acute attack into pareses and severe CNS involvement and a fatal outcome in two patients. Nerve conduction studies and needle EMG were performed in eleven AIP patients during an acute attack and/or in remission. Nine patients had severe PNP and two patients had an acute encephalopathy but no clinically evident PNP. In addition to axonopathy, features suggestive of demyelination could be demonstrated in patients with severe PNP during an acute attack. PNP with a moderate muscle weakness was mainly pure axonal. Sensory involvement was common in acute PNP and could be subclinical. Decreased conduction velocities with normal amplitudes of evoked potentials during acute attacks with no clinically evident PNP indicated subclinical polyneuropathy. Reversible symmetrical lesions comparable with posterior reversible encephalopathy syndrome (PRES) were revealed in two patients' brain CT or MRI during an acute attack. In other five patients brain MRI during or soon after the symptoms was normal. The frequency of reversible brain oedema in AIP is probably under-estimated since it may be short-lasting and often indistinguishable on CT or MRI. In the present study, nine different mutations were identified in the HMBS gene in 11 unrelated Russian AIP patients from North Western Russia and their 32 relatives. AIP was diagnosed in nine symptom-free relatives. The majority of the mutations were family-specific and confirmed allelic heterogeneity also among Russian AIP patients. Three mutations, c.825+5G>C, c.825+3_825+6del and c.770T>C, were novel. Six mutations, c.77G>A (p.R26H), c.517C>T (p.R173W), c.583C>T (p.R195C), c.673C>T (p.R225X), c.739T>C (p.C247R) and c.748G>C (p.E250A), have previously been identified in AIP patients from Western and other Eastern European populations. The effects of novel mutations were studied by amplification and sequencing of the reverse-transcribed total RNA obtained from the patients' lymphoblastoid or fibroblast cell lines. The mutations c.825+5G>C and c.770T>C resulted in varyable amounts of abnormal transcripts, r.822_825del (p.C275fsX2) and [r.770u>c, r.652_771del, r.613_771del (p.L257P, p.G218_L257del, p.I205_L257del)]. All mutations demonstrated low residual activities (0.1-1.3 %) when expressed in COS-1 cells confirming the causality of the mutations and the enzymatic defect of the disease. The clinical outcome, prognosis and correlation between the HMBS genotype and phenotype were studied in 143 Finnish and Russian AIP patients with ten mutations (c.33G>T, c.97delA, InsAlu333, p.R149X, p.R167W, p.R173W, p.R173Q, p.R225G, p.R225X, c.1073delA) and more than six patients in each group. The patients were selected from the pool of 287 Finnish AIP patients presented in a Finnish Porphyria Register (1966-2003) and 23 Russian AIP patients (diagnosed 1995-2003). Patients with the p.R167W and p.R225G mutations showed lower penetrance (19% and 11%) and the recurrence rate (33% and 0%) in comparison to the patients with other mutations (range 36 to 67% and 0 to 66%, respectively), as well as milder biochemical abnormalities [urinary porphobilinogen 47±10 vs. 163±21 mol/L, p<0.001; uroporphyrin 130±40 vs. 942±183 nmol/L, p<0.001] suggesting a milder form of AIP in these patients. Erythrocyte HMBS activity did not correlate with the porphobilinogen excretion in remission or the clinical of the disease. In all AIP severity patients, normal PBG excretion predicted freedom from acute attacks. Urinary PBG excretion together with gender, age at the time of diagnosis and mutation type could predict the likelihood of acute attacks in AIP patients.

Characterization of LRRTM and NGR gene families: Expression and functions

Some leucine-rich repeat (LRR) -containing membrane proteins are known regulators of neuronal growth and synapse formation. In this work I characterize two gene families encoding neuronal LRR membrane proteins, namely the LRRTM (leucine-rich repeat, transmembrane neuronal) and NGR (Nogo-66 receptor) families. I studied LRRTM and NGR family member's mRNA tissue distribution by RT-PCR and by in situ hybridization. Subcellular localization of LRRTM1 protein was studied in neurons and in non-neuronal cells. I discovered that LRRTM and NGR family mRNAs are predominantly expressed in the nervous system, and that each gene possesses a specific expression pattern. I also established that LRRTM and NGR family mRNAs are expressed by neurons, and not by glial cells. Within neurons, LRRTM1 protein is not transported to the plasma membrane; rather it localizes to endoplasmic reticulum. Nogo-A (RTN4), MAG, and OMgp are myelin-associated proteins that bind to NgR1 to limit axonal regeneration after central nervous system injury. To better understand the functions of NgR2 and NgR3, and to explore the possible redundancy in the signaling of myelin inhibitors of neurite growth, I mapped the interactions between NgR family and the known and candidate NgR1 ligands. I identified high-affinity interactions between RTN2-66, RTN3-66 and NgR1. I also demonstrate that Rtn3 mRNA is expressed in the same glial cell population of mouse spinal cord white matter as Nogo-A mRNA, and thus it could have a role in myelin inhibition of axonal growth. To understand how NgR1 interacts with multiple structurally divergent ligands, I aimed first to map in more detail the nature of Nogo-A:NgR1 interactions, and then to systematically map the binding sites of multiple myelin ligands in NgR1 by using a library of NgR1 expression constructs encoding proteins with one or multiple surface residues mutated to alanine. My analysis of the Nogo-A:NgR1 -interactions revealed a novel interaction site between the proteins, suggesting a trivalent Nogo-A:NgR1-interaction. Our analysis also defined a central binding region on the concave side of NgR1's LRR domain that is required for the binding of all known ligands, and a surrounding region critical for binding MAG and OMgp. To better understand the biological role of LRRTMs, I generated Lrrtm1 and Lrrtm3 knock out mice. I show here that reporter genes expressed from the targeted loci can be used for maping the neuronal connections of Lrrtm1 and Lrrtm3 expressing neurons in finer detail. With regard to LRRTM1's role in humans, we found a strong association between a 70 kb-spanning haplotype in the proposed promoter region of LRRTM1 gene and two possibly related phenotypes: left-handedness and schizophrenia. Interestingly, the responsible haplotype was linked to phenotypic variability only when paternally inherited. In summary, I identified two families of neuronal receptor-like proteins, and mapped their expression and certain protein-protein interactions. The identification of a central binding region in NgR1 shared by multiple ligands may facilitate the design and development of small molecule therapeutics blocking binding of all NgR1 ligands. Additionally, the genetic association data suggests that allelic variation upstream of LRRTM1 may play a role in the development of left-right brain asymmetry in humans. Lrrtm1 and Lrrtm3 knock out mice developed as a part of this study will likely be useful for schizophrenia and Alzheimer s disease research.