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.

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...

Impact of the Human Bacterial Environment on Mycobacteriosis and Allergy

My work describes two sectors of the human bacterial environment: 1. The sources of exposure to infectious non-tuberculous mycobacteria. 2. Bacteria in dust, reflecting the airborne bacterial exposure in environments protecting from or predisposing to allergic disorders. Non-tuberculous mycobacteria (NTM) transmit to humans and animals from the environment. Infection by NTM in Finland has increased during the past decade beyond that by Mycobacterium tuberculosis. Among the farm animals, porcine mycobacteriosis is the predominant NTM disease in Finland. Symptoms of mycobacteriosis are found in 0.34 % of slaughtered pigs. Soil and drinking water are suspected as sources for humans and bedding materials for pigs. To achieve quantitative data on the sources of human and porcine NTM exposure, methods for quantitation of environmental NTM are needed. We developed a quantitative real-time PCR method, utilizing primers targeted at the 16S rRNA gene of the genus of Mycobacterium. With this method, I found in Finnish sphagnum peat, sandy soils and mud high contents of mycobacterial DNA, 106 to 107 genome equivalents per gram. A similar result was obtained by a method based on the Mycobacterium-specific hybridization of 16S rRNA. Since rRNA is found mainly in live cells, this result shows that the DNA detected by qPCR mainly represented live mycobacteria. Next, I investigated the occurrence of environmental mycobacteria in the bedding materials obtained from 5 pig farms with high prevalence (>4 %) of mycobacteriosis. When I used for quantification the same qPCR methods as for the soils, I found that piggery samples contained non-mycobacterial DNA that was amplified in spite of several mismatches with the primers. I therefore improved the qPCR assay by designing Mycobacterium-specific detection probes. Using the probe qPCR assay, I found 105 to 107 genome equivalents of mycobacterial DNA in unused bedding materials and up to 1000 fold more in the bedding collected after use in the piggery. This result shows that there was a source of mycobacteria in the bedding materials purchased by the piggery and that mycobacteria increased in the bedding materials during use in the piggery. Allergic diseases have reached epidemic proportions in urbanized countries. At the same time, childhood in rural environment or simple living conditions appears to protect against allergic disorders. Exposure to immunoreactive microbial components in rural environments seems to prevent allergies. I searched for differences in the bacterial communities of two indoor dusts, an urban house dust shown to possess immunoreactivity of the TH2-type and a farm barn dust with TH1-activity. The immunoreactivities of the dusts were revealed by my collaborators, in vitro in human dendritic cells and in vivo in mouse. The dusts accumulated >10 years in the respiratory zone (>1.5 m above floor), thus reflecting the long-term content of airborne bacteria at the two sites. I investigated these dusts by cloning and sequencing of bacterial 16S rRNA genes from dust contained DNA. From the TH2-active urban house dust, I isolated 139 16S rRNA gene clones. The most prevalent genera among the clones were Corynebacterium (5 species, 34 clones), Streptococcus (8 species, 33 clones), Staphylococcus (5 species, 9 clones) and Finegoldia (1 species, 9 clones). Almost all of these species are known as colonizers of the human skin and oral cavity. Species of Corynebacterium and Streptococcus have been reported to contain anti-inflammatory lipoarabinomannans and immunmoreactive beta-glucans respectively. Streptococcus mitis, found in the urban house dust is known as an inducer of TH2 polarized immunity, characteristic of allergic disorders. I isolated 152 DNA clones from the TH1-active farm barn dust and found species quite different from those found from the urban house dust. Among others, I found DNA clones representing Bacillus licheniformis, Acinetobacter lwoffii and Lactobacillus each of which was recently reported to possess anti-allergy immunoreactivity. Moreover, the farm barn dust contained dramatically higher bacterial diversity than the urban house dust. Exposure to this dust thus stimulated the human dendritic cells by multiple microbial components. Such stimulation was reported to promote TH1 immunity. The biodiversity in dust may thus be connected to its immunoreactivity. Furthermore, the bacterial biomass in the farm barn dust consisted of live intact bacteria mainly. In the urban house dust only ~1 % of the biomass appeared as intact bacteria, as judged by microscoping. Fragmented microbes may possess bioactivity different from that of intact cells. This was recently shown for moulds. If this is also valid for bacteria, the different immunoreactivities of the two dusts may be explained by the intactness of dustborne bacteria. Based on these results, we offer three factors potentially contributing to the polarized immunoreactivities of the two dusts: (i) the species-composition, (ii) the biodiversity and (iii) the intactness of the dustborne bacterial biomass. The risk of childhood atopic diseases is 4-fold lower in the Russian compared with the Finnish Karelia. This difference across the country border is not explainable by different geo-climatic factors or genetic susceptibilities of the two populations. Instead, the explanation must be lifestyle-related. It has already been reported that the microbiological quality of drinking water differs on the two sides of the borders. In collaboration with allergists, I investigated dusts collected from homes in the Russian Karelia and in the Finnish Karelia. I found that bacterial 16S rRNA genes cloned from the Russian Karelian dusts (10 homes, 234 clones) predominantly represented Gram-positive taxa (the phyla Actinobacteria and Firmicutes, 67%). The Russian Karelian dusts contained nine-fold more of muramic acid (60 to 70 ng mg-1) than the Finnish Karelian dusts (3 to 11 ng mg-1). Among the DNA clones isolated from the Finnish side (n=231), Gram-negative taxa (40%) outnumbered the Gram-positives (34%). Out of the 465 DNA clones isolated from the Karelian dusts, 242 were assigned to cultured validly described bacterial species. In Russian Karelia, animal-associated species e.g. Staphylococcus and Macrococcus were numerous (27 clones, 14 unique species). This finding may connect to the difference in the prevalence of allergy, as childhood contacts with pets and farm animals have been connected with low allergy risk. Plant-associated bacteria and plant-borne 16S rRNA genes (chloroplast) were frequent among the DNA clones isolated from the Finnish Karelia, indicating components originating from plants. In conclusion, my work revealed three major differences between the bacterial communtites in the Russian and in the Finnish Karelian homes: (i) the high prevalence of Gram-positive bacteria on the Russian side and of Gram-negative bacteria on the Finnish side and (ii) the rich presence of animal-associated bacteria on the Russian side whereas (iii) plant-associated bacteria prevailed on the Finnish side. One or several of these factors may connect to the differences in the prevalence of allergy.