Search

Dingoes are culturally and ecologically important free-living canids whose ancestors arrived in Australia over 3,000 B.P., likely transported by seafaring people. However, the early history of dingoes in Australia-including the number of founding populations and their routes of introduction-remains uncertain. This uncertainty arises partly from the complex and poorly understood relationship between modern dingoes and New Guinea singing dogs, and suspicions that post-Colonial hybridization has introduced recent domestic dog ancestry into the genomes of many wild dingo populations. In this study, we analyzed genome-wide data from nine ancient dingo specimens ranging in age from 400 to 2,746 y old, predating the introduction of domestic dogs to Australia by European colonists. We uncovered evidence that the continent-wide population structure observed in modern dingo populations had already emerged several thousand years ago. We also detected excess allele sharing between New Guinea singing dogs and ancient dingoes from coastal New South Wales (NSW) compared to ancient dingoes from southern Australia, irrespective of any post-Colonial hybrid ancestry in the genomes of modern individuals. Our results are consistent with several demographic scenarios, including a scenario where the ancestry of dingoes from the east coast of Australia results from at least two waves of migration from source populations with varying affinities to New Guinea singing dogs. We also contribute to the growing body of evidence that modern dingoes derive little genomic ancestry from post-Colonial hybridization with other domestic dog lineages, instead descending primarily from ancient canids introduced to Sahul thousands of years ago., Competing Interests: Competing interests statement:The authors declare no competing interest.

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., (Copyright © 2023 International Society of Nephrology. All rights reserved.)

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 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., (© 2022. The Author(s).)

Competing Interests: Declaration of interests We declare no competing interests.

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 ., Competing Interests: The authors declare no conflict of interest., (© 2021 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.)

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.

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., (© The Author(s) 2021. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

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., (© The Author(s) 2021. Published by Oxford University Press.)

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., Competing Interests: The authors declare no competing interests., (© 2021 The Author(s).)

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., (Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

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.

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., (Copyright © 2020 American Society of Human Genetics. All rights reserved.)

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 Author(s) 2020. Published by Oxford University Press.)

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., Competing Interests: No competing interests were disclosed., (Copyright: © 2019 Llamas B et al.)

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 computing environments, and for the nurturing of technical expertise in workflow languages and containerization technologies. Building on the network's Standard Operating Procedures (SOPs) for common genomic analyses, H3ABioNet arranged its first Cloud Computing and Reproducible Workflows Hackathon in 2016, with the purpose of translating those SOPs into 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., Competing Interests: Competing interests: MRC in his role as CWL Community Engineer, has had his salary supported in the past by grants from Seven Bridges Genomics, Inc to its employers. The other authors declare that they have no competing interests., (Copyright: © 2019 Ahmed AE 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.

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.

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., Competing Interests: Competing interests: MRC in his role as CWL Community Engineer, has had his salary supported in the past by grants from Seven Bridges Genomics, Inc to its employers. The other authors declare that they have no competing interests., (Copyright: © 2018 Ahmed AE et al.)

Background: Next Generation Genome sequencing techniques became affordable for massive sequencing efforts devoted to clinical characterization of human diseases. However, the cost of providing cloud-based data analysis of the mounting datasets remains a concerning bottleneck for providing cost-effective clinical services. To address this computational problem, it is important to optimize the variant analysis workflow and the used analysis tools to reduce the overall computational processing time, and concomitantly reduce the processing cost. Furthermore, it is important to capitalize on the use of the recent development in the cloud computing market, which have witnessed more providers competing in terms of products and prices., Results: In this paper, we present a new package called MC-GenomeKey (Multi-Cloud GenomeKey) that efficiently executes the variant analysis workflow for detecting and annotating mutations using cloud resources from different commercial cloud providers. Our package supports Amazon, Google, and Azure clouds, as well as, any other cloud platform based on OpenStack. Our package allows different scenarios of execution with different levels of sophistication, up to the one where a workflow can be executed using a cluster whose nodes come from different clouds. MC-GenomeKey also supports scenarios to exploit the spot instance model of Amazon in combination with the use of other cloud platforms to provide significant cost reduction. To the best of our knowledge, this is the first solution that optimizes the execution of the workflow using computational resources from different cloud providers., Conclusions: MC-GenomeKey provides an efficient multicloud based solution to detect and annotate mutations. The package can run in different commercial cloud platforms, which enables the user to seize the best offers. The package also provides a reliable means to make use of the low-cost spot instance model of Amazon, as it provides an efficient solution to the sudden termination of spot machines as a result of a sudden price increase. The package has a web-interface and it is available for free for academic use.

Background: While next-generation sequencing (NGS) costs have plummeted in recent years, cost and complexity of computation remain substantial barriers to the use of NGS in routine clinical care. The clinical potential of NGS will not be realized until robust and routine whole genome sequencing data can be accurately rendered to medically actionable reports within a time window of hours and at scales of economy in the 10's of dollars., Results: We take a step towards addressing this challenge, by using COSMOS, a cloud-enabled workflow management system, to develop GenomeKey, an NGS whole genome analysis workflow. COSMOS implements complex workflows making optimal use of high-performance compute clusters. Here we show that the Amazon Web Service (AWS) implementation of GenomeKey via COSMOS provides a fast, scalable, and cost-effective analysis of both public benchmarking and large-scale heterogeneous clinical NGS datasets., Conclusions: Our systematic benchmarking reveals important new insights and considerations to produce clinical turn-around of whole genome analysis optimization and workflow management including strategic batching of individual genomes and efficient cluster resource configuration.

This is a summary of the activities and scientific content of the first International Society for Computational Biology Student Council symposium in Africa. This meeting organized by the students for the students took place 8th of March 2015 in Dar Es Salaam, Tanzania.

Summary: Efficient workflows to shepherd clinically generated genomic data through the multiple stages of a next-generation sequencing pipeline are of critical importance in translational biomedical science. Here we present COSMOS, a Python library for workflow management that allows formal description of pipelines and partitioning of jobs. In addition, it includes a user interface for tracking the progress of jobs, abstraction of the queuing system and fine-grained control over the workflow. Workflows can be created on traditional computing clusters as well as cloud-based services., Availability and Implementation: Source code is available for academic non-commercial research purposes. Links to code and documentation are provided at http://lpm.hms.harvard.edu and http://wall-lab.stanford.edu., Contact: dpwall@stanford.edu or peter_tonellato@hms.harvard.edu., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author 2014. Published by Oxford University Press.)