The genome-wide structure of the Jewish people

Enriching iPSC research diversity: Harnessing human biobank collections for improved ethnic representation

Biobanks of human biosamples and cell lines are indispensable for biomedical research on human health and disease and for drug development projects. Many human cell line biobanks worldwide hold collections of lymphoblastoid cell lines (LCLs), representing thousands of affected and control donors from diverse ethnic/ancestry groups. In recent years, induced human pluripotent stem cells (iPSCs) and differentiated human cells derived from these iPSCs have become indispensable for applied biomedical research. Establishing iPSCs remains a laborious and costly step towards generating differentiated human cells. To address this research need, several non-profit and commercial biobanks have established iPSC collections for distribution to researchers, thereby serving as a resource for generating differentiated human cells. The most common starting materials for generation of iPSCs are a skin biopsy for harvesting fibroblasts, or a blood sample for collection of peripheral blood mononuclear cells. However untapped resources include the large established collections of biobanked human LCLs which can be reprogrammed to iPSCs using a variety of published protocols including the use of non-integrating episomal vectors. Many biobanks curate LCLs from diverse ethnic/ancestry populations, an aspect largely absent in most established iPSC biobanks which tend to primarily reflect populations from developed countries. Here, we call upon researchers across the breadth of iPSC research to tap the unique resource of existing and diverse human LCL collections for establishing biobanked iPSC panels that better represent the varied human ethnic (and hence genomic) diversity, thereby benefiting precision medicine and drug development research on a global scale.

Genome wide association study and genomic risk prediction of age related macular degeneration in Israel

The risk of developing age-related macular degeneration (AMD) is influenced by genetic background. In 2016, the International AMD Genomics Consortium (IAMDGC) identified 52 risk variants in 34 loci, and a polygenic risk score (PRS) from these variants was associated with AMD. The Israeli population has a unique genetic composition: Ashkenazi Jewish (AJ), Jewish non-Ashkenazi, and Arab sub-populations. We aimed to perform a genome-wide association study (GWAS) for AMD in Israel, and to evaluate PRSs for AMD. Our discovery set recruited 403 AMD patients and 256 controls at Hadassah Medical Center. We genotyped individuals via custom exome chip. We imputed non-typed variants using cosmopolitan and AJ reference panels. We recruited additional 155 cases and 69 controls for validation. To evaluate predictive power of PRSs for AMD, we used IAMDGC summary-statistics excluding our study and developed PRSs via clumping/thresholding or LDpred2. In our discovery set, 31/34 loci reported by IAMDGC were AMD-associated (P < 0.05). Of those, all effects were directionally consistent with IAMDGC and 11 loci had a P-value under Bonferroni-corrected threshold (0.05/34 = 0.0015). At a 5 × 10-5 threshold, we discovered four suggestive associations in FAM189A1, IGDCC4, C7orf50, and CNTNAP4. Only the FAM189A1 variant was AMD-associated in the replication cohort after Bonferroni-correction. A prediction model including LDpred2-based PRS + covariates had an AUC of 0.82 (95% CI 0.79-0.85) and performed better than covariates-only model (P = 5.1 × 10-9). Therefore, previously reported AMD-associated loci were nominally associated with AMD in Israel. A PRS developed based on a large international study is predictive in Israeli populations.

Evidence for dynastic succession among early Celtic elites in Central Europe

The early Iron Age (800 to 450 BCE) in France, Germany and Switzerland, known as the 'West-Hallstattkreis', stands out as featuring the earliest evidence for supra-regional organization north of the Alps. Often referred to as 'early Celtic', suggesting tentative connections to later cultural phenomena, its societal and population structure remain enigmatic. Here we present genomic and isotope data from 31 individuals from this context in southern Germany, dating between 616 and 200 BCE. We identify multiple biologically related groups spanning three elite burials as far as 100 km apart, supported by trans-regional individual mobility inferred from isotope data. These include a close biological relationship between two of the richest burial mounds of the Hallstatt culture. Bayesian modelling points to an avuncular relationship between the two individuals, which may suggest a practice of matrilineal dynastic succession in early Celtic elites. We show that their ancestry is shared on a broad geographic scale from Iberia throughout Central-Eastern Europe, undergoing a decline after the late Iron Age (450 BCE to ~50 CE).

Smoking changes adaptive immunity with persistent effects

Individuals differ widely in their immune responses, with age, sex and genetic factors having major roles in this inherent variability1-6. However, the variables that drive such differences in cytokine secretion-a crucial component of the host response to immune challenges-remain poorly defined. Here we investigated 136 variables and identified smoking, cytomegalovirus latent infection and body mass index as major contributors to variability in cytokine response, with effects of comparable magnitudes with age, sex and genetics. We find that smoking influences both innate and adaptive immune responses. Notably, its effect on innate responses is quickly lost after smoking cessation and is specifically associated with plasma levels of CEACAM6, whereas its effect on adaptive responses persists long after individuals quit smoking and is associated with epigenetic memory. This is supported by the association of the past smoking effect on cytokine responses with DNA methylation at specific signal trans-activators and regulators of metabolism. Our findings identify three novel variables associated with cytokine secretion variability and reveal roles for smoking in the short- and long-term regulation of immune responses. These results have potential clinical implications for the risk of developing infections, cancers or autoimmune diseases.

Body mass index and visual impairment in Israeli adolescents: A nationwide study

BACKGROUND

Previous research on the association between body mass index (BMI) and visual impairment (VI) in youth has reported inconsistent findings. We aimed to investigate this association in a national cohort of Israeli adolescents.

METHODS

This retrospective, population-based, cross-sectional study included 1 697 060 adolescents (56.4% men; mean age 17 years) who underwent mandatory pre-military service assessments from 1993 to 2017. BMI was classified based on the US age- and sex-matched percentiles. Unilateral or bilateral VI was classified as best-corrected visual acuity (BCVA) worse than 6/9 in either or both eyes, respectively. Sex-stratified regression models adjusted for sociodemographic variables were used to analyse the BMI-VI relationship.

RESULTS

Overall, 17 871 (1.05%) and 5148 (0.30%) adolescents had unilateral and bilateral VI, respectively. Compared with high-normal BMI (50th to 85th percentile), adjusted odds ratios (ORs) for unilateral and bilateral VI gradually increased with higher BMI, reaching 1.33 (1.13-1.55) and 1.80 (1.37-2.35) in men with severe obesity, and 1.51 (1.24-1.84) and 1.52 (1.08-2.14) in women with severe obesity, respectively. Men with underweight also had increased ORs for unilateral and bilateral VI (1.23; 1.14-1.33 and 1.59; 1.37-1.84, respectively), a pattern not observed in women (0.96; 0.86-1.07 and 1.02; 0.83-1.25, respectively). Results were maintained when the outcome was restricted to mild VI, as well as in subgroups of adolescents with unimpaired health and those without moderate-to-severe myopia.

CONCLUSIONS

Abnormal BMI, and particularly obesity, is associated with increased OR for VI in late adolescence.

Genetic Affinities and Adaptation of the South-West Coast Populations of India

Evolutionary event has not only altered the genetic structure of human populations but also associated with social and cultural transformation. South Asian populations were the result of migration and admixture of genetically and culturally diverse groups. Most of the genetic studies pointed to large-scale admixture events between Ancestral North Indian (ANI) and Ancestral South Indian (ASI) groups, also additional layers of recent admixture. In the present study, we have analyzed 213 individuals inhabited in South-west coast India with traditional warriors and feudal lord status and historically associated with migratory events from North/North West India and possible admixture with West Eurasian populations, whose genetic links are still missing. Analysis of autosomal Single Nucleotide Polymorphism (SNP) markers suggests that these groups possibly derived their ancestry from some groups of North West India having additional Middle Eastern genetic components. Higher distribution of West Eurasian mitochondrial haplogroups also points to female-mediated admixture. Estimation of Effective Migration Surface (EEMS) analysis indicates Central India and Godavari basin as a crucial transition zone for population migration from North and North West India to South-west coastal India. Selection screen using 3 distinct outlier-based approaches revealed genetic signatures related to Immunity and protection from Viral infections. Thus, our study suggests that the South-west coastal groups with traditional warriors and feudal lords' status are of a distinct lineage compared to Dravidian and Gangetic plain Indo-Europeans and are remnants of very early migrations from North West India following the Godavari basin to Karnataka and Kerala.

Evaluation of European-based polygenic risk score for breast cancer in Ashkenazi Jewish women in Israel

BACKGROUND

Polygenic risk score (PRS), calculated based on genome-wide association studies (GWASs), can improve breast cancer (BC) risk assessment. To date, most BC GWASs have been performed in individuals of European (EUR) ancestry, and the generalisation of EUR-based PRS to other populations is a major challenge. In this study, we examined the performance of EUR-based BC PRS models in Ashkenazi Jewish (AJ) women.

METHODS

We generated PRSs based on data on EUR women from the Breast Cancer Association Consortium (BCAC). We tested the performance of the PRSs in a cohort of 2161 AJ women from Israel (1437 cases and 724 controls) from BCAC (BCAC cohort from Israel (BCAC-IL)). In addition, we tested the performance of these EUR-based BC PRSs, as well as the established 313-SNP EUR BC PRS, in an independent cohort of 181 AJ women from Hadassah Medical Center (HMC) in Israel.

RESULTS

In the BCAC-IL cohort, the highest OR per 1 SD was 1.56 (±0.09). The OR for AJ women at the top 10% of the PRS distribution compared with the middle quintile was 2.10 (±0.24). In the HMC cohort, the OR per 1 SD of the EUR-based PRS that performed best in the BCAC-IL cohort was 1.58±0.27. The OR per 1 SD of the commonly used 313-SNP BC PRS was 1.64 (±0.28).

CONCLUSIONS

Extant EUR GWAS data can be used for generating PRSs that identify AJ women with markedly elevated risk of BC and therefore hold promise for improving BC risk assessment in AJ women.

Genome-wide association study and genomic risk prediction of age-related macular degeneration in Israel

Purpose

The risk of developing age-related macular degeneration(AMD) is influenced by genetic background. In 2016, International AMD Genomics Consortium(IAMDGC) identified 52 risk variants in 34 loci, and a polygenic risk score(PRS) based on these variants was associated with AMD. The Israeli population has a unique genetic composition: Ashkenazi Jewish(AJ), Jewish non-Ashkenazi, and Arab sub-populations. We aimed to perform a genome-wide association study(GWAS) for AMD in Israel, and to evaluate PRSs for AMD.

Methods

For our discovery set, we recruited 403 AMD patients and 256 controls at Hadassah Medical Center. We genotyped all individuals via custom exome chip. We imputed non-typed variants using cosmopolitan and AJ reference panels. We recruited additional 155 cases and 69 controls for validation. To evaluate predictive power of PRSs for AMD, we used IAMDGC summary statistics excluding our study and developed PRSs via either clumping/thresholding or LDpred2.

Results

In our discovery set, 31/34 loci previously reported by the IAMDGC were AMD associated with P<0.05. Of those, all effects were directionally consistent with the IAMDGC and 11 loci had a p-value under Bonferroni-corrected threshold(0.05/34=0.0015). At a threshold of 5x10 -5 , we discovered four suggestive associations in FAM189A1 , IGDCC4 , C7orf50 , and CNTNAP4 . However, only the FAM189A1 variant was AMD associated in the replication cohort after Bonferroni-correction. A prediction model including LDpred2-based PRS and other covariates had an AUC of 0.82(95%CI:0.79-0.85) and performed better than a covariates-only model(P=5.1x10 -9 ).

Conclusions

Previously reported AMD-associated loci were nominally associated with AMD in Israel. A PRS developed based on a large international study is predictive in Israeli populations.

The Qatari population's genetic structure and gene flow as revealed by the Y chromosome

The Y-chromosome has been widely used in forensic genetic applications and human population genetic studies due to its uniparental origins. A large database on the Qatari population was created for comparison with other databases from the Arabian Peninsula, the Middle East, and Africa. We provide a study of 23 Y-STR loci included in PowerPlex Y23 (Promega, USA) that were genotyped to produce haplotypes in 379 unrelated males from Qatar, a country at the crossroads of migration patterns. Overall, the most polymorphic locus provided by the Promega kit was DYS458, with a genetic diversity value of 0.85 and a haplotype diversity of 0.998924. Athey's Haplogroup Predictor tool was used to predict haplogroups from Y-STR haplotypes in the Qatari population. In a median-joining network, the haplogroup J1 predominance (49%) in Qatar generated a star-like expansion cluster. The graph of population Q-matrix was developed using Y-STR data from 38 Middle Eastern and 97 African populations (11,305 individuals), and it demonstrated a stronger sub-grouping of countries within each ethnic group and showed the effect of Arabs on the indigenous Berbers of North Africa. The estimated migration rate between the Qatari and other Arabian populations was inferred using Bayesian coalescence theory in the Migrate-n program. According to the Gene Flow study, the main migration route was from Yemen to Kuwait through Qatar. Our research, using the PowerPlex Y23 database, shows the importance of gene diversity, as well as regional and social structuring, in determining the utility of demographic and forensic databases.

A pathogenic variant in the uncharacterized RNF212B gene results in severe aneuploidy male infertility and repeated IVF failure

Quantitative and qualitative spermatogenic impairments are major causes of men's infertility. Although in vitro fertilization (IVF) is effective, some couples persistently fail to conceive. To identify causal variants in patients with severe male infertility factor and repeated IVF failures, we sequenced the exome of two consanguineous family members who underwent several failed IVF cycles and were diagnosed with low sperm count and motility. We identified a rare homozygous nonsense mutation in a previously uncharacterized gene, RNF212B, as the causative variant. Recurrence was identified in another unrelated, infertile patient who also faced repeated failed IVF treatments. scRNA-seq demonstrated meiosis-specific expression of RNF212B. Sequence analysis located a protein domain known to be associated with aneuploidy, which can explain multiple IVF failures. Accordingly, FISH analysis revealed a high aneuploidy rate in the patients' sperm cells and their IVF embryos. Finally, inactivation of the Drosophila orthologs significantly reduced male fertility. Given that members of the evolutionary conserved RNF212 gene family are involved in meiotic recombination and crossover maturation, our findings indicate a critical role of RNF212B in meiosis, genome stability, and in human fertility. Since recombination is completely absent in Drosophila males, our findings may indicate an additional unrelated role for the RNF212-like paralogs in spermatogenesis.

Glucose intolerance in pregnancy and risk of early-onset type 2 diabetes: a population-based cohort study

BACKGROUND

The risk of type 2 diabetes among women with glucose intolerance during pregnancy that does not meet gestational diabetes criteria requires further investigation. We aimed to explore the associations between various degrees of gestational glucose intolerance and the risk of type 2 diabetes in young adulthood.

METHODS

For this population-based cohort study, the national Israeli conscription database was linked to Maccabi Healthcare Services (MHS), the second-largest state-mandated health provider in Israel. We included 177 241 women who underwent a pre-recruitment evaluation at adolescence (age 16-20 years), 1 year before mandatory military service, and later underwent, from Jan 1, 2001, to Dec 31, 2019, two-step gestational diabetes screening with a 50 g glucose challenge test (GCT) based on a threshold of 140 mg/dL (7·8 mmol/L), followed as needed by a 100 g oral glucose tolerance test (OGTT). Abnormal OGTT values were defined according to the Carpenter-Coustan thresholds: 95 mg/dL (5·3 mmol/L) or higher in the fasting state; 180 mg/dL (10·0 mmol/L) or higher at 1 h; 155 mg/dL (8·6 mmol/L) or higher at 2 h; and 140 mg/dL (7·8 mmol/L) or higher at 3 h. The primary outcome was incident type 2 diabetes in the MHS diabetes registry. Cox proportional hazards models were applied to estimate adjusted hazard ratios (HRs) with 95% CIs for incident type 2 diabetes.

FINDINGS

During a cumulative follow-up of 1 882 647 person-years, and with a median follow-up of 10·8 (IQR 5·2-16·4) years, 1262 women were diagnosed with type 2 diabetes. Crude incidence rates of type 2 diabetes were 2·6 (95% CI 2·4-2·9) per 10 000 person-years in women with gestational normoglycaemia, 8·9 (7·4-10·6) per 10 000 person-years in women with an abnormal GCT and normal OGTT, 26·1 (22·4-30·1) per 10 000 person-years in women with one abnormal OGTT value (in the fasting state or 1 h, 2 h, or 3 h post-challenge), and 71·9 (66·0-78·3) per 10 000 person-years in women with gestational diabetes. After adjustment for sociodemographic characteristics, adolescent BMI, and age at gestational screening, the risk of type 2 diabetes was higher, compared to the gestational normoglycaemia group, in women with an abnormal GCT and normal OGTT (adjusted hazard ratio [HR] 3·39 [95% CI 2·77-4·16]; p<0·0001), in women with one abnormal OGTT value (9·11 [7·64-10·86]; p<0·0001), and in women with gestational diabetes (24·84 [21·78-28·34]; p<0·0001). The risk of type 2 diabetes was modestly increased in women with isolated elevated fasting glucose (adjusted HR 11·81 [95% CI 8·58-16·25]; p<0·0001), and in women with gestational diabetes and an abnormal fasting glucose (38·02 [32·41-44·61]; p<0·0001).

INTERPRETATION

Gestational glucose intolerance, including conditions not meeting gestational diabetes criteria of the two-step strategy, confers a high risk of type 2 diabetes in young adulthood. These conditions should be recognised as risk factors for type 2 diabetes, especially among women with abnormal fasting glucose concentrations during pregnancy.

FUNDING

None.

Genomic analyses of hair from Ludwig van Beethoven

Ludwig van Beethoven (1770-1827) remains among the most influential and popular classical music composers. Health problems significantly impacted his career as a composer and pianist, including progressive hearing loss, recurring gastrointestinal complaints, and liver disease. In 1802, Beethoven requested that following his death, his disease be described and made public. Medical biographers have since proposed numerous hypotheses, including many substantially heritable conditions. Here we attempt a genomic analysis of Beethoven in order to elucidate potential underlying genetic and infectious causes of his illnesses. We incorporated improvements in ancient DNA methods into existing protocols for ancient hair samples, enabling the sequencing of high-coverage genomes from small quantities of historical hair. We analyzed eight independently sourced locks of hair attributed to Beethoven, five of which originated from a single European male. We deemed these matching samples to be almost certainly authentic and sequenced Beethoven's genome to 24-fold genomic coverage. Although we could not identify a genetic explanation for Beethoven's hearing disorder or gastrointestinal problems, we found that Beethoven had a genetic predisposition for liver disease. Metagenomic analyses revealed furthermore that Beethoven had a hepatitis B infection during at least the months prior to his death. Together with the genetic predisposition and his broadly accepted alcohol consumption, these present plausible explanations for Beethoven's severe liver disease, which culminated in his death. Unexpectedly, an analysis of Y chromosomes sequenced from five living members of the Van Beethoven patrilineage revealed the occurrence of an extra-pair paternity event in Ludwig van Beethoven's patrilineal ancestry.

Exploring regional aspects of 3D facial variation within European individuals

Facial ancestry can be described as variation that exists in facial features that are shared amongst members of a population due to environmental and genetic effects. Even within Europe, faces vary among subregions and may lead to confounding in genetic association studies if unaccounted for. Genetic studies use genetic principal components (PCs) to describe facial ancestry to circumvent this issue. Yet the phenotypic effect of these genetic PCs on the face has yet to be described, and phenotype-based alternatives compared. In anthropological studies, consensus faces are utilized as they depict a phenotypic, not genetic, ancestry effect. In this study, we explored the effects of regional differences on facial ancestry in 744 Europeans using genetic and anthropological approaches. Both showed similar ancestry effects between subgroups, localized mainly to the forehead, nose, and chin. Consensus faces explained the variation seen in only the first three genetic PCs, differing more in magnitude than shape change. Here we show only minor differences between the two methods and discuss a combined approach as a possible alternative for facial scan correction that is less cohort dependent, more replicable, non-linear, and can be made open access for use across research groups, enhancing future studies in this field.

The genetic history of Scandinavia from the Roman Iron Age to the present

We investigate a 2,000-year genetic transect through Scandinavia spanning the Iron Age to the present, based on 48 new and 249 published ancient genomes and genotypes from 16,638 modern individuals. We find regional variation in the timing and magnitude of gene flow from three sources: the eastern Baltic, the British-Irish Isles, and southern Europe. British-Irish ancestry was widespread in Scandinavia from the Viking period, whereas eastern Baltic ancestry is more localized to Gotland and central Sweden. In some regions, a drop in current levels of external ancestry suggests that ancient immigrants contributed proportionately less to the modern Scandinavian gene pool than indicated by the ancestry of genomes from the Viking and Medieval periods. Finally, we show that a north-south genetic cline that characterizes modern Scandinavians is mainly due to the differential levels of Uralic ancestry and that this cline existed in the Viking Age and possibly earlier.

Genome-wide data from medieval German Jews show that the Ashkenazi founder event pre-dated the 14th century

We report genome-wide data from 33 Ashkenazi Jews (AJ), dated to the 14^th century, obtained following a salvage excavation at the medieval Jewish cemetery of Erfurt, Germany. The Erfurt individuals are genetically similar to modern AJ, but they show more variability in Eastern European-related ancestry than modern AJ. A third of the Erfurt individuals carried a mitochondrial lineage common in modern AJ and eight carried pathogenic variants known to affect AJ today. These observations, together with high levels of runs of homozygosity, suggest that the Erfurt community had already experienced the major reduction in size that affected modern AJ. The Erfurt bottleneck was more severe, implying substructure in medieval AJ. Overall, our results suggest that the AJ founder event and the acquisition of the main sources of ancestry pre-dated the 14^th century and highlight late medieval genetic heterogeneity no longer present in modern AJ.

Genomes from a medieval mass burial show Ashkenazi-associated hereditary diseases pre-date the 12th century

We report genome sequence data from six individuals excavated from the base of a medieval well at a site in Norwich, UK. A revised radiocarbon analysis of the assemblage is consistent with these individuals being part of a historically attested episode of antisemitic violence on 6 February 1190 CE. We find that four of these individuals were closely related and all six have strong genetic affinities with modern Ashkenazi Jews. We identify four alleles associated with genetic disease in Ashkenazi Jewish populations and infer variation in pigmentation traits, including the presence of red hair. Simulations indicate that Ashkenazi-associated genetic disease alleles were already at appreciable frequencies, centuries earlier than previously hypothesized. These findings provide new insights into a significant historical crime, into Ashkenazi population history, and into the origins of genetic diseases associated with modern Jewish populations.

The Anglo-Saxon migration and the formation of the early English gene pool

The history of the British Isles and Ireland is characterized by multiple periods of major cultural change, including the influential transformation after the end of Roman rule, which precipitated shifts in language, settlement patterns and material culture1. The extent to which migration from continental Europe mediated these transitions is a matter of long-standing debate2-4. Here we study genome-wide ancient DNA from 460 medieval northwestern Europeans-including 278 individuals from England-alongside archaeological data, to infer contemporary population dynamics. We identify a substantial increase of continental northern European ancestry in early medieval England, which is closely related to the early medieval and present-day inhabitants of Germany and Denmark, implying large-scale substantial migration across the North Sea into Britain during the Early Middle Ages. As a result, the individuals who we analysed from eastern England derived up to 76% of their ancestry from the continental North Sea zone, albeit with substantial regional variation and heterogeneity within sites. We show that women with immigrant ancestry were more often furnished with grave goods than women with local ancestry, whereas men with weapons were as likely not to be of immigrant ancestry. A comparison with present-day Britain indicates that subsequent demographic events reduced the fraction of continental northern European ancestry while introducing further ancestry components into the English gene pool, including substantial southwestern European ancestry most closely related to that seen in Iron Age France5,6.

Alzheimer’s Disease Polygenic Risk Score Is Not Associated With Cognitive Decline Among Older Adults With Type 2 Diabetes

Objectives

Multiple risk loci for late-onset Alzheimer’s disease (LOAD) have been identified. Type 2 diabetes (T2D) is a risk factor for cognitive decline, dementia and Alzheimer’s disease (AD). We investigated the association of polygenic risk score (PRS) for LOAD with overall cognitive functioning and longitudinal decline, among older adults with T2D.

Methods

The study included 1046 Jewish participants from the Israel Diabetes and Cognitive Decline (IDCD) study, aged ≥ 65 years, diagnosed with T2D, and cognitively normal at baseline. The PRS included variants from 26 LOAD associated loci (at genome-wide significance level), and was calculated with and without APOE. Outcome measures, assessed in 18 months intervals, were global cognition and the specific domains of episodic memory, attention/working memory, executive functions, and language/semantic categorization. Random coefficient models were used for analysis, adjusting for demographic variables, T2D-related characteristics, and cardiovascular factors. Additionally, in a subsample of 202 individuals, we analyzed the association of PRS with the volumes of total gray matter, frontal lobe, hippocampus, amygdala, and white matter hyperintensities. Last, the association of PRS with amyloid beta (Aβ) burden was examined in 44 participants who underwent an 18F-flutemetamol PET scan.

Results

The PRS was not significantly associated with overall functioning or decline in global cognition or any of the specific cognitive domains. Similarly, following correction for multiple testing, there was no association with Aβ burden and other brain imaging phenotypes.

Conclusion

Our results suggest that the cumulative effect of LOAD susceptibility loci is not associated with a greater rate of cognitive decline in older adults with T2D, and other pathways may underlie this link.

Principal Component Analyses (PCA)-based findings in population genetic studies are highly biased and must be reevaluated

Principal Component Analysis (PCA) is a multivariate analysis that reduces the complexity of datasets while preserving data covariance. The outcome can be visualized on colorful scatterplots, ideally with only a minimal loss of information. PCA applications, implemented in well-cited packages like EIGENSOFT and PLINK, are extensively used as the foremost analyses in population genetics and related fields (e.g., animal and plant or medical genetics). PCA outcomes are used to shape study design, identify, and characterize individuals and populations, and draw historical and ethnobiological conclusions on origins, evolution, dispersion, and relatedness. The replicability crisis in science has prompted us to evaluate whether PCA results are reliable, robust, and replicable. We analyzed twelve common test cases using an intuitive color-based model alongside human population data. We demonstrate that PCA results can be artifacts of the data and can be easily manipulated to generate desired outcomes. PCA adjustment also yielded unfavorable outcomes in association studies. PCA results may not be reliable, robust, or replicable as the field assumes. Our findings raise concerns about the validity of results reported in the population genetics literature and related fields that place a disproportionate reliance upon PCA outcomes and the insights derived from them. We conclude that PCA may have a biasing role in genetic investigations and that 32,000-216,000 genetic studies should be reevaluated. An alternative mixed-admixture population genetic model is discussed.

Glucose Intolerance in Pregnancy and Offspring Obesity in Late Adolescence

OBJECTIVE

Gestational hyperglycemia is associated with deleterious neonatal outcomes, but long-term risks for offspring obesity are less clear. We estimated the odds for offspring adolescent overweight and obesity among mothers with gestational glucose intolerance.

RESEARCH DESIGN AND METHODS

In a mother-offspring historical cohort, the Israel military conscription data set was linked to a large health maintenance organization. Included were women who were evaluated at adolescence and underwent two-step gestational diabetes screening (mean age, 31 years) with a 50-g glucose challenge test (GCT), followed by a 100-g oral glucose tolerance test (OGTT) if the result was abnormal. Glucose tolerance categories included gestational normoglycemia, abnormal GCT with normal OGTT, impaired glucose tolerance (IGT; one abnormal OGTT value), and gestational diabetes. The primary outcome was offspring overweight/obesity (BMI ≥85th percentile) at adolescence, measured prior to military conscription. Logistic regression models were applied.

RESULTS

Of 33,482 mother-offspring pairs, overweight and obesity were observed in 6,516 offspring. Across increasing categories of pregnancy glycemia, the proportions of offspring with adolescent overweight/obesity increased: normoglycemia, 19%; abnormal GCT with normal OGTT, 22%; gestational IGT, 24%; and gestational diabetes, 25% (P < 0.0001). Corresponding odds ratios after adjustment for the mother's late adolescent characteristics (sociodemographic confounders and BMI) and pregnancy age were 1.2 (95% CI 1.1-1.4), 1.3 (1.2-1.5), and 1.4 (1.3-1.6), respectively. Further adjustment for offspring birth weight percentile and sociodemographic variables did not materially change results. Associations were more pronounced with increasing obesity severity.

CONCLUSIONS

Gestational glucose intolerance, including categories not meeting the gestational diabetes threshold, was associated with increased odds for offspring overweight/obesity at late adolescence.

An efficient method to identify, date, and describe admixture events using haplotype information

We present fastGLOBETROTTER, an efficient new haplotype-based technique to identify, date, and describe admixture events using genome-wide autosomal data. With simulations, we show how fastGLOBETROTTER reduces computation time by an order of magnitude relative to the related technique GLOBETROTTER without suffering loss of accuracy. We apply fastGLOBETROTTER to a cohort of more than 6000 Europeans from 10 countries, revealing previously unreported admixture signals. In particular, we infer multiple periods of admixture related to East Asian or Siberian-like sources, starting >2000 yr ago, in people living in countries north of the Baltic Sea. In contrast, we infer admixture related to West Asian, North African, and/or Southern European sources in populations south of the Baltic Sea, including admixture dated to ∼300-700 CE, overlapping the fall of the Roman Empire, in people from Belgium, France, and parts of Germany. Our new approach scales to analyzing hundreds to thousands of individuals from a putatively admixed population and, hence, is applicable to emerging large-scale cohorts of genetically homogeneous populations.

Invasive Management in Older Adults (≥80 Years) With Non-ST Elevation Myocardial Infarction

OBJECTIVE

To evaluate the association of invasive management (coronary angiogram) with all-cause mortality among older adult (≥80 years of age) patients presenting with non-ST elevation myocardial infarction (NSTEMI) by frailty status.

PATIENTS AND METHODS

This study used a retrospective cohort of consecutive older adult patients who were hospitalized with NSTEMI as their primary clinical diagnosis between August 1, 2008, and December 31, 2019. Cox regression models were applied with stratification by frailty status (low, medium, and high). Extensive sensitivity analyses were conducted including propensity score matching and inverse probability treatment weighting models.

RESULTS

The study population included 2317 patients with median age of 86 years (IQR, 83-90 years) of whom 1243 (53.6%) were men. Patients who were managed invasively (n=581 [25%]) were less likely to be frail (7% vs 44%, P<.001). During the follow-up (median, 19 months, IQR, 4-41 months), 1599 (69%) patients died. In a multivariable Cox model, invasive approach was associated with adjusted hazard ratio (HR) of 0.61 (95% CI, 0.53 to 0.71) for the risk of death. The benefit of invasive approach was consistent among low, medium, and high frailty subgroups with adjusted HRs of 0.74 (95% CI, 0.58 to 0.93), 0.65 (95% CI, 0.50 to 0.85), and 0.52 (95% CI, 0.34 to 0.78), respectively (P for interaction = 0.48). Results were consistent with propensity score matching and inverse probability treatment weighting analyses (HR, 0.6; 95% CI, 0.50 to 0.71 and HR, 0.67; 95% CI, 0.55 to 0.82, respectively). Sensitivity analysis addressing potential immortal time bias and residual confounding yielded similar results.

CONCLUSION

Invasive approach is associated with improved survival among older adults with NSTEMI irrespective of frailty status.

Assessing temporal and geographic contacts across the Adriatic Sea through the analysis of genome-wide data from Southern Italy

Southern Italy was characterised by a complex prehistory that started with different Palaeolithic cultures, later followed by the Neolithization and the demic dispersal from the Pontic-Caspian Steppe during the Bronze Age. Archaeological and historical evidences point to a link between Southern Italians and the Balkans still present in modern times. To shed light on these dynamics, we analysed around 700 South Mediterranean genomes combined with informative ancient DNAs. Our findings revealed high affinities of South-Eastern Italians with modern Eastern Peloponnesians, and a closer affinity of ancient Greek genomes with those from specific regions of South Italy than modern Greek genomes. The higher similarity could be associated with a Bronze Age component ultimately originating from the Caucasus with high Iranian and Anatolian Neolithic ancestries. Furthermore, extremely differentiated allele frequencies among Northern and Southern Italy revealed putatively adapted SNPs in genes involved in alcohol metabolism, nevi features and immunological traits.

The spatiotemporal patterns of major human admixture events during the European Holocene

Recent studies have shown that admixture has been pervasive throughout human history. While several methods exist for dating admixture in contemporary populations, they are not suitable for sparse, low coverage ancient genomic data. Thus, we developed DATES (Distribution of Ancestry Tracts of Evolutionary Signals) that leverages ancestry covariance patterns across the genome of a single individual to infer the timing of admixture. DATES provides reliable estimates under various demographic scenarios and outperforms available methods for ancient DNA applications. Using DATES on~1100 ancient genomes from sixteen regions in Europe and west Asia, we reconstruct the chronology of the formation of the ancestral populations and the fine-scale details of the spread of Neolithic farming and Steppe pastoralist-related ancestry across Europe. By studying the genetic formation of Anatolian farmers, we infer that gene flow related to Iranian Neolithic farmers occurred before 9600 BCE, predating the advent of agriculture in Anatolia. Contrary to the archaeological evidence, we estimate that early Steppe pastoralist groups (Yamnaya and Afanasievo) were genetically formed more than a millennium before the start of Steppe pastoralism. Our analyses provide new insights on the origins and spread of farming and Indo-European languages, highlighting the power of genomic dating methods to elucidate the legacy of human migrations.

Genome-Wide Marker Data-Based Comparative Population Analysis of Szeklers From Korond, Transylvania, and From Transylvania Living Non-Szekler Hungarians

Genome-wide genotype data from 48 carefully selected population samples of Transylvania-living Szeklers and non-Szekler Hungarians were analyzed by comparative analysis. Our analyses involved contemporary Hungarians living in Hungary, other Europeans, and Eurasian samples counting 530 individuals altogether. The source of the Szekler samples was the commune of Korond, Transylvania. The analyzed non-Szekler Hungarian samples were collected from villages with a history dating back to the era of the Árpád Dynasty. Population structure by principal component analysis and ancestry analysis also revealed a great within-group similarity of the analyzed Szeklers and non-Szekler Transylvanian Hungarians. These groups also showed similar genetic patterns with each other. Haplotype analyses using identity-by-descent segment discovering tools showed that average pairwise identity-by-descent sharing is similar in the investigated populations, but the Korond Szekler samples had higher average sharing with the Hungarians from Hungary than non-Szekler Transylvanian Hungarians. Average sharing results showed that both groups are isolated compared to other Europeans, and pointed out that the non-Szekler Transylvanian Hungarian inhabitants of the investigated Árpád Age villages are more isolated than investigated Szeklers from Korond. This was confirmed by our autozygosity analysis as well. Identity-by-descent segment analyses and 4-population tests also confirmed that these Hungarian-speaking Transylvanian ethnic groups are strongly related to Hungarians living in Hungary.

Fine-Scale Genetic Structure in the United Arab Emirates Reflects Endogamous and Consanguineous Culture, Population History, and Geography

The indigenous population of the United Arab Emirates (UAE) has a unique demographic and cultural history. Its tradition of endogamy and consanguinity is expected to produce genetic homogeneity and partitioning of gene pools while population movements and intercontinental trade are likely to have contributed to genetic diversity. Emiratis and neighboring populations of the Middle East have been underrepresented in the population genetics literature with few studies covering the broader genetic history of the Arabian Peninsula. Here, we genotyped 1,198 individuals from the seven Emirates using 1.7 million markers and by employing haplotype-based algorithms and admixture analyses, we reveal the fine-scale genetic structure of the Emirati population. Shared ancestry and gene flow with neighboring populations display their unique geographic position while increased intra- versus inter-Emirati kinship and sharing of uniparental haplogroups, reflect the endogamous and consanguineous cultural traditions of the Emirates and their tribes.

Recurring pathogenic variants in the BRCA2 gene in the Ethiopian Jewish population. Founder mutations?

Mutations in the BRCA1 and BRCA2 genes increase the risk for various cancers including breast, ovarian, prostate, pancreas and melanoma. Identifying BRCA1/2 mutation carriers enables risk assessment, surveillance, early detection and risk reduction. In certain Israeli sub-populations recurring and founder mutations have been identified and for these, testing for founder mutations is simple, efficient and cost-effective. Founder mutations in the Jewish Ethiopian population have not been described. We report here the identification of a recurring BRCA2 mutation in the Ethiopian Jewish population; c.5159C>A; p.Ser1720Ter, which has only been described once before in this population. In addition, in another family of the same origin we found the BRCA2 c.7579delG; p.Val2527Ter mutation that has been previously described in two different Jewish Ethiopian families. In Israel genetic testing is performed in a sequential stepwise manner, first testing a panel of predominant mutations and if negative further testing by gene sequencing is offered. Recently it has been decided to expand the founder mutation panel to include mutations which have been found in two or more separate families. This new panel will include the BRCA2 c.7579delG; p.Val2527Ter mutation, and we recommend that the BRCA2 c.5159C>A; p.Ser1720Ter mutation should also be added to the new predominant mutation panel.

The opposing trends of body mass index and blood pressure during 1977-2020; nationwide registry of 2.8 million male and female adolescents

Background

Elevated blood pressure among adolescents has been shown to be associated with future adverse cardiovascular outcomes and early onset diabetes. Most data regarding systolic and diastolic blood pressure trends are based on surveys of selected populations within 10–20-year periods. The goal of this study was to characterize the secular trend of blood pressure given the rising prevalence of adolescent obesity.

Methods

This nationwide population-based study included 2,785,515 Israeli adolescents (41.6% females, mean age 17.4 years) who were medically evaluated and whose weight, height and blood pressure were measured, prior to mandatory military service between 1977 and 2020. The study period was divided into 5-year intervals. Linear regression models were used to describe the P for trend along the time intervals. Analysis of covariance was used to calculate means of blood pressure adjusted for body mass index.

Results

During the study period, the mean body mass index increased by 2.1 and 1.6 kg/m² in males and females, respectively (P for trend  < 0.001 in both sexes). The mean diastolic blood pressure decreased by 3.6 mmHg in males and by 2.9 mmHg in females (P  < 0.001 in both sexes). The mean systolic blood pressure increased by 1.6 mmHg in males and decreased by 1.9 mmHg in females. These trends were also consistent when blood pressure values were adjusted to body mass index.

Conclusion

Despite the increase in body mass index over the last four decades, diastolic blood pressure decreased in both sexes while systolic blood pressure increased slightly in males and decreased in females.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12933-021-01433-0.

A common founder effect of the splice site variant c.-23 + 1G > A in GJB2 gene causing autosomal recessive deafness 1A (DFNB1A) in Eurasia

Mutations in the GJB2 gene are known to be a major cause of autosomal recessive deafness 1A (OMIM 220290). The most common pathogenic variants of the GJB2 gene have a high ethno-geographic specificity in their distribution, being attributed to a founder effect related to the Neolithic migration routes of Homo sapiens. The c.-23 + 1G > A splice site variant is frequently found among deaf patients of both Caucasian and Asian origins. It is currently unknown whether the spread of this mutation across Eurasia is a result of the founder effect or if it could have multiple local centers of origin. To determine the origin of c.-23 + 1G > A, we reconstructed haplotypes by genotyping SNPs on an Illumina OmniExpress 730 K platform of 23 deaf individuals homozygous for this variant from different populations of Eurasia. The analyses revealed the presence of common regions of homozygosity in different individual genomes in the sample. These data support the hypothesis of the common founder effect in the distribution of the c.-23 + 1G > A variant of the GJB2 gene. Based on the published data on the c.-23 + 1G > A prevalence among 16,177 deaf people and the calculation of the TMRCA of the modified f2-haplotypes carrying this variant, we reconstructed the potential migration routes of the carriers of this mutation around the world. This analysis indicates that the c.-23 + 1G > A variant in the GJB2 gene may have originated approximately 6000 years ago in the territory of the Caucasus or the Middle East then spread throughout Europe, South and Central Asia and other regions of the world.

Asthma in Youth and Early-onset Type 2 Diabetes: A Nationwide Study of 1.72 Million Israeli Adolescents

BACKGROUND

The prevalence of both asthma and early-onset diabetes is on the rise; however, the association between them remains unclear. We examined a possible association of asthma at adolescence with type 2 diabetes in young adulthood.

METHODS

This is a nationwide, population-based study of 1 718 541 Israeli adolescents (57% males; mean age 17.3 years; range 16-19 years), examined before compulsory military service between 1992 and 2016, with data linked to the Israeli National Diabetes Registry. Asthma diagnosis and severity were determined by a board-certified pulmonologist and based on spirometry tests.

RESULTS

Type 2 diabetes developed in 58/9090 (0.64%), 507/97 059 (0.52%), 114/23 332 (0.49%), and 7095/1 589 060 (0.44%) persons with moderate-to-severe, mild, inactive, and no history of asthma, respectively, during a mean follow-up >13 years. The respective odds ratios (ORs) were 1.33 (95% CI, 1.02-1.74), 1.17 (1.06-1.28), and 1.09 (0.9-1.31), considering those without asthma history as the reference, in a model adjusted for birth year, sex, body mass index, and other sociodemographic variables. The association persisted when the analysis accounted for coexisting morbidities, and when persons without asthma, individually matched by age, sex, birth year, and body mass index were the reference. Both mild and moderate-to-severe asthma were associated with type 2 diabetes before age 35 years: ORs 1.18 (1.05-1.34) and 1.44 (1.05-2.00), respectively. The strength of the association was accentuated over time. The effect was unchanged when adjusted for oral and inhaled glucocorticoid use.

CONCLUSION

Adolescents with active asthma have higher risk to develop type 2 diabetes. This seems related to disease severity, independent of adolescent obesity status, apparent before age 35 years, and more pronounced in recent years.

Ancestral patterns of recessive dystrophic epidermolysis bullosa mutations in Hispanic populations suggest sephardic ancestry

Recessive dystrophic epidermolysis bullosa (RDEB) is a rare genodermatosis caused by mutations in the gene coding for type VII collagen (COL7A1). More than 800 different pathogenic mutations in COL7A1 have been described to date; however, the ancestral origins of many of these mutations have not been precisely identified. In this study, 32 RDEB patient samples from the Southwestern United States, Mexico, Chile, and Colombia carrying common mutations in the COL7A1 gene were investigated to determine the origins of these mutations and the extent to which shared ancestry contributes to disease prevalence. The results demonstrate both shared European and American origins of RDEB mutations in distinct populations in the Americas and suggest the influence of Sephardic ancestry in at least some RDEB mutations of European origins. Knowledge of ancestry and relatedness among RDEB patient populations will be crucial for the development of future clinical trials and the advancement of novel therapeutics.

Dissecting the genetic history of the Roman Catholic populations of West Coast India

Migration and admixture history of populations have always been curious and an interesting theme. The West Coast of India harbours a rich diversity, bestowing various ethno-linguistic groups, with many of them having well-documented history of migrations. The Roman Catholic is one such distinct group, whose origin was much debated. While some historians and anthropologists relating them to ancient group of Gaud Saraswat Brahmins, others relating them for being members of the Jews Lost Tribes in the first Century migration to India. Historical records suggests that this community was later forcibly converted to Christianity by the Portuguese in Goa during the Sixteenth Century. Till date, no genetic study was done on this group to infer their origin and genetic affinity. Hence, we analysed 110 Roman Catholics from three different locations of West Coast of India including Goa, Kumta and Mangalore using both uniparental and autosomal markers to understand their genetic history. We found that the Roman Catholics have close affinity with the Indo-European linguistic groups, particularly Brahmins. Additionally, we detected genetic signal of Jews in the linkage disequilibrium-based admixture analysis, which was absent in other Indo-European populations, who are inhabited in the same geographical regions. Haplotype-based analysis suggests that the Roman Catholics consist of South Asian-specific ancestry and showed high drift. Ancestry-specific historical population size estimation points to a possible bottleneck around the time of Goan inquisition (fifteenth century). Analysis of the Roman Catholics data along with ancient DNA data of Neolithic and bronze age revealed that the Roman Catholics fits well in a basic model of ancient ancestral composition, typical of most of the Indo-European caste groups of India. Mitochondrial DNA (mtDNA) analysis suggests that most of the Roman Catholics have aboriginal Indian maternal genetic ancestry; while the Y chromosomal DNA analysis indicates high frequency of R1a lineage, which is predominant in groups with higher ancestral North Indian (ANI) component. Therefore, we conclude that the Roman Catholics of Goa, Kumta and Mangalore regions are the remnants of very early lineages of Brahmin community of India, having Indo-Europeans genetic affinity along with cryptic Jewish admixture, which needs to be explored further.

Evidence of the interplay of genetics and culture in Ethiopia

The rich linguistic, ethnic and cultural diversity of Ethiopia provides an unprecedented opportunity to understand the level to which cultural factors correlate with-and shape-genetic structure in human populations. Using primarily new genetic variation data covering 1,214 Ethiopians representing 68 different ethnic groups, together with information on individuals' birthplaces, linguistic/religious practices and 31 cultural practices, we disentangle the effects of geographic distance, elevation, and social factors on the genetic structure of Ethiopians today. We provide evidence of associations between social behaviours and genetic differences among present-day peoples. We show that genetic similarity is broadly associated with linguistic affiliation, but also identify pronounced genetic similarity among groups from disparate language classifications that may in part be attributable to recent intermixing. We also illustrate how groups reporting the same culture traits are more genetically similar on average and show evidence of recent intermixing, suggesting that shared cultural traits may promote admixture. In addition to providing insights into the genetic structure and history of Ethiopia, we identify the most important cultural and geographic predictors of genetic differentiation and provide a resource for designing sampling protocols for future genetic studies involving Ethiopians.

Integrating Linguistics, Social Structure, and Geography to Model Genetic Diversity within India

India represents an intricate tapestry of population substructure shaped by geography, language, culture, and social stratification. Although geography closely correlates with genetic structure in other parts of the world, the strict endogamy imposed by the Indian caste system and the large number of spoken languages add further levels of complexity to understand Indian population structure. To date, no study has attempted to model and evaluate how these factors have interacted to shape the patterns of genetic diversity within India. We merged all publicly available data from the Indian subcontinent into a data set of 891 individuals from 90 well-defined groups. Bringing together geography, genetics, and demographic factors, we developed Correlation Optimization of Genetics and Geodemographics to build a model that explains the observed population genetic substructure. We show that shared language along with social structure have been the most powerful forces in creating paths of gene flow in the subcontinent. Furthermore, we discover the ethnic groups that best capture the diverse genetic substructure using a ridge leverage score statistic. Integrating data from India with a data set of additional 1,323 individuals from 50 Eurasian populations, we find that Indo-European and Dravidian speakers of India show shared genetic drift with Europeans, whereas the Tibeto-Burman speaking tribal groups have maximum shared genetic drift with East Asians.

The forensic landscape and the population genetic analyses of Hainan Li based on massively parallel sequencing DNA profiling

Due to the formation of the Qiongzhou Strait by climate change and marine transition, Hainan island was isolated from the mainland southern China during the Last Glacial Maximum. Hainan island, located at the southernmost part of China and separated from the Leizhou Peninsula by the Qiongzhou Strait, laid on one of the modern human northward migration routes from Southeast Asia to East Asia. The Hlai language-speaking Li minority, the second largest population after Han Chinese in Hainan island, is the direct descendants of the initial migrants in Hainan island and has unique ethnic properties and derived characteristics; however, the forensic-associated studies on Hainan Li population are still insufficient. Hence, 136 Hainan Li individuals were genotyped in this study using the MPS-based ForenSeq™ DNA Signature Prep Kit (DNA Primer Set A, DPMA) to characterize the forensic genetic polymorphism landscape, and DNA profiles were obtained from 152 different molecular genetic markers (27 autosomal STRs, 24 Y-STRs, 7 X-STRs, and 94 iiSNPs). A total of 419 distinct length variants and 586 repeat sequence sub-variants, with 31 novel alleles (at 17 loci), were identified across the 58 STR loci from the DNA profiles of Hainan Li population. We evaluated the forensic characteristics and efficiencies of DPMA, demonstrating that the STRs and iiSNPs in DPMA were highly polymorphic in Hainan Li population and could be employed in forensic applications. In addition, we set up three datasets, which included the genetic data of (i) iiSNPs (27 populations, 2640 individuals), (ii) Y-STRs (42 populations, 8281 individuals), and (iii) Y haplogroups (123 populations, 4837 individuals) along with the population ancestries and language families, to perform population genetic analyses separately from different perspectives. In conclusion, the phylogenetic analyses indicated that Hainan Li, with a southern East Asia origin and Tai-Kadai language-speaking language, is an isolated population relatively. But the genetic pool of Hainan Li influenced by the limited gene flows from other Tai-Kadai populations and Hainan populations. Furthermore, the establishment of isolated population models will be beneficial to clarify the exquisite population structures and develop specific genetic markers for subpopulations in forensic genetic fields.

European genetic ancestry associated with risk of childhood ependymoma

BACKGROUND

Ependymoma is a histologically defined central nervous system tumor most commonly occurring in childhood. Population-level incidence differences by race/ethnicity are observed, with individuals of European ancestry at highest risk. We aimed to determine whether extent of European genetic ancestry is associated with ependymoma risk in US populations.

METHODS

In a multi-ethnic study of Californian children (327 cases, 1970 controls), we estimated the proportions of European, African, and Native American ancestry among recently admixed Hispanic and African American subjects and estimated European admixture among non-Hispanic white subjects using genome-wide data. We tested whether genome-wide ancestry differences were associated with ependymoma risk and performed admixture mapping to identify associations with local ancestry. We also evaluated race/ethnicity-stratified ependymoma incidence data from the Central Brain Tumor Registry of the United States (CBTRUS).

RESULTS

CBTRUS data revealed that African American and Native American children have 33% and 36%, respectively, reduced incidence of ependymoma compared with non-Hispanic whites. In genetic analyses, a 20% increase in European ancestry was associated with a 1.31-fold higher odds of ependymoma among self-reported Hispanics and African Americans (95% CI: 1.08-1.59, Pmeta = 6.7 × 10-3). Additionally, eastern European ancestral substructure was associated with increased ependymoma risk in non-Hispanic whites (P = 0.030) and in Hispanics (P = 0.043). Admixture mapping revealed a peak at 20p13 associated with increased local European ancestry, and targeted fine-mapping identified a lead variant at rs6039499 near RSPO4 (odds ratio = 1.99; 95% CI: 1.45-2.73; P = 2.2 × 10-5) but which was not validated in an independent set of posterior fossa type A patients.

CONCLUSIONS

Interethnic differences in ependymoma risk are recapitulated in the genomic ancestry of ependymoma patients, implicating regions to target in future association studies.

Dissecting dynamics and differences of selective pressures in the evolution of human pigmentation

Human pigmentation is a highly diverse and complex trait among populations and has drawn particular attention from both academic and non-academic investigators for thousands of years. Previous studies detected selection signals in several human pigmentation genes, but few studies have integrated contribution from multiple genes to the evolution of human pigmentation. Moreover, none has quantified selective pressures on human pigmentation over epochs and between populations. Here, we dissect dynamics and differences of selective pressures during different periods and between distinct populations with new approaches. We use genotype data of 19 genes associated with human pigmentation from 17 publicly available datasets and obtain data for 2346 individuals of six representative population groups from across the world. Our results quantify the strength of natural selection on light pigmentation not only in modern Europeans (0.0259/generation) but also in proto-Eurasians (0.00650/generation). Our results also suggest that several derived alleles associated with human dark pigmentation may be under positive directional selection in some African populations. Our study provides the first attempt to quantitatively investigate the dynamics of selective pressures during different time periods in the evolution of human pigmentation.This article has an associated First Person interview with the first author of the article.

Adolescent Nonalcoholic Fatty Liver Disease and Type 2 Diabetes in Young Adulthood

CONTEXT

The long-term risk of type 2 diabetes in adolescents with nonalcoholic fatty liver disease (NAFLD) is unclear.

OBJECTIVE

To assess type 2 diabetes risk among adolescents with NAFLD.

DESIGN AND SETTING

A nationwide, population-based study of Israeli adolescents who were examined before military service during 1997-2011 and were followed until December 31, 2016.

PARTICIPANTS

A total of 1 025 796 normoglycemic adolescents were included.

INTERVENTIONS

Biopsy or radiographic tests were prerequisite for NAFLD diagnosis. Data were linked to the Israeli National Diabetes Registry.

MAIN OUTCOME MEASURES

Type 2 diabetes incidence.

RESULTS

During a mean follow-up of 13.3 years, 12 of 633 adolescents with NAFLD (1.9%; all with high body mass index [BMI] at baseline) were diagnosed with type 2 diabetes compared with 2917 (0.3%) adolescents without NAFLD. The hazard ratio (HR) for type 2 diabetes was 2.59 (95% confidence interval [CI], 1.47-4.58) for the NAFLD vs. the non-NAFLD group after adjustment for BMI and sociodemographic confounders. The elevated risk persisted in several sensitivity analyses. These included an analysis of persons without other metabolic comorbidities (adjusted HR, 2.75 [95% CI, 1.48-5.14]) and of persons with high BMI; and an analysis whose outcome was type 2 diabetes by age 30 years (adjusted HR, 2.14 [95% CI, 1.02-4.52]). The results remained significant when a sex-, birth year-, and BMI-matched control group was the reference (adjusted HR, 2.98 [95% CI, 1.54-5.74]).

CONCLUSIONS

Among normoglycemic adolescents, NAFLD was associated with an increased adjusted risk for type 2 diabetes, which may be apparent before age 30 years.

Cartographic atlas of frequency variation for 45 pharmacogenetic markers in populations of Russia and its neighbor states

The lack of information about the frequency of pharmacogenetic markers in Russia impedes the adoption of personalized treatment algorithms originally developed for West European populations. The aim of this paper was to study the distribution of some clinically significant pharmacogenetic markers across Russia. A total of 45 pharmacogenetic markers were selected from a few population genetic datasets, including ADME, drug target and hemostasis-controlling genes. The total number of donors genotyped for these markers was 2,197. The frequencies of these markers were determined for 50 different populations, comprised of 137 ethnic and subethnic groups. A comprehensive pharmacogenetic atlas was created, i.e. a systematic collection of gene geographic maps of frequency variation for 45 pharmacogenetic DNA markers in Russia and its neighbor states. The maps revealed 3 patterns of geographic variation. Clinal variation (a gradient change in frequency along the East-West axis) is observed in the pharmacogenetic markers that follow the main pattern of variation for North Eurasia (13% of the maps). Uniform distribution singles out a group of markers that occur at average frequency in most Russian regions (27% of the maps). Focal variation is observed in the markers that are specific to a certain group of populations and are absent in other regions (60% of the maps). The atlas reveals that the average frequency of the marker and its frequency in individual populations do not indicate the type of its distribution in Russia: a gene geographic map is needed to uncover the pattern of its variation.

Middle eastern genetic legacy in the paternal and maternal gene pools of Chuetas

Chuetas are a group of descendants of Majorcan Crypto-Jews (Balearic Islands, Spain) who were socially stigmatized and segregated by their Majorcan neighbours until recently; generating a community that, although after the seventeenth century no longer contained Judaic religious elements, maintained strong group cohesion, Jewishness consciousness, and endogamy. Collective memory fixed 15 surnames as a most important defining element of Chueta families. Previous studies demonstrated Chuetas were a differentiated population, with a considerable proportion of their original genetic make-up. Genetic data of Y-chromosome polymorphism and mtDNA control region showed, in Chuetas’ paternal lineages, high prevalence of haplogroups J2-M172 (33%) and J1-M267 (18%). In maternal lineages, the Chuetas hallmark is the presence of a new sub-branching of the rare haplogroup R0a2m as their modal haplogroup (21%). Genetic diversity in both Y-chromosome and mtDNA indicates the Chueta community has managed to avoid the expected heterogeneity decrease in their gene pool after centuries of isolation and inbreeding. Moreover, the composition of their uniparentally transmitted lineages demonstrates a remarkable signature of Middle Eastern ancestry—despite some degree of host admixture—confirming Chuetas have retained over the centuries a considerable degree of ancestral genetic signature along with the cultural memory of their Jewish origin.

Janus-faced race: Is race biological, social, or mythical?

As belief in the reality of race as a biological category among U.S. anthropologists has fallen, belief in the reality of race as a social category has risen in its place. The view that race simply does not exist-that it is a myth-is treated with suspicion. While racial classification is linked to many of the worst evils of recent history, it is now widely believed to be necessary to fight back against racism. In this article, I argue that race is indeed a biological fiction, but I critique the claim that race is socially real. I defend a form of anti-realist reconstructionism about race, which says that there are no races, only racialized groups-groups mistakenly believed to be races. I argue that this is the most attractive position about race from a metaphysical perspective, and that it is also the position most conductive to public understanding and social justice.

Investigating the genetic characteristics of the Csangos, a traditionally Hungarian speaking ethnic group residing in Romania

Csango people are an East-Central European ethnographic group living mostly in the historical region of Moldavia, Romania. Their traditional language, the Csango is an old Hungarian dialect, which is a severely endangered language due to language shift. Their origin is still disputed among experts and there are many hypotheses since the 19th century. Previous genetic studies found connection with ethnic groups living in Hungary and provided evidence which might support their Hungarian origin. Another study found Inner Asian Altaic ancestry in their genetic makeup. The goal of this study was to analyze the genetic characteristics of the Csango people by comparing their genetic characteristics to contemporary Eurasian populations based on genome-wide autosomal marker data. Our findings suggest that genetic affinity of Csangos to Hungarians is more significant than to Romanians. They also have a detectable connection with Central-Asian and Siberian Turkic ethnic groups. Besides the presumable Middle Eastern/Central-Asian Turkic ancestry, Csangos show ~4% Turkic ancestry from Central Asia/Siberia, which makes them unique in comparison to all other East-Central European populations investigated in this study. The admixture that resulted in this Turkic ancestry could have occurred 30-40 generations ago, which date interval corresponds to Hungarian historical events regarding their migration and the conquest of the Carpathian basin.

Recent Demographic History Inferred by High-Resolution Analysis of Linkage Disequilibrium

Inferring changes in effective population size (Ne) in the recent past is of special interest for conservation of endangered species and for human history research. Current methods for estimating the very recent historical Ne are unable to detect complex demographic trajectories involving multiple episodes of bottlenecks, drops, and expansions. We develop a theoretical and computational framework to infer the demographic history of a population within the past 100 generations from the observed spectrum of linkage disequilibrium (LD) of pairs of loci over a wide range of recombination rates in a sample of contemporary individuals. The cumulative contributions of all of the previous generations to the observed LD are included in our model, and a genetic algorithm is used to search for the sequence of historical Ne values that best explains the observed LD spectrum. The method can be applied from large samples to samples of fewer than ten individuals using a variety of genotyping and DNA sequencing data: haploid, diploid with phased or unphased genotypes and pseudohaploid data from low-coverage sequencing. The method was tested by computer simulation for sensitivity to genotyping errors, temporal heterogeneity of samples, population admixture, and structural division into subpopulations, showing high tolerance to deviations from the assumptions of the model. Computer simulations also show that the proposed method outperforms other leading approaches when the inference concerns recent timeframes. Analysis of data from a variety of human and animal populations gave results in agreement with previous estimations by other methods or with records of historical events.

Adolescent Obesity and Early-Onset Type 2 Diabetes

OBJECTIVE

Type 2 diabetes (T2D) is increasingly diagnosed at younger ages. We investigated the association of adolescent obesity with incident T2D at early adulthood.

RESEARCH DESIGN AND METHODS

A nationwide, population-based study evaluated 1,462,362 adolescents (59% men, mean age 17.4 years) during 1996-2016. Data were linked to the Israeli National Diabetes Registry. Weight and height were measured at study entry. Cox proportional models were applied.

RESULTS

During 15,810,751 person-years, 2,177 people (69% men) developed T2D (mean age at diagnosis 27 years). There was an interaction among BMI, sex, and incident T2D (P _interaction = 0.023). In a model adjusted for sociodemographic variables, the hazard ratios for diabetes diagnosis were 1.7 (95% CI 1.4-2.0), 2.8 (2.3-3.5), 5.8 (4.9-6.9), 13.4 (11.5-15.7), and 25.8 (21.0-31.6) among men in the 50th-74th percentile, 75th-84th percentile, overweight, mild obesity, and severe obesity groups, respectively, and 2.2 (1.6-2.9), 3.4 (2.5-4.6), 10.6 (8.3-13.6), 21.1 (16.0-27.8), and 44.7 (32.4-61.5), respectively, in women. An inverse graded relationship was observed between baseline BMI and mean age of T2D diagnosis: 27.8 and 25.9 years among men and women with severe obesity, respectively, and 29.5 and 28.5 years among low-normal BMI (5th-49th percentile; reference), respectively. The projected fractions of adult-onset T2D that were attributed to high BMI (≥85th percentile) at adolescence were 56.9% (53.8-59.9%) and 61.1% (56.8-65.2%) in men and women, respectively.

CONCLUSIONS

Severe obesity significantly increases the risk for incidence of T2D in early adulthood in both sexes. The rise in adolescent severe obesity is likely to increase diabetes incidence in young adults in coming decades.

Populations dynamics in Northern Eurasian forests: a long-term perspective from Northeast Asia

The 'Northern Eurasian Greenbelt' (NEG) is the northern forest zone stretching from the Japanese Archipelago to Northern Europe. The NEG has created highly productive biomes for humanity to exploit since the end of the Pleistocene. This research explores how the ecological conditions in northern Eurasia contributed to and affected human migrations and cultural trajectories by synthesizing the complimentary viewpoints of environmental archaeology, Geographic Information Science (GIS), genetics and linguistics. First, the environmental archaeology perspective raises the possibility that the NEG functioned as a vessel fostering people to develop diverse cultures and engage in extensive cross-cultural exchanges. Second, geographical analysis of genomic data on mitochondrial DNA using GIS reveals the high probability that population dynamics in the southeastern NEG promoted the peopling of the Americas at the end of the Pleistocene. Finally, a linguistic examination of environmental- and landscape-related vocabulary of the proto-Turkic language groups enables the outline of their original cultural landscape and natural conditions, demonstrating significant cultural spheres, i.e. from southern Siberia to eastern Inner Mongolia during Neolithization. All of these results combine to suggest that the ecological complex in the southern edge of the NEG in northeast Asia played a significant role in peopling across the continents during prehistory.

Association of a Variant in VWA3A with Response to Anti-Vascular Endothelial Growth Factor Treatment in Neovascular AMD

Purpose

Anti–vascular endothelial growth factor (VEGF) therapy for neovascular AMD (nvAMD) obtains a variable outcome. We performed a genome-wide association study for anti-VEGF treatment response in nvAMD to identify variants potentially underlying such a variable outcome.

Methods

Israeli patients with nvAMD who underwent anti-VEGF treatment (n = 187) were genotyped on a whole exome chip containing approximately 500,000 variants. Genotyping was correlated with delta visual acuity (deltaVA) between baseline and after three injections of anti-VEGF. Top principal components, age, and baseline VA were included in the analysis. Two lead associated variants were genotyped in an independent validation set of patients with nvAMD (n = 108).

Results

Linear regression analysis on 5,353,842 variants revealed five exonic variants with an association P value of less than 6 × 10⁻⁵. The top variant in the gene VWA3A (P = 1.77 × 10⁻⁶) was tested in the validation cohort. The minor allele of the VWA3A variant was associated with worse response to treatment (P = 0.02). The average deltaVA of discovery plus validation was –0.214 logMAR (≈ a gain of 10.7 Early Treatment Diabetic Retinopathy Study letters) for homozygote for the major allele, 0.172 logMAR for heterozygotes (≈ a loss of 8.6 Early Treatment Diabetic Retinopathy Study letters), and 0.21 logMAR for homozygote for the minor allele (≈ a loss of 10.5 Early Treatment Diabetic Retinopathy Study letters). Minor allele carriers had a higher frequency of macular hemorrhage at baseline.

Conclusions

An VWA3A gene variant was associated with worse response to anti-VEGF treatment in Israeli patients with nvAMD. The VWA3A protein is a precursor of the multimeric von Willebrand factor which is involved in blood coagulation, a system previously associated with nvAMD.

Ancient Rome: A genetic crossroads of Europe and the Mediterranean

Ancient Rome was the capital of an empire of ~70 million inhabitants, but little is known about the genetics of ancient Romans. Here we present 127 genomes from 29 archaeological sites in and around Rome, spanning the past 12,000 years. We observe two major prehistoric ancestry transitions: one with the introduction of farming and another prior to the Iron Age. By the founding of Rome, the genetic composition of the region approximated that of modern Mediterranean populations. During the Imperial period, Rome's population received net immigration from the Near East, followed by an increase in genetic contributions from Europe. These ancestry shifts mirrored the geopolitical affiliations of Rome and were accompanied by marked interindividual diversity, reflecting gene flow from across the Mediterranean, Europe, and North Africa.

Loss-of-Function Variants in SERPINA12 Underlie Autosomal Recessive Palmoplantar Keratoderma

Inherited palmoplantar keratodermas refer to a large and heterogeneous group of conditions resulting from abnormal epidermal differentiation and featuring thickening of the skin of the palms and soles. Here, we aimed at delineating the genetic basis of an autosomal recessive form of palmoplantar keratodermas manifesting with erythematous hyperkeratotic plaques over the palms and soles, extending to non-palmoplantar areas. Whole-exome sequencing in affected individuals revealed homozygous nonsense variants in the SERPINA12 gene. SERPINA12 encodes the visceral adipose tissue-derived serpin A12, a serine protease inhibitor. The pathogenic variants were found to result in reduced visceral adipose tissue-derived serpin A12 expression in patients' skin biopsies in comparison to healthy controls. In addition, SERPINA12 downregulation in three-dimensional skin equivalents was associated with marked epidermal acanthosis and hyperkeratosis, replicating the human phenotype. Moreover, decreased SERPINA12 expression resulted in reduced visceral adipose tissue-derived serpin A12-mediated inhibition of kallikrein 7 activity as well as decreased levels of desmoglein-1 and corneodesmosin, two known kallikrein 7 substrates, which are required for normal epidermal differentiation. The present data, taken collectively, demarcate a unique type of autosomal recessive palmoplantar keratodermas, attribute to visceral adipose tissue-derived serpin A12 a role in skin biology, and emphasize the importance of mechanisms regulating proteolytic activity for normal epidermal differentiation.

Genetic Landscapes Reveal How Human Genetic Diversity Aligns with Geography

Geographic patterns in human genetic diversity carry footprints of population history and provide insights for genetic medicine and its application across human populations. Summarizing and visually representing these patterns of diversity has been a persistent goal for human geneticists, and has revealed that genetic differentiation is frequently correlated with geographic distance. However, most analytical methods to represent population structure do not incorporate geography directly, and it must be considered post hoc alongside a visual summary of the genetic structure. Here, we estimate "effective migration" surfaces to visualize how human genetic diversity is geographically structured. The results reveal local patterns of differentiation in detail and emphasize that while genetic similarity generally decays with geographic distance, the relationship is often subtly distorted. Overall, the visualizations provide a new perspective on genetics and geography in humans and insight to the geographic distribution of human genetic variation.