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 Myopathic Protein Myotilin in Developing Mouse and in Muscle Function

Skeletal muscle cells are highly specialised in order to accomplish their function. During development, the fusion of hundreds of immature myoblasts creates large syncytial myofibres with a highly ordered cytoplasm filled with packed myofibrils. The assembly and organisation of contractile myofibrils must be tightly controlled. Indeed, the number of proteins involved in sarcomere building is impressive, and the role of many of them has only recently begun to be elucidated. Myotilin was originally identified as a high affinity a-actinin binding protein in yeast twohybrid screen. It was then found to interact also with filamin C, actin, ZASP and FATZ-1. Human myotilin is mainly expressed in striated muscle and induces efficient actin bundling in vitro and in cells. Moreover, mutations in myotilin cause different forms of muscle disease, now collectively known as myotilinopathies. In this thesis, consisting of three publications, the work on the mouse orthologue is presented. First, the cloning and molecular characterisation of the mouse myotilin gene showed that human and mouse myotilin share high sequence homology and a similar expression pattern and gene regulation. Functional analysis of the mouse promoter revealed the myogenic factor-binding elements that are required for myotilin gene transcription. Secondly, expression of myotilin was studied during mouse embryogenesis. Surprisingly, myotilin was expressed in a wide array of tissues at some stages of development; its expression pattern became more restricted at perinatal stages and in adult life. Immunostaining of human embryos confirmed broader myotilin expression compared to the sarcomeric marker titin. Finally, in the third article, targeted deletion of myotilin gene in mice revealed that it is not essential for muscle development and function. These data altogether indicate that the mouse can be used as a model for human myotilinopathy and that loss of myotilin does not alter significantly muscle structure and function. Therefore, disease-associated mutant myotilin may act as a dominant myopathic factor.

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.

Human Torque teno virus: epidemiology, cell biology and immunology

Torque teno virus (TTV) was discovered in 1997 in the serum of a Japanese patient who had a post-transfusion hepatitis of unknown etiology. It is a small virus containing a circular single-stranded DNA genome which is unique among human viruses. Within a few years after its discovery, the TTVs were noted to form a large family of viruses with numerous genotypes. TTV is highly prevalent among the general population throughout the world, and persistent infections and co-infections with several genotypes occur frequently. However, the pathogenicity and the mechanism for the sustained occurrence of the virus in blood are at present unclear. To determine the prevalence of TTV in Finland, we set up PCR methods and examined the sera of asymptomatic subjects for the presence of TTV DNA and for genotype-6 DNA. TTV was found to be highly prevalent also in Finland; 85% of adults harbored TTV in their blood, and 4% were infected with genotype-6. In addition, TTV DNA was detected in a number of different tissues, with no tissue-type or symptom specificity. Most cell-biological events during TTV infections are at the moment unknown. Replicating TTV DNA has, however, been detected in liver and the hematopoietic compartment, and three mRNAs are known to be generated. To characterize TTV cell biology in more detail, we cloned in full length the genome of TTV genotype 6. We showed that in human kidney-derived cells TTV produces altogether six proteins with distinct subcellular localizations. TTV mRNA transcription was detected in all cell lines transfected with the full-length clone, and TTV DNA replicated in several of them, including those of erythroid, kidney, and hepatic origin. Furthermore, the viral DNA replication was shown to utilize the cellular DNA polymerases. Diagnoses of TTV infections have been based almost solely on PCR, whereas serological tests, measuring antibody responses, would give more information on many aspects of these infections. To investigate the TTV immunology in more detail, we produced all six TTV proteins for use as antigens in serological tests. We detected in human sera IgM and IgG antibodies to occur simultaneously with TTV DNA, and observed appearance of TTV DNA regardless of pre-existing antibodies, and disappearance of TTV DNA after antibody appearance. The genotype-6 nucleotide sequence remained stable for years within the infected subjects, suggesting that some mechanism other than mutations is used by this minute virus to evade our immune system and to establish chronic infections in immunocompetent subjects.

Human papillomavirus infection and vaccination: Awareness and knowledge of HPV and acceptability of HPV vaccine among mothers of teenage daughters in Weihai, Shandong, China

In preparation for the introduction of human papillomavirus (HPV) vaccine, we investigated awareness and knowledge of HPV/HPV vaccine and potential acceptability to HPV vaccine among mothers with a teenage daughter in Weihai, Shandong, China. A cross-sectional survey was conducted in 2013 with a sample of 1850 mothers who had a daughter (aged 9–17 years) attending primary, junior and senior high schools. In the final sample (N = 1578, response rate 85.30%), awareness of HPV was reported by 305 (19.32%) mothers. Awareness varied significantly by daughter’s age (P<0.01), mother’s education level (P<0.01), mother’s occupation (P<0.01), household income (P<0.01) and residence type (P<0.01). Knowledge about HPV/HPV vaccine was poor with a mean total score of 3.56 (SD = 2.40) out of a possible score of 13. Mothers with a higher education level reported higher levels of knowledge (P = 0.02). Slightly more than one-fourth (26.49%) of mothers expressed their potential acceptability of HPV vaccine for their daughters. Acceptability increased along with increased daughters’ age (P<0.01), household income (P<0.01) and knowledge level (P<0.01). House wives and unemployed mothers had the highest acceptability (P<0.01). The most common reasons for not accepting HPV vaccination were “My daughter is too young to have risk of cervical cancer (30.95%)”, “The vaccine has not been widely used, and the decision will be made after it is widely used (24.91%)”, “Worry about the safety of the vaccine (22.85%)”. Awareness and knowledge of HPV/HPV vaccines are poor and HPV vaccine acceptability is low among these Chinese mothers. These results may help inform appropriate health education programs in this population.

Exclusion of angiotensinogen gene in molecular basis of human hypertension: Sibpair linkage and association analyses in Australian anglo-caucasians

Linkage with essential hypertension has been claimed for a microsatellite marker near the angiotensinogen gene (AGT; chromosome 1q42), as has association for the AGT variants M235T, G(-6)A and A(-20)C. To more rigorously evaluate AGT as a candidate gene for hypertension we performed sibpair analysis with multiple microsatellite markers surrounding this locus and using more sophisticated analysis programs. We also performed an association study of the AGT variants in unrelated subjects with a strong family history (two affected parents). For the linkage study, single and multiplex polymerase chain reaction (PCRs) and automated genescan analysis were conducted on DNA from 175 Australian Anglo-Celtic Caucasian hypertensives for the following markers: D1S2880-(2.1 cM)-D1S213-(2.8 cM)-D1S251-(6.5 cM)-AGT-(2.0 cM) -D1S235. Statistical evaluation of genotype data by nonparametric methods resulted in the following scores: Single-point analysis - SPLINK, P > 0.18; APM method, P > 0.25; ASPEX, MLOD < 0.28; SIB-PAIR, P > 0. 24; Multipoint analysis - MAPMAKER/SIBS, MLOD < 0.24; GENEHUNTER, P > 0.35. Exclusion scores of Lod -4.1 to -5.1 were obtained for these markers using MAPMAKER/SIBS for a lambda(s) of 1.6. The association study of G(-6)A, A(-20)C and M235T variants in 111 hypertensives with strong family history and 190 normotensives with no family history showed significant linkage disequilibrium between particular haplotypes, but we could find no association with hypertension. The present study therefore excludes AGT in the etiology of hypertension, at least in the population of Australian Anglo-Celtic Caucasians studied.

Editorial: Measuring refractive index distribution in the human eye lens

Non-invasive measurements of the age dependence of refractive index distribution in human eye lenses in vitro using a novel X-ray Talbot Interferometry method. In their paper, the authors make frequent reference to our own work in which we employed magnetic resonance imaging (MRI) to make similar non-invasive measurements of the refractive index distribution in the human eye lens [2, 3]. Prior to the current work, ours was the only method for making such measurements both non-invasively and without prior assumptions about the shape of the refractive index distribution. For this reason, the latest work is to be welcomed. However at several points in the paper, Pierscionek et al. [1] make statements about our technique which are factually incorrect...