Fragmented Genome Assemblies Research Articles

Patchwork: Alignment-Based Retrieval and Concatenation of Phylogenetic Markers from Genomic Data.

Low-coverage whole-genome sequencing (also known as "genome skimming") is becoming an increasingly affordable approach to large-scale phylogenetic analyses. While already routinely used to recover organellar genomes, genome skimming is rather rarely utilized for recovering single-copy nuclear markers. One reason might be that only few tools exist to work with this data type within a phylogenomic context, especially to deal with fragmented genome assemblies. We here present a new software tool called Patchwork for mining phylogenetic markers from highly fragmented short-read assemblies as well as directly from sequence reads. Patchwork is an alignment-based tool that utilizes the sequence aligner DIAMOND and is written in the programming language Julia. Homologous regions are obtained via a sequence similarity search, followed by a "hit stitching" phase, in which adjacent or overlapping regions are merged into a single unit. The novel sliding window algorithm trims away any noncoding regions from the resulting sequence. We demonstrate the utility of Patchwork by recovering near-universal single-copy orthologs within a benchmarking study, and we additionally assess the performance of Patchwork in comparison with other programs. We find that Patchwork allows for accurate retrieval of (putatively) single-copy genes from genome skimming data sets at different sequencing depths with high computational speed, outperforming existing software targeting similar tasks. Patchwork is released under the GNU General Public License version 3. Installation instructions, additional documentation, and the source code itself are all available via GitHub at https://github.com/fethalen/Patchwork.

Syntenic relationship of chromosomes in Strongyloides species and Rhabditophanes diutinus based on the chromosome-level genome assemblies.

The Strongyloides clade, to which the parasitic nematode genus Strongyloides belongs, contains taxa with diverse lifestyles, ranging from free-living to obligate vertebrate parasites. Reproductive strategies are also diverse in this group of nematodes, employing not only sexual reproduction but also parthenogenesis, making it an attractive group to study genome adaptation to specific conditions. An in-depth understanding of genome evolution, however, has been hampered by fragmented genome assemblies. In this study, we generated chromosome-level genome assemblies for two Strongyloides species and the outgroup species Rhabditophanes diutinus using long-read sequencing and high-throughput chromosome conformation capture (Hi-C). Our synteny analyses revealed a clearer picture of chromosome evolution in this group, suggesting that a functional sex chromosome has been maintained throughout the group. We further investigated sex chromosome dynamics in the lifecycle of Strongyloides ratti and found that bivalent formation in oocytes appears to be important for male production in the mitotic parthenogenesis. This article is part of the Theo Murphy meeting issue 'Strongyloides: omics to worm-free populations'.

Adversarial and variational autoencoders improve metagenomic binning

Assembly of reads from metagenomic samples is a hard problem, often resulting in highly fragmented genome assemblies. Metagenomic binning allows us to reconstruct genomes by re-grouping the sequences by their organism of origin, thus representing a crucial processing step when exploring the biological diversity of metagenomic samples. Here we present Adversarial Autoencoders for Metagenomics Binning (AAMB), an ensemble deep learning approach that integrates sequence co-abundances and tetranucleotide frequencies into a common denoised space that enables precise clustering of sequences into microbial genomes. When benchmarked, AAMB presented similar or better results compared with the state-of-the-art reference-free binner VAMB, reconstructing ~7% more near-complete (NC) genomes across simulated and real data. In addition, genomes reconstructed using AAMB had higher completeness and greater taxonomic diversity compared with VAMB. Finally, we implemented a pipeline Integrating VAMB and AAMB that enabled improved binning, recovering 20% and 29% more simulated and real NC genomes, respectively, compared to VAMB, with moderate additional runtime.

PlASgraph2: using graph neural networks to detect plasmid contigs from an assembly graph.

Identification of plasmids from sequencing data is an important and challenging problem related to antimicrobial resistance spread and other One-Health issues. We provide a new architecture for identifying plasmid contigs in fragmented genome assemblies built from short-read data. We employ graph neural networks (GNNs) and the assembly graph to propagate the information from nearby nodes, which leads to more accurate classification, especially for short contigs that are difficult to classify based on sequence features or database searches alone. We trained plASgraph2 on a data set of samples from the ESKAPEE group of pathogens. plASgraph2 either outperforms or performs on par with a wide range of state-of-the-art methods on testing sets of independent ESKAPEE samples and samples from related pathogens. On one hand, our study provides a new accurate and easy to use tool for contig classification in bacterial isolates; on the other hand, it serves as a proof-of-concept for the use of GNNs in genomics. Our software is available at https://github.com/cchauve/plasgraph2 and the training and testing data sets are available at https://github.com/fmfi-compbio/plasgraph2-datasets.

Improved Genome Assembly Resource of the Plant Pathogen Fusarium avenaceum

Fusarium avenaceum is a generalist plant pathogen of concern due to its potential to produce mycotoxins on plant products. Previous research efforts have sequenced and assembled genomes of F. avenaceum. However, those works relied on limited next-generation sequencing technologies that resulted in fragmented and incomplete genome assemblies. To address this, we utilized high-depth third-generation long-read sequencing and several different genome assembly software to generate a new, highly contiguous genome of F. avenaceum. Moreover, we conducted a thorough annotation of the genome using a mix of long-read cDNA and short-read RNAseq data. Our genome was more contiguous than the current reference genome Fustri1, matched the estimated genome size and chromosome number, and contained a similar number of annotated genes to other F. avenaceum genome assemblies. Lastly, we conducted a secondary metabolite (SM) cluster analysis that identified 60 gene clusters associated with SM biosynthesis, five more than the reference F. avenaceum genome. In conclusion, our genome and associated annotation information will help advance research on plant-fungal interactions, food safety, and Fusarium spp. diversity. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license .

Chromosome-level assembly of the Phytophthora agathidicida genome reveals adaptation in effector gene families.

Phytophthora species are notorious plant pathogens, with some causing devastating tree diseases that threaten the survival of their host species. One such example is Phytophthora agathidicida, the causal agent of kauri dieback - a root and trunk rot disease that kills the ancient, iconic and culturally significant tree species, Agathis australis (New Zealand kauri). A deeper understanding of how Phytophthora pathogens infect their hosts and cause disease is critical for the development of effective treatments. Such an understanding can be gained by interrogating pathogen genomes for effector genes, which are involved in virulence or pathogenicity. Although genome sequencing has become more affordable, the complete assembly of Phytophthora genomes has been problematic, particularly for those with a high abundance of repetitive sequences. Therefore, effector genes located in repetitive regions could be truncated or missed in a fragmented genome assembly. Using a combination of long-read PacBio sequences, chromatin conformation capture (Hi-C) and Illumina short reads, we assembled the P. agathidicida genome into ten complete chromosomes, with a genome size of 57 Mb including 34% repeats. This is the first Phytophthora genome assembled to chromosome level and it reveals a high level of syntenic conservation with the complete genome of Peronospora effusa, the only other completely assembled genome sequence of an oomycete. All P. agathidicida chromosomes have clearly defined centromeres and contain candidate effector genes such as RXLRs and CRNs, but in different proportions, reflecting the presence of gene family clusters. Candidate effector genes are predominantly found in gene-poor, repeat-rich regions of the genome, and in some cases showed a high degree of duplication. Analysis of candidate RXLR effector genes that occur in multicopy gene families indicated half of them were not expressed in planta. Candidate CRN effector gene families showed evidence of transposon-mediated recombination leading to new combinations of protein domains, both within and between chromosomes. Further analysis of this complete genome assembly will help inform new methods of disease control against P. agathidicida and other Phytophthora species, ultimately helping decipher how Phytophthora pathogens have evolved to shape their effector repertoires and how they might adapt in the future.

A chromosome-scale reference genome for Spironucleus salmonicida

Spironucleus salmonicida is a diplomonad causing systemic infection in salmon. The first S. salmonicida genome assembly was published 2014 and has been a valuable reference genome in protist research. However, the genome assembly is fragmented without assignment of the sequences to chromosomes. In our previous Giardia genome study, we have shown how a fragmented genome assembly can be improved with long-read sequencing technology complemented with optical maps. Combining Pacbio long-read sequencing technology and optical maps, we are presenting here this new S. salmonicida genome assembly in nine near-complete chromosomes with only three internal gaps at long repeats. This new genome assembly is not only more complete sequence-wise but also more complete at annotation level, providing more details into gene families, gene organizations and chromosomal structure. This near-complete reference genome will aid comparative genomics at chromosomal level, and serve as a valuable resource for the diplomonad community and protist research.

Novel Multiplex PCR Method and Genome Sequence-Based Analog for High-Resolution Subclonal Assignment and Characterization of Escherichia coli Sequence Type 131 Isolates.

ABSTRACTEscherichia coli sequence type 131 (ST131) is a pandemic, multidrug-resistant extraintestinal pathogen. The multiple distinctive ST131 subclones differ for rfb and fliC alleles (O and H antigens), fimH allele (type-1 fimbriae adhesin), resistance phenotype and genotype, clinical correlates, and host predilection. Current PCR assays for detecting ST131 and its main subclones offer limited sub-ST characterization. Here we combined 22 novel and 14 published primers for a multiplex PCR assay to detect and extensively characterize ST131 isolates. The primers target mdh36, gyrB47, trpA72, sbmA, plsB, nupC, rmuC, kefC, ybbW, the O16 and O25b rfb variants, five fimH alleles (fimH22, fimH27, fimH30, fimH35, and fimH41), two fliC alleles (H4 and H5), a (subclone-specific) fluoroquinolone resistance-associated parC allele, and a (subclone-specific) prophage marker. The resulting amplicons resolve 15 molecular subsets within ST131, including 3 within clade A (H41 subclone), 5 within clade B (H22 subclone), and 7 within clade C (H30 subclone), which includes subclones C0 (H30S: 2 subsets), C1 and C1-M27 (H30R1: 2 subsets), and C2 (H30Rx: 3 subsets). Validation in three laboratories showed that this assay provides a rapid, accurate, and portable method for rapidly detecting and characterizing E. coli ST131 and its key subsets. Additionally, for users with whole genome sequencing (WGS) capability, we developed a command-line executable called ST131Typer, an in silico version of the extended multiplex PCR assay. Its accuracy was 87.8%, with most issues due to incomplete or fragmented input genome assemblies. These two novel assays should facilitate detailed ST131 subtyping using either endpoint PCR or WGS.IMPORTANCE These novel assays provide greater subclonal resolution and characterization of E. coli ST131 isolates than do the available comparable PCR assays, plus offer a novel sequence-based alternative to PCR. They may prove useful for molecular epidemiological studies, surveillance, and, potentially, clinical management.

ExceS-A: an exon-centric split aligner.

Spliced alignments are a key step in the construction of high-quality homology-based annotations of protein sequences. The exon/intron structure, which is computed as part of spliced alignment procedures, often conveys important information for the distinguishing paralogous members of gene families. Here we present an exon-centric pipeline for spliced alignment that is intended in particular for applications that involve exon-by-exon comparisons of coding sequences. We show that the simple, blat-based approach has advantages over established tools in particular for genes with very large introns and applications to fragmented genome assemblies.

Genome-Guided Discovery of Natural Products through Multiplexed Low-Coverage Whole-Genome Sequencing of Soil Actinomycetes on Oxford Nanopore Flongle.

ABSTRACTGenome mining is an important tool for discovery of new natural products; however, the number of publicly available genomes for natural product-rich microbes such as actinomycetes, relative to human pathogens with smaller genomes, is small. To obtain contiguous DNA assemblies and identify large (ca. 10 to greater than 100 kb) biosynthetic gene clusters (BGCs) with high GC (>70%) and high-repeat content, it is necessary to use long-read sequencing methods when sequencing actinomycete genomes. One of the hurdles to long-read sequencing is the higher cost. In the current study, we assessed Flongle, a recently launched platform by Oxford Nanopore Technologies, as a low-cost DNA sequencing option to obtain contiguous DNA assemblies and analyze BGCs. To make the workflow more cost-effective, we multiplexed up to four samples in a single Flongle sequencing experiment while expecting low-sequencing coverage per sample. We hypothesized that contiguous DNA assemblies might enable analysis of BGCs even at low sequencing depth. To assess the value of these assemblies, we collected high-resolution mass spectrometry data and conducted a multi-omics analysis to connect BGCs to secondary metabolites. In total, we assembled genomes for 20 distinct strains across seven sequencing experiments. In each experiment, 50% of the bases were in reads longer than 10 kb, which facilitated the assembly of reads into contigs with an average N50 value of 3.5 Mb. The programs antiSMASH and PRISM predicted 629 and 295 BGCs, respectively. We connected BGCs to metabolites for N,N-dimethyl cyclic-di-tryptophan, two novel lasso peptides, and three known actinomycete-associated siderophores, namely, mirubactin, heterobactin, and salinichelin.IMPORTANCE Short-read sequencing of GC-rich genomes such as those from actinomycetes results in a fragmented genome assembly and truncated biosynthetic gene clusters (often 10 to >100 kb long), which hinders our ability to understand the biosynthetic potential of a given strain and predict the molecules that can be produced. The current study demonstrates that contiguous DNA assemblies, suitable for analysis of BGCs, can be obtained through low-coverage, multiplexed sequencing on Flongle, which provides a new low-cost workflow ($30 to 40 per strain) for sequencing actinomycete strain libraries.

A linkage disequilibrium-based approach to position unmapped SNPs in crop species

BackgroundHigh-density SNP arrays are now available for a wide range of crop species. Despite the development of many tools for generating genetic maps, the genome position of many SNPs from these arrays is unknown. Here we propose a linkage disequilibrium (LD)-based algorithm to allocate unassigned SNPs to chromosome regions from sparse genetic maps. This algorithm was tested on sugarcane, wheat, and barley data sets. We calculated the algorithm’s efficiency by masking SNPs with known locations, then assigning their position to the map with the algorithm, and finally comparing the assigned and true positions.ResultsIn the 20-fold cross-validation, the mean proportion of masked mapped SNPs that were placed by the algorithm to a chromosome was 89.53, 94.25, and 97.23% for sugarcane, wheat, and barley, respectively. Of the markers that were placed in the genome, 98.73, 96.45 and 98.53% of the SNPs were positioned on the correct chromosome. The mean correlations between known and new estimated SNP positions were 0.97, 0.98, and 0.97 for sugarcane, wheat, and barley. The LD-based algorithm was used to assign 5920 out of 21,251 unpositioned markers to the current Q208 sugarcane genetic map, representing the highest density genetic map for this species to date.ConclusionsOur LD-based approach can be used to accurately assign unpositioned SNPs to existing genetic maps, improving genome-wide association studies and genomic prediction in crop species with fragmented and incomplete genome assemblies. This approach will facilitate genomic-assisted breeding for many orphan crops that lack genetic and genomic resources.

Genome survey of sago palm (Metroxylon sagu Rottboll)

The sago palm (Metroxylon sagu Rottboll) is a halophytic, underrated starch-producing palm that is both disease resistant and capable of producing a high yield of storage starch. It is deemed the future tree of life as it has a high potential to be cultivated for its starch to solve global food security issues in the future. Due to low scientific interest of this palm, scarcity of its genomic information has halted our comprehension on the mechanism behind its starch biosynthesis as well as salt and disease resistance pathways. In this study, we report on the genome survey and microsatellite analysis of sago palm. Based on the k-mer statistical approach, the sago palm has an estimated genome size of 464 Mbp to 616 Mbp, consisting of 31.9 to 35.7% repetitive sequence. The heterozygosity of the individual sequenced, a true sago palm clonal breed, was 0.63 to 0.66%. A sum of 123,430,321,878 clean bases were yielded from the 124,258,713,300 raw bases. BUSCO analysis revealed that only around 22% of the 1375 single copy plant genes are complete in both sago palm genome scaffolds and contigs whereas the remaining are missing (~45%) and fragmented (~32%). The fragmented genome assembly generated from short-read Illumina data is likely due to the high repetitiveness of the genome. It is recommended that a long-read sequencing to be carried out to improve the genome information of sago palm. The Illumina data generated in this study will be instrumental for polishing the genome assembly generated from a long-read assembly that is known to be error-prone.

Hybrid <em>De Novo</em> Genome Assembly for the Generation of Complete Genomes of Urinary Bacteria using Short- and Long-read Sequencing Technologies

Complete genome sequences provide valuable data for the understanding of genetic diversity and unique colonization factors of urinary microbes. These data may include mobile genetic elements, such as plasmids and extrachromosomal phage, that contribute to the dissemination of antimicrobial resistance and further complicate treatment of urinary tract infection (UTI). In addition to providing fine resolution of genome structure, complete, closed genomes allow for the detailed comparative genomics and evolutionary analyses. The generation of complete genomes de novo has long been a challenging task due to limitations of available sequencing technology. Paired-end Next Generation Sequencing (NGS) produces high quality short reads often resulting in accurate but fragmented genome assemblies. On the contrary, Nanopore sequencing provides long reads of lower quality normally leading to error-prone complete assemblies. Such errors may hamper genome-wide association studies or provide misleading variant analysis results. Therefore, hybrid approaches combining both short and long reads have emerged as reliable methods to achieve highly accurate closed bacterial genomes. Reported herein is a comprehensive method for the culture of diverse urinary bacteria, species identification by 16S rRNA gene sequencing, extraction of genomic DNA (gDNA), and generation of short and long reads by NGS and Nanopore platforms, respectively. Additionally, this method describes a bioinformatic pipeline of quality control, assembly, and gene prediction algorithms for the generation of annotated complete genome sequences. Combination of bioinformatic tools enables the selection of high quality read data for hybrid genome assembly and downstream analysis. The streamlined approach for the hybrid de novo genome assembly described in this protocol may be adapted for the use in any culturable bacteria.

Hybrid De Novo Genome Assembly for the Generation of Complete Genomes of Urinary Bacteria using Short- and Long-read Sequencing Technologies.

Complete genome sequences provide valuable data for the understanding of genetic diversity and unique colonization factors of urinary microbes. These data may include mobile genetic elements, such as plasmids and extrachromosomal phage, that contribute to the dissemination of antimicrobial resistance and further complicate treatment of urinary tract infection (UTI). In addition to providing fine resolution of genome structure, complete, closed genomes allow for the detailed comparative genomics and evolutionary analyses. The generation of complete genomes de novo has long been a challenging task due to limitations of available sequencing technology. Paired-end Next Generation Sequencing (NGS) produces high quality short reads often resulting in accurate but fragmented genome assemblies. On the contrary, Nanopore sequencing provides long reads of lower quality normally leading to error-prone complete assemblies. Such errors may hamper genome-wide association studies or provide misleading variant analysis results. Therefore, hybrid approaches combining both short and long reads have emerged as reliable methods to achieve highly accurate closed bacterial genomes. Reported herein is a comprehensive method for the culture of diverse urinary bacteria, species identification by 16S rRNA gene sequencing, extraction of genomic DNA (gDNA), and generation of short and long reads by NGS and Nanopore platforms, respectively. Additionally, this method describes a bioinformatic pipeline of quality control, assembly, and gene prediction algorithms for the generation of annotated complete genome sequences. Combination of bioinformatic tools enables the selection of high quality read data for hybrid genome assembly and downstream analysis. The streamlined approach for the hybrid de novo genome assembly described in this protocol may be adapted for the use in any culturable bacteria.

CRISPRclassify: Repeat-Based Classification of CRISPR Loci.

Detection and classification of CRISPR-Cas systems in metagenomic data have become increasingly prevalent in recent years due to their potential for diverse applications in genome editing. Traditionally, CRISPR-Cas systems are classified through reference-based identification of proximate cas genes. Here, we present a machine learning approach for the detection and classification of CRISPR loci using repeat sequences in a cas-independent context, enabling identification of unclassified loci missed by traditional cas-based approaches. Using biological attributes of the CRISPR repeat, the core element in CRISPR arrays, and leveraging methods from natural language processing, we developed a machine learning model capable of accurate classification of CRISPR loci in an extensive set of metagenomes, resulting in an F1 measure of 0.82 across all predictions and an F1 measure of 0.97 when limiting to classifications with probabilities >0.85. Furthermore, assessing performance on novel repeats yielded an F1 measure of 0.96. Although the performance of cas-based identification will exceed that of a repeat-based approach in many cases, CRISPRclassify provides an efficient approach to classification of CRISPR loci for cases in which cas gene information is unavailable, such as metagenomes and fragmented genome assemblies.

Chromosomal assembly of the nuclear genome of the endosymbiont-bearing trypanosomatid Angomonas deanei.

Angomonas deanei is an endosymbiont-bearing trypanosomatid with several highly fragmented genome assemblies and unknown chromosome number. We present an assembly of the A. deanei nuclear genome based on Oxford Nanopore sequence that resolves into 29 complete or close-to-complete chromosomes. The assembly has several previously unknown special features; it has a supernumerary chromosome, a chromosome with a 340-kb inversion, and there is a translocation between two chromosomes. We also present an updated annotation of the chromosomal genome with 10,365 protein-coding genes, 59 transfer RNAs, 26 ribosomal RNAs, and 62 noncoding RNAs.

High-resolution three-dimensional chromatin profiling of the Chinese hamster ovary cell genome.

Chinese hamster ovary (CHO) cell lines are the pillars of a multibillion‐dollar biopharmaceutical industry producing recombinant therapeutic proteins. The effects of local chromatin organization and epigenetic repression within these cell lines result in unpredictable and unstable transgene expression following random integration. Limited knowledge of the CHO genome and its higher order chromatin organization has thus far impeded functional genomics approaches required to tackle these issues. Here, we present an integrative three‐dimensional (3D) map of genome organization within the CHOK1SV® 10E9 cell line in conjunction with an improved, less fragmented CHOK1SV 10E9 genome assembly. Using our high‐resolution chromatin conformation datasets, we have assigned ≈90% of sequence to a chromosome‐scale genome assembly. Our genome‐wide 3D map identifies higher order chromatin structures such as topologically associated domains, incorporates our chromatin accessibility data to enhance the identification of active cis‐regulatory elements, and importantly links these cis‐regulatory elements to target promoters in a 3D promoter interactome. We demonstrate the power of our improved functional annotation by evaluating the 3D landscape of a transgene integration site and two phenotypically different cell lines. Our work opens up further novel genome engineering targets, has the potential to inform vital improvements for industrial biotherapeutic production, and represents a significant advancement for CHO cell line development.

GenAPI: a tool for gene absence-presence identification in fragmented bacterial genome sequences

BackgroundBacterial gene loss and acquisition is a well-known phenomenon which contributes to bacterial adaptation through changes in important phenotypes such as virulence, antibiotic resistance and metabolic capability. While advances in DNA sequencing have accelerated our ability to generate short genome sequence reads to disentangle phenotypic changes caused by gene loss and acquisition, the short-read genome sequencing often results in fragmented genome assemblies as a basis for identification of gene loss and acquisition events. However, sensitive and precise determination of gene content change for fragmented genome assemblies remains challenging as analysis needs to account for cases when only a fragment of the gene is assembled or when the gene assembly is split in more than one contig.ResultsWe developed GenAPI, a command-line tool that is designed to compare the gene content of bacterial genomes for which only fragmented genome assemblies are available. GenAPI, unlike other available tools of similar purpose, accounts for imperfections in sequencing and assembly, and aims to compensate for them. We tested the performance of GenAPI on three different datasets to show that GenAPI has a high sensitivity while it maintains precision when dealing with partly assembled genes in both simulated and real datasets. Furthermore, we benchmarked the performance of GenAPI with six popular tools for gene presence-absence identification.ConclusionsOur developed bioinformatics tool, called GenAPI, has the same precision and recall rates when analyzing complete genome sequences as the other tools of the same purpose; however, GenAPI’s performance is markedly better on fragmented genome assemblies.

Spatial inter-centromeric interactions facilitated the emergence of evolutionary new centromeres.

Centromeres of Candida albicans form on unique and different DNA sequences but a closely related species, Candida tropicalis, possesses homogenized inverted repeat (HIR)-associated centromeres. To investigate the mechanism of centromere type transition, we improved the fragmented genome assembly and constructed a chromosome-level genome assembly of C. tropicalis by employing PacBio sequencing, chromosome conformation capture sequencing (3C-seq), chromoblot, and genetic analysis of engineered aneuploid strains. Further, we analyzed the 3D genome organization using 3C-seq data, which revealed spatial proximity among the centromeres as well as telomeres of seven chromosomes in C. tropicalis. Intriguingly, we observed evidence of inter-centromeric translocations in the common ancestor of C. albicans and C. tropicalis. Identification of putative centromeres in closely related Candida sojae, Candida viswanathii and Candida parapsilosis indicates loss of ancestral HIR-associated centromeres and establishment of evolutionary new centromeres (ENCs) in C. albicans. We propose that spatial proximity of the homologous centromere DNA sequences facilitated karyotype rearrangements and centromere type transitions in human pathogenic yeasts of the CUG-Ser1 clade.

A chromosome-level assembly of the Atlantic herring genome-detection of a supergene and other signals of selection.

The Atlantic herring is a model species for exploring the genetic basis for ecological adaptation, due to its huge population size and extremely low genetic differentiation at selectively neutral loci. However, such studies have so far been hampered because of a highly fragmented genome assembly. Here, we deliver a chromosome-level genome assembly based on a hybrid approach combining a de novo Pacific Biosciences (PacBio) assembly with Hi-C-supported scaffolding. The assembly comprises 26 autosomes with sizes ranging from 12.4 to 33.1 Mb and a total size, in chromosomes, of 726 Mb, which has been corroborated by a high-resolution linkage map. A comparison between the herring genome assembly with other high-quality assemblies from bony fishes revealed few inter-chromosomal but frequent intra-chromosomal rearrangements. The improved assembly facilitates analysis of previously intractable large-scale structural variation, allowing, for example, the detection of a 7.8-Mb inversion on Chromosome 12 underlying ecological adaptation. This supergene shows strong genetic differentiation between populations. The chromosome-based assembly also markedly improves the interpretation of previously detected signals of selection, allowing us to reveal hundreds of independent loci associated with ecological adaptation.