The GenomeAsia 100K Project enables genetic discoveries across Asia

Association between single-nucleotide polymorphisms in miRNA and breast cancer risk: an updated review

Breast cancer (BC), a heterogeneous, aggressive illness with high mortality, is essentially a genomic disease. While the high-penetrance genes BRCA1 and BRCA2 play important roles in tumorigenesis, moderate- and low-penetrance genes are also involved. Single-nucleotide polymorphisms (SNPs) in microRNA (miRNA) genes have recently been identified as BC risk factors. miRNA genes are currently classified as low-penetrance. SNPs are the most common variations in the human genome. While the role of miRNA SNPs in BC susceptibility has been studied extensively, results have been inconsistent. This review analyzes the results of association studies between miRNA SNPs and BC risk from countries around the world. We conclude that: (a) By continent, the largest proportion of studies to date were conducted in Asia (65.0 %) and the smallest proportion in Africa (1.8 %); (b) Association studies have been completed for 67 different SNPs; (c) 146a, 196a2, 499, 27a, and 423 are the most-studied miRNAs; (d) The SNPs rs2910164 (miRNA-146a), rs11614913 (miRNA-196a2), rs3746444 (miRNA-499) and rs6505162 (miRNA-423) were the most widely associated with increased BC risk; (e) The majority of studies had small samples, which may affect the precision and power of the results; and (f) The effect of an SNP on BC risk depends on the ethnicity of the population. This review also discusses potential explanations for controversial findings.

The genomic history of the Middle East

The Middle East region is important to understand human evolution and migrations but is underrepresented in genomic studies. Here, we generated 137 high-coverage physically phased genome sequences from eight Middle Eastern populations using linked-read sequencing. We found no genetic traces of early expansions out-of-Africa in present-day populations but found Arabians have elevated Basal Eurasian ancestry that dilutes their Neanderthal ancestry. Population sizes within the region started diverging 15–20 kya, when Levantines expanded while Arabians maintained smaller populations that derived ancestry from local hunter-gatherers. Arabians suffered a population bottleneck around the aridification of Arabia 6 kya, while Levantines had a distinct bottleneck overlapping the 4.2 kya aridification event. We found an association between movement and admixture of populations in the region and the spread of Semitic languages. Finally, we identify variants that show evidence of selection, including polygenic selection. Our results provide detailed insights into the genomic and selective histories of the Middle East.

•

Middle Easterners do not have ancestry from an early out-of-Africa expansion

•

Basal Eurasian and African ancestry in Arabians deplete their Neanderthal ancestry

•

Populations experienced bottlenecks overlapping aridification events

•

Identification of recent single and polygenic signals of selection in Arabia

Philippine Ayta possess the highest level of Denisovan ancestry in the world

Multiple lines of evidence show that modern humans interbred with archaic Denisovans. Here, we report an account of shared demographic history between Australasians and Denisovans distinctively in Island Southeast Asia. Our analyses are based on ∼2.3 million genotypes from 118 ethnic groups of the Philippines, including 25 diverse self-identified Negrito populations, along with high-coverage genomes of Australopapuans and Ayta Magbukon Negritos. We show that Ayta Magbukon possess the highest level of Denisovan ancestry in the world-∼30%-40% greater than that of Australians and Papuans-consistent with an independent admixture event into Negritos from Denisovans. Together with the recently described Homo luzonensis, we suggest that there were multiple archaic species that inhabited the Philippines prior to the arrival of modern humans and that these archaic groups may have been genetically related. Altogether, our findings unveil a complex intertwined history of modern and archaic humans in the Asia-Pacific region, where distinct Islander Denisovan populations differentially admixed with incoming Australasians across multiple locations and at various points in time.

Blood groups of Neandertals and Denisova decrypted

Blood group systems were the first phenotypic markers used in anthropology to decipher the origin of populations, their migratory movements, and their admixture. The recent emergence of new technologies based on the decoding of nucleic acids from an individual’s entire genome has relegated them to their primary application, blood transfusion. Thus, despite the finer mapping of the modern human genome in relation to Neanderthal and Denisova populations, little is known about red cell blood groups in these archaic populations. Here we analyze the available high-quality sequences of three Neanderthals and one Denisovan individuals for 7 blood group systems that are used today in transfusion (ABO including H/Se, Rh (Rhesus), Kell, Duffy, Kidd, MNS, Diego). We show that Neanderthal and Denisova were polymorphic for ABO and shared blood group alleles recurrent in modern Sub-Saharan populations. Furthermore, we found ABO-related alleles currently preventing from viral gut infection and Neanderthal RHD and RHCE alleles nowadays associated with a high risk of hemolytic disease of the fetus and newborn. Such a common blood group pattern across time and space is coherent with a Neanderthal population of low genetic diversity exposed to low reproductive success and with their inevitable demise. Lastly, we connect a Neanderthal RHD allele to two present-day Aboriginal Australian and Papuan, suggesting that a segment of archaic genome was introgressed in this gene in non-Eurasian populations. While contributing to both the origin and late evolutionary history of Neanderthal and Denisova, our results further illustrate that blood group systems are a relevant piece of the puzzle helping to decipher it.

Discovery of genomic variation across a generation

Over the past 30 years (the timespan of a generation), advances in genomics technologies have revealed tremendous and unexpected variation in the human genome and have provided increasingly accurate answers to long-standing questions of how much genetic variation exists in human populations and to what degree the DNA complement changes between parents and offspring. Tracking the characteristics of these inherited and spontaneous (or de novo) variations has been the basis of the study of human genetic disease. From genome-wide microarray and next-generation sequencing scans, we now know that each human genome contains over 3 million single nucleotide variants when compared with the ~ 3 billion base pairs in the human reference genome, along with roughly an order of magnitude more DNA—approximately 30 megabase pairs (Mb)—being ‘structurally variable’, mostly in the form of indels and copy number changes. Additional large-scale variations include balanced inversions (average of 18 Mb) and complex, difficult-to-resolve alterations. Collectively, ~1% of an individual’s genome will differ from the human reference sequence. When comparing across a generation, fewer than 100 new genetic variants are typically detected in the euchromatic portion of a child’s genome. Driven by increasingly higher-resolution and higher-throughput sequencing technologies, newer and more accurate databases of genetic variation (for instance, more comprehensive structural variation data and phasing of combinations of variants along chromosomes) of worldwide populations will emerge to underpin the next era of discovery in human molecular genetics.

Accelerated deciphering of the genetic architecture of agricultural economic traits in pigs using a low-coverage whole-genome sequencing strategy

BACKGROUND

Uncovering the genetic architecture of economic traits in pigs is important for agricultural breeding. However, high-density haplotype reference panels are unavailable in most agricultural species, limiting accurate genotype imputation in large populations. Moreover, the infinitesimal model of quantitative traits implies that weak association signals tend to be spread across most of the genome, further complicating the genetic analysis. Hence, there is a need to develop new methods for sequencing large cohorts without large reference panels.

RESULTS

We describe a Tn5-based highly accurate, cost- and time-efficient, low-coverage sequencing method to obtain 11.3 million whole-genome single-nucleotide polymorphisms in 2,869 Duroc boars at a mean depth of 0.73×. On the basis of these single-nucleotide polymorphisms, a genome-wide association study was performed, resulting in 14 quantitative trait loci (QTLs) for 7 of 21 important agricultural traits in pigs. These QTLs harbour genes, such as ABCD4 for total teat number and HMGA1 for back fat thickness, and provided a starting point for further investigation. The inheritance models of the different traits varied greatly. Most follow the minor-polygene model, but this can be attributed to different reasons, such as the shaping of genetic architecture by artificial selection for this population and sufficiently interconnected minor gene regulatory networks.

CONCLUSIONS

Genome-wide association study results for 21 important agricultural traits identified 14 QTLs/genes and showed their genetic architectures, providing guidance for genetic improvement harnessing genomic features. The Tn5-based low-coverage sequencing method can be applied to large-scale genome studies for any species without a good reference panel and can be used for agricultural breeding.

Haplotype-based membership inference from summary genomic data

Motivation

The availability of human genomic data, together with the enhanced capacity to process them, is leading to transformative technological advances in biomedical science and engineering. However, the public dissemination of such data has been difficult due to privacy concerns. Specifically, it has been shown that the presence of a human subject in a case group can be inferred from the shared summary statistics of the group, e.g. the allele frequencies, or even the presence/absence of genetic variants (e.g. shared by the Beacon project) in the group. These methods rely on the availability of the target’s genome, i.e. the DNA profile of a target human subject, and thus are often referred to as the membership inference method.

Results

In this article, we demonstrate the haplotypes, i.e. the sequence of single nucleotide variations (SNVs) showing strong genetic linkages in human genome databases, may be inferred from the summary of genomic data without using a target’s genome. Furthermore, novel haplotypes that did not appear in the database may be reconstructed solely from the allele frequencies from genomic datasets. These reconstructed haplotypes can be used for a haplotype-based membership inference algorithm to identify target subjects in a case group with greater power than existing methods based on SNVs.

Availability and implementation

The implementation of the membership inference algorithms is available at https://github.com/diybu/Haplotype-based-membership-inferences.

Quantifying and Understanding the Higher Risk of Atherosclerotic Cardiovascular Disease Among South Asian Individuals: Results From the UK Biobank Prospective Cohort Study

BACKGROUND

Individuals of South Asian ancestry represent 23% of the global population, corresponding to 1.8 billion people, and have substantially higher risk of atherosclerotic cardiovascular disease compared with most other ethnicities. US practice guidelines now recognize South Asian ancestry as an important risk-enhancing factor. The magnitude of enhanced risk within the context of contemporary clinical care, the extent to which it is captured by existing risk estimators, and its potential mechanisms warrant additional study.

METHODS

Within the UK Biobank prospective cohort study, 8124 middle-aged participants of South Asian ancestry and 449 349 participants of European ancestry who were free of atherosclerotic cardiovascular disease at the time of enrollment were examined. The relationship of ancestry to risk of incident atherosclerotic cardiovascular disease-defined as myocardial infarction, coronary revascularization, or ischemic stroke-was assessed with Cox proportional hazards regression, along with examination of a broad range of clinical, anthropometric, and lifestyle mediators.

RESULTS

The mean age at study enrollment was 57 years, and 202 405 (44%) were male. Over a median follow-up of 11 years, 554 of 8124 (6.8%) individuals of South Asian ancestry experienced an atherosclerotic cardiovascular disease event compared with 19 756 of 449 349 (4.4%) individuals of European ancestry, corresponding to an adjusted hazard ratio of 2.03 (95% CI, 1.86-2.22; P<0.001). This higher relative risk was largely consistent across a range of age, sex, and clinical subgroups. Despite the >2-fold higher observed risk, the predicted 10-year risk of cardiovascular disease according to the American Heart Association/American College of Cardiology Pooled Cohort equations and QRISK3 equations was nearly identical for individuals of South Asian and European ancestry. Adjustment for a broad range of clinical, anthropometric, and lifestyle risk factors led to only modest attenuation of the observed hazard ratio to 1.45 (95% CI, 1.28-1.65, P<0.001). Assessment of variance explained by 18 candidate risk factors suggested greater importance of hypertension, diabetes, and central adiposity in South Asian individuals.

CONCLUSIONS

Within a large prospective study, South Asian individuals had substantially higher risk of atherosclerotic cardiovascular disease compared with individuals of European ancestry, and this risk was not captured by the Pooled Cohort Equations.

A recurrent, homozygous EMC10 frameshift variant is associated with a syndrome of developmental delay with variable seizures and dysmorphic features

PURPOSE

The endoplasmic reticulum membrane complex (EMC) is a highly conserved, multifunctional 10-protein complex related to membrane protein biology. In seven families, we identified 13 individuals with highly overlapping phenotypes who harbor a single identical homozygous frameshift variant in EMC10.

METHODS

Using exome, genome, and Sanger sequencing, a recurrent frameshift EMC10 variant was identified in affected individuals in an international cohort of consanguineous families. Multiple families were independently identified and connected via Matchmaker Exchange and internal databases. We assessed the effect of the frameshift variant on EMC10 RNA and protein expression and evaluated EMC10 expression in normal human brain tissue using immunohistochemistry.

RESULTS

A homozygous variant EMC10 c.287delG (Refseq NM_206538.3, p.Gly96Alafs*9) segregated with affected individuals in each family, who exhibited a phenotypic spectrum of intellectual disability (ID) and global developmental delay (GDD), variable seizures and variable dysmorphic features (elongated face, curly hair, cubitus valgus, and arachnodactyly). The variant arose on two founder haplotypes and results in significantly reduced EMC10 RNA expression and an unstable truncated EMC10 protein.

CONCLUSION

We propose that a homozygous loss-of-function variant in EMC10 causes a novel syndromic neurodevelopmental phenotype. Remarkably, the recurrent variant is likely the result of a hypermutable site and arose on distinct founder haplotypes.

Cultural Evolution of Genetic Heritability

Behavioral genetics and cultural evolution have both revolutionized our understanding of human behavior-largely independent of each other. Here we reconcile these two fields under a dual inheritance framework, offering a more nuanced understanding of the interaction between genes and culture. Going beyond typical analyses of gene-environment interactions, we describe the cultural dynamics that shape these interactions by shaping the environment and population structure. A cultural evolutionary approach can explain, for example, how factors such as rates of innovation and diffusion, density of cultural sub-groups, and tolerance for behavioral diversity impact heritability estimates, thus yielding predictions for different social contexts. Moreover, when cumulative culture functionally overlaps with genes, genetic effects become masked, unmasked, or even reversed, and the causal effects of an identified gene become confounded with features of the cultural environment. The manner of confounding is specific to a particular society at a particular time, but a WEIRD (Western, educated, industrialized, rich, democratic) sampling problem obscures this boundedness. Cultural evolutionary dynamics are typically missing from models of gene-to-phenotype causality, hindering generalizability of genetic effects across societies and across time. We lay out a reconciled framework and use it to predict the ways in which heritability should differ between societies, between socioeconomic levels and other groupings within some societies but not others, and over the life course. An integrated cultural evolutionary behavioral genetic approach cuts through the nature-nurture debate and helps resolve controversies in topics such as IQ.

Refinement of the clinical variant interpretation framework by statistical evidence and machine learning

BACKGROUND

Although the American College of Medical Genetics and Genomics/Association for Molecular Pathology (ACMG/AMP) guidelines for variant interpretation are used widely in clinical genetics, there is room for improvement of these knowledge-based guidelines.

METHODS

Statistical assessment of average deleteriousness of start-lost, stop-lost, and in-frame insertion and deletion (indel) variants and extraction of deleterious subsets was performed, being informed by proportions of rare variants in the general population of the Genome Aggregation Database (gnomAD). A machine learning-based model scoring the pathogenicity of start-lost variants (the PoStaL model) was constructed by predicting possible translation initiation sites on transcripts by deep learning and training a random forest on known pathogenic and likely benign variants.

FINDINGS

The proportion of rare variants was highest in stop-lost variants, followed by in-frame indels and start-lost variants, suggesting that the criteria in the ACMG/AMP guidelines assigning PVS (pathogenic very strong) to start-lost variants and PM (pathogenic moderate) to stop-lost and in-frame indel variants would not be appropriate. Regarding deleterious subsets, stop-lost variants introducing extensions of more than 30 amino acids and in-frame indels computationally predicted to be damaging are enriched for rare and known pathogenic variants. For start-lost variants, we developed the PoStaL model, which outperforms existing tools. We also provide comprehensive lists of the PoStaL scores for start-lost variants and the length of extended amino acids by stop-lost variants.

CONCLUSIONS

Our study could contribute to refinement of the ACMG/AMP guidelines, provides resources for future investigation, and provides an example of how to improve knowledge-based frameworks by data-driven approaches.

FUNDING

The study was supported by grants from the Japan Agency for Medical Research and Development (AMED) and the Japan Society for the Promotion of Science (JSPS).

Genetic risk factors for colorectal cancer in multiethnic Indonesians

Colorectal cancer is a common cancer in Indonesia, yet it has been understudied in this resource-constrained setting. We conducted a genome-wide association study focused on evaluation and preliminary discovery of colorectal cancer risk factors in Indonesians. We administered detailed questionnaires and collecting blood samples from 162 colorectal cancer cases throughout Makassar, Indonesia. We also established a control set of 193 healthy individuals frequency matched by age, sex, and ethnicity. A genome-wide association analysis was performed on 84 cases and 89 controls passing quality control. We evaluated known colorectal cancer genetic variants using logistic regression and established a genome-wide polygenic risk model using a Bayesian variable selection technique. We replicate associations for rs9497673, rs6936461 and rs7758229 on chromosome 6; rs11255841 on chromosome 10; and rs4779584, rs11632715, and rs73376930 on chromosome 15. Polygenic modeling identified 10 SNP associated with colorectal cancer risk. This work helps characterize the relationship between variants in the SCL22A3, SCG5, GREM1, and STXBP5-AS1 genes and colorectal cancer in a diverse Indonesian population. With further biobanking and international research collaborations, variants specific to colorectal cancer risk in Indonesians will be identified.

Advances in Genomic Discovery and Implications for Personalized Prevention and Medicine: Estonia as Example

The current paradigm of personalized medicine envisages the use of genomic data to provide predictive information on the health course of an individual with the aim of prevention and individualized care. However, substantial efforts are required to realize the concept: enhanced genetic discoveries, translation into intervention strategies, and a systematic implementation in healthcare. Here we review how further genetic discoveries are improving personalized prediction and advance functional insights into the link between genetics and disease. In the second part we give our perspective on the way these advances in genomic research will transform the future of personalized prevention and medicine using Estonia as a primer.

How understudied populations have contributed to our understanding of Alzheimer's disease genetics

The majority of genome-wide association studies have been conducted using samples with a broadly European genetic background. As a field, we acknowledge this limitation and the need to increase the diversity of populations studied. A major challenge when designing and conducting such studies is to assimilate large samples sizes so that we attain enough statistical power to detect variants associated with disease, particularly when trying to identify variants with low and rare minor allele frequencies. In this review, we aimed to illustrate the benefits to genetic characterization of Alzheimer's disease, in researching currently understudied populations. This is important for both fair representation of world populations and the translatability of findings. To that end, we conducted a literature search to understand the contributions of studies, on different populations, to Alzheimer's disease genetics. Using both PubMed and Alzforum Mutation Database, we systematically quantified the number of studies reporting variants in known disease-causing genes, in a worldwide manner, and discuss the contributions of research in understudied populations to the identification of novel genetic factors in this disease. Additionally, we compared the effects of genome-wide significant single nucleotide polymorphisms across populations by focusing on loci that show different association profiles between populations (a key example being APOE). Reports of variants in APP, PSEN1 and PSEN2 can initially determine whether patients from a country have been studied for Alzheimer's disease genetics. Most genome-wide significant associations in non-Hispanic white genome-wide association studies do not reach genome-wide significance in such studies of other populations, with some suggesting an opposite effect direction; this is likely due to much smaller sample sizes attained. There are, however, genome-wide significant associations first identified in understudied populations which have yet to be replicated. Familial studies in understudied populations have identified rare, high effect variants, which have been replicated in other populations. This work functions to both highlight how understudied populations have furthered our understanding of Alzheimer's disease genetics, and to help us gauge our progress in understanding the genetic architecture of this disease in all populations.

Genetic Modifiers of Age at Onset for Parkinson's Disease in Asians: A Genome-Wide Association Study

BACKGROUND

Age at onset (AAO) is an essential feature of Parkinson's disease (PD) and can help predict disease progression and mortality. Identification of genetic variants influencing AAO of PD could lead to a better understanding of the disease's biological mechanism and provide clinical guidance. However, genetic determinants for AAO of PD remain mostly unknown, especially in the Asian population.

OBJECTIVES

To identify genetic determinants for AAO of PD in the Asian population.

METHODS

We performed a genome-wide association meta-analysis on AAO of PD in 5166 Chinese patients with PD (N_discovery = 3628, N_replication = 1538). We then conducted a further cross-ethnic meta-analysis using our results and summary statistics for the AAO of PD from the European population.

RESULTS

The total heritability of AAO of PD was around 0.10 ~ 0.14, similar to that (~0.11) estimated in populations of European ancestry. One novel significant intergenic locus rs9783733 (NDN; PWRN4) was identified (P = 3.14E-09, beta = 2.30, SE = 0.39). Remarkably, this variant could delay AAO of PD by ~2.43 years, with a more considerable effect on males (~3.18 years) than females (~1.45 years). The variant was suggestively significant in the cross-ethnic meta-analysis and suggested a positive selection in the East Asian population. Additionally, cross-ethnic meta-analysis identified a significant locus rs356203 in SNCA (P = 2.35E-11, beta = -0.71, SE = 0.01).

CONCLUSIONS

These findings improve the current understanding of the genetic etiology of AAO of PD in different ethnic groups, and provide a new target for further research on PD pathogenesis and potential therapeutic options. © 2021 International Parkinson and Movement Disorder Society.

Validation of whole genome sequencing from dried blood spots

BACKGROUND

Dried blood spots (DBS) are a relatively inexpensive source of nucleic acids and are easy to collect, transport, and store in large-scale field surveys, especially in resource-limited settings. However, their performance in whole-genome sequencing (WGS) relative to that of venous blood DNA has not been analyzed for various downstream applications.

METHODS

This study compares the WGS performance of DBS paired with venous blood samples collected from 12 subjects.

RESULTS

Results of standard quality checks of coverage, base quality, and mapping quality were found to be near identical between DBS and venous blood. Concordance for single-nucleotide variants, insertions and deletions, and copy number variants was high between these two sample types. Additionally, downstream analyses typical of population-based studies were performed, such as mitochondrial heteroplasmy detection, haplotype analysis, mitochondrial copy number changes, and determination of telomere lengths. The absolute mitochondrial copy number values were higher for DBS than for venous blood, though the trend in sample-to-sample variation was similar between DBS and blood. Telomere length estimates in most DBS samples were on par with those from venous blood.

CONCLUSION

DBS samples can serve as a robust and feasible alternative to venous blood for studies requiring WGS analysis.

Three Signatures of Adaptive Polymorphism Exemplified by Malaria-Associated Genes

Malaria has been one of the strongest selective pressures on our species. Many of the best-characterized cases of adaptive evolution in humans are in genes tied to malaria resistance. However, the complex evolutionary patterns at these genes are poorly captured by standard scans for nonneutral evolution. Here, we present three new statistical tests for selection based on population genetic patterns that are observed more than once among key malaria resistance loci. We assess these tests using forward-time evolutionary simulations and apply them to global whole-genome sequencing data from humans, and thus we show that they are effective at distinguishing selection from neutrality. Each test captures a distinct evolutionary pattern, here called Divergent Haplotypes, Repeated Shifts, and Arrested Sweeps, associated with a particular period of human prehistory. We clarify the selective signatures at known malaria-relevant genes and identify additional genes showing similar adaptive evolutionary patterns. Among our top outliers, we see a particular enrichment for genes involved in erythropoiesis and for genes previously associated with malaria resistance, consistent with a major role for malaria in shaping these patterns of genetic diversity. Polymorphisms at these genes are likely to impact resistance to malaria infection and contribute to ongoing host-parasite coevolutionary dynamics.

Identification of Novel Candidate Genes and Variants for Hearing Loss and Temporal Bone Anomalies

Background: Hearing loss remains an important global health problem that is potentially addressed through early identification of a genetic etiology, which helps to predict outcomes of hearing rehabilitation such as cochlear implantation and also to mitigate the long-term effects of comorbidities. The identification of variants for hearing loss and detailed descriptions of clinical phenotypes in patients from various populations are needed to improve the utility of clinical genetic screening for hearing loss. Methods: Clinical and exome data from 15 children with hearing loss were reviewed. Standard tools for annotating variants were used and rare, putatively deleterious variants were selected from the exome data. Results: In 15 children, 21 rare damaging variants in 17 genes were identified, including: 14 known hearing loss or neurodevelopmental genes, 11 of which had novel variants; and three candidate genes IST1, CBLN3 and GDPD5, two of which were identified in children with both hearing loss and enlarged vestibular aqueducts. Patients with variants within IST1 and MYO18B had poorer outcomes after cochlear implantation. Conclusion: Our findings highlight the importance of identifying novel variants and genes in ethnic groups that are understudied for hearing loss.

Human ACE2 receptor polymorphisms and altered susceptibility to SARS-CoV-2

COVID-19 is a respiratory illness caused by a novel coronavirus called SARS-CoV-2. The viral spike (S) protein engages the human angiotensin-converting enzyme 2 (ACE2) receptor to invade host cells with ~10-15-fold higher affinity compared to SARS-CoV S-protein, making it highly infectious. Here, we assessed if ACE2 polymorphisms can alter host susceptibility to SARS-CoV-2 by affecting this interaction. We analyzed over 290,000 samples representing >400 population groups from public genomic datasets and identified multiple ACE2 protein-altering variants. Using reported structural data, we identified natural ACE2 variants that could potentially affect virus-host interaction and thereby alter host susceptibility. These include variants S19P, I21V, E23K, K26R, T27A, N64K, T92I, Q102P and H378R that were predicted to increase susceptibility, while variants K31R, N33I, H34R, E35K, E37K, D38V, Y50F, N51S, M62V, K68E, F72V, Y83H, G326E, G352V, D355N, Q388L and D509Y were predicted to be protective variants that show decreased binding to S-protein. Using biochemical assays, we confirmed that K31R and E37K had decreased affinity, and K26R and T92I variants showed increased affinity for S-protein when compared to wildtype ACE2. Consistent with this, soluble ACE2 K26R and T92I were more effective in blocking entry of S-protein pseudotyped virus suggesting that ACE2 variants can modulate susceptibility to SARS-CoV-2.

A novel homozygous missense mutation p.P388S in TULP1 causes protein instability and retinitis pigmentosa

Purpose

Retinitis pigmentosa (RP) is an inherited retinal disorder that results in the degeneration of photoreceptor cells, ultimately leading to severe visual impairment. We characterized a consanguineous family from Southern India wherein a 25 year old individual presented with night blindness since childhood. The purpose of this study was to identify the causative mutation for RP in this individual as well as characterize how the mutation may ultimately affect protein function.

Methods

We performed a complete ophthalmologic examination of the proband followed by exome sequencing. The likely causative mutation was identified and modeled in cultured cells, evaluating its expression, solubility (both with western blotting), subcellular distribution, (confocal microscopy), and testing whether this variant induced endoplasmic reticulum (ER) stress (quantitative PCR [qPCR] and western blotting).

Results

The proband presented with generalized and parafoveal retinal pigmented epithelium (RPE) atrophy with bone spicule-like pigmentation in the midperiphery and arteriolar attenuation. Optical coherence tomography scans through the macula of both eyes showed atrophy of the outer retinal layers with loss of the ellipsoid zone, whereas the systemic examination of this individual was normal. The proband's parents and sibling were asymptomatic and had normal funduscopic examinations. We discovered a novel homozygous p.Pro388Ser mutation in the tubby-like protein 1 (TULP1) gene in the individual with RP. In cultured cells, the P388S mutation does not alter the subcellular distribution of TULP1 or induce ER stress when compared to wild-type TULP1, but instead significantly lowers protein stability as indicated with steady-state and cycloheximide-chase experiments.

Conclusions

These results add to the list of known mutations in TULP1 identified in individuals with RP and suggest a possible unique pathogenic mechanism in TULP1-induced RP, which may be shared among select mutations in TULP1.

Genomic insights into population history and biological adaptation in Oceania

The Pacific region is of major importance for addressing questions regarding human dispersals, interactions with archaic hominins and natural selection processes1. However, the demographic and adaptive history of Oceanian populations remains largely uncharacterized. Here we report high-coverage genomes of 317 individuals from 20 populations from the Pacific region. We find that the ancestors of Papuan-related ('Near Oceanian') groups underwent a strong bottleneck before the settlement of the region, and separated around 20,000-40,000 years ago. We infer that the East Asian ancestors of Pacific populations may have diverged from Taiwanese Indigenous peoples before the Neolithic expansion, which is thought to have started from Taiwan around 5,000 years ago2-4. Additionally, this dispersal was not followed by an immediate, single admixture event with Near Oceanian populations, but involved recurrent episodes of genetic interactions. Our analyses reveal marked differences in the proportion and nature of Denisovan heritage among Pacific groups, suggesting that independent interbreeding with highly structured archaic populations occurred. Furthermore, whereas introgression of Neanderthal genetic information facilitated the adaptation of modern humans related to multiple phenotypes (for example, metabolism, pigmentation and neuronal development), Denisovan introgression was primarily beneficial for immune-related functions. Finally, we report evidence of selective sweeps and polygenic adaptation associated with pathogen exposure and lipid metabolism in the Pacific region, increasing our understanding of the mechanisms of biological adaptation to island environments.

Multiple migrations to the Philippines during the last 50,000 years

A key link to understand human history in Island Southeast Asia is the Philippine archipelago and its poorly investigated genetic diversity. We analyzed the most comprehensive set of population-genomic data for the Philippines: 1,028 individuals covering 115 indigenous communities. We demonstrate that the Philippines were populated by at least five waves of human migration. The Cordillerans migrated into the Philippines prior to the arrival of rice agriculture, where some remain as the least admixed East Asians carrying an ancestry shared by all Austronesian-speaking populations, thereby challenging an exclusive out-of-Taiwan model of joint farming–language–people dispersal. Altogether, our findings portray the Philippines as a crucial gateway, with a multilayered history, that ultimately changed the genetic landscape of the Asia-Pacific region.

Island Southeast Asia has recently produced several surprises regarding human history, but the region’s complex demography remains poorly understood. Here, we report ∼2.3 million genotypes from 1,028 individuals representing 115 indigenous Philippine populations and genome-sequence data from two ∼8,000-y-old individuals from Liangdao in the Taiwan Strait. We show that the Philippine islands were populated by at least five waves of human migration: initially by Northern and Southern Negritos (distantly related to Australian and Papuan groups), followed by Manobo, Sama, Papuan, and Cordilleran-related populations. The ancestors of Cordillerans diverged from indigenous peoples of Taiwan at least ∼8,000 y ago, prior to the arrival of paddy field rice agriculture in the Philippines ∼2,500 y ago, where some of their descendants remain to be the least admixed East Asian groups carrying an ancestry shared by all Austronesian-speaking populations. These observations contradict an exclusive “out-of-Taiwan” model of farming–language–people dispersal within the last four millennia for the Philippines and Island Southeast Asia. Sama-related ethnic groups of southwestern Philippines additionally experienced some minimal South Asian gene flow starting ∼1,000 y ago. Lastly, only a few lowlanders, accounting for <1% of all individuals, presented a low level of West Eurasian admixture, indicating a limited genetic legacy of Spanish colonization in the Philippines. Altogether, our findings reveal a multilayered history of the Philippines, which served as a crucial gateway for the movement of people that ultimately changed the genetic landscape of the Asia-Pacific region.

A data set of variants derived from 1455 clinical and research exomes is efficient in variant prioritization for early-onset monogenic disorders in Indians

Given the genomic uniqueness, a local data set is most desired for Indians, who are underrepresented in existing public databases. We hypothesize patients with rare monogenic disorders and their family members can provide a reliable source of common variants in the population. Exome sequencing (ES) data from families with rare Mendelian disorders was aggregated from five centers in India. The dataset was refined by excluding related individuals and removing the disease-causing variants (refined cohort). The efficiency of these data sets was assessed in a new set of 50 exomes against gnomAD and GenomeAsia. Our original cohort comprised 1455 individuals from 1203 families. The refined cohort had 836 unrelated individuals that retained 1,251,064 variants with 181,125 population-specific and 489,618 common variants. The allele frequencies from our cohort helped to define 97,609 rare variants in gnomAD and 44,520 rare variants in GenomeAsia as common variants in our population. Our variant dataset provided an additional 1.7% and 0.1% efficiency for prioritizing heterozygous and homozygous variants respectively for rare monogenic disorders. We observed additional 19 genes/human knockouts. We list carrier frequency for 142 recessive disorders. This is a large and useful resource of exonic variants for Indians. Despite limitations, datasets from patients are efficient tools for variant prioritization in a resource-limited setting.

Scanning the human genome for "signatures" of positive selection: Transformative opportunities and ethical obligations

The relationship history of evolutionary anthropology and genetics is complex. At best, genetics is a beautifully integrative part of the discipline. Yet this integration has also been fraught, with punctuated, disruptive challenges to dogma, periodic reluctance by some members of the field to embrace results from analyses of genetic data, and occasional over-assertions of genetic definitiveness by geneticists. At worst, evolutionary genetics has been a tool for reinforcing racism and colonialism. While a number of genetics/genomics papers have disproportionately impacted evolutionary anthropology, here we highlight the 2002 presentation of an elegantly powerful approach for identifying "signatures" of past positive selection from haplotype-based patterns of genetic variation. Together with technological advances in genotyping methods, this article transformed our field by facilitating genome-wide "scans" for signatures of past positive selection in human populations. This approach helped researchers test longstanding evolutionary anthropology hypotheses while simultaneously providing opportunities to develop entirely new ones. Genome-wide scans for signatures of positive selection have since been conducted in diverse worldwide populations, with striking findings of local adaptation and convergent evolution. Yet there are ethical considerations with respect to the ubiquity of these studies and the cross-application of the genome-wide scan approach to existing datasets, which we also discuss.

Founder variants and population genomes-Toward precision medicine

Human migration and community specific cultural practices have contributed to founder events and enrichment of the variants associated with genetic diseases. While many founder events in isolated populations have remained uncharacterized, the application of genomics in clinical settings as well as for population scale studies in the recent years have provided an unprecedented push towards identification of founder variants associated with human health and disease. The discovery and characterization of founder variants could have far reaching implications not only in understanding the history or genealogy of the disease, but also in implementing evidence based policies and genetic testing frameworks. This further enables precise diagnosis and prevention in an attempt towards precision medicine. This review provides an overview of founder variants along with methods and resources cataloging them. We have also discussed the public health implications and examples of prevalent disease associated founder variants in specific populations.

Genetic profiles of 103,106 individuals in the Taiwan Biobank provide insights into the health and history of Han Chinese

Personalized medical care focuses on prediction of disease risk and response to medications. To build the risk models, access to both large-scale genomic resources and human genetic studies is required. The Taiwan Biobank (TWB) has generated high-coverage, whole-genome sequencing data from 1492 individuals and genome-wide SNP data from 103,106 individuals of Han Chinese ancestry using custom SNP arrays. Principal components analysis of the genotyping data showed that the full range of Han Chinese genetic variation was found in the cohort. The arrays also include thousands of known functional variants, allowing for simultaneous ascertainment of Mendelian disease-causing mutations and variants that affect drug metabolism. We found that 21.2% of the population are mutation carriers of autosomal recessive diseases, 3.1% have mutations in cancer-predisposing genes, and 87.3% carry variants that affect drug response. We highlight how TWB data provide insight into both population history and disease burden, while showing how widespread genetic testing can be used to improve clinical care.

The impact of global and local Polynesian genetic ancestry on complex traits in Native Hawaiians

Epidemiological studies of obesity, Type-2 diabetes (T2D), cardiovascular diseases and several common cancers have revealed an increased risk in Native Hawaiians compared to European- or Asian-Americans living in the Hawaiian islands. However, there remains a gap in our understanding of the genetic factors that affect the health of Native Hawaiians. To fill this gap, we studied the genetic risk factors at both the chromosomal and sub-chromosomal scales using genome-wide SNP array data on ~4,000 Native Hawaiians from the Multiethnic Cohort. We estimated the genomic proportion of Native Hawaiian ancestry ("global ancestry," which we presumed to be Polynesian in origin), as well as this ancestral component along each chromosome ("local ancestry") and tested their respective association with binary and quantitative cardiometabolic traits. After attempting to adjust for non-genetic covariates evaluated through questionnaires, we found that per 10% increase in global Polynesian genetic ancestry, there is a respective 8.6%, and 11.0% increase in the odds of being diabetic (P = 1.65×10-4) and having heart failure (P = 2.18×10-4), as well as a 0.059 s.d. increase in BMI (P = 1.04×10-10). When testing the association of local Polynesian ancestry with risk of disease or biomarkers, we identified a chr6 region associated with T2D. This association was driven by an uniquely prevalent variant in Polynesian ancestry individuals. However, we could not replicate this finding in an independent Polynesian cohort from Samoa due to the small sample size of the replication cohort. In conclusion, we showed that Polynesian ancestry, which likely capture both genetic and lifestyle risk factors, is associated with an increased risk of obesity, Type-2 diabetes, and heart failure, and that larger cohorts of Polynesian ancestry individuals will be needed to replicate the putative association on chr6 with T2D.

Origins of modern human ancestry

New finds in the palaeoanthropological and genomic records have changed our view of the origins of modern human ancestry. Here we review our current understanding of how the ancestry of modern humans around the globe can be traced into the deep past, and which ancestors it passes through during our journey back in time. We identify three key phases that are surrounded by major questions, and which will be at the frontiers of future research. The most recent phase comprises the worldwide expansion of modern humans between 40 and 60 thousand years ago (ka) and their last known contacts with archaic groups such as Neanderthals and Denisovans. The second phase is associated with a broadly construed African origin of modern human diversity between 60 and 300 ka. The oldest phase comprises the complex separation of modern human ancestors from archaic human groups from 0.3 to 1 million years ago. We argue that no specific point in time can currently be identified at which modern human ancestry was confined to a limited birthplace, and that patterns of the first appearance of anatomical or behavioural traits that are used to define Homo sapiens are consistent with a range of evolutionary histories.

CYP2C9 Variations and Their Pharmacogenetic Implications Among Diverse South Asian Populations

Introduction

Allelic frequency distribution of drug metabolizing enzyme genes among populations is important to identify risk groups for adverse drug reaction and to select representative populations for clinical trials. Although India emerged as an important hub for clinical trials, information about the pharmacogenetic diversity for this region is still lacking. Here, we investigated genetic diversity of cytochrome-P450-2C9 (CYP2C9) gene which metabolizes wide range of drugs and is highly expressed in the human liver.

Methods

In total, 1278 individuals from 36 diverse Indian populations, 210 individuals from in-house data-repository and 489 other South Asian samples from the 1000 Genomes Project were selected. Variants observed in CYP2C9 gene were subjected to various statistical analyses.

Results

High frequency of CYP2C9*3 (~13%) and CYP2C9*3/*3 (~1%) was observed among South Asians, compared to 21 populations living outside the Indian subcontinent. The allelic/genotypic frequency does not correlate with geographical location or linguistic affiliation, except populations speaking Tibeto-Burmans language, who have lower frequency of CYP2C9*3 and CYP2C9*3/*3. Since, South Asians practice strict endogamy, presence of unique mutation and high frequency of homozygous genotypes not surprising. CYP2C9*3 has been associated with therapeutic response.The effect of CYP2C9*3/*3 is more pronounced compared to heterozygous and wild type homozygous genotypes as evident in many in vitro studies. As South Asians have high frequency, it would be interesting to explore potential of CYP2C9*3 as a marker for personalized therapy. Our study revealed several rare functional variants, which form eight novel and rare haplotypes of CYP2C9 (CYP2C9*63–*70). Of which, CYP2C9*64, *65, *66, *68, *69 and *70 haplotypes are South Asian-specific.

Conclusion

Overall, we find high genetic heterogeneity within South Asians and identified South Asian-specific putative functional CYP2C9 haplotypes. High frequency of CYP2C9*3 and CYP2C9*3/*3 was observed in South Asian populations. Taken together, current study greatly enriches the knowledge of naturally occurring CYP2C9 variants and its diversity in South Asia, which are relevant to further CYP2C9-related functional research and for personalized medicine.

V, Jain S, Dash D, Kumar NS, Vanlallawma A, Sarma R, Chhakchhuak L, Kalyanaraman S, Mahadevan R, Kandasamy S, B. M. P, Rajagopal RE, J. ER, P. ND, Bajaj A, Gupta V, Mathew S, Goswami S, Mangla M, Prakash S, Joshi K, S. S, Gajjar D, Soraisham R, Yadav R, Devi YS, Gupta A, Mukerji M, Ramalingam S, B. K. B, Scaria V, Sivasubbu S. IndiGenomes: a comprehensive resource of genetic variants from over 1000 Indian genomes

With the advent of next-generation sequencing, large-scale initiatives for mining whole genomes and exomes have been employed to better understand global or population-level genetic architecture. India encompasses more than 17% of the world population with extensive genetic diversity, but is under-represented in the global sequencing datasets. This gave us the impetus to perform and analyze the whole genome sequencing of 1029 healthy Indian individuals under the pilot phase of the 'IndiGen' program. We generated a compendium of 55,898,122 single allelic genetic variants from geographically distinct Indian genomes and calculated the allele frequency, allele count, allele number, along with the number of heterozygous or homozygous individuals. In the present study, these variants were systematically annotated using publicly available population databases and can be accessed through a browsable online database named as 'IndiGenomes' http://clingen.igib.res.in/indigen/. The IndiGenomes database will help clinicians and researchers in exploring the genetic component underlying medical conditions. Till date, this is the most comprehensive genetic variant resource for the Indian population and is made freely available for academic utility. The resource has also been accessed extensively by the worldwide community since it's launch.

Genome of a middle Holocene hunter-gatherer from Wallacea

Much remains unknown about the population history of early modern humans in southeast Asia, where the archaeological record is sparse and the tropical climate is inimical to the preservation of ancient human DNA¹. So far, only two low-coverage pre-Neolithic human genomes have been sequenced from this region. Both are from mainland Hòabìnhian hunter-gatherer sites: Pha Faen in Laos, dated to 7939-7751 calibrated years before present (yr cal BP; present taken as AD 1950), and Gua Cha in Malaysia (4.4-4.2 kyr cal BP)¹. Here we report, to our knowledge, the first ancient human genome from Wallacea, the oceanic island zone between the Sunda Shelf (comprising mainland southeast Asia and the continental islands of western Indonesia) and Pleistocene Sahul (Australia-New Guinea). We extracted DNA from the petrous bone of a young female hunter-gatherer buried 7.3-7.2 kyr cal BP at the limestone cave of Leang Panninge² in South Sulawesi, Indonesia. Genetic analyses show that this pre-Neolithic forager, who is associated with the 'Toalean' technocomplex^3,4, shares most genetic drift and morphological similarities with present-day Papuan and Indigenous Australian groups, yet represents a previously unknown divergent human lineage that branched off around the time of the split between these populations approximately 37,000 years ago⁵. We also describe Denisovan and deep Asian-related ancestries in the Leang Panninge genome, and infer their large-scale displacement from the region today.

Practical guide for managing large-scale human genome data in research

Studies in human genetics deal with a plethora of human genome sequencing data that are generated from specimens as well as available on public domains. With the development of various bioinformatics applications, maintaining the productivity of research, managing human genome data, and analyzing downstream data is essential. This review aims to guide struggling researchers to process and analyze these large-scale genomic data to extract relevant information for improved downstream analyses. Here, we discuss worldwide human genome projects that could be integrated into any data for improved analysis. Obtaining human whole-genome sequencing data from both data stores and processes is costly; therefore, we focus on the development of data format and software that manipulate whole-genome sequencing. Once the sequencing is complete and its format and data processing tools are selected, a computational platform is required. For the platform, we describe a multi-cloud strategy that balances between cost, performance, and customizability. A good quality published research relies on data reproducibility to ensure quality results, reusability for applications to other datasets, as well as scalability for the future increase of datasets. To solve these, we describe several key technologies developed in computer science, including workflow engine. We also discuss the ethical guidelines inevitable for human genomic data analysis that differ from model organisms. Finally, the future ideal perspective of data processing and analysis is summarized.

Exome sequencing identifies ARID2 as a novel tumor suppressor in early-onset sporadic rectal cancer

Early-onset sporadic rectal cancer (EOSRC) is a unique and predominant colorectal cancer (CRC) subtype in India. In order to understand the tumorigenic process in EOSRC, we performed whole-exome sequencing of 47 microsatellite stable EOSRC samples. Signature 1 was the predominant mutational signature in EOSRC, as previously shown in other CRC exome studies. More importantly, we identified TP53, KRAS, APC, PIK3R1, SMAD4 and ZNF880 as significantly mutated (q < 0.1) and ARID1A and ARID2 as near-significantly mutated (restricted hypothesis testing; q < 0.1) candidate drivers. Unlike the other candidates, the tumorigenic potential of ARID2, encoding a component of the SWI/SNF chromatin remodeling complex, is largely unexplored in CRC. shRNA-mediated ARID2 knockdown performed in different CRC cell lines resulted in significant alterations in transcript levels of cancer-related target genes. More importantly, ARID2 knockdown promoted several tumorigenic features including cell viability, proliferation, ability to override contact inhibition of growth, and migration besides significantly increasing tumor formation ability in nude mice. The observed gain in tumorigenic features was rescued upon ectopic expression of wild type but not mutant ARID2. Analyses of the TCGA pan-cancer dataset revealed several modes of ARID2 inactivation and of the CRC dataset revealed poorer survival in patients with ARID2 alterations. We therefore propose ARID2 as a novel tumor suppressor in CRC.

Whole-exome analyses of congenital muscular dystrophy and congenital myopathy patients from India reveal a wide spectrum of known and novel mutations

BACKGROUND AND PURPOSE

Congenital muscular dystrophies (CMDs) and congenital myopathies (CMs) are a group of genetically and clinically heterogeneous degenerative primary muscle disorders with onset at birth or during infancy. Due to vast heterogeneity, clinical examination and protein-based analyses often fail to identify the genetic causes of these diseases. The aim of this study was to genetically diagnose a cohort of 36 difficult-to-diagnose CMD and CM cases of Indian origin using next-generation sequencing methods.

METHODS

Whole-exome sequencing (WES) was performed to identify pathogenic mutations in previously reported CMD and CM-related genes using variant calling and stringent variant filtration process. Subsequently, in silico homology modelling and molecular dynamics simulations (MDS) studies were undertaken for a number of novel and missense variants.

RESULTS

A total of 33 and 21 rare and deleterious mutations were identified in 28 genes previously reported in CMD and CM based on OMIM, ClinVar and Orphanet, respectively. We could accurately diagnose 54% patients (n = 12/22) in the CMD group and 35% patients (n = 5/14) in the CM group. Furthermore, MDS studies for mutations located in LMNA, LAMA2 and RYR1 suggest that the wild-type proteins are more stable than their mutant counterparts, implying a potential mechanism of pathogenesis.

CONCLUSION

The WES findings led us to identify reported as well as novel variants for the first time in Indian patients with CMD and CM. This allowed us to achieve an accurate genetic diagnosis, which was difficult using conventional diagnostic tools. Transferring these WES findings to clinical practice will help guide clinical care of the affected patients and inform genetic counselling.

Dysbiotic Lesional Microbiome With Filaggrin Missense Variants Associate With Atopic Dermatitis in India

Background: Atopic Dermatitis (AD) has been associated with the loss of function (LoF) mutations in Filaggrin (FLG) gene and increase in relative abundance of specific microbes in the lesional skin, predominantly in Caucasians. Our study aims to determine, in Indian AD patients, (a) the prevalence of FLG LoF and missense mutations, and (b) the nature and extent of dysbiosis and altered microbial pathways with and without mutations in FLG. AD patients (n = 34) and healthy controls (n = 54) were recruited from India in this study and shotgun sequencing was carried out in a subset of samples with adequate microbiome DNA concentration. Host DNA from the same subset of samples was subjected to FLG coding region sequencing and host-microbiome association was estimated. Results: The prevalence of FLG LoFs that are associated with AD globally were significantly lesser in our cases and controls (8.6%, 0%) than those reported in Europeans (27%, 2.6%). Staphylococcus aureus was present only on AD skin [abundance in Pediatric AD: 32.86%; Adult AD: 22.17%], but not on healthy skin on which Staphylococcus hominis (Adult controls: 16.43%, Adult AD: 0.20%; p = 0.002), Cutibacterium acnes (Adult controls:10.84%, Adult AD: 0.90%; p = 0.02), and Malassezia globosa (Adult controls: 8.89%, Adult AD: 0.005%; p = 0.001) were significantly more abundant. Microbial pathways mostly associated with skin barrier permeability, ammonia production and inflammation (Arginine and Proline metabolism, Histidine Metabolism and Staphylococcus aureus infection) were significantly enriched on AD skin metagenome. These pathways are also reported to impair antimicrobial peptide activity. Among AD patients with missense single nucleotide polymorphisms harboring "potentially damaging" alleles in FLG gene, damaging allele dosage was significantly (p < 0.02) positively correlated with relative abundance of phylum_Proteobacteria up to order_Pseudomonadales and negatively correlated with phylum_Firmicutes up to species_Staphylococcus aureus. Conclusion: Our study has provided evidence that host DNA profile is significantly associated with microbiome composition in the development of AD. Species and strain level analysis showed that the microbial pathways enriched in AD cases were mostly found in MRSA strains. These evidences can be harnessed to control AD by modulating the microbiome using a personalized strategy. Our findings on the association of FLG genotypes with the microbiome dysbiosis may pave the way for a personalized strategy to provide a more effective control of AD.

Global Picture of Genetic Relatedness and the Evolution of Humankind

Simple Summary

The intricacies of human ancestry are buried deep within our DNA. For years, scientists have been working to piece together a vast picture of our genetic lineage. The purpose of this study was to further reveal this global picture of human genetic relatedness using identical-by-descent (IBD) genomic fragments. We processed over 65 million very rare single nucleotide polymorphic (SNP) alleles and detected over 17 million shared IBD fragments, including very short IBD fragments that allowed us to trace common ancestors back to 200,000 years ago. We also determined nine geographical regions representing nine unique genetic components for mankind: East and West Africa, Northern Europe, Arctica, East Asia, Oceania, South Asia, Middle East, and South America. The levels of admixture in every studied population could be assigned to one of these regions and long-term neighboring populations are strikingly similar, despite any political, religious, and cultural differences. Additionally, we observed the topmost admixture to be in central Eurasia. The entire picture of relatedness of all the studied populations presents itself in the form of shared number/size of IBDs, providing novel insights into geographical admixtures and genetic contributions that shaped human ancestry into what it is today.

Abstract

We performed an exhaustive pairwise comparison of whole-genome sequences of 3120 individuals, representing 232 populations from all continents and seven prehistoric people including archaic and modern humans. In order to reveal an intricate picture of worldwide human genetic relatedness, 65 million very rare single nucleotide polymorphic (SNP) alleles have been bioinformatically processed. The number and size of shared identical-by-descent (IBD) genomic fragments for every pair of 3127 individuals have been revealed. Over 17 million shared IBD fragments have been described. Our approach allowed detection of very short IBD fragments (<20 kb) that trace common ancestors who lived up to 200,000 years ago. We detected nine distinct geographical regions within which individuals had strong genetic relatedness, but with negligible relatedness between the populations of these regions. The regions, comprising nine unique genetic components for mankind, are the following: East and West Africa, Northern Europe, Arctica, East Asia, Oceania, South Asia, Middle East, and South America. The level of admixture in every studied population has been apportioned among these nine genetic components. Genetically, long-term neighboring populations are strikingly similar to each other in spite of any political, religious, and cultural differences. The topmost admixture has been observed at the center of Eurasia. These admixed populations (including Uyghurs, Azerbaijanis, Uzbeks, and Iranians) have roughly equal genetic contributions from the Middle East, Europe, China, and India, with additional significant traces from Africa and Arctic. The entire picture of relatedness of all the studied populations unfolds and presents itself in the form of shared number/size of IBDs.

Phenotype prediction and characterization of 25 pharmacogenes in Thais from whole genome sequencing for clinical implementation

Publicly available pharmacogenomics (PGx) databases enable translation of genotype data into clinically actionable information. As variation within pharmacogenes is population-specific, this study investigated the spectrum of 25 clinically relevant pharmacogenes in the Thai population (n = 291) from whole genome sequencing. The bioinformatics tool Stargazer was used for phenotype prediction, through assignment of alleles and detection of structural variation. Known and unreported potentially deleterious PGx variants were identified. Over 25% of Thais carried a high-risk diplotype in CYP3A5, CYP2C19, CYP2D6, NAT2, SLCO1B1, and UGT1A1. CYP2D6 structural variants accounted for 83.8% of all high-risk diplotypes. Of 39 known PGx variants identified, six variants associated with adverse drug reactions were common. Allele frequencies of CYP3A5*3 (rs776746), CYP2B6*6 (rs2279343), and NAT2 (rs1041983) were significantly higher in Thais than East-Asian and global populations. 121 unreported variants had potential to exert clinical impact, majority were rare and population-specific, with 60.3% of variants absent from gnomAD database. This study demonstrates the population-specific variation in clinically relevant pharmacogenes, the importance of CYP2D6 structural variation detection in the Thai population, and potential of unreported variants in explaining drug response. These findings are essential in development of dosing guidelines, PGx testing, clinical trials, and drugs.

Ethnic-specific BRCA1/2 variation within Asia population: evidence from over 78 000 cancer and 40 000 non-cancer cases of Indian, Chinese, Korean and Japanese populations

BACKGROUND

Germline mutation in BRCA1 and BRCA2 (BRCA) is genetic predisposition for breast and ovarian cancer. Identification of mutation carriers is a critical step to prevent and treat the cancer in the mutation carriers. Human BRCA variation has been well determined as ethnic-specific by studies in Ashkenazi Jewish, Polish and Icelandic populations in the 1990s. However, sufficient evidence is lacking to determine if ethnic-specific BRCA variation is also present in Asia population, which is the largest and the most diversified in modern humans. Our current study aims to investigate ethnic-specific BRCA variation in Asian population.

METHODS

We performed a comprehensive data mining to collect BRCA variation data in Indian, Chinese, Korean and Japanese populations derived from over 78 000 cancer and 40 000 non-cancer cases. We standardised all BRCA variation data following the international standard. We made a systematic comparison between the datasets including variant composition, variation spectrum, variant type, clinical class, founder mutation and high-frequent variants.

RESULTS

Our analysis showed that over half of the Asian BRCA variants were Asian-specific, and significant differences were present between the four Asia populations in each category analysed.

CONCLUSION

Data from our study reveal that ethnic-specific BRCA variation is commonly present in Asia population as existing in non-Asian populations. Our study indicates that ethnicity should be an important factor to consider in prevention and treatment of BRCA mutation-related cancer in the Asia population. We recommend that the current BRCA variation databases should include ethnic variation information in order to function as true global BRCA references.

An integrated personal and population-based Egyptian genome reference

A small number of de novo assembled human genomes have been reported to date, and few have been complemented with population-based genetic variation, which is particularly important for North Africa, a region underrepresented in current genome-wide references. Here, we combine long- and short-read whole-genome sequencing data with recent assembly approaches into a de novo assembly of an Egyptian genome. The assembly demonstrates well-balanced quality metrics and is complemented with variant phasing via linked reads into haploblocks, which we associate with gene expression changes in blood. To construct an Egyptian genome reference, we identify genome-wide genetic variation within a cohort of 110 Egyptian individuals. We show that differences in allele frequencies and linkage disequilibrium between Egyptians and Europeans may compromise the transferability of European ancestry-based genetic disease risk and polygenic scores, substantiating the need for multi-ethnic genome references. Thus, the Egyptian genome reference will be a valuable resource for precision medicine.

Frequency spectrum of rare and clinically relevant markers in multiethnic Indian populations (ClinIndb): A resource for genomic medicine in India

There have been concerted efforts toward cataloging rare and deleterious variants in different world populations using high-throughput genotyping and sequencing-based methods. The Indian population is underrepresented or its information with respect to clinically relevant variants is sparse in public data sets. The aim of this study was to estimate the burden of monogenic disease-causing variants in Indian populations. Toward this, we have assessed the frequency profile of monogenic phenotype-associated ClinVar variants. The study utilized a genotype data set (global screening array, Illumina) from 2795 individuals (multiple in-house genomics cohorts) representing diverse ethnic and geographically distinct Indian populations. Of the analyzed variants from Global Screening Array, ~9% were found to be informative and were either not known earlier or underrepresented in public databases in terms of their frequencies. These variants were linked to disorders, namely inborn errors of metabolism, monogenic diabetes, hereditary cancers, and various other hereditary conditions. We have also shown that our study cohort is genetically a better representative of the Indian population than its representation in the 1000 Genome Project (South Asians). We have created a database, ClinIndb, linked to the Leiden Open Variation Database, to help clinicians and researchers in diagnosis, counseling, and development of appropriate genetic screening tools relevant to the Indian populations and Indians living abroad.

Prospective avenues for human population genomics and disease mapping in southern Africa

Population substructure within human populations is globally evident and a well-known confounding factor in many genetic studies. In contrast, admixture mapping exploits population stratification to detect genotype-phenotype correlations in admixed populations. Southern Africa has untapped potential for disease mapping of ancestry-specific disease risk alleles due to the distinct genetic diversity in its populations compared to other populations worldwide. This diversity contributes to a number of phenotypes, including ancestry-specific disease risk and response to pathogens. Although the 1000 Genomes Project significantly improved our understanding of genetic variation globally, southern African populations are still severely underrepresented in biomedical and human genetic studies due to insufficient large-scale publicly available data. In addition to a lack of genetic data in public repositories, existing software, algorithms and resources used for imputation and phasing of genotypic data (amongst others) are largely ineffective for populations with a complex genetic architecture such as that seen in southern Africa. This review article, therefore, aims to summarise the current limitations of conducting genetic studies on populations with a complex genetic architecture to identify potential areas for further research and development.

The road ahead in genetics and genomics

In celebration of the 20th anniversary of Nature Reviews Genetics, we asked 12 leading researchers to reflect on the key challenges and opportunities faced by the field of genetics and genomics. Keeping their particular research area in mind, they take stock of the current state of play and emphasize the work that remains to be done over the next few years so that, ultimately, the benefits of genetic and genomic research can be felt by everyone.

To celebrate the first 20 years of Nature Reviews Genetics, we asked 12 leading scientists to reflect on the key challenges and opportunities faced by the field of genetics and genomics.

Amy L. McGuire is the Leon Jaworski Professor of Biomedical Ethics and Director of the Center for Medical Ethics and Health Policy at Baylor College of Medicine. She has received numerous teaching awards at Baylor College of Medicine, was recognized by the Texas Executive Women as a Woman on the Move in 2016 and was invited to give a TedMed talk titled “There is No Genome for the Human Spirit” in 2014. In 2020, she was elected as a Hastings Center Fellow. Her research focuses on ethical and policy issues related to emerging technologies, with a particular focus on genomic research, personalized medicine and the clinical integration of novel neurotechnologies.

Stacey Gabriel is the Senior Director of the Genomics Platform at the Broad Institute since 2012 and has led platform development, execution and operation since its founding. She is Chair of Institute Scientists and serves on the institute’s executive leadership team. She is widely recognized as a leader in genomic technology and project execution. She has led the Broad’s contributions to numerous flagship projects in human genetics, including the International HapMap Project, the 1000 Genomes Project, The Cancer Genome Atlas, the National Heart, Lung, and Blood Institute’s Exome Sequencing Project and the TOPMed programme. She is Principal Investigator of the Broad’s All of Us (AoU) Genomics Center and serves on the AoU Program Steering Committee.

Sarah A. Tishkoff is the David and Lyn Silfen University Associate Professor in Genetics and Biology at the University of Pennsylvania, Philadelphia, USA, and holds appointments in the School of Medicine and the School of Arts and Sciences. She is a member of the US National Academy of Sciences and a recipient of an NIH Pioneer Award, a David and Lucile Packard Career Award, a Burroughs/Wellcome Fund Career Award and an American Society of Human Genetics Curt Stern Award. Her work focuses on genomic variation in Africa, human evolutionary history, the genetic basis of adaptation and phenotypic variation in Africa, and the genetic basis of susceptibility to infectious disease in Africa.

Ambroise Wonkam is Professor of Medical Genetics, Director of GeneMAP (Genetic Medicine of African Populations Research Centre) and Deputy Dean Research in the Faculty of Health Sciences, University of Cape Town, South Africa. He has successfully led numerous NIH- and Wellcome Trust-funded projects over the past decade to investigate clinical variability in sickle cell disease, hearing impairment genetics and the return of individual findings in genetic research in Africa. He won the competitive Clinical Genetics Society International Award for 2014 from the British Society of Genetic Medicine. He is president of the African Society of Human Genetics.

Aravinda Chakravarti is Director of the Center for Human Genetics and Genomics, the Muriel G. and George W. Singer Professor of Neuroscience and Physiology, and Professor of Medicine at New York University School of Medicine. He is an elected member of the US National Academy of Sciences, the US National Academy of Medicine and the Indian National Science Academy. He has been a key participant in the Human Genome Project, the International HapMap Project and the 1000 Genomes Project. His research attempts to understand the molecular basis of multifactorial disease. He was awarded the 2013 William Allan Award by the American Society of Human Genetics and the 2018 Chen Award by the Human Genome Organization.

Eileen E. M. Furlong is Head of the Genome Biology Department at the European Molecular Biology Laboratory (EMBL) and a member of the EMBL Directorate. She is an elected member of the European Molecular Biology Organization (EMBO) and the Academia Europaea, and a European Research Council (ERC) advanced investigator. Her group dissects fundamental principles of how the genome is regulated and how it drives cell fate decisions during embryonic development, including how developmental enhancers are organized and function within the 3D nucleus. Her work combines genetics, (single-cell) genomics, imaging and computational approaches to understand these processes. Her research has advanced the development of genomic methods for use in complex multicellular organisms.

Barbara Treutlein is Associate Professor of Quantitative Developmental Biology in the Department of Biosystems Science and Engineering of ETH Zurich in Basel, Switzerland. Her group uses and develops single-cell genomics approaches in combination with stem cell-based 2D and 3D culture systems to study how human organs develop and regenerate and how cell fate is regulated. For her work, Barbara has received multiple awards, including the Friedmund Neumann Prize of the Schering Foundation, the Dr. Susan Lim Award for Outstanding Young Investigator of the International Society of Stem Cell Research and the EMBO Young Investigator Award.

Alexander Meissner is a scientific member of the Max Planck Society and currently Managing Director of the Max Planck Institute (MPI) for Molecular Genetics in Berlin, Germany. He heads the Department of Genome Regulation and is a visiting scientist in the Department of Stem Cell and Regenerative Biology at Harvard University. Before his move to the MPI, he was a tenured professor at Harvard University and a senior associate member of the Broad Institute, where he co-directed the epigenomics programme. In 2018, he was elected as an EMBO member. His laboratory uses genomic tools to study developmental and disease biology with a particular focus on epigenetic regulation.

Howard Y. Chang is the Virginia and D. K. Ludwig Professor of Cancer Genomics at Stanford University and an investigator at the Howard Hughes Medical Institute. He is a physician–scientist who has focused on deciphering the hidden information in the non-coding genome. His laboratory is best known for studies of long non-coding RNAs in gene regulation and development of new epigenomic technologies. He is an elected member of the US National Academy of Sciences, the US National Academy of Medicine, and the American Academy of Arts and Sciences.

Núria López-Bigas is ICREA research Professor at the Institute for Research in Biomedicine and Associate Professor at the University Pompeu Fabra. She obtained an ERC Consolidator Grant in 2015 and was elected as an EMBO member in 2016. Her work has been recognized with the prestigious Banc de Sabadell Award for Research in Biomedicine, the Catalan National Award for Young Research Talent and the Career Development Award from the Human Frontier Science Program. Her research focuses on the identification of cancer driver mutations, genes and pathways across tumour types and in understanding the mutational processes that lead to the accumulation of mutations in cancer cells.

Eran Segal is Professor in the Department of Computer Science and Applied Mathematics at the Weizmann Institute of Science, heading a multidisciplinary laboratory with extensive experience in machine learning, computational biology and analysis of heterogeneous high-throughput genomic data. His research focuses on the microbiome, nutrition and genetics, and their effect on health and disease and aims to develop personalized medicine based on big data from human cohorts. He has published more than 150 publications and received several awards and honours for his work, including the Overton and the Michael Bruno awards. He was recently elected as an EMBO member and as a member of the Israel Young Academy.

Jin-Soo Kim is Director of the Center for Genome Engineering in the Institute for Basic Science in Daejon, South Korea. He has received numerous awards, including the 2017 Asan Award in Medicine, the 2017 Yumin Award in Science and the 2019 Research Excellence Award (Federation of Asian and Oceanian Biochemists and Molecular Biologists). He was featured as one of ten Science Stars of East Asia in Nature (558, 502–510 (2018)) and has been recognized as a highly cited researcher by Clarivate Analytics since 2018. His work focuses on developing tools for genome editing in biomedical research.

Polygenic background modifies penetrance of monogenic variants for tier 1 genomic conditions

Genetic variation can predispose to disease both through (i) monogenic risk variants that disrupt a physiologic pathway with large effect on disease and (ii) polygenic risk that involves many variants of small effect in different pathways. Few studies have explored the interplay between monogenic and polygenic risk. Here, we study 80,928 individuals to examine whether polygenic background can modify penetrance of disease in tier 1 genomic conditions — familial hypercholesterolemia, hereditary breast and ovarian cancer, and Lynch syndrome. Among carriers of a monogenic risk variant, we estimate substantial gradients in disease risk based on polygenic background — the probability of disease by age 75 years ranged from 17% to 78% for coronary artery disease, 13% to 76% for breast cancer, and 11% to 80% for colon cancer. We propose that accounting for polygenic background is likely to increase accuracy of risk estimation for individuals who inherit a monogenic risk variant.

Genetic variation predisposes to disease via monogenic and polygenic risk variants. Here, the authors assess the interplay between these types of variation on disease penetrance in 80,928 individuals. In carriers of monogenic variants, they show that disease risk is a gradient influenced by polygenic background.

Transcriptional adaptation: a mechanism underlying genetic robustness

Mutations play a crucial role in evolution as they provide the genetic variation that allows evolutionary change. Although some mutations in regulatory elements or coding regions can be beneficial, a large number of them disrupt gene function and reduce fitness. Organisms utilize several mechanisms to compensate for the damaging consequences of genetic perturbations. One such mechanism is the recently identified process of transcriptional adaptation (TA): during this event, mutations that cause mutant mRNA degradation trigger the transcriptional modulation of so-called adapting genes. In some cases, for example when one (or more) of the upregulated genes is functionally redundant with the mutated gene, this process compensates for the loss of the mutated gene's product. Notably, unlike other mechanisms underlying genetic robustness, TA is not triggered by the loss of protein function, an observation that has prompted studies into the machinery of TA and the contexts in which it functions. Here, we review the discovery and current understanding of TA, and discuss how its main features appear to be conserved across species. In light of these findings, we also speculate on the importance of TA in the context of human disease, and provide some recommendations for genome-editing strategies that should be more effective.

Validation of a Genome-Wide Polygenic Score for Coronary Artery Disease in South Asians

BACKGROUND

Genome-wide polygenic scores (GPS) integrate information from many common DNA variants into a single number. Because rates of coronary artery disease (CAD) are substantially higher among South Asians, a GPS to identify high-risk individuals may be particularly useful in this population.

OBJECTIVES

This analysis used summary statistics from a prior genome-wide association study to derive a new GPSCAD for South Asians.

METHODS

This GPSCAD was validated in 7,244 South Asian UK Biobank participants and tested in 491 individuals from a case-control study in Bangladesh. Next, a static ancestry and GPSCAD reference distribution was built using whole-genome sequencing from 1,522 Indian individuals, and a framework was tested for projecting individuals onto this static ancestry and GPSCAD reference distribution using 1,800 CAD cases and 1,163 control subjects newly recruited in India.

RESULTS

The GPSCAD, containing 6,630,150 common DNA variants, had an odds ratio (OR) per SD of 1.58 in South Asian UK Biobank participants and 1.60 in the Bangladeshi study (p < 0.001 for each). Next, individuals of the Indian case-control study were projected onto static reference distributions, observing an OR/SD of 1.66 (p < 0.001). Compared with the middle quintile, risk for CAD was most pronounced for those in the top 5% of the GPSCAD distribution-ORs of 4.16, 2.46, and 3.22 in the South Asian UK Biobank, Bangladeshi, and Indian studies, respectively (p < 0.05 for each).

CONCLUSIONS

The new GPSCAD has been developed and tested using 3 distinct South Asian studies, and provides a generalizable framework for ancestry-specific GPS assessment.

Current Status on Population Genome Catalogues in different Countries

Genomics has become indispensable for research in the last two decades. Completed and ongoing genome projects such as the UK Biobank, ICGC, TCGA and GenomeAsia100K have helped to understand several life-threatening diseases like cancer. Such initiatives from different countries have offered genomics-based diagnostics along with glues for therapies towards personalized healthcare. The Indian Agencies has started initiatives to catalogue the genome sequences of 20,000 individuals. The Department of Biotechnology (DBT) along with other scientific agencies has plans to sequence 10,000 healthy individuals and 10,000 diseased individuals. The Council of Scientific and Industrial Research (CSIR) also developed the "IndiGen" genome project where genome sequences for 1008 individuals are made available in Phase I. This will enable the development of a genome catalogue to introduce novel genomics-based clinical applications in future healthcare plan.

Making neurogenetics a global endeavour

This scientific commentary refers to ‘The role of genetics in Parkinson’s disease: a large cohort study in Chinese mainland population’, by Zhao et al. (doi:10.1093/brain/awaa167).

A virtual assessment of the suprainiac depressions on the Eyasi I (Tanzania) and Aduma ADU-VP-1/3 (Ethiopia) Pleistocene hominin crania

Despite a steady increase in our understanding of the phenotypic variation of Pleistocene Homo, debate continues over phylogenetically informative features. One such trait is the suprainiac fossa, a depression on the occipital bone above inion that is commonly considered an autapomorphy of the Neanderthal lineage. Challenging this convention, depressions in the suprainiac region have also been described for two Pleistocene hominin crania from sub-Saharan Africa: Eyasi I (Tanzania) and ADU-VP-1/3 (Ethiopia). Here, we use a combined quantitative and qualitative approach, using μCT imaging, to investigate the occipital depressions on these specimens. The results show that neither the external nor the internal morphologies of these depressions bear any resemblance to the Neanderthal condition. A principal component analysis based on multiple thickness measurements along the occipital squama demonstrates that the relative thickness values for the internal structures in Eyasi I and ADU-VP-1/3 are within the range of Homo sapiens. Thus, our results support the autapomorphic status of the Neanderthal suprainiac fossa and highlight the need to use nuanced approaches and multiple lines of evidence.

Methods for detecting introgressed archaic sequences

Analysis of genome sequences from archaic and modern humans have revealed multiple episodes of admixture between highly-diverged population groups. Statistical methods that attempt to localize DNA segments introduced by these events offer a powerful tool to investigate recent human evolution. We review recent advances in methods for detecting introgressed sequences.