A High-Quality Blue Whale Genome, Segmental Duplications, and Historical Demography.

The blue whale, Balaenoptera musculus, is the largest animal known to have ever existed, making it an important case study in longevity and resistance to cancer. To further this and other blue whale-related research, we report a reference-quality, long-read-based genome assembly of this fascinating species. We assembled the genome from PacBio long reads and utilized Illumina/10×, optical maps, and Hi-C data for scaffolding, polishing, and manual curation. We also provided long read RNA-seq data to facilitate the annotation of the assembly by NCBI and Ensembl. Additionally, we annotated both haplotypes using TOGA and measured the genome size by flow cytometry. We then compared the blue whale genome with other cetaceans and artiodactyls, including vaquita (Phocoena sinus), the world's smallest cetacean, to investigate blue whale's unique biological traits. We found a dramatic amplification of several genes in the blue whale genome resulting from a recent burst in segmental duplications, though the possible connection between this amplification and giant body size requires further study. We also discovered sites in the insulin-like growth factor-1 gene correlated with body size in cetaceans. Finally, using our assembly to examine the heterozygosity and historical demography of Pacific and Atlantic blue whale populations, we found that the genomes of both populations are highly heterozygous and that their genetic isolation dates to the last interglacial period. Taken together, these results indicate how a high-quality, annotated blue whale genome will serve as an important resource for biology, evolution, and conservation research.

Near chromosome-level and highly repetitive genome assembly of the snake pipefish Entelurus aequoreus (Syngnathiformes: Syngnathidae).

The snake pipefish, Entelurus aequoreus (Linnaeus, 1758), is a northern Atlantic fish inhabiting open seagrass environments that recently expanded its distribution range. Here, we present a highly contiguous, near chromosome-scale genome of E. aequoreus. The final assembly spans 1.6 Gbp in 7,391 scaffolds, with a scaffold N50 of 62.3 Mbp and L50 of 12. The 28 largest scaffolds (>21 Mbp) span 89.7% of the assembly length. A BUSCO completeness score of 94.1% and a mapping rate above 98% suggest a high assembly completeness. Repetitive elements cover 74.93% of the genome, one of the highest proportions identified in vertebrates. Our demographic modeling identified a peak in population size during the last interglacial period, suggesting the species might benefit from warmer water conditions. Our updated snake pipefish assembly is essential for future analyses of the morphological and molecular changes unique to the Syngnathidae.

The de novo genome of the Black-necked Snakefly (Venustoraphidia nigricollis Albarda, 1891): A resource to study the evolution of living fossils.

Snakeflies (Raphidioptera) are the smallest order of holometabolous insects that have kept their distinct and name-giving appearance since the Mesozoic, probably since the Jurassic, and possibly even since their emergence in the Carboniferous, more than 300 million years ago. Despite their interesting nature and numerous publications on their morphology, taxonomy, systematics, and biogeography, snakeflies have never received much attention from the general public, and only a few studies were devoted to their molecular biology. Due to this lack of molecular data, it is therefore unknown, if the conserved morphological nature of these living fossils translates to conserved genomic structures. Here, we present the first genome of the species and of the entire order of Raphidioptera. The final genome assembly has a total length of 669 Mbp and reached a high continuity with an N50 of 5.07 Mbp. Further quality controls also indicate a high completeness and no meaningful contamination. The newly generated data was used in a large-scaled phylogenetic analysis of snakeflies using shared orthologous sequences. Quartet score and gene concordance analyses revealed high amounts of conflicting signals within this group that might speak for substantial incomplete lineage sorting and introgression after their presumed re-radiation after the asteroid impact 66 million years ago. Overall, this reference genome will be a door-opening dataset for many future research applications, and we demonstrated its utility in a phylogenetic analysis that provides new insights into the evolution of this group of living fossils.

Acylated- and unacylated ghrelin during an oral glucose tolerance test in humans at risk for type 2 diabetes mellitus.

BACKGROUND/OBJECTIVES: The orexigenic peptide hormone ghrelin has been implicated in the pathophysiology of obesity and type 2 diabetes mellitus through its effects on nutrient homeostasis. Ghrelin is subject to a unique post-translational acyl modification regulating its biochemical activity. SUBJECTS/METHODS: In this study we aimed to investigate the relation of acylated (AcG) as well as unacylated ghrelin (UnG) with body weight and insulin resistance in the fasting (n = 545) and post-oral glucose tolerance test (oGTT) state (n = 245) in a metabolically well characterized cohort covering a broad range of BMI (17.95 kg/m²-76.25 kg/m²). RESULTS: Fasting AcG (median 94.2 pg/ml) and UnG (median 175.3 pg/ml) were negatively and the AcG/UnG ratio was positively correlated with BMI (all p < 0.0001). Insulin sensitivity (ISI) correlated positively with AcG (p = 0.0014) and UnG (p = 0.0004) but not with the AcG/UnG ratio. In a multivariate analysis, including ISI and BMI, only BMI, but not ISI was independently associated with AcG and UnG concentrations. Significant changes of AcG and UnG concentrations were detectable after oGTT stimulation, with slight decreases after 30 min and increases after 90-120 min. Subject stratification into BMI-divergent groups revealed more pronounced AcG increases in the two groups with BMI < 40 kg/m². CONCLUSION: Our data demonstrate lower concentrations for both AcG and UnG with increasing BMI as well as an increased proportion of the biologically active, acylated form of ghrelin giving point to pharmacologic intervention in ghrelin acylation and/or increase in UnG for treatment of obesity despite decreased absolute AcG levels.

The genome of the pygmy right whale illuminates the evolution of rorquals.

BACKGROUND: Baleen whales are a clade of gigantic and highly specialized marine mammals. Their genomes have been used to investigate their complex evolutionary history and to decipher the molecular mechanisms that allowed them to reach these dimensions. However, many unanswered questions remain, especially about the early radiation of rorquals and how cancer resistance interplays with their huge number of cells. The pygmy right whale is the smallest and most elusive among the baleen whales. It reaches only a fraction of the body length compared to its relatives and it is the only living member of an otherwise extinct family. This placement makes the pygmy right whale genome an interesting target to update the complex phylogenetic past of baleen whales, because it splits up an otherwise long branch that leads to the radiation of rorquals. Apart from that, genomic data of this species might help to investigate cancer resistance in large whales, since these mechanisms are not as important for the pygmy right whale as in other giant rorquals and right whales. RESULTS: Here, we present a first de novo genome of the species and test its potential in phylogenomics and cancer research. To do so, we constructed a multi-species coalescent tree from fragments of a whole-genome alignment and quantified the amount of introgression in the early evolution of rorquals. Furthermore, a genome-wide comparison of selection rates between large and small-bodied baleen whales revealed a small set of conserved candidate genes with potential connections to cancer resistance. CONCLUSIONS: Our results suggest that the evolution of rorquals is best described as a hard polytomy with a rapid radiation and high levels of introgression. The lack of shared positive selected genes between different large-bodied whale species supports a previously proposed convergent evolution of gigantism and hence cancer resistance in baleen whales.

Range-wide whole-genome resequencing of the brown bear reveals drivers of intraspecies divergence.

Population-genomic studies can shed new light on the effect of past demographic processes on contemporary population structure. We reassessed phylogeographical patterns of a classic model species of postglacial recolonisation, the brown bear (Ursus arctos), using a range-wide resequencing dataset of 128 nuclear genomes. In sharp contrast to the erratic geographical distribution of mtDNA and Y-chromosomal haplotypes, autosomal and X-chromosomal multi-locus datasets indicate that brown bear population structure is largely explained by recent population connectivity. Multispecies coalescent based analyses reveal cases where mtDNA haplotype sharing between distant populations, such as between Iberian and southern Scandinavian bears, likely results from incomplete lineage sorting, not from ancestral population structure (i.e., postglacial recolonisation). However, we also argue, using forward-in-time simulations, that gene flow and recombination can rapidly erase genomic evidence of former population structure (such as an ancestral population in Beringia), while this signal is retained by Y-chromosomal and mtDNA, albeit likely distorted. We further suggest that if gene flow is male-mediated, the information loss proceeds faster in autosomes than in X chromosomes. Our findings emphasise that contemporary autosomal genetic structure may reflect recent population dynamics rather than postglacial recolonisation routes, which could contribute to mtDNA and Y-chromosomal discordances.

A chromosome-scale reference genome assembly of the great sand eel, Hyperoplus lanceolatus.

Despite increasing sequencing efforts, numerous fish families still lack a reference genome, which complicates genetic research. One such understudied family is the sand lances (Ammodytidae, literally: "sand burrower"), a globally distributed clade of over 30 fish species that tend to avoid tidal currents by burrowing into the sand. Here, we present the first annotated chromosome-level genome assembly of the great sand eel (Hyperoplus lanceolatus). The genome assembly was generated using Oxford Nanopore Technologies long sequencing reads and Illumina short reads for polishing. The final assembly has a total length of 808.5 Mbp, of which 97.1% were anchored into 24 chromosome-scale scaffolds using proximity-ligation scaffolding. It is highly contiguous with a scaffold and contig N50 of 33.7 and 31.3 Mbp, respectively, and has a BUSCO completeness score of 96.9%. The presented genome assembly is a valuable resource for future studies of sand lances, as this family is of great ecological and commercial importance and may also contribute to studies aiming to resolve the suprafamiliar taxonomy of bony fishes.

Genomic Impact of Whaling in North Atlantic Fin Whales.

It is generally recognized that large-scale whaling in the 19th and 20th century led to a substantial reduction of the size of many cetacean populations, particularly those of the baleen whales (Mysticeti). The impact of these operations on genomic diversity of one of the most hunted whales, the fin whale (Balaenoptera physalus), has remained largely unaddressed because of the paucity of adequate samples and the limitation of applicable techniques. Here, we have examined the effect of whaling on the North Atlantic fin whale based on genomes of 51 individuals from Icelandic waters, representing three temporally separated intervals, 1989, 2009 and 2018 and provide a reference genome for the species. Demographic models suggest a noticeable drop of the effective population size of the North Atlantic fin whale around a century ago. The present results suggest that the genome-wide heterozygosity is not markedly reduced and has remained comparable with other baleen whale species. Similarly, there are no signs of apparent inbreeding, as measured by the proportion of long runs of homozygosity, or of a distinctively increased mutational load, as measured by the amount of putative deleterious mutations. Compared with other baleen whales, the North Atlantic fin whale appears to be less affected by anthropogenic influences than other whales such as the North Atlantic right whale, consistent with the presence of long runs of homozygosity and higher levels of mutational load in an otherwise more heterozygous genome. Thus, genome-wide assessments of other species and populations are essential for future, more specific, conservation efforts.

Two high-quality de novo genomes from single ethanol-preserved specimens of tiny metazoans (Collembola).

BACKGROUND: Genome sequencing of all known eukaryotes on Earth promises unprecedented advances in biological sciences and in biodiversity-related applied fields such as environmental management and natural product research. Advances in long-read DNA sequencing make it feasible to generate high-quality genomes for many non-genetic model species. However, long-read sequencing today relies on sizable quantities of high-quality, high molecular weight DNA, which is mostly obtained from fresh tissues. This is a challenge for biodiversity genomics of most metazoan species, which are tiny and need to be preserved immediately after collection. Here we present de novo genomes of 2 species of submillimeter Collembola. For each, we prepared the sequencing library from high molecular weight DNA extracted from a single specimen and using a novel ultra-low input protocol from Pacific Biosciences. This protocol requires a DNA input of only 5 ng, permitted by a whole-genome amplification step. RESULTS: The 2 assembled genomes have N50 values >5.5 and 8.5 Mb, respectively, and both contain ∼96% of BUSCO genes. Thus, they are highly contiguous and complete. The genomes are supported by an integrative taxonomy approach including placement in a genome-based phylogeny of Collembola and designation of a neotype for 1 of the species. Higher heterozygosity values are recorded in the more mobile species. Both species are devoid of the biosynthetic pathway for ß-lactam antibiotics known in several Collembola, confirming the tight correlation of antibiotic synthesis with the species way of life. CONCLUSIONS: It is now possible to generate high-quality genomes from single specimens of minute, field-preserved metazoans, exceeding the minimum contig N50 (1 Mb) required by the Earth BioGenome Project.

Education in the genomics era: Generating high-quality genome assemblies in university courses.

Recent advances in genome sequencing technologies have simplified the generation of genome data and reduced the costs for genome assemblies, even for complex genomes like those of vertebrates. More practically oriented genomic courses can prepare university students for the increasing importance of genomic data used in biological and medical research. Low-cost third-generation sequencing technology, along with publicly available data, can be used to teach students how to process genomic data, assemble full chromosome-level genomes, and publish the results in peer-reviewed journals, or preprint servers. Here we outline experiences gained from 2 master's-level courses and discuss practical considerations for teaching hands-on genome assembly courses.

Chromosome-level genome assembly of a benthic associated Syngnathiformes species: the common dragonet, Callionymus lyra.

Background: The common dragonet, Callionymus lyra, is one of three Callionymus species inhabiting the North Sea. All three species show strong sexual dimorphism. The males show strong morphological differentiation, e.g., species-specific colouration and size relations, while the females of different species have few distinguishing characters. Callionymus belongs to the 'benthic associated clade' of the order Syngnathiformes. The 'benthic associated clade' so far is not represented by genome data and serves as an important outgroup to understand the morphological transformation in 'long-snouted' syngnatiformes such as seahorses and pipefishes. Findings: Here, we present the chromosome-level genome assembly of C. lyra. We applied Oxford Nanopore Technologies' long-read sequencing, short-read DNBseq, and proximity-ligation-based scaffolding to generate a high-quality genome assembly. The resulting assembly has a contig N50 of 2.2 Mbp and a scaffold N50 of 26.7 Mbp. The total assembly length is 568.7 Mbp, of which over 538 Mbp were scaffolded into 19 chromosome-length scaffolds. The identification of 94.5% complete BUSCO genes indicates high assembly completeness. Additionally, we sequenced and assembled a multi-tissue transcriptome with a total length of 255.5 Mbp that was used to aid the annotation of the genome assembly. The annotation resulted in 19,849 annotated transcripts and identified a repeat content of 27.7%. Conclusions: The chromosome-level assembly of C. lyra provides a high-quality reference genome for future population genomic, phylogenomic, and phylogeographic analyses.

Comparison of the automated KRYPTOR chromogranin A assay with the DAKO ELISA.

BACKGROUND: Currently available assays for the widely used marker for neuroendocrine tumors chromogranin A have a weak comparability, involve manual work, and require batch processing of the samples. In this study, we evaluated the automated chromogranin A KRYPTOR assay compared with the widely used DAKO ELISA. METHODS: 83 samples were measured with the DAKO ELISA (EDTA plasma) and the KRYPTOR assay (serum), since different sample materials are recommended. Furthermore, different sample materials were compared. RESULTS: The results between the two assays were highly correlated (KRYPTOR serum (ng/mL) = 2.887 ELISA EDTA (U/L) + 5.028; r = 0.99). The inter-assay variation for the KRYPTOR assay was determined as 5.56% at low (95 ng/mL; n = 34) and 6.21% at high (530 ng/mL; n = 33) chromogranin A concentrations. EDTA plasma samples revealed significantly lower results than serum in the KRYPTOR assays and could not be used. CONCLUSIONS: In conclusion, the new chromogranin A KRYPTOR assay is fast, reliable, and compares well to an established test after adaptation of sample material and reported units.

Production and release of acylcarnitines by primary myotubes reflect the differences in fasting fat oxidation of the donors.

CONTEXT: Acylcarnitines are biomarkers of incomplete ß-oxidation and mitochondrial lipid overload but indicate also high rates of mitochondrial fatty acid oxidation. It is unknown whether the production of acylcarnitines in primary human myotubes obtained from lean, metabolically healthy subjects reflects the fat oxidation in vivo. OBJECTIVE: Our objective was to quantify the acylcarnitine production in myotubes obtained from subjects with low and high fasting respiratory quotient (RQ). METHODS: Fasting RQ was determined by indirect calorimetry. Muscle biopsies from the vastus lateralis muscle were taken from 6 subjects with low fasting RQ (mean 0.79 ± 0.03) and 6 with high fasting RQ (0.90 ± 0.03), and satellite cells were isolated, cultured, and differentiated to myotubes. Myotubes were cultivated with 125 µM (13)C-labeled palmitate for 30 minutes and 4 and 24 hours. Quantitative profiling of 42 intracellular and 31 extracellular acylcarnitines was performed by stable isotope dilution-based metabolomics analysis by liquid chromatography coupled to mass spectrometry. RESULTS: Myotubes from donors with high fasting RQ produced and released significant higher amounts of medium-chain acylcarnitines. High (13)C8 and (13)C10 acylcarnitine levels in the extracellular compartment correlated with high fasting RQ. The decreased expression of medium-chain acyl-coenzyme A dehydrogenase (MCAD) in these myotubes can explain the higher rate of incomplete fatty acid oxidation. A lower intracellular [(13)C]acetylcarnitine to carnitine and lower intracellular (13)C16/(13)C18 acylcarnitine to carnitine ratio indicate reduced fatty acid oxidation capacity in these myotubes. Mitochondrial DNA content was not different. CONCLUSION: Acylcarnitine production and release from primary human myotubes of donors with high fasting RQ indicate a reduced fatty acid oxidation capacity and a higher rate of incomplete fatty acid oxidation. Thus, quantitative profiling of acylcarnitine production in human myotubes can be a suitable tool to identify muscular determinants of fat oxidation in vivo.