Linked-read Sequencing Research Articles

Comparative genome analysis and the genome-shaping role of long terminal repeat retrotransposons in the evolutionary divergence of fungal pathogens Blastomyces dermatitidis and Blastomyces gilchristii.

Blastomyces dermatitidis and Blastomyces gilchristii are cryptic species of fungi that cause blastomycosis, an often severe disease involving pulmonary infection capable of systemic dissemination. While these species appear morphologically identical, differences exist in the genetic makeup, geographical range, and possibly the clinical presentation of infection. Here, we show genetic divergence between the cryptic species through both a Blastomyces species tree constructed from orthologous protein sequences and whole genome single-nucleotide variant phylogenomic analysis. Following linked-read sequencing and de novo genome assembly, we characterized and compared the genomes of 3 B. dermatitidis and 3 B. gilchristii isolates. The B. gilchristii genomes (73.25-75.4 Mb) were ∼8 Mb larger than the B. dermatitidis genomes (64.88-66.61 Mb). Average nucleotide identity was lower between genomes of different species than genomes of the same species, yet functional classification of genes suggested similar proteomes. The most striking difference involved long terminal repeat retrotransposons. Although the same retrotransposon elements were detected in the genomes, the quantity of elements differed between the 2 species. Gypsy retrotransposon content was significantly higher in B. gilchristii (38.04-39.26 Mb) than in B. dermatitidis (30.85-32.40 Mb), accounting for the majority of genome size difference between species. Age estimation and phylogenetic analysis of the reverse transcriptase domains suggested that these retrotransposons are relatively ancient, with genome insertion predating the speciation of B. dermatitidis and B. gilchristii. We postulate that different trajectories of genome contraction led to genetic incompatibility, reproductive isolation, and speciation, highlighting the role of transposable elements in fungal evolution.

Targeted Linked-Read Sequencing for Direct Haplotype Phasing of Parental GJB2/SLC26A4 Alleles: A Universal and Dependable Noninvasive Prenatal Diagnosis Method Applied to Autosomal Recessive Nonsyndromic Hearing Loss in At-Risk Families

Noninvasive prenatal diagnosis (NIPD) for autosomal recessive nonsyndromic hearing loss (ARNSHL) has been rarely reported until recent years. Additionally, the existing method can not be used for challenging genome loci (eg, copy number variations, deletions, inversions, or gene recombinants) or on families without proband genotype. This study assessed the performance of relative haplotype dosage analysis (RHDO)-based NIPD for identifying fetal genotyping in pregnancies at risk of ARNSHL. Fifty couples carrying pathogenic variants associated with ARNSHL in either GJB2 or SLC26A4 were recruited. The RHDO-based targeted linked-read sequencing combined with whole gene coverage probes was used to genotype the fetal cell-free DNA of 49 families who met the quality control standard. Fetal amniocyte samples were genotyped using invasive prenatal diagnosis (IPD) to assess the performance of NIPD. The NIPD results showed 100% (49/49) concordance with those obtained through IPD. Two families with copy number variation and recombination were also successfully identified. Sufficient specific informative single-nucleotide polymorphisms for haplotyping, as well as the fetal cell-free DNA concentration and sequencing depth, are prerequisites for RHDO-based NIPD. This method has the merits of covering the entire genes of GJB2 and SLC26A4, qualifying for copy number variation and recombination analysis with remarkable sensitivity and specificity. Therefore, it has clinical potential as an alternative to traditional IPD for ARNSHL.

Transposable Element Insertions Are Associated with Batesian Mimicry in the Pantropical Butterfly Hypolimnas misippus.

Hypolimnas misippus is a Batesian mimic of the toxic African Queen butterfly (Danaus chrysippus). Female H. misippus butterflies use two major wing patterning loci (M and A) to imitate three color morphs of D. chrysippus found in different regions of Africa. In this study, we examine the evolution of the M locus and identify it as an example of adaptive atavism. This phenomenon involves a morphological reversion to an ancestral character that results in an adaptive phenotype. We show that H. misippus has re-evolved an ancestral wing pattern present in other Hypolimnas species, repurposing it for Batesian mimicry of a D. chrysippus morph. Using haplotagging, a linked-read sequencing technology, and our new analytical tool, Wrath, we discover two large transposable element insertions located at the M locus and establish that these insertions are present in the dominant allele responsible for producing mimetic phenotype. By conducting a comparative analysis involving additional Hypolimnas species, we demonstrate that the dominant allele is derived. This suggests that, in the derived allele, the transposable elements disrupt a cis-regulatory element, leading to the reversion to an ancestral phenotype that is then utilized for Batesian mimicry of a distinct model, a different morph of D. chrysippus. Our findings present a compelling instance of convergent evolution and adaptive atavism, in which the same pattern element has independently evolved multiple times in Hypolimnas butterflies, repeatedly playing a role in Batesian mimicry of diverse model species.

Haplotype information of large neuromuscular disease genes provided by linked-read sequencing has a potential to increase diagnostic yield.

Rare or novel missense variants in large genes such as TTN and NEB are frequent in the general population, which hampers the interpretation of putative disease-causing biallelic variants in patients with sporadic neuromuscular disorders. Often, when the first initial genetic analysis is performed, the reconstructed haplotype, i.e. phasing information of the variants is missing. Segregation analysis increases the diagnostic turnaround time and is not always possible if samples from family members are lacking. To overcome this difficulty, we investigated how well the linked-read technology succeeded to phase variants in these large genes, and whether it improved the identification of structural variants. Linked-read sequencing data of nemaline myopathy, distal myopathy, and proximal myopathy patients were analyzed for phasing, single nucleotide variants, and structural variants. Variant phasing was successful in the large muscle genes studied. The longest continuous phase blocks were gained using high-quality DNA samples with long DNA fragments. Homozygosity increased the number of phase blocks, especially in exome sequencing samples lacking intronic variation. In our cohort, linked-read sequencing added more information about the structural variation but did not lead to a molecular genetic diagnosis. The linked-read technology can support the clinical diagnosis of neuromuscular and other genetic disorders.

Divergence and gene flow history at two large chromosomal inversions underlying ecotype differentiation in the long-snouted seahorse.

Chromosomal inversions can play an important role in divergence and reproductive isolation by building and maintaining distinct allelic combinations between evolutionary lineages. Alternatively, they can take the form of balanced polymorphisms that segregate within populations until one arrangement becomes fixed. Many questions remain about how inversion polymorphisms arise, how they are maintained over the long term, and ultimately, whether and how they contribute to speciation. The long-snouted seahorse (Hippocampus guttulatus) is genetically subdivided into geographic lineages and marine-lagoon ecotypes, with shared structural variation underlying lineage and ecotype divergence. Here, we aim to characterize structural variants and to reconstruct their history and suspected role in ecotype formation. We generated a near chromosome-level genome assembly and described genome-wide patterns of diversity and divergence through the analysis of 112 whole-genome sequences from Atlantic, Mediterranean, and Black Sea populations. By also analysing linked-read sequencing data, we found evidence for two chromosomal inversions that were several megabases in length and showed contrasting allele frequency patterns between lineages and ecotypes across the species range. We reveal that these inversions represent ancient intraspecific polymorphisms, one likely being maintained by divergent selection and the other by pseudo-overdominance. A possible selective coupling between the two inversions was further supported by the absence of specific haplotype combinations and a putative functional interaction between the two inversions in reproduction. Lastly, we detected gene flux eroding divergence between inverted alleles at varying levels for the two inversions, with a likely impact on their dynamics and contribution to divergence and speciation.

LRTK: a platform agnostic toolkit for linked-read analysis of both human genome and metagenome.

Linked-read sequencing technologies generate high-base quality short reads that contain extrapolative information on long-range DNA connectedness. These advantages of linked-read technologies are well known and have been demonstrated in many human genomic and metagenomic studies. However, existing linked-read analysis pipelines (e.g., Long Ranger) were primarily developed to process sequencing data from the human genome and are not suited for analyzing metagenomic sequencing data. Moreover, linked-read analysis pipelines are typically limited to 1 specific sequencing platform. To address these limitations, we present the Linked-Read ToolKit (LRTK), a unified and versatile toolkit for platform agnostic processing of linked-read sequencing data from both human genome and metagenome. LRTK provides functions to perform linked-read simulation, barcode sequencing error correction, barcode-aware read alignment and metagenome assembly, reconstruction of long DNA fragments, taxonomic classification and quantification, and barcode-assisted genomic variant calling and phasing. LRTK has the ability to process multiple samples automatically and provides users with the option to generate reproducible reports during processing of raw sequencing data and at multiple checkpoints throughout downstream analysis. We applied LRTK on linked reads from simulation, mock community, and real datasets for both human genome and metagenome. We showcased LRTK's ability to generate comparative performance results from preceding benchmark studies and to report these results in publication-ready HTML document plots. LRTK provides comprehensive and flexible modules along with an easy-to-use Python-based workflow for processing linked-read sequencing datasets, thereby filling the current gap in the field caused by platform-centric genome-specific linked-read data analysis tools.

Chromosome-scale genome of the human blood fluke Schistosoma mekongi and its implications for public health.

Schistosoma mekongi is a human blood fluke causing schistosomiasis that threatens approximately 1.5 million humans in the world. Nonetheless, the limited available S. mekongi genomic resources have hindered understanding of its biology and parasite-host interactions for disease management and pathogen control. The aim of our study was to integrate multiple technologies to construct a high-quality chromosome-level assembly of the S. mekongi genome. The reference genome for S. mekongi was generated through integrating Illumina, PacBio sequencing, 10 × Genomics linked-read sequencing, and high-throughput chromosome conformation capture (Hi-C) methods. In this study, we conducted de novo assembly, alignment, and gene prediction to assemble and annotate the genome. Comparative genomics allowed us to compare genomes across different species, shedding light on conserved regions and evolutionary relationships. Additionally, our transcriptomic analysis focused on genes associated with parasite-snail interactions in S. mekongi infection. We employed gene ontology (GO) enrichment analysis for functional annotation of these genes. In the present study, the S. mekongi genome was both assembled into 8 pseudochromosomes with a length of 404Mb, with contig N50 and scaffold N50 lengths of 1168kb and 46,759kb, respectively. We detected that 43% of the genome consists of repeat sequences and predicted 9103 protein-coding genes. We also focused on proteases, particularly leishmanolysin-like metalloproteases (M8), which are crucial in the invasion of hosts by 12 flatworm species. Through phylogenetic analysis, it was discovered that the M8 gene exhibits lineage-specific amplification among the genus Schistosoma. Lineage-specific expansion of M8 was observed in blood flukes. Additionally, the results of the RNA-seq revealed that a mass of genes related to metabolic and biosynthetic processes were up-regulated, which might be beneficial for cercaria production. This study delivers a high-quality, chromosome-scale reference genome of S. mekongi, enhancing our understanding of the divergence and evolution of Schistosoma. The molecular research conducted here also plays a pivotal role in drug discovery and vaccine development. Furthermore, our work greatly advances the understanding of host-parasite interactions, providing crucial insights for schistosomiasis intervention strategies.

BLR: a flexible pipeline for haplotype analysis of multiple linked-read technologies.

Linked-read sequencing promises a one-method approach for genome-wide insights including single nucleotide variants (SNVs), structural variants, and haplotyping. We introduce Barcode Linked Reads (BLR), an open-source haplotyping pipeline capable of handling millions of barcodes and data from multiple linked-read technologies including DBS, 10× Genomics, TELL-seq and stLFR. Running BLR on DBS linked-reads yielded megabase-scale phasing with low (<0.2%) switch error rates. Of 13616 protein-coding genes phased in the GIAB benchmark set (v4.2.1), 98.6% matched the BLR phasing. In addition, large structural variants showed concordance with HPRC-HG002 reference assembly calls. Compared to diploid assembly with PacBio HiFi reads, BLR phasing was more continuous when considering switch errors. We further show that integrating long reads at low coverage (∼10×) can improve phasing contiguity and reduce switch errors in tandem repeats. When compared to Long Ranger on 10× Genomics data, BLR showed an increase in phase block N50 with low switch-error rates. For TELL-Seq and stLFR linked reads, BLR generated longer or similar phase block lengths and low switch error rates compared to results presented in the original publications. In conclusion, BLR provides a flexible workflow for comprehensive haplotype analysis of linked reads from multiple platforms.

Whole genomic analysis reveals atypical non-homologous off-target large structural variants induced by CRISPR-Cas9-mediated genome editing

CRISPR-Cas9 genome editing has promising therapeutic potential for genetic diseases and cancers, but safety could be a concern. Here we use whole genomic analysis by 10x linked-read sequencing and optical genome mapping to interrogate the genome integrity after editing and in comparison to four parental cell lines. In addition to the previously reported large structural variants at on-target sites, we identify heretofore unexpected large chromosomal deletions (91.2 and 136 Kb) at atypical non-homologous off-target sites without sequence similarity to the sgRNA in two edited lines. The observed large structural variants induced by CRISPR-Cas9 editing in dividing cells may result in pathogenic consequences and thus limit the usefulness of the CRISPR-Cas9 editing system for disease modeling and gene therapy. In this work, our whole genomic analysis may provide a valuable strategy to ensure genome integrity after genomic editing to minimize the risk of unintended effects in research and clinical applications.

Linked-read based analysis of the medulloblastoma genome.

Medulloblastoma is the most common type of malignant pediatric brain tumor with group 4 medulloblastomas (G4 MBs) accounting for 40% of cases. However, the molecular mechanisms that underlie this subgroup are still poorly understood. Point mutations are detected in a large number of genes at low incidence per gene while the detection of complex structural variants in recurrently affected genes typically requires the application of long-read technologies. Here, we applied linked-read sequencing, which combines the long-range genome information of long-read sequencing with the high base pair accuracy of short read sequencing and very low sample input requirements. We demonstrate the detection of complex structural variants and point mutations in these tumors, and, for the first time, the detection of extrachromosomal DNA (ecDNA) with linked-reads. We provide further evidence for the high heterogeneity of somatic mutations in G4 MBs and add new complex events associated with it. We detected several enhancer-hijacking events, an ecDNA containing the MYCN gene, and rare structural rearrangements, such a chromothripsis in a G4 medulloblastoma, chromoplexy involving 8 different chromosomes, a TERT gene rearrangement, and a PRDM6 duplication.

MTG-Link: leveraging barcode information from linked-reads to assemble specific loci

BackgroundLocal assembly with short and long reads has proven to be very useful in many applications: reconstruction of the sequence of a locus of interest, gap-filling in draft assemblies, as well as alternative allele reconstruction of large Structural Variants. Whereas linked-read technologies have a great potential to assemble specific loci as they provide long-range information while maintaining the power and accuracy of short-read sequencing, there is a lack of local assembly tools for linked-read data.ResultsWe present MTG-Link, a novel local assembly tool dedicated to linked-reads. The originality of the method lies in its read subsampling step which takes advantage of the barcode information contained in linked-reads mapped in flanking regions. We validated our approach on several datasets from different linked-read technologies. We show that MTG-Link is able to assemble successfully large sequences, up to dozens of Kb. We also demonstrate that the read subsampling step of MTG-Link considerably improves the local assembly of specific loci compared to other existing short-read local assembly tools. Furthermore, MTG-Link was able to fully characterize large insertion variants and deletion breakpoints in a human genome and to reconstruct dark regions in clinically-relevant human genes. It also improved the contiguity of a 1.3 Mb locus of biological interest in several individual genomes of the mimetic butterfly Heliconius numata.ConclusionsMTG-Link is an efficient local assembly tool designed for different linked-read sequencing technologies. MTG-Link source code is available at https://github.com/anne-gcd/MTG-Link and as a Bioconda package.

Chromosome-length genome assembly and karyotype of the endangered black-footed ferret (Mustela nigripes).

The black-footed ferret (Mustela nigripes) narrowly avoided extinction to become an oft-cited example of the benefits of intensive management, research, and collaboration to save a species through ex situ conservation breeding and reintroduction into its former range. However, the species remains at risk due to possible inbreeding, disease susceptibility, and multiple fertility challenges. Here, we report the de novo genome assembly of a male black-footed ferret generated through a combination of linked-read sequencing, optical mapping, and Hi-C proximity ligation. In addition, we report the karyotype for this species, which was used to anchor and assign chromosome numbers to the chromosome-length scaffolds. The draft assembly was ~2.5 Gb in length, with 95.6% of it anchored to 19 chromosome-length scaffolds, corresponding to the 2n = 38 chromosomes revealed by the karyotype. The assembly has contig and scaffold N50 values of 148.8 kbp and 145.4 Mbp, respectively, and is up to 96% complete based on BUSCO analyses. Annotation of the assembly, including evidence from RNA-seq data, identified 21,406 protein-coding genes and a repeat content of 37.35%. Phylogenomic analyses indicated that the black-footed ferret diverged from the European polecat/domestic ferret lineage 1.6 million yr ago. This assembly will enable research on the conservation genomics of black-footed ferrets and thereby aid in the further restoration of this endangered species.

Abstract LB048: AR-targeted linked-read sequencing reveals complex AR structural variants in CRPC tumors

Abstract Inhibition of the androgen receptor (AR) transcription factor is the mainstay therapeutic strategy for systemic prostate cancer. However, development of a resistance phenotype termed castration resistant prostate cancer (CRPC) enables disease progression that underlies virtually all prostate cancer mortality. In the majority of CRPC, AR is transcriptionally re-activated in part due to accumulation of numerous alterations in the AR gene. These genomic alterations include missense mutations in exons that encode the ligand binding domain of AR, amplification of the AR gene body and/or an upstream transcriptional enhancer, and AR structural variants (SVs) resulting from inaccurate repair of DNA double strand breaks. In CRPC patient samples, these AR genomic alterations co-occur within tumors, which makes elucidating the fitness contribution of each alteration difficult. While short-read DNA sequencing has allowed for characterization of these events, they are underpowered to preserve linkage information. To address this challenge, we performed AR-targeted linked-read DNA-sequencing of CRPC cell line and patient derived xenograft (PDX) models with the goal of resolving sub-clonal and complex AR SVs. This approach confirmed the heterogeneity of AR SVs in these PDXs, including phased AR SVs that can be traced back to a single high molecular weight DNA molecule. These phased AR SVs were consistently found to co-occur in CRPC specimens with amplification of AR. Our results define a new class of AR SVs where multiple rearrangements co-occur on amplified DNA molecules to yield a complex array of AR gene structures in CRPC. Citation Format: Andrej Zivanovic, Sarah A. Munro, Jeffrey T. Miller, Yingming Li, Courtney N. Passow, Scott M. Dehm. AR-targeted linked-read sequencing reveals complex AR structural variants in CRPC tumors [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 2 (Clinical Trials and Late-Breaking Research); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(8_Suppl):Abstract nr LB048.

Benchmarking genome assembly methods on metagenomic sequencing data.

Metagenome assembly is an efficient approach to reconstruct microbial genomes from metagenomic sequencing data. Although short-read sequencing has been widely used for metagenome assembly, linked- and long-read sequencing have shown their advancements in assembly by providing long-range DNA connectedness. Many metagenome assembly tools were developed to simplify the assembly graphs and resolve the repeats in microbial genomes. However, there remains no comprehensive evaluation of metagenomic sequencing technologies, and there is a lack of practical guidance on selecting the appropriate metagenome assembly tools. This paper presents a comprehensive benchmark of 19 commonly used assembly tools applied to metagenomic sequencing datasets obtained from simulation, mock communities or human gut microbiomes. These datasets were generated using mainstream sequencing platforms, such as Illumina and BGISEQ short-read sequencing, 10x Genomics linked-read sequencing, and PacBio and Oxford Nanopore long-read sequencing. The assembly tools were extensively evaluated against many criteria, which revealed that long-read assemblers generated high contig contiguity but failed to reveal some medium- and high-quality metagenome-assembled genomes (MAGs). Linked-read assemblers obtained the highest number of overall near-complete MAGs from the human gut microbiomes. Hybrid assemblers using both short- and long-read sequencing were promising methods to improve both total assembly length and the number of near-complete MAGs. This paper also discussed the running time and peak memory consumption of these assembly tools and provided practical guidance on selecting them.

Benchmarking of structural variant detection in the tetraploid potato genome using linked-read sequencing

It has recently been shown that structural variants (SV) can have a higher impact on gene expression variation compared to single nucleotide variants (SNV) in different plant species. Additionally, SV were associated with phenotypic variation in several crops. However, compared to the established SV detection based on short-read sequencing, less approaches were described for linked-read based SV calling. We therefore evaluated the performance of six linked-read SV callers compared to an established short-read SV caller based on simulated linked-reads in tetraploid potato. The objectives of our study were to i) compare the performance of SV callers based on linked-read sequencing to short-read sequencing, ii) examine the influence of SV type, SV length, haplotype incidence (HI), as well as sequencing coverage on the SV calling performance in the tetraploid potato genome, and iii) evaluate the accuracy of detecting insertions by linked-read compared to short-read sequencing. We observed high break point resolutions (BPR) detecting short SV and slightly lower BPR for large SV. Our observations highlighted the importance of short-read signals provided by Manta and LinkedSV to detect short SV. Manta and NAIBR performed well for detecting larger deletions, inversions, and duplications. Detected large SV were weakly influenced by the HI. Furthermore, we illustrated that large insertions can be assembled by Novel-X. Our results suggest the usage of the short-read and linked-read SV callers Manta, NAIBR, LinkedSV, and Novel-X based on at least 90x linked-read sequencing coverage to ensure the detection of a broad range of SV in the tetraploid potato genome.

High-quality read-based phasing of cystic fibrosis cohort informs genetic understanding of disease modification.

Phasing of heterozygous alleles is critical for interpretation of cis-effects of disease-relevant variation. We sequenced 477 individuals with cystic fibrosis (CF) using linked-read sequencing, which display an average phase block N50 of 4.39 Mb. We use these samples to construct a graph representation of CFTR haplotypes, demonstrating its utility for understanding complex CF alleles. These are visualized in a Web app, CFTbaRcodes, that enables interactive exploration of CFTR haplotypes present in this cohort. We perform fine-mapping and phasing of the chr7q35 trypsinogen locus associated with CF meconium ileus, an intestinal obstruction at birth associated with more severe CF outcomes and pancreatic disease. A 20-kb deletion polymorphism and a PRSS2 missense variant p.Thr8Ile (rs62473563) are shown to independently contribute to meconium ileus risk (p= 0.0028, p= 0.011, respectively) and are PRSS2 pancreas eQTLs (p=9.5×10-7 and p= 1.4×10-4, respectively), suggesting the mechanism by which these polymorphisms contribute to CF. The phase information from linked reads provides a putative causal explanation for variation at a CF-relevant locus, which also has implications for the genetic basis of non-CF pancreatitis, to which this locus has been reported to contribute.

On the origin and structure of haplotype blocks.

The term "haplotype block" is commonly used in the developing field of haplotype-based inference methods. We argue that the term should be defined based on the structure of the Ancestral Recombination Graph (ARG), which contains complete information on the ancestry of a sample. We use simulated examples to demonstrate key features of the relationship between haplotype blocks and ancestral structure, emphasizing the stochasticity of the processes that generate them. Even the simplest cases of neutrality or of a "hard" selective sweep produce a rich structure, often missed by commonly used statistics. We highlight a number of novel methods for inferring haplotype structure, based on the full ARG, or on a sequence of trees, and illustrate how they can be used to define haplotype blocks using an empirical data set. While the advent of new, computationally efficient methods makes it possible to apply these concepts broadly, they (and additional new methods) could benefit from adding features to explore haplotype blocks, as we define them. Understanding and applying the concept of the haplotype block will be essential to fully exploit long and linked-read sequencing technologies.

RF27 | LBSUN342 Ar-targeted Linked-read Sequencing Reveals Complex Ar Structural Variants In Castration Resistant Prostate Cancer Models

Abstract Androgen deprivation therapies are employed for the treatment of advanced prostate cancer to suppress the transcriptional activity of the androgen receptor (AR). However, over the course of therapy, AR transcriptional activity re-emerges, and the disease progresses to castration resistant prostate cancer (CRPC). During CRPC progression, AR becomes heavily altered by genomic aberrations that include amplification, missense mutations, and structural variants (SVs). AR SVs are heterogeneous across CRPC specimens, with variability noted in SV class, breakpoint locations, and co-occurrence with amplification. Further variability has been observed for sub-clonal enrichment of AR SVs within tumors. A barrier to understanding the clinical significance of AR SVs is the challenge of deciphering this high degree of variability from short-read DNA-sequencing data. To address this challenge, we performed AR-targeted linked-read DNA-sequencing of CRPC cell line and patient-derived xenograft (PDX) models with the goal of resolving sub-clonal and complex AR SVs. We benchmarked this approach using prostate cancer cell lines and LuCaP86.2 PDX tissue harboring known AR SVs. Using this benchmarking approach, we identified linked AR SVs that co-occurred on single DNA molecules in LuCaP35CR and LuCaP105CR PDX models. We also identified a novel AR SV in the LuCaP77 PDX model and a castration-resistant derivative, LuCaP77CR. These results have revealed heterogeneous trajectories of AR SVs in CRPC and define a new class of complex AR SVs where multiple rearrangements co-occur on single DNA molecules to yield highly rearranged AR gene structures. Presentation: Monday, June 13, 2022 1:05 p.m. - 1:10 p.m.

QuickDeconvolution: fast and scalable deconvolution of linked-read sequencing data

QuickDeconvolution: fast and scalable deconvolution of linked-read sequencing data

Genomic novelty within a "great speciator" revealed by a high-quality reference genome of the collared kingfisher (Todiramphus chloris collaris).

Islands are natural laboratories for studying patterns and processes of evolution. Research on island endemic birds has revealed elevated speciation rates and rapid phenotypic evolution in several groups (e.g. white-eyes, Darwin's finches). However, understanding the evolutionary processes behind these patterns requires an understanding of how genotypes map to novel phenotypes. To date, there are few high-quality reference genomes for species found on islands. Here, we sequence the genome of one of Ernst Mayr's "great speciators," the collared kingfisher (Todiramphus chloris collaris). Utilizing high molecular weight DNA and linked-read sequencing technology, we assembled a draft high-quality genome with highly contiguous scaffolds (scaffold N50 = 19 Mb). Based on universal single-copy orthologs, we estimated a gene space completeness of 96.6% for the draft genome assembly. The population demographic history analyses reveal a distinct pattern of contraction and expansion in population size throughout the Pleistocene. Comparative genomic analysis of gene family evolution revealed that species-specific and rapidly expanding gene families in the collared kingfisher (relative to other Coraciiformes) are mainly involved in the ErbB signaling pathway and focal adhesion. Todiramphus kingfishers are a species-rich group that has become a focus of speciation research. This draft genome will be a platform for future taxonomic, phylogeographic, and speciation research in the group. For example, target genes will enable testing of changes in sensory structures associated with changes in vision and taste genes across kingfishers.