The role of genetic selection and climatic factors in the dispersal of anatomically modern humans out of Africa.

The evolutionarily recent dispersal of anatomically modern humans (AMH) out of Africa (OoA) and across Eurasia provides a unique opportunity to examine the impacts of genetic selection as humans adapted to multiple new environments. Analysis of ancient Eurasian genomic datasets (~1,000 to 45,000 y old) reveals signatures of strong selection, including at least 57 hard sweeps after the initial AMH movement OoA, which have been obscured in modern populations by extensive admixture during the Holocene. The spatiotemporal patterns of these hard sweeps provide a means to reconstruct early AMH population dispersals OoA. We identify a previously unsuspected extended period of genetic adaptation lasting ~30,000 y, potentially in the Arabian Peninsula area, prior to a major Neandertal genetic introgression and subsequent rapid dispersal across Eurasia as far as Australia. Consistent functional targets of selection initiated during this period, which we term the Arabian Standstill, include loci involved in the regulation of fat storage, neural development, skin physiology, and cilia function. Similar adaptive signatures are also evident in introgressed archaic hominin loci and modern Arctic human groups, and we suggest that this signal represents selection for cold adaptation. Surprisingly, many of the candidate selected loci across these groups appear to directly interact and coordinately regulate biological processes, with a number associated with major modern diseases including the ciliopathies, metabolic syndrome, and neurodegenerative disorders. This expands the potential for ancestral human adaptation to directly impact modern diseases, providing a platform for evolutionary medicine.

The impact of the cytoplasmic ubiquitin ligase TNFAIP3 gene variation on transcription factor NF-κB activation in acute kidney injury.

Nuclear factor κB (NF-κB) activation is a deleterious molecular mechanism that drives acute kidney injury (AKI) and manifests in transplanted kidneys as delayed graft function. The TNFAIP3 gene encodes A20, a cytoplasmic ubiquitin ligase and a master negative regulator of the NF- κB signaling pathway. Common population-specific TNFAIP3 coding variants that reduce A20's enzyme function and increase NF- κB activation have been linked to heightened protective immunity and autoimmune disease, but have not been investigated in AKI. Here, we functionally identified a series of unique human TNFAIP3 coding variants linked to the autoimmune genome-wide association studies single nucleotide polymorphisms of F127C; namely F127C;R22Q, F127C;G281E, F127C;W448C and F127C;N449K that reduce A20's anti-inflammatory function in an NF- κB reporter assay. To investigate the impact of TNFAIP3 hypomorphic coding variants in AKI we tested a mouse Tnfaip3 hypomorph in a model of ischemia reperfusion injury (IRI). The mouse Tnfaip3 coding variant I325N increases NF- κB activation without overt inflammatory disease, providing an immune boost as I325N mice exhibit enhanced innate immunity to a bacterial challenge. Surprisingly, despite exhibiting increased intra-kidney NF- κB activation with inflammation in IRI, the kidney of I325N mice was protected. The I325N variant influenced the outcome of IRI by changing the dynamic expression of multiple cytoprotective mechanisms, particularly by increasing NF- κB-dependent anti-apoptotic factors BCL-2, BCL-XL, c-FLIP and A20, altering the active redox state of the kidney with a reduction of superoxide levels and the enzyme super oxide dismutase-1, and enhancing cellular protective mechanisms including increased Foxp3+ T cells. Thus, TNFAIP3 gene variants represent a kidney and population-specific molecular factor that can dictate the course of IRI.

Admixture has obscured signals of historical hard sweeps in humans.

The role of natural selection in shaping biological diversity is an area of intense interest in modern biology. To date, studies of positive selection have primarily relied on genomic datasets from contemporary populations, which are susceptible to confounding factors associated with complex and often unknown aspects of population history. In particular, admixture between diverged populations can distort or hide prior selection events in modern genomes, though this process is not explicitly accounted for in most selection studies despite its apparent ubiquity in humans and other species. Through analyses of ancient and modern human genomes, we show that previously reported Holocene-era admixture has masked more than 50 historic hard sweeps in modern European genomes. Our results imply that this canonical mode of selection has probably been underappreciated in the evolutionary history of humans and suggest that our current understanding of the tempo and mode of selection in natural populations may be inaccurate.

Response to Comment on "A global environmental crisis 42,000 years ago".

Our study on the exact timing and the potential climatic, environmental, and evolutionary consequences of the Laschamps Geomagnetic Excursion has generated the hypothesis that geomagnetism represents an unrecognized driver in environmental and evolutionary change. It is important for this hypothesis to be tested with new data, and encouragingly, none of the studies presented by Picin et al. undermine our model.

Response to Comment on "A global environmental crisis 42,000 years ago".

Our paper about the impacts of the Laschamps Geomagnetic Excursion 42,000 years ago has provoked considerable scientific and public interest, particularly in the so-called Adams Event associated with the initial transition of the magnetic poles. Although we welcome the opportunity to discuss our new ideas, Hawks' assertions of misrepresentation are especially disappointing given his limited examination of the material.

Ancient mitochondrial genomes from the Argentinian Pampas inform the early peopling of the Southern Cone of South America.

The Southern Cone of South America (SCSA) is a key region for investigations about the peopling of the Americas. However, little is known about the eastern sector, the Argentinian Pampas. We analyzed 18 mitochondrial genomes-7 of which are novel-from human skeletal remains from 3 Early to Late Holocene archaeological sites. The Pampas present a distinctive genetic makeup compared to other Middle to Late Holocene pre-Columbian SCSA populations. We also report the earliest individuals carrying SCSA-specific mitochondrial haplogroups D1j and D1g from Early and Middle Holocene, respectively. Using these deep calibration time points in Bayesian phylogenetic reconstructions, we suggest that the first settlers of the Pampas were part of a single and rapid dispersal ∼15,600 years ago. Finally, we propose that present-day genetic differences between the Pampas and the rest of the SCSA are due to founder effects, genetic drift, and a partial population replacement ∼9,000 years ago.

An ancient viral epidemic involving host coronavirus interacting genes more than 20,000 years ago in East Asia.

The current severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has emphasized the vulnerability of human populations to novel viral pressures, despite the vast array of epidemiological and biomedical tools now available. Notably, modern human genomes contain evolutionary information tracing back tens of thousands of years, which may help identify the viruses that have impacted our ancestors-pointing to which viruses have future pandemic potential. Here, we apply evolutionary analyses to human genomic datasets to recover selection events involving tens of human genes that interact with coronaviruses, including SARS-CoV-2, that likely started more than 20,000 years ago. These adaptive events were limited to the population ancestral to East Asian populations. Multiple lines of functional evidence support an ancient viral selective pressure, and East Asia is the geographical origin of several modern coronavirus epidemics. An arms race with an ancient coronavirus, or with a different virus that happened to use similar interactions as coronaviruses with human hosts, may thus have taken place in ancestral East Asian populations. By learning more about our ancient viral foes, our study highlights the promise of evolutionary information to better predict the pandemics of the future. Importantly, adaptation to ancient viral epidemics in specific human populations does not necessarily imply any difference in genetic susceptibility between different human populations, and the current evidence points toward an overwhelming impact of socioeconomic factors in the case of coronavirus disease 2019 (COVID-19).

Systematic benchmark of ancient DNA read mapping.

The current standard practice for assembling individual genomes involves mapping millions of short DNA sequences (also known as DNA 'reads') against a pre-constructed reference genome. Mapping vast amounts of short reads in a timely manner is a computationally challenging task that inevitably produces artefacts, including biases against alleles not found in the reference genome. This reference bias and other mapping artefacts are expected to be exacerbated in ancient DNA (aDNA) studies, which rely on the analysis of low quantities of damaged and very short DNA fragments (~30-80 bp). Nevertheless, the current gold-standard mapping strategies for aDNA studies have effectively remained unchanged for nearly a decade, during which time new software has emerged. In this study, we used simulated aDNA reads from three different human populations to benchmark the performance of 30 distinct mapping strategies implemented across four different read mapping software-BWA-aln, BWA-mem, NovoAlign and Bowtie2-and quantified the impact of reference bias in downstream population genetic analyses. We show that specific NovoAlign, BWA-aln and BWA-mem parameterizations achieve high mapping precision with low levels of reference bias, particularly after filtering out reads with low mapping qualities. However, unbiased NovoAlign results required the use of an IUPAC reference genome. While relevant only to aDNA projects where reference population data are available, the benefit of using an IUPAC reference demonstrates the value of incorporating population genetic information into the aDNA mapping process, echoing recent results based on graph genome representations.

Genozip: a universal extensible genomic data compressor.

We present Genozip, a universal and fully featured compression software for genomic data. Genozip is designed to be a general-purpose software and a development framework for genomic compression by providing five core capabilities-universality (support for all common genomic file formats), high compression ratios, speed, feature-richness and extensibility. Genozip delivers high-performance compression for widelyused genomic data formats in genomics research, namely FASTQ, SAM/BAM/CRAM, VCF, GVF, FASTA, PHYLIP and 23andMe formats. Our test results show that Genozip is fast and achieves greatly improved compression ratios, even when the files are already compressed. Further, Genozip is architected with a separation of the Genozip Framework from file-format-specific Segmenters and data-type-specific Codecs. With this, we intend for Genozip to be a general-purpose compression platform where researchers can implement compression for additional file formats, as well as new codecs for data types or fields within files, in the future. We anticipate that this will ultimately increase the visibility and adoption of these algorithms by the user community, thereby accelerating further innovation in this space. AVAILABILITY AND IMPLEMENTATION: Genozip is written in C. The code is open-source and available on http://www.genozip.com. The package is free for non-commercial use. It is distributed through the Conda package manager, github, and as a Docker container on DockerHub. Genozip is tested on Linux, Mac and Windows. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A global environmental crisis 42,000 years ago.

Geological archives record multiple reversals of Earth's magnetic poles, but the global impacts of these events, if any, remain unclear. Uncertain radiocarbon calibration has limited investigation of the potential effects of the last major magnetic inversion, known as the Laschamps Excursion [41 to 42 thousand years ago (ka)]. We use ancient New Zealand kauri trees (Agathis australis) to develop a detailed record of atmospheric radiocarbon levels across the Laschamps Excursion. We precisely characterize the geomagnetic reversal and perform global chemistry-climate modeling and detailed radiocarbon dating of paleoenvironmental records to investigate impacts. We find that geomagnetic field minima ~42 ka, in combination with Grand Solar Minima, caused substantial changes in atmospheric ozone concentration and circulation, driving synchronous global climate shifts that caused major environmental changes, extinction events, and transformations in the archaeological record.

Additional evaluations show that specific BWA-aln settings still outperform BWA-mem for ancient DNA data alignment.

Xu et al. (2021) recently recommended a new parameterization of BWA-mem as a superior alternative to the widely-used BWA-aln algorithm to map ancient DNA sequencing data. Here, we compare the BWA-mem parameterization recommended by Xu et al. with the best-performing alignment methods determined in the recent benchmarks of Oliva and colleagues (2021), demonstrating that BWA-aln is still the gold-standard for ancient DNA read alignment .

Ancient DNA Studies in Pre-Columbian Mesoamerica.

Mesoamerica is a historically and culturally defined geographic area comprising current central and south Mexico, Belize, Guatemala, El Salvador, and border regions of Honduras, western Nicaragua, and northwestern Costa Rica. The permanent settling of Mesoamerica was accompanied by the development of agriculture and pottery manufacturing (2500 BCE-150 CE), which led to the rise of several cultures connected by commerce and farming. Hence, Mesoamericans probably carried an invaluable genetic diversity partly lost during the Spanish conquest and the subsequent colonial period. Mesoamerican ancient DNA (aDNA) research has mainly focused on the study of mitochondrial DNA in the Basin of Mexico and the Yucatán Peninsula and its nearby territories, particularly during the Postclassic period (900-1519 CE). Despite limitations associated with the poor preservation of samples in tropical areas, recent methodological improvements pave the way for a deeper analysis of Mesoamerica. Here, we review how aDNA research has helped discern population dynamics patterns in the pre-Columbian Mesoamerican context, how it supports archaeological, linguistic, and anthropological conclusions, and finally, how it offers new working hypotheses.

Equitable Expanded Carrier Screening Needs Indigenous Clinical and Population Genomic Data.

Expanded carrier screening (ECS) for recessive monogenic diseases requires prior knowledge of genomic variation, including DNA variants that cause disease. The composition of pathogenic variants differs greatly among human populations, but historically, research about monogenic diseases has focused mainly on people with European ancestry. By comparison, less is known about pathogenic DNA variants in people from other parts of the world. Consequently, inclusion of currently underrepresented Indigenous and other minority population groups in genomic research is essential to enable equitable outcomes in ECS and other areas of genomic medicine. Here, we discuss this issue in relation to the implementation of ECS in Australia, which is currently being evaluated as part of the national Government's Genomics Health Futures Mission. We argue that significant effort is required to build an evidence base and genomic reference data so that ECS can bring significant clinical benefit for many Aboriginal and/or Torres Strait Islander Australians. These efforts are essential steps to achieving the Australian Government's objectives and its commitment "to leveraging the benefits of genomics in the health system for all Australians." They require culturally safe, community-led research and community involvement embedded within national health and medical genomics programs to ensure that new knowledge is integrated into medicine and health services in ways that address the specific and articulated cultural and health needs of Indigenous people. Until this occurs, people who do not have European ancestry are at risk of being, in relative terms, further disadvantaged.

genozip: a fast and efficient compression tool for VCF files.

MOTIVATION: genozip is a new lossless compression tool for Variant Call Format (VCF) files. By applying field-specific algorithms and fully utilizing the available computational hardware, genozip achieves the highest compression ratios amongst existing lossless compression tools known to the authors, at speeds comparable with the fastest multi-threaded compressors. AVAILABILITY AND IMPLEMENTATION: genozip is freely available to non-commercial users. It can be installed via conda-forge, Docker Hub, or downloaded from github.com/divonlan/genozip. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The Dogma of Dingoes-Taxonomic status of the dingo: A reply to Smith et al.

Adopting the name Canis dingo for the Dingo to explicitly denote a species-level taxon separate from other canids was suggested by Crowther et al.  (2014) as a means to eliminate taxonomic instability and contention. However, Jackson et al.  (2017), using standard taxonomic and nomenclatural approaches and principles, called instead for continued use of the nomen C. familiaris for all domestic dogs and their derivatives, including the Dingo. (This name, C. familiaris, is applied to all dogs that derive from the domesticated version of the Gray Wolf, Canis lupus, based on nomenclatural convention.) The primary reasons for this call by Jackson et al.  (2017) were: (1) a lack of evidence to show that recognizing multiple species amongst the dog, including the Dingo and New Guinea Singing Dog, was necessary taxonomically, and (2) the principle of nomenclatural priority (the name familiaris Linnaeus, 1758, antedates dingo Meyer, 1793). Overwhelming current evidence from archaeology and genomics indicates that the Dingo is of recent origin in Australia and shares immediate ancestry with other domestic dogs as evidenced by patterns of genetic and morphological variation. Accordingly, for Smith et al.  (2019) to recognise Canis dingo as a distinct species, the onus was on them to overturn current interpretations of available archaeological, genomic, and morphological datasets and instead show that Dingoes have a deeply divergent evolutionary history that distinguishes them from other named forms of Canis (including C. lupus and its domesticated version, C. familiaris). A recent paper by Koepfli et al.  (2015) demonstrates exactly how this can be done in a compelling way within the genus Canis-by demonstrating deep evolutionary divergence between taxa, on the order of hundreds of thousands of years, using data from multiple genetic systems. Smith et al.  (2019) have not done this; instead they have misrepresented the content and conclusions of Jackson et al.  (2017), and contributed extraneous arguments that are not relevant to taxonomic decisions. Here we dissect Smith et al.  (2019), identifying misrepresentations, to show that ecological, behavioural and morphological evidence is insufficient to recognise Dingoes as a separate species from other domestic dogs. We reiterate: the correct binomial name for the taxon derived from Gray Wolves (C. lupus) by passive and active domestication, including Dingoes and other domestic dogs, is Canis familiaris. We are strongly sympathetic to arguments about the historical, ecological, cultural, or other significance of the Dingo, but these are issues that will have to be considered outside of the more narrow scope of taxonomy and nomenclature.

A strategy for building and using a human reference pangenome.

In March 2019, 45 scientists and software engineers from around the world converged at the University of California, Santa Cruz for the first pangenomics codeathon. The purpose of the meeting was to propose technical specifications and standards for a usable human pangenome as well as to build relevant tools for genome graph infrastructures. During the meeting, the group held several intense and productive discussions covering a diverse set of topics, including advantages of graph genomes over a linear reference representation, design of new methods that can leverage graph-based data structures, and novel visualization and annotation approaches for pangenomes. Additionally, the participants self-organized themselves into teams that worked intensely over a three-day period to build a set of pipelines and tools for specific pangenomic applications. A summary of the questions raised and the tools developed are reported in this manuscript.

Developing reproducible bioinformatics analysis workflows for heterogeneous computing environments to support African genomics.

BACKGROUND: The Pan-African bioinformatics network, H3ABioNet, comprises 27 research institutions in 17 African countries. H3ABioNet is part of the Human Health and Heredity in Africa program (H3Africa), an African-led research consortium funded by the US National Institutes of Health and the UK Wellcome Trust, aimed at using genomics to study and improve the health of Africans. A key role of H3ABioNet is to support H3Africa projects by building bioinformatics infrastructure such as portable and reproducible bioinformatics workflows for use on heterogeneous African computing environments. Processing and analysis of genomic data is an example of a big data application requiring complex interdependent data analysis workflows. Such bioinformatics workflows take the primary and secondary input data through several computationally-intensive processing steps using different software packages, where some of the outputs form inputs for other steps. Implementing scalable, reproducible, portable and easy-to-use workflows is particularly challenging. RESULTS: H3ABioNet has built four workflows to support (1) the calling of variants from high-throughput sequencing data; (2) the analysis of microbial populations from 16S rDNA sequence data; (3) genotyping and genome-wide association studies; and (4) single nucleotide polymorphism imputation. A week-long hackathon was organized in August 2016 with participants from six African bioinformatics groups, and US and European collaborators. Two of the workflows are built using the Common Workflow Language framework (CWL) and two using Nextflow. All the workflows are containerized for improved portability and reproducibility using Docker, and are publicly available for use by members of the H3Africa consortium and the international research community. CONCLUSION: The H3ABioNet workflows have been implemented in view of offering ease of use for the end user and high levels of reproducibility and portability, all while following modern state of the art bioinformatics data processing protocols. The H3ABioNet workflows will service the H3Africa consortium projects and are currently in use. All four workflows are also publicly available for research scientists worldwide to use and adapt for their respective needs. The H3ABioNet workflows will help develop bioinformatics capacity and assist genomics research within Africa and serve to increase the scientific output of H3Africa and its Pan-African Bioinformatics Network.

Organizing and running bioinformatics hackathons within Africa: The H3ABioNet cloud computing experience.

The need for portable and reproducible genomics analysis pipelines is growing globally as well as in Africa, especially with the growth of collaborative projects like the Human Health and Heredity in Africa Consortium (H3Africa). The Pan-African H3Africa Bioinformatics Network (H3ABioNet) recognized the need for portable, reproducible pipelines adapted to heterogeneous compute environments, and for the nurturing of technical expertise in workflow languages and containerization technologies. To address this need, in 2016 H3ABioNet arranged its first Cloud Computing and Reproducible Workflows Hackathon, with the purpose of building key genomics analysis pipelines able to run on heterogeneous computing environments and meeting the needs of H3Africa research projects. This paper describes the preparations for this hackathon and reflects upon the lessons learned about its impact on building the technical and scientific expertise of African researchers. The workflows developed were made publicly available in GitHub repositories and deposited as container images on quay.io.