PGG.SNV: understanding the evolutionary and medical implications of human single nucleotide variations in diverse populations

Novel genetic association of the Furin gene polymorphism rs1981458 with COVID-19 severity among Indian populations

SARS CoV-2, the causative agent for the ongoing COVID-19 pandemic, it enters the host cell by activating the ACE2 receptor with the help of two proteasesi.e., Furin and TMPRSS2. Therefore, variations in these genes may account for differential susceptibility and severity between populations. Previous studies have shown that the role of ACE2 and TMPRSS2 gene variants in understanding COVID-19 susceptibility among Indian populations. Nevertheless, a knowledge gap exists concerning the COVID-19 susceptibility of Furin gene variants among diverse South Asian ethnic groups. Investigating the role of Furin gene variants and their global phylogeographic structure is essential to comprehensively understanding COVID-19 susceptibility in these populations. We have used 450 samples from diverse Indian states and performed linear regression to analyse the Furin gene variant's with COVID-19 Case Fatality Rate (CFR) that could be epidemiologically associated with disease severity outcomes. Associated genetic variants were further evaluated for their expression and regulatory potential through various Insilco analyses. Additionally, we examined the Furin gene using next-generation sequencing (NGS) data from 393 diverse global samples, with a particular emphasis on South Asia, to investigate its Phylogeographic structure among diverse world populations. We found a significant positive association for the SNP rs1981458 with COVID-19 CFR (p < 0.05) among diverse Indian populations at different timelines of the first and second waves. Further, QTL and other regulatory analyses showed various significant associations for positive regulatory roles of rs1981458 and Furin gene, mainly in Immune cells and virus infection process, highlighting their role in host immunity and viral assembly and processing. The Furin protein-protein interaction suggested that COVID-19 may contribute to Pulmonary arterial hypertension via a typical inflammation mechanism. The phylogeographic architecture of the Furin gene demonstrated a closer genetic affinity of South Asia with West Eurasian populations. Therefore, it is worth proposing that for the Furin gene, the COVID-19 susceptibility of South Asians will be more similar to the West Eurasian population. Our previous studies on the ACE2 and TMPRSS2 genes showed genetic affinity of South Asian with East Eurasians and West Eurasians, respectively. Therefore, with the collective information from these three important genes (ACE2, TMPRSS2 and Furin) we modelled COVID-19 susceptibilityof South Asia in between these two major ancestries with an inclination towards West Eurasia. In conclusion, this study, for the first time, concluded the role of rs1981458 in COVID-19 severity among the Indian population and outlined its regulatory potential.This study also highlights that the genetic structure for COVID-19 susceptibilityof South Asia is distinct, however, inclined to the West Eurasian population. We believe this insight may be utilised as a genetic biomarker to identify vulnerable populations, which might be directly relevant for developing policies and allocating resources more effectively during an epidemic.

Reconstructing the ancestral gene pool to uncover the origins and genetic links of Hmong-Mien speakers

BACKGROUND

Hmong-Mien (HM) speakers are linguistically related and live primarily in China, but little is known about their ancestral origins or the evolutionary mechanism shaping their genomic diversity. In particular, the lack of whole-genome sequencing data on the Yao population has prevented a full investigation of the origins and evolutionary history of HM speakers. As such, their origins are debatable.

RESULTS

Here, we made a deep sequencing effort of 80 Yao genomes, and our analysis together with 28 East Asian populations and 968 ancient Asian genomes suggested that there is a strong genetic basis for the formation of the HM language family. We estimated that the most recent common ancestor dates to 5800 years ago, while the genetic divergence between the HM and Tai-Kadai speakers was estimated to be 8200 years ago. We proposed that HM speakers originated from the Yangtze River Basin and spread with agricultural civilization. We identified highly differentiated variants between HM and Han Chinese, in particular, a deafness-related missense variant (rs72474224) in the GJB2 gene is in a higher frequency in HM speakers than in others.

CONCLUSIONS

Our results indicated complex gene flow and medically relevant variants involved in the HM speakers' evolution history.

HervD Atlas: a curated knowledgebase of associations between human endogenous retroviruses and diseases

Human endogenous retroviruses (HERVs), as remnants of ancient exogenous retrovirus infected and integrated into germ cells, comprise ∼8% of the human genome. These HERVs have been implicated in numerous diseases, and extensive research has been conducted to uncover their specific roles. Despite these efforts, a comprehensive source of HERV-disease association still needs to be added. To address this gap, we introduce the HervD Atlas (https://ngdc.cncb.ac.cn/hervd/), an integrated knowledgebase of HERV-disease associations manually curated from all related published literature. In the current version, HervD Atlas collects 60 726 HERV-disease associations from 254 publications (out of 4692 screened literature), covering 21 790 HERVs (21 049 HERV-Terms and 741 HERV-Elements) belonging to six types, 149 diseases and 610 related/affected genes. Notably, an interactive knowledge graph that systematically integrates all the HERV-disease associations and corresponding affected genes into a comprehensive network provides a powerful tool to uncover and deduce the complex interplay between HERVs and diseases. The HervD Atlas also features a user-friendly web interface that allows efficient browsing, searching, and downloading of all association information, research metadata, and annotation information. Overall, the HervD Atlas is an essential resource for comprehensive, up-to-date knowledge on HERV-disease research, potentially facilitating the development of novel HERV-associated diagnostic and therapeutic strategies.

Genetic and cultural adaptations underlie the establishment of dairy pastoralism in the Tibetan Plateau

BACKGROUND

Domestication and introduction of dairy animals facilitated the permanent human occupation of the Tibetan Plateau. Yet the history of dairy pastoralism in the Tibetan Plateau remains poorly understood. Little is known how Tibetans adapted to milk and dairy products.

RESULTS

We integrated archeological evidence and genetic analysis to show the picture that the dairy ruminants, together with dogs, were introduced from West Eurasia into the Tibetan Plateau since ~ 3600 years ago. The genetic admixture between the exotic and indigenous dogs enriched the candidate lactase persistence (LP) allele 10974A > G of West Eurasian origin in Tibetan dogs. In vitro experiments demonstrate that - 13838G > A functions as a LP allele in Tibetans. Unlike multiple LP alleles presenting selective signatures in West Eurasians and South Asians, the de novo origin of Tibetan-specific LP allele - 13838G > A with low frequency (~ 6-7%) and absence of selection corresponds - 13910C > T in pastoralists across eastern Eurasia steppe.

CONCLUSIONS

Results depict a novel scenario of genetic and cultural adaptations to diet and expand current understanding of the establishment of dairy pastoralism in the Tibetan Plateau.

Historical Pathogen-Driven Selection May Contribute to Contemporary Ethnic Difference in Bladder Cancer Susceptibility

BACKGROUND

The rationale for ethnic differences in bladder cancer (BCa) susceptibility is an important open question. In this study, we raised the hypothesis that the APOBEC3-rs1014971 variant associated with BCa risk and APOBEC-mutagenesis probably contribute to ethnic differences.

METHODS

We calculated the ethnicity-stratified 5-year age-adjusted incidence rates of BCa using the US SEER database. We performed somatic mutational-signature analyses and compared the APOBEC-related mutational contribution across BCa tumors in patients of different ethnicities. We analyzed the allele frequency distribution of APOBEC3-related rs1014971 in contemporary populations of different ethnicities and in ancient human genomes. We also analyzed the natural selection profiles and ages of the investigated SNPs.

RESULTS

We validated the ethnic difference in BCa risk using US SEER data, revealing Caucasians to be at >2-fold greater risk than Asians / Pacific islanders. In contemporary populations, we observed a coherent ethnic distribution in terms not only of the allele frequency of APOBEC3-related rs1014971, but also the mutational contribution of APOBEC-mediated mutagenesis in BCa tumors. Population genetics and ancient genome analyses further suggested that the diverse ethnic distribution of rs1014971 could be rooted in human evolution.

CONCLUSIONS

It is possible that APOBEC3-related rs1014971 is involved in the different BCa incidence across ethnic groups, and this difference is potentially derived from human evolution. Our findings suggested an evolutionary link between contemporary population-level variations in malignancy susceptibility and pathogen-driven selection in the past, not unlike previously reported cases of certain autoimmune and metabolic disorders.

Genome-wide association study on coordination and agility in 461 Chinese Han males

There is growing evidence that genetic factors can influence human athletic performance. In many sports performances, excellent coordination and agility are the keys to mastery. However, few studies have been devoted to identifying genetic influences on athletic performance. Methods: We generated a derived measure of coordination and agility from the data of hexagonal jumps and T-runs and conducted genome-wide association and meta-analysis studies focused on coordination and agility. Results: The phenotypic correlation and genetic covariance analysis indicated that hexagonal jumps and T-runs were possibly influenced by the same set of genetic factors (R = 0.27, genetic covariance = 0.59). Meta-analysis identified rs117047321 genome-wide significant association (N = 143, P < 10E-5) with coordination and agility, and this association was replicated in the replication group (N = 318, P < 0.05). The CG genotype samples of this single nucleotide polymorphism (SNP) required a longer average movement time than the CC genotype samples, and the CG genotype only exists in Asia, which may belong to the East Asia-specific variation. This SNP is located on MYO5B, which is highly expressed in tissues such as the brain, heart, and muscle, suggesting that this locus might be a genetic factor related to human energy metabolism. Conclusion: Our study indicated that genetic factors can affect the athletic performance of coordination and agility. These findings may provide valuable insights for using genetic factors to evaluate sports characteristics.

A comprehensive meta-analysis comprising 149 case-control studies to investigate the association between IL-6 gene rs1800795 polymorphism and multiple disease risk

BACKGROUND

Individual genome-wide association studies (GWAS) or single case-specific meta-analyses may not be sufficient evidence to take action against a specific gene function. Thus, we tried to determine a consensus association between the IL-6 gene rs1800795 polymorphism and multiple disease risks through an updated statistical meta-analysis.

METHOD

After systematically searching online databases, we found 149 case-control relevant datasets with a sample size of 96,153 (cases: 38,291 and controls: 57862) and conducted the meta-analysis using updated statistical models.

RESULTS

The analyses of this comprehensive meta-analysis revealed a significant association between IL-6 -174G/C polymorphism and overall disorder risk under all genetic models (C vs G: OR = 1.11, 95% CI = 1.08-1.13; p-value = 4.8E-17; CC vs GG: OR = 1.19, 95% CI = 1.13-1.26; p-value = 9.4E-12; CG vs GG: OR = 1.10, 95% CI = 1.06-1.14; p-value = 1.1E-07; CC + CG vs GG: OR = 1.13, 95% CI = 1.10-1.17; p-value = 1.1E-13; CC vs CG + GG: OR = 1.18, 95% CI = 1.06-1.31; p-value = 0.0019) and (OR > 1) with Asian ethnicity. The subgroup analyses based on the diseases revealed that the polymorphism was highly significantly increasing the risk of coronary artery disease (CAD) under all genetic models. Likewise, a significant association was observed with increased risk under three genetic models of inflammatory diseases (C vs G; CC vs GG; and CC vs CG + GG), and rheumatoid arthritis (C vs G; CG vs GG; and CC + CG vs GG). Conversely, the -174G/C SNP significantly decreased the risk of ischemic stroke under the two genetic models (C vs G; and CG vs GG). However, the other diseases included in this study showed no significant association with IL-6 (-174G/C) polymorphism.

CONCLUSION

This meta-analysis provided strong evidence for the association between IL-6 gene rs1800795 polymorphism and multiple disease risks. The IL-6 gene could be a useful prognostic biomarker for CAD, inflammatory disease, ischemic stroke, and rheumatoid arthritis.

Expression profiles of east-west highly differentiated genes in Uyghur genomes

It remains unknown and debatable how European-Asian-differentiated alleles affect individual phenotypes. Here, we made the first effort to analyze the expression profiles of highly differentiated genes with eastern and western origins in 90 Uyghurs using whole-genome (30× to 60×) and transcriptome data. We screened 921 872 east-west highly differentiated genetic variants, of which ∼4.32% were expression quantitative trait loci (eQTLs), ∼0.12% were alternative splicing quantitative trait loci (sQTLs), and ∼0.12% showed allele-specific expression (ASE). The 8305 highly differentiated eQTLs of strong effects appear to have undergone natural selection, associated with immunity and metabolism. European-origin alleles tend to be more biasedly expressed; highly differentiated ASEs were enriched in diabetes-associated genes, likely affecting the diabetes susceptibility in the Uyghurs. We proposed an admixture-induced expression model to dissect the highly differentiated expression profiles. We provide new insights into the genetic basis of phenotypic differentiation between Western and Eastern populations, advancing our understanding of the impact of genetic admixture.

PGG.MHC: toward understanding the diversity of major histocompatibility complexes in human populations

The human leukocyte antigen (HLA) system, or the human version of the major histocompatibility complex (MHC), is known for its extreme polymorphic nature and high heterogeneity. Taking advantage of whole-genome and whole-exome sequencing data, we developed PGG.MHC to provide a platform to explore the diversity of the MHC in Asia as well as in global populations. PGG.MHC currently archives high-resolution HLA alleles of 53 254 samples representing 190 populations spanning 66 countries. PGG.MHC provides: (i) high-quality allele frequencies for eight classical HLA loci (HLA-A, -B, -C, -DQA1, -DQB1, -DRB1, -DPA1 and -DPB1); (ii) visualization of population prevalence of HLA alleles on global, regional, and country-wide levels; (iii) haplotype structure of 134 populations; (iv) two online analysis tools including 'HLA imputation' for inferring HLA alleles from SNP genotyping data and 'HLA association' to perform case/control studies for HLA-related phenotypes and (v) East Asian-specific reference panels for HLA imputation. Equipped with high-quality frequency data and user-friendly computer tools, we expect that the PGG.MHC database can advance the understanding and facilitate applications of MHC genomic diversity in both evolutionary and medical studies. The PGG.MHC database is freely accessible via https://pog.fudan.edu.cn/pggmhc or https://www.pggmhc.org/pggmhc.

PGG.SV: a whole-genome-sequencing-based structural variant resource and data analysis platform

Structural variations (SVs) play important roles in human evolution and diseases, but there is a lack of data resources concerning representative samples, especially for East Asians. Taking advantage of both next-generation sequencing and third-generation sequencing data at the whole-genome level, we developed the database PGG.SV to provide a practical platform for both regionally and globally representative structural variants. In its current version, PGG.SV archives 584 277 SVs obtained from whole-genome sequencing data of 6048 samples, including 1030 long-read sequencing genomes representing 177 global populations. PGG.SV provides (i) high-quality SVs with fine-scale and precise genomic locations in both GRCh37 and GRCh38, covering underrepresented SVs in existing sequencing and microarray data; (ii) hierarchical estimation of SV prevalence in geographical populations; (iii) informative annotations of SV-related genes, potential functions and clinical effects; (iv) an analysis platform to facilitate SV-based case-control association studies and (v) various visualization tools for understanding the SV structures in the human genome. Taken together, PGG.SV provides a user-friendly online interface, easy-to-use analysis tools and a detailed presentation of results. PGG.SV is freely accessible via https://www.biosino.org/pggsv.

Genomic data integration and user-defined sample-set extraction for population variant analysis

BACKGROUND

Population variant analysis is of great importance for gathering insights into the links between human genotype and phenotype. The 1000 Genomes Project established a valuable reference for human genetic variation; however, the integrative use of the corresponding data with other datasets within existing repositories and pipelines is not fully supported. Particularly, there is a pressing need for flexible and fast selection of population partitions based on their variant and metadata-related characteristics.

RESULTS

Here, we target general germline or somatic mutation data sources for their seamless inclusion within an interoperable-format repository, supporting integration among them and with other genomic data, as well as their integrated use within bioinformatic workflows. In addition, we provide VarSum, a data summarization service working on sub-populations of interest selected using filters on population metadata and/or variant characteristics. The service is developed as an optimized computational framework with an Application Programming Interface (API) that can be called from within any existing computing pipeline or programming script. Provided example use cases of biological interest show the relevance, power and ease of use of the API functionalities.

CONCLUSIONS

The proposed data integration pipeline and data set extraction and summarization API pave the way for solid computational infrastructures that quickly process cumbersome variation data, and allow biologists and bioinformaticians to easily perform scalable analysis on user-defined partitions of large cohorts from increasingly available genetic variation studies. With the current tendency to large (cross)nation-wide sequencing and variation initiatives, we expect an ever growing need for the kind of computational support hereby proposed.

Tracing Bai-Yue Ancestry in Aboriginal Li People on Hainan Island

As the most prevalent aboriginal group on Hainan Island located between South China and the mainland of Southeast Asia, the Li people are believed to preserve some unique genetic information due to their isolated circumstances, although this has been largely uninvestigated. We performed the first whole-genome sequencing of 55 Hainan Li (HNL) individuals with high coverage (∼30-50×) to gain insight into their genetic history and potential adaptations. We identified the ancestry enriched in HNL (∼85%) is well preserved in present-day Tai-Kadai speakers residing in South China and North Vietnam, that is, Bai-Yue populations. A lack of admixture signature due to the geographical restriction exacerbated the bottleneck in the present-day HNL. The genetic divergence among Bai-Yue populations began ∼4,000-3,000 years ago when the proto-HNL underwent migration and the settling of Hainan Island. Finally, we identified signatures of positive selection in the HNL, some outstanding examples included FADS1 and FADS2 related to a diet rich in polyunsaturated fatty acids. In addition, we observed that malaria-driven selection had occurred in the HNL, with population-specific variants of malaria-related genes (e.g., CR1) present. Interestingly, HNL harbors a high prevalence of malaria leveraged gene variants related to hematopoietic function (e.g., CD3G) that may explain the high incidence of blood disorders such as B-cell lymphomas in the present-day HNL. The results have advanced our understanding of the genetic history of the Bai-Yue populations and have provided new insights into the adaptive scenarios of the Li people.

Genomic diversity and post-admixture adaptation in the Uyghurs

Population admixture results in genome-wide combinations of genetic variants derived from different ancestral populations of distinct ancestry, thus providing a unique opportunity for understanding the genetic determinants of phenotypic variation in humans. Here, we used whole-genome sequencing of 92 individuals with high coverage (30–60×) to systematically investigate genomic diversity in the Uyghurs living in Xinjiang, China (XJU), an admixed population of both European-like and East-Asian-like ancestry. The XJU population shows greater genetic diversity, especially a higher proportion of rare variants, compared with their ancestral source populations, corresponding to greater phenotypic diversity of XJU. Admixture-induced functional variants in EDAR were associated with the diversity of facial morphology in XJU. Interestingly, the interaction of functional variants between SLC24A5 and OCA2 likely influences the diversity of skin pigmentation. Notably, selection has seemingly been relaxed or canceled in several genes with significantly biased ancestry, such as HERC2–OCA2. Moreover, signatures of post-admixture adaptation in XJU were identified, including genes related to metabolism (e.g. CYP2D6), digestion (e.g. COL11A1), olfactory perception (e.g. ANO2) and immunity (e.g. HLA). Our results demonstrated population admixture as a driving force, locally or globally, in shaping human genetic and phenotypic diversity as well as in adaptive evolution.

Genetic association of TMPRSS2 rs2070788 polymorphism with COVID-19 case fatality rate among Indian populations

SARS-CoV-2, the causative agent for COVID-19, an ongoing pandemic, engages the ACE2 receptor to enter the host cell through S protein priming by a serine protease, TMPRSS2. Variation in the TMPRSS2 gene may account for the disparity in disease susceptibility between populations. Therefore, in the present study, we have used next-generation sequencing (NGS) data of world populations from 393 individuals and analyzed the TMPRSS2 gene using a haplotype-based approach with a major focus on South Asia to study its phylogenetic structure and their haplotype sharing among various populations worldwide. Our analysis of phylogenetic relatedness showed a closer affinity of South Asians with the West Eurasian populations therefore, host disease susceptibility and severity particularly in the context of TMPRSS2 will be more akin to West Eurasian instead of East Eurasian. This is in contrast to our prior study on the ACE2 gene which shows South Asian haplotypes have a strong affinity towards West Eurasians. Thus ACE2 and TMPRSS2 have an antagonistic genetic relatedness among South Asians. Considering the significance of the TMPRSS2 gene in the SARS-CoV-2 pathogenicity, COVID-19 infection and intensity trends could be directly associated with increased expression therefore, we have also tested the SNPs frequencies of this gene among various Indian state populations with respect to the case fatality rate (CFR). Interestingly, we found a significant positive association between the rs2070788 SNP (G Allele) and the CFR among Indian populations. Further our cis eQTL analysis of rs2070788 shows that the GG genotype of the rs2070788 tends to have a significantly higher expression of TMPRSS2 gene in the lung compared to the AG and AA genotypes thus validating the previous observation and therefore it might play a vital part in determining differential disease vulnerability. We trust that this information will be useful in understanding the role of the TMPRSS2 variant in COVID-19 susceptibility and using it as a biomarker may help to predict populations at risk.

Advances in integrative African genomics

There has been a rapid increase in human genome sequencing in the past two decades, resulting in the identification of millions of previously unknown genetic variants. However, African populations are under-represented in sequencing efforts. Additional sequencing from diverse African populations and the construction of African-specific reference genomes is needed to better characterize the full spectrum of variation in humans. However, sequencing alone is insufficient to address the molecular and cellular mechanisms underlying variable phenotypes and disease risks. Determining functional consequences of genetic variation using multi-omics approaches is a fundamental post-genomic challenge. We discuss approaches to close the knowledge gaps about African genomic diversity and review advances in African integrative genomic studies and their implications for precision medicine.

Regional TMPRSS2 V197M Allele Frequencies Are Correlated with COVID-19 Case Fatality Rates

Coronavirus disease, COVID-19 (coronavirus disease 2019), caused by SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2), has a higher case fatality rate in European countries than in others, especially East Asian ones. One potential explanation for this regional difference is the diversity of the viral infection efficiency. Here, we analyzed the allele frequencies of a nonsynonymous variant rs12329760 (V197M) in the TMPRSS2 gene, a key enzyme essential for viral infection and found a significant association between the COVID-19 case fatality rate and the V197M allele frequencies, using over 200,000 present-day and ancient genomic samples. East Asian countries have higher V197M allele frequencies than other regions, including European countries which correlates to their lower case fatality rates. Structural and energy calculation analysis of the V197M amino acid change showed that it destabilizes the TMPRSS2 protein, possibly negatively affecting its ACE2 and viral spike protein processing.

Fine-Scale Genetic Structure and Natural Selection Signatures of Southwestern Hans Inferred From Patterns of Genome-Wide Allele, Haplotype, and Haplogroup Lineages

The evolutionary and admixture history of Han Chinese have been widely discussed via traditional autosomal and uniparental genetic markers [e.g., short tandem repeats, low-density single nucleotide polymorphisms). However, their fine-scale genetic landscapes (admixture scenarios and natural selection signatures) based on the high-density allele/haplotype sharing patterns have not been deeply characterized. Here, we collected and generated genome-wide data of 50 Han Chinese individuals from four populations in Guizhou Province, one of the most ethnolinguistically diverse regions, and merged it with over 3,000 publicly available modern and ancient Eurasians to describe the genetic origin and population admixture history of Guizhou Hans and their neighbors. PCA and ADMIXTURE results showed that the studied four populations were homogeneous and grouped closely to central East Asians. Genetic homogeneity within Guizhou populations was further confirmed via the observed strong genetic affinity with inland Hmong-Mien people through the observed genetic clade in Fst and outgroup f₃/f₄-statistics. qpGraph-based phylogenies and f₄-based demographic models illuminated that Guizhou Hans were well fitted via the admixture of ancient Yellow River Millet farmers related to Lajia people and southern Yangtze River farmers related to Hanben people. Further ChromoPainter-based chromosome painting profiles and GLOBETROTTER-based admixture signatures confirmed the two best source matches for southwestern Hans, respectively, from northern Shaanxi Hans and southern indigenes with variable mixture proportions in the historical period. Further three-way admixture models revealed larger genetic contributions from coastal southern East Asians into Guizhou Hans compared with the proposed inland ancient source from mainland Southeast Asia. We also identified candidate loci (e.g., MTUS2, NOTCH4, EDAR, ADH1B, and ABCG2) with strong natural selection signatures in Guizhou Hans via iHS, nSL, and ihh, which were associated with the susceptibility of the multiple complex diseases, morphology formation, alcohol and lipid metabolism. Generally, we provided a case and ideal strategy to reconstruct the detailed demographic evolutionary history of Guizhou Hans, which provided new insights into the fine-scale genomic formation of one ethnolinguistically specific targeted population from the comprehensive perspectives of the shared unlinked alleles, linked haplotypes, and paternal and maternal lineages.

False Alarms in Consumer Genomics Add to Public Fear and Potential Health Care Burden

We have entered an era of direct-to-consumer (DTC) genomics. Patients have relayed many success stories of DTC genomics about finding causal mutations of genetic diseases before showing any symptoms and taking precautions. However, consumers may also take unnecessary medical actions based on false alarms of “pathogenic alleles”. The severity of this problem is not well known. Using publicly available data, we compared DTC microarray genotyping data with deep-sequencing data of 5 individuals and manually checked each inconsistently reported single nucleotide variants (SNVs). We estimated that, on average, a person would have ~5 “pathogenic” alleles reported due to wrongly reported genotypes if using a 23andMe genotyping microarray. We also found that the number of wrongly classified “pathogenic” alleles per person is at least as significant as those due to wrongly reported genotypes. We show that the scale of the false alarm problem could be large enough that the medical costs will become a burden to public health.

Triose Kinase Controls the Lipogenic Potential of Fructose and Dietary Tolerance

The surge in fructose consumption is a major factor behind the rapid rise of nonalcoholic fatty liver disease in modern society. Through flux and genetic analyses, we demonstrate that fructose is catabolized at a much higher rate than glucose, and triose kinase (TK) couples fructolysis with lipogenesis metabolically and transcriptionally. In the absence of TK, fructose oxidation is accelerated through the activation of aldehyde dehydrogenase (ALDH) and serine biosynthesis, accompanied by increased oxidative stress and fructose aversion. TK is also required by the endogenous fructolysis pathway to drive lipogenesis and hepatic triglyceride accumulation under high-fat diet and leptin-deficient conditions. Intriguingly, a nonsynonymous TK allele (rs2260655_A) segregated during human migration out of Africa behaves as TK null for its inability to rescue fructose toxicity and increase hepatic triglyceride accumulation. Therefore, we posit TK as a metabolic switch controlling the lipogenic potential of fructose and its dietary tolerance.

Most frequent South Asian haplotypes of ACE2 share identity by descent with East Eurasian populations

It was shown that the human Angiotensin-converting enzyme 2 (ACE2) is the receptor of recent coronavirus SARS-CoV-2, and variation in this gene may affect the susceptibility of a population. Therefore, we have analysed the sequence data of ACE2 among 393 samples worldwide, focusing on South Asia. Genetically, South Asians are more related to West Eurasian populations rather than to East Eurasians. In the present analyses of ACE2, we observed that the majority of South Asian haplotypes are closer to East Eurasians rather than to West Eurasians. The phylogenetic analysis suggested that the South Asian haplotypes shared with East Eurasians involved two unique event polymorphisms (rs4646120 and rs2285666). In contrast with the European/American populations, both of the SNPs have largely similar frequencies for East Eurasians and South Asians, Therefore, it is likely that among the South Asians, host susceptibility to the novel coronavirus SARS-CoV-2 will be more similar to that of East Eurasians rather than to that of Europeans.

Indigenous Genomic Databases: Pragmatic Considerations and Cultural Contexts

The potential to grow genomic knowledge and harness the subsequent clinical benefits has escalated the building of background variant databases (BVDs) for genetic diagnosis across the globe. Alongside the upsurge of this precision medicine, potential benefits have been highlighted for both rare genetic conditions and other diagnoses. However, with the ever-present “genomic divide,” Indigenous peoples globally have valid concerns as they endure comparatively greater health disparities but stand to benefit the least from these novel scientific discoveries and progress in healthcare. The paucity of Indigenous healthcare providers and researchers in these fields contributes to this genomic divide both in access to, and availability of culturally safe, relevant and respectful healthcare using this genetic knowledge. The vital quest to provide equitable clinical research, and provision and use of genomic services and technologies provides a strong rationale for building BVDs for Indigenous peoples. Such tools would ground their representation and participation in accompanying genomic health research and benefit acquisition. We describe two, independent but highly similar initiatives–the “Silent Genomes” in Canada and the “Aotearoa Variome” in New Zealand–as exemplars that have had to address the aforementioned issues and work to create Indigenous BVDs with these populations. Taking into account the baseline inequities in genomic medicine for Indigenous populations and the ongoing challenges of implementing genomic research with Indigenous communities, we provide a rationale for multiple changes required that will assure communities represented in BVDs, as well as Indigenous researchers, that their participation will maximize benefits and minimize risk.

Database Resources of the National Genomics Data Center in 2020

The National Genomics Data Center (NGDC) provides a suite of database resources to support worldwide research activities in both academia and industry. With the rapid advancements in higher-throughput and lower-cost sequencing technologies and accordingly the huge volume of multi-omics data generated at exponential scales and rates, NGDC is continually expanding, updating and enriching its core database resources through big data integration and value-added curation. In the past year, efforts for update have been mainly devoted to BioProject, BioSample, GSA, GWH, GVM, NONCODE, LncBook, EWAS Atlas and IC4R. Newly released resources include three human genome databases (PGG.SNV, PGG.Han and CGVD), eLMSG, EWAS Data Hub, GWAS Atlas, iSheep and PADS Arsenal. In addition, four web services, namely, eGPS Cloud, BIG Search, BIG Submission and BIG SSO, have been significantly improved and enhanced. All of these resources along with their services are publicly accessible at https://bigd.big.ac.cn.