Systematic benchmark of ancient DNA read mapping

Identification of the 18 World War II executed citizens of Adele, Rethymnon, Crete using an ancient DNA approach and low coverage genomes

In the Battle of Crete during the World War II occupation of Greece, the German forces faced substantial civilian resistance. To retribute the numerous German losses, a series of mass executions took place in numerous places in Crete; a common practice reported from Greece and elsewhere. In Adele, a village in the regional unit of Rethymnon, 18 male civilians were executed and buried in a burial pit at the Sarakina site. In this study, the first one conducted for a conflict that occurred in Greece, we identified for humanitarian purposes the 18 skulls of the Sarakina victims, following a request from the local community of Adele. The molecular identification of historical human remains via ancient DNA approaches and low coverage whole genome sequencing has only recently been introduced. Here, we performed genome skimming on the living relatives of the victims, as well as high throughput historical DNA analysis on the skulls to infer the kinship degrees among the victims via genetic relatedness analyses. We also conducted targeted anthropological analysis to successfully complete the identification of all Sarakina victims. We demonstrate that our methodological approach constitutes a potentially highly informative forensic tool to identify war victims. It can hence be applied to analogous studies on degraded DNA, thus, paving the path for systematic war victim identification in Greece and beyond.

Disregarding multimappers leads to biases in the functional assessment of NGS data

BACKGROUND

Standard ChIP-seq and RNA-seq processing pipelines typically disregard sequencing reads whose origin is ambiguous ("multimappers"). This usual practice has potentially important consequences for the functional interpretation of the data: genomic elements belonging to clusters composed of highly similar members are left unexplored.

RESULTS

In particular, disregarding multimappers leads to the underrepresentation in epigenetic studies of recently active transposable elements, such as AluYa5, L1HS and SVAs. Furthermore, this common strategy also has implications for transcriptomic analysis: members of repetitive gene families, such the ones including major histocompatibility complex (MHC) class I and II genes, are under-quantified.

CONCLUSION

Revealing inherent biases that permeate routine tasks such as functional enrichment analysis, our results underscore the urgency of broadly adopting multimapper-aware bioinformatic pipelines -currently restricted to specific contexts or communities- to ensure the reliability of genomic and transcriptomic studies.

Colonial-driven extinction of the blue antelope despite genomic adaptation to low population size

Low genomic diversity is generally indicative of small population size and is considered detrimental by decreasing long-term adaptability.1,2,3,4,5,6 Moreover, small population size may promote gene flow with congeners and outbreeding depression.7,8,9,10,11,12,13 Here, we examine the connection between habitat availability, effective population size (Ne), and extinction by generating a 40× nuclear genome from the extinct blue antelope (Hippotragus leucophaeus). Historically endemic to the relatively small Cape Floristic Region in southernmost Africa,14,15 populations were thought to have expanded and contracted across glacial-interglacial cycles, tracking suitable habitat.16,17,18 However, we found long-term low Ne, unaffected by glacial cycles, suggesting persistence with low genomic diversity for many millennia prior to extinction in ∼AD 1800. A lack of inbreeding, alongside high levels of genetic purging, suggests adaptation to this long-term low Ne and that human impacts during the colonial era (e.g., hunting and landscape transformation), rather than longer-term ecological processes, were central to its extinction. Phylogenomic analyses uncovered gene flow between roan (H. equinus) and blue antelope, as well as between roan and sable antelope (H. niger), approximately at the time of divergence of blue and sable antelope (∼1.9 Ma). Finally, we identified the LYST and ASIP genes as candidates for the eponymous bluish pelt color of the blue antelope. Our results revise numerous aspects of our understanding of the interplay between genomic diversity and evolutionary history and provide the resources for uncovering the genetic basis of this extinct species' unique traits.

Shotgun metagenomics reveals the gut microbial diversity and functions in Vespa mandarinia (Hymenoptera: Vespidae) at multiple life stages

Wasps play important roles as predators and pollinators in the ecosystem. The Jingpo minority residing in Yunnan Province, China, has a traditional practice of using wine infused with mature wasps as a customary remedy for managing rheumatoid arthritis. The larva of the wasp is also a tasteful folk dish that has created a tremendous market. There is a paucity of survival knowledge, which has greatly restricted their potential applications in food and healthcare. Recent research has highlighted the importance of gut microbiota in insect growth. Nevertheless, there is still a lack of understanding regarding the composition, changes, and functions of the gut microbiota in Vespa mandarinia during development. In this research, the gut microbiota were investigated across three growth stages of Vespa mandarinia using a metagenomic technology. The result revealed that there are significant variations in the proportion of main gut microbes during the metamorphosis of Vespa mandarinia. Tenericutes were found to dominate during the larval stage, while Proteobacteria emerged as the dominant group post-pupation. Through a comprehensive analysis of the gut microbiota metagenome, this study revealed functional differences in the wasp gut microbiota at various growth stages. During the larval stage, the gut microbiota plays a central role in promoting metabolism. Following pupation, the gut microbiota exhibited diversified functions, likely due to the complex environments and diverse food sources encountered after metamorphosis. These functions included amino acid metabolism, compound degradation, and defense mechanisms. This research provides an extensive dataset on the gut microbiota during the metamorphosis of Vespa mandarinia, contributing to a deeper understanding of the influence of gut microbiota on wasp growth. Furthermore, this study uncovers a unique microbial treasure within insect guts, which is important for advancing the application of wasps in the fields of food and medicine.

Benchmarking a targeted 16S ribosomal RNA gene enrichment approach to reconstruct ancient microbial communities

The taxonomic characterization of ancient microbiomes is a key step in the rapidly growing field of paleomicrobiology. While PCR amplification of the 16S ribosomal RNA (rRNA) gene is a widely used technique in modern microbiota studies, this method has systematic biases when applied to ancient microbial DNA. Shotgun metagenomic sequencing has proven to be the most effective method in reconstructing taxonomic profiles of ancient dental calculus samples. Nevertheless, shotgun sequencing approaches come with inherent limitations that could be addressed through hybridization enrichment capture. When employed together, shotgun sequencing and hybridization capture have the potential to enhance the characterization of ancient microbial communities. Here, we develop, test, and apply a hybridization enrichment capture technique to selectively target 16S rRNA gene fragments from the libraries of ancient dental calculus samples generated with shotgun techniques. We simulated data sets generated from hybridization enrichment capture, indicating that taxonomic identification of fragmented and damaged 16S rRNA gene sequences was feasible. Applying this enrichment approach to 15 previously published ancient calculus samples, we observed a 334-fold increase of ancient 16S rRNA gene fragments in the enriched samples when compared to unenriched libraries. Our results suggest that 16S hybridization capture is less prone to the effects of background contamination than 16S rRNA amplification, yielding a higher percentage of on-target recovery. While our enrichment technique detected low abundant and rare taxa within a given sample, these assignments may not achieve the same level of specificity as those achieved by unenriched methods.

Allelic bias when performing in-solution enrichment of ancient human DNA

In-solution hybridisation enrichment of genetic variation is a valuable methodology in human paleogenomics. It allows enrichment of endogenous DNA by targeting genetic markers that are comparable between sequencing libraries. Many studies have used the 1240k reagent-which enriches 1,237,207 genome-wide SNPs-since 2015, though access was restricted. In 2021, Twist Biosciences and Daicel Arbor Biosciences independently released commercial kits that enabled all researchers to perform enrichments for the same 1240 k SNPs. We used the Daicel Arbor Biosciences Prime Plus kit to enrich 132 ancient samples from three continents. We identified a systematic assay bias that increases genetic similarity between enriched samples and that cannot be explained by batch effects. We present the impact of the bias on population genetics inferences (e.g. Principal Components Analysis, ƒ-statistics) and genetic relatedness (READ). We compare the Prime Plus bias to that previously reported of the legacy 1240k enrichment assay. In ƒ-statistics, we find that all Prime-Plus-generated data exhibit artefactual excess shared drift, such that within-continent relationships cannot be correctly determined. The bias is more subtle in READ, though interpretation of the results can still be misleading in specific contexts. We expect the bias may affect analyses we have not yet tested. Our observations support previously reported concerns for the integration of different data types in paleogenomics. We also caution that technological solutions to generate 1240k data necessitate a thorough validation process before their adoption in the paleogenomic community.

The diurnal salivary glands transcriptome of Dermacentor nuttalli from the first four days of blood feeding

The ixodid tick Dermacentor nuttalli is distributed from southern Siberia to North China and is a vector of many pathogens. This species can have severe impacts on animal husbandry and human health. To date, the control of D. nuttalli is limited to the use of acaricides such as organophosphorus, synthetic pyrethroids and amidine pesticides. There are no environmentally friendly or reliable prevention and control measures, and little is known regarding key antigens involved in blood feeding. Salivary glands are major tissues involved in the blood feeding and pathogen transmission of ticks. Therefore, this study focused on salivary glands tissue to identify the dominant antigens of D. nuttalli involved in tick feeding. For this, high-throughput RNA sequencing (RNA-seq) was used for analysis. The transcriptome of female D. nuttalli ticks was assembled and characterized, and differentially expressed genes (DEGs) were identified in the salivary glands of ticks that had not fed (0 h) and of ticks after 24, 48, 72 and 96 h of feeding. There were 22,802,784, 22,275,013, 26,629,453, 24,982,389, and 22,596,230 high-quality clean reads obtained from salivary glands tissues at the five different blood feeding time points. The total number of annotated unigenes was 100,347. The differences in gene expression between different time points were compared, and functional enrichment was performed. Quantitative reverse transcription PCR (RT‒qPCR) was used to validate the RNA-seq results, the results of which showed that the differences in expressed transcripts presented similar trends. Among the identified DEGs, the most numerous were those with catalytic and binding activities and those involved in diverse metabolic pathways and cellular processes. The expression patterns of homologous and family-member proteins throughout the blood feeding period exhibited significant differences, strongly suggesting that the transcriptome composition is highly dynamic and likely subjected to important variation throughout the life cycle. Studies of gene sequences in D. nuttalli will greatly increase the information on tick protective antigens, which could potentially function as effective vaccine candidates or drug targets for the development of environmentally friendly acaricides.

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 .

Mining museums for historical DNA: advances and challenges in museomics

Historical DNA (hDNA), obtained from museum and herbarium specimens, has yielded spectacular new insights into the history of organisms. This includes documenting historical genetic erosion and extinction, discovering species new to science, resolving evolutionary relationships, investigating epigenetic effects, and determining origins of infectious diseases. However, the development of best-practices in isolating, processing, and analyzing hDNA remain under-explored, due to the substantial diversity of specimen preparation types, tissue sources, archival ages, and collecting histories. Thus, for hDNA to reach its full potential, and justify the destructive sampling of the rarest specimens, more experimental work using time-series collections, and the development of improved methods to correct for data asymmetries and biases due to DNA degradation are required.