HyRAD-X Exome Capture Museomics Unravels Giant Ground Beetle Evolution

Molecular phylogeny of mega-diverse Carabus attests late Miocene evolution of alpine environments in the Himalayan-Tibetan Orogen

The timing, sequence, and scale of uplift of the Himalayan-Tibetan Orogen (HTO) are controversially debated. Many geoscientific studies assume paleoelevations close to present-day elevations and the existence of alpine environments across the HTO already in the late Paleogene, contradicting fossil data. Using molecular genetic data of ground beetles, we aim to reconstruct the paleoenvironmental history of the HTO, focusing on its southern margin (Himalayas, South Tibet). Based on a comprehensive sampling of extratropical Carabus, and ~ 10,000 bp of mitochondrial and nuclear DNA we applied Bayesian and Maximum likelihood methods to infer the phylogenetic relationships. We show that Carabus arrived in the HTO at the Oligocene-Miocene boundary. During the early Miocene, five lineages diversified in different parts of the HTO, initially in its southern center and on its eastern margin. Evolution of alpine taxa occurred during the late Miocene. There were apparently no habitats for Carabus before the late Oligocene. Until the Late Oligocene elevations must have been low throughout the HTO. Temperate forests emerged in South Tibet in the late Oligocene at the earliest. Alpine environments developed in the HTO from the late Miocene and, in large scale, during the Pliocene-Quaternary. Findings are consistent with fossil records but contrast with uplift models recovered from stable isotope paleoaltimetry.

Dynamic evolution of size and colour in the highly specialized Zodarion ant-eating spiders

Ecological specialists constitute relevant case studies for understanding the mechanisms, potential and limitations of evolution. The species-rich and strictly myrmecophagous spiders of the genus Zodarion show diversified defence mechanisms, including myrmecomorphy of different ant species and nocturnality. Through Hybridization Capture Using RAD Probes (hyRAD), a phylogenomic technique designed for sequencing poorly preserved specimens, we reconstructed a phylogeny of Zodarion using 52 (approx. a third of the nominal) species that cover its phylogenetic and distributional diversity. We then estimated the evolution of body size and colour, traits that have diversified noticeably and are linked to defence mechanisms, across the group. Our genomic matrix of 300 loci led to a well-supported phylogenetic hypothesis that uncovered two main clades inside Zodarion. Ancestral state estimation revealed the highly dynamic evolution of body size and colour across the group, with multiple transitions and convergences in both traits, which we propose is likely indicative of multiple transitions in ant specialization across the genus. Our study will allow the informed targeted selection of Zodarion taxa of special interest for research into the group's remarkable adaptations to ant specialization. It also exemplifies the utility of hyRAD for phylogenetic studies using museum material.

The identity of Argyria lacteella (Fabricius, 1794) (Lepidoptera, Pyraloidea, Crambinae), synonyms, and related species revealed by morphology and DNA capture in type specimens

In this study the aim was to resolve the taxonomy of several species of Argyria Hübner (Pyraloidea, Crambinae) with previously unrecognised morphological variation. By analysing the DNA barcode (COI-5P) in numerous specimens, the aim was to reconstruct phylogenetic relationships between species, to provide better evidence for synonymies, and to circumscribe their geographical distribution. Using an innovative DNA hybridisation capture protocol, the DNA barcode of the lectotype of Argyria lacteella (Fabricius, 1794) was partially recovered for comparison with the 229 DNA barcode sequences of Argyria specimens available in the Barcode of Life Datasystems, and this firmly establishes the identity of the species. The same protocol was used for the following type specimens: the Argyria abronalis (Walker, 1859) holotype, thus confirming the synonymy of this name with A. lacteella, the holotype of A. lusella (Zeller, 1863), syn. rev., the holotype of A. multifacta Dyar, 1914, syn. nov. newly synonymised with A. lacteella, and a specimen of Argyria diplomochalis Dyar, 1913, collected in 1992. In addition, nine specimens of A. lacteella, A. diplomochalis, A. centrifugens Dyar, 1914 and A. gonogramma Dyar, 1915, from North to South America were sampled using classical COI amplification and Sanger sequencing. Argyria gonogramma Dyar, described from Bermuda, is the name to be applied to the more widespread North American species formerly identified as A. lacteella. Following morphological study of its holotype, Argyria vestalis Butler, 1878, syn. nov. is also synonymised with A. lacteella. The name A. pusillalis Hübner, 1818, is considered a nomen dubium associated with A. gonogramma. The adult morphology is diagnosed and illustrated, and distributions are plotted for A. lacteella, A. diplomochalis, A. centrifugens, and A. gonogramma based on slightly more than 800 specimens. For the first time, DNA barcode sequences are provided for the Antillean A. diplomochalis. This work provides a modified, improved protocol for the efficient hybrid capture enrichment of DNA barcodes from 18th and 19th century type specimens in order to solve taxonomic issues in Lepidoptera.

Predictors of sequence capture in a large-scale anchored phylogenomics project

Next-generation sequencing (NGS) technologies have revolutionized phylogenomics by decreasing the cost and time required to generate sequence data from multiple markers or whole genomes. Further, the fragmented DNA of biological specimens collected decades ago can be sequenced with NGS, reducing the need for collecting fresh specimens. Sequence capture, also known as anchored hybrid enrichment, is a method to produce reduced representation libraries for NGS sequencing. The technique uses single-stranded oligonucleotide probes that hybridize with pre-selected regions of the genome that are sequenced via NGS, culminating in a dataset of numerous orthologous loci from multiple taxa. Phylogenetic analyses using these sequences have the potential to resolve deep and shallow phylogenetic relationships. Identifying the factors that affect sequence capture success could save time, money, and valuable specimens that might be destructively sampled despite low likelihood of sequencing success. We investigated the impacts of specimen age, preservation method, and DNA concentration on sequence capture (number of captured sequences and sequence quality) while accounting for taxonomy and extracted tissue type in a large-scale butterfly phylogenomics project. This project used two probe sets to extract 391 loci or a subset of 13 loci from over 6,000 butterfly specimens. We found that sequence capture is a resilient method capable of amplifying loci in samples of varying age (0–111 years), preservation method (alcohol, papered, pinned), and DNA concentration (0.020 ng/μl - 316 ng/ul). Regression analyses demonstrate that sequence capture is positively correlated with DNA concentration. However, sequence capture and DNA concentration are negatively correlated with sample age and preservation method. Our findings suggest that sequence capture projects should prioritize the use of alcohol-preserved samples younger than 20 years old when available. In the absence of such specimens, dried samples of any age can yield sequence data, albeit with returns that diminish with increasing age.

Life-history traits drive spatial genetic structuring in Dinaric cave spiders

The subterranean ecosystem exerts strong selection pressures on the organisms that thrive in it. In response, obligate cave-dwellers have developed a series of morphological, physiological, and behavioral adaptations, such as eye reduction, appendage elongation, low metabolic rates or intermittent activity patterns, collectively referred to as troglomorphism. Traditionally, studies on cave organisms have been hampered by the difficulty of sampling (i.e., small population sizes, temporal heterogeneity in specimen occurrence, challenges imposed by the difficult-to-access nature of caves). Here, we circumvent this limitation by implementing a museomics approach. Specifically, we aim at comparing the genetic population structures of five cave spider species demonstrating contrasting life histories and levels of troglomorphism across different caves in the northern Dinarides (Balkans, Europe). We applied a genome-wide hybridization-capture approach (i.e., HyRAD) to capture DNA from 117 historical samples. By comparing the population genetic structures among five species and by studying isolation by distance, we identified deeper population structuring and more pronounced patterns of isolation by distance in the highly troglomorphic Parastalita stygia and Stalita pretneri ground dwellers, while the three web-building Troglohyphantes species, two of which can occasionally be found in surface habitats, showed less structured populations compatible with higher dispersal ability. The spatial distribution of genetic groups revealed common phylogeographic breaks among lineages across the studied species, which hint at the importance of environmental features in driving dispersal potential and shaping underground diversity.

Revised phylogeny from complete mitochondrial genomes of phyllostomid bats resolves subfamilial classification

Classically, molecular phylogenetic trees of Phyllostomidae have been inferred using a combination of a few mitochondrial and nuclear markers. However, there is still uncertainty in the relationships, especially among deep clades within the family. In this study, we provide newly sequenced complete mitochondrial genomes from 26 bat species, including genomes of 23 species reported here for the first time. By carefully analysing these genomes using maximum likelihood and Bayesian methods and different ingroup and outgroup samples, partition schemes and data types, we investigated the robustness and sensitivity of our phylogenetic results. The optimal topologies were those inferred from the complete data matrix of nucleotides, with complex and highly parameterized substitution models and partition schemes. Our results show a statistically robust picture of the evolutionary relationships between phyllostomid subfamilies and clarify hitherto uncertain relationships of Lonchorhininae and Macrotinae.

The untapped potential of macrofossils in ancient plant DNA research

The rapid development of ancient DNA analysis in the last decades has induced a paradigm shift in ecology and evolution. Driven by a combination of breakthroughs in DNA isolation techniques, high-throughput sequencing, and bioinformatics, ancient genome-scale data for a rapidly growing variety of taxa are now available, allowing researchers to directly observe demographic and evolutionary processes over time. However, the vast majority of paleogenomic studies still focus on human or animal remains. In this article, we make the case for a vast untapped resource of ancient plant material that is ideally suited for paleogenomic analyses: plant remains, such as needles, leaves, wood, seeds, or fruits, that are deposited in natural archives, such as lake sediments, permafrost, or even ice caves. Such plant remains are commonly found in large numbers and in stratigraphic sequence through time and have so far been used primarily to reconstruct past local species presences and abundances. However, they are also unique repositories of genetic information with the potential to revolutionize the fields of ecology and evolution by directly studying microevolutionary processes over time. Here, we give an overview of the current state-of-the-art, address important challenges, and highlight new research avenues to inspire future research.

Performance and automation of ancient DNA capture with RNA hyRAD probes

DNA hybridization-capture techniques allow researchers to focus their sequencing efforts on preselected genomic regions. This feature is especially useful when analysing ancient DNA (aDNA) extracts, which are often dominated by exogenous environmental sources. Here, we assessed, for the first time, the performance of hyRAD as an inexpensive and design-free alternative to commercial capture protocols to obtain authentic aDNA data from osseous remains. HyRAD relies on double enzymatic restriction of fresh DNA extracts to produce RNA probes that cover only a fraction of the genome and can serve as baits for capturing homologous fragments from aDNA libraries. We found that this approach could retrieve sequence data from horse remains coming from a range of preservation environments, including beyond radiocarbon range, yielding up to 146.5-fold on-target enrichment for aDNA extracts showing extremely low endogenous content (<1%). Performance was, however, more limited for those samples already characterized by good DNA preservation (>20%-30%), while the fraction of endogenous reads mapping on- and off-target was relatively insensitive to the original endogenous DNA content. Procedures based on two instead of a single round of capture increased on-target coverage up to 3.6-fold. Additionally, we used methylation-sensitive restriction enzymes to produce probes targeting hypomethylated regions, which improved data quality by reducing post-mortem DNA damage and mapping within multicopy regions. Finally, we developed a fully automated hyRAD protocol utilizing inexpensive robotic platforms to facilitate capture processing. Overall, our work establishes hyRAD as a cost-effective strategy to recover a set of shared orthologous variants across multiple ancient samples.