Ancient genomes reveal origin and rapid trans-Eurasian migration of 7th century Avar elites.

The Avars settled the Carpathian Basin in 567/68 CE, establishing an empire lasting over 200 years. Who they were and where they came from is highly debated. Contemporaries have disagreed about whether they were, as they claimed, the direct successors of the Mongolian Steppe Rouran empire that was destroyed by the Turks in ∼550 CE. Here, we analyze new genome-wide data from 66 pre-Avar and Avar-period Carpathian Basin individuals, including the 8 richest Avar-period burials and further elite sites from Avar's empire core region. Our results provide support for a rapid long-distance trans-Eurasian migration of Avar-period elites. These individuals carried Northeast Asian ancestry matching the profile of preceding Mongolian Steppe populations, particularly a genome available from the Rouran period. Some of the later elite individuals carried an additional non-local ancestry component broadly matching the steppe, which could point to a later migration or reflect greater genetic diversity within the initial migrant population.

Underrepresentation of Asian awardees of United States biomedical research prizes.

The lack of racial diversity among the winners of United States biomedical research prizes reflects a chronic problem of the underappreciation of certain groups of biomedical scientists. Asians continue to be severely underrepresented as awardees of United States biomedical research prizes, a trend that shows no obvious recent improvement.

Limb development genes underlie variation in human fingerprint patterns.

Fingerprints are of long-standing practical and cultural interest, but little is known about the mechanisms that underlie their variation. Using genome-wide scans in Han Chinese cohorts, we identified 18 loci associated with fingerprint type across the digits, including a genetic basis for the long-recognized "pattern-block" correlations among the middle three digits. In particular, we identified a variant near EVI1 that alters regulatory activity and established a role for EVI1 in dermatoglyph patterning in mice. Dynamic EVI1 expression during human development supports its role in shaping the limbs and digits, rather than influencing skin patterning directly. Trans-ethnic meta-analysis identified 43 fingerprint-associated loci, with nearby genes being strongly enriched for general limb development pathways. We also found that fingerprint patterns were genetically correlated with hand proportions. Taken together, these findings support the key role of limb development genes in influencing the outcome of fingerprint patterning.

TMEM41B Is a Pan-flavivirus Host Factor.

Flaviviruses pose a constant threat to human health. These RNA viruses are transmitted by the bite of infected mosquitoes and ticks and regularly cause outbreaks. To identify host factors required for flavivirus infection, we performed full-genome loss of function CRISPR-Cas9 screens. Based on these results, we focused our efforts on characterizing the roles that TMEM41B and VMP1 play in the virus replication cycle. Our mechanistic studies on TMEM41B revealed that all members of the Flaviviridae family that we tested require TMEM41B. We tested 12 additional virus families and found that SARS-CoV-2 of the Coronaviridae also required TMEM41B for infection. Remarkably, single nucleotide polymorphisms present at nearly 20% in East Asian populations reduce flavivirus infection. Based on our mechanistic studies, we propose that TMEM41B is recruited to flavivirus RNA replication complexes to facilitate membrane curvature, which creates a protected environment for viral genome replication.

Trans-ethnic and Ancestry-Specific Blood-Cell Genetics in 746,667 Individuals from 5 Global Populations.

Most loci identified by GWASs have been found in populations of European ancestry (EUR). In trans-ethnic meta-analyses for 15 hematological traits in 746,667 participants, including 184,535 non-EUR individuals, we identified 5,552 trait-variant associations at p < 5 × 10-9, including 71 novel associations not found in EUR populations. We also identified 28 additional novel variants in ancestry-specific, non-EUR meta-analyses, including an IL7 missense variant in South Asians associated with lymphocyte count in vivo and IL-7 secretion levels in vitro. Fine-mapping prioritized variants annotated as functional and generated 95% credible sets that were 30% smaller when using the trans-ethnic as opposed to the EUR-only results. We explored the clinical significance and predictive value of trans-ethnic variants in multiple populations and compared genetic architecture and the effect of natural selection on these blood phenotypes between populations. Altogether, our results for hematological traits highlight the value of a more global representation of populations in genetic studies.

Integrative Proteomic Characterization of Human Lung Adenocarcinoma.

Genomic studies of lung adenocarcinoma (LUAD) have advanced our understanding of the disease's biology and accelerated targeted therapy. However, the proteomic characteristics of LUAD remain poorly understood. We carried out a comprehensive proteomics analysis of 103 cases of LUAD in Chinese patients. Integrative analysis of proteome, phosphoproteome, transcriptome, and whole-exome sequencing data revealed cancer-associated characteristics, such as tumor-associated protein variants, distinct proteomics features, and clinical outcomes in patients at an early stage or with EGFR and TP53 mutations. Proteome-based stratification of LUAD revealed three subtypes (S-I, S-II, and S-III) related to different clinical and molecular features. Further, we nominated potential drug targets and validated the plasma protein level of HSP 90ß as a potential prognostic biomarker for LUAD in an independent cohort. Our integrative proteomics analysis enables a more comprehensive understanding of the molecular landscape of LUAD and offers an opportunity for more precise diagnosis and treatment.

An Ancient Harappan Genome Lacks Ancestry from Steppe Pastoralists or Iranian Farmers.

We report an ancient genome from the Indus Valley Civilization (IVC). The individual we sequenced fits as a mixture of people related to ancient Iranians (the largest component) and Southeast Asian hunter-gatherers, a unique profile that matches ancient DNA from 11 genetic outliers from sites in Iran and Turkmenistan in cultural communication with the IVC. These individuals had little if any Steppe pastoralist-derived ancestry, showing that it was not ubiquitous in northwest South Asia during the IVC as it is today. The Iranian-related ancestry in the IVC derives from a lineage leading to early Iranian farmers, herders, and hunter-gatherers before their ancestors separated, contradicting the hypothesis that the shared ancestry between early Iranians and South Asians reflects a large-scale spread of western Iranian farmers east. Instead, sampled ancient genomes from the Iranian plateau and IVC descend from different groups of hunter-gatherers who began farming without being connected by substantial movement of people.

Large-Scale Whole-Genome Sequencing of Three Diverse Asian Populations in Singapore.

Underrepresentation of Asian genomes has hindered population and medical genetics research on Asians, leading to population disparities in precision medicine. By whole-genome sequencing of 4,810 Singapore Chinese, Malays, and Indians, we found 98.3 million SNPs and small insertions or deletions, over half of which are novel. Population structure analysis demonstrated great representation of Asian genetic diversity by three ethnicities in Singapore and revealed a Malay-related novel ancestry component. Furthermore, demographic inference suggested that Malays split from Chinese ∼24,800 years ago and experienced significant admixture with East Asians ∼1,700 years ago, coinciding with the Austronesian expansion. Additionally, we identified 20 candidate loci for natural selection, 14 of which harbored robust associations with complex traits and diseases. Finally, we show that our data can substantially improve genotype imputation in diverse Asian and Oceanian populations. These results highlight the value of our data as a resource to empower human genetics discovery across broad geographic regions.

Multiple Deeply Divergent Denisovan Ancestries in Papuans.

Genome sequences are known for two archaic hominins-Neanderthals and Denisovans-which interbred with anatomically modern humans as they dispersed out of Africa. We identified high-confidence archaic haplotypes in 161 new genomes spanning 14 island groups in Island Southeast Asia and New Guinea and found large stretches of DNA that are inconsistent with a single introgressing Denisovan origin. Instead, modern Papuans carry hundreds of gene variants from two deeply divergent Denisovan lineages that separated over 350 thousand years ago. Spatial and temporal structure among these lineages suggest that introgression from one of these Denisovan groups predominantly took place east of the Wallace line and continued until near the end of the Pleistocene. A third Denisovan lineage occurs in modern East Asians. This regional mosaic suggests considerable complexity in archaic contact, with modern humans interbreeding with multiple Denisovan groups that were geographically isolated from each other over deep evolutionary time.

Analysis of Human Sequence Data Reveals Two Pulses of Archaic Denisovan Admixture.

Anatomically modern humans interbred with Neanderthals and with a related archaic population known as Denisovans. Genomes of several Neanderthals and one Denisovan have been sequenced, and these reference genomes have been used to detect introgressed genetic material in present-day human genomes. Segments of introgression also can be detected without use of reference genomes, and doing so can be advantageous for finding introgressed segments that are less closely related to the sequenced archaic genomes. We apply a new reference-free method for detecting archaic introgression to 5,639 whole-genome sequences from Eurasia and Oceania. We find Denisovan ancestry in populations from East and South Asia and Papuans. Denisovan ancestry comprises two components with differing similarity to the sequenced Altai Denisovan individual. This indicates that at least two distinct instances of Denisovan admixture into modern humans occurred, involving Denisovan populations that had different levels of relatedness to the sequenced Altai Denisovan. VIDEO ABSTRACT.

US Immigration Westernizes the Human Gut Microbiome.

Many US immigrant populations develop metabolic diseases post immigration, but the causes are not well understood. Although the microbiome plays a role in metabolic disease, there have been no studies measuring the effects of US immigration on the gut microbiome. We collected stool, dietary recalls, and anthropometrics from 514 Hmong and Karen individuals living in Thailand and the United States, including first- and second-generation immigrants and 19 Karen individuals sampled before and after immigration, as well as from 36 US-born European American individuals. Using 16S and deep shotgun metagenomic DNA sequencing, we found that migration from a non-Western country to the United States is associated with immediate loss of gut microbiome diversity and function in which US-associated strains and functions displace native strains and functions. These effects increase with duration of US residence and are compounded by obesity and across generations.

Genomic Analyses from Non-invasive Prenatal Testing Reveal Genetic Associations, Patterns of Viral Infections, and Chinese Population History.

We analyze whole-genome sequencing data from 141,431 Chinese women generated for non-invasive prenatal testing (NIPT). We use these data to characterize the population genetic structure and to investigate genetic associations with maternal and infectious traits. We show that the present day distribution of alleles is a function of both ancient migration and very recent population movements. We reveal novel phenotype-genotype associations, including several replicated associations with height and BMI, an association between maternal age and EMB, and between twin pregnancy and NRG1. Finally, we identify a unique pattern of circulating viral DNA in plasma with high prevalence of hepatitis B and other clinically relevant maternal infections. A GWAS for viral infections identifies an exceptionally strong association between integrated herpesvirus 6 and MOV10L1, which affects piwi-interacting RNA (piRNA) processing and PIWI protein function. These findings demonstrate the great value and potential of accumulating NIPT data for worldwide medical and genetic analyses.

Physiological and Genetic Adaptations to Diving in Sea Nomads.

Understanding the physiology and genetics of human hypoxia tolerance has important medical implications, but this phenomenon has thus far only been investigated in high-altitude human populations. Another system, yet to be explored, is humans who engage in breath-hold diving. The indigenous Bajau people ("Sea Nomads") of Southeast Asia live a subsistence lifestyle based on breath-hold diving and are renowned for their extraordinary breath-holding abilities. However, it is unknown whether this has a genetic basis. Using a comparative genomic study, we show that natural selection on genetic variants in the PDE10A gene have increased spleen size in the Bajau, providing them with a larger reservoir of oxygenated red blood cells. We also find evidence of strong selection specific to the Bajau on BDKRB2, a gene affecting the human diving reflex. Thus, the Bajau, and possibly other diving populations, provide a new opportunity to study human adaptation to hypoxia tolerance. VIDEO ABSTRACT.

African-specific molecular taxonomy of prostate cancer.

Prostate cancer is characterized by considerable geo-ethnic disparity. African ancestry is a significant risk factor, with mortality rates across sub-Saharan Africa of 2.7-fold higher than global averages1. The contributing genetic and non-genetic factors, and associated mutational processes, are unknown2,3. Here, through whole-genome sequencing of treatment-naive prostate cancer samples from 183 ancestrally (African versus European) and globally distinct patients, we generate a large cancer genomics resource for sub-Saharan Africa, identifying around 2 million somatic variants. Significant African-ancestry-specific findings include an elevated tumour mutational burden, increased percentage of genome alteration, a greater number of predicted damaging mutations and a higher total of mutational signatures, and the driver genes NCOA2, STK19, DDX11L1, PCAT1 and SETBP1. Examining all somatic mutational types, we describe a molecular taxonomy for prostate cancer differentiated by ancestry and defined as global mutational subtypes (GMS). By further including Chinese Asian data, we confirm that GMS-B (copy-number gain) and GMS-D (mutationally noisy) are specific to African populations, GMS-A (mutationally quiet) is universal (all ethnicities) and the African-European-restricted subtype GMS-C (copy-number losses) predicts poor clinical outcomes. In addition to the clinical benefit of including individuals of African ancestry, our GMS subtypes reveal different evolutionary trajectories and mutational processes suggesting that both common genetic and environmental factors contribute to the disparity between ethnicities. Analogous to gene-environment interaction-defined here as a different effect of an environmental surrounding in people with different ancestries or vice versa-we anticipate that GMS subtypes act as a proxy for intrinsic and extrinsic mutational processes in cancers, promoting global inclusion in landmark studies.

Genes with differential expression across ancestries are enriched in ancestry-specific disease effects likely due to gene-by-environment interactions.

Multi-ancestry genome-wide association studies (GWASs) have highlighted the existence of variants with ancestry-specific effect sizes. Understanding where and why these ancestry-specific effects occur is fundamental to understanding the genetic basis of human diseases and complex traits. Here, we characterized genes differentially expressed across ancestries (ancDE genes) at the cell-type level by leveraging single-cell RNA-sequencing data in peripheral blood mononuclear cells for 21 individuals with East Asian (EAS) ancestry and 23 individuals with European (EUR) ancestry (172,385 cells); then, we tested whether variants surrounding those genes were enriched in disease variants with ancestry-specific effect sizes by leveraging ancestry-matched GWASs of 31 diseases and complex traits (average n ∼ 90,000 and ∼ 267,000 in EAS and EUR, respectively). We observed that ancDE genes tended to be cell-type specific and enriched in genes interacting with the environment and in variants with ancestry-specific disease effect sizes, which suggests cell-type-specific, gene-by-environment interactions shared between regulatory and disease architectures. Finally, we illustrated how different environments might have led to ancestry-specific myeloid cell leukemia 1 (MCL1) expression in B cells and ancestry-specific allele effect sizes in lymphocyte count GWASs for variants surrounding MCL1. Our results imply that large single-cell and GWAS datasets from diverse ancestries are required to improve our understanding of human diseases.

Concurrent Asian monsoon strengthening and early modern human dispersal to East Asia during the last interglacial.

The relationship between initial Homo sapiens dispersal from Africa to East Asia and the orbitally paced evolution of the Asian summer monsoon (ASM)-currently the largest monsoon system-remains underexplored due to lack of coordinated synthesis of both Asian paleoanthropological and paleoclimatic data. Here, we investigate orbital-scale ASM dynamics during the last 280 thousand years (kyr) and their likely influences on early H. sapiens dispersal to East Asia, through a unique integration of i) new centennial-resolution ASM records from the Chinese Loess Plateau, ii) model-based East Asian hydroclimatic reconstructions, iii) paleoanthropological data compilations, and iv) global H. sapiens habitat suitability simulations. Our combined proxy- and model-based reconstructions suggest that ASM precipitation responded to a combination of Northern Hemisphere ice volume, greenhouse gas, and regional summer insolation forcing, with cooccurring primary orbital cycles of ~100-kyr, 41-kyr, and ~20-kyr. Between ~125 and 70 kyr ago, summer monsoon rains and temperatures increased in vast areas across Asia. This episode coincides with the earliest H. sapiens fossil occurrence at multiple localities in East Asia. Following the transcontinental increase in simulated habitat suitability, we suggest that ASM strengthening together with Southeast African climate deterioration may have promoted the initial H. sapiens dispersal from their African homeland to remote East Asia during the last interglacial.

Genetic variation in CCDC93 is associated with elevated central systolic blood pressure, impaired arterial relaxation, and mitochondrial dysfunction.

Genetic studies of blood pressure (BP) traits to date have been performed on conventional measures by brachial cuff sphygmomanometer for systolic BP (SBP) and diastolic BP, integrating several physiologic occurrences. Genetic associations with central SBP (cSBP) have not been well-studied. Genetic discovery studies of BP have been most often performed in European-ancestry samples. Here, we investigated genetic associations with cSBP in a Chinese population and functionally validated the impact of a novel associated coiled-coil domain containing 93 (CCDC93) gene on BP regulation. An exome-wide association study (EWAS) was performed using a mixed linear model of non-invasive cSBP and peripheral BP traits in a Han Chinese population (N = 5,954) from Beijing, China genotyped with a customized Illumina ExomeChip array. We identified four SNP-trait associations with three SNPs, including two novel associations (rs2165468-SBP and rs33975708-cSBP). rs33975708 is a coding variant in the CCDC93 gene, c.535C>T, p.Arg179Cys (MAF = 0.15%), and was associated with increased cSBP (ß = 29.3 mmHg, P = 1.23x10-7). CRISPR/Cas9 genome editing was used to model the effect of Ccdc93 loss in mice. Homozygous Ccdc93 deletion was lethal prior to day 10.5 of embryonic development. Ccdc93+/- heterozygous mice were viable and morphologically normal, with 1.3-fold lower aortic Ccdc93 protein expression (P = 0.0041) and elevated SBP as compared to littermate Ccdc93+/+ controls (110±8 mmHg vs 125±10 mmHg, P = 0.016). Wire myography of Ccdc93+/- aortae showed impaired acetylcholine-induced relaxation and enhanced phenylephrine-induced contraction. RNA-Seq transcriptome analysis of Ccdc93+/- mouse thoracic aortae identified significantly enriched pathways altered in fatty acid metabolism and mitochondrial metabolism. Plasma free fatty acid levels were elevated in Ccdc93+/- mice (96±7mM vs 124±13mM, P = 0.0031) and aortic mitochondrial dysfunction was observed through aberrant Parkin and Nix protein expression. Together, our genetic and functional studies support a novel role of CCDC93 in the regulation of BP through its effects on vascular mitochondrial function and endothelial function.

Arranged cohabitation among Chinese Muslims.

While modern family-related ideas and behaviors have become more widely accepted in contemporary China, Chinese Muslim minorities continue to hold on to traditional religious practices. Surprisingly, data from our survey conducted in Gansu province in China's northwestern borderlands reveal that Muslims of the Hui and Dongxiang ethnicities reported much higher rates of cohabitation experience than the secular majority Han. Based on follow-up qualitative interviews, we found the answer to lie in the interplay between the highly interventionist Chinese state and the robust cultural resilience of local Islamic communities. While the state maintains a high minimum legal age of marriage, the early marriage norm remains strong in Chinese Muslim communities, where religion constitutes an alternative and often more powerful source of legitimacy-at least in the private sphere of life. Using the 2000 census data, we further show that women in almost all 10 Muslim ethnic groups have higher percentages of underage births and premarital births than Han women, both nationally and in the northwest where most Chinese Muslims live. As the once-outlawed behavior of cohabitation became more socially acceptable during the reform and opening-up era, young Muslim Chinese often found themselves in "arranged cohabitations" as de facto marriages formed at younger-than-legal ages. In doing so, Chinese Muslim communities have reinvented the meaning of cohabitation. Rather than liberal intimate relationship based on individual autonomy, cohabitation has served as a coping strategy by which Islamic patriarchs circumvent the Chinese state's aggressive regulations aimed at "modernizing" the Muslim family.

The genetic admixture and assimilation of Ahom: a historic migrant from Thailand to India.

The Northeastern region of India is considered a gateway for modern humans' dispersal throughout Asia. This region is a mixture of various ethnic and indigenous populations amalgamating multiple ancestries. One reason for such amalgamation is that, South Asia experienced multiple historic migrations from various parts of the world. A few examples explored genetically are Jews, Parsis and Siddis. Ahom is a dynasty that historically migrated to India during the 12th century. However, this putative migration has not been studied genetically at high resolution. Therefore, to validate this historical evidence, we genotyped autosomal data of the Modern Ahom population residing in seven sister states of India. Principal Component and Admixture analyses haave suggested a substantial admixture of the Ahom population with the local Tibeto-Burman populations. Moreover, the haplotype-based analysis has linked these Ahom individuals mainly with the Kusunda (a language isolated from Nepal) and Khasi (an Austroasiatic population of Meghalaya). Such unexpected presence of widespread population affinities suggests that Ahom mixed and assimilated a wide variety of Trans-Himalayan populations inhabiting this region after the migration. In summary, we observed a significant deviation of Ahom from their ancestral homeland (Thailand) and extensive admixture and assimilation with the local South Asian populations.

Power of inclusion: Enhancing polygenic prediction with admixed individuals.

Admixed individuals offer unique opportunities for addressing limited transferability in polygenic scores (PGSs), given the substantial trans-ancestry genetic correlation in many complex traits. However, they are rarely considered in PGS training, given the challenges in representing ancestry-matched linkage-disequilibrium reference panels for admixed individuals. Here we present inclusive PGS (iPGS), which captures ancestry-shared genetic effects by finding the exact solution for penalized regression on individual-level data and is thus naturally applicable to admixed individuals. We validate our approach in a simulation study across 33 configurations with varying heritability, polygenicity, and ancestry composition in the training set. When iPGS is applied to n = 237,055 ancestry-diverse individuals in the UK Biobank, it shows the greatest improvements in Africans by 48.9% on average across 60 quantitative traits and up to 50-fold improvements for some traits (neutrophil count, R2 = 0.058) over the baseline model trained on the same number of European individuals. When we allowed iPGS to use n = 284,661 individuals, we observed an average improvement of 60.8% for African, 11.6% for South Asian, 7.3% for non-British White, 4.8% for White British, and 17.8% for the other individuals. We further developed iPGS+refit to jointly model the ancestry-shared and -dependent genetic effects when heterogeneous genetic associations were present. For neutrophil count, for example, iPGS+refit showed the highest predictive performance in the African group (R2 = 0.115), which exceeds the best predictive performance for the White British group (R2 = 0.090 in the iPGS model), even though only 1.49% of individuals used in the iPGS training are of African ancestry. Our results indicate the power of including diverse individuals for developing more equitable PGS models.