#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Human genome diversity data reveal that L564P is the predominant TPC2 variant and a prerequisite for the blond hair associated M484L gain-of-function effect


Authors: Julia Böck aff001;  Einar Krogsaeter aff001;  Marcel Passon aff001;  Yu-Kai Chao aff001;  Sapna Sharma aff002;  Harald Grallert aff002;  Annette Peters aff002;  Christian Grimm aff001
Authors place of work: Walther Straub Institute of Pharmacology and Toxicology, Faculty of Medicine, Ludwig-Maximilians-Universität, Munich, Germany aff001;  Helmholtz Zentrum–Deutsches Forschungszentrum für Gesundheit und Umwelt (GmbH), Institute of Epidemiology, Neuherberg, Germany aff002
Published in the journal: Human genome diversity data reveal that L564P is the predominant TPC2 variant and a prerequisite for the blond hair associated M484L gain-of-function effect. PLoS Genet 17(1): e1009236. doi:10.1371/journal.pgen.1009236
Category: Research Article
doi: https://doi.org/10.1371/journal.pgen.1009236

Summary

The endo-lysosomal two-pore channel (TPC2) has been established as an intracellular cation channel of significant physiological and pathophysiological relevance in recent years. For example, TPC2-/- mice show defects in cholesterol degradation, leading to hypercholesterinemia; TPC2 absence also results in mature-onset obesity, and a role in glucagon secretion and diabetes has been proposed. Infections with bacterial toxins or viruses e.g., cholera toxin or Ebola virus result in reduced infectivity rates in the absence of TPC2 or after pharmacological blockage, and TPC2-/- cancer cells lose their ability to migrate and metastasize efficiently. Finally, melanin production is affected by changes in hTPC2 activity, resulting in pigmentation defects and hair color variation. Here, we analyzed several publicly available genome variation data sets and identified multiple variations in the TPC2 protein in distinct human populations. Surprisingly, one variation, L564P, was found to be the predominant TPC2 isoform on a global scale. By applying endo-lysosomal patch-clamp electrophysiology, we found that L564P is a prerequisite for the previously described M484L gain-of-function effect that is associated with blond hair. Additionally, other gain-of-function variants with distinct geographical and ethnic distribution were discovered and functionally characterized. A meta-analysis of genome-wide association studies was performed, finding the polymorphisms to be associated with both distinct and overlapping traits. In sum, we present the first systematic analysis of variations in TPC2. We functionally characterized the most common variations and assessed their association with various disease traits. With TPC2 emerging as a novel drug target for the treatment of various diseases, this study provides valuable insights into ethnic and geographical distribution of TPC2 polymorphisms and their effects on channel activity.

Keywords:

Bone density – Fibroblasts – Genome-wide association studies – Homozygosity – Human genomics – Mammalian genomics – Single nucleotide polymorphisms

Introduction

DNA variants are traditionally divided into DNA mutations and single nucleotide polymorphisms (SNPs), based on their occurrence in <1% and >1% of the population, respectively. Variations with a proven effect on gene function are also termed functional polymorphisms. Such functional polymorphisms may impact development of disease or response to pathogens, chemicals, vaccines or drugs. This is particularly important in the context of personalized medicine, as gene variation can affect, e.g. drug metabolism or side effects. Biomedical research is increasingly focusing on a better understanding of the functional relevance of such variations and their association with diseases and therapies. For instance, variations associated with a certain disease in one population may appear harmless or are irrelevant in other populations. One famous example is malaria resistance in people carrying certain variations in hemoglobin, which on the one hand results in sickle-cell anemia, but on the other hand is a selection advantage, at least for heterozygous carriers, in geographic regions where malaria is prevalent [1,2]. As with other genetic variations, functional variations in TPC2 may also impact health and survival under certain environmental conditions and in certain geographical areas, while in others they might not.

The two-pore cation channel (TPC2) has recently emerged as an intracellular ion channel of significant physiological and pathophysiological relevance. TPC2 has been associated with various essential functions in endo-lysosomes, such as trafficking, autophagy, exocytosis, and lysosomal cation/pH homeostasis [36]. Accordingly, changes to its function bear consequences for numerous diseases, e.g. infectious diseases in which viruses or bacterial toxins are trafficked through the endo-lysosomal system [79], metabolic diseases caused by defects in endo-lysosomal trafficking or cargo degradation[10,11], cancer [1216], or pathologies involving lysosome-related organelles such as platelet dense granules [17], melanosomes [18,19], or cytolytic granules [20]. Recently, two variations in TPC2 associated with hair pigmentation defects were found to be gain-of-function (GOF) variants: rs35264875 (encoding M484L) results in an increased sensitivity to the endogenous TPC2 ligand PI(3,5)P2, while rs3829241 (encoding G734E) results in reduced channel inhibition by ATP [21]. Given their direct functional relevance to channel activity, and the link between channel function and various diseases, we predict TPC2 polymorphisms confer traits which can be acted upon by various selection pressures. For example, TPC2 GOF polymorphisms could represent risk factors, enhancing endosomal trafficking and predisposing carriers to infections [79,22]. Furthermore, since TPC2 activity has been implicated in glucose homeostasis, altered TPC2 activity may correlate with the development of diabetes mellitus [2325].

In an attempt to assess the global distribution of TPC2 variations in humans, we analyzed several human genome datasets: the Simons Genome Diversity Project (SGDP) dataset, comprising genomes of 279 individuals from 142 indigenous populations in Africa, Europe, Asia, Australia, Oceania, and America [26], the 1000 Genome Project dataset (1000GP, 2504 genomes) [27], and the gnomAD data set (~140.000 genomes) [28]. We screened these datasets for polymorphisms in TPC2 and other endo-lysosomal cation channels, namely TPC1, TRPML1, TRPML2, and TRPML3. We found that TPC2 variations are more frequent than variations in the other channels. One variation in TPC2, L564P, occurs much more frequently on a global scale than the nominal “wild-type” TPC2 isoform. Other variations were found to occur in more geographically restricted manners. Unexpectedly, we found L564P to be a prerequisite for the blond hair-associated M484L variant to exert its GOF effect, as demonstrated by endo-lysosomal patch-clamp electrophysiology. Furthermore, additional variants were found to occur frequently in a homozygous manner with distinct ethnic or geographical distribution. While some of these variations had no apparent effect on channel activity, others showed increased sensitivity to the endogenous TPC2 ligand PI(3,5)P2, akin to the previously characterized M484L polymorphism. These results reveal a functional diversity of TPC2 between distinct populations, with potential relevance in manifestation and progression of TPC2-associated diseases.

Results

L564P is the predominant TPC2 variant on a global scale

We used the data provided by SGDP, 1000GP, and gnomAD to assess the occurrence of variations in the sequence of TPC2 and its related endo-lysosomal cation channels TPC1, TRPML1, TRPML2, and TRPML3 across different populations and ethnicities. In the SGDP dataset, 22 variations were identified for TPC2, while only 7 were found for TPC1 and TRPML1 each, 10 for TRPML2, and 5 for TRPML3 (Figs 1A and 1B and S1). One of the identified variations in TPC2, L564P, was found with particularly high frequency (Figs 1A, 1C, 1D, 2E and 2G). Other homozygous variations found to occur with increased frequency in TPC2 in the SGDP data set were V219I, K376R, G387D, M484L, and G734E. The polymorphisms V219I, M484L and G734E are frequently present in European populations (Fig 2A, 2D, 2F and 2G). G387D showed a low frequency across most populations (Fig 2C and 2G), while K376R (Fig 2B and 2G) showed comparably high frequencies across all populations with the exception of Native Americans, who over all showed low TPC2 variation rates in their genomes (Fig 2G). L564P occurs at a very high frequency in a homozygous manner on a global scale except for Native Americans and some African and Asian populations. Nearly all European SGDP samples were found to be homozygous for L564P, with only one sampled Hungarian being heterozygous for L564P. Similarly, in most Asian, African, Australian and Oceanian populations, the homozygous L564P variant predominates (Fig 2G).

Fig. 1. Analysis of the sequencing data derived from the SGDP data set and the 1000GP data set.
Analysis of the sequencing data derived from the SGDP data set and the 1000GP data set.
(A) Mean allele frequency in the Simons Genome Diversity Project (SGDP) data set of all homozygous (colored) and heterozygous SNPs in the TPC2 gene leading to a missense mutation. (B) Cartoon showing the estimated positions of the polymorphisms in the human TPC2 protein. (C) Map of the worldwide distribution of TPC2L564P in the SGDP data set. Each circle represents one individual from an indigenous population, respectively. White circles represent individuals homozygous for TPC2P564L, red circles represent individuals homozygous for TPC2L564P, and pink circles represent heterozygous individuals. (D) Map of the worldwide distribution of TPC2L564P in the 1000GP data set. Each circle represents one of the 26 populations (indigenous or non-indigenous (other descent) as indicated). Color coding scheme corresponds to the one in (C).
Fig. 2. Geographical distribution of all homozygous TPC2 SNPs found in the SGDP data set.
Geographical distribution of all homozygous TPC2 SNPs found in the SGDP data set.
The 297 samples were grouped according to their geographic origin. (A-F) Mean allele frequency of each SNP (mean ± SEM) grouped by geographic region, each symbol (circle, square, triangle) represents one sample (G) Mean allele frequency of each SNP grouped by continent.

Analysis of the 1000GP and the gnomAD datasets confirmed the SGDP dataset findings, including the very high frequency of the L564P variation across all populations (Figs 3 and S2 and S3). Like the SGDP dataset, the 1000GP dataset revealed a higher variation frequency in TPC2 (63 SNPs) compared to TPC1 (33 SNPs), TRPML1 (34 SNPs), TRPML2 (34 SNPs), and TRPML3 (31 SNPs) (S2 Fig). Accordingly, the gnomAD dataset supported these observations, revealing higher variation frequencies for TPC2 (412 SNPs) compared to TPC1 (381), TRPML1 (272), TRPML2 (278), and TRPML3 (260) (S3 Fig). In general, variations were more frequent in TPC2 compared to the other endo-lysosomal cation channels (S1S3 Figs). In addition to the homozygously occurring variations in TPC2 in the SGDP dataset (V219I, K376R, G387D, M484L, L564P, and G734E), another variation, S681L, was found to occur in a homozygous manner in the 1000GP and gnomAD datasets. This variation occurs mainly in Han Chinese and Japanese samples (Figs 3 and S2 and S3). Homozygously occurring variations were less frequent in other endolysosomal ion channels. In the 1000GP dataset, only one homozygously occurring variation was found in TPC1 (G803W) and three in TRPML2 (V20I, M365V, and K370Q), but none in TRPML1 and TRPML3 (S2 Fig). Only in the gnomAD dataset, homozygously occurring variations were found for TRPML1 and TRPML3 (S3 Fig). Neither of the homozygously occurring TRPML1 variations are known to cause mucolipidosis type IV, a rare neurodegenerative lysosomal storage disease caused by certain TRPML1 mutations [29,30]. For TRPML2 and TRPML3, homozygous premature stop codon variants were detected, revealing the existence of people with complete loss of either TRPML2 or TRPML3 channels.

Fig. 3. Geographical distribution of all homozygous TPC2 SNPs found in the 1000GP data set.
Geographical distribution of all homozygous TPC2 SNPs found in the 1000GP data set.
The 2054 samples were grouped according to geographic origin. (A-G) Mean allele frequency of each SNP (mean ± SEM) in each population grouped by geographic region, each symbol (circle, square, triangle, rhombus) represents one sample. (H) Mean allele frequency of each SNP grouped by continent.

TPC2 in ancient human genomes and in evolution

Next, we analyzed published genome sequencing data from ancient samples. In Neanderthal, Denisovan, and ancient modern humans, both L564 and 564P were present. Again, 564P occurred with higher frequency, suggesting that at least in the last 40.000 to 50.000 years, the occurrence and geographical distribution of L564 versus 564P was not very different from contemporary samples. An analysis of TPC2 variations in non-human primates revealed that all assessed primate reference genomes have the L564P variation, suggesting an evolutionary conserved high prevalence of 564P (S4 Fig).

TPC2 variations have predominantly developed on the background of L564P

The structure of the human TPC2 channel was recently described, pinpointing amino acid 564 to be present on the five-residue 2TMD4/2TMD5 (IIS4-S5) loop, which is of importance in channel gating [31]. As position 564 is located on a critical part of the protein, we assessed channel function of the TPC2L564P in comparison to the “wildtype” variant, TPC2P564L by using the endo-lysosomal patch-clamp technique. Stimulation of both variants with the endogenous TPC2 agonist, PI(3,5)P2, revealed no significant difference in channel conductance (Figs 4A, 4B and 5H). However, the genomic data analysis demonstrated that individuals homozygous for 564P have an increased likelihood of carrying additional homozygous variations in TPC2 (Fig 4C–4H). Only one SNP, S681L, predominantly found in East Asian populations, appeared associated with the L564 background (Fig 4H). Taking this into consideration, we re-evaluated the human donor fibroblast samples used by Chao et al. (2017) for studying the GOF variations of TPC2. We found that all fibroblasts were homozygous for the L564P polymorphism, and that the cDNA originally isolated from HEK293 cells used by Chao et al. (2017) was likewise homozygous for L564P [21]. Therefore, the question arose whether amino acid 564 may affect the recently described GOF variation of TPC2, namely TPC2M484L.

Fig. 4. TPC2L564P in association with other TPC2 SNPs.
TPC2<sup>L564P</sup> in association with other TPC2 SNPs.
(A, B) Representative PI(3,5)P2 (1 μM) activated current densities in vacuolin-enlarged lysosomal vesicles expressing TPC2L564P or TPC2P564L and the corresponding currents after block with ATP (1 mM). (C-F) Frequency of other homozygous and heterozygous TPC2 SNPs on either L564 or 564P background. (C, D) Results from the SGDP data set. (E, F) Results from the 1000GP data set. (G-H) Frequency of homozygous SNPs on either L564 or 564P background in the SGDP data set or the 1000GP data set.

L564P is essential for the blond hair associated GOF effect of M484L

Surprisingly, in endo-lysosomal patch-clamp experiments we found that M484L only acts as a GOF variant on the 564P, but not on the L564, background (Fig 5A, 5B and 5H). Thus, PI(3,5)P2-elicited currents in the TPC2M484L/P564L variant were not significantly different from TPC2P564L (“wild type”) or TPC2L564P, while TPC2M484L/L564P showed a GOF effect as reported previously [21]. To further elucidate the surprising finding that M484L only results in a GOF in combination with proline at position 564, we applied site-directed mutagenesis and found that introduction of prolines up- or downstream of 564, i.e., at position 563 or 565, did not result in GOF in combination with M484L (S7A and S7B Fig). Furthermore, we tested 564A and 564G in combination with M484L. As expected, alanine showed no GOF effect in combination with M484L. Also glycine, which may act as a helix breaker, albeit milder than proline, showed no GOF effect (Fig 5H–5J). These data suggest that only a strong helix breaking amino acid (P) but not A, L or G is capable of mediating the observed GOF effect in combination with 484L. Based on the recently resolved hTPC2 (M484/564P) cryo-EM structure [31] it is evident that L564P is situated directly between M484L and the channel pore. We therefore assume the M484L-associated GOF to be transduced through the L564P-encoding IIS4-S5 linker, dilating the channel pore further upon PI(3,5)P2 activation. The IIS4-S5 linker has already been implicated in channel gating, as an extension of the IIS4-S5 linker appears necessary to provide space for pore dilation and channel opening. Substituting the helix-initiating 564P with a leucine would dramatically affect this linker helix extension. We therefore propose a model where M484L is amplifying the effect of PI(3,5)P2 activation, requiring signal transduction through L564P to result in pore dilation. However, to confirm this, comparative cryo-EM of polymorphic channels and/or molecular dynamics simulations would be required (S7C Fig).

Fig. 5. Effect of PI(3,5)P2 on different variations of TPC2.
Effect of PI(3,5)P<sub>2</sub> on different variations of TPC2.
(A-I) Representative PI(3,5)P2 (1 μM) activated current densities in vacuolin-enlarged lysosomal vesicles isolated from HEK293 cells overexpressing different human TPC2 variants and the corresponding currents after block with ATP (1 mM). (J) Statistical summary of data as shown in A-I and in Fig 4A and 4B. Shown are average current densities (mean ± SEM) at -100 mV. Unpaired t-test were applied. *p < 0.05 and **p < 0.01.

In addition to the GOF variant TPC2M484L/P564L, endo-lysosomal patch-clamp experimentation revealed two further variants that showed moderately increased activation levels after stimulation with PI(3,5)P2 compared to control: TPC2K376R/L564P and TPC2G387D/L564P (Fig 5C, 5D and 5J). Both combinations occur in several populations with similar frequency, but are the least common in African, and Native American populations (Figs 2 and 3). In contrast, the variant TPC2V219I/L564P, mainly found in European populations, showed no significant difference compared to TPC2L564P, nor did the V219I polymorphism appear to affect the GOF effect of its associated SNP M484L (Fig 5E, 5F and 5J). As noted previously, the SNP S681L was only found on the L564 background in contrast to all other homozygously occurring variations. TPC2P564L/S681L is predominantly found in East Asian populations. Albeit showing a trend towards increased basal activity, the activity level of the TPC2P564L/S681L variant upon stimulation with PI(3,5)P2 was not significantly different from control, i.e., TPC2P564L (Fig 5G and 5J). By contrast, the blond-hair associated GOF variant TPC2M484L/L564P is barely found among East Asians and mainly present in European, Caucasian and West Asian populations (S5 Fig). All variants showed similar subcellular localization, highly correlating with lysotracker localization. In addition, no significant expression differences in western blot analysis were found (S6 Fig).

Functional analysis of TPC2 SNPs in endogenously expressing human donor fibroblast samples

As a next step we wanted to confirm these findings in endogenously expressing cells. 136 human donor samples were genotyped [21] out of which, nine donor fibroblast samples were collected for subculture and further experimentation. All fibroblast donors were of Caucasian/European origin. In accordance with the interrogated genome databases, all donor samples were homozygous for L564P (S8A Fig). Likewise, as expected, neither of the donors carried the East Asian variant S681L. Four donors (donor 1, 2, 5, and 8) were homozygous for M484L. All four were simultaneously either homo- or heterozygous for V219I, a SNP not significantly affecting basal or PI(3,5)P2-stimulated channel activity. Of the four M484L donors, three were blond, and one brown-haired, in accordance with previously published data [21]. Another blond-haired donor reported previously [21] was homozygous for G734E (donor 6). Donors 3 and 9 were both heterozygous for K376R and G387D. Donor 4 was homozygous for K376R and donor 7 homozygous for G387D (S8A Fig). Endo-lysosomal patch-clamp experiments revealed significantly higher PI(3,5)P2-stimulated TPC2 activity for all four M484L fibroblast samples compared to donor 3, which does not encode the M484L variation, nor any of the other variations in a homozygous manner (S8B and S8C Fig). Donor 4 (homozygous for K376R) and donor 7 (homozygous for G387D) both showed significantly higher PI(3,5)P2-stimulated activity compared to control, but lower activity compared to the M484L samples. This mild GOF activity of both K376R and G387D is in accordance with overexpression data described above (Fig 5). In sum, these data reproduce key findings obtained from overexpressing HEK293 cells. The data further show complex TPC2 variation patterns in different individuals with several SNPs occurring at the same time. As expected, the combination of M484L and P564L was not found among the 136 samples, as it does not seem to occur naturally.

Analysis of GWAS data

Publicly available GWAS databases were assessed and their inclusion of TPC2 polymorphisms interrogated [32,33]. Analysis of the NHGRI-EBI GWAS Catalog revealed associations of V219I, K376R, M484L, and G734E variants with hair color, which is in line with previous studies on the matter (see also S1 Table) [21,34]. Our meta-analysis revealed additional associations beyond hair pigmentation, showing V219I and M484L to be associated with a decreased risk for type 2 diabetes mellitus (T2DM). Furthermore, the G734E polymorphism appeared to be associated with bone mineral density. We investigated these findings further by interrogating the Type 2 Diabetes Knowledge Portal, a publicly available web-based GWAS tool providing access to human genetic information linked to T2DM and various other traits [33]. The top 5 significant associations of 7 TPC2 SNPs with phenotypes were collected (S2 Table) and are depicted in Fig 6A–6G. The data presented were collected from different consortia, such as the UK biobank eBMD, GIANT, GoT2D, and DIAMANTE, with sample sizes ranging from 1997 to 898130. V219I, M484L, and G734E appeared most strongly associated with hair pigmentation. V219I and M484L appeared protective against T2DM, and G734E was associated with increased bone mineral density and decreased height. Of interest, we found K376R to be associated with decreased bone mineral density and increased height, contrasting the associations of G734E. It should be further noted that the associations of V219I could be confused with associations due to the M484L polymorphism, as the two polymorphisms appear to be linked (Fig 6H). This observation is also supported by the donor fibroblast genotypes (S8 Fig). Taken together, TPC2 polymorphisms appear to be associated with various phenotypes, most consistently with hair color, T2DM, and bone mineral density.

Fig. 6. GWAS analysis to assess TPC2 SNP disease association.
GWAS analysis to assess TPC2 SNP disease association.
(A-G) The Type 2 Diabetes Knowledge Portal was employed to assess phenotypes associated with high-frequency TPC2 SNPs [33]. Odds ratios were converted into effect sizes to render the data format uniform, and standard error of the mean estimated from P values, sample sizes, and effect sizes. Forest plots were subsequently plotted for the top five associated traits, centered around an effect size of 0. Decreasing effect sizes would suggest inverse associations between the trait and the SNP (such as for M484L and type 2 diabetes), while increasing effect sizes suggest direct associations (such as for G734E and bone mineral density). Colors indicate statistical significance, with dark red data-points indicating genome-wide significance (P < 10−8). (H) We further investigated the association between V219I and other SNPs. We found V219I and M484L SNPs to be linked, with less linkage between M484L and other SNPs. Since M484L affects channel function more dramatically than V219I, we assume shared phenotypes associated with V219I to be largely attributable to the linked M484L.

Discussion

We provide here the first in-depth analysis of variations of the endo-lysosomal cation channels TPC1, TPC2, TRPML1, TRPML2, and TRPML3, across different human populations including ancient human genomes. We analyzed genomes from three different genome datasets: The Simons genome data set, the 1000 genome project dataset, and the gnomAD dataset. All three datasets revealed higher occurrence of variations in TPC2 compared to the other endo-lysosomal, non-selective cation channels TPC1, TRPML1, TRPML2, and TRPML3. Surprisingly, the TPC2 variant L564P occurs with such high frequency that the variant described hitherto as “wild-type” seems to be a rather rare variant. In addition, most other variations were found to occur on the 564P background rather than the L564 background, apart from S681L. Finally, all publicly available primate reference genomes carry proline instead of leucine at position 564. Functionally, by using endo-lysosomal patch-clamp experimentation we found that 564P has a surprising effect on the recently described M484L GOF variant of TPC2: without 564P, the M484L variation shows no GOF effect. Therefore, 564P appears to be a prerequisite for the M484L GOF effect, which is associated with a higher likelihood of light pigmentation and blond hair color, as previously reported [21,34]. The data further suggest that M484L has developed on the background of 564P, as it does not occur in combination with L564 in any of the analyzed genome samples. Furthermore, additional GOF variants of TPC2 were identified.

Given the many physiological and pathophysiological roles that TPC2 plays as outlined above, our data suggest that individuals with different TPC2 variations will likely show different susceptibilities to various diseases. Indeed, certain TPC2 polymorphisms appear to be correlated with bone mineral density (G734E, loss-of-inhibition) [35] and inversely correlate with type 2 diabetes mellitus (M484L, gain-of-function) [36], conditions previously associated with TPC2 function [2325,37]. Beyond its association with diseases, the functional spectra of TPC2 also appear of interest pharmacologically. TPC2 was recently shown to be druggable, with two groups independently developing channel agonists either by drug repurposing [38] or drug screening for high-efficacy agonists [39]. Genetically informed drug design has been estimated to double the success rate of clinical drug development [40], for example the development of SLC30A8 antagonists as a potential therapy for type 2 diabetes [41,42]. Given the relevance of TPC2 in health and disease alongside its potential translational value, the presented findings suggest the channel to be a likely contributor to multigenic traits, with implications in basic research and pharmacogenomics. Additionally, the functional variability of TPC2 appears intriguing, and warrants further investigation into how individual TPC2 polymorphisms historically have been selected for by various selection pressures.

Materials and methods

Analysis of human genome variation data sets

The Simons Genome Project (SGDP) was accessed through the Cancer Genomics Cloud, and visualized using the Integrative Genomics Viewer through loading.BAM and.BAI URLs. Available metadata was recorded, identifying sample gender, ancestry, and continental affiliation. 1000 genome dataset and primate genomes were accessed through the ncbi.nlm.nih.gov website, the gnomAD dataset through gnomad.broadinstitute.org. All missense, nonsense or stop gained SNPs of each channel were included in the analysis. Due to their indirect affiliation with the channels, intronic variants were excluded from the analysis. Wild-type samples were assigned with the number 0, heterozygous samples with 0.5 und samples with homozygous polymorphisms with the number 1, and mean allele frequencies (MAF) calculated. Published genomes of ancient samples were accessed as described in the respective publications [4362]. Only genomes that had sequencing coverage of 5 or more for the amino acid 564 of TPCN2 were included.

Human donor gDNA purification and sequencing

Genomic DNA was isolated from previously described, cultured primary human fibroblasts [21]. gDNA extraction was performed using the PureLink Genomic DNA Kit (Invitrogen) according to manufacturer’s instructions. Loci of interest were amplified by PCR using Q5 polymerase (NEB) according to manufacturer’s instructions, and primers suitable for polymorphic site amplification. Primers for amplifying and sequencing M484L and G734E loci were previously described [21]. The following primers for amplifying and sequencing the newly described TPC2 SNPs from gDNA were used: GCCCCCAGGGTTTCATTTGT (V219I forward), GAGTCCGCACATGGGTTAGG (V219I reverse, fragment size 427 bp), GGTGAAGTCAGTTTGCGCC (K376R forward), AAGGTGAGCGCCCTCGAA (K376R reverse, fragment size 444 bp), GCCATGATGGAGGTACCCG (G387D forward), GTTCCCCAGGTGATCAAGGG (G387D reverse, fragment size 381 bp), GCCTGACAGGCTGTGTGG (L564P forward), ACTCCAGTGCATAACCCGCC (L564P reverse, fragment size 303 bp), GTGCTCTTGCTTTGCTCATC (S681L forward), CGACTCTCCCATCAAAGTTCC (S681L reverse, fragment size 246 bp). Amplified PCR products were purified by agarose gel electrophoresis, band excision, and gel extraction using the QIAquick Gel Extraction Kit (Qiagen) according to manufacturer’s instructions. Subsequently, samples were sequenced using forward and reverse amplification primers to identify TPC2 genotypes.

Endo-lysosomal patch-clamp experiments

For whole-endolysosome manual patch-clamp recordings, cells were treated with vacuolin (1 μM) overnight [63]. Currents were recorded using an EPC-10 patch-clamp amplifier (HEKA, Lambrecht, Germany) and PatchMaster acquisition software (HEKA). Data were digitized at 40 kHz and filtered at 2.8 kHz. Fast and slow capacitive transients were cancelled by the compensation circuit of the EPC-10 amplifier. Recording glass pipettes were polished and had a resistance of 4–8 MΩ. Cytoplasmic solution contained 140 mM K-MSA, 5 mM KOH, 4 mM NaCl, 0.39 mM CaCl2, 1 mM EGTA and 10 mM HEPES (pH was adjusted with KOH to 7.2). Luminal solution contained 140 mM Na-MSA, 5 mM K-MSA, 2 mM Ca-MSA, 1 mM CaCl2, 10 mM HEPES and 10 mM MES (pH was adjusted with NaOH to 4.6). In all experiments, 500-ms voltage ramps from -100 to +100 mV were applied every 5 s, holding potential at -60 mV. The current amplitudes at -100 mV were extracted from individual ramp current recordings. All statistical analyses were performed using Origin8 software.

Cell culture and mutagenesis

HEK293 cells were maintained in DMEM supplemented with 10% FBS, 100 U penicillin/mL, and 100 μg streptomycin/mL at 37°C and 5% CO2. For patch clamp experiments, cells were plated on poly-L-lysine (0,1%)-coated glass coverslips 60–96 hours before experimentation, and transiently transfected with Turbofect (ThermoFisher) according to the manufacturer’s protocols and used for patch clamp experiments 36 hours after transfection.

Human genomic DNA sampling and human fibroblast isolation were approved by the Ludwig-Maximilians-Universität Ethics Committee (headed by Prof. Dr. Eisenmenger; reference no. 254–16). Acquisition of human material was performed after obtaining written informed consent by the donors as described previously [21]. Primary fibroblasts were isolated by Prof. Dr. Carola Berking and colleagues (Department of Dermatology, Ludwig-Maximilians-Universität München) from the skin of healthy adult donors. Epidermis was separated from dermis using dispase II (10 mg/mL in PBS, pH 7.2–7.4, D4693; Sigma) and dermis was digested in collagenase (1 mg/mL in DMEM, C0130; Sigma) for 22 h at room temperature. Fibroblasts were cultured in DMEM with glutamine (Life Technologies, Inc.) and 10% FBS (FBS Superior, S0615; Biochrom). Human TPC2 SNPs were generated by site-directed mutagenesis as previously described [21]. The following oligonucleotide primers were used for site-directed mutagenesis PCRs to generate newly identified TPC2 SNP variants or other mutants: V219I forward: CGGAAATGGCCAGCATCGGGCTGCTGCTGGCC; V219I reverse: GGCCAGCAGC AGCCCGATGCTGGCCATTTCCG; K376R forward: GCTGGACAGCTCCCACAGACAGGCCATGATGGAG; K376R reverse: CTCCATCATG GCCTGTCTGTGGGAGCTGTCCAGC; G387D forward: AAGGTGCGTTCCTACGACAGTGTTCTGCTGTCAGC; G387D reverse: GCTGACAGCA GAACACTGTCGTAGGAACGCACCTT; P564L forward: GCGTATCATCCCCAGCATGAAGCTGATGGCCGTGGTGGCC; P564L reverse: GGCCACCACGGCCATCAGCTTCATGCTGGGGATGATACGC; S681L forward: GTGGTGGCTGGTGTTGTCTGTCATCTGGG; S681L reverse: CCCAGATGACAGACAACACCAGCCACCAC. L564A forward: ATCATCCCCAGCATGAAGGCGATGGCCGTGGTG; L564A reverse: CACCACGGCCATCGCCTTCATGCTGGGGATGAT; L564G forward: ATCATCCCCAGCATGAAGGGGATGGCCGTGGTG; L564G reverse: CACCACGGCCATCCCCTTCATGCTGGGGATGAT; M562P forward: GCGTATCATCCCCAGCCCGAAGCTGATGGCCGTG; M562P reverse: CACGGCCATCAGCTTCGGGCTGGGGATGATACGC; K563P forward: TATCATCCCCAGCATGCCGCTGATGGCCGTGGTG; K563P reverse: CACCACGGCCATCAGCGGCATGCTGGGGATGATA; M565P forward: CCAGCATGAAGCTGCCGGCCGTGGTGGCCA; M565P reverse: TGGCCACCACGGCCGGCAGCTTCATGCTGG.

TPC2 genotype-phenotype linkage analysis

The high-frequency TPC2 polymorphisms were investigated for associations with phenotypes using the following identifiers: rs72928978 (V219I), rs3750965 (K376R), rs61746574 (G387D), rs35264875 (M484L), rs2376558 (L564P), rs78034812 (S681L), and rs3829241 (G734E). GWAS results obtained from the NHGRI-EBI GWAS Catalog were compiled in S1 Table for polymorphisms showing phenotype associations (V219I, K376R, M484L, and G734E) [32]. TPC2 genotypes were next investigated in the Type 2 Diabetes Knowledge Portal, which includes numerous phenotypes beyond type 2 diabetes mellitus (T2DM) [33]. Where necessary, odds ratios were converted into effect sizes as described previously [64]. For plotting forest plots, SEM values were extracted using R version 1.2 as follows: Traits, effect sizes, sample sizes, were extracted, reformatted, and stored in R dataframes for statistical analysis. T-statistics were calculated using the qt() function, providing reported P values and degrees of freedom as inputs. SEM was subsequently calculated upon dividing effect sizes with the respective T-statistic. Forest plots of the top five associations for each polymorphism were plotted using GraphPad Prism v8, and data-points colored according to P values reported by the Type 2 Diabetes Knowledge Portal.

Western blotting

For Western blot experiments transiently transfected HEK293 cells were washed twice with 1x PBS and pellets were collected. Total cell lysates were obtained by solubilizing in TRIS HCl 10 mM pH 8.0 and 0.2% SDS supplemented with protease and phosphatase inhibitors (Sigma). Protein concentrations were quantified via Bradford assay. Proteins were separated via a 7% sodium dodecyl sulphate polyacrylamide gel electrophoresis (SDS-PAGE; BioRad) and transferred to polyvinylidene difluoride (PVDF; BioRad) membranes. Membranes were blocked with 5% bovine serum albumin (Sigma) diluted in Tris Buffered Saline supplemented with 0.5% Tween-20 (TBS-T) for 1 h at room temperature (RT), then incubated with primary antibody at 4°C overnight. Then, membranes were washed with TBS-T and incubated with horseradish peroxidase (HRP) conjugated anti-mouse or anti-rabbit secondary antibody at RT for 1 h. Membranes were then washed and developed by incubation with Immobilon Crescendo Western HRP substrate (Merck) and by using an Odyssey imaging system (LI-COR Biosciences). Quantification was carried out using unsaturated images on ImageJ 1.52a software. The following primary antibodies were used: Rabbit polyclonal Anti-GFP, abcam, Cat# ab6556 (1:5000 in TBS-T with 5% BSA), mouse monoclonal Anti-β-tubulin, Cell Signaling, Cat# 86298S (1:1000 in TBS-T with 5% BSA). The following secondary antibodies were used: Anti-rabbit IgG, HRP-linked, Cell Signaling, Cat# 7074S (1:5000 in TBS-T with 5% BSA) and anti-mouse IgG, HRP-linked, Cell Signaling, Cat# 7076S (1:5000 in TBS-T with 5% BSA).

Statistical analysis

Details of statistical analyses and n values are provided in the Materials and Methods or the figures or figure legends. Statistical analyses were carried out using Origin 8 and GraphPad Prism 8. All error bars are depicted as mean ± SEM. Statistical significance is denoted on figures as outlined in the legends. Statistics carried out on datasets upon obtaining three independent biological and technical replicates.

Supporting information

S1 Fig [eps]
Mean allele frequency of SNPs in other endo-lysosomal cation channels in the SGDP data set.

S2 Fig [f]
Mean allele frequency of SNPs in TPC2 and other endo-lysosomal cation channels in the 1000GP data.

S3 Fig [f]
Mean allele frequency of SNPs in TPC2 and other endo-lysosomal cation channels in the gnomAD data set.

S4 Fig [a]
Analysis of ancient genome and primate samples.

S5 Fig [a]
Maps of the worldwide distribution of TPC2 in the 1000GP and the gnomAD data set.

S6 Fig [a]
Analysis of expression and subcellular localization of TPC2 variants.

S7 Fig [a]
Analysis of TPC2 mutations in the IIS4-S5 linker region.

S8 Fig [a]
TPC2 SNPs in human donor fibroblast samples.

S1 Table [xlsx]
TPC2 polymorphisms reaching genome-wide significance in GWAS studies.

S2 Table [xlsx]
Data obtained from the GWAS Catalog and Type 2 Diabetes Knowledge Portal.


Zdroje

1. Malowany JI, Butany J. Pathology of sickle cell disease. Semin Diagn Pathol. 2012;29: 49–55. doi: 10.1053/j.semdp.2011.07.005 22372205

2. Ferreira A, Marguti I, Bechmann I, Jeney V, Chora Â, Palha NR, et al. Sickle hemoglobin confers tolerance to plasmodium infection. Cell. 2011;145: 398–409. doi: 10.1016/j.cell.2011.03.049 21529713

3. Grimm C, Butz E, Chen C, Wahl-schott C, Biel M. From mucolipidosis type IV to Ebola: TRPML and two-pore channels at the crossroads of endo-lysosomal trafficking and disease. Cell Calcium. 2017;67: 148–155. doi: 10.1016/j.ceca.2017.04.003 28457591

4. Cang C, Zhou Y, Navarro B, Seo YJ, Aranda K, Shi L, et al. MTOR regulates lysosomal ATP-sensitive two-pore Na+channels to adapt to metabolic state. Cell. 2013;152: 778–790. doi: 10.1016/j.cell.2013.01.023 23394946

5. Sun W, Yue J. TPC2 mediates autophagy progression and extracellular vesicle secretion in cancer cells. Exp Cell Res. 2018;370: 478–489. doi: 10.1016/j.yexcr.2018.07.013 29990474

6. Wang X, Zhang X, Dong X, Samie M, Li X, Clapham DE, et al. TPC Proteins Are Phosphoinositide-activated Sodium-selective Ion Channels in Endosomes and Lysosomes. Cell. 2012;151: 372–383. doi: 10.1016/j.cell.2012.08.036 23063126

7. Sakurai Y, Kolokoltsov AA, Chen C-CC, Tidwell MW, Bauta WE, Klugbauer N, et al. Two-pore channels control Ebola virus host cell entry and are drug targets for disease treatment. Science. 2015;347: 995–998. doi: 10.1126/science.1258758 25722412

8. Gunaratne GS, Yang Y, Li F, Walseth TF, Marchant JS. NAADP-dependent Ca2+ signaling regulates Middle East Respiratory Syndrome-Coronavirus pseudovirus translocation through the endolysosomal system. Cell Calcium. 2018;75: 30–41. doi: 10.1016/j.ceca.2018.08.003 30121440

9. Ou X, Liu Y, Lei X, Li P, Mi D, Ren L, et al. Characterization of spike glycoprotein of SARS-CoV-2 on virus entry and its immune cross-reactivity with SARS-CoV. Nat Commun. 2020;11: 1620. doi: 10.1038/s41467-020-15562-9 32221306

10. Grimm C, Holdt LM, Chen C-C, Hassan S, Muller C, Jors S, et al. High susceptibility to fatty liver disease in two-pore channel 2-deficient mice. Nat Comms. 2014;5: 4699–4712. doi: 10.1038/ncomms5699 25144390

11. Lear P V., González-Touceda D, Couto BP, Viaño P, Guymer V, Remzova E, et al. Absence of intracellular ion channels TPC1 and TPC2 leads to mature-onset obesity in male mice, due to impaired lipid availability for thermogenesis in brown adipose tissue. Endocrinology. 2015;156: 975–986. doi: 10.1210/en.2014-1766 25545384

12. Favia A, Desideri M, Gambara G, D’Alessio A, Ruas M, Esposito B, et al. VEGF-induced neoangiogenesis is mediated by NAADP and two-pore channel-2-dependent Ca2+ signaling. Proc Natl Acad Sci. 2014;111: 4706–4715. doi: 10.1073/pnas.1406029111 25331892

13. Parrington J, Lear P, Hachem A. Calcium signals regulated by NAADP and two-pore channels -their role in development, differentiation and cancer. Int J Dev Biol. 2015;59: 327–340. doi: 10.1387/ijdb.150197ml 26198142

14. Nguyen ONP, Grimm C, Schneider LS, Chao Y-K, Atzberger C, Bartel K, et al. Two-Pore Channel Function Is Crucial for the Migration of Invasive Cancer Cells. Cancer Res. 2017;77: 1427–1438. doi: 10.1158/0008-5472.CAN-16-0852 28108508

15. Grimm C, Bartel K, Vollmar AM, Biel M. Endolysosomal cation channels and cancer—a link with great potential. Pharmaceuticals. 2018;11: 1–8. doi: 10.3390/ph11010004 29303993

16. Alharbi AF, Parrington J. Endolysosomal Ca2+ Signaling in Cancer: The Role of TPC2, From Tumorigenesis to Metastasis. Front Cell Dev Biol. 2019;7: 1–7. doi: 10.3389/fcell.2019.00001 30733944

17. Ambrosio AL, Boyle JA, Di Pietro SM. TPC2 mediates new mechanisms of platelet dense granule membrane dynamics through regulation of Ca2+ release. Mol Biol Cell. 2015;26: 3263–3274. doi: 10.1091/mbc.E15-01-0058 26202466

18. Ambrosio AL, Boyle JA, Aradi AE, Christian KA, Di Pietro SM. TPC2 controls pigmentation by regulating melanosome pH and size. Proc Natl Acad Sci U S A. 2016;113: 5622–5627. doi: 10.1073/pnas.1600108113 27140606

19. Bellono NW, Escobar IE, Oancea E. A melanosomal two-pore sodium channel regulates pigmentation. Sci Rep. 2016;6: 1–11. doi: 10.1038/s41598-016-0001-8 28442746

20. Davis LC, Morgan AJ, Chen JL, Snead CM, Bloor-Young D, Shenderov E, et al. NAADP Activates two-pore channels on t cell cytolytic granules to stimulate exocytosis and killing. Curr Biol. 2012;22: 2331–2337. doi: 10.1016/j.cub.2012.10.035 23177477

21. Chao Y-K, Schludi V, Chen C-C, Butz E, Nguyen ONP, Müller M, et al. TPC2 polymorphisms associated with a hair pigmentation phenotype in humans result in gain of channel function by independent mechanisms. Proc Natl Acad Sci. 2017;114: E8595–E8602. doi: 10.1073/pnas.1705739114 28923947

22. Grimm C, Butz E, Chen CC, Wahl-Schott C, Biel M. From mucolipidosis type IV to Ebola: TRPML and two-pore channels at the crossroads of endo-lysosomal trafficking and disease. Cell Calcium. 2017;67: 148–155. doi: 10.1016/j.ceca.2017.04.003 28457591

23. Arredouani A, Ruas M, Collins SC, Parkesh R, Clough F, Pillinger T, et al. Nicotinic Acid Adenine Dinucleotide Phosphate (NAADP) and Endolysosomal Two-pore Channels Modulate Membrane Excitability and Stimulus-Secretion Coupling in Mouse Pancreatic beta-Cells. J Biol Chem. 2015;290: 21376–21392. doi: 10.1074/jbc.M115.671248 26152717

24. Cane MC, Parrington J, Rorsman P, Galione A, Rutter GA. Cell Calcium The two pore channel TPC2 is dispensable in pancreatic ␤ -cells for normal Ca 2 + dynamics and insulin secretion. Cell Calcium. 2016;59: 32–40. doi: 10.1016/j.ceca.2015.12.004 26769314

25. Hamilton A, Zhang Q, Salehi A, Willems M, Knudsen JG, Stephen S, et al. Adrenaline stimulates glucagon secretion by Tpc2-dependent Ca2+ mobilization from acidic stores in pancreatic α-cells. Diabetes. 2018;67: 1128–1139. doi: 10.2337/db17-1102 29563152

26. Mallick S, Li H, Lipson M, Mathieson I, Gymrek M, Racimo F, et al. The Simons Genome Diversity Project: 300 genomes from 142 diverse populations. Nature. 2016/09/21. 2016;538: 201–206. doi: 10.1038/nature18964 27654912

27. Auton A, Abecasis GR, Altshuler DM, Durbin RM, Bentley DR, Chakravarti A, et al. A global reference for human genetic variation. Nature. 2015;526: 68–74. doi: 10.1038/nature15393 26432245

28. Karczewski KJ, Francioli LC, Tiao G, Cummings BB, Wang Q, Collins RL, et al. Variation across 141,456 human exomes and genomes reveals the spectrum of loss-of- function intolerance across human protein-coding genes. bioRxiv. https://doi.org/10.1101/531210

29. Sun M, Goldin E, Stahl S, Falardeau JL, Kennedy JC, Acierno JS, et al. Mucolipidosis type IV is caused by mutations in a gene encoding a novel transient receptor potential channel. Hum Mol Genet. 2000;9: 2471–2478. doi: 10.1093/hmg/9.17.2471 11030752

30. Chen C-C, Keller M, Hess M, Schiffmann R, Urban N, Wolfgardt A, et al. A small molecule restores function to TRPML1 mutant isoforms responsible for mucolipidosis type IV. Nat Commun. 2014;5: 4681–4691. doi: 10.1038/ncomms5681 25119295

31. She J, Zeng W, Guo J, Chen Q. Structural mechanisms of phospholipid activation of the human TPC2 channel. Elife. 2019;8: e45222. doi: 10.7554/eLife.45222 30860481

32. Buniello A, Macarthur JAL, Cerezo M, Harris LW, Hayhurst J, Malangone C, et al. The NHGRI-EBI GWAS Catalog of published genome-wide association studies, targeted arrays and summary statistics 2019. Nucleic Acids Res. 2019;47: D1005–D1012. doi: 10.1093/nar/gky1120 30445434

33. AMP-T2D Consortium. Type 2 Diabetes Knowledge Portal. [cited 18 Mar 2020]. Available: http://www.type2diabetesgenetics.org/

34. Sulem P, Gudbjartsson DF, Stacey SN, Helgason A, Rafnar T, Jakobsdottir M, et al. Two newly identified genetic determinants of pigmentation in Europeans. Nat Genet. 2008;40: 835–837. doi: 10.1038/ng.160 18488028

35. Morris JA, Kemp JP, Youlten SE, Laurent L, Logan JG, Chai RC, et al. An atlas of genetic influences on osteoporosis in humans and mice. Nat Genet. 2019;51: 258–266. doi: 10.1038/s41588-018-0302-x 30598549

36. Mahajan A, Wessel J, Willems SM, Zhao W, Robertson NR, Chu AY, et al. Refining the accuracy of validated target identification through coding variant fine-mapping in type 2 diabetes. Nat Genet. 2018;50: 559–571. doi: 10.1038/s41588-018-0084-1 29632382

37. Notomi T, Kuno M, Hiyama A, Nozaki T, Ohura K, Ezura Y, et al. Role of lysosomal channel protein TPC2 in osteoclast differentiation and bone remodeling under normal and low-magnesium conditions. J Biol Chem. 2017;292: 20998–21010. doi: 10.1074/jbc.M117.780072 29084844

38. Zhang X, Chen W, Li P, Calvo R, Southall N, Hu X, et al. Agonist-specific voltage-dependent gating of lysosomal two-pore Na+ channels. Elife. 2019;8: e51423. doi: 10.7554/eLife.51423 31825310

39. Gerndt S, Chen C-C, Chao Y-K, Yuan Y, Rosato AS, Krogsaeter E, et al. Agonist-mediated switching of ion selectivity in TPC2 differentially promotes lysosomal function. Elife. 2020;9: e54712. doi: 10.7554/eLife.54712 32167471

40. Nelson MR, Tipney H, Painter JL, Shen J, Nicoletti P, Shen Y, et al. The support of human genetic evidence for approved drug indications. Nat Genet. 2015;47. doi: 10.1038/ng.3164 25485836

41. Flannick J, Thorleifsson G, Beer NL, Jacobs SBR, Grarup N, Burtt NP, et al. Loss-of-function mutations in SLC30A8 protect against type 2 diabetes. Nat Publ Gr. 2014;46: 357–364. doi: 10.1038/ng.2915 24584071

42. Flannick J, Mercader JM, Fuchsberger C, Udler MS, Mahajan A, Wessel J, et al. Exome sequencing of 20,791 cases of type 2 diabetes and 24,440 controls. Nature. 2019;570: 71–76. doi: 10.1038/s41586-019-1231-2 31118516

43. Prüfer K, De Filippo C, Grote S, Mafessoni F, Korlević P, Hajdinjak M, et al. A high-coverage Neandertal genome from Vindija Cave in Croatia. Science. 2017;358: 655–658. doi: 10.1126/science.aao1887 28982794

44. Pääbo S. A high-quality Neandertal genome sequence.

45. Broushaki F, Thomas MG, Link V, López S, Dorp L Van, Hofmanová Z, et al. Early Neolithic genomes from the eastern Fertile Crescent. Science. 2016;353: 499–503. doi: 10.1126/science.aaf7943 27417496

46. Gamba C, Jones ER, Teasdale MD, McLaughlin RL, Gonzalez-Fortes G, Mattiangeli V, et al. Genome flux and stasis in a five millennium transect of European prehistory. Nat Commun. 2014;5: 1–9. doi: 10.1038/ncomms6257 25334030

47. Allentoft ME, Sikora M, Sjögren KG, Rasmussen S, Rasmussen M, Stenderup J, et al. Population genomics of Bronze Age Eurasia. Nature. 2015;522: 167–172. doi: 10.1038/nature14507 26062507

48. Schlebusch CM, Malmström H, Günther T, Sjödin P, Coutinho A, Edlund H, et al. Southern African ancient genomes estimate modern human divergence to 350,000 to 260,000 years ago. Science. 2017;358: 652–655. doi: 10.1126/science.aao6266 28971970

49. Moreno-Mayar JV, Potter BA, Vinner L, Steinrücken M, Rasmussen S, Terhorst J, et al. Terminal Pleistocene Alaskan genome reveals first founding population of Native Americans. Nature. 2018;553: 203–207. doi: 10.1038/nature25173 29323294

50. Kılınç GM, Omrak A, Özer F, Günther T, Büyükkarakaya AM, Bıçakçı E, et al. The Demographic Development of the First Farmers in Anatolia. Curr Biol. 2016/08/04. 2016;26: 2659–2666. doi: 10.1016/j.cub.2016.07.057 27498567

51. Rasmussen M, Anzick SL, Waters MR, Skoglund P, Degiorgio M, Stafford TW, et al. The genome of a Late Pleistocene human from a Clovis burial site in western Montana. Nature. 2014;506: 225–229. doi: 10.1038/nature13025 24522598

52. Fu Q, Li H, Moorjani P, Jay F, Slepchenko SM, Bondarev AA, et al. Genome sequence of a 45,000-year-old modern human from western Siberia. Nature. 2014;514: 445–449. doi: 10.1038/nature13810 25341783

53. Malaspinas A-S, Lao O, Schroeder H, Rasmussen M, Raghavan M, Moltke I, et al. Two ancient human genomes reveal Polynesian ancestry among the indigenous Botocudos of Brazil. Curr Biol. 2014/10/23. 2014;24: R1035–R1037. doi: 10.1016/j.cub.2014.09.078 25455029

54. Jones ER, Gonzalez-Fortes G, Connell S, Siska V, Eriksson A, Martiniano R, et al. Upper Palaeolithic genomes reveal deep roots of modern Eurasians. Nat Commun. 2015;6: 1–8. doi: 10.1038/ncomms9912 26567969

55. Pääbo S, Mafessoni F. High coverage Chagyrskaya Neandertal.

56. Meyer M, Kircher M, Gansauge M-T, Li H, Racimo F, Mallick S, et al. A High-Coverage Genome Sequence from an Archaic Denisovan Individual. Science. 2012;338: 222–226. doi: 10.1126/science.1224344 22936568

57. Rasmussen M, Guo X, Wang Y, Lohmueller KE, Rasmussen S, Albrechtsen A, et al. An aboriginal Australian genome reveals separate human dispersals into Asia. Science. 2011;334: 94–98. doi: 10.1126/science.1211177 21940856

58. Gallego Llorente M, Jones ER, Eriksson A, Siska V, Arthur KW, Arthur JW, et al. Ancient Ethiopian genome reveals extensive Eurasian admixture in Eastern Africa. Science. 2015;350: 821–822.

59. Sikora M, Seguin-Orlando A, Sousa VC, Albrechtsen A, Korneliussen T, Ko A, et al. Ancient genomes show social and reproductive behavior of early Upper Paleolithic foragers. Science. 2017;662: 659–662. doi: 10.1126/science.aao1807 28982795

60. Lazaridis I, Patterson N, Mittnik A, Renaud G, Mallick S, Kirsanow K, et al. Ancient human genomes suggest three ancestral populations for present-day Europeans. Nature. 2014;513: 409–413. doi: 10.1038/nature13673 25230663

61. Olalde I, Allentoft ME, Sánchez-Quinto F, Santpere G, Chiang CWK, DeGiorgio M, et al. Derived immune and ancestral pigmentation alleles in a 7,000-year-old Mesolithic European. Nature. 2014;507: 225–228. doi: 10.1038/nature12960 24463515

62. Hofmanová Z, Kreutzer S, Hellenthal G, Sell C, Diekmann Y, Díez-Del-Molino D, et al. Early farmers from across Europe directly descended from Neolithic Aegeans. Proc Natl Acad Sci. 2016;113: 6886–6891. doi: 10.1073/pnas.1523951113 27274049

63. Chen C-C, Cang C, Fenske S, Butz E, Chao Y-K, Biel M, et al. Patch-clamp technique to characterize ion channels in enlarged individual endolysosomes. Nat Protoc. 2017;12: 1639–1658. doi: 10.1038/nprot.2017.036 28726848

64. Chinn S. A simple method for converting an odds ratio to effect size for use in meta-analysis. Stat Med. 2000;19: 3127–3131. doi: 10.1002/1097-0258(20001130)19:22<3127::aid-sim784>3.0.co;2-m 11113947

65. Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35: 1547–1549. doi: 10.1093/molbev/msy096 29722887


Článek vyšel v časopise

PLOS Genetics


2021 Číslo 1
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy

Zvyšte si kvalifikaci online z pohodlí domova

Důležitost adherence při depresivním onemocnění
nový kurz
Autoři: MUDr. Eliška Bartečková, Ph.D.

Koncepce osteologické péče pro gynekology a praktické lékaře
Autoři: MUDr. František Šenk

Sekvenční léčba schizofrenie
Autoři: MUDr. Jana Hořínková, Ph.D.

Hypertenze a hypercholesterolémie – synergický efekt léčby
Autoři: prof. MUDr. Hana Rosolová, DrSc.

Multidisciplinární zkušenosti u pacientů s diabetem
Autoři: Prof. MUDr. Martin Haluzík, DrSc., prof. MUDr. Vojtěch Melenovský, CSc., prof. MUDr. Vladimír Tesař, DrSc.

Všechny kurzy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#