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

Inferring human phenotypes using ancient DNA: from molecules to populations

The increasing availability of ancient DNA (aDNA) from human groups across space and time has yielded deep insights into the movements of our species. However, given the challenges of mapping from genotype to phenotype, aDNA has revealed less about the phenotypes of ancient individuals. In this review, we highlight recent advances in inferring ancient phenotypes - from the molecular to population scale - with a focus on applications enabled by new machine learning approaches. The genetic architecture of complex traits across human groups suggests that the prediction of individual-level complex traits, like behavior or disease risk, is often challenging across the relevant evolutionary distances. Thus, we propose an approach that integrates predictions of molecular phenotypes, whose mechanisms are more conserved, with nongenetic data.

Social and genetic diversity in first farmers of central Europe

The Linearbandkeramik (LBK) Neolithic communities were the first to spread farming across large parts of Europe. We report genome-wide data for 250 individuals: 178 individuals from whole-cemetery surveys of the Alföld Linearbankeramik Culture eastern LBK site of Polgár-Ferenci-hát, the western LBK site of Nitra Horné Krškany and the western LBK settlement and massacre site of Asparn-Schletz, as well as 48 LBK individuals from 16 other sites and 24 earlier Körös and Starčevo individuals from 17 more sites. Here we show a systematically higher percentage of western hunter-gatherer ancestry in eastern than in western LBK sites, showing that these two distinct LBK groups had different genetic trajectories. We find evidence for patrilocality, with more structure across sites in the male than in the female lines and a higher rate of within-site relatives for males. At Asparn-Schletz we find almost no relatives, showing that the massacred individuals were from a large population, not a small community.

Genetic Polymorphism of Y-Chromosome in Turkmen Population from Turkmenistan

This study investigates the Y-chromosome genetic diversity of the Turkmen population in Turkmenistan, analyzing 23 Y-STR loci for the first time in a sample of 100 individuals. Combined with comparative data from Turkmen populations in Afghanistan, Iran, Iraq, Russia, and Uzbekistan, this analysis offers insights into the genetic structure and relationships among Turkmen populations across regions across Central Asia and the Near East. High haplotype diversity in the Turkmen of Turkmenistan is shaped by founder effects (lineage expansions) from distinct haplogroups, with haplogroups Q and R1a predominating. Subhaplogroups Q1a and Q1b identified in Turkmenistan trace back to ancient Y-chromosome lineages from the Bronze Age. Comparative analyses, including genetic distance (RST), median-joining network, and multidimensional scaling (MDS), highlight the genetic proximity of the Turkmen in Turkmenistan to those in Afghanistan and Iran, while Iraqi Turkmen display unique characteristics, aligning with Near Eastern populations. This study underscores the Central Asian genetic affinity across most Turkmen populations. It demonstrates the value of deep-sequencing Y-chromosome data in tracing the patrilineal history of Central Asia for future studies. These findings contribute to a more comprehensive understanding of Turkmen genetic ancestry and add new data to the ongoing study of Central Asian population genetics.

Unveiling ancestral threads: Exploring CCR5 ∆32 mutation frequencies in Colombian populations for HIV/AIDS therapeutics

AIDS remains a significant global health challenge since its emergence in 1981, with millions of deaths and new cases every year. The CCR5 ∆32 genetic deletion confers immunity to HIV infection by altering a cell membrane protein crucial for viral entry. Stem cell transplants from homozygous carriers of this mutation to HIV-infected individuals have resulted in viral load reduction and disease remission, suggesting a potential therapeutic avenue. This study aims to investigate the relationship between genetic ancestry and the frequency of the CCR5 ∆32 mutation in Colombian populations, exploring the feasibility of targeted donor searches based on ancestry composition. Utilizing genomic data from the CÓDIGO-Colombia consortium, comprising 532 individuals, the study assessed the presence of the CCR5 ∆32 mutation and examined if the population was on Hardy-Weinberg equilibrium. Individuals were stratified into clusters based on African, American, and European ancestry percentages, with logistic regression analysis performed to evaluate the association between ancestry and mutation frequency. Additionally, global genomic databases were utilized to visualize the worldwide distribution of the mutation. The findings revealed a significant positive association between European ancestry and the CCR5 ∆32 mutation frequency, underscoring its relevance in donor selection. African and American ancestry showed negative but non-significant associations with CCR5 ∆32 frequency, which may be attributed to the study's limitations. These results emphasize the potential importance of considering ancestry in donor selection strategies, reveal the scarcity of potential donors in Colombia, and underscore the need to consider donors from other populations with mainly European ancestry if the CCR5 ∆32 stem cell transplant becomes a routine treatment for HIV/AIDS in Colombia.

Pervasive findings of directional selection realize the promise of ancient DNA to elucidate human adaptation

We present a method for detecting evidence of natural selection in ancient DNA time-series data that leverages an opportunity not utilized in previous scans: testing for a consistent trend in allele frequency change over time. By applying this to 8433 West Eurasians who lived over the past 14000 years and 6510 contemporary people, we find an order of magnitude more genome-wide significant signals than previous studies: 347 independent loci with >99% probability of selection. Previous work showed that classic hard sweeps driving advantageous mutations to fixation have been rare over the broad span of human evolution, but in the last ten millennia, many hundreds of alleles have been affected by strong directional selection. Discoveries include an increase from ~0% to ~20% in 4000 years for the major risk factor for celiac disease at HLA-DQB1; a rise from ~0% to ~8% in 6000 years of blood type B; and fluctuating selection at the TYK2 tuberculosis risk allele rising from ~2% to ~9% from ~5500 to ~3000 years ago before dropping to ~3%. We identify instances of coordinated selection on alleles affecting the same trait, with the polygenic score today predictive of body fat percentage decreasing by around a standard deviation over ten millennia, consistent with the "Thrifty Gene" hypothesis that a genetic predisposition to store energy during food scarcity became disadvantageous after farming. We also identify selection for combinations of alleles that are today associated with lighter skin color, lower risk for schizophrenia and bipolar disease, slower health decline, and increased measures related to cognitive performance (scores on intelligence tests, household income, and years of schooling). These traits are measured in modern industrialized societies, so what phenotypes were adaptive in the past is unclear. We estimate selection coefficients at 9.9 million variants, enabling study of how Darwinian forces couple to allelic effects and shape the genetic architecture of complex traits.

Improving data archiving practices in ancient genomics

Ancient DNA is producing a rich record of past genetic diversity in humans and other species. However, unless the primary data is appropriately archived, its long-term value will not be fully realised. I surveyed publicly archived data from 42 recent ancient genomics studies. Half of the studies archived incomplete datasets, preventing accurate replication and representing a loss of data of potential future use. No studies met all criteria that could be considered best practice. Based on these results, I make six recommendations for data producers: (1) archive all sequencing reads, not just those that aligned to a reference genome, (2) archive read alignments too, but as secondary analysis files, (3) provide correct experiment metadata on samples, libraries and sequencing runs, (4) provide informative sample metadata, (5) archive data from low-coverage and negative experiments, and (6) document archiving choices in papers, and peer review these. Given the reliance on destructive sampling of finite material, ancient genomics studies have a particularly strong responsibility to ensure the longevity and reusability of generated data.

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.

Genetic analysis and natural history of Parkinson's disease due to the LRRK2 G2019S variant

The LRRK2 G2019S variant is the most common cause of monogenic Parkinson's disease (PD); however, questions remain regarding the penetrance, clinical phenotype and natural history of carriers. We performed a 3.5-year prospective longitudinal online study in a large number of 1286 genotyped LRRK2 G2019S carriers and 109 154 controls, with and without PD, recruited from the 23andMe Research Cohort. We collected self-reported motor and non-motor symptoms every 6 months, as well as demographics, family histories and environmental risk factors. Incident cases of PD (phenoconverters) were identified at follow-up. We determined lifetime risk of PD using accelerated failure time modelling and explored the impact of polygenic risk on penetrance. We also computed the genetic ancestry of all LRRK2 G2019S carriers in the 23andMe database and identified regions of the world where carrier frequencies are highest. We observed that despite a 1 year longer disease duration (P = 0.016), LRRK2 G2019S carriers with PD had similar burden of motor symptoms, yet significantly fewer non-motor symptoms including cognitive difficulties, REM sleep behaviour disorder (RBD) and hyposmia (all P-values ≤ 0.0002). The cumulative incidence of PD in G2019S carriers by age 80 was 49%. G2019S carriers had a 10-fold risk of developing PD versus non-carriers. This rose to a 27-fold risk in G2019S carriers with a PD polygenic risk score in the top 25% versus non-carriers in the bottom 25%. In addition to identifying ancient founding events in people of North African and Ashkenazi descent, our genetic ancestry analyses infer that the G2019S variant was later introduced to Spanish colonial territories in the Americas. Our results suggest LRRK2 G2019S PD appears to be a slowly progressive predominantly motor subtype of PD with a lower prevalence of hyposmia, RBD and cognitive impairment. This suggests that the current prodromal criteria, which are based on idiopathic PD, may lack sensitivity to detect the early phases of LRRK2 PD in G2019S carriers. We show that polygenic burden may contribute to the development of PD in the LRRK2 G2019S carrier population. Collectively, the results should help support screening programmes and candidate enrichment strategies for upcoming trials of LRRK2 inhibitors in early-stage disease.

The Allen Ancient DNA Resource (AADR) a curated compendium of ancient human genomes

More than two hundred papers have reported genome-wide data from ancient humans. While the raw data for the vast majority are fully publicly available testifying to the commitment of the paleogenomics community to open data, formats for both raw data and meta-data differ. There is thus a need for uniform curation and a centralized, version-controlled compendium that researchers can download, analyze, and reference. Since 2019, we have been maintaining the Allen Ancient DNA Resource (AADR), which aims to provide an up-to-date, curated version of the world's published ancient human DNA data, represented at more than a million single nucleotide polymorphisms (SNPs) at which almost all ancient individuals have been assayed. The AADR has gone through six public releases at the time of writing and review of this manuscript, and crossed the threshold of >10,000 individuals with published genome-wide ancient DNA data at the end of 2022. This note is intended as a citable descriptor of the AADR.

An empirical test of the estimation of historical effective population size using Drosophila melanogaster

The availability of a large number of high-density markers (SNPs) allows the estimation of historical effective population size (Ne ) from linkage disequilibrium between loci. A recent refinement of methods to estimate historical Ne from the recent past has been shown to be rather accurate with simulation data. The method has also been applied to real data for numerous species. However, the simulation data cannot encompass all the complexities of real genomes, and the performance of any estimation method with real data is always uncertain, as the true demography of the populations is not known. Here, we carried out an experimental design with Drosophila melanogaster to test the method with real data following a known demographic history. We used a population maintained in the laboratory with a constant census size of about 2800 individuals and subjected the population to a drastic decline to a size of 100 individuals. After a few generations, the population was expanded back to the previous size and after a few further generations again expanded to twice the initial size. Estimates of historical Ne were obtained with the software GONE both for autosomal and X chromosomes from samples of 17 individuals sequenced for the whole genome. Estimates of the historical effective size were able to infer the patterns of changes that occurred in the populations showing generally good performance of the method. We discuss the limitations of the method and the application of the software carried out so far.

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.

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.

Familial dysautonomia

Familial dysautonomia (FD) is an autosomal recessive hereditary sensory and autonomic neuropathy (HSAN, type 3) expressed at birth with profound sensory loss and early death. The FD founder mutation in the ELP1 gene arose within the Ashkenazi Jews in the sixteenth century and is present in 1:30 Jews of European ancestry. The mutation yield a tissue-specific skipping of exon 20 and a loss of function of the elongator-1 protein (ELP1), which is essential for the development and survival of neurons. Patients with FD produce variable amounts of ELP1 in different tissues, with the brain producing mostly mutant transcripts. Patients have excessive blood pressure variability due to the failure of the IXth and Xth cranial nerves to carry baroreceptor signals. Neurogenic dysphagia causes frequent aspiration leading to chronic pulmonary disease. Characteristic hyperadrenergic "autonomic crises" consisting of brisk episodes of severe hypertension, tachycardia, skin blotching, retching, and vomiting occur in all patients. Progressive features of the disease include retinal nerve fiber loss and blindness, and proprioceptive ataxia with severe gait impairment. Chemoreflex failure may explain the high frequency of sudden death in sleep. Although 99.5% of patients are homozygous for the founder mutation, phenotypic severity varies, suggesting that modifier genes impact expression. Medical management is currently symptomatic and preventive. Disease-modifying therapies are close to clinical testing. Endpoints to measure efficacy have been developed, and the ELP1 levels are a good surrogate endpoint for target engagement. Early intervention may be critical for treatment to be successful.