Genetics of Cardiovascular Disease: A Candidate Gene Study of USF1

Cardiovascular diseases (CVD) are major contributors to morbidity and mortality worldwide. Several interacting environmental, biochemical, and genetic risk factors can increase disease susceptibility. While some of the genes involved in the etiology of CVD are known, many are yet to be discovered. During the last few decades, scientists have searched for these genes with genome-wide linkage and association methods, and with more targeted candidate gene studies. This thesis investigates variation within the upstream transcription factor 1 (USF1) gene locus in relation to CVD risk factors, atherosclerosis, and incidence and prevalence of CVD. This candidate gene was first identified in Finnish families ascertained for familial combined hyperlipidemia, a common dyslipidemia predisposing to coronary heart disease. The gene is a ubiquitously expressed transcription factor regulating expression of several genes from lipid and glucose metabolism, inflammation, and endothelial function. First, we examined association between USF1 variants and several CVD risk factors, such as lipid phenotypes, body composition measures, and metabolic syndrome, in two prospective population cohorts. Our data suggested that USF1 contributes to these CVD risk factors at the population level. Notably, the associations with quantitative measurements were mostly detected among study subjects with CVD or metabolic syndrome, suggesting complex interactions between USF1 effects and the pathophysiological state of an individual. Second, we investigated how variation at the USF1 locus contributes to atherosclerotic lesions of the coronary arteries and abdominal aorta. For this, we used two study samples of middle-aged men with detailed measurements of atherosclerosis obtained in autopsy. USF1 variation significantly associated with areas of several types of lesions, especially with calcification of the arteries. Next, we tested what effect the USF1 risk variants have on sudden cardiac death and incidence of CVD. The atherosclerosis-associated risk variant increased the risk of sudden cardiac death of the same study subjects. Furthermore, USF1 alleles associated with incidence of CVD in the Finnish population follow-up cohorts. These associations were especially prominent among women, suggesting a sex specific effect, which has also been detected in subsequent studies. Finally, as some of the low-yield DNA samples of the Finnish follow-up study cohort needed to be whole-genome amplified (WGA) prior to genotyping, we evaluated whether the produced WGA genotypes were of good quality. Although the samples giving genotype discrepancies could not be detected before genotyping with standard laboratory quality control methods, our results suggested that enhanced quality control at the time of the genotyping could identify such samples. In addition, combining two WGA reactions into one pooled DNA sample for genotyping markedly reduced the number of discrepancies and samples showing them. In conclusion, USF1 seems to have a role in the etiology of CVD. Additional studies are warranted to identify functional variants and to study interactions between USF1 and other genetic or environmental factors. This USF1 study, and other studies with low DNA yield of some samples, can benefit from whole genome amplification of the low-yield samples prior to genotyping. Careful quality control procedures are, however, needed in WGA genotyping.

Identification of the semaphorin receptor PLXNA2 as a candidate for susceptibility to schizophrenia

The discovery of genetic factors that contribute to schizophrenia susceptibility is a key challenge in understanding the etiology of this disease. Here, we report the identification of a novel schizophrenia candidate gene on chromosome 1q32, plexin A2 (PLXNA2), in a genome-wide association study using 320 patients with schizophrenia of European descent and 325 matched controls. Over 25,000 single-nucleotide polymorphisms (SNPs) located within approximately 14,000 genes were tested. Out of 62 markers found to be associated with disease status, the most consistent finding was observed for a candidate locus on chromosome 1q32. The marker SNP rs752016 showed suggestive association with schizophrenia (odds ratio (OR) = 1.49, P = 0.006). This result was confirmed in an independent case-control sample of European Americans (combined OR = 1.38, P = 0.035) and similar genetic effects were observed in smaller subsets of Latin Americans (OR = 1.26) and Asian Americans (OR = 1.37). Supporting evidence was also obtained from two family-based collections, one of which reached statistical significance (OR = 2.2, P = 0.02). High-density SNP mapping showed that the region of association spans approximately 60 kb of the PLXNA2 gene. Eight out of 14 SNPs genotyped showed statistically significant differences between cases and controls. These results are in accordance with previous genetic findings that identified chromosome 1q32 as a candidate region for schizophrenia. PLXNA2 is a member of the transmembrane semaphorin receptor family that is involved in axonal guidance during development and may modulate neuronal plasticity and regeneration. The PLXNA2 ligand semaphorin 3A has been shown to be upregulated in the cerebellum of individuals with schizophrenia. These observations, together with the genetic results, make PLXNA2 a likely candidate for the 1q32 schizophrenia susceptibility locus.

Rhipicephalus (Boophilus) microplus tick in vitro feeding methods for functional (dsRNA) and vaccine candidate (antibody) screening

Rhipicephalus (Boophilus) microplus (Acari: Ixodidae) ticks cause economic losses for cattle industries throughout tropical and subtropical regions of the world estimated at $US2.5 billion annually. Lack of access to efficacious long-lasting vaccination regimes and increases in tick acaricide resistance have led to the investigation of targets for the development of novel tick vaccines and treatments. In vitro tick feeding has been used for many tick species to study the effect of new acaricides on the transmission of tick-borne pathogens. Few studies have reported the use of in vitro feeding for functional genomic studies using RNA interference and/or the effect of specific anti-tick antibodies. In particular, in vitro feeding reports for the cattle tick are limited due to its relatively short hypostome. Previously published methods were further modified to broaden optimal tick sizes/weights, feeding sources including bovine and ovine serum, optimisation of commercially available blood anti-coagulant tubes, and IgG concentrations for effective antibody delivery. Ticks are fed overnight and monitored for ∼5–6 weeks to determine egg output and success of larval emergence using a humidified incubator. Lithium-heparin blood tubes provided the most reliable anti-coagulant for bovine blood feeding compared with commercial citrated (CPDA) and EDTA tubes. Although >30 mg semi-engorged ticks fed more reliably, ticks as small as 15 mg also fed to repletion to lay viable eggs. Ticks which gained less than ∼10 mg during in vitro feeding typically did not lay eggs. One mg/ml IgG from Bm86-vaccinated cattle produced a potent anti-tick effect in vitro (83% efficacy) similar to that observed in vivo. Alternatively, feeding of dsRNA targeting Bm86 did not demonstrate anti-tick effects (11% efficacy) compared with the potent effects of ubiquitin dsRNA. This study optimises R. microplus tick in vitro feeding methods which support the development of cattle tick vaccines and treatments.

Evaluation of potential candidate genes involved in salinity tolerance in striped catfish (Pangasianodon hypophthalmus) using an RNA-Seq approach

Increasing salinity levels in freshwater and coastal environments caused by sea level rise linked to climate change is now recognized to be a major factor that can impact fish growth negatively, especially for freshwater teleost species. Striped catfish (Pangasianodon hypophthalmus) is an important freshwater teleost that is now widely farmed across the Mekong River Delta in Vietnam. Understanding the basis for tolerance and adaptation to raised environmental salinity conditions can assist the regional culture industry to mitigate predicted impacts of climate change across this region. Attempt of next generation sequencing using the ion proton platform results in more than 174 million raw reads from three tissue libraries (gill, kidney and intestine). Reads were filtered and de novo assembled using a variety of assemblers and then clustered together to generate a combined reference transcriptome. Downstream analysis resulted in a final reference transcriptome that contained 60,585 transcripts with an N50 of 683 bp. This resource was further annotated using a variety of bioinformatics databases, followed by differential gene expression analysis that resulted in 3062 transcripts that were differentially expressed in catfish samples raised under two experimental conditions (0 and 15 ppt). A number of transcripts with a potential role in salinity tolerance were then classified into six different functional gene categories based on their gene ontology assignments. These included; energy metabolism, ion transportation, detoxification, signal transduction, structural organization and detoxification. Finally, we combined the data on functional salinity tolerance genes into a hypothetical schematic model that attempted to describe potential relationships and interactions among target genes to explain the molecular pathways that control adaptive salinity responses in P. hypophthalmus. Our results indicate that P. hypophthalmus exhibit predictable plastic regulatory responses to elevated salinity by means of characteristic gene expression patterns, providing numerous candidate genes for future investigations.

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.

Prediction of Candidate Primary Immunodeficiency Disease Genes Using a Support Vector Machine Learning Approach

Screening and early identification of primary immunodeficiency disease (PID) genes is a major challenge for physicians. Many resources have catalogued molecular alterations in known PID genes along with their associated clinical and immunological phenotypes. However, these resources do not assist in identifying candidate PID genes. We have recently developed a platform designated Resource of Asian PDIs, which hosts information pertaining to molecular alterations, protein-protein interaction networks, mouse studies and microarray gene expression profiling of all known PID genes. Using this resource as a discovery tool, we describe the development of an algorithm for prediction of candidate PID genes. Using a support vector machine learning approach, we have predicted 1442 candidate PID genes using 69 binary features of 148 known PID genes and 3162 non-PID genes as a training data set. The power of this approach is illustrated by the fact that six of the predicted genes have recently been experimentally confirmed to be PID genes. The remaining genes in this predicted data set represent attractive candidates for testing in patients where the etiology cannot be ascribed to any of the known PID genes.

Complex trait of acute pancreatitis : Studies of candidate genes and adipokines

Acute pancreatitis (AP) is a common disease. Mild disease resolves spontaneously in a few days. Severe forms of the disease can lead to local complications, necrosis, and abscesses in and around the pancreas. Systemic inflammation in severe AP is associated with distant organ failures. The aim of this study is to identify genetically determined prognostic factors involved in the clinical features of AP. The study employs a candidate-gene approach, and the genes are involved in trysinogen activation in the initiation phase of the disease, as well as in the systemic inflammation as the disease proceeds. The last study examines adipokines, fat-derived hormones characterized with the capacity to modify inflammation. SPINK 1 is a gene coding trypsin activation inhibitor. Mutations N34S and P55N were determined by minisequencing methods in 371 AP patients and in 459 controls. The mutation N34S was more common in AP patients (7.8%) than in controls (2.6%). This suggests that SPINK 1 gene mutation N34S is a risk factor for AP. In the fourth study, in 12 matched pairs of patients with severe and mild AP, levels of adipokines, adiponectin, and leptin were evaluated. Plasma adipokine levels did not differ between patients with mild and severe AP. The results suggest that in AP, adipokine plasma levels are not factors predisposing to organ failures. This study identified the SPINK 1 mutation N34S to be a risk factor for AP in the general population. As AP is a multifactorial disease, and extensive genetic heterogeneity is likely, further identification of genetic factors in the disease requires larger future studies with more advanced genetic study models. Further identification of the patient characteristics associated with organ failures offers another direction of the study to achieve more detailed understanding of the severe form of AP.

On gravitational fields of stars in f(R) gravity theories

Einstein's general relativity is a classical theory of gravitation: it is a postulate on the coupling between the four-dimensional, continuos spacetime and the matter fields in the universe, and it yields their dynamical evolution. It is believed that general relativity must be replaced by a quantum theory of gravity at least at extremely high energies of the early universe and at regions of strong curvature of spacetime, cf. black holes. Various attempts to quantize gravity, including conceptually new models such as string theory, have suggested that modification to general relativity might show up even at lower energy scales. On the other hand, also the late time acceleration of the expansion of the universe, known as the dark energy problem, might originate from new gravitational physics. Thus, although there has been no direct experimental evidence contradicting general relativity so far - on the contrary, it has passed a variety of observational tests - it is a question worth asking, why should the effective theory of gravity be of the exact form of general relativity? If general relativity is modified, how do the predictions of the theory change? Furthermore, how far can we go with the changes before we are face with contradictions with the experiments? Along with the changes, could there be new phenomena, which we could measure to find hints of the form of the quantum theory of gravity? This thesis is on a class of modified gravity theories called f(R) models, and in particular on the effects of changing the theory of gravity on stellar solutions. It is discussed how experimental constraints from the measurements in the Solar System restrict the form of f(R) theories. Moreover, it is shown that models, which do not differ from general relativity at the weak field scale of the Solar System, can produce very different predictions for dense stars like neutron stars. Due to the nature of f(R) models, the role of independent connection of the spacetime is emphasized throughout the thesis.

Molecular line and continuum studies of the early stages of star formation

New stars form in dense interstellar clouds of gas and dust called molecular clouds. The actual sites where the process of star formation takes place are the dense clumps and cores deeply embedded in molecular clouds. The details of the star formation process are complex and not completely understood. Thus, determining the physical and chemical properties of molecular cloud cores is necessary for a better understanding of how stars are formed. Some of the main features of the origin of low-mass stars, like the Sun, are already relatively well-known, though many details of the process are still under debate. The mechanism through which high-mass stars form, on the other hand, is poorly understood. Although it is likely that the formation of high-mass stars shares many properties similar to those of low-mass stars, the very first steps of the evolutionary sequence are unclear. Observational studies of star formation are carried out particularly at infrared, submillimetre, millimetre, and radio wavelengths. Much of our knowledge about the early stages of star formation in our Milky Way galaxy is obtained through molecular spectral line and dust continuum observations. The continuum emission of cold dust is one of the best tracers of the column density of molecular hydrogen, the main constituent of molecular clouds. Consequently, dust continuum observations provide a powerful tool to map large portions across molecular clouds, and to identify the dense star-forming sites within them. Molecular line observations, on the other hand, provide information on the gas kinematics and temperature. Together, these two observational tools provide an efficient way to study the dense interstellar gas and the associated dust that form new stars. The properties of highly obscured young stars can be further examined through radio continuum observations at centimetre wavelengths. For example, radio continuum emission carries useful information on conditions in the protostar+disk interaction region where protostellar jets are launched. In this PhD thesis, we study the physical and chemical properties of dense clumps and cores in both low- and high-mass star-forming regions. The sources are mainly studied in a statistical sense, but also in more detail. In this way, we are able to examine the general characteristics of the early stages of star formation, cloud properties on large scales (such as fragmentation), and some of the initial conditions of the collapse process that leads to the formation of a star. The studies presented in this thesis are mainly based on molecular line and dust continuum observations. These are combined with archival observations at infrared wavelengths in order to study the protostellar content of the cloud cores. In addition, centimetre radio continuum emission from young stellar objects (YSOs; i.e., protostars and pre-main sequence stars) is studied in this thesis to determine their evolutionary stages. The main results of this thesis are as follows: i) filamentary and sheet-like molecular cloud structures, such as infrared dark clouds (IRDCs), are likely to be caused by supersonic turbulence but their fragmentation at the scale of cores could be due to gravo-thermal instability; ii) the core evolution in the Orion B9 star-forming region appears to be dynamic and the role played by slow ambipolar diffusion in the formation and collapse of the cores may not be significant; iii) the study of the R CrA star-forming region suggests that the centimetre radio emission properties of a YSO are likely to change with its evolutionary stage; iv) the IRDC G304.74+01.32 contains candidate high-mass starless cores which may represent the very first steps of high-mass star and star cluster formation; v) SiO outflow signatures are seen in several high-mass star-forming regions which suggest that high-mass stars form in a similar way as their low-mass counterparts, i.e., via disk accretion. The results presented in this thesis provide constraints on the initial conditions and early stages of both low- and high-mass star formation. In particular, this thesis presents several observational results on the early stages of clustered star formation, which is the dominant mode of star formation in our Galaxy.

Ultra high throughput sequencing excludes MDH1 as candidate gene for RP28-linked retinitis pigmentosa

Purpose: Mutations in IDH3B, an enzyme participating in the Krebs cycle, have recently been found to cause autosomal recessive retinitis pigmentosa (arRP). The MDH1 gene maps within the RP28 arRP linkage interval and encodes cytoplasmic malate dehydrogenase, an enzyme functionally related to IDH3B. As a proof of concept for candidate gene screening to be routinely performed by ultra high throughput sequencing (UHTs), we analyzed MDH1 in a patient from each of the two families described so far to show linkage between arRP and RP28. Methods: With genomic long-range PCR, we amplified all introns and exons of the MDH1 gene (23.4 kb). PCR products were then sequenced by short-read UHTs with no further processing. Computer-based mapping of the reads and mutation detection were performed by three independent software packages. Results: Despite the intrinsic complexity of human genome sequences, reads were easily mapped and analyzed, and all algorithms used provided the same results. The two patients were homozygous for all DNA variants identified in the region, which confirms previous linkage and homozygosity mapping results, but had different haplotypes, indicating genetic or allelic heterogeneity. None of the DNA changes detected could be associated with the disease. Conclusions: The MDH1 gene is not the cause of RP28-linked arRP. Our experimental strategy shows that long-range genomic PCR followed by UHTs provides an excellent system to perform a thorough screening of candidate genes for hereditary retinal degeneration.

Ultra high throughput sequencing excludes MDH1 as candidate gene for RP28-linked retinitis pigmentosa

Purpose: Mutations in IDH3B, an enzyme participating in the Krebs cycle, have recently been found to cause autosomal recessive retinitis pigmentosa (arRP). The MDH1 gene maps within the RP28 arRP linkage interval and encodes cytoplasmic malate dehydrogenase, an enzyme functionally related to IDH3B. As a proof of concept for candidate gene screening to be routinely performed by ultra high throughput sequencing (UHTs), we analyzed MDH1 in a patient from each of the two families described so far to show linkage between arRP and RP28. Methods: With genomic long-range PCR, we amplified all introns and exons of the MDH1 gene (23.4 kb). PCR products were then sequenced by short-read UHTs with no further processing. Computer-based mapping of the reads and mutation detection were performed by three independent software packages. Results: Despite the intrinsic complexity of human genome sequences, reads were easily mapped and analyzed, and all algorithms used provided the same results. The two patients were homozygous for all DNA variants identified in the region, which confirms previous linkage and homozygosity mapping results, but had different haplotypes, indicating genetic or allelic heterogeneity. None of the DNA changes detected could be associated with the disease.

Heat Shock Protein 90 as a Drug Target against Protozoan Infections Biochemical characterization of HSP90 From plasmodium falciparum and trypanosoma evansi and evaluation of its inhibitor as a candidate drug

Using a pharmacological inhibitor of Hsp90 in cultured malarial parasite, we have previously implicated Plasmodium falciparum Hsp90 (PfHsp90) as a drug target against malaria. In this study, we have biochemically characterized PfHsp90 in terms of its ATPase activity and interaction with its inhibitor geldanamycin (GA) and evaluated its potential as a drug target in a preclinical mouse model of malaria. In addition, we have explored the potential of Hsp90 inhibitors as drugs for the treatment of Trypanosoma infection in animals. Our studies with full-length PfHsp90 showed it to have the highest ATPase activity of all known Hsp90s; its ATPase activity was 6 times higher than that of human Hsp90. Also, GA brought about more robust inhibition of PfHsp90 ATPase activity as compared with human Hsp90. Mass spectrometric analysis of PfHsp90 expressed in P. falciparum identified a site of acetylation that overlapped with Aha1 and p23 binding domain, suggesting its role in modulating Hsp90 multichaperone complex assembly. Indeed, treatment of P. falciparum cultures with a histone deacetylase inhibitor resulted in a partial dissociation of PfHsp90 complex. Furthermore, we found a well known, semisynthetic Hsp90 inhibitor, namely 17-(allylamino)-17-demethoxygeldanamycin, to be effective in attenuating parasite growth and prolonging survival in a mouse model of malaria. We also characterized GA binding to Hsp90 from another protozoan parasite, namely Trypanosoma evansi. We found 17-(allylamino)-17-demethoxygeldanamycin to potently inhibit T. evansi growth in a mouse model of trypanosomiasis. In all, our biochemical characterization, drug interaction, and animal studies supported Hsp90 as a drug target and its inhibitor as a potential drug against protozoan diseases.

3D MHD simulations of subsurface convection in OB stars

During their main sequence evolution, massive stars can develop convective regions very close to their surface. These regions are caused by an opacity peak associated with iron ionization. Cantiello et al. (2009) found a possible connection between the presence of sub-photospheric convective motions and small scale stochastic velocities in the photosphere of early-type stars. This supports a physical mechanism where microturbulence is caused by waves that are triggered by subsurface convection zones. They further suggest that clumping in the inner parts of the winds of OB stars could be related to subsurface convection, and that the convective layers may also be responsible for stochastic excitation of non-radial pulsations. Furthermore, magnetic fields produced in the iron convection zone could appear at the surface of such massive stars. Therefore subsurface convection could be responsible for the occurrence of observable phenomena such as line profile variability and discrete absorption components. These phenomena have been observed for decades, but still evade a clear theoretical explanation. Here we present preliminary results from 3D MHD simulations of such subsurface convection.