Sequence capture by hybridization to explore modern and ancient genomic diversity in model and nonmodel organisms

Targeting FGFR for cancer therapy

The FGFR signaling pathway is integral to cellular activities, including proliferation, differentiation, and survival. Dysregulation of this pathway is implicated in numerous human cancers, positioning FGFR as a prominent therapeutic target. Here, we conduct a comprehensive review of the function, signaling pathways and abnormal alterations of FGFR, as well as its role in tumorigenesis and development. Additionally, we provide an in-depth analysis of pivotal phase 2 and 3 clinical trials evaluating the performance and safety of FGFR inhibitors in oncology, thereby shedding light on the current state of clinical research in this field. Then, we highlight four drugs that have been approved for marketing by the FDA, offering insights into their molecular mechanisms and clinical achievements. Our discussion encompasses the intricate landscape of FGFR-driven tumorigenesis, current techniques for pinpointing FGFR anomalies, and clinical experiences with FGFR inhibitor regimens. Furthermore, we discuss the inherent challenges of targeting the FGFR pathway, encompassing resistance mechanisms such as activation by gatekeeper mutations, alternative pathways, and potential adverse reactions. By synthesizing the current evidence, we underscore the potential of FGFR-centric therapies to enhance patient prognosis, while emphasizing the imperative need for continued research to surmount resistance and optimize treatment modalities.

Computational deconvolution of DNA methylation data from mixed DNA samples

In this review, we provide a comprehensive overview of the different computational tools that have been published for the deconvolution of bulk DNA methylation (DNAm) data. Here, deconvolution refers to the estimation of cell-type proportions that constitute a mixed sample. The paper reviews and compares 25 deconvolution methods (supervised, unsupervised or hybrid) developed between 2012 and 2023 and compares the strengths and limitations of each approach. Moreover, in this study, we describe the impact of the platform used for the generation of methylation data (including microarrays and sequencing), the applied data pre-processing steps and the used reference dataset on the deconvolution performance. Next to reference-based methods, we also examine methods that require only partial reference datasets or require no reference set at all. In this review, we provide guidelines for the use of specific methods dependent on the DNA methylation data type and data availability.

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.

Applying targeted gene hybridization capture to viruses with a focus to SARS-CoV-2

Although next-generation sequencing technologies are advancing rapidly, many research topics often require selective sequencing of genomic regions of interest. In addition, sequencing low-titre viruses is challenging, especially for coronaviruses, which are the largest RNA viruses. Prior to sequencing, enrichment of viral particles can help to significantly increase target sequence information as well as avoid large sequencing efforts and, consequently, can increase sensitivity and reduce sequencing costs. Targeting nucleic acids using capture by hybridization is another efficient method that can be performed by applying complementary probes (DNA or RNA baits) to directly enrich genetic information of interest while removing background non-target material. In studies where sequence capture by hybridization has been applied to the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus, most authors agree that this technique is useful to easily access sequence targets in complex samples. Furthermore, this approach allows for complete or near-complete sequencing of the viral genome, even in samples with low viral load or poor nucleic acid integrity. In addition, this strategy is highly efficient at discovering new variants by facilitating downstream investigations, such as phylogenetics, epidemiology, and evolution. Commercial kits, as well as in-house protocols, have been developed for enrichment of viral sequences. However, these kits have multiple variations in procedure, with differences in performance. This review compiles and describes studies in which hybridization capture has been applied to SARS-CoV-2 variant genomes.

Applying molecular and genetic methods to trees and their fungal communities

Forests provide invaluable economic, ecological, and social services. At the same time, they are exposed to several threats, such as fragmentation, changing climatic conditions, or increasingly destructive pests and pathogens. Trees, the inherent species of forests, cannot be viewed as isolated organisms. Manifold (micro)organisms are associated with trees playing a pivotal role in forest ecosystems. Of these organisms, fungi may have the greatest impact on the life of trees. A multitude of molecular and genetic methods are now available to investigate tree species and their associated organisms. Due to their smaller genome sizes compared to tree species, whole genomes of different fungi are routinely compared. Such studies have only recently started in forest tree species. Here, we summarize the application of molecular and genetic methods in forest conservation genetics, tree breeding, and association genetics as well as for the investigation of fungal communities and their interrelated ecological functions. These techniques provide valuable insights into the molecular basis of adaptive traits, the impacts of forest management, and changing environmental conditions on tree species and fungal communities and can enhance tree-breeding cycles due to reduced time for field testing. It becomes clear that there are multifaceted interactions among microbial species as well as between these organisms and trees. We demonstrate the versatility of the different approaches based on case studies on trees and fungi. KEY POINTS: • Current knowledge of genetic methods applied to forest trees and associated fungi. • Genomic methods are essential in conservation, breeding, management, and research. • Important role of phytobiomes for trees and their ecosystems.

Myeloma immunoglobulin rearrangement and translocation detection through targeted capture sequencing

Multiple myeloma is a plasma cell neoplasm characterized by clonal immunoglobulin V(D)J signatures and oncogenic immunoglobulin gene translocations. Additional subclonal genomic changes are acquired with myeloma progression and therapeutic selection. PCR-based methods to detect V(D)J rearrangements can have biases introduced by highly multiplexed reactions and primers undermined by somatic hypermutation, and are not readily extended to include mutation detection. Here, we report a hybrid-capture approach (CapIG-seq) targeting the 3' and 5' ends of the V and J segments of all immunoglobulin loci that enable the efficient detection of V(D)J rearrangements. We also included baits for oncogenic translocations and mutation detection. We demonstrate complete concordance with matched whole-genome sequencing and/or PCR clonotyping of 24 cell lines and report the clonal sequences for 41 uncharacterized cell lines. We also demonstrate the application to patient specimens, including 29 bone marrow and 39 cell-free DNA samples. CapIG-seq shows concordance between bone marrow and cfDNA blood samples (both contemporaneous and follow-up) with regard to the somatic variant, V(D)J, and translocation detection. CapIG-seq is a novel, efficient approach to examining genomic alterations in myeloma.

Clinically relevant fusion oncogenes: detection and practical implications

Mechanistically, chimeric genes result from DNA rearrangements and include parts of preexisting normal genes combined at the genomic junction site. Some rearranged genes encode pathological proteins with altered molecular functions. Those which can aberrantly promote carcinogenesis are called fusion oncogenes. Their formation is not a rare event in human cancers, and many of them were documented in numerous study reports and in specific databases. They may have various molecular peculiarities like increased stability of an oncogenic part, self-activation of tyrosine kinase receptor moiety, and altered transcriptional regulation activities. Currently, tens of low molecular mass inhibitors are approved in cancers as the drugs targeting receptor tyrosine kinase (RTK) oncogenic fusion proteins, that is, including ALK, ABL, EGFR, FGFR1-3, NTRK1-3, MET, RET, ROS1 moieties. Therein, the presence of the respective RTK fusion in the cancer genome is the diagnostic biomarker for drug prescription. However, identification of such fusion oncogenes is challenging as the breakpoint may arise in multiple sites within the gene, and the exact fusion partner is generally unknown. There is no gold standard method for RTK fusion detection, and many alternative experimental techniques are employed nowadays to solve this issue. Among them, RNA-seq-based methods offer an advantage of unbiased high-throughput analysis of only transcribed RTK fusion genes, and of simultaneous finding both fusion partners in a single RNA-seq read. Here we focus on current knowledge of biology and clinical aspects of RTK fusion genes, related databases, and laboratory detection methods.

Targeted 16S rRNA Gene Capture by Hybridization and Bioinformatic Analysis

Next-generation sequencing technologies have impressively unlocked capacities to depict the complexity of microbial communities. Microbial community structure is for now routinely monitored by sequencing of 16S rRNA gene, a phylogenetic marker almost conserved among bacteria and archaea. Nevertheless, amplicon sequencing, the most popular used approach, suffers from several biases impacting the picture of microbial communities. Here, we describe an innovative method based on gene capture by hybridization for the targeted enrichment of 16S rDNA biomarker from metagenomic samples. Coupled to near full-length 16S rDNA reconstruction, this approach enables an exhaustive and accurate description of microbial communities by enhancing taxonomic and phylogenetic resolutions. Furthermore, access of captured 16S flanking regions opens link between structure and function in microbial communities.

RiboTaxa: combined approaches for rRNA genes taxonomic resolution down to the species level from metagenomics data revealing novelties

Metagenomic classifiers are widely used for the taxonomic profiling of metagenomics data and estimation of taxa relative abundance. Small subunit rRNA genes are a gold standard for phylogenetic resolution of microbiota, although the power of this marker comes down to its use as full-length. We aimed at identifying the tools that can efficiently lead to taxonomic resolution down to the species level. To reach this goal, we benchmarked the performance and accuracy of rRNA-specialized versus general-purpose read mappers, reference-targeted assemblers and taxonomic classifiers. We then compiled the best tools (BBTools, FastQC, SortMeRNA, MetaRib, EMIRGE, VSEARCH, BBMap and QIIME 2’s Sklearn classifier) to build a pipeline called RiboTaxa. Using metagenomics datasets, RiboTaxa gave the best results compared to other tools (i.e. Kraken2, Centrifuge, METAXA2, phyloFlash, SPINGO, BLCA, MEGAN) with precise taxonomic identification and relative abundance description without false positive detection (F-measure of 100% and 83.7% at genus level and species level, respectively). Using real datasets from various environments (i.e. ocean, soil, human gut) and from different approaches (e.g. metagenomics and gene capture by hybridization), RiboTaxa revealed microbial novelties not discerned by current bioinformatics analysis opening new biological perspectives in human and environmental health.

Conservation genomics and systematics of a near-extinct island radiation

Conservation benefits from incorporating genomics to explore the impacts of population declines, inbreeding, loss of genetic variation and hybridization. Here we use the near-extinct Mariana Islands reedwarbler radiation to showcase how ancient DNA approaches can allow insights into the population dynamics of extinct species and threatened populations for which historical museum specimens or material with low DNA yield (e.g., scats, feathers) are the only sources for DNA. Despite their having paraphyletic mtDNA, nuclear SNPs support the distinctiveness of critically endangered Acrocephalus hiwae and the other three species in the radiation that went extinct between the 1960s and 1990s. Two extinct species, A. yamashinae and A. luscinius, were deeply divergent from each other and from a third less differentiated lineage containing A. hiwae and extinct A. nijoi. Both mtDNA and SNPs suggest that the two isolated populations of A. hiwae from Saipan and Alamagan Islands are sufficiently distinct to warrant subspecies recognition and separate conservation management. We detected no significant differences in genetic diversity or inbreeding between Saipan and Alamagan, nor strong signatures of geographic structuring within either island. However, the implications of possible signatures of inbreeding in both Saipan and Alamagan, and long-term population declines in A. hiwae that predate modern anthropogenic threats require further study with denser population sampling. Our study highlights the value conservation genomics studies of island radiations have as windows onto the possible future for the world's biota as climate change and habitat destruction increasingly fragments their ranges and contributes to rapid declines in population abundances.

Revealing microbial species diversity using sequence capture by hybridization

Targeting small parts of the 16S rDNA phylogenetic marker by metabarcoding reveals microorganisms of interest but cannot achieve a taxonomic resolution at the species level, precluding further precise characterizations. To identify species behind operational taxonomic units (OTUs) of interest, even in the rare biosphere, we developed an innovative strategy using gene capture by hybridization. From three OTU sequences detected upon polyphenol supplementation and belonging to the rare biosphere of the human gut microbiota, we revealed 59 nearly full-length 16S rRNA genes, highlighting high bacterial diversity hidden behind OTUs while evidencing novel taxa. Inside each OTU, revealed 16S rDNA sequences could be highly distant from each other with similarities down to 85 %. We identified one new family belonging to the order Clostridiales, 39 new genera and 52 novel species. Related bacteria potentially involved in polyphenol degradation have also been identified through genome mining and our results suggest that the human gut microbiota could be much more diverse than previously thought.

Capture by hybridization for full-length barcode-based eukaryotic and prokaryotic biodiversity inventories of deep sea ecosystems

Biodiversity inventory of marine systems remains limited due to unbalanced access to the three ocean dimensions. The use of environmental DNA (eDNA) for metabarcoding allows fast and effective biodiversity inventory and is forecast as a future biodiversity research and biomonitoring tool. However, in poorly understood ecosystems, eDNA results remain difficult to interpret due to large gaps in reference databases and PCR bias limiting the detection of some major phyla. Here, we aimed to circumvent these limitations by avoiding PCR and recollecting larger DNA fragments to improve assignment of detected taxa through phylogenetic reconstruction. We applied capture by hybridization (CBH) to enrich DNA from deep-sea sediment samples and compared the results with those obtained through an up-to-date metabarcoding PCR-based approach (MTB). Originally developed for bacterial communities and targeting 16S rDNA, the CBH approach was applied to 18S rDNA to improve the detection of species forming benthic communities of eukaryotes, with a particular focus on metazoans. The results confirmed the possibility of extending CBH to metazoans with two major advantages: (i) CBH revealed a broader spectrum of prokaryotic, eukaryotic, and particularly metazoan diversity, and (ii) CBH allowed much more robust phylogenetic reconstructions of full-length barcodes with up to 1900 base pairs. This is particularly important for taxa whose assignment is hampered by gaps in reference databases. This study provides a database and probes to apply 18S CBH to diverse marine systems, confirming this promising new tool to improve biodiversity assessments in data-poor ecosystems such as those in the deep sea.

Enrichment of low abundance DNA/RNA by oligonucleotide-clicked iron oxide nanoparticles

Detection of low abundance target DNA/RNA for clinical or research purposes is challenging because the target sequences can be hidden under a large background of human genomic or non-human metagenomic sequences. We describe a probe-based capture method to enrich for target sequences with DNA-clicked iron oxide nanoparticles. Our method was tested against commercial capture assays using streptavidin beads, on a set of probes derived from a common genotype of the hepatitis C virus. We showed that our method is more specific and sensitive, most likely due to the combination of an inert silica coating and a high density of DNA probes clicked to the nanoparticles. This facilitates target capture below the limits of detection for TaqMan qPCR, and we believe that this method has the potential to transform management of infectious diseases.

Phylogenomic Resolution of the Cetacean Tree of Life Using Target Sequence Capture

The evolution of cetaceans, from their early transition to an aquatic lifestyle to their subsequent diversification, has been the subject of numerous studies. However, although the higher-level relationships among cetacean families have been largely settled, several aspects of the systematics within these groups remain unresolved. Problematic clades include the oceanic dolphins (37 spp.), which have experienced a recent rapid radiation, and the beaked whales (22 spp.), which have not been investigated in detail using nuclear loci. The combined application of high-throughput sequencing with techniques that target specific genomic sequences provide a powerful means of rapidly generating large volumes of orthologous sequence data for use in phylogenomic studies. To elucidate the phylogenetic relationships within the Cetacea, we combined sequence capture with Illumina sequencing to generate data for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\sim $\end{document}3200 protein-coding genes for 68 cetacean species and their close relatives including the pygmy hippopotamus. By combining data from \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$>$\end{document}38,000 exons with existing sequences from 11 cetaceans and seven outgroup taxa, we produced the first comprehensive comparative genomic data set for cetaceans, spanning 6,527,596 aligned base pairs (bp) and 89 taxa. Phylogenetic trees reconstructed with maximum likelihood and Bayesian inference of concatenated loci, as well as with coalescence analyses of individual gene trees, produced mostly concordant and well-supported trees. Our results completely resolve the relationships among beaked whales as well as the contentious relationships among oceanic dolphins, especially the problematic subfamily Delphinidae. We carried out Bayesian estimation of species divergence times using MCMCTree and compared our complete data set to a subset of clocklike genes. Analyses using the complete data set consistently showed less variance in divergence times than the reduced data set. In addition, integration of new fossils (e.g., Mystacodon selenensis) indicates that the diversification of Crown Cetacea began before the Late Eocene and the divergence of Crown Delphinidae as early as the Middle Miocene. [Cetaceans; phylogenomics; Delphinidae; Ziphiidae; dolphins; whales.]

Assessing Genome-Wide Diversity in European Hantaviruses through Sequence Capture from Natural Host Samples

Research on the ecology and evolution of viruses is often hampered by the limitation of sequence information to short parts of the genomes or single genomes derived from cultures. In this study, we use hybrid sequence capture enrichment in combination with high-throughput sequencing to provide efficient access to full genomes of European hantaviruses from rodent samples obtained in the field. We applied this methodology to Tula (TULV) and Puumala (PUUV) orthohantaviruses for which analyses from natural host samples are typically restricted to partial sequences of their tri-segmented RNA genome. We assembled a total of ten novel hantavirus genomes de novo with very high coverage (on average >99%) and sequencing depth (average >247×). A comparison with partial Sanger sequences indicated an accuracy of >99.9% for the assemblies. An analysis of two common vole (Microtus arvalis) samples infected with two TULV strains each allowed for the de novo assembly of all four TULV genomes. Combining the novel sequences with all available TULV and PUUV genomes revealed very similar patterns of sequence diversity along the genomes, except for remarkably higher diversity in the non-coding region of the S-segment in PUUV. The genomic distribution of polymorphisms in the coding sequence was similar between the species, but differed between the segments with the highest sequence divergence of 0.274 for the M-segment, 0.265 for the S-segment, and 0.248 for the L-segment (overall 0.258). Phylogenetic analyses showed the clustering of genome sequences consistent with their geographic distribution within each species. Genome-wide data yielded extremely high node support values, despite the impact of strong mutational saturation that is expected for hantavirus sequences obtained over large spatial distances. We conclude that genome sequencing based on capture enrichment protocols provides an efficient means for ecological and evolutionary investigations of hantaviruses at an unprecedented completeness and depth.

A Guide to Carrying Out a Phylogenomic Target Sequence Capture Project

High-throughput DNA sequencing techniques enable time- and cost-effective sequencing of large portions of the genome. Instead of sequencing and annotating whole genomes, many phylogenetic studies focus sequencing effort on large sets of pre-selected loci, which further reduces costs and bioinformatic challenges while increasing coverage. One common approach that enriches loci before sequencing is often referred to as target sequence capture. This technique has been shown to be applicable to phylogenetic studies of greatly varying evolutionary depth. Moreover, it has proven to produce powerful, large multi-locus DNA sequence datasets suitable for phylogenetic analyses. However, target capture requires careful considerations, which may greatly affect the success of experiments. Here we provide a simple flowchart for designing phylogenomic target capture experiments. We discuss necessary decisions from the identification of target loci to the final bioinformatic processing of sequence data. We outline challenges and solutions related to the taxonomic scope, sample quality, and available genomic resources of target capture projects. We hope this review will serve as a useful roadmap for designing and carrying out successful phylogenetic target capture studies.

Assessing changes in genomic divergence following a century of human-mediated secondary contact among wild and captive-bred ducks

Along with manipulating habitat, the direct release of domesticated individuals into the wild is a practice used worldwide to augment wildlife populations. We test between possible outcomes of human-mediated secondary contact using genomic techniques at both historical and contemporary timescales for two iconic duck species. First, we sequence several thousand ddRAD-seq loci for contemporary mallards (Anas platyrhynchos) throughout North America and two domestic mallard types (i.e., known game-farm mallards and feral Khaki Campbell's). We show that North American mallards may well be becoming a hybrid swarm due to interbreeding with domesticated game-farm mallards released for hunting. Next, to attain a historical perspective, we applied a bait-capture array targeting thousands of loci in century-old (1842-1915) and contemporary (2009-2010) mallard and American black duck (Anas rubripes) specimens. We conclude that American black ducks and mallards have always been closely related, with a divergence time of ~600,000 years before present, and likely evolved through prolonged isolation followed by limited bouts of gene flow (i.e., secondary contact). They continue to maintain genetic separation, a finding that overturns decades of prior research and speculation suggesting the genetic extinction of the American black duck due to contemporary interbreeding with mallards. Thus, despite having high rates of hybridization, actual gene flow is limited between mallards and American black ducks. Conversely, our historical and contemporary data confirm that the intensive stocking of game-farm mallards during the last ~100 years has fundamentally changed the genetic integrity of North America's wild mallard population, especially in the east. It thus becomes of great interest to ask whether the iconic North American mallard is declining in the wild due to introgression of maladaptive traits from domesticated forms. Moreover, we hypothesize that differential gene flow from domestic game-farm mallards into the wild mallard population may explain the overall temporal increase in differentiation between wild black ducks and mallards, as well as the uncoupling of genetic diversity and effective population size estimates across time in our results. Finally, our findings highlight how genomic methods can recover complex population histories by capturing DNA preserved in traditional museum specimens.

Amand P, Le HQ, Su Z, Bai G. Multiplex restriction amplicon sequencing: a novel next-generation sequencing-based marker platform for high-throughput genotyping

To enable rapid selection of traits in marker-assisted breeding, markers must be technically simple, low-cost, high-throughput and randomly distributed in a genome. We developed such a technology, designated as Multiplex Restriction Amplicon Sequencing (MRASeq), which reduces genome complexity by polymerase chain reaction (PCR) amplification of amplicons flanked by restriction sites. The first PCR primers contain restriction site sequences at 3'-ends, preceded by 6-10 bases of specific or degenerate nucleotide sequences and then by a unique M13-tail sequence which serves as a binding site for a second PCR that adds sequencing primers and barcodes to allow sample multiplexing for sequencing. The sequences of restriction sites and adjacent nucleotides can be altered to suit different species. Physical mapping of MRASeq SNPs from a biparental population of allohexaploid wheat (Triticum aestivum L.) showed a random distribution of SNPs across the genome. MRASeq generated thousands of SNPs from a wheat biparental population and natural populations of wheat and barley (Hordeum vulgare L.). This novel, next-generation sequencing-based genotyping platform can be used for linkage mapping to screen quantitative trait loci (QTL), background selection in breeding and many other genetics and breeding applications of various species.

A nuclease-polymerase chain reaction enables amplification of probes used for capture-based DNA target enrichment

DNA target enrichment via hybridization capture is a commonly adopted approach which remains expensive due in-part to using biotinylated-probe panels. Here we provide a novel isothermal amplification reaction to amplify rapidly existing probe panels without knowledge of the sequences involved, thereby decreasing a major portion of the overall sample preparation cost. The reaction employs two thermostable enzymes, BST-polymerase and duplex-specific nuclease DSN. DSN initiates random ‘nicks’ on double-stranded-DNA which enable BST to polymerize DNA by displacing the nicked-strand. Displaced strands re-hybridize and the process leads to an exponential chain-reaction generating biotinylated DNA fragments within minutes. When starting from single-stranded-DNA, DNA is first converted to double-stranded-DNA via terminal-deoxynucleotidyl-transferase (TdT) prior to initiation of BST–DSN reaction. Biotinylated probes generated by TdT–BST–DSN (TBD) reactions using panels of 33, 190 or 7186 DNA targets are used for hybrid-capture-based target enrichment from amplified circulating-DNA, followed by targeted re-sequencing. Polymerase-nuclease isothermal-chain-reactions generate random amplified probes with no apparent sequence dependence. One round of target-capture using TBD probes generates a modest on-target sequencing ratio, while two successive rounds of capture generate >80% on-target reads with good sequencing uniformity. TBD-reactions generate enough capture-probes to increase by approximately two to three orders-of-magnitude the target-enrichment experiments possible from an initial set of probes.

Quaternary DNA: A Multidisciplinary Research Field

The purpose of this Milankovitch review is to explain the significance of Quaternary DNA studies and the importance of the recent methodological advances that have enabled the study of late Quaternary remains in more detail, and the testing of new assumptions in evolutionary biology and phylogeography to reconstruct the past. The topic is wide, and this review is not intended to be an exhaustive account of all the aDNA work performed in the last three decades on late-Quaternary remains. Instead, it is a selection of relevant studies aimed at illustrating how aDNA has been used to reconstruct not only environments of the past, but also the history of many species including our own.

Off-target capture data, endosymbiont genes and morphology reveal a relict lineage that is sister to all other singing cicadas

Phylogenetic asymmetry is common throughout the tree of life and results from contrasting patterns of speciation and extinction in the paired descendant lineages of ancestral nodes. On the depauperate side of a node, we find extant ‘relict’ taxa that sit atop long, unbranched lineages. Here, we show that a tiny, pale green, inconspicuous and poorly known cicada in the genus Derotettix, endemic to degraded salt-plain habitats in arid regions of central Argentina, is a relict lineage that is sister to all other modern cicadas. Nuclear and mitochondrial phylogenies of cicadas inferred from probe-based genomic hybrid capture data of both target and non-target loci and a morphological cladogram support this hypothesis. We strengthen this conclusion with genomic data from one of the cicada nutritional bacterial endosymbionts, Sulcia, an ancient and obligate endosymbiont of the larger plant-sucking bugs (Auchenorrhyncha) and an important source of maternally inherited phylogenetic data. We establish Derotettiginae subfam. nov. as a new, monogeneric, fifth cicada subfamily, and compile existing and new data on the distribution, ecology and diet of Derotettix. Our consideration of the palaeoenvironmental literature and host-plant phylogenetics allows us to predict what might have led to the relict status of Derotettix over 100 Myr of habitat change in South America.

Terrestrial mammal surveillance using hybridization capture of environmental DNA from African waterholes

Determining species distributions can be extremely challenging but is crucial to ecological and conservation research. Environmental DNA (eDNA) approaches have shown particular promise in aquatic systems for several vertebrate and invertebrate species. For terrestrial animals, however, eDNA-based surveys are considerably more difficult due to the lack of or difficulty in obtaining appropriate sampling substrate. In water-limited ecosystem where terrestrial mammals are often forced to congregate at waterholes, water and sediment from shared water sources may be a suitable substrate for noninvasive eDNA approaches. We characterized mitochondrial DNA sequences from a broad range of terrestrial mammal species in two different African ecosystems (in Namibia and Tanzania) using eDNA isolated from native water, sediment and water filtered through glass fibre filters. A hybridization capture enrichment with RNA probes targeting the mitochondrial genomes of 38 mammal species representing the genera/families expected at the respective ecosystems was employed, and 16 species were identified, with a maximum mitogenome coverage of 99.8%. Conventional genus-specific PCRs were tested on environmental samples for two genera producing fewer positive results than hybridization capture enrichment. An experiment with mock samples using DNA from non-African mammals showed that baits covering 30% of nontarget mitogenomes produced 91% mitogenome coverage after capture. In the mock samples, over-representation of DNA of one species still allowed for the detection of DNA of other species that was at a 100-fold lower concentration. Hybridization capture enrichment of eDNA is therefore an effective method for monitoring terrestrial mammal species from shared water sources.

Nonfunctionalized PNAs as Beacons for Nucleic Acid Detection in a Nanopore System

In this work, single-channel current recordings were used to selectively detect individual ssDNA strands in the vestibule of the α-hemolysin (α-HL) protein nanopore. The sensing mechanism was based on the detection of the intrinsic topological change of target ssDNA molecules after the hybridization with complementary PNA fragments. The readily distinguishable current signatures of PNA-DNA duplexes reversible association with the α-HL's vestibule, in terms of blockade amplitudes and kinetic features, allows specific detection of nucleic acid hybridization.

Optimized DNA extraction and library preparation for minute arthropods: Application to target enrichment in chalcid wasps used for biocontrol

Target enrichment is increasingly used for genotyping of plant and animal species or to better understand the evolutionary history of important lineages through the inference of statistically robust phylogenies. Limitations to routine target enrichment are both the complexity of current protocols and low input DNA quantity. Thus, working with tiny organisms such as microarthropods can be challenging. Here, we propose easy to set up optimizations for DNA extraction and library preparation prior to target enrichment. Prepared libraries were used to capture 1,432 ultraconserved elements (UCEs) from microhymenoptera (Chalcidoidea), which are among the tiniest insects on Earth and the most commercialized worldwide for biological control purposes. Results show no correlation between input DNA quantities (1.8-250 ng, 0.4 ng with an extra whole genome amplification step) and the number of sequenced UCEs on an Illumina MiSeq. Phylogenetic inferences highlight the potential of UCEs to solve relationships within the families of chalcid wasps, which has not been achieved so far. The protocol (library preparation + target enrichment) allows processing 96 specimens in five working days, by a single person, without requiring the use of expensive robotic molecular biology platforms, which could help to generalize the use of target enrichment for minute specimens.

Cinacalcet therapy in an infant with an R185Q calcium-sensing receptor mutation causing hyperparathyroidism: a case report and review of the literature

Background Neonatal severe hyperparathyroidism (NSHPT) is commonly treated with either parathyroidectomy or pharmacologic agents with varying efficacy and numerous side effects. Reports of using cinacalcet for NSHPT have increased, however, the effective dose for pediatric patients from the onset of symptoms through infancy has not been established. Case presentation We describe the clinical course of a newborn with a de novo R185Q mutation in the calcium-sensing receptor (CASR) gene, causing NSHPT. The infant received cinacalcet from the first days of life until 1 year of age. Conclusions Cinacalcet therapy effectively controlled the patient's serum calcium, phosphorus, and parathyroid hormone (PTH) levels without side effects.

Design and evaluation of a sequence capture system for genome-wide SNP genotyping in highly heterozygous plant genomes: a case study with a keystone Neotropical hardwood tree genome

Targeted sequence capture coupled to high-throughput sequencing has become a powerful method for the study of genome-wide sequence variation. Following our recent development of a genome assembly for the Pink Ipê tree (Handroanthus impetiginosus), a widely distributed Neotropical timber species, we now report the development of a set of 24,751 capture probes for single-nucleotide polymorphisms (SNPs) characterization and genotyping across 18,216 distinct loci, sampling more than 10 Mbp of the species genome. This system identifies nearly 200,000 SNPs located inside or in close proximity to almost 14,000 annotated protein-coding genes, generating quality genotypic data in populations spanning wide geographic distances across the species native range. To provide recommendations for future developments of similar systems for highly heterozygous plant genomes we investigated issues such as probe design, sequencing coverage and bioinformatics, including the evaluation of the capture efficiency and a reassessment of the technical reproducibility of the assay for SNPs recall and genotyping precision. Our results highlight the value of a detailed probe screening on a preliminary genome assembly to produce reliable data for downstream genetic studies. This work should inspire and assist the development of similar genomic resources for other orphan crops and forest trees with highly heterozygous genomes.

MATAM: reconstruction of phylogenetic marker genes from short sequencing reads in metagenomes

Motivation

Advances in the sequencing of uncultured environmental samples, dubbed metagenomics, raise a growing need for accurate taxonomic assignment. Accurate identification of organisms present within a community is essential to understanding even the most elementary ecosystems. However, current high-throughput sequencing technologies generate short reads which partially cover full-length marker genes and this poses difficult bioinformatic challenges for taxonomy identification at high resolution.

Results

We designed MATAM, a software dedicated to the fast and accurate targeted assembly of short reads sequenced from a genomic marker of interest. The method implements a stepwise process based on construction and analysis of a read overlap graph. It is applied to the assembly of 16S rRNA markers and is validated on simulated, synthetic and genuine metagenomes. We show that MATAM outperforms other available methods in terms of low error rates and recovered fractions and is suitable to provide improved assemblies for precise taxonomic assignments.

Availability and implementation

https://github.com/bonsai-team/matam.

Contact

pierre.pericard@gmail.com or helene.touzet@univ-lille1.fr.

Supplementary information

Supplementary data are available at Bioinformatics online.

Hybrid Capture-Based Next Generation Sequencing and Its Application to Human Infectious Diseases

This review describes target-enrichment approaches followed by next generation sequencing and their recent application to the research and diagnostic field of modern and past infectious human diseases caused by viruses, bacteria, parasites and fungi.

Chloromethane formation and degradation in the fern phyllosphere

Chloromethane (CH₃Cl) is the most abundant halogenated trace gas in the atmosphere. It plays an important role in natural stratospheric ozone destruction. Current estimates of the global CH₃Cl budget are approximate. The strength of the CH₃Cl global sink by microbial degradation in soils and plants is under discussion. Some plants, particularly ferns, have been identified as substantial emitters of CH₃Cl. Their ability to degrade CH₃Cl remains uncertain. In this study, we investigated the potential of leaves from 3 abundant ferns (Osmunda regalis, Cyathea cooperi, Dryopteris filix-mas) to produce and degrade CH₃Cl by measuring their production and consumption rates and their stable carbon and hydrogen isotope signatures. Investigated ferns are able to degrade CH₃Cl at rates from 2.1 to 17 and 0.3 to 0.9μgg_dw^-1day^-1 for C. cooperi and D. filix-mas respectively, depending on CH₃Cl supplementation and temperature. The stable carbon isotope enrichment factor of remaining CH₃Cl was -39±13‰, whereas negligible isotope fractionation was observed for hydrogen (-8±19‰). In contrast, O. regalis did not consume CH₃Cl, but produced it at rates ranging from 0.6 to 128μgg_dw^-1day^-1, with stable isotope values of -97±8‰ for carbon and -202±10‰ for hydrogen, respectively. Even though the 3 ferns showed clearly different formation and consumption patterns, their leaf-associated bacterial diversity was not notably different. Moreover, we did not detect genes associated with the only known chloromethane utilization pathway "cmu" in the microbial phyllosphere of the investigated ferns. Our study suggests that still unknown CH₃Cl biodegradation processes on plants play an important role in global cycling of atmospheric CH₃Cl.

Comparing 16S rDNA amplicon sequencing and hybridization capture for pea aphid microbiota diversity analysis

OBJECTIVE

Targeted sequencing of 16S rDNA amplicons is routinely used for microbial community profiling but this method suffers several limitations such as bias affinity of universal primers and short read size. Gene capture by hybridization represents a promising alternative. Here we used a metagenomic extract from the pea aphid Acyrthosiphon pisum to compare the performances of two widely used PCR primer pairs with DNA capture, based on solution hybrid selection.

RESULTS

All methods produced an exhaustive description of the 8 bacterial taxa known to be present in this sample. In addition, the methods yielded similar quantitative results, with the number of reads strongly correlating with quantitative PCR controls. Both methods can thus be considered as qualitatively and quantitatively robust on such a sample with low microbial complexity.

Hybridization capture reveals microbial diversity missed using current profiling methods

BACKGROUND

Microorganisms comprise the majority of living organisms on our planet. For many years, exploration of the composition of microbial communities has been performed through the PCR-based study of the small subunit rRNA gene due to its high conservation across the domains of life. The application of this method has resulted in the discovery of many unexpected evolutionary lineages. However, amplicon sequencing is subject to numerous biases, with some taxa being missed, and is limited by the read length of second-generation sequencing platforms, which drastically reduces the phylogenetic resolution.

RESULTS

Here, we describe a hybridization capture strategy that allows the enrichment of 16S rRNA genes from metagenomic samples and enables an exhaustive identification and a complete reconstruction of the biomarker. Applying this approach to a microbial mock community and a soil sample, we demonstrated that hybridization capture is able to reveal greater microbial diversity than 16S rDNA amplicon sequencing and shotgun sequencing. The reconstruction of full-length 16S rRNA genes facilitated the improvement of phylogenetic resolution and the discovery of novel prokaryotic taxa.

CONCLUSIONS

Our results demonstrate that hybridization capture can lead to major breakthroughs in our understanding of microbial diversity, overcoming the limitations of conventional 16S rRNA gene studies. If applied to a broad range of environmental samples, this innovative approach could reveal the undescribed diversity of the still underexplored microbial communities and could provide a better understanding of ecosystem function.

Efficient genome-wide genotyping strategies and data integration in crop plants

Next-generation sequencing (NGS) has revolutionized plant and animal research by providing powerful genotyping methods. This review describes and discusses the advantages, challenges and, most importantly, solutions to facilitate data processing, the handling of missing data, and cross-platform data integration. Next-generation sequencing technologies provide powerful and flexible genotyping methods to plant breeders and researchers. These methods offer a wide range of applications from genome-wide analysis to routine screening with a high level of accuracy and reproducibility. Furthermore, they provide a straightforward workflow to identify, validate, and screen genetic variants in a short time with a low cost. NGS-based genotyping methods include whole-genome re-sequencing, SNP arrays, and reduced representation sequencing, which are widely applied in crops. The main challenges facing breeders and geneticists today is how to choose an appropriate genotyping method and how to integrate genotyping data sets obtained from various sources. Here, we review and discuss the advantages and challenges of several NGS methods for genome-wide genetic marker development and genotyping in crop plants. We also discuss how imputation methods can be used to both fill in missing data in genotypic data sets and to integrate data sets obtained using different genotyping tools. It is our hope that this synthetic view of genotyping methods will help geneticists and breeders to integrate these NGS-based methods in crop plant breeding and research.

Capture Hybridization of Long-Range DNA Fragments for High-Throughput Sequencing

Capture hybridization coupled with high-throughput sequencing (HTS) has become one of the most popular approaches to address some scientific problems not only for fundamental evolution but also for ecology and human disease in recent years. However, the technical problem of limited probe capture ability affects its widespread application. Here, we propose to capture hybridize long-range DNA fragments for HTS (termed LR-LCH). We provide a case of three amphibian samples to examine LR-LCH with 2 kb libraries and comparison of standard capture hybridization with 480 bp libraries. Capture sensitivity increased from an average 13.57% of standard capture hybridization to an average 19.80% of LR-LCH; capture efficiency also increased from an average 72.56% of standard capture hybridization to an average 97.71% of LR-LCH. These indicate that longer fragments in the library generally contain both relatively variable regions and relatively conservative regions. The divergent parts of target DNA are enriched along with conservative parts of DNA sequence that effectively captured during hybridization. We present a protocol that allows users to overcome the low capture sensitivity problem for high divergent regions.

Comparative scaffolding and gap filling of ancient bacterial genomes applied to two ancient Yersinia pestis genomes

Yersinia pestis is the causative agent of the bubonic plague, a disease responsible for several dramatic historical pandemics. Progress in ancient DNA (aDNA) sequencing rendered possible the sequencing of whole genomes of important human pathogens, including the ancient Y. pestis strains responsible for outbreaks of the bubonic plague in London in the 14th century and in Marseille in the 18th century, among others. However, aDNA sequencing data are still characterized by short reads and non-uniform coverage, so assembling ancient pathogen genomes remains challenging and often prevents a detailed study of genome rearrangements. It has recently been shown that comparative scaffolding approaches can improve the assembly of ancient Y. pestis genomes at a chromosome level. In the present work, we address the last step of genome assembly, the gap-filling stage. We describe an optimization-based method AGapEs (ancestral gap estimation) to fill in inter-contig gaps using a combination of a template obtained from related extant genomes and aDNA reads. We show how this approach can be used to refine comparative scaffolding by selecting contig adjacencies supported by a mix of unassembled aDNA reads and comparative signal. We applied our method to two Y. pestis data sets from the London and Marseilles outbreaks, for which we obtained highly improved genome assemblies for both genomes, comprised of, respectively, five and six scaffolds with 95 % of the assemblies supported by ancient reads. We analysed the genome evolution between both ancient genomes in terms of genome rearrangements, and observed a high level of synteny conservation between these strains.

Investigation of orthologous pathogen recognition gene-rich regions in solanaceous species

Pathogen receptor proteins such as receptor-like protein (RLP), receptor-like kinase (RLK), and nucleotide-binding leucine-rich repeat (NLR) play a leading role in plant immunity activation. The genome architecture of such genes has been extensively investigated in several plant species. However, we still know little about their elaborate reorganization that arose during the plant speciation process. Using recently released pepper and eggplant genome sequences, we were able to identify 1097 pathogen recognition genes (PRGs) in the cultivated pepper Zunla-1 and 775 in the eggplant line Nakate-Shinkuro. The retrieved genes were analysed for their tendency to cluster, using different methods to infer the means of grouping. Orthologous relationships among clustering loci were found, and interesting reshuffling within given loci was observed for each analysed species. The information obtained was integrated into a comparative map to highlight the evolutionary dynamics in which the PRG loci were involved. Diversification of 14 selected PRG-rich regions was also explored using a DNA target-enrichment approach. A large number of gene variants were found as well as rearrangements of sequences encoding single protein domain and changes in chromosome gene order among species. Gene duplication and transposition activity have clearly influenced plant genome R-gene architecture and diversification. Our findings contribute to addressing several biological questions concerning the parallel evolution that occurred between genomes of the family Solanaceae. Moreover, the integration of different methods proved a powerful approach to reconstruct the evolutionary history in plant families and to transfer important biology findings among plant genomes.

Revealing large metagenomic regions through long DNA fragment hybridization capture

BACKGROUND

High-throughput DNA sequencing technologies have revolutionized genomic analysis, including the de novo assembly of whole genomes from single organisms or metagenomic samples. However, due to the limited capacity of short-read sequence data to assemble complex or low coverage regions, genomes are typically fragmented, leading to draft genomes with numerous underexplored large genomic regions. Revealing these missing sequences is a major goal to resolve concerns in numerous biological studies.

METHODS

To overcome these limitations, we developed an innovative target enrichment method for the reconstruction of large unknown genomic regions. Based on a hybridization capture strategy, this approach enables the enrichment of large genomic regions allowing the reconstruction of tens of kilobase pairs flanking a short, targeted DNA sequence.

RESULTS

Applied to a metagenomic soil sample targeting the linA gene, the biomarker of hexachlorocyclohexane (HCH) degradation, our method permitted the enrichment of the gene and its flanking regions leading to the reconstruction of several contigs and complete plasmids exceeding tens of kilobase pairs surrounding linA. Thus, through gene association and genome reconstruction, we identified microbial species involved in HCH degradation which constitute targets to improve biostimulation treatments.

CONCLUSIONS

This new hybridization capture strategy makes surveying and deconvoluting complex genomic regions possible through large genomic regions enrichment and allows the efficient exploration of metagenomic diversity. Indeed, this approach enables to assign identity and function to microorganisms in natural environments, one of the ultimate goals of microbial ecology.

Discovering Protein-Coding Genes from the Environment: Time for the Eukaryotes?

Eukaryotic microorganisms from diverse environments encompass a large number of taxa, many of them still unknown to science. One strategy to mine these organisms for genes of biotechnological relevance is to use a pool of eukaryotic mRNA directly extracted from environmental samples. Recent reports demonstrate that the resulting metatranscriptomic cDNA libraries can be screened by expression in yeast for a wide range of genes and functions from many of the different eukaryotic taxa. In combination with novel emerging high-throughput technologies, we anticipate that this approach should contribute to exploring the functional diversity of the eukaryotic microbiota.

The sequence capture by hybridization: a new approach for revealing the potential of mono-aromatic hydrocarbons bioattenuation in a deep oligotrophic aquifer

The formation water of a deep aquifer (853 m of depth) used for geological storage of natural gas was sampled to assess the mono-aromatic hydrocarbons attenuation potential of the indigenous microbiota. The study of bacterial diversity suggests that Firmicutes and, in particular, sulphate-reducing bacteria (Peptococcaceae) predominate in this microbial community. The capacity of the microbial community to biodegrade toluene and m- and p-xylenes was demonstrated using a culture-based approach after several hundred days of incubation. In order to reveal the potential for biodegradation of these compounds within a shorter time frame, an innovative approach named the solution hybrid selection method, which combines sequence capture by hybridization and next-generation sequencing, was applied to the same original water sample. The bssA and bssA-like genes were investigated as they are considered good biomarkers for the potential of toluene and xylene biodegradation. Unlike a PCR approach which failed to detect these genes directly from formation water, this innovative strategy demonstrated the presence of the bssA and bssA-like genes in this oligotrophic ecosystem, probably harboured by Peptococcaceae. The sequence capture by hybridization shows significant potential to reveal the presence of genes of functional interest which have low-level representation in the biosphere.

OCaPPI-Db: an oligonucleotide probe database for pathogen identification through hybridization capture

The detection and identification of bacterial pathogens involved in acts of bio- and agroterrorism are essential to avoid pathogen dispersal in the environment and propagation within the population. Conventional molecular methods, such as PCR amplification, DNA microarrays or shotgun sequencing, are subject to various limitations when assessing environmental samples, which can lead to inaccurate findings. We developed a hybridization capture strategy that uses a set of oligonucleotide probes to target and enrich biomarkers of interest in environmental samples. Here, we present Oligonucleotide Capture Probes for Pathogen Identification Database (OCaPPI-Db), an online capture probe database containing a set of 1,685 oligonucleotide probes allowing for the detection and identification of 30 biothreat agents up to the species level. This probe set can be used in its entirety as a comprehensive diagnostic tool or can be restricted to a set of probes targeting a specific pathogen or virulence factor according to the user’s needs.

Database URL: http://ocappidb.uca.works

From Genomes to GENE-omes: Exome Sequencing Concept and Applications in Crop Improvement

Exome sequencing represents targeted capture and sequencing of 1-2% of 'high-value genomic regions' (subset of the genome) which are enriched for functional variants and harbors low level of repetitive regions. We discuss here an overview of exome sequencing, ways to approach plant exomes, and advantages and applicability of this powerful approach in deciphering functional regions of genomes. Though initially this approach was developed as an alternative to whole genome sequencing (WGS), but the multitude of benefits conferred by sequence capture via hybridization approaches created a niche for itself to solve many of biological riddles, particularly for resolving phylogenetic distances. The technique has also proved to be successful in understanding the basis of natural and induced molecular variation, marker development and developing genomic resources for complex, wild and non-model species, which are still intractable for WGS efforts. Thus, with profound applications of this powerful sequencing strategy, near future is expected to witness a collective expansion of both techniques, i.e., sequence capture via hybridization for evolutionary and ecological research and WGS approaches for its universal accessibility.

Taming the Past: Ancient DNA and the Study of Animal Domestication

During the last decade, ancient DNA research has been revolutionized by the availability of increasingly powerful DNA sequencing and ancillary genomics technologies, giving rise to the new field of paleogenomics. In this review, we show how our understanding of the genetic basis of animal domestication and the origins and dispersal of livestock and companion animals during the Upper Paleolithic and Neolithic periods is being rapidly transformed through new scientific knowledge generated with paleogenomic methods. These techniques have been particularly informative in revealing high-resolution patterns of artificial and natural selection and evidence for significant admixture between early domestic animal populations and their wild congeners.