Aiming off the target: recycling target capture sequencing reads for investigating repetitive DNA

Maximising informativeness for target capture-based phylogenomics in Erica (Ericaceae)

Plant phylogenetics has been revolutionised in the genomic era, with target capture acting as the primary workhorse of most recent research in the new field of phylogenomics. Target capture (aka Hyb-Seq) allows researchers to sequence hundreds of genomic regions (loci) of their choosing, at relatively low cost per sample, from which to derive phylogenetically informative data. Although this highly flexible and widely applicable method has rightly earned its place as the field's de facto standard, it does not come without its challenges. In particular, users have to specify which loci to sequence-a surprisingly difficult task, especially when working with non-model groups, as it requires pre-existing genomic resources in the form of assembled genomes and/or transcriptomes. In the absence of taxon-specific genomic resources, target sets exist that are designed to work across broad taxonomic scales. However, the highly conserved loci that they target may lack informativeness for difficult phylogenetic problems, such as that presented by the rapid radiation of Erica in southern Africa. We designed a target set for Erica phylogenomics intended to maximise informativeness and minimise paralogy while maintaining universality by including genes from the widely used Angiosperms353 set. Comprising just over 300 genes, the targets had excellent recovery rates in roughly 90 Erica species as well as outgroups from Calluna, Daboecia, and Rhododendron, and had high information content as measured by parsimony informative sites and Quartet Internode Resolution Probability (QIRP) at shallow nodes. Notably, QIRP was positively correlated with intron content, while including introns in targets-rather than recovering them via exon-flanking "bycatch"-substantially improved intron recovery. Overall, our results show the value of building a custom target set, and we provide a suite of open-source tools that can be used to replicate our approach in other groups (https://github.com/SethMusker/TargetVet).

Holocentric repeat landscapes: From micro-evolutionary patterns to macro-evolutionary associations with karyotype evolution

Repetitive elements can cause large-scale chromosomal rearrangements, for example through ectopic recombination, potentially promoting reproductive isolation and speciation. Species with holocentric chromosomes, that lack a localized centromere, might be more likely to retain chromosomal rearrangements that lead to karyotype changes such as fusions and fissions. This is because chromosome segregation during cell division should be less affected than in organisms with a localized centromere. The relationships between repetitive elements and chromosomal rearrangements and how they may translate to patterns of speciation in holocentric organisms are though poorly understood. Here, we use a reference-free approach based on low-coverage short-read sequencing data to characterize the repeat landscape of two independently evolved holocentric groups: Erebia butterflies and Carex sedges. We consider both micro- and macro-evolutionary scales to investigate the repeat landscape differentiation between Erebia populations and the association between repeats and karyotype changes in a phylogenetic framework for both Erebia and Carex. At a micro-evolutionary scale, we found population differentiation in repeat landscape that increases with overall intraspecific genetic differentiation among four Erebia species. At a macro-evolutionary scale, we found indications for an association between repetitive elements and karyotype changes along both Erebia and Carex phylogenies. Altogether, our results suggest that repetitive elements are associated with the level of population differentiation and chromosomal rearrangements in holocentric clades and therefore likely play a role in adaptation and potentially species diversification.

The effect of missing data on evolutionary analysis of sequence capture bycatch, with application to an agricultural pest

Sequence capture is a genomic technique that selectively enriches target sequences before high throughput next-generation sequencing, to generate specific sequences of interest. Off-target or 'bycatch' data are often discarded from capture experiments, but can be leveraged to address evolutionary questions under some circumstances. Here, we investigated the effects of missing data on a variety of evolutionary analyses using bycatch from an exon capture experiment on the global pest moth, Helicoverpa armigera. We added > 200 new samples from across Australia in the form of mitogenomes obtained as bycatch from targeted sequence capture, and combined these into an additional larger dataset to total > 1000 mitochondrial cytochrome c oxidase subunit I (COI) sequences across the species' global distribution. Using discriminant analysis of principal components and Bayesian coalescent analyses, we showed that mitogenomes assembled from bycatch with up to 75% missing data were able to return evolutionary inferences consistent with higher coverage datasets and the broader literature surrounding H. armigera. For example, low-coverage sequences broadly supported the delineation of two H. armigera subspecies and also provided new insights into the potential for geographic turnover among these subspecies. However, we also identified key effects of dataset coverage and composition on our results. Thus, low-coverage bycatch data can offer valuable information for population genetic and phylodynamic analyses, but caution is required to ensure the reduced information does not introduce confounding factors, such as sampling biases, that drive inference. We encourage more researchers to consider maximizing the potential of the targeted sequence approach by examining evolutionary questions with their off-target bycatch where possible-especially in cases where no previous mitochondrial data exists-but recommend stratifying data at different genome coverage thresholds to separate sampling effects from genuine genomic signals, and to understand their implications for evolutionary research.

Oligo-barcode illuminates holocentric karyotype evolution in Rhynchospora (Cyperaceae)

Holocentric karyotypes are assumed to rapidly evolve through chromosome fusions and fissions due to the diffuse nature of their centromeres. Here, we took advantage of the recent availability of a chromosome-scale reference genome for Rhynchospora breviuscula, a model species of this holocentric genus, and developed the first set of oligo-based barcode probes for a holocentric plant. These probes were applied to 13 additional species of the genus, aiming to investigate the evolutionary dynamics driving the karyotype evolution in Rhynchospora. The two sets of probes were composed of 27,392 (green) and 23,968 (magenta) oligonucleotides (45-nt long), and generated 15 distinct FISH signals as a unique barcode pattern for the identification of all five chromosome pairs of the R. breviuscula karyotype. Oligo-FISH comparative analyzes revealed different types of rearrangements, such as fusions, fissions, putative inversions and translocations, as well as genomic duplications among the analyzed species. Two rounds of whole genome duplication (WGD) were demonstrated in R. pubera, but both analyzed accessions differed in the complex chain of events that gave rise to its large, structurally diploidized karyotypes with 2n = 10 or 12. Considering the phylogenetic relationships and divergence time of the species, the specificity and synteny of the probes were maintained up to species with a divergence time of ~25 My. However, karyotype divergence in more distant species hindered chromosome mapping and the inference of specific events. This barcoding system is a powerful tool to study chromosomal variations and genomic evolution in holocentric chromosomes of Rhynchospora species.

Repeat-based phylogenomics shed light on unclear relationships in the monocentric genus Juncus L. (Juncaceae)

The repetitive fraction (repeatome) of eukaryotic genomes is diverse and usually fast evolving, being an important tool for clarify plant systematics. The genus Juncus L. comprises 332 species, karyotypically recognized by having holocentric chromosomes. However, four species were recently described as monocentric, yet our understanding of their genome evolution is largely masked by unclear phylogenetic relationships. Here, we reassess the current Juncus systematics using low-coverage genome skimming data of 33 taxa to construct repeats, nuclear rDNA and plastome-based phylogenetic hypothesis. Furthermore, we characterize the repeatome and chromosomal distribution of Juncus-specific centromeric repeats/CENH3 protein to test the monocentricity reach in the genus. Repeat-base phylogenies revealed topologies congruent with the rDNA tree, but not with the plastome tree. The incongruence between nuclear and plastome chloroplast dataset suggest an ancient hybridization in the divergence of Juncotypus and Tenageia sections 40 Myr ago. The phylogenetic resolution at section level was better fitted with the rDNA/repeat-based approaches, with the recognition of two monophyletic sections (Stygiopsis and Tenageia). We found specific repeatome trends for the main lineages, such as the higher abundances of TEs in the Caespitosi and Iridifolii + Ozophyllum clades. CENH3 immunostaining confirmed the monocentricity of Juncus, which can be a generic synapomorphy for the genus. The heterogeneity of the repeatomes, with high phylogenetic informativeness, identified here may be correlated with their ancient origin (56 Mya) and reveals the potential of comparative genomic analyses for understanding plant systematics and evolution.

Target capture and genome skimming for plant diversity studies

Recent technological advances in long-read high-throughput sequencing and assembly methods have facilitated the generation of annotated chromosome-scale whole-genome sequence data for evolutionary studies; however, generating such data can still be difficult for many plant species. For example, obtaining high-molecular-weight DNA is typically impossible for samples in historical herbarium collections, which often have degraded DNA. The need to fast-freeze newly collected living samples to conserve high-quality DNA can be complicated when plants are only found in remote areas. Therefore, short-read reduced-genome representations, such as target capture and genome skimming, remain important for evolutionary studies. Here, we review the pros and cons of each technique for non-model plant taxa. We provide guidance related to logistics, budget, the genomic resources previously available for the target clade, and the nature of the study. Furthermore, we assess the available bioinformatic analyses, detailing best practices and pitfalls, and suggest pathways to combine newly generated data with legacy data. Finally, we explore the possible downstream analyses allowed by the type of data generated using each technique. We provide a practical guide to help researchers make the best-informed choice regarding reduced genome representation for evolutionary studies of non-model plants in cases where whole-genome sequencing remains impractical.

Investigating the diversification of holocentromeric satellite DNA Tyba in Rhynchospora (Cyperaceae)

BACKGROUND AND AIMS

Satellite DNAs (satDNAs) are repetitive sequences composed by tandemly arranged, often highly homogenized units called monomers. Although satDNAs are usually fast evolving, some satDNA families can be conserved across species separated by several millions of years, probably because of their functional roles in the genomes. Tyba was the first centromere-specific satDNA described for a holocentric organism, until now being characterized for only eight species of the genus Rhynchospora Vahl. (Cyperaceae). Here, we characterized Tyba across a broad sampling of the genus, analysing and comparing its evolutionary patterns with other satDNAs.

METHODS

We characterized the structure and sequence evolution of satDNAs across a robust dadated phylogeny based on Hybrid Target-Capture Sequencing (hyb-seq) of 70 species. We mined the repetitive fraction for Tyba-like satellites to compare its features with other satDNAs and to construct a Tyba-based phylogeny for the genus.

KEY RESULTS

Our results show that Tyba is present in the majority of examined species of the genus, spanning four of the five major clades and maintaining intrafamily pairwise identity of 70.9% over 31 Myr. In comparison, other satellite families presented higher intrafamily pairwise identity but are phylogenetically restricted. Furthermore, Tyba sequences could be divided into 12 variants grouped into three different clade-specific subfamilies, showing evidence of traditional models of satDNA evolution, such as the concerted evolution and library models. Besides, a Tyba-based phylogeny showed high congruence with the hyb-seq topology. Our results show structural indications of a possible relationship of Tyba with nucleosomes, given its high curvature peaks over conserved regions and overall high bendability values compared with other non-centromeric satellites.

CONCLUSIONS

Overall, Tyba shows a remarkable sequence conservation and phylogenetic significance across the genus Rhynchospora, which suggests that functional roles might lead to long-term stability and conservation for satDNAs in the genome.

Fishing for DNA? Designing baits for population genetics in target enrichment experiments: Guidelines, considerations and the new tool supeRbaits

Targeted sequencing is an increasingly popular next-generation sequencing (NGS) approach for studying populations that involves focusing sequencing efforts on specific parts of the genome of a species of interest. Methodologies and tools for designing targeted baits are scarce but in high demand. Here, we present specific guidelines and considerations for designing capture sequencing experiments for population genetics for both neutral genomic regions and regions subject to selection. We describe the bait design process for three diverse fish species: Atlantic salmon, Atlantic cod and tiger shark, which was carried out in our research group, and provide an evaluation of the performance of our approach across both historical and modern samples. The workflow used for designing these three bait sets has been implemented in the R-package supeRbaits, which encompasses our considerations and guidelines for bait design for the benefit of researchers and practitioners. The supeRbaits R-package is user-friendly and versatile. It is written in C++ and implemented in R. supeRbaits and its manual are available from Github: https://github.com/BelenJM/supeRbaits.

Larix species range dynamics in Siberia since the Last Glacial captured from sedimentary ancient DNA

Climate change is expected to cause major shifts in boreal forests which are in vast areas of Siberia dominated by two species of the deciduous needle tree larch (Larix). The species differ markedly in their ecosystem functions, thus shifts in their respective ranges are of global relevance. However, drivers of species distribution are not well understood, in part because paleoecological data at species level are lacking. This study tracks Larix species distribution in time and space using target enrichment on sedimentary ancient DNA extracts from eight lakes across Siberia. We discovered that Larix sibirica, presently dominating in western Siberia, likely migrated to its northern distribution area only in the Holocene at around 10,000 years before present (ka BP), and had a much wider eastern distribution around 33 ka BP. Samples dated to the Last Glacial Maximum (around 21 ka BP), consistently show genotypes of L. gmelinii. Our results suggest climate as a strong determinant of species distribution in Larix and provide temporal and spatial data for species projection in a changing climate.

Using ancient sedimentary DNA from up to 50 kya, dramatic distributional shifts are documented in two dominant boreal larch species, likely guided by environmental changes suggesting climate as a strong determinant of species distribution.

Can we have it all? Repurposing target capture for repeat genomics. A commentary on: 'Aiming off the target: recycling target capture sequencing reads for investigating repetitive DNA'

This article comments on:

Lucas Costa, André Marques, Chris Buddenhagen, William Wayt Thomas, Bruno Huettel, Veit Schubert, Steven Dodsworth, Andreas Houben, Gustavo Souza and Andrea Pedrosa-Harand, Aiming off the target: recycling target capture sequencing reads for investigating repetitive DNA, Annals of Botany, Volume 128, Issue 7, 2 December 2021, Pages 835–848 https://doi.org/10.1093/aob/mcab063