Short-read Technologies Research Articles

Evaluation of commercial RNA extraction kits for long-read metatranscriptomics in soil

Metatranscriptomic analysis of the soil microbiome has the potential to reveal molecular mechanisms that drive soil processes regulated by the microbial community. Therefore, RNA samples must be of sufficient yield and quality to robustly quantify differential gene expression. While short-read sequencing technology is often favoured for metatranscriptomics, long-read sequencing has the potential to provide several benefits over short-read technologies. The ability to resolve complete transcripts on a portable sequencing platform for a relatively low capital expenditure makes Oxford Nanopore Technology an attractive prospect for addressing many of the challenges of soil metatranscriptomics. To fully enable long-read metatranscriptomic analysis of the functional molecular pathways expressed in these diverse habitats, RNA purification methods from soil must be optimised for long-read sequencing. Here we compare RNA samples purified using five commercially available extraction kits designed for use with soil. We found that the Qiagen RNeasy PowerSoil Total RNA Kit performed the best across RNA yield, quality and purity and was robust across different soil types. We found that sufficient sequencing depth can be achieved to characterise the active community for total RNA samples using Oxford Nanopore Technology, and discuss its current limitations for differential gene expression analysis in soil studies.

Decomposing a San Francisco estuary microbiome using long-read metagenomics reveals species- and strain-level dominance from picoeukaryotes to viruses.

Ocean and estuarine microbiomes play critical roles in global element cycling and ecosystem function. Despite the importance of these microbial communities, many species still have not been cultured in the lab. Environmental sequencing is the primary way the function and population dynamics of these communities can be studied. Long-read sequencing provides an avenue to overcome limitations of short-read technologies to obtain complete microbial genomes but comes with its own technical challenges, such as needed sequencing depth and obtaining high-quality DNA. We present here new sampling and bioinformatics methods to attempt decomposing an estuarine microbiome into its constituent genomes. Our results suggest there are only a few strains that comprise most of the species abundances from viruses to picoeukaryotes, and to fully decompose a metagenome of this diversity requires 1 Tbp of long-read sequencing. We anticipate that as long-read sequencing technologies continue to improve, less sequencing will be needed.

Assembly and analysis of sequence from a spring and winter type Camelina sativa by whole genome PacBio HiFi technologies

Camelina sativa is an emerging oilseed crop with potential for multi-cropping with more traditional cash crops. Until very recently, the only reference genome for camelina was produced using short read technologies and was considerably fragmented. To facilitate the mapping of traits associated with vernalization requirement and freezing tolerance in an F7 recombinant inbred line (RIL) population developed by crossing a spring- and winter-biotype of camelina, we performed long read sequencing using PacBio HiFi technology. Here we report on the genome assembly and gene annotation obtained from the long read sequencing of the two parental lines. Both assemblies formed 20 chromosomal units with a genome size of 667,897,934 and 714,129,675 bases for C046 and Joelle, respectively. Assessment of completeness by BUSCO analysis indicated the C046 and Joelle genomes were both 99.5 complete with 98.1 % of the conserved genes being duplicated. Skim-sequencing of the F6:7 RILs resulted in 823,142 markers that mapped to variants identified by comparisons of the parental genome sequences. These sequences provide a valuable resource for breeders seeking to improve the food and industrial attributes of camelina.

Maptcha: an efficient parallel workflow for hybrid genome scaffolding

BackgroundGenome assembly, which involves reconstructing a target genome, relies on scaffolding methods to organize and link partially assembled fragments. The rapid evolution of long read sequencing technologies toward more accurate long reads, coupled with the continued use of short read technologies, has created a unique need for hybrid assembly workflows. The construction of accurate genomic scaffolds in hybrid workflows is complicated due to scale, sequencing technology diversity (e.g., short vs. long reads, contigs or partial assemblies), and repetitive regions within a target genome.ResultsIn this paper, we present a new parallel workflow for hybrid genome scaffolding that would allow combining pre-constructed partial assemblies with newly sequenced long reads toward an improved assembly. More specifically, the workflow, called Maptcha, is aimed at generating long scaffolds of a target genome, from two sets of input sequences—an already constructed partial assembly of contigs, and a set of newly sequenced long reads. Our scaffolding approach internally uses an alignment-free mapping step to build a ⟨\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\langle $$\\end{document}contig,contig⟩\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\rangle $$\\end{document} graph using long reads as linking information. Subsequently, this graph is used to generate scaffolds. We present and evaluate a graph-theoretic “wiring” heuristic to perform this scaffolding step. To enable efficient workload management in a parallel setting, we use a batching technique that partitions the scaffolding tasks so that the more expensive alignment-based assembly step at the end can be efficiently parallelized. This step also allows the use of any standalone assembler for generating the final scaffolds.ConclusionsOur experiments with Maptcha on a variety of input genomes, and comparison against two state-of-the-art hybrid scaffolders demonstrate that Maptcha is able to generate longer and more accurate scaffolds substantially faster. In almost all cases, the scaffolds produced by Maptcha are at least an order of magnitude longer (in some cases two orders) than the scaffolds produced by state-of-the-art tools. Maptcha runs significantly faster too, reducing time-to-solution from hours to minutes for most input cases. We also performed a coverage experiment by varying the sequencing coverage depth for long reads, which demonstrated the potential of Maptcha to generate significantly longer scaffolds in low coverage settings (1×\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$1\ imes $$\\end{document}–10×\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$10\ imes $$\\end{document}).

Chromosome-Scale Assembly and Annotation of Eight Arabidopsis thaliana Ecotypes.

The plant Arabidopsis thaliana is a model system used by researchers through much of plant research. Recent efforts have focused on discovering the genomic variation found in naturally occurring ecotypes isolated from around the world. These ecotypes have come from diverse climates and therefore have faced and adapted to a variety of abiotic and biotic stressors. The sequencing and comparative analysis of these genomes can offer insight into the adaptive strategies of plants. While there are a large number of ecotype genome sequences available, the majority were created using short-read technology. Mapping of short-reads containing structural variation to a reference genome bereft of that variation leads to incorrect mapping of those reads, resulting in a loss of genetic information and introduction of false heterozygosity. For this reason, long-read de novo sequencing of genomes is required to resolve structural variation events. In this article, we sequenced the genomes of eight natural variants of A. thaliana using nanopore sequencing. This resulted in highly contiguous assemblies with >95% of the genome contained within five contigs. The sequencing results from this study include five ecotypes from relict and African populations, an area of untapped genetic diversity. With this study, we increase the knowledge of diversity we have across A. thaliana ecotypes and contribute to ongoing production of an A. thaliana pan-genome.

Whole-genome sequence of Mucor velutinosus NIH1002, a strain associated with disseminated disease.

The draft genome of Mucor velutinosus NIH1002, a 2011 isolate from a case of disseminated disease, was sequenced using PacBio long-read and HiSeq short-read technologies. The genome has 43 contigs, an N50 of 2.65 Mb, and 13,295 protein-coding genes. It is the most complete M. velutinosus genome to date.

CoRAL accurately resolves extrachromosomal DNA genome structures with long-read sequencing.

Extrachromosomal DNA (ecDNA) is a central mechanism for focal oncogene amplification in cancer, occurring in ∼15% of early-stage cancers and ∼30% of late-stage cancers. ecDNAs drive tumor formation, evolution, and drug resistance by dynamically modulating oncogene copy number and rewiring gene-regulatory networks. Elucidating the genomic architecture of ecDNA amplifications is critical for understanding tumor pathology and developing more effective therapies. Paired-end short-read (Illumina) sequencing and mapping have been utilized to represent ecDNA amplifications using a breakpoint graph, in which the inferred architecture of ecDNA is encoded as a cycle in the graph. Traversals of breakpoint graphs have been used to successfully predict ecDNA presence in cancer samples. However, short-read technologies are intrinsically limited in the identification of breakpoints, phasing together complex rearrangements and internal duplications, and deconvolution of cell-to-cell heterogeneity of ecDNA structures. Long-read technologies, such as from Oxford Nanopore Technologies, have the potential to improve inference as the longer reads are better at mapping structural variants and are more likely to span rearranged or duplicated regions. Here, we propose Complete Reconstruction of Amplifications with Long reads (CoRAL) for reconstructing ecDNA architectures using long-read data. CoRAL reconstructs likely cyclic architectures using quadratic programming that simultaneously optimizes parsimony of reconstruction, explained copy number, and consistency of long-read mapping. CoRAL substantially improves reconstructions in extensive simulations and 10 data sets from previously characterized cell lines compared with previous short- and long-read-based tools. As long-read usage becomes widespread, we anticipate that CoRAL will be a valuable tool for profiling the landscape and evolution of focal amplifications in tumors.

Effects of different assembly strategies on gene annotation in activated sludge

Activated sludge comprises diverse bacteria, fungi, and other microorganisms, featuring a rich repertoire of genes involved in antibiotic resistance, pollutant degradation, and elemental cycling. In this regard, hybrid assembly technology can revolutionize metagenomics by detecting greater gene diversity in environmental samples. Nonetheless, the optimal utilization and comparability of genomic information between hybrid assembly and short- or long-read technology remain unclear. To address this gap, we compared the performance of the hybrid assembly, short- and long-read technologies, abundance and diversity of annotated genes, and taxonomic diversity by analysing 46, 161, and 45 activated sludge metagenomic datasets, respectively. The results revealed that hybrid assembly technology exhibited the best performance, generating the most contiguous and longest contigs but with a lower proportion of high-quality metagenome-assembled genomes than short-read technology. Compared with short- or long-read technologies, hybrid assembly technology can detect a greater diversity of microbiota and antibiotic resistance genes, as well as a wider range of potential hosts. However, this approach may yield lower gene abundance and pathogen detection. Our study revealed the specific advantages and disadvantages of hybrid assembly and short- and long-read applications in wastewater treatment plants, and our approach could serve as a blueprint to be extended to terrestrial environments.

Mibianto: ultra-efficient online microbiome analysis through k-mer based metagenomics.

Quantifying microbiome species and composition from metagenomic assays is often challenging due to its time-consuming nature and computational complexity. In Bioinformatics, k-mer-based approaches were long established to expedite the analysis of large sequencing data and are now widely used to annotate metagenomic data. We make use of k-mer counting techniques for efficient and accurate compositional analysis of microbiota from whole metagenome sequencing. Mibianto solves this problem by operating directly on read files, without manual preprocessing or complete data exchange. It handles diverse sequencing platforms, including short single-end, paired-end, and long read technologies. Our sketch-based workflow significantly reduces the data volume transferred from the user to the server (up to 99.59% size reduction) to subsequently perform taxonomic profiling with enhanced efficiency and privacy. Mibianto offers functionality beyond k-mer quantification; it supports advanced community composition estimation, including diversity, ordination, and differential abundance analysis. Our tool aids in the standardization of computational workflows, thus supporting reproducibility of scientific sequencing studies. It is adaptable to small- and large-scale experimental designs and offers a user-friendly interface, thus making it an invaluable tool for both clinical and research-oriented metagenomic studies. Mibianto is freely available without the need for a login at: https://www.ccb.uni-saarland.de/mibianto.

Nanopore and Illumina sequencing reveal different viral populations from human gut samples.

The advent of viral metagenomics, or viromics, has improved our knowledge and understanding of global viral diversity. High-throughput sequencing technologies enable explorations of the ecological roles, contributions to host metabolism, and the influence of viruses in various environments, including the human intestinal microbiome. However, bacterial metagenomic studies frequently have the advantage. The adoption of advanced technologies like long-read sequencing has the potential to be transformative in refining viromics and metagenomics. Here, we examined the effectiveness of long-read and hybrid sequencing by comparing Illumina short-read and Oxford Nanopore Technology (ONT) long-read sequencing technologies and different assembly strategies on recovering viral genomes from human faecal samples. Our findings showed that if a single sequencing technology is to be chosen for virome analysis, Illumina is preferable due to its superior ability to recover fully resolved viral genomes and minimise erroneous genomes. While ONT assemblies were effective in recovering viral diversity, the challenges related to input requirements and the necessity for amplification made it less ideal as a standalone solution. However, using a combined, hybrid approach enabled a more authentic representation of viral diversity to be obtained within samples.

Analysis and benchmarking of small and large genomic variants across tandem repeats.

Tandem repeats (TRs) are highly polymorphic in the human genome, have thousands of associated molecular traits and are linked to over 60 disease phenotypes. However, they are often excluded from at-scale studies because of challenges with variant calling and representation, as well as a lack of a genome-wide standard. Here, to promote the development of TR methods, we created a catalog of TR regions and explored TR properties across 86 haplotype-resolved long-read human assemblies. We curated variants from the Genome in a Bottle (GIAB) HG002 individual to create a TR dataset to benchmark existing and future TR analysis methods. We also present an improved variant comparison method that handles variants greater than 4 bp in length and varying allelic representation. The 8.1% of the genome covered by the TR catalog holds ~24.9% of variants per individual, including 124,728 small and 17,988 large variants for the GIAB HG002 'truth-set' TR benchmark. We demonstrate the utility of this pipeline across short-read and long-read technologies.

The fly route of extended-spectrum-β-lactamase-producing Enterobacteriaceae dissemination in a cattle farm: from the ecosystem to the molecular scale

IntroductionThis study aimed to understand the origin and to explain the maintenance of extended-spectrum β-lactamase (ESBL) Enterobacteriaceae isolated from food-producing animals in a third-generation cephalosporin (3GC)-free farm.MethodsCulture and molecular approaches were used to test molecules other than 3GC such as antibiotics (tetracycline and oxytetracycline), antiparasitics (ivermectin, flumethrin, fenbendazol, and amitraz), heavy metal [arsenic, HNO3, aluminum, HNO3, cadmium (CdSO4), zinc (ZnCl2), copper (CuSO4), iron (FeCl3), and aluminum (Al2SO4)], and antioxidant (butylated hydroxytoluene) as sources of selective pressure. Whole-genome sequencing using short read (Illumina™) and long read (Nanopore™) technologies was performed on 34 genomes. In silico gene screening and comparative analyses were used to characterize the genetic determinants of resistance, their mobility, and the genomic relatedness among isolates.ResultsOur analysis unveiled a low diversity among the animal ESBL-producing strains. Notably, E. coli ST3268 was recurrently isolated from both flies (n = 9) and cattle (n = 5). These E. coli ST3268/blaCTX-M-15/blaTEM-1B have accumulated multiple plasmids and genes, thereby representing a reservoir of resistance and virulence factors. Our findings suggest that flies could act as effective mechanical vectors for antimicrobial gene transfer and are capable of transporting resistant bacteria across different environments and to multiple hosts, facilitating the spread of pathogenic traits. A significantly higher mean minimum inhibitory concentration of oxytetracycline (841.4 ± 323.5 mg/L vs. 36.0 ± 52.6 mg/L, p = 0.0022) in ESBL E. coli than in non-ESBL E. coli and blaCTX-M-15 gene overexpression in oxytetracycline-treated vs. untreated ESBL E. coli (RQOxy = 3.593, p = 0.024) confirmed oxytetracycline as a source of selective pressure in ESBL E. coli.DiscussionThe occurrence of ESBL E. coli in a farm without 3GC use is probably due to an as yet undefined human origin of Enterobacteriaceae blaCTX-M-15 gene transmission to animals in close contact with cattle farm workers and the maintenance of the local ESBL E. coli reservoir by a high fly diversity and oxytetracycline selective pressure. These findings highlight the critical need for stringent vector control to mitigate antimicrobial resistance spread for preserving public health. Addressing this issue necessitates a multifaceted approach combining microbial genetics, vector ecology, and farm management practices.

Chromosome-level reference genome and resequencing of 322 accessions reveal evolution, genomic imprint and key agronomic traits in adzuki bean.

Adzuki bean (Vigna angularis) is an important legume crop cultivated in over 30 countries worldwide. We developed a high-quality chromosome-level reference genome of adzuki bean cultivar Jingnong6 by combining PacBio Sequel long-read sequencing with short-read and Hi-C technologies. The assembled genome covers 97.8% of the adzuki bean genome with a contig N50 of approximately 16 Mb and a total of 32 738 protein-coding genes. We also generated a comprehensive genome variation map of adzuki bean by whole-genome resequencing (WGRS) of 322 diverse adzuki beans accessions including both wild and cultivated. Furthermore, we have conducted comparative genomics and a genome-wide association study (GWAS) on key agricultural traits to investigate the evolution and domestication. GWAS identified several candidate genes, including VaCycA3;1, VaHB15, VaANR1 and VaBm, that exhibited significant associations with domestication traits. Furthermore, we conducted functional analyses on the roles of VaANR1 and VaBm in regulating seed coat colour. We provided evidence for the highest genetic diversity of wild adzuki (Vigna angularis var. nipponensis) in China with the presence of the most original wild adzuki bean, and the occurrence of domestication process facilitating transition from wild to cultigen. The present study elucidates the genetic basis of adzuki bean domestication traits and provides crucial genomic resources to support future breeding efforts in adzuki bean.

The long and short of it: benchmarking viromics using Illumina, Nanopore and PacBio sequencing technologies.

Viral metagenomics has fuelled a rapid change in our understanding of global viral diversity and ecology. Long-read sequencing and hybrid assembly approaches that combine long- and short-read technologies are now being widely implemented in bacterial genomics and metagenomics. However, the use of long-read sequencing to investigate viral communities is still in its infancy. While Nanopore and PacBio technologies have been applied to viral metagenomics, it is not known to what extent different technologies will impact the reconstruction of the viral community. Thus, we constructed a mock bacteriophage community of previously sequenced phage genomes and sequenced them using Illumina, Nanopore and PacBio sequencing technologies and tested a number of different assembly approaches. When using a single sequencing technology, Illumina assemblies were the best at recovering phage genomes. Nanopore- and PacBio-only assemblies performed poorly in comparison to Illumina in both genome recovery and error rates, which both varied with the assembler used. The best Nanopore assembly had errors that manifested as SNPs and INDELs at frequencies 41 and 157 % higher than found in Illumina only assemblies, respectively. While the best PacBio assemblies had SNPs at frequencies 12 and 78 % higher than found in Illumina-only assemblies, respectively. Despite high-read coverage, long-read-only assemblies recovered a maximum of one complete genome from any assembly, unless reads were down-sampled prior to assembly. Overall the best approach was assembly by a combination of Illumina and Nanopore reads, which reduced error rates to levels comparable with short-read-only assemblies. When using a single technology, Illumina only was the best approach. The differences in genome recovery and error rates between technology and assembler had downstream impacts on gene prediction, viral prediction, and subsequent estimates of diversity within a sample. These findings will provide a starting point for others in the choice of reads and assembly algorithms for the analysis of viromes.

Long-Read Sequencing and De Novo Genome Assembly Pipeline of Two Plasmodium falciparum Clones (Pf3D7, PfW2) Using Only the PromethION Sequencer from Oxford Nanopore Technologies without Whole-Genome Amplification.

Antimalarial drug resistance has become a real public health problem despite WHO measures. New sequencing technologies make it possible to investigate genomic variations associated with resistant phenotypes at the genome-wide scale. Based on the use of hemisynthetic nanopores, the PromethION technology from Oxford Nanopore Technologies can produce long-read sequences, in contrast to previous short-read technologies used as the gold standard to sequence Plasmodium. Two clones of P. falciparum (Pf3D7 and PfW2) were sequenced in long-read using the PromethION sequencer from Oxford Nanopore Technologies without genomic amplification. This made it possible to create a processing analysis pipeline for human Plasmodium with ONT Fastq only. De novo assembly revealed N50 lengths of 18,488 kb and 17,502 kb for the Pf3D7 and PfW2, respectively. The genome size was estimated at 23,235,407 base pairs for the Pf3D7 clone and 21,712,038 base pairs for the PfW2 clone. The average genome coverage depth was estimated at 787X and 653X for the Pf3D7 and PfW2 clones, respectively. This study proposes an assembly processing pipeline for the human Plasmodium genome using software adapted to large ONT data and the high AT percentage of Plasmodium. This search provides all the parameters which were optimized for use with the software selected in the pipeline.

Full-length transcriptome and RNA-Seq analyses reveal the resistance mechanism of sesame in response to Corynespora cassiicola

BackgroundCorynespora leaf spot is a common leaf disease occurring in sesame, and the disease causes leaf yellowing and even shedding, which affects the growth quality of sesame. At present, the mechanism of sesame resistance to this disease is still unclear. Understanding the resistance mechanism of sesame to Corynespora leaf spot is highly important for the control of infection. In this study, the leaves of the sesame resistant variety (R) and the sesame susceptible variety (S) were collected at 0–48 hpi for transcriptome sequencing, and used a combined third-generation long-read and next-generation short-read technology approach to identify some key genes and main pathways related to resistance.ResultsThe gene expression levels of the two sesame varieties were significantly different at 0, 6, 12, 24, 36 and 48 hpi, indicating that the up-regulation of differentially expressed genes in the R might enhanced the resistance. Moreover, combined with the phenotypic observations of sesame leaves inoculated at different time points, we found that 12 hpi was the key time point leading to the resistance difference between the two sesame varieties at the molecular level. The WGCNA identified two modules significantly associated with disease resistance, and screened out 10 key genes that were highly expressed in R but low expressed in S, which belonged to transcription factors (WRKY, AP2/ERF-ERF, and NAC types) and protein kinases (RLK-Pelle_DLSV, RLK-Pelle_SD-2b, and RLK-Pelle_WAK types). These genes could be the key response factors in the response of sesame to infection by Corynespora cassiicola. GO and KEGG enrichment analysis showed that specific modules could be enriched, which manifested as enrichment in biologically important pathways, such as plant signalling hormone transduction, plant-pathogen interaction, carbon metabolism, phenylpropanoid biosynthesis, glutathione metabolism, MAPK and other stress-related pathways.ConclusionsThis study provides an important resource of genes contributing to disease resistance and will deepen our understanding of the regulation of disease resistance, paving the way for further molecular breeding of sesame.

2476. Universal detection of healthcare-related bacterial transmissions using de novo assembly with short-read whole genome sequencing

Abstract Background Whole genome sequencing (WGS) is the gold standard for molecular evidence of pathogen transmission. To date, most WGS investigations within healthcare settings have used low cost short-read technologies (e.g., Illumina), aligning reads to a high quality reference genome. However, this relies on either public repositories having closely-related high quality references, which is not guaranteed for arbitrary clinical isolates, or additional long-read sequencing. De novo assembly avoids reliance on references, but the utility of fragmented assemblies of short reads is not well-established in hospital outbreak investigations. Methods 24 clusters of clinical isolates were identified by automated and routine methods at 17 U.S. hospitals and underwent paired-end WGS on Illumina platforms, followed by two parallel analyses. For a reference-free analysis, de novo assembly was performed with Shovill, assemblies were clustered with Mash, core-genome alignments were created with Harvest, and finally single nucleotide polymorphism (SNP) distances were visualized using the PathoSPOT phylogenomics toolkit to highlight likely transmission events. For a reference-guided analysis, references were selected for each isolate using StrainGST, and SNPs called by read alignment were used to define phylogenomic distances. Finally, conclusions from both analyses were compared. Results 92 bacterial isolates underwent WGS, of which 86 (93%) completed de novo assembly and passed quality control, comprising 10 isolates of Acinetobacter baumanii, 4 of Klebsiella oxytoca, 46 of Klebsiella pneumoniae, and 26 of Staphylococcus aureus. Reference-free analysis with PathoSPOT identified close relatedness consistent with transmission in 13 of 24 clusters (54%), while the reference-guided analysis identified the same in 12 clusters (50%). Overall, 21 clusters (88%) showed concordance between reference-free and reference-guided methods. Cohen’s kappa indicated substantial agreement (0.75; 95% CI, 0.49–1.0). Conclusion De novo assembly of short reads is tenable for WGS investigation of hospital outbreaks, elides the need for reference genomes or long-read sequencing, and compares favorably with SNP calling via read alignment for detecting likely healthcare-related transmissions. Disclosures Laura E. McLean, M.Ed., HCA Healthcare: Stocks/Bonds Susan S. Huang, MD MPH, Medline Industries, Inc: Conducted studies whereby participating nursing homes and hospital patients received cleaning & antiseptic products|Xttrium Laboratories: Conducted studies where participating nursing homes and hospital patients received antiseptic products Yonatan H. Grad, MD, PhD, Day Zero Diagnostics: Board Member|GSK: Advisor/Consultant

Panmixia in the American eel extends to its tropical range of distribution: Biological implications and policymaking challenges.

The American eel (Anguilla rostrata) has long been regarded as a panmictic fish and has been confirmed as such in the northern part of its range. In this paper, we tested for the first time whether panmixia extends to the tropical range of the species. To do so, we first assembled a reference genome (975 Mbp, 19 chromosomes) combining long (PacBio and Nanopore and short (Illumina paired-end) reads technologies to support both this study and future research. To test for population structure, we estimated genotype likelihoods from low-coverage whole-genome sequencing of 460 American eels, collected at 21 sampling sites (in seven geographic regions) ranging from Canada to Trinidad and Tobago. We estimated genetic distance between regions, performed ADMIXTURE-like clustering analysis and multivariate analysis, and found no evidence of population structure, thus confirming that panmixia extends to the tropical range of the species. In addition, two genomic regions with putative inversions were observed, both geographically widespread and present at similar frequencies in all regions. We discuss the implications of lack of genetic population structure for the species. Our results are key for the future genomic research in the American eel and the implementation of conservation measures throughout its geographic range. Additionally, our results can be applied to fisheries management and aquaculture of the species.

Emerging Opportunities to Study Mobile Element Insertions and Their Source Elements in an Expanding Universe of Sequenced Human Genomes.

Three mobile element classes, namely Alu, LINE-1 (L1), and SVA elements, remain actively mobile in human genomes and continue to produce new mobile element insertions (MEIs). Historically, MEIs have been discovered and studied using several methods, including: (1) Southern blots, (2) PCR (including PCR display), and (3) the detection of MEI copies from young subfamilies. We are now entering a new phase of MEI discovery where these methods are being replaced by whole genome sequencing and bioinformatics analysis to discover novel MEIs. We expect that the universe of sequenced human genomes will continue to expand rapidly over the next several years, both with short-read and long-read technologies. These resources will provide unprecedented opportunities to discover MEIs and study their impact on human traits and diseases. They also will allow the MEI community to discover and study the source elements that produce these new MEIs, which will facilitate our ability to study source element regulation in various tissue contexts and disease states. This, in turn, will allow us to better understand MEI mutagenesis in humans and the impact of this mutagenesis on human biology.

Benchmarking long-read RNA-sequencing analysis tools using in silico mixtures.

The lack of benchmark data sets with inbuilt ground-truth makes it challenging to compare the performance of existing long-read isoform detection and differential expression analysis workflows. Here, we present a benchmark experiment using two human lung adenocarcinoma cell lines that were each profiled in triplicate together with synthetic, spliced, spike-in RNAs (sequins). Samples were deeply sequenced on both Illumina short-read and Oxford Nanopore Technologies long-read platforms. Alongside the ground-truth available via the sequins, we created in silico mixture samples to allow performance assessment in the absence of true positives or true negatives. Our results show that StringTie2 and bambu outperformed other tools from the six isoform detection tools tested, DESeq2, edgeR and limma-voom were best among the five differential transcript expression tools tested and there was no clear front-runner for performing differential transcript usage analysis between the five tools compared, which suggests further methods development is needed for this application.