Phylogenetic assessment within a species complex of a subterranean rodent (Geomys bursarius) with conservation implications for isolated subspecies.

Range contraction and expansion from glaciation has led to genetic divergence that may be particularly pronounced in fossorial species with low dispersal. The plains pocket gopher (Geomys bursarius) is a fossorial species that ranges widely across North America but has a poorly understood phylogeny. We used mitogenomes (14,996 base pairs) from 56 individuals across seven subspecies, plus two outgroup species, to assess genetic divergence from minimum spanning trees, measure genetic distances, and infer phylogenetic trees using BEAST. We found G. b. wisconsinensis was monophyletic with recent divergence. Further assessment is needed for G. b. major because it was paraphyletic and exhibited inconsistent groupings with other clades. Importantly, we identified G. b. illinoensis as being genetically distinct and monophyletic likely due to a unique colonization event eastward across the Mississippi River. Because G. b. illinoensis faces continued pressures from niche reduction and habitat loss, we recommend that G. b. illinoensis be considered an evolutionary significant unit warranting conservation actions to promote connectivity and restore suitable habitat. Such conservation efforts should benefit other grassland species including those originating from clades west of the Mississippi River that may also be evolutionary significant units.

Genomic analyses correspond with deep persistence of peoples of Blackfoot Confederacy from glacial times.

Mutually beneficial partnerships between genomics researchers and North American Indigenous Nations are rare yet becoming more common. Here, we present one such partnership that provides insight into the peopling of the Americas and furnishes another line of evidence that can be used to further treaty and Indigenous rights. We show that the genomics of sampled individuals from the Blackfoot Confederacy belong to a previously undescribed ancient lineage that diverged from other genomic lineages in the Americas in Late Pleistocene times. Using multiple complementary forms of knowledge, we provide a scenario for Blackfoot population history that fits with oral tradition and provides a plausible model for the evolutionary process of the peopling of the Americas.

Genomic analysis supports Cape Lion population connectivity prior to colonial eradication and extinction.

Cape lions (Panthera leo melanochaitus) formerly ranged throughout the grassland plains of the "Cape Flats" in what is today known as the Western Cape Province, South Africa. Cape lions were likely eradicated because of overhunting and habitat loss after European colonization. European naturalists originally described Cape lions as "black-maned lions" and claimed that they were phenotypically distinct. However, other depictions and historical descriptions of lions from the Cape report mixed or light coloration and without black or extensively developed manes. These findings suggest that, rather than forming a distinct population, Cape lions may have had phenotypic and genotypic variation similar to other African lions. Here we investigate Cape lion genome characteristics, population dynamics, and genetic distinctiveness prior to their extinction. We generated genomic data from 2 historic Cape lions to compare to 118 existing high-coverage mitogenomes, and low-coverage nuclear genomes of 53 lions from 13 African countries. We show that, before their eradication, lions from the Cape Flats had diverse mitogenomes and nuclear genomes that clustered with lions from both southern and eastern Africa. Cape lions had high genome-wide heterozygosity and low inbreeding coefficients, indicating that populations in the Cape Flats went extinct so rapidly that genomic effects associated with long-term small population size and isolation were not detectable. Our findings do not support the characterization of Cape lions as phylogeographically distinct, as originally put forth by some European naturalists, and illustrates how alternative knowledge systems, for example, Indigenous perspectives, could potentially further inform interpretations of species histories.

Association between gene methylation and experiences of historical trauma in Alaska Native peoples.

BACKGROUND: Historical trauma experienced by Indigenous peoples of North America is correlated with health disparities and is hypothesized to be associated with DNA methylation. Massive group traumas such as genocide, loss of land and foodways, and forced conversion to Western lifeways may be embodied and affect individuals, families, communities, cultures, and health. This study approaches research with Alaska Native people using a community-engaged approach designed to create mutually-beneficial partnerships, including intentional relationship development, capacity building, and sample and data care. METHODS: A total of 117 Alaska Native individuals from two regions of Alaska joined the research study. Participants completed surveys on cultural identification, historical trauma (historical loss scale and historical loss associated symptoms scale), and general wellbeing. Participants provided a blood sample which was used to assess DNA methylation with the Illumina Infinium MethylationEPIC array. RESULTS: We report an association between historical loss associated symptoms and DNA methylation at five CpG sites, evidencing the embodiment of historical trauma. We further report an association between cultural identification and general wellbeing, complementing evidence from oral narratives and additional studies that multiple aspects of cultural connection may buffer the effects of and/or aid in the healing process from historical trauma. CONCLUSION: A community-engaged approach emphasizes balanced partnerships between communities and researchers. Here, this approach helps better understand embodiment of historical trauma in Alaska Native peoples. This analysis reveals links between the historical trauma response and DNA methylation. Indigenous communities have been stigmatized for public health issues instead caused by systemic inequalities, social disparities, and discrimination, and we argue that the social determinants of health model in Alaska Native peoples must include the vast impact of historical trauma and ongoing colonial violence.

"Where do I even start?" Recommendations for faculty diversifying syllabi in ecology, evolution, and the life sciences.

Diversifying curricula is of increasing interest in higher education, including in ecology and evolution and allied fields. Yet, many educators may not know where to start. Here we provide a framework for meeting standard curriculum goals while enacting anti-racist and anti-colonial syllabi that is grounded in the development of a sustainable network of educators. In addition to highlighting this professional learning process and sharing the list of resources our group has developed, we provide suggestions to help educators highlight contributions of minoritized groups, explore multiple ways of knowing, and perform critical assessments of foundational views of life and environmental science fields. We further discuss the key classroom dynamics that affect the success of such anti-racist and anti-colonial initiatives. The retention and success of minoritized students in ecology and evolution depends on whether we address injustices in our fields. Our hope is that our fellow educators will use this paper to catalyze their own efforts to diversify their courses.

Numt Parser: Automated identification and removal of nuclear mitochondrial pseudogenes (numts) for accurate mitochondrial genome reconstruction in Panthera.

Nuclear mitochondrial pseudogenes (numts) may hinder the reconstruction of mtDNA genomes and affect the reliability of mtDNA datasets for phylogenetic and population genetic comparisons. Here, we present the program Numt Parser, which allows for the identification of DNA sequences that likely originate from numt pseudogene DNA. Sequencing reads are classified as originating from either numt or true cytoplasmic mitochondrial (cymt) DNA by direct comparison against cymt and numt reference sequences. Classified reads can then be parsed into cymt or numt datasets. We tested this program using whole genome shotgun-sequenced data from 2 ancient Cape lions (Panthera leo), because mtDNA is often the marker of choice for ancient DNA studies and the genus Panthera is known to have numt pseudogenes. Numt Parser decreased sequence disagreements that were likely due to numt pseudogene contamination and equalized read coverage across the mitogenome by removing reads that likely originated from numts. We compared the efficacy of Numt Parser to 2 other bioinformatic approaches that can be used to account for numt contamination. We found that Numt Parser outperformed approaches that rely only on read alignment or Basic Local Alignment Search Tool (BLAST) properties, and was effective at identifying sequences that likely originated from numts while having minimal impacts on the recovery of cymt reads. Numt Parser therefore improves the reconstruction of true mitogenomes, allowing for more accurate and robust biological inferences.

Ten simple rules for improving communication among scientists.

Communication is a fundamental part of scientific development and methodology. With the advancement of the internet and social networks, communication has become rapid and sometimes overwhelming, especially in science. It is important to provide scientists with useful, effective, and dynamic tools to establish and build a fluid communication framework that allows for scientific advancement. Therefore, in this article, we present advice and recommendations that can help promote and improve science communication while respecting an adequate balance in the degree of commitment toward collaborative work. We have developed 10 rules shown in increasing order of commitment that are grouped into 3 key categories: (1) speak (based on active participation); (2) join (based on joining scientific groups); and (3) assess (based on the analysis and retrospective consideration of the weaknesses and strengths). We include examples and resources that provide actionable strategies for involvement and engagement with science communication, from basic steps to more advanced, introspective, and long-term commitments. Overall, we aim to help spread science from within and encourage and engage scientists to become involved in science communication effectively and dynamically.

Ancient and modern genomics of the Ohlone Indigenous population of California.

SignificanceCalifornia supports a high cultural and linguistic diversity of Indigenous peoples. In a partnership of researchers with the Muwekma Ohlone tribe, we studied genomes of eight present-day tribal members and 12 ancient individuals from two archaeological sites in the San Francisco Bay Area, spanning ∼2,000 y. We find that compared to genomes of Indigenous individuals from throughout the Americas, the 12 ancient individuals are most genetically similar to ancient individuals from Southern California, and that despite spanning a large time period, they share distinctive ancestry. This ancestry is also shared with present-day tribal members, providing evidence of genetic continuity between past and present Indigenous individuals in the region, in contrast to some popular reconstructions based on archaeological and linguistic information.

Combining methods for non-invasive fecal DNA enables whole genome and metagenomic analyses in wildlife biology.

Non-invasive biological samples benefit studies that investigate rare, elusive, endangered, or dangerous species. Integrating genomic techniques that use non-invasive biological sampling with advances in computational approaches can benefit and inform wildlife conservation and management. Here, we used non-invasive fecal DNA samples to generate low- to medium-coverage genomes (e.g., >90% of the complete nuclear genome at six X-fold coverage) and metagenomic sequences, combining widely available and accessible DNA collection cards with commonly used DNA extraction and library building approaches. DNA preservation cards are easy to transport and can be stored non-refrigerated, avoiding cumbersome or costly sample methods. The genomic library construction and shotgun sequencing approach did not require enrichment or targeted DNA amplification. The utility and potential of the data generated was demonstrated through genome scale and metagenomic analyses of zoo and free-ranging African savanna elephants (Loxodonta africana). Fecal samples collected from free-ranging individuals contained an average of 12.41% (5.54-21.65%) endogenous elephant DNA. Clustering of these elephants with others from the same geographic region was demonstrated by a principal component analysis of genetic variation using nuclear genome-wide SNPs. Metagenomic analyses identified taxa that included Loxodonta, green plants, fungi, arthropods, bacteria, viruses and archaea, showcasing the utility of this approach for addressing complementary questions based on host-associated DNA, e.g., pathogen and parasite identification. The molecular and bioinformatic analyses presented here contributes towards the expansion and application of genomic techniques to conservation science and practice.

Sourcing Elephant Ivory from a Sixteenth-Century Portuguese Shipwreck.

The oldest known shipwreck in southern Africa was found in Namibia in 2008.1-4 Forty tons of cargo, including gold and silver coins, helped identify the ship as the Bom Jesus, a Portuguese nau (trading vessel) lost in 1533 while headed to India.4-6 The cargo included >100 elephant tusks,7 which we examined using paleogenomic and stable isotope analyses. Nuclear DNA identified the ivory source as African forest (Loxodonta cyclotis) rather than savanna (Loxodonta africana) elephants. Mitochondrial sequences traced them to West and not Central Africa and from ≥17 herds with distinct haplotypes. Four of the haplotypes are known from modern populations; others were potentially lost to subsequent hunting of elephants for ivory. Stable isotope analyses (δ13C and δ15N) indicated that the elephants were not from deep rainforests but from savanna and mixed habitats. Such habitats surround the Guinean forest block of West Africa8 and accord with the locations of major historic Portuguese trading ports.9,10 West African forest elephants currently range into savanna habitats;11-13 our findings suggest that this was not consequent to regional decimation of savanna elephants for their ivory in the 19th and 20th centuries. During the time of the Bom Jesus, ivory was a central driver in the formation of maritime trading systems connecting Europe, Africa, and Asia. Our integration of paleogenomic, archeological, and historical methods to analyze the Bom Jesus ivory provides a framework for examining vast collections of archaeological ivories around the world, in shipwrecks and other contexts.

A comparison of proteomic, genomic, and osteological methods of archaeological sex estimation.

Sex estimation of skeletons is fundamental to many archaeological studies. Currently, three approaches are available to estimate sex-osteology, genomics, or proteomics, but little is known about the relative reliability of these methods in applied settings. We present matching osteological, shotgun-genomic, and proteomic data to estimate the sex of 55 individuals, each with an independent radiocarbon date between 2,440 and 100 cal BP, from two ancestral Ohlone sites in Central California. Sex estimation was possible in 100% of this burial sample using proteomics, in 91% using genomics, and in 51% using osteology. Agreement between the methods was high, however conflicts did occur. Genomic sex estimates were 100% consistent with proteomic and osteological estimates when DNA reads were above 100,000 total sequences. However, more than half the samples had DNA read numbers below this threshold, producing high rates of conflict with osteological and proteomic data where nine out of twenty conditional DNA sex estimates conflicted with proteomics. While the DNA signal decreased by an order of magnitude in the older burial samples, there was no decrease in proteomic signal. We conclude that proteomics provides an important complement to osteological and shotgun-genomic sex estimation.

Alternative Quantifications of Landscape Complementation to Model Gene Flow in Banded Longhorn Beetles [Typocerus v. velutinus (Olivier)].

Rapid progression of human socio-economic activities has altered the structure and function of natural landscapes. Species that rely on multiple, complementary habitat types (i.e., landscape complementation) to complete their life cycle may be especially at risk. However, such landscape complementation has received little attention in the context of landscape connectivity modeling. A previous study on flower longhorn beetles (Cerambycidae: Lepturinae) integrated landscape complementation into a continuous habitat suitability 'surface', which was then used to quantify landscape connectivity between pairs of sampling sites using gradient-surface metrics. This connectivity model was validated with molecular genetic data collected for the banded longhorn beetle (Typocerus v. velutinus) in Indiana, United States. However, this approach has not been compared to alternative models in a landscape genetics context. Here, we used a discrete land use/land cover map to calculate landscape metrics related to landscape complementation based on a patch mosaic model (PMM) as an alternative to the previously published, continuous habitat suitability model (HSM). We evaluated the HSM surface with gradient surface metrics (GSM) and with two resistance-based models (RBM) based on least cost path (LCP) and commute distance (CD), in addition to an isolation-by-distance (IBD) model based on Euclidean distance. We compared the ability of these competing models of connectivity to explain pairwise genetic distances (R ST) previously calculated from ten microsatellite genotypes of 454 beetles collected from 17 sites across Indiana, United States. Model selection with maximum likelihood population effects (MLPE) models found that GSM were most effective at explaining pairwise genetic distances as a proxy for gene flow across the landscape, followed by the landscape metrics calculated from the PMM, whereas the LCP model performed worse than both the CD and the isolation by distance model. We argue that the analysis of a continuous HSM with GSM might perform better because of their combined ability to effectively represent and quantify the continuous degree of landscape complementation (i.e., availability of complementary habitats in vicinity) found at and in-between sites, on which these beetles depend. Our findings may inform future studies that seek to model habitat connectivity in complex heterogeneous landscapes as natural habitats continue to become more fragmented in the Anthropocene.

Accurate Sex Identification of Ancient Elephant and Other Animal Remains Using Low-Coverage DNA Shotgun Sequencing Data.

Sex identification of ancient animal biological remains can benefit our understanding of historical population structure, demography and social behavior. Traditional methods for sex identification (e.g., osteological and morphometric comparisons) may be ineffective when animal remains are not well preserved, when sex distinguishing characteristics have not yet developed, or where organisms do not exhibit sex-associated phenotypic dimorphisms. Here we adapt a method developed for human sex determination so that it can be used to identify the sex of ancient and modern animal taxa. The method identifies sex by calculating the ratio of DNA reads aligning to the X chromosome to DNA reads aligning to autosomes (termed the Rx ratio). We tested the accuracy of this method using low coverage genomes from 15 modern elephants (Loxodonta africana) for which sex was known. We then applied this method to ancient elephant ivory samples for which sex was unknown, and describe how this method can be further adapted to the genomes of other taxa. This method may be especially useful when only low-coverage genomic data are obtainable. Furthermore, because this method relies on only the X and not the Y chromosome, it can be used to determine the sex of organisms for which a reference genome was obtained from a female or for which only the X chromosome is reported. Such taxa include the domestic cat, sheep, goat, and horse; and non-domesticated animals such as the Sumatran orangutan, western lowland gorilla and meerkat.

Loxodonta Localizer: A Software Tool for Inferring the Provenance of African Elephants and Their Ivory Using Mitochondrial DNA.

Illegal hunting is a major threat to the elephants of Africa, with more elephants killed by poachers than die from natural causes. DNA from tusks has been used to infer the source populations for confiscated ivory, relying on nuclear genetic markers. However, mitochondrial DNA (mtDNA) sequences can also provide information on the geographic origins of elephants due to female elephant philopatry. Here, we introduce the Loxodonta Localizer (LL; www.loxodontalocalizer.org), an interactive software tool that uses a database of mtDNA sequences compiled from previously published studies to provide information on the potential provenance of confiscated ivory. A 316 bp control region sequence, which can be readily generated from DNA extracted from ivory, is used as a query. The software generates a listing of haplotypes reported among 1917 African elephants in 24 range countries, sorted in order of similarity to the query sequence. The African locations from which haplotype sequences have been previously reported are shown on a map. We demonstrate examples of haplotypes reported from only a single locality or country, examine the utility of the program in identifying elephants from countries with varying degrees of sampling, and analyze batches of confiscated ivory. The LL allows for the source of confiscated ivory to be assessed within days, using widely available molecular methods that do not depend on a particular platform or laboratory. The program enables identification of potential regions or localities from which elephants are being poached, with capacity for rapid identification of populations newly or consistently targeted by poachers.

Species identification and mitochondrial genomes of ancient fish bones from the Riverine Kachemak tradition of the Kenai Peninsula, Alaska.

Seven fish vertebrae were chosen for analysis from the 49-KEN-147 archaeological site in the Kenai Peninsula, Alaska. Mitochondrial DNA analysis of the ancient fish bones revealed that they were from sockeye and coho salmon. Here, we report the ancient mitochondrial genomes for three sockeye salmon and one coho salmon fish bone.