Extremely sparse models of linkage disequilibrium in ancestrally diverse association studies

Linkage disequilibrium (LD) is the correlation among nearby genetic variants. In genetic association studies, LD is often modeled using large correlation matrices, but this approach is inefficient, especially in ancestrally diverse studies. In the present study, we introduce LD graphical models (LDGMs), which are an extremely sparse and efficient representation of LD. LDGMs are derived from genome-wide genealogies; statistical relationships among alleles in the LDGM correspond to genealogical relationships among haplotypes. We published LDGMs and ancestry-specific LDGM precision matrices for 18 million common variants (minor allele frequency >1%) in five ancestry groups, validated their accuracy and demonstrated order-of-magnitude improvements in runtime for commonly used LD matrix computations. We implemented an extremely fast multiancestry polygenic prediction method, BLUPx-ldgm, which performs better than a similar method based on the reference LD correlation matrix. LDGMs will enable sophisticated methods that scale to ancestrally diverse genetic association data across millions of variants and individuals.

Expanding the stdpopsim species catalog, and lessons learned for realistic genome simulations

Simulation is a key tool in population genetics for both methods development and empirical research, but producing simulations that recapitulate the main features of genomic datasets remains a major obstacle. Today, more realistic simulations are possible thanks to large increases in the quantity and quality of available genetic data, and the sophistication of inference and simulation software. However, implementing these simulations still requires substantial time and specialized knowledge. These challenges are especially pronounced for simulating genomes for species that are not well-studied, since it is not always clear what information is required to produce simulations with a level of realism sufficient to confidently answer a given question. The community-developed framework stdpopsim seeks to lower this barrier by facilitating the simulation of complex population genetic models using up-to-date information. The initial version of stdpopsim focused on establishing this framework using six well-characterized model species (Adrion et al., 2020). Here, we report on major improvements made in the new release of stdpopsim (version 0.2), which includes a significant expansion of the species catalog and substantial additions to simulation capabilities. Features added to improve the realism of the simulated genomes include non-crossover recombination and provision of species-specific genomic annotations. Through community-driven efforts, we expanded the number of species in the catalog more than threefold and broadened coverage across the tree of life. During the process of expanding the catalog, we have identified common sticking points and developed the best practices for setting up genome-scale simulations. We describe the input data required for generating a realistic simulation, suggest good practices for obtaining the relevant information from the literature, and discuss common pitfalls and major considerations. These improvements to stdpopsim aim to further promote the use of realistic whole-genome population genetic simulations, especially in non-model organisms, making them available, transparent, and accessible to everyone.

A unified genealogy of modern and ancient genomes

The sequencing of modern and ancient genomes from around the world has revolutionized our understanding of human history and evolution. However, the problem of how best to characterize ancestral relationships from the totality of human genomic variation remains unsolved. Here, we address this challenge with nonparametric methods that enable us to infer a unified genealogy of modern and ancient humans. This compact representation of multiple datasets explores the challenges of missing and erroneous data and uses ancient samples to constrain and date relationships. We demonstrate the power of the method to recover relationships between individuals and populations as well as to identify descendants of ancient samples. Finally, we introduce a simple nonparametric estimator of the geographical location of ancestors that recapitulates key events in human history.

East Anglian early Neolithic monument burial linked to contemporary Megaliths

In the fourth millennium BCE a cultural phenomenon of monumental burial structures spread along the Atlantic façade. Megalithic burials have been targeted for aDNA analyses, but a gap remains in East Anglia, where Neolithic structures were generally earthen or timber. An early Neolithic (3762–3648 cal. BCE) burial monument at the site of Trumpington Meadows, Cambridgeshire, UK, contained the partially articulated remains of at least three individuals. To determine whether this monument fits a pattern present in megalithic burials regarding sex bias, kinship, diet and relationship to modern populations, teeth and ribs were analysed for DNA and carbon and nitrogen isotopic values, respectively. Whole ancient genomes were sequenced from two individuals to a mean genomic coverage of 1.6 and 1.2X and genotypes imputed. Results show that they were brothers from a small population genetically and isotopically similar to previously published British Neolithic individuals, with a level of genome-wide homozygosity consistent with a small island population sourced from continental Europe, but bearing no signs of recent inbreeding. The first Neolithic whole genomes from a monumental burial in East Anglia confirm that this region was connected with the larger pattern of Neolithic megaliths in the British Isles and the Atlantic façade.

Evaluating macroscopic sex estimation methods using genetically sexed archaeological material: The medieval skeletal collection from St John's Divinity School, Cambridge

OBJECTIVES

In tests on known individuals macroscopic sex estimation has between 70% and 98% accuracy. However, materials used to create and test these methods are overwhelming modern. As sexual dimorphism is dependent on multiple factors, it is unclear whether macroscopic methods have similar success on earlier materials, which differ in lifestyle and nutrition. This research aims to assess the accuracy of commonly used traits by comparing macroscopic sex estimates to genetic sex in medieval English material.

MATERIALS AND METHODS

Sixty-six individuals from the 13th to 16th century Hospital of St John the Evangelist, Cambridge, were assessed. Genetic sex was determined using a shotgun approach. Eighteen skeletal traits were examined, and macroscopic sex estimates were derived from the os coxae, skull, and os coxae and skull combined. Each trait was tested for accuracy to explore sex estimates errors.

RESULTS

The combined estimate (97.7%) outperformed the os coxae only estimate (95.7%), which outperformed the skull only estimate (90.4%). Accuracy rates for individual traits varied: Phenice traits were most accurate, whereas supraorbital margins, frontal bossing, and gonial flaring were least accurate. The preauricular sulcus and arc compose showed a bias in accuracy between sexes.

DISCUSSION

Macroscopic sex estimates are accurate when applied to medieval material from Cambridge. However, low trait accuracy rates may relate to differences in dimorphism between the method derivative sample and the St John's collection. Given the sex bias, the preauricular sulcus, frontal bossing, and arc compose should be reconsidered as appropriate traits for sex estimation for this group.