An improved genome assembly uncovers prolific tandem repeats in Atlantic cod

Abstract

BackgroundThe first Atlantic cod (Gadus morhua) genome assembly published in 2011 was one of the early genome assemblies exclusively based on high-throughput 454 pyrosequencing. Since then, rapid advances in sequencing technologies have led to a multitude of assemblies generated for complex genomes, although many of these are of a fragmented nature with a significant fraction of bases in gaps. The development of long-read sequencing and improved software now enable the generation of more contiguous genome assemblies.ResultsBy combining data from Illumina, 454 and the longer PacBio sequencing technologies, as well as integrating the results of multiple assembly programs, we have created a substantially improved version of the Atlantic cod genome assembly. The sequence contiguity of this assembly is increased fifty-fold and the proportion of gap-bases has been reduced fifteen-fold. Compared to other vertebrates, the assembly contains an unusual high density of tandem repeats (TRs). Indeed, retrospective analyses reveal that gaps in the first genome assembly were largely associated with these TRs. We show that 21% of the TRs across the assembly, 19% in the promoter regions and 12% in the coding sequences are heterozygous in the sequenced individual.ConclusionsThe inclusion of PacBio reads combined with the use of multiple assembly programs drastically improved the Atlantic cod genome assembly by successfully resolving long TRs. The high frequency of heterozygous TRs within or in the vicinity of genes in the genome indicate a considerable standing genomic variation in Atlantic cod populations, which is likely of evolutionary importance.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-3448-x) contains supplementary material, which is available to authorized users.

Switching on the light: using metagenomic shotgun sequencing to characterize the intestinal microbiome of Atlantic cod

“…0.20.4) (Schmieder and Edwards, 2011) (Supporting Information Table S15). PhiX, host and human sequences were removed by mapping reads to the phiX reference genome [GenBank: J02482.1], the Atlantic cod genome assembly (gadMor 2), (Tørresen et al, 2017) and a masked version of the human genome (HG19) (Genome Reference Consortium, 2009) using BWA (ver. 0.7.13) (Li and Durbin, 2009) or BBMap (ver.…”

A nanopore based chromosome-level assembly representing Atlantic cod from the Celtic Sea

“…Following this, sequence flanking each marker was 171 used to precisely position all genetic markers to contigs in the gadMor_Celtic 172 assembly using megablast (Altschul et al 1990), and thereby associate sequence 173 with linkage groups. This analysis revealed 2 chimeric contigs containing at 174 markers from each of different linkage groups that were selectively 'broken' 175 using alignments with the gadMor2 assembly (Torresen et al 2017 have directly compared reference genomes from different cod ecotypes to define 304 and confirm the underlying mechanism, or to locate the genomic regions 305 containing the inversion breakpoints or to define the exact complement of genes 306 they contain. We aligned the recently released gadMor3 assembly (NCBI 307 accession ID: GCF_902167405.1) constructed from a NEAC individual to our 308 gadMor_Celtic assembly using LASTZ (Harris 2007).…”

Single cell transcriptome profiling of the Atlantic cod immune system

“…We used the Drop-seq Core Computational Protocol (Namesh 2015) using STAR alignment to map the raw sequencing data to the most recent version of the Atlantic cod genome, gadMor3 (RefSeq accession GCF_902167405.1). A gene of interest, GATA-3, was present in gadMor2 (Torresen, Star et al 2017) but missing in gadMor3, so the GATA-3 gene sequence was manually added to the gadMor3 assembly fasta file. Reads were then grouped by cell barcode and the unique molecular identifiers (UMIs) for each gene counted, resulting in a digital expression matrix showing the number of transcripts per gene per cell.…”