Genetic Sharing with Cardiovascular Disease Risk Factors and Diabetes Reveals Novel Bone Mineral Density Loci

Bone Mineral Density (BMD) is a highly heritable trait, but genome-wide association studies have identified few genetic risk factors. Epidemiological studies suggest associations between BMD and several traits and diseases, but the nature of the suggestive comorbidity is still unknown. We used a novel genetic pleiotropy-informed conditional False Discovery Rate (FDR) method to identify single nucleotide polymorphisms (SNPs) associated with BMD by leveraging cardiovascular disease (CVD) associated disorders and metabolic traits. By conditioning on SNPs associated with the CVD-related phenotypes, type 1 diabetes, type 2 diabetes, systolic blood pressure, diastolic blood pressure, high density lipoprotein, low density lipoprotein, triglycerides and waist hip ratio, we identified 65 novel independent BMD loci (26 with femoral neck BMD and 47 with lumbar spine BMD) at conditional FDR < 0.01. Many of the loci were confirmed in genetic expression studies. Genes validated at the mRNA levels were characteristic for the osteoblast/osteocyte lineage, Wnt signaling pathway and bone metabolism. The results provide new insight into genetic mechanisms of variability in BMD, and a better understanding of the genetic underpinnings of clinical comorbidity.

Detection of HPV in mouth floor squamous cell carcinoma and its correlation with clinicopathologic variables, risk factors and survival

The human papillomavirus (HPV) has been historically associated with head and neck cancers, although its role in oral carcinogenesis remains poorly defined. The purpose of this study was to investigate the prevalence of HPV in mouth floor squamous cell carcinoma and correlate it with clinicopathologic variables, risk factors and survival. HPV presence was evaluated by nested polymerase chain reaction (nPCR) in 29 paraffin-embedded specimens of mouth floor squamous cell carcinoma. HPV DNA was detected in 17.2% (5 of 29) of the specimens; the highest prevalence was observed in non-smoking patients over the age of 60 years. All HPV DNA positive specimens were detected in men with clinical stage III and IV lesions, being most of which were moderately differentiated. Despite this correlation there were no statistically significant differences observed among the analyzed variables, including patients' survival. The relatively low incidence of HPV DNA present in these tumors suggests that this virus does not, by itself, have a significant role in the development of mouth floor squamous cell carcinoma. J Oral Pathol Med (2008) 37: 593-598

Osteoporosis case finding in the general practice: phalangeal radiographic absorptiometry with and without risk factors for osteoporosis to select postmenopausal women eligible for lumbar spine and hip densitometry

Mass screening for osteoporosis using DXA measurements at the spine and hip is presently not recommended by health authorities. Instead, risk factor questionnaires and peripheral bone measurements may facilitate the selection of women eligible for axial bone densitometry. The aim of this study was to validate a case finding strategy for postmenopausal women who would benefit most from subsequent DXA measurement by using phalangeal radiographic absorptiometry (RA) alone or in combination with risk factors in a general practice setting. The sensitivity and specificity of this strategy in detecting osteoporosis (T-score < or =2.5 SD at the spine and/or the hip) were compared with those of the current reimbursement criteria for DXA measurements in Switzerland. Four hundred and twenty-three postmenopausal women with one or more risk factors for osteoporosis were recruited by 90 primary care physicians who also performed the phalangeal RA measurements. All women underwent subsequent DXA measurement of the spine and the hip at the Osteoporosis Policlinic of the University Hospital of Berne. They were allocated to one of two groups depending on whether they matched with the Swiss reimbursement conditions for DXA measurement or not. Logistic regression models were used to predict the likelihood of osteoporosis versus "no osteoporosis" and to derive ROC curves for the various strategies. Differences in the areas under the ROC curves (AUC) were tested for significance. In women lacking reimbursement criteria, RA achieved a significantly larger AUC (0.81; 95% CI 0.72-0.89) than the risk factors associated with patients' age, height and weight (0.71; 95% C.I. 0.62-0.80). Furthermore, in this study, RA provided a better sensitivity and specificity in identifying women with underlying osteoporosis than the currently accepted criteria for reimbursement of DXA measurement. In the Swiss environment, RA is a valid case finding tool for patients with risk factors for osteoporosis, especially for those who do not qualify for DXA reimbursement.

Spatiotemporal dynamics of Puumala hantavirus associated with its rodent host, Myodes glareolus

Many viruses significantly impact human and animal health. Understanding the population dynamics of these viruses and their hosts can provide important insights for epidemiology and virus evolution. Puumala virus (PUUV) is a European hantavirus that may cause regional outbreaks of hemorrhagic fever with renal syndrome in humans. Here, we analyzed the spatiotemporal dynamics of PUUV circulating in local populations of its rodent reservoir host, the bank vole (Myodes glareolus) during eight years. Phylogenetic and population genetic analyses of all three genome segments of PUUV showed strong geographical structuring at a very local scale. There was a high temporal turnover of virus strains in the local bank vole populations, but several virus strains persisted through multiple years. Phylodynamic analyses showed no significant changes in the local effective population sizes of PUUV, although vole numbers and virus prevalence fluctuated widely. Microsatellite data demonstrated also a temporally persisting subdivision between local vole populations, but these groups did not correspond to the subdivision in the virus strains. We conclude that restricted transmission between vole populations and genetic drift play important roles in shaping the genetic structure and temporal dynamics of PUUV in its natural host which has several implications for zoonotic risks of the human population.

Ross River virus disease in Australia, 1886-1998, with analysis of risk factors associated with outbreaks

Ross River virus (RE) is a mosquito-borne arbovirus responsible for outbreaks of polyarthritic disease throughout Australia. To better understand human and environmental factors driving such events, 57 historical reports oil RR Outbreaks between 1896 and 1998 were examined collectively. The magnitude, regularity, seasonality, and locality of outbreaks were found to be wide ranging; however, analysis of climatic and tidal data highlighted that environmental conditions let differently ill tropical, arid, and temperate regions. Overall, rainfall seems to be the single most important risk factor, with over 90% of major outbreak locations receiving higher than average rainfall in preceding mouths. Many temperatures were close to average, particularly in tropical populations; however, in arid regions, below average maximum temperatures predominated, and ill southeast temperate regions, above average minimum temperatures predominated. High spring tides preceded coastal Outbreaks, both in the presence and absence of rainfall, and the relationship between rainfall and the Southern Oscillation Index and Lit Nina episodes suggest they may be useful predictive tools, but only ill southeast temperate regions. Such heterogeneity predisposing outbreaks supports the notion that there are different RE epidemiologies throughout Australia but also Suggests that generic parameters for the prediction and control of outbreaks are of limited use at a local level.

The Structure and Functions of Hantavirus Nucleocapsid Protein

Hantaviruses, members of the genus Hantavirus in the Bunyaviridae family, are enveloped single-stranded RNA viruses with tri-segmented genome of negative polarity. In humans, hantaviruses cause two diseases, hemorrhagic fever with renal syndrome (HFRS) and hantavirus pulmonary syndrome (HPS), which vary in severity depending on the causative agent. Each hantavirus is carried by a specific rodent host and is transmitted to humans through excreta of infected rodents. The genome of hantaviruses encodes four structural proteins: the nucleocapsid protein (N), the glycoproteins (Gn and Gc), and the polymerase (L) and also the nonstructural protein (NSs). This thesis deals with the functional characterization of hantavirus N protein with regard to its structure. Structural studies of the N protein have progressed slowly and the crystal structure of the whole protein is still not available, therefore biochemical assays coupled with bioinformatical modeling proved essential for studying N protein structure and functions. Presumably, during RNA encapsidation, the N protein first forms intermediate trimers and then oligomers. First, we investigated the role of N-terminal domain in the N protein oligomerization. The results suggested that the N-terminal region of the N protein forms a coiled-coil, in which two antiparallel alpha helices interact via their hydrophobic seams. Hydrophobic residues L4, I11, L18, L25 and V32 in the first helix and L44, V51, L58 and L65 in the second helix were crucial for stabilizing the structure. The results were consistent with the head-to-head, tail-to-tail model for hantavirus N protein trimerization. We demonstrated that an intact coiled-coil structure of the N terminus is crucial for the oligomerization capacity of the N protein. We also added new details to the head-to-head, tail-to-tail model of trimerization by suggesting that the initial step is based on interaction(s) between intact intra-molecular coiled-coils of the monomers. We further analyzed the importance of charged aa residues located within the coiled-coil for the N protein oligomerization. To predict the interacting surfaces of the monomers we used an upgraded in silico model of the coiled-coil domain that was docked into a trimer. Next the predicted target residues were mutated. The results obtained using the mammalian two-hybrid assay suggested that conserved charged aa residues within the coiled-coil make a substantial contribution to the N protein oligomerization. This contribution probably involves the formation of interacting surfaces of the N monomers and also stabilization of the coiled-coil via intramolecular ionic bridging. We proposed that the tips of the coiled-coils are the first to come into direct contact and thus initiate tight packing of the three monomers into a compact structure. This was in agreement with the previous results showing that an increase in ionic strength abolished the interaction between N protein molecules. We also showed that residues having the strongest effect on the N protein oligomerization are not scattered randomly throughout the coiled-coil 3D model structure, but form clusters. Next we found evidence for the hantaviral N protein interaction with the cytoplasmic tail of the glycoprotein Gn. In order to study this interaction we used the GST pull-down assay in combination with mutagenesis technique. The results demonstrated that intact, properly folded zinc fingers of the Gn protein cytoplasmic tail as well as the middle domain of the N protein (that includes aa residues 80 248 and supposedly carries the RNA-binding domain) are essential for the interaction. Since hantaviruses do not have a matrix protein that mediates the packaging of the viral RNA in other negatve stranded viruses (NSRV), hantaviral RNPs should be involved in a direct interaction with the intraviral domains of the envelope-embedded glycoproteins. By showing the N-Gn interaction we provided the evidence for one of the crucial steps in the virus replication at which RNPs are directed to the site of the virus assembly. Finally we started analysis of the N protein RNA-binding region, which is supposedly located in the middle domain of the N protein molecule. We developed a model for the initial step of RNA-binding by the hantaviral N protein. We hypothesized that the hantaviral N protein possesses two secondary structure elements that initiate the RNA encapsidation. The results suggest that amino acid residues (172-176) presumably act as a hook to catch vRNA and that the positively charged interaction surface (aa residues 144-160) enhances the initial N-RNA interacation. In conclusion, we elucidated new functions of hantavirus N protein. Using in silico modeling we predicted the domain structure of the protein and using experimental techniques showed that each domain is responsible for executing certain function(s). We showed that intact N terminal coiled-coil domain is crucial for oligomerization and charged residues located on its surface form a interaction surface for the N monomers. The middle domain is essential for interaction with the cytoplasmic tail of the Gn protein and RNA binding.

Hantavirus glycoprotein GN in virion assembly and regulation of interferon response

Hantaviruses have a tri-segmented negative-stranded RNA genome. The S segment encodes the nucleocapsid protein (N), M segment two glycoproteins, Gn and Gc, and the L segment the RNA polymerase. Gn and Gc are co-translationally cleaved from a precursor and targeted to the cis-Golgi compartment. The Gn glycoprotein consists of an external domain, a transmembrane domain and a C-terminal cytoplasmic domain. In addition, the S segment of some hantaviruses, including Tula and Puumala virus, have an open reading frame (ORF) encoding a nonstructural potein NSs that can function as a weak interferon antagonist. The mechanisms of hantavirus-induced pathogenesis are not fully understood but it is known that both hemorrhagic fever with renal syndrome (HFRS) and hantavirus (cardio) pulmonary syndrome (HCPS) share various features such as increased capillary permeability, thrombocytopenia and upregulation of TNF-. Several hantaviruses have been reported to induce programmed cell death (apoptosis), such as TULV-infected Vero E6 cells which is known to be defective in interferon signaling. Recently reports describing properties of the hantavirus Gn cytoplasmic tail (Gn-CT) have appeared. The Gn-CT of hantaviruses contains animmunoreceptor tyrosine-based activation motif (ITAM) which directs receptor signaling in immune and endothelial cells; and contain highly conserved classical zinc finger domains which may have a role in the interaction with N protein. More functions of Gn protein have been discovered, but much still remains unknown. Our aim was to study the functions of Gn protein from several aspects: synthesis, degradation and interaction with N protein. Gn protein was reported to inhibit interferon induction and amplication. For this reason, we also carried out projects studying the mechanisms of IFN induction and evasion by hantavirus. We first showed degradation and aggresome formation of the Gn-CT of the apathogenic TULV. It was reported earlier that the degradation of Gn-CT is related to the pathogenicity of hantavirus. We found that the Gn-CT of the apathogenic hantaviruses (TULV, Prospect Hill virus) was degraded through the ubiquitin-proteasome pathway, and TULV Gn-CT formed aggresomes upon treatment with proteasomal inhibitor. Thus the results suggest that degradation and aggregation of the Gn-CT may be a general property of most hantaviruses, unrelated to pathogenicity. Second, we investigated the interaction of TULV N protein and the TULV Gn-CT. The Gn protein is located on the Golgi membrane and its interaction with N protein has been thought to determine the cargo of the hantaviral ribonucleoprotein which is an important step in virus assembly, but direct evidence has not been reported. We found that TULV Gn-CT fused with GST tag expressed in bacteria can pull-down the N protein expressed in mammalian cells; a mutagenesis assay was carried out, in which we found that the zinc finger motif in Gn-CT and RNA-binding motif in N protein are indispensable for the interaction. For the study of mechanisms of IFN induction and evasion by Old World hantavirus, we found that Old World hantaviruses do not produce detectable amounts of dsRNA in infected cells and the 5 -termini of their genomic RNAs are monophosphorylated. DsRNA and tri-phosphorylated RNA are considered to be critical activators of innate immnity response by interacting with PRRs (pattern recognition receptors). We examined systematically the 5´-termini of hantavirus genomic RNAs and the dsRNA production by different species of hantaviruses. We found that no detectable dsRNA was produced in cells infected by the two groups of the old world hantaviruses: Seoul, Dobrava, Saaremaa, Puumala and Tula. We also found that the genomic RNAs of these Old World hantaviruses carry 5´-monophosphate and are unable to trigger interferon induction. The antiviral response is mainly mediated by alpha/beta interferon. Recently the glycoproteins of the pathogenic hantaviruses Sin Nombre and New York-1 viruses were reported to regulate cellular interferon. We found that Gn-CT can inhibit the induction of IFN activation through Toll-like receptor (TLR) and retinoic acid-inducible gene I-like RNA helicases (RLH) pathway and that the inhibition target lies at the level of TANK-binding kinase 1 (TBK-1)/ IKK epislon complex and myeloid differentiation primary response gene (88) (MyD88) / interferon regulatory factor 7 (IRF-7) complex.

Glycoprotein Interactions in the Assembly of Hantaviruses

Hantaviruses are one of the five genera of the vector-borne virus family Bunyaviridae. While other members of the family are transmitted via arthropods, hantaviruses are carried and transmitted by rodents and insectivores. Occasional transmission to humans occurs via inhalation of aerosolized rodent excreta. When transmitted to man hantaviruses cause hemorrhagic fever with renal syndrome (HFRS, in Eurasia, mortality ~10%) and hantavirus cardiopulmonary syndrome (HCPS, in the Americas, mortality ~40%). The single-stranded, negative-sense RNA genome of hantaviruses is in segments S, M and L that respectively encode for nucleocapsid (N), glycoproteins Gn and Gc, and RNA-dependent RNA-polymerase (RdRp or L protein). The genome segments, encapsidated by N protein to form ribonucleoprotein (RNP), are enclosed inside a lipid envelope decorated by spikes formed of Gn and Gc. The focus of this study was to understand the mechanisms and interactions through which the virion is formed and maintained. We observed that when extracted from virions both Gn and Gc favor homo- over hetero-oligomerization. The minimal glycoprotein complexes extracted from virion by detergent were observed, by using ultracentrifugation and gel filtration, to be tetrameric Gn and homodimeric Gc. These results led us to suggest a model where tetrameric Gn complexes are interconnected through homodimeric Gc units to form the grid-like surface architecture described for hantaviruses. This model was found to correlate with the three-dimensional (3D) reconstruction of virion surface created using cryo-electron tomography (cryo-ET). The 3D-density map showed the spike complex formed of Gn and Gc to be 10 nm high and to display a four-fold symmetry with dimensions of 15 nm times 15 nm. This unique square-shaped complex on a roughly round virion creates a hitch for the assembly, since a sphere cannot be broken into rectangles. Thus additional interactions are likely required for the virion assembly. In cryo-ET we observed that the RNP makes occasional contacts to the viral membrane, suggesting an interaction between the spike and RNP. We were able to demonstrate this interaction using various techniques, and showed that both Gn and Gc contribute to the interaction. This led us to suggest that in addition to the interactions between Gn and Gc, also the interaction between spike and RNP is required for assembly. We found galectin-3 binding protein (referred to as 90K) to co-purify with the virions and showed an interaction between 90K and the virion. Analysis of plasma samples taken from patients hospitalized for Puumala virus infection showed increased concentrations of 90K in the acute phase and the increased 90K level was found to correlate with several parameters that reflect the severity of acute HFRS. The results of these studies confirmed, but also challenged some of the dogmas on the structure and assembly of hantaviruses. We confirmed that Gn and RNP do interact, as long assumed. On the other hand we demonstrated that the glycoproteins Gn and Gc exist as homo-oligomers or appear in large hetero-oligomeric complexes, rather than form primarily heterodimers as was previously assumed. This work provided new insight into the structure and assembly of hantaviruses.

Hantavirus infection: Insights into entry, assembly and pathogenesis

Hantaviruses (family Bunyaviridae, genus Hantavirus) are enveloped viruses incorporating a segmented, negative-sense RNA genome. Each hantavirus is carried by its specific host, either a rodent or an insectivore (shrew), in which the infection is asymptomatic and persistent. In humans, hantaviruses cause Hemorrhagic fever with renal syndrome (HFRS) in Eurasia and Hantavirus cardiopulmonary syndrome (HCPS) in the Americas. In Finland, Puumala virus (genus Hantavirus) is the causative agent of NE, a mild form of HFRS. The HFRS-type diseases are often associated with renal failure and proteinuria that might be mechanistically explained by infected kidney tubular cell degeneration in patients. Previously, it has been shown that non-pathogenic hantavirus, Tula virus (TULV), could cause programmed cell death, apoptosis, in cell cultures. This suggested that the infected kidney tubular degeneration could be caused directly by virus replication. In the first paper of this thesis the molecular mechanisms involved in TULV-induced apoptosis was further elucidated. A virus replication-dependent down-regulation of ERK1/2, concomitantly with the induced apoptosis, was identified. In addition, this phenomenon was not restricted to TULV or to non-pathogenic hantaviruses in general since also a pathogenic hantavirus, Seoul virus, could inhibit ERK1/2 activity. Hantaviruses consist of membrane-spanning glycoproteins Gn and Gc, RNA-dependent RNA polymerase (L protein) and nucleocapsid protein N, which encapsidates the viral genome, and thus forms the ribonucleoprotein (RNP). Interaction between the cytoplasmic tails of viral glycoproteins and RNP is assumed to be the only means how viral genetic material is incorporated into infectious virions. In the second paper of this thesis, it was shown by immunoprecipitation that viral glycoproteins and RNP interact in the purified virions. It was further shown that peptides derived from the cytoplasmic tails (CTs) of both Gn and Gc could bind RNP and recombinant N protein. In the fourth paper the cytoplamic tail of Gn but not Gc was shown to interact with genomic RNA. This interaction was probably rather unspecific since binding of Gn-CT with unrelated RNA and even single-stranded DNA were also observed. However, since the RNP consists of both N protein and N protein-encapsidated genomic RNA, it is possible that the viral genome plays a role in packaging of RNPs into virions. On the other hand, the nucleic acid-binding activity of Gn may have importance in the synthesis of viral RNA. Binding sites of Gn-CT with N protein or nucleic acids were also determined by peptide arrays, and they were largely found to overlap. The Gn-CT of hantaviruses contain a conserved zinc finger (ZF) domain with an unknown function. Some viruses need ZFs in entry or post-entry steps of the viral life cycle. Cysteine residues are required for the folding of ZFs by coordinating zinc-ions, and alkylation of these residues can affect virus infectivity. In the third paper, it was shown that purified hantavirions could be inactivated by treatment with cysteine-alkylating reagents, especially N-ethyl maleimide. However, the effect could not be pin-pointed to the ZF of Gn-CT since also other viral proteins reacted with maleimides, and it was, therefore, impossible to exclude the possibility that other cysteines besides those that were essential in the formation of ZF are required for hantavirus infectivity.

Two problems with variational expectation maximisation for time-series models

Variational methods are a key component of the approximate inference and learning toolbox. These methods fill an important middle ground, retaining distributional information about uncertainty in latent variables, unlike maximum a posteriori methods (MAP), and yet generally requiring less computational time than Monte Carlo Markov Chain methods. In particular the variational Expectation Maximisation (vEM) and variational Bayes algorithms, both involving variational optimisation of a free-energy, are widely used in time-series modelling. Here, we investigate the success of vEM in simple probabilistic time-series models. First we consider the inference step of vEM, and show that a consequence of the well-known compactness property of variational inference is a failure to propagate uncertainty in time, thus limiting the usefulness of the retained distributional information. In particular, the uncertainty may appear to be smallest precisely when the approximation is poorest. Second, we consider parameter learning and analytically reveal systematic biases in the parameters found by vEM. Surprisingly, simpler variational approximations (such a mean-field) can lead to less bias than more complicated structured approximations.

A critical value function approach, with an application to persistent time-series

Researchers often rely on the t-statistic to make inference on parameters in statistical models. It is common practice to obtain critical values by simulation techniques. This paper proposes a novel numerical method to obtain an approximately similar test. This test rejects the null hypothesis when the test statistic islarger than a critical value function (CVF) of the data. We illustrate this procedure when regressors are highly persistent, a case in which commonly-used simulation methods encounter dificulties controlling size uniformly. Our approach works satisfactorily, controls size, and yields a test which outperforms the two other known similar tests.