The application of molecular genetic approaches to the study of human evolution

Polygenic Indices (a.k.a. Polygenic Scores) in Social Science: A Guide for Interpretation and Evaluation

Polygenic indices (PGI)-the new recommended label for polygenic scores (PGS) in social science-are genetic summary scales often used to represent an individual's liability for a disease, trait, or behavior based on the additive effects of measured genetic variants. Enthusiasm for linking genetic data with social outcomes and the inclusion of premade PGIs in social science datasets have facilitated increased uptake of PGIs in social science research-a trend that will likely continue. Yet, most social scientists lack the expertise to interpret and evaluate PGIs in social science research. Here, we provide a primer on PGIs for social scientists focusing on key concepts, unique statistical genetic considerations, and best practices in calculation, estimation, reporting, and interpretation. We summarize our recommended best practices as a checklist to aid social scientists in evaluating and interpreting studies with PGIs. We conclude by discussing the similarities between PGIs and standard social science scales and unique interpretative considerations.

mtDNA analysis using Mitopore

Mitochondrial DNA (mtDNA) analysis is crucial for the diagnosis of mitochondrial disorders, forensic investigations, and basic research. Existing pipelines are complex, expensive, and require specialized personnel. In many cases, including the diagnosis of detrimental single nucleotide variants (SNVs), mtDNA analysis is still carried out using Sanger sequencing. Here, we developed a simple workflow and a publicly available webserver named Mitopore that allows the detection of mtDNA SNVs, indels, and haplogroups. To simplify mtDNA analysis, we tailored our workflow to process noisy long-read sequencing data for mtDNA analysis, focusing on sequence alignment and parameter optimization. We implemented Mitopore with eliBQ (eliminate bad quality reads), an innovative quality enhancement that permits the increase of per-base quality of over 20% for low-quality data. The whole Mitopore workflow and webserver were validated using patient-derived and induced pluripotent stem cells harboring mtDNA mutations. Mitopore streamlines mtDNA analysis as an easy-to-use fast, reliable, and cost-effective analysis method for both long- and short-read sequencing data. This significantly enhances the accessibility of mtDNA analysis and reduces the cost per sample, contributing to the progress of mtDNA-related research and diagnosis.

Analysis of the Overlength Main Noncoding Region in Metacarcinus magister (Decapoda: Brachyura) and a Phylogenetic Study of the Cancroidea Species

Complete mitochondrial genomes (mitogenomes) can provide important information regarding the molecular evolution and phylogenetic relationships of marine invertebrates, especially in Brachyura. Only one Cancroidea species of mitogenomes has been sequenced before; in this research, the mitogenomic characteristics of Metacarcinus magister (Cancridae: Cancroidea) are newly studied. The length of the M. magister mitogenome was 48,820 bp, and it contained the typical 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. We performed a series of analyses on the characteristics of the mNCR of M. magister. The phylogenetics, life circumstances, and selective pressures were all analyzed to explain the formation of this length, which revealed the length of the M. magister mitogenome to be approximately three times greater than the normal length of Brachyuran mitogenomes. Phylogenetic analyses based on a dataset of 215 Decapodan mitogenomes indicated that all Eriphioidea crabs were clustered together as a group. Moreover, the rearrangement mechanism of the Cancroidea species was predicted to provide stronger evidence for the phylogenetic analysis. In general, the results obtained in this study will contribute to a better understanding of the cause of the unusual length of the M. magister mitogenome and provide new insights into the phylogeny of Brachyura.

SNV-FEAST: microbial source tracking with single nucleotide variants

Elucidating the sources of a microbiome can provide insight into the ecological dynamics responsible for the formation of these communities. Source tracking approaches to date leverage species abundance information; however, single nucleotide variants (SNVs) may be more informative because of their high specificity to certain sources. To overcome the computational burden of utilizing all SNVs for a given sample, we introduce a novel method to identify signature SNVs for source tracking. Signature SNVs used as input into a previously designed source tracking algorithm, FEAST, can more accurately estimate contributions than species and provide novel insights, demonstrated in three case studies.

Global Distribution of Founder Variants Associated with Non-Syndromic Hearing Impairment

The genetic etiology of non-syndromic hearing impairment (NSHI) is highly heterogeneous with over 124 distinct genes identified. The wide spectrum of implicated genes has challenged the implementation of molecular diagnosis with equal clinical validity in all settings. Differential frequencies of allelic variants in the most common NSHI causal gene, gap junction beta 2 (GJB2), has been described as stemming from the segregation of a founder variant and/or spontaneous germline variant hot spots. We aimed to systematically review the global distribution and provenance of founder variants associated with NSHI. The study protocol was registered on PROSPERO, the International Prospective Register of Systematic Reviews, with the registration number "CRD42020198573". Data from 52 reports, involving 27,959 study participants from 24 countries, reporting 56 founder pathogenic or likely pathogenic (P/LP) variants in 14 genes (GJB2, GJB6, GSDME, TMC1, TMIE, TMPRSS3, KCNQ4, PJVK, OTOF, EYA4, MYO15A, PDZD7, CLDN14, and CDH23), were reviewed. Varied number short tandem repeats (STRs) and single nucleotide polymorphisms (SNPs) were used for haplotype analysis to identify the shared ancestral informative markers in a linkage disequilibrium and variants' origins, age estimates, and common ancestry computations in the reviewed reports. Asia recorded the highest number of NSHI founder variants (85.7%; 48/56), with variants in all 14 genes, followed by Europe (16.1%; 9/56). GJB2 had the highest number of ethnic-specific P/LP founder variants. This review reports on the global distribution of NSHI founder variants and relates their evolution to population migration history, bottleneck events, and demographic changes in populations linked with the early evolution of deleterious founder alleles. International migration and regional and cultural intermarriage, coupled to rapid population growth, may have contributed to re-shaping the genetic architecture and structural dynamics of populations segregating these pathogenic founder variants. We have highlighted and showed the paucity of data on hearing impairment (HI) variants in Africa, establishing unexplored opportunities in genetic traits.

Evolutionary origin of pathogenic GJB2 alleles in China

The frequency of the pathogenic allele of the autosomal recessive deafness gene GJB2 varies among different populations in the world, and accumulates to a sufficiently high frequency in certain population. The purpose of this study is to investigate the origin and evolution of GJB2 pathogenic alleles in Chinese deaf patients. Children with non-syndromic hearing loss, and their parents, from 295 families were recruited. Customized capture probes targeted at 943 SNPs related to GJB2 gene were designed for sequencing of genomic DNA in blood samples. Haplotypes carrying pathogenic allele were analyzed through linkage disequilibrium block building, ancestry tracing, and extended haplotype heterozygosity calculation. Two pathogenic GJB2 alleles, c.235delC (18.41%) and c.109G>A (15.57%), were observed in 867 donors. For c.235delC allele, 3 different core haplotypes with one major haplotype (97.32%) were found, and their core SNPs were 100% conserved. For c.109G>A allele, 6 different haplotypes with one major haplotype (93.28%) were found and the major c.109G>A allele evolved from a specific ancestral haplotype. Geographical origins of donors carrying GJB2 c.109G>A and c.235delC core haplotypes centered between Qinghai and Neimenggu. GJB2 c.235delC has long-range linkage disequilibrium. No positive selection signature was found for GJB2 c.235delC or c.109G>A in the studied population. In conclusion, we discovered a single origin of GJB2 c.235delC allele and multiple independent origins of GJB2 c.109G>A allele. Alternative to positive selection or multiple independent recurrent mutation event, population bottleneck effect might account for the observed high population frequency of these pathogenic alleles.

Genomic Insights Into the Unique Demographic History and Genetic Structure of Five Hmong-Mien-Speaking Miao and Yao Populations in Southwest China

Southern China was the original center of multiple ancestral populations related to modern Hmong-Mien, Tai-Kadai, Austroasiatic, and Austronesian people. More recent genetic surveys have focused on the fine-scale genetic structure and admixture history of southern Chinese populations, but the genetic formation and diversification of Hmong-Mien speakers are far from clear due to the sparse genetic sampling. Here, we reported nearly 700,000 single-nucleotide polymorphisms (SNPs) data from 130 Guizhou Miao and Yao individuals. We used principal component analysis, ADMIXTURE, f-statistics, qpAdm, phylogenetic tree, fineSTRUCTURE, and ALDER to explore the fine-scale population genetic structure and admixture pattern of Hmong-Mien people. The sharing allele patterns showed that our studied populations had a strong genetic affinity with ancient and modern groups from southern and southeastern East Asia. We identified one unique ancestry component maximized in Yao people, which widely existed in other Hmong-Mien-speaking populations in southern China and Southeast Asia and ancient samples of Guangxi. Guizhou Hmong-Mien speakers harbored the dominant proportions of ancestry related to southern indigenous East Asians and minor proportions of northern ancestry related to Yellow River farmers, suggesting the possibility of genetic admixture between Hmong-Mien people and recent southward Sino-Tibetan-related populations. Furthermore, we found a genetic substructure among geographically different Miao and Yao people in Leishan and Songtao. The Yao and Miao people in Leishan harbored more southern East Asian ancestry, but Miao in Songtao received more northern East Asian genetic influence. We observed high mtDNA but low Y-chromosome diversity in studied Hmong-Mien groups, supporting the role of sex-specific residence in influencing human genetic variation. Our data provide valuable clues for further exploring population dynamics in southern China.

Diversity and its causes: Lewontin on racism, biological determinism and the adaptationist programme

Lewontin's 1972 paper (RC Lewontin, 1972 The apportionment of human diversity, in Evolutionary biology, vol. 6 (eds T Dobzhansky, MK Hecht, WC Steere), pp. 381–398) can be viewed as one foray in his battle against biological determinism. Our paper shows where Lewontin, The apportionment of human diversity, fits in the debate over human classification that it stimulated. We outline three assumptions inherent in the biological deterministic view of human phenotypic diversity and show how the 1972 paper, as well as Lewontin's papers in 1970 and 1974 on the problems with the heritability statistic and his 1979 criticism of naive pan-selectionism, invalidate these assumptions. These papers were crucial components of his campaign against biological determinism and the racism with which it was associated. In the current climate of widespread racism and the rise of sociogenomics, it is important to revisit Lewontin's writings and to disseminate the messages they contain.

This article is part of the theme issue ‘Celebrating 50 years since Lewontin's apportionment of human diversity’.

A unified genealogy of modern and ancient genomes

The sequencing of modern and ancient genomes from around the world has revolutionized our understanding of human history and evolution. However, the problem of how best to characterize ancestral relationships from the totality of human genomic variation remains unsolved. Here, we address this challenge with nonparametric methods that enable us to infer a unified genealogy of modern and ancient humans. This compact representation of multiple datasets explores the challenges of missing and erroneous data and uses ancient samples to constrain and date relationships. We demonstrate the power of the method to recover relationships between individuals and populations as well as to identify descendants of ancient samples. Finally, we introduce a simple nonparametric estimator of the geographical location of ancestors that recapitulates key events in human history.

Review of the Forensic Applicability of Biostatistical Methods for Inferring Ancestry from Autosomal Genetic Markers

The inference of ancestry has become a part of the services many forensic genetic laboratories provide. Interest in ancestry may be to provide investigative leads or identify the region of origin in cases of unidentified missing persons. There exist many biostatistical methods developed for the study of population structure in the area of population genetics. However, the challenges and questions are slightly different in the context of forensic genetics, where the origin of a specific sample is of interest compared to the understanding of population histories and genealogies. In this paper, the methodologies for modelling population admixture and inferring ancestral populations are reviewed with a focus on their strengths and weaknesses in relation to ancestry inference in the forensic context.

De novo mutation rates at the single-mutation resolution in a human HBB gene-region associated with adaptation and genetic disease

Although it is known that the mutation rate varies across the genome, previous estimates were based on averaging across various numbers of positions. Here, we describe a method to measure the origination rates of target mutations at target base positions and apply it to a 6-bp region in the human hemoglobin subunit beta (HBB) gene and to the identical, paralogous hemoglobin subunit delta (HBD) region in sperm cells from both African and European donors. The HBB region of interest (ROI) includes the site of the hemoglobin S (HbS) mutation, which protects against malaria, is common in Africa, and has served as a classic example of adaptation by random mutation and natural selection. We found a significant correspondence between de novo mutation rates and past observations of alleles in carriers, showing that mutation rates vary substantially in a mutation-specific manner that contributes to the site frequency spectrum. We also found that the overall point mutation rate is significantly higher in Africans than in Europeans in the HBB region studied. Finally, the rate of the 20A→T mutation, called the “HbS mutation” when it appears in HBB, is significantly higher than expected from the genome-wide average for this mutation type. Nine instances were observed in the African HBB ROI, where it is of adaptive significance, representing at least three independent originations; no instances were observed elsewhere. Further studies will be needed to examine mutation rates at the single-mutation resolution across these and other loci and organisms and to uncover the molecular mechanisms responsible.

Population affinity and variation of sexual dimorphism in three-dimensional facial forms: comparisons between Turkish and Japanese populations

Examining the extent to which sex differences in three-dimensional (3D) facial soft tissue configurations are similar across diverse populations could suggest the source of the indirect evolutionary benefits of facial sexual dimorphism traits. To explore this idea, we selected two geographically distinct populations. Three-dimensional model faces were derived from 272 Turkish and Japanese men and women; their facial morphologies were evaluated using landmark and surface-based analyses. We found four common facial features related to sexual dimorphism. Both Turkish and Japanese females had a shorter lower face height, a flatter forehead, greater sagittal cheek protrusion in the infraorbital region but less prominence of the cheek in the parotid-masseteric region, and an antero-posteriorly smaller nose when compared with their male counterparts. The results indicated the possible phylogenetic contribution of the masticatory organ function and morphogenesis on sexual dimorphism of the human face in addition to previously reported biological and psychological characteristics, including sexual maturity, reproductive potential, mating success, general health, immune response, age, and personality.

Genetic characterisation of the North-West Indian populations: analysis of mitochondrial DNA control region variations

BACKGROUND

Human mitochondrial DNA presents several interesting characteristics, making it a favourable tool in the field of molecular anthropology, medical genetics, population history, and forensic science.

AIM

The present study investigated the mitochondrial DNA (mtDNA) control region variations in diverse ethnic groups of North-West India for which population data is insufficient.

SUBJECTS AND METHODS

The complete mtDNA control regions of 197 unrelated (for up to three generations) healthy individuals belonging to different ethnic groups of North-West India were sequenced. The haplotype frequencies, haplogroup distribution, and pairwise F_ST values between the studied and other worldwide populations were generated to study patterns of variation in human mtDNA.

RESULTS

The results ascertained high gene diversity (0.998) in the studied maternal lineages, identifying 166 distinct haplotypes, of which 158 were unique and characterised by 117 variable sites. Three haplogroups: M3, M30, and U7 were observed to be the most prevalent, and phylogeographically a total of 55.86% of sequences were characterised into South Asian, followed by West Eurasian (40.18%) and East Asian (3.96%), ancestry haplogroups. Pairwise genetic differentiation comparisons revealed maternal homogeneity in the studied groups. No population substructure was detected within the North-West Indian populations.

CONCLUSION

The results of this preliminary study will contribute to an existing database of mtDNA variations of the Indian population and facilitate prospective studies investigating population genetics and human diseases.

Preliminary study on the function of TMEM50A and its correlation with the RH genes

OBJECTIVES

The purpose of this study was to investigate the association and impact of TMEM50A on RH genes activity and function.

BACKGROUND

SMP1 is located on chromosome 1p36.11 in the RH gene locus, between the RHD and RHCE gene, where its position may be linked to RH haplotypes and contribute to selective pressures regarding certain RH haplotypes. TMEM50A is encoded by the SMP1 located in the intergenic region of RH, its influence on the function of the RH genes remains unclear.

METHODS

The expression of TMEM50A was regulated by transfection of plasmid and siRNA in K562 cell model. Western blot and real-time PCR were used to detect possible expression changes in the RH. The ammonium transport function of cells was monitored using pH-sensitive dye, while transcriptome sequencing was used to predict the potential function of TMEM50A.

RESULTS

The overexpression of TMEM50A significantly up-regulated RHCE gene activity (63.56%). The inhibition of TMEM50A resulted in significantly decreased RHCE (41.82%) and RHD expression (27.35%). Compared to control group, there was no significant change in the NH₄ ⁺ transport function of cells in the overexpressed TMEM50A group. Transcriptome analysis showed that TMEM50A not only affected the transcription of target gene through splicing activities, but also played a role in the development of embryonic nervous system.

CONCLUSIONS

TMEM50A may regulate the expression of RH gene by affecting the stability of RH mRNA through splicing function. It speculates that TMEM50A may play an important role in the development of embryonic nervous system.

Archetypes and the 'Impoverished Genome' argument: updates from evolutionary genetics

Throughout his career, Jung felt the psyche had 'ancestral layers' that contained elements of an individual's species history, and clinical experience has shown that this idea can be an aid to psychological healing and emotional well-being. Thus, some later thinkers have attempted to link such theoretical constructs to the genome, as Jung had little knowledge of genetics in his day. But in the early 2000s, genome studies suggested that the genome might contain too little content to be capable of encoding symbolic information. This opinion gave rise to an oft-repeated 'impoverished genome' argument, i.e. that the genome could not provide a significant contribution to the collective unconscious, prompting theorists to propose other sources for it, or to argue that it doesn't exist. Today, however, developments in evolutionary neurogenetics calls the impoverished genome argument into question for a number of independent reasons. These developments re-open the idea that the genome may be worth reconsidering as the biological substrate for the collective unconscious.

The Y Chromosome: A Complex Locus for Genetic Analyses of Complex Human Traits

The Human Y chromosome (ChrY) has been demonstrated to be a powerful tool for phylogenetics, population genetics, genetic genealogy and forensics. However, the importance of ChrY genetic variation in relation to human complex traits is less clear. In this review, we summarise existing evidence about the inherent complexities of ChrY variation and their use in association studies of human complex traits. We present and discuss the specific particularities of ChrY genetic variation, including Y chromosomal haplogroups, that need to be considered in the design and interpretation of genetic epidemiological studies involving ChrY.

You Will Never Walk Alone: Codispersal of JC Polyomavirus with Human Populations

JC polyomavirus (JCPyV) is one of the most prevalent human viruses. Findings based on the geographic distribution of viral subtypes suggested that JCPyV codiverged with human populations. This view was however challenged by data reporting a much more recent origin and expansion of JCPyV. We collected information on ∼1,100 worldwide strains and we show that their geographic distribution roughly corresponds to major human migratory routes. Bayesian phylogeographic analysis inferred a Subsaharan origin for JCPyV, although with low posterior probability. High confidence inference at internal nodes provided strong support for a long-standing association between the virus and human populations. In line with these data, pairwise FST values for JCPyV and human mtDNA sampled from the same areas showed a positive and significant correlation. Likewise, very strong relationships were found when node ages in the JCPyV phylogeny were correlated with human population genetic distances (nuclear-marker based FST). Reconciliation analysis detected a significant cophylogenetic signal for the human population and JCPyV trees. Notably, JCPyV also traced some relatively recent migration events such as the expansion of people from the Philippines/Taiwan area into Remote Oceania, the gene flow between North-Eastern Siberian and Ainus, and the Koryak contribution to Circum-Arctic Americans. Finally, different molecular dating approaches dated the origin of JCPyV in a time frame that precedes human out-of-Africa migration. Thus, JCPyV infected early human populations and accompanied our species during worldwide dispersal. JCPyV typing can provide reliable geographic information and the virus most likely adapted to the genetic background of human populations.

Identification of single nucleotide variants in the Moroccan population by whole-genome sequencing

BACKGROUND

Large-scale human sequencing projects have described around a hundred-million single nucleotide variants (SNVs). These studies have predominately involved individuals with European ancestry despite the fact that genetic diversity is expected to be highest in Africa where Homo sapiens evolved and has maintained a large population for the longest time. The African Genome Variation Project examined several African populations but these were all located south of the Sahara. Morocco is on the northwest coast of Africa and mostly lies north of the Sahara, which makes it very attractive for studying genetic diversity. The ancestry of present-day Moroccans is unknown and may be substantially different from Africans found South of the Sahara desert, Recent genomic data of Taforalt individuals in Eastern Morocco revealed 15,000-year-old modern humans and suggested that North African individuals may be genetically distinct from previously studied African populations.

RESULTS

We present SNVs discovered by whole genome sequencing (WGS) of three Moroccans. From a total of 5.9 million SNVs detected, over 200,000 were not identified by 1000G and were not in the extensive gnomAD database. We summarise the SNVs by genomic position, type of sequence gene context and effect on proteins encoded by the sequence. Analysis of the overall genomic information of the Moroccan individuals to individuals from 1000G supports the Moroccan population being distinct from both sub-Saharan African and European populations.

CONCLUSIONS

We conclude that Moroccan samples are genetically distinct and lie in the middle of the previously observed cline between populations of European and African ancestry. WGS of Moroccan individuals can identify a large number of novel SNVs and aid in functional characterisation of the genome.

A comprehensive survey of models for dissecting local ancestry deconvolution in human genome

Over the past decade, studies of admixed populations have increasingly gained interest in both medical and population genetics. These studies have so far shed light on the patterns of genetic variation throughout modern human evolution and have improved our understanding of the demographics and adaptive processes of human populations. To date, there exist about 20 methods or tools to deconvolve local ancestry. These methods have merits and drawbacks in estimating local ancestry in multiway admixed populations. In this article, we survey existing ancestry deconvolution methods, with special emphasis on multiway admixture, and compare these methods based on simulation results reported by different studies, computational approaches used, including mathematical and statistical models, and biological challenges related to each method. This should orient users on the choice of an appropriate method or tool for given population admixture characteristics and update researchers on current advances, challenges and opportunities behind existing ancestry deconvolution methods.

Comparative study between Helicobacter pylori and host human genetics in the Dominican Republic

Background

Helicobacter pylori, a bacterium that infects the human stomach, has high genetic diversity. Because its evolution is parallel to human, H. pylori is used as a tool to trace human migration. However, there are few studies about the relationship between phylogeography of H. pylori and its host human.

Methods

We examined both H. pylori DNA and the host mitochondrial DNA and Y-chromosome DNA obtained from a total 119 patients in the Dominican Republic, where human demography consists of various ancestries. DNA extracted from cultured H. pylori were analyzed by multi locus sequence typing. Mitochondrial DNA and Y-chromosome DNA were evaluated by haplogroup analyses.

Results

H. pylori strains were divided into 2 populations; 68 strains with African group (hpAfrica1) and 51 strains with European group (hpEurope). In Y-chromosomal haplogroup, European origin was dominant, whereas African origin was dominant both in H. pylori and in mtDNA haplogroup. These results supported the hypothesis that mother-to-child infection is predominant in H. pylori infection. The Amerindian type of mtDNA haplogroup was observed in 11.8% of the patients; however, Amerindian type (hspAmerind) of H. pylori was not observed. Although subpopulation type of most hpAfrica1 strains in Central America and South America were hybrid (hspWAfrica/hpEurope), most Dominican Republic hpAfrica1 strains were similar to those of African continent.

Conclusions

Genetic features of H. pylori, mtDNA, and Y haplogroups reflect the history of colonial migration and slave trade in the Dominican Republic. Discrepancy between H. pylori and the host human genotypes support the hypothesis that adaptability of hspAmerind H. pylori strains are weaker than hpEurope strains. H. pylori strains in the Dominican Republic seem to contain larger proportion of African ancestry compared to other American continent strains.

The Association Between Schizophrenia Risk Variants and Creativity in Healthy Han Chinese Subjects

Although previous evidence has suggested that there is a genetic link between schizophrenia and creativity, the specific genetic variants that underlie the link are still largely unknown. To further explore the potential genetic link between schizophrenia and creativity, in a sample of 580 healthy Han Chinese subjects, this study aimed to (1) validate the role of Neuregulin 1 (NRG1) rs6994992 (one schizophrenia risk variant that has been previously linked to creativity in the European population) in the relationship between schizophrenia and creativity and (2) explore the associations between 10 other schizophrenia risk variants and creativity. For NRG1 rs6994992, the result validated its association with creativity measures. However, since NRG1 rs6994992 is not a schizophrenia risk variant in the Han Chinese population, the validated association suggested that ethnic difference may exist in the relationship between NRG1 rs6994992, schizophrenia and creativity. For other schizophrenia risk variants, the result only demonstrated a nominal association between ZNF536 rs2053079 and creativity measures which would not survive correction for multiple testing. No association between polygenic risk score for these 10 schizophrenia risk variants and creativity measures was observed. In conclusion, this study provides limited evidence for the associations between these schizophrenia risk variants and creativity in healthy Han Chinese subjects. Future studies are warranted to better understand the potential genetic link between schizophrenia and creativity.

MtDNA polymorphism analyses in the Chinese Mongolian group: Efficiency evaluation and further matrilineal genetic structure exploration

BACKGROUND

Profiling of mitochondrial DNA is surely to provide valuable investigative clues for forensic cases involving highly degraded specimens or complex maternal lineage kinship determination. But traditionally used hypervariable region sequencing of mitochondrial DNA is less frequently suggested by the forensic community for insufficient informativeness. Genome-wide sequencing of mitochondrial DNA can provide considerable amount of variant information but can be high cost at the same time.

METHODS

Efficiency of the 60 mitochondrial DNA polymorphic sites dispersing across the control region and coding region of mitochondrial DNA genome was evaluated with 106 Mongolians recruited from the Xinjiang Uyghur Autonomous Region, China, and allele-specific PCR technique was employed for mitochondrial DNA typing.

RESULTS

Altogether 58 haplotypes were observed and the haplotypic diversity, discrimination power and random match probability were calculated to be 0.981, 0.972, and 0.028, respectively. Mitochondrial DNA haplogroup affiliation exhibited an exceeding percentage (12.26%) of west Eurasian lineage (H haplogroup) in the studied Mongolian group, which needed to be further verified with more samples. Furthermore, the genetic relationships between the Xinjiang Mongolian group and the comparison populations were also investigated and the genetic affinity was discovered between the Xinjiang Mongolian group and the Xinjiang Kazak group in this study.

CONCLUSION

It was indicated that the panel was potentially enough to be used as a supplementary tool for forensic applications. And the matrilineal genetic structure analyses based on mitochondrial DNA variants in the Xinjiang Mongolian group could be helpful for subsequent anthropological studies.

The inference of sex-biased human demography from whole-genome data

Sex-biased demographic events ("sex-bias") involve unequal numbers of females and males. These events are typically inferred from the relative amount of X-chromosomal to autosomal genetic variation and have led to conflicting conclusions about human demographic history. Though population size changes alter the relative amount of X-chromosomal to autosomal genetic diversity even in the absence of sex-bias, this has generally not been accounted for in sex-bias estimators to date. Here, we present a novel method to identify sex-bias from genetic sequence data that models population size changes and estimates the female fraction of the effective population size during each time epoch. Compared to recent sex-bias inference methods, our approach can detect sex-bias that changes on a single population branch without requiring data from an outgroup or knowledge of divergence events. When applied to simulated data, conventional sex-bias estimators are biased by population size changes, especially recent growth or bottlenecks, while our estimator is unbiased. We next apply our method to high-coverage exome data from the 1000 Genomes Project and estimate a male bias in Yorubans (47% female) and Europeans (44%), possibly due to stronger background selection on the X chromosome than on the autosomes. Finally, we apply our method to the 1000 Genomes Project Phase 3 high-coverage Complete Genomics whole-genome data and estimate a female bias in Yorubans (63% female), Europeans (84%), Punjabis (82%), as well as Peruvians (56%), and a male bias in the Southern Han Chinese (45%). Our method additionally identifies a male-biased migration out of Africa based on data from Europeans (20% female). Our results demonstrate that modeling population size change is necessary to estimate sex-bias parameters accurately. Our approach gives insight into signatures of sex-bias in sexual species, and the demographic models it produces can serve as more accurate null models for tests of selection.

Worldwide human mitochondrial haplogroup distribution from urban sewage

Community level genetic information can be essential to direct health measures and study demographic tendencies but is subject to considerable ethical and legal challenges. These concerns become less pronounced when analyzing urban sewage samples, which are ab ovo anonymous by their pooled nature. We were able to detect traces of the human mitochondrial DNA (mtDNA) in urban sewage samples and to estimate the distribution of human mtDNA haplogroups. An expectation maximization approach was used to determine mtDNA haplogroup mixture proportions for samples collected at each different geographic location. Our results show reasonable agreement with both previous studies of ancient evolution or migration and current US census data; and are also readily reproducible and highly robust. Our approach presents a promising alternative for sample collection in studies focusing on the ethnic and genetic composition of populations or diseases associated with different mtDNA haplogroups and genotypes.

Crossbred evaluations using single-step genomic BLUP and algorithm for proven and young with different sources of data1

Genomic selection (GS) is routinely applied to many purebreds and lines of farm species. However, this method can be extended to predictions across purebreds as well as for crossbreds. This is useful for swine and poultry, for which selection in nucleus herds is typically performed on purebred animals, whereas the commercial product is the crossbred animal. Single-step genomic BLUP (ssGBLUP) is a widely applied method that can explore the recently developed algorithm for proven and young (APY). The APY allows for greater efficiency in computing BLUP solutions by exploiting the theory of limited dimensionality of genomic information and chromosome segments (Me). This study investigates the predictivity as a proxy for accuracy across and within 2 purebred pig lines and their crosses, under the application of APY in ssGBLUP setup, and different levels of Me overlapping across populations. The data consisted of approximately 210k phenotypic records for 2 traits (T1 and T2) with moderate heritability. Genotypes for 43k SNP markers were available for approximately 46k animals, from which 26k and 16k belong to 2 pure lines (L1 and L2), and approximately 4k are crossbreds. The complete pedigree had more than 720k animals. Different multivariate ssGBLUP models were applied, either with the regular or APY inverse of the genomic relationship matrix (G). The models included a standard bivariate animal model with 3 lines evaluated as 1 joint line, and for each trait individually, a 3-trait animal model with each line treated as a separate trait. Both models provided the same predictivity across and within the lines. Using either of the pure lines data as a training set resulted in a similar predictivity for the crossbreed animals (0.18 to 0.21). Across-line predictive ability was limited to less than half of the maximum predictivity for each line (L1T1 0.33, L1T2 0.25, L2T1 0.35, L2T2 0.36). For crossbred predictions, APY performed equivalently to regular G inverse when the number of core animals was equal to the number of eigenvalues explaining between 98% and 99% of the variance of G (4k to 8k) including all lines. Predictivity across the lines is achievable because of the shared Me between them. The number of those shared segments can be obtained via eigenvalue decomposition of genomic information available for each line.

Genomic analyses in African populations identify novel risk loci for cleft palate

Orofacial clefts are common developmental disorders that pose significant clinical, economical and psychological problems. We conducted genome-wide association analyses for cleft palate only (CPO) and cleft lip with or without palate (CL/P) with ~17 million markers in sub-Saharan Africans. After replication and combined analyses, we identified novel loci for CPO at or near genome-wide significance on chromosomes 2 (near CTNNA2) and 19 (near SULT2A1). In situ hybridization of Sult2a1 in mice showed expression of SULT2A1 in mesenchymal cells in palate, palatal rugae and palatal epithelium in the fused palate. The previously reported 8q24 was the most significant locus for CL/P in our study, and we replicated several previously reported loci including PAX7 and VAX1.

Genetic diversity of the Thao people of Taiwan using Y-chromosome, mitochondrial DNA and HLA gene systems

BACKGROUND

Despite attempts in retracing the history of the Thao people in Taiwan using folktales, linguistics, physical anthropology, and ethnic studies, their history remains incomplete. The heritage of Thao has been associated with the Pazeh Western plains peoples and several other mountain peoples of Taiwan. In the last 400 years, their culture and genetic profile have been reshaped by East Asian migrants. They were displaced by the Japanese and the construction of a dam and almost faced extinction. In this paper, genetic information from mitochondrial DNA (mtDNA), Histoleucocyte antigens (HLA), and the non-recombining Y chromosome of 30 Thao individuals are compared to 836 other Taiwan Mountain and Plains Aborigines (TwrIP & TwPp), 384 Non-Aboriginal Taiwanese (non-TwA) and 149 Continental East Asians.

RESULTS

The phylogeographic analyses of mtDNA haplogroups F4b and B4b1a2 indicated gene flow between Thao, Bunun, and Tsou, and suggested a common ancestry from 10,000 to 3000 years ago. A claim of close contact with the heavily Sinicized Pazeh of the plains was not rejected and suggests that the plains and mountain peoples most likely shared the same Austronesian agriculturist gene pool in the Neolithic.

CONCLUSIONS

Having been moving repeatedly since their arrival in Taiwan between 6000 and 4500 years ago, the Thao finally settled in the central mountain range. They represent the last plains people whose strong bonds with their original culture allowed them to preserve their genetic heritage, despite significant gene flow from the mainland of Asia. Representing a considerable contribution to the genealogical history of the Thao people, the findings of this study bear on ongoing anthropological and linguistic debates on their origin.

Immune Gene Diversity in Archaic and Present-day Humans

Genome-wide analyses of two Neandertals and a Denisovan have shown that these archaic humans had lower genetic heterozygosity than present-day people. A similar reduction in genetic diversity of protein-coding genes (gene diversity) was found in exome sequences of three Neandertals. Reduced gene diversity, particularly in genes involved in immunity, may have important functional consequences. In fact, it has been suggested that reduced diversity in immune genes may have contributed to Neandertal extinction. We therefore explored gene diversity in different human groups, and at different time points on the Neandertal lineage, with a particular focus on the diversity of genes involved in innate immunity and genes of the Major Histocompatibility Complex (MHC).

We find that the two Neandertals and a Denisovan have similar gene diversity, all significantly lower than any present-day human. This is true across gene categories, with no gene set showing an excess decrease in diversity compared with the genome-wide average. Innate immune-related genes show a similar reduction in diversity to other genes, both in present-day and archaic humans. There is also no observable decrease in gene diversity over time in Neandertals, suggesting that there may have been no ongoing reduction in gene diversity in later Neandertals, although this needs confirmation with a larger sample size. In both archaic and present-day humans, genes with the highest levels of diversity are enriched for MHC-related functions. In fact, in archaic humans the MHC genes show evidence of having retained more diversity than genes involved only in the innate immune system.

Identification of CACNA1D variants associated with sinoatrial node dysfunction and deafness in additional Pakistani families reveals a clinical significance

Sinoatrial node dysfunction and deafness (SANDD) syndrome is rare and characterized by a low heart beat and severe-to-profound deafness. Additional features include fatigue, dizziness, and episodic syncope. The sinoatrial node (SAN) drives heart automaticity and continuously regulates heart rate. The CACNA1D gene encoding the Ca_v1.3 protein expressed in inner hair cells, atria and SAN, induces loss-of-function in channel activity and underlies SANDD. To date, only one variant c.1208_1209insGGG:p.(G403_V404insG) has been reported for SANDD syndrome. We studied five Pakistani families with SANDD and characterized a new missense variant p.(A376V) in CACNA1D in one family, and further characterized the founder variant p.(G403_V404insG) in four additional pedigrees. We show that affected individuals in the four families which segregate p.(G403_V404insG) share a 1.03 MB haplotype on 3p21.1 suggesting they share a common distant ancestor. In conclusion, we identified new and known variants in CACNA1D in five Pakistani families with SANDD. This study is of clinical importance as the CACNA1D founder variant is only observed in families from the Khyber Pakhtunkhwa (KPK) province, in Pakistan. Therefore, screening patients with congenital deafness for SAN dysfunction in this province could ensure adequate follow-up and prevent cardiac failure associated with SAN.

Patterns of Genetic Coding Variation in a Native American Population before and after European Contact

The effects of European colonization on the genomes of Native Americans may have produced excesses of potentially deleterious features, mainly due to the severe reductions in population size and corresponding losses of genetic diversity. This assumption, however, neither considers actual genomic patterns that existed before colonization nor does it adequately capture the effects of admixture. In this study, we analyze the whole-exome sequences of modern and ancient individuals from a Northwest Coast First Nation, with a demographic history similar to other indigenous populations from the Americas. We show that in approximately ten generations from initial European contact, the modern individuals exhibit reduced levels of novel and low-frequency variants, a lower proportion of potentially deleterious alleles, and decreased heterozygosity when compared to their ancestors. This pattern can be explained by a dramatic population decline, resulting in the loss of potentially damaging low-frequency variants, and subsequent admixture. We also find evidence that the indigenous population was on a steady decline in effective population size for several thousand years before contact, which emphasizes regional demography over the common conception of a uniform expansion after entry into the Americas. This study examines the genomic consequences of colonialism on an indigenous group and describes the continuing role of gene flow among modern populations.

Blood Groups Distribution and Gene Diversity of the ABO and Rh (D) Loci in the Mexican Population

Objective

To determine the frequency and distribution of ABO and Rh (D) antigens and, additionally, investigate gene diversity and the structure of Mexican populations.

Materials and Methods

Blood groups were tested in 271,164 subjects from 2014 to 2016. The ABO blood group was determined by agglutination using the antibodies anti-A, Anti-B, and Anti-D for the Rh factor, respectively.

Results

The overall distribution of ABO and Rh (D) groups in the population studied was as follows: O: 61.82%; A: 27.44%; B: 8.93%; and AB: 1.81%. For the Rh group, 95.58% of people were Rh (D), and 4.42% were Rh (d). Different distributions of blood groups across regions were found; additionally, genetic analysis revealed that the I^O and I^D allele showed an increasing trend from the north to the center, while the I^A and I^d allele tended to increase from the center to the north. Also, we found more gene diversity in both loci in the north compared with the center, suggesting population structure in Mexico.

Conclusion

This work could help health institutions to identify where they can obtain blood products necessary for medical interventions. Moreover, this piece of information contributes to the knowledge of the genetic structure of the Mexican populations which could have significant implications in different fields of biomedicine.

Human bony labyrinth is an indicator of population history and dispersal from Africa

The cavity system of the inner ear—the so-called bony labyrinth—houses the senses of balance and hearing. This structure is embedded in dense petrous bone, fully formed by birth and generally well preserved in human skeletal remains, thus providing a rich source of morphological information about past populations. Here we show that labyrinthine morphology tracks genetic distances and geography in an isolation-by-distance model with dispersal from Africa. Because petrous bones have become prime targets of ancient DNA recovery, we propose that all destructive studies first acquire high-resolution 3D computed-tomography data prior to any invasive sampling. Such data will constitute an important archive of morphological variation in past and present populations, and will permit individual-based genotype–phenotype comparisons.

The dispersal of modern humans from Africa is now well documented with genetic data that track population history, as well as gene flow between populations. Phenetic skeletal data, such as cranial and pelvic morphologies, also exhibit a dispersal-from-Africa signal, which, however, tends to be blurred by the effects of local adaptation and in vivo phenotypic plasticity, and that is often deteriorated by postmortem damage to skeletal remains. These complexities raise the question of which skeletal structures most effectively track neutral population history. The cavity system of the inner ear (the so-called bony labyrinth) is a good candidate structure for such analyses. It is already fully formed by birth, which minimizes postnatal phenotypic plasticity, and it is generally well preserved in archaeological samples. Here we use morphometric data of the bony labyrinth to show that it is a surprisingly good marker of the global dispersal of modern humans from Africa. Labyrinthine morphology tracks genetic distances and geography in accordance with an isolation-by-distance model with dispersal from Africa. Our data further indicate that the neutral-like pattern of variation is compatible with stabilizing selection on labyrinth morphology. Given the increasingly important role of the petrous bone for ancient DNA recovery from archaeological specimens, we encourage researchers to acquire 3D morphological data of the inner ear structures before any invasive sampling. Such data will constitute an important archive of phenotypic variation in present and past populations, and will permit individual-based genotype–phenotype comparisons.

A common copy-number variant within SIRPB1 correlates with human Out-of-Africa migration after genetic drift correction

Previous reports have proposed that personality may have played a role on human Out-Of-Africa migration, pinpointing some genetic variants that were positively selected in the migrating populations. In this work, we discuss the role of a common copy-number variant within the SIRPB1 gene, recently associated with impulsive behavior, in the human Out-Of-Africa migration. With the analysis of the variant distribution across forty-two different populations, we found that the SIRPB1 haplotype containing duplicated allele significantly correlated with human migratory distance, being one of the few examples of positively selected loci found across the human world colonization. Circular Chromosome Conformation Capture (4C-seq) experiments from the SIRPB1 promoter revealed important 3D modifications in the locus depending on the presence or absence of the duplication variant. In addition, a 3’ enhancer showed neural activity in transgenic models, suggesting that the presence of the CNV may compromise the expression of SIRPB1 in the central nervous system, paving the way to construct a molecular explanation of the SIRPB1 variants role in human migration.

Genetic and phenotypic variability of iris color in Buenos Aires population

The aim of this work was to describe the phenotypic and genotypic variability related to iris color for the population of Buenos Aires province (Argentina), and to assess the usefulness of current methods of analysis for this country. We studied five Single Nucleotide Polymorphisms (SNPs) included in the IrisPlex kit, in 118 individuals, and we quantified eye color with Digital Iris Analysis Tool. The markers fit Hardy-Weinberg equilibrium for the whole sample, but not for rs12913832 within the group of brown eyes (LR=8.429; p=0.004). We found a remarkable association of HERC2 rs12913832 GG with blue color (p < 0.01) but the other markers did not show any association with iris color. The results for the Buenos Aires population differ from those of other populations of the world for these polymorphisms (p < 0,01). The differences we found might respond to the admixed ethnic composition of Argentina; therefore, methods of analysis used in European populations should be carefully applied when studying the population of Argentina. These findings reaffirm the importance of this investigation in the Argentinian population for people identification based on iris color.

Genetic and linguistic histories in Central Asia inferred using approximate Bayesian computations

Linguistic and genetic data have been widely compared, but the histories underlying these descriptions are rarely jointly inferred. We developed a unique methodological framework for analysing jointly language diversity and genetic polymorphism data, to infer the past history of separation, exchange and admixture events among human populations. This method relies on approximate Bayesian computations that enable the identification of the most probable historical scenario underlying each type of data, and to infer the parameters of these scenarios. For this purpose, we developed a new computer program PopLingSim that simulates the evolution of linguistic diversity, which we coupled with an existing coalescent-based genetic simulation program, to simulate both linguistic and genetic data within a set of populations. Applying this new program to a wide linguistic and genetic dataset of Central Asia, we found several differences between linguistic and genetic histories. In particular, we showed how genetic and linguistic exchanges differed in the past in this area: some cultural exchanges were maintained without genetic exchanges. The methodological framework and the linguistic simulation tool developed here can be used in future work for disentangling complex linguistic and genetic evolutions underlying human biological and cultural histories.

Assessment of mitochondrial DNA heteroplasmy detected on commercial panel using MPS system with artificial mixture samples

Mitochondrial DNA (mtDNA) heteroplasmy is a potential genetic marker for forensic mtDNA analysis as well as phylogenic studies. Frequency of mtDNA heteroplasmy has been investigated in different populations through massively parallel sequencing (MPS) analysis, revealing various levels of frequency based on different MPS systems. For accurate heteroplasmy identification, it is essential to explore reliable detection threshold on various MPS systems. In addition, software solutions and pipelines need to be evaluated to analyze data effectively. In this study, heteroplasmy analysis was conducted on a commercially available mtDNA analysis system developed for forensic caseworks with artificially mixed DNA samples known for ratios and variant positions for assessment. mtDNA heteroplasmy > 10% was detectable with Torrent Variant Caller (TVC) while lower levels were identified using GeneMarker® HTS specialized software for minor variant detection. This study implies that analytical parameters and tools need to be optimized and evaluated for low-level heteroplasmy identification. Automated system with simple and efficient workflow is needed for forensic caseworks.

Worldwide patterns of human epigenetic variation

DNA methylation is an epigenetic modification, influenced by both genetic and environmental variation, that plays a key role in transcriptional regulation and many organismal phenotypes. Although patterns of DNA methylation have been shown to differ between human populations, it remains to be determined how epigenetic diversity relates to the patterns of genetic and gene expression variation at a global scale. Here we measured DNA methylation at 485,000 CpG sites in five diverse human populations, and analysed these data together with genome-wide genotype and gene expression data. We found that population-specific DNA methylation mirrors genetic variation, and has greater local genetic control than mRNA levels. We estimated the rate of epigenetic divergence between populations, which indicates far greater evolutionary stability of DNA methylation in humans than has been observed in plants. This study provides a deeper understanding of worldwide patterns of human epigenetic diversity, as well as initial estimates of the rate of epigenetic divergence in recent human evolution.

Mutation Spectrum of STAR and a Founder Effect of the p.Q258* in Korean Patients with Congenital Lipoid Adrenal Hyperplasia

Congenital lipoid adrenal hyperplasia (CLAH) is the most severe form of congenital adrenal hyperplasia, caused by defects in the steroidogenic acute regulatory protein (STAR). The STAR p.Q258* mutation is the most common mutation in China, Japan, and Korea, suggesting a founder effect. This study aimed to investigate the phenotypic and mutation spectrum of STAR defects and identify a founder effect of the p.Q258* mutation in Korean patients with CLAH. For 45 patients from 42 independent pedigrees, haplotype analysis was performed in 10 unrelated trio families, including patients with the p.Q258* mutation whose DNA samples were available, using 1,972 single nucleotide polymorphism (SNP) and six short tandem repeat (STR) markers. An Illumina Infinium® Human Omni2.5-8 v1.3 performed the SNP genotyping. Among 84 alleles from 42 unrelated families, mutation p.Q258* was found in 74 alleles (88.1%) from 41 families. A shared haplotype was identified in 17 of 20 alleles from 10 patients (size, 198 kb). The age of the founder mutation was estimated as 4,875 years (95% credible set: 3,575-7,925 years) assuming an intergenerational time interval of 25 years. The STAR p.Q258* mutation is the most common in Korean patients with CLAH, suggesting a founder effect. The age of the mutation corresponded with the date when the Korean people settled in the Korean peninsula. The high prevalence of p.Q258* in Japan and China also suggests a founder effect in Asian countries.

Intracellular Nucleic Acid Detection in Autoimmunity

Protective immune responses to viral infection are initiated by innate immune sensors that survey extracellular and intracellular space for foreign nucleic acids. The existence of these sensors raises fundamental questions about self/nonself discrimination because of the abundance of self-DNA and self-RNA that occupy these same compartments. Recent advances have revealed that enzymes that metabolize or modify endogenous nucleic acids are essential for preventing inappropriate activation of the innate antiviral response. In this review, we discuss rare human diseases caused by dysregulated nucleic acid sensing, focusing primarily on intracellular sensors of nucleic acids. We summarize lessons learned from these disorders, we rationalize the existence of these diseases in the context of evolution, and we propose that this framework may also apply to a number of more common autoimmune diseases for which the underlying genetics and mechanisms are not yet fully understood.