Search

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 data sets remains a major obstacle. Today, more realistic simulations are possible thanks to large increases in the quantity and quality of available genetic data, and to 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 frameworkstdpopsimseeks to lower this barrier by facilitating the simulation of complex population genetic models using up-to-date information. The initial version ofstdpopsimfocused 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 ofstdpopsim(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 three-fold and broadened coverage across the tree of life. During the process of expanding the catalog, we have identified common sticking points and developed 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 tostdpopsimaim 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.

Linkage disequilibrium (LD) is the correlation among nearby genetic variants. In genetic association studies, LD is often modeled using massive local correlation matrices, but this approach is slow, especially in ancestrally diverse studies. Here, 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 publish LDGMs and ancestry specific LDGM precision matrices for 18 million common SNPs (MAF>1%) in five ancestry groups, validate their accuracy, and demonstrate order-of-magnitude improvements in runtime for commonly used LD matrix computations. We implement an extremely fast multi-ancestry 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 genetic association data across millions of variants and individuals.

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.

Stochastic simulation is a key tool in population genetics, since the models involved are often analytically intractable and simulation is usually the only way of obtaining ground-truth data to evaluate inferences. Because of this necessity, a large number of specialised simulation programs have been developed, each filling a particular niche, but with largely overlapping functionality and a substantial duplication of effort. Here, we introducemsprimeversion 1.0, which efficiently implements ancestry and mutation simulations based on the succinct tree sequence data structure andtskitlibrary. We summarisemsprime’s many features, and show that its performance is excellent, often many times faster and more memory efficient than specialised alternatives. These high-performance features have been thoroughly tested and validated, and built using a collaborative, open source development model, which reduces duplication of effort and promotes software quality via community engagement.

Inferring the full genealogical history of a set of DNA sequences is a core problem in evolutionary biology, because this history encodes information about the events and forces that have influenced a species. However, current methods are limited, and the most accurate techniques are able to process no more than a hundred samples. As datasets that consist of millions of genomes are now being collected, there is a need for scalable and efficient inference methods to fully utilize these resources. Here we introduce an algorithm that is able to not only infer whole-genome histories with comparable accuracy to the state-of-the-art but also process four orders of magnitude more sequences. The approach also provides an 'evolutionary encoding' of the data, enabling efficient calculation of relevant statistics. We apply the method to human data from the 1000 Genomes Project, Simons Genome Diversity Project and UK Biobank, showing that the inferred genealogies are rich in biological signal and efficient to process.

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. ispartof: AMERICAN JOURNAL OF PHYSICAL ANTHROPOLOGY vol:168 issue:2 pages:340-351 ispartof: location:United States status: published

A central problem in evolutionary biology is to infer the full genealogical history of a set of DNA sequences. This history contains rich information about the forces that have influenced a sexually reproducing species. However, existing methods are limited: the most accurate is unable to cope with more than a few dozen samples. With modern genetic data sets rapidly approaching millions of genomes, there is an urgent need for efficient inference methods to exploit such rich resources. We introduce an algorithm to infer whole-genome history which has comparable accuracy to the state-of-the-art but can process around four orders of magnitude more sequences. Additionally, our method results in an “evolutionary encoding” of the original sequence data, enabling efficient access to genealogies and calculation of genetic statistics over the data. We apply this technique to human data from the 1000 Genomes Project, Simons Genome Diversity Project and UK Biobank, showing that the genealogies we estimate are both rich in biological signal and efficient to process.