Search for genetic variants influencing human height

Human growth and attained height are determined by a combination of genetic and environmental effects and in modern Western societies > 80% of the observed variation in height is determined by genetic factors. Height is a fundamental human trait that is associated with many socioeconomic and psychosocial factors and health measures, however little is known of the identity of the specific genes that influence height variation in the general population. This thesis work aimed to identify the genetic variants that influence height in the general population by genome-wide linkage analysis utilizing large family samples. The study focused on analysis of three separate sets of families consisting of: 1) 1,417 individuals from 277 Finnish families (FinnHeight), 2) 8,450 individuals from 3,817 families from Australia and Europe (EUHeight) and 3) 9,306 individuals from 3,302 families from the United States (USHeight). The most significant finding in this study was found in the Finnish family sample where we a locus in the chromosomal region 1p21 was linked to adult height. Several regions showed evidence for linkage in the Australian, European and US families with 8q21 and 15q25 being the most significant. The region on 1p21 was followed up with further studies and we were able to show that the collagen 11-alpha-1 gene (COL11A1) residing at this location was associated with adult height. This association was also confirmed in an independent Finnish population cohort (Health 2000) consisting of 6,542 individuals. From this population sample, we estimated that homozygous males and females for this gene variant were 1.1 and 0.6 cm taller than the respective controls. In this thesis work we identified a gene variant in the COL11A1 gene that influences human height, although this variant alone explains only 0.1% of height variation in the Finnish population. We also demonstrated in this study that special stratification strategies such as performing sex-limited analyses, focusing on dizygous twin pairs, analyzing ethnic groups within a population separately and utilizing homogenous populations such as the Finns can improve the statistical power of finding QTL significantly. Also, we concluded from the results of this study that even though genetic effects explain a great proportion of height variance, it is likely that there are tens or even hundreds of genes with small individual effects underlying the genetic architecture of height.

Human herpesvirus 6 in multiple sclerosis and encephalitis

Human herpesvirus 6 (HHV-6) was identified from patients with HIV and lymphoproliferative diseases in 1986. It is a β-herpesvirus and is divided into two subgroups, variants A and B. HHV-6 variant B is the cause of exanthema subitum, while variant A has not yet definitely proven to cause any disease. HHV-6, especially variant A, is a highly neurotropic virus and has been associated with many diseases of the central nervous system (CNS) such as encephalitis and multiple sclerosis (MS). The present studies were aimed to elucidate the role of HHV-6 and its two variants in neurological infections. Special attention was given to study the possible role of HHV-6 in the pathogenesis of MS. We studied the expression of HHV-6 antigens using immunohistochemistry in brain autopsy samples from patients with MS and controls. HHV-6 antigen was identified in 70% of MS specimens whereas 30% of control specimens expressed HHV-6 antigen. Serum and cerebrospinal fluid (CSF) samples were collected from patients with MS and patients with other neurological diseases (OND) from patients visiting Helsinki University Central Hospital Neurological Outpatient Clinic during the years 2003 and 2004. In addition, we studied 53 children with suspected encephalitis. We developed an immunofluorescence IgG-avidity assay for the detection of primary HHV-6A and HHV-6B infection. For HHV-6B antibodies, no differences were observed between patients with MS and OND. For HHV-6A both seroprevalence and mean titers were significantly higher in MS compared to OND. HHV-6A low-avidity IgG antibodies, suggestive of primary infection, were found in serum of two, three and one patient with definite MS, possible MS and OND, respectively. From pediatric patients with suspected encephalitis, six serum samples (11.3%) contained low-avidity antibodies, indicating a temporal association between HHV-6A infection and onset of encephalitis. Three out of 26 patients with CDMS and four out of 19 patients with CPMS had HHV-6 antibodies in their CSF compared to none of the patients with OND (p=0.06 and p=0.01, respectively). Two patients with CDMS and three patients with CPMS appeared to have specific intrathecal synthesis of HHV-6A antibodies. In addition, oligoclonal bands (OCB) were observed in the CSF of five out of nine MS patients tested, and in two the OCBs reacted specifically with HHV-6 antigen, which is a novel finding. These results indicate HHV-6 specific antibody production in the CNS and suggest that there is a subset of MS patients with an active or chronic HHV-6A infection in the CNS that might be involved in the pathogenesis of MS. Our studies suggest that HHV-6 is an important causative or associated virus in some neurological infections, such as encephalitis and it might contribute to the development of MS, at least in some cases. In conclusion, HHV-6 is a neurotropic virus that should be taken into consideration when studying acute and chronic CNS diseases of unknown origin.

In vitro development of islets from human adult pancreatic tissues

Transplantation of isolated islets from cadaver pancreas is a promising possibility for the optimal treatment of type 1 diabetes. The lack of islets is a major problem. Here we have investigated the possibility of generating islets in tissue culture of human pancreatic cells. We first reproduced a previously reported method of in vitro generation of endocrine cells from human adult pancreatic tissue. By tracing the bromodeoxyuridine-labeled cells in differentiated islet buds, we found that the pancreatic progenitor cells represented a subpopulation of cytokeratin 19 (CK19)-positive ductal cells. Serum-free medium and Matrigel overlay were essential for the endocrine differentiation. We then examined the involvement of preexisting islet cells in islet neogenesis. About 6-10% of endocrine cells dedifferentiated and acquired a transitional phenotype by coexpressing CK19. Significant cell proliferation was only observed in CK19-positive cells, but not in chromogranin A-positive endocrine cells. The in vitro-derived human islets were morphologically and functionally immature when compared with normal islets. Their insulin mRNA levels were only 4-5% of that found in fresh human islets, and glucose-stimulated insulin release was 3 times lower than that of control islets. Moreover, some immature endocrine cells coexpressed insulin and glucagon. After transplantation in nude mice, the in vitro-generated islets became mature with one type of hormone per endocrine cell. In addition, we also found that also in both fresh islet transplants many cells coexpressed endocrine markers and ductal marker CK19 as a sign of ductal to endocrine cell transition. Finally, we studied the effects of clinically used immunosuppressive drugs on precursor cell proliferation and differentiation. Mycophenolate mofetil (MMF) severely hampered duct-cell proliferation, and significantly reduced the total DNA content indicating its antiproliferative effect on the precursors. Tacrolimus mainly affected differentiated beta cells by decreasing the insulin content per DNA as well as the proportion of insulin-positive cells. Sirolimus and daclizumab did not show any individual or synergistic side effects suggesting that these drugs are amenable for use in clinical islet transplantation. In summary, we confirm the capacity of endocrine differentiation from progenitors present in the adult human pancreas. The plasticity of differentiated cell types of human pancreas may be a potential mechanism of human pancreas regeneration. Ductal cell differentiation into endocrine cells in transplanted islets may be an important factor in sustaining the long-term function of islet transplants. The immunosuppressive protocol is likely to be an important determinant of long-term clinical islet graft function. Moreover, these results provide new information on the mechanisms of pancreatic islet regeneration and provide the basis for the development of new strategies for the treatment of insulin deficient diabetes mellitus.

Common SNPs explain a large proportion of the heritability for human height

SNPs discovered by genome-wide association studies (GWASs) account for only a small fraction of the genetic variation of complex traits in human populations. Where is the remaining heritability? We estimated the proportion of variance for human height explained by 294,831 SNPs genotyped on 3,925 unrelated individuals using a linear model analysis, and validated the estimation method with simulations based on the observed genotype data. We show that 45% of variance can be explained by considering all SNPs simultaneously. Thus, most of the heritability is not missing but has not previously been detected because the individual effects are too small to pass stringent significance tests. We provide evidence that the remaining heritability is due to incomplete linkage disequilibrium between causal variants and genotyped SNPs, exacerbated by causal variants having lower minor allele frequency than the SNPs explored to date.

Hundreds of variants clustered in genomic loci and biological pathways affect human height

Most common human traits and diseases have a polygenic pattern of inheritance: DNA sequence variants at many genetic loci influence the phenotype. Genome-wide association (GWA) studies have identified more than 600 variants associated with human traits, but these typically explain small fractions of phenotypic variation, raising questions about the use of further studies. Here, using 183,727 individuals, we show that hundreds of genetic variants, in at least 180 loci, influence adult height, a highly heritable and classic polygenic trait. The large number of loci reveals patterns with important implications for genetic studies of common human diseases and traits. First, the 180 loci are not random, but instead are enriched for genes that are connected in biological pathways (P = 0.016) and that underlie skeletal growth defects (P < 0.001). Second, the likely causal gene is often located near the most strongly associated variant: in 13 of 21 loci containing a known skeletal growth gene, that gene was closest to the associated variant. Third, at least 19 loci have multiple independently associated variants, suggesting that allelic heterogeneity is a frequent feature of polygenic traits, that comprehensive explorations of already-discovered loci should discover additional variants and that an appreciable fraction of associated loci may have been identified. Fourth, associated variants are enriched for likely functional effects on genes, being over-represented among variants that alter amino-acid structure of proteins and expression levels of nearby genes. Our data explain approximately 10% of the phenotypic variation in height, and we estimate that unidentified common variants of similar effect sizes would increase this figure to approximately 16% of phenotypic variation (approximately 20% of heritable variation). Although additional approaches are needed to dissect the genetic architecture of polygenic human traits fully, our findings indicate that GWA studies can identify large numbers of loci that implicate biologically relevant genes and pathways.

Hierarchical Bayesian modelling of plant pest invasions with human-mediated dispersal

Hierarchical Bayesian models can assimilate surveillance and ecological information to estimate both invasion extent and model parameters for invading plant pests spread by people. A reliability analysis framework that can accommodate multiple dispersal modes is developed to estimate human-mediated dispersal parameters for an invasive species. Uncertainty in the observation process is modelled by accounting for local natural spread and population growth within spatial units. Broad scale incursion dynamics are based on a mechanistic gravity model with a Weibull distribution modification to incorporate a local pest build-up phase. The model uses Markov chain Monte Carlo simulations to infer the probability of colonisation times for discrete spatial units and to estimate connectivity parameters between these units. The hierarchical Bayesian model with observational and ecological components is applied to a surveillance dataset for a spiralling whitefly (Aleurodicus dispersus) invasion in Queensland, Australia. The model structure provides a useful application that draws on surveillance data and ecological knowledge that can be used to manage the risk of pest movement.

Molecular pathogenesis of human parechovirus 1 infection

The Parechoviruses (HPEV) belong to the family Picornaviridae of positive-stranded RNA viruses. Although the parechovirus genome shares the general properties of other picornaviruses, the genus has several unique features when compared to other family members. We found that HPEV1 attaches to αv integrins on the cell surface and is internalized through the clathrin-mediated endocytic pathway. During he course of the infection, the Golgi was found to disintegrate and the ER membranes to swell and loose their ribosomes. The replication of HPEV1 was found to take place on small clusters of vesicles which contained the trans-Golgi marker GalT as well as the viral non-structural 2C protein. 2C was additionally found on stretches of modified ER-membranes, seemingly not involved in RNA replication. The viral non-structural 2A and 2C proteins were studied in further detail and were found to display several interesting features. The 2A protein was found to be a RNA-binding protein that preferably binds to positive sense 3 UTR RNA. It was found to bind also duplex RNA containing 3 UTR(+)-3 UTR(-), but not other dsRNA molecules studied. Mutagenesis revealed that the N-terminal basic-rich region as well as the C-terminus, are important for RNA-binding. The 2C protein on the other hand, was found to have both ATP-diphosphohydrolase and AMP kinase activities. Neither dATP nor other NTP:s were suitable substrates. Furthermore, we found that as a result of theses activities the protein is autophosphorylated. The intracellular changes brought about by the individual HPEV1 non-structural proteins were studied through the expression of fusion proteins. None of the proteins expressed were able to induce membrane changes similar to those seen during HPEV1 infection. However, the 2C protein, which could be found on the surface of lipid droplets but also on diverse intracellular membranes, was partly relocated to viral replication complexes in transfected, superinfected cells. Although Golgi to ER traffic was arrested in HPEV1-infected cells, none of the individually expressed non-structural proteins had any visible effect on the anterograde membrane traffic. Our results suggest that the HPEV1 replication strategy is different from that of many other picornaviruses. Furthermore, this study shows how relatively small differences in genome sequence result in very different intracellular pathology.

The Role of Genes on Chromosome Locus 4q12 in Human Central Nervous System Tumors

Human central nervous system (CNS) tumors are a heterogeneous group of tumors occurring in brain, brainstem and spinal cord. Malignant gliomas (astrocytic and oligodendroglial tumors), which arise from the neuroepithelial cells are the most common CNS neoplasms in human. Malignant gliomas are highly aggressive and invasive tumors, and have a very poor prognosis. The development and progression of gliomas involve a stepwise accumulation of genetic alterations that generally affect either signal transduction pathways activated by receptor tyrosine kinases (RTKs), or cell cycle arrest pathways. Constitutive activation or deregulated signaling by RTKs is caused by gene amplification, overexpression or mutations. The aberrant RTK signaling results in turn in the activation of several downstream pathways, which ultimately lead to malignant transformation and tumor proliferation. Many genetic abnormalities implicated in nervous system tumors involve the genes located at the chromosomal region 4q12. This locus harbors the receptor tyrosine kinases KIT, PDGFRA and VEGFR2, and other genes (REST, LNX1) with neural function. Gene amplification and protein expression of KIT, PDGFRA, and VEGFR2 was studied using clinical tumor material. REST and LNX1, as well as NUMBL, the interaction partner of LNX1, were studied for their gene mutations and amplifications. In our studies, amplification of LNX1 was associated with KIT and PDGFRA amplification in glioblastomas, and coamplification of KIT, PDGFRA and VEGFR2 was detected in medulloblastomas and CNS primitive neuroectodermal tumors. PDGFRA amplification was also correlated with poor overall survival. Coamplification of KIT, PDGFRA and VEGFR2 was observed in a subset of human astrocytic and oligodendroglial tumors. We suggest that genes at 4q12 could be a part of a larger amplified region, which is deregulated in gliomas, and could be used as a prognostic marker of tumorigenic process. The signaling pathways activated due to gene amplifications, activating gene mutations, and overexpressed proteins may be useful as therapeutic targets for glioma treatment. This study also includes the characterization of KIT overexpressing astrocytes, analyzed by various in vitro functional assays. Our results show that overexpression of KIT in mouse astrocytes promotes cell proliferation and anchorage-independent growth, as well as phenotypic changes in the cells. Furthermore, the increased proliferation is partly inhibited by imatinib, a small molecule inhibitor of KIT. These results suggest that KIT may play a role in astrocyte growth regulation, and might have an oncogenic role in brain tumorigenesis. Elucidation of the altered signaling pathways due to specific gene amplifications, activating gene mutations, and overexpressed proteins may be useful as therapeutic targets for glioma treatment.