A linguistically informed autosomal STR survey of human populations residing in the greater Himalayan region

Advancing FDSTools by integrating STRNaming 1.1

The introduction of massively parallel sequencing in forensic analysis has been facilitated with typing kits, analysis software and allele naming tools such as the ForenSeq DNA Signature Prep (DSP) kit, FDSTools and STRNaming respectively. Here we describe how FDSTools 2.0 with integrated and refined STRNaming nomenclature was validated for implementation under ISO 17025 accreditation for the ForenSeq DSP kit. Newly-added options result in efficient automatic allele calling for the majority of markers while specific settings are applied for 'novel' sequence variants to avoid the calling of remaining variable noise observed in samples sequenced with the ForenSeq DSP kit that seem to arise in the PCR. Genome-wide built-in reference data allows for greatly simplified configuration of allele naming for human targets.

An evaluation of inter and intra population structure of Uttar Pradesh, inferred from 24 autosomal STRs

Aim: Present study was designed to explore the STR diversity and genomic history of the inhabitants of the most populous subdivision of the country. A set of 24 hyper variable autosomal STRs was used to estimate the genetic diversity within the studied population. A panel of 15 autosomal STRs, which is most common in the previously reported data sets, was used to estimate the genetic diversity between the studied population and obtained unique relations were reported here.Method: The genetic diversity and polymorphism among 636 individuals of different ethnic groups, residing in Bareilly, Pilibhit, Shahjahanpur, Gorakhpur, Jhansi and Varanasi regions of Uttar Pradesh, India was investigated. This investigation was carried out via 24 autosomal STRs.Result: Studied 24 loci showed the highest value of combined power of discrimination (CPD =1), combined power of exclusion (CPE =0.99999999985), combined paternity index (CPI =6.10x10⁹) and lowest combined matching probability (CPM =7.90x10^-31).Conclusion: Studied population showed the genetic closeness with the population of Uttarakhand, Jats of Delhi, Jat Sikh (Punjab) and population of Rajasthan. The locus SE33 and Penta E was found to be most useful in the terms of highest discrimination power, lowest matching probability, highest power of exclusion and highest polymorphism information content for Uttar Pradesh population among the tested loci.

Genetic diversity, structure, and admixture in Mayans from Guatemala and Mexico based on 15 short tandem repeats

OBJECTIVE

To analyze the genetic origin, relationships, structure, and admixture in Mayan Native American groups from Guatemala and Mexico based on 15 autosomal short tandem repeats (STRs) loci commonly used in human identification (HID).

METHODS

We genotyped 513 unrelated Mayan samples from Guatemala based on 15 STR loci (AmpFlSTR® Identifiler kit). Moreover, we included 4408 genotypes previously reported, as following: Mayas from Guatemala and Mexico (n = 1666) and from Latin American, European, and African (n = 2742) populations. Forensic parameters, genetic distances, admixture, and population structure were assessed.

RESULTS

Forensic parameters of the 15 STRs in different Mayan groups from Guatemala were reported. Low (Fst = 0.78%; p = 0.000) and non-significant differentiation (Fst = 1.8%; p = 0.108) were observed in Mayas from Guatemala and Mexico, respectively. The relative homogeneity observed among Mayan groups supported theories of extensive pre-Columbian gene flow and trade throughout the Mayan Empire. The distribution of the three Native American ancestries among these Mayan groups did not support the presumable Guatemalan origin of Tojolabal and Lacandon people (South, Mexico). The nonsignificant differentiation between Ladinos and Mayas suggests a relative panmixia in Guatemala. Mestizos from southeastern Mexico and Guatemala constitute a core of Native American ancestry in Latin America related to the Mayan Empire in Central America.

CONCLUSIONS

The higher European admixture and homogeneity in Mexican Mayas of the Yucatan Peninsula suggest more intensive post-Columbian gene flow in this region than in Guatemalan Mayas.

STRNaming: Generating simple, informative names for sequenced STR alleles in a standardised and automated manner

The introduction of Massively Parallel Sequencing in the forensic domain has exposed the need for comprehensive nomenclature of sequenced Short Tandem Repeat (STR) alleles. In general, three strategies are at hand: 1) the full sequence mapped to the human genome reference sequence, which ensures exact data exchange; 2) shortened, human-readable formats for forensic reporting and data presentation and 3) very short codes that enable compact figures and tables but do not convey any sequence information. Here, we describe an algorithm of the second type: STRNaming, which generates human-readable names for sequenced STR alleles. STRNaming is guided by a reference sequence at each locus and then functions independently to automatically assign a unique, sequence-descriptive name that also includes the capillary electrophoresis allele number. STRNaming settings were established based on preferences that were surveyed internationally in the forensic community. These settings ensure that a small change in the sequence corresponds to a small change in the allele name, which is helpful for recognising for instance stutter products. Sequence variants outside of the repeat units are indicated as simple variant calls. Since the STR name is sequence-descriptive, the sequence can be traced back from the allele name. Because STRNaming is fully guided by an assignable reference sequence, no central coordination or configuration is required and the method will work for any STR locus, be it autosomal, Y-, X-chromosomal in current or future use. The algorithm is publicly available online and offline.

Dated phylogeny suggests early Neolithic origin of Sino-Tibetan languages

An accurate reconstruction of Sino-Tibetan language evolution would greatly advance our understanding of East Asian population history. Two recent phylogenetic studies attempted to do so but several of their conclusions are different from each other. Here we reconstruct the phylogeny of the Sino-Tibetan language family, using Bayesian computational methods applied to a larger and linguistically more diverse sample. Our results confirm previous work in finding that the ancestral Sino-Tibetans first split into Sinitic and Tibeto-Burman clades, and support the existence of key internal relationships. But we find that the initial divergence of this group occurred earlier than previously suggested, at approximately 8000 years before the present, coinciding with the onset of millet-based agriculture and significant environmental changes in the Yellow River region. Our findings illustrate that key aspects of phylogenetic history can be replicated in this complex language family, and calls for a more nuanced understanding of the first Sino-Tibetan speakers in relation to the "early farming dispersal" theory of language evolution.

Demographic History and Genetic Adaptation in the Himalayan Region Inferred from Genome-Wide SNP Genotypes of 49 Populations

We genotyped 738 individuals belonging to 49 populations from Nepal, Bhutan, North India, or Tibet at over 500,000 SNPs, and analyzed the genotypes in the context of available worldwide population data in order to investigate the demographic history of the region and the genetic adaptations to the harsh environment. The Himalayan populations resembled other South and East Asians, but in addition displayed their own specific ancestral component and showed strong population structure and genetic drift. We also found evidence for multiple admixture events involving Himalayan populations and South/East Asians between 200 and 2,000 years ago. In comparisons with available ancient genomes, the Himalayans, like other East and South Asian populations, showed similar genetic affinity to Eurasian hunter-gatherers (a 24,000-year-old Upper Palaeolithic Siberian), and the related Bronze Age Yamnaya. The high-altitude Himalayan populations all shared a specific ancestral component, suggesting that genetic adaptation to life at high altitude originated only once in this region and subsequently spread. Combining four approaches to identifying specific positively selected loci, we confirmed that the strongest signals of high-altitude adaptation were located near the Endothelial PAS domain-containing protein 1 and Egl-9 Family Hypoxia Inducible Factor 1 loci, and discovered eight additional robust signals of high-altitude adaptation, five of which have strong biological functional links to such adaptation. In conclusion, the demographic history of Himalayan populations is complex, with strong local differentiation, reflecting both genetic and cultural factors; these populations also display evidence of multiple genetic adaptations to high-altitude environments.

Genetic structure in the Sherpa and neighboring Nepalese populations

Background

We set out to describe the fine-scale population structure across the Eastern region of Nepal. To date there is relatively little known about the genetic structure of the Sherpa residing in Nepal and their genetic relationship with the Nepalese. We assembled dense genotype data from a total of 1245 individuals representing Nepal and a variety of different populations resident across the greater Himalayan region including Tibet, China, India, Pakistan, Kazakhstan, Uzbekistan, Tajikistan and Kirghizstan. We performed analysis of principal components, admixture and homozygosity.

Results

We identified clear substructure across populations resident in the Himalayan arc, with genetic structure broadly mirroring geographical features of the region. Ethnic subgroups within Nepal show distinct genetic structure, on both admixture and principal component analysis. We detected differential proportions of ancestry from northern Himalayan populations across Nepalese subgroups, with the Nepalese Rai, Magar and Tamang carrying the greatest proportions of Tibetan ancestry.

Conclusions

We show that populations dwelling on the Himalayan plateau have had a clear impact on the Northern Indian gene pool. We illustrate how the Sherpa are a remarkably isolated population, with little gene flow from surrounding Nepalese populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3469-5) contains supplementary material, which is available to authorized users.

Genetic Diversity in Gorkhas: an Autosomal STR Study

Genotyping of highly polymorphic autosomal short tandem repeat (STR) markers is a potent tool for elucidating genetic diversity. In the present study, fifteen autosomal STR markers were analyzed in unrelated healthy male Gorkha individuals (n = 98) serving in the Indian Army by using AmpFlSTR Identifiler Plus PCR Amplification Kit. In total, 138 alleles were observed with corresponding allele frequencies ranging from 0.005 to 0.469. The studied loci were in Hardy-Weinberg Equilibrium (HWE). Heterozygosity ranged from 0.602 to 0.867. The most polymorphic locus was Fibrinogen Alpha (FGA) chain which was also the most discriminating locus as expected. Neighbor Joining (NJ) tree and principal component analysis (PCA) plot clustered the Gorkhas with those of Nepal and other Tibeto-Burman population while lowlander Indian population formed separate cluster substantiating the closeness of the Gorkhas with the Tibeto-Burman linguistic phyla. Furthermore, the dataset of STR markers obtained in the study presents a valuable information source of STR DNA profiles from personnel for usage in disaster victim identification in military exigencies and adds to the Indian database of military soldiers and military hospital repository.

Importance of the geographic barriers to promote gene drift and avoid pre- and post-Columbian gene flow in Mexican native groups: Evidence from forensic STR Loci

OBJECTIVE

To analyze the origin, structure, relationships, and recent admixture in Mexican Native groups based on 15 STRs commonly used in human identification.

METHODS

We analyzed 39 Mexican Native population samples using STR databases based on the AmpFlSTR® Identifiler kit (n = 3,135), including Mexican-Mestizos (admixed), European and African populations, as reference.

RESULTS

Based upon effective population size (Ne) differences, Native groups were clustered into three regions: i) Center-Southeast groups, characterized by larger Ne, migration rate (Nm), genetic diversity (He), and relative homogeneity principally in the Yucatan Peninsula; ii) Isolated southern groups from Chiapas and Oaxaca, characterized by lower Ne, Nm, and He (i.e. higher isolation and genetic differentiation); iii) North-Northwest groups, which are similar to the previous group but are characterized by generating the widest gene flow barrier in the Pre-Hispanic Mexican territory, and currently by elevated admixture in some northern Native groups. Despite the relative congruence between genetic relationships with cultural, linguistic, geographic criteria, these factors do not explain the present-day population structure of Native groups, excepting in those linguistically related to the Mayan that show higher homogeneity. The Isolation by distance model was demonstrated at long distances (>1,500 km), whereas geographic isolation stands as a determining factor to avoid both non-indigenous admixture and bottleneck processes.

CONCLUSIONS

Different dynamics of gene flow and drift were observed among Mexican Native groups, highlighting the geographic barriers (mountains, canyons and jungle regions) as the main factor differentiating Pre-Hispanic populations, and eventually helping to avoid Post-European contact admixture and population bottleneck. Am J Phys Anthropol 160:298-316, 2016. © 2016 Wiley Periodicals, Inc.

Wide distribution and altitude correlation of an archaic high-altitude-adaptive EPAS1 haplotype in the Himalayas

High-altitude adaptation in Tibetans is influenced by introgression of a 32.7-kb haplotype from the Denisovans, an extinct branch of archaic humans, lying within the endothelial PAS domain protein 1 (EPAS1), and has also been reported in Sherpa. We genotyped 19 variants in this genomic region in 1507 Eurasian individuals, including 1188 from Bhutan and Nepal residing at altitudes between 86 and 4550 m above sea level. Derived alleles for five SNPs characterizing the core Denisovan haplotype (AGGAA) were present at high frequency not only in Tibetans and Sherpa, but also among many populations from the Himalayas, showing a significant correlation with altitude (Spearman’s correlation coefficient = 0.75, p value 3.9 × 10⁻¹¹). Seven East- and South-Asian 1000 Genomes Project individuals shared the Denisovan haplotype extending beyond the 32-kb region, enabling us to refine the haplotype structure and identify a candidate regulatory variant (rs370299814) that might be interacting in an additive manner with the derived G allele of rs150877473, the variant previously associated with high-altitude adaptation in Tibetans. Denisovan-derived alleles were also observed at frequencies of 3–14 % in the 1000 Genomes Project African samples. The closest African haplotype is, however, separated from the Asian high-altitude haplotype by 22 mutations whereas only three mutations, including rs150877473, separate the Asians from the Denisovan, consistent with distant shared ancestry for African and Asian haplotypes and Denisovan adaptive introgression.