Highly Continuous Genome Assembly of Eurasian Perch (Perca fluviatilis) Using Linked-Read Sequencing

Multi-genome comparisons reveal gain-and-loss evolution of anti-Mullerian hormone receptor type 2 as a candidate master sex-determining gene in Percidae

BACKGROUND

The Percidae family comprises many fish species of major importance for aquaculture and fisheries. Based on three new chromosome-scale assemblies in Perca fluviatilis, Perca schrenkii, and Sander vitreus along with additional percid fish reference genomes, we provide an evolutionary and comparative genomic analysis of their sex-determination systems.

RESULTS

We explored the fate of a duplicated anti-Mullerian hormone receptor type-2 gene (amhr2bY), previously suggested to be the master sex-determining (MSD) gene in P. flavescens. Phylogenetically related and structurally similar amhr2 duplicates (amhr2b) were found in P. schrenkii and Sander lucioperca, potentially dating this duplication event to their last common ancestor around 19-27 Mya. In P. fluviatilis and S. vitreus, this amhr2b duplicate has been likely lost while it was subject to amplification in S. lucioperca. Analyses of the amhr2b locus in P. schrenkii suggest that this duplication could be also male-specific as it is in P. flavescens. In P. fluviatilis, a relatively small (100 kb) non-recombinant sex-determining region (SDR) was characterized on chromosome 18 using population-genomics approaches. This SDR is characterized by many male-specific single-nucleotide variations (SNVs) and no large duplication/insertion event, suggesting that P. fluviatilis has a male heterogametic sex-determination system (XX/XY), generated by allelic diversification. This SDR contains six annotated genes, including three (c18h1orf198, hsdl1, tbc1d32) with higher expression in the testis than in the ovary.

CONCLUSIONS

Together, our results provide a new example of the highly dynamic sex chromosome turnover in teleosts and provide new genomic resources for Percidae, including sex-genotyping tools for all three known Perca species.

Multi-genome comparisons reveal gain-and-loss evolution of the anti-Mullerian hormone receptor type 2 gene, an old master sex determining gene, in Percidae

The Percidae family comprises many fish species of major importance for aquaculture and fisheries. Based on three new chromosome-scale assemblies in Perca fluviatilis, Perca schrenkii and Sander vitreus along with additional percid fish reference genomes, we provide an evolutionary and comparative genomic analysis of their sex-determination systems. We explored the fate of a duplicated anti-Mullerian hormone receptor type-2 gene (amhr2bY), previously suggested to be the master sex determining (MSD) gene in P. flavescens. Phylogenetically related and structurally similar amhr2 duplications (amhr2b) were found in P. schrenkii and Sander lucioperca, potentially dating this duplication event to their last common ancestor around 19-27 Mya. In P. fluviatilis and S. vitreus, this amhr2b duplicate has been lost while it was subject to amplification in S. lucioperca. Analyses of the amhr2b locus in P. schrenkii suggest that this duplication could be also male-specific as it is in P. flavescens. In P. fluviatilis, a relatively small (100 kb) non-recombinant sex-determining region (SDR) was characterized on chromosome-18 using population-genomics approaches. This SDR is characterized by many male-specific single-nucleotide variants (SNVs) and no large duplication/insertion event, suggesting that P. fluviatilis has a male heterogametic sex determination system (XX/XY), generated by allelic diversification. This SDR contains six annotated genes, including three (c18h1orf198, hsdl1, tbc1d32) with higher expression in testis than ovary. Together, our results provide a new example of the highly dynamic sex chromosome turnover in teleosts and provide new genomic resources for Percidae, including sex-genotyping tools for all three known Perca species.

Chromosome-level genome of butterflyfish unveils genomic features of unique colour patterns and morphological traits

Chaetodontidae, known as butterflyfishes, are typical fish in coral ecosystems, exhibiting remarkable interspecific differences including body colour patterns and feeding ecology. In this study, we report genomes of three butterflyfish species (Chelmon rostratus, Chaetodon trifasciatus and Chaetodon auriga) and a closely related species from the Pomacanthidae family, Centropyge bicolour, with an average genome size of 65,611 Mb. Chelmon rostratus, comprising 24 chromosomes assembled to the chromosome level, could be served as a reference genome for butterflyfish. By conducting a collinearity analysis between butterflyfishes and several fishes, we elucidated the specific and conserved genomic features of butterflyfish, with particular emphasis on novel genes arising from tandem duplications and their potential functions. In addition to the two melanocyte-specific tyr genes commonly found in fish, we found the gene tyrp3, a new tyrosinase-related proteins gene in the reef fish, including butterflyfish and clownfish, implicating their involvement in the pigmentation diversity of fish. Additionally, we observed a tandem duplication expansion of three copies of nell1 gene in C. rostratus genome, which likely contribute to its unique jaw development and distinctive morphology of its sharp mouth. These results provided valuable genomic resources for further investigations into the genetic diversity and evolutionary adaptations of reef fish.

Quality of DNA extracted from freshwater fish scales and mucus and its application in genetic diversity studies of Perca fluviatilis and Rutilus rutilus

Studies on genetic diversity require biological material containing a reliable source of DNA that can be extracted and analyzed. Recently, non-invasive sampling has become a preferred sampling method of biological material. The suitability of a less invasive approach that involves obtaining samples by swabbing the fish skin (including live, non-anesthetized fish) should be considered. In this study, we compared the efficiency of DNA extraction, amplification, and sequencing of mtDNA fragments of two fish species Perca fluviatilis and Rutilus rutilus based on DNA collected from the scales and mucus using the modified Aljanabi and Martinez method. The results revealed a higher quality of DNA extracted from the mucus; however, the mean DNA concentration obtained from the scales of both fish species was higher. We verified the method suitable for amplification and sequencing of mtDNA fragments of both fish species using newly designed markers (D-loop, ATP6) and examined the potential risk of intraspecific cross-contamination. The DNA sequence alignment analysis revealed identical sequences attributed to the same individual when DNA, extracted from two different sources (scales and mucus), was used. We demonstrated that the quantity and quality of DNA extracted from the scales and mucus using the proposed method were high enough to carry out genetic diversity studies based on sampling of live fish with the possibility to release it after collecting samples.

Near-Complete Genome Sequence of Lötschberg Virus (Mononegavirales: Filoviridae) Identified in European Perch (Perca fluviatilis Linnaeus, 1758)

We obtained the near-complete genome sequence of a novel virus, Lötschberg virus (LTBV), from a European perch metatranscriptome. Genome organization and pairwise sequence comparison indicated that LTBV represents a tentative new species and genus of the mononegaviral family Filoviridae.

Draft genome assembly for the colombian freshwater bocachico fish, Prochilodus magdalenae

We report the first draft genome assembly for Prochilodus magdalenae, the leading representative species of the Prochilodontidae family in Colombia. This 1.2-Gb assembly, with a GC content of 42.0% and a repetitive content of around 31.0%, is in the range of previously reported characid species genomes. Annotation identified 34,725 nuclear genes, and BUSCO completeness value was 94.9%. Gene ontology and primary metabolic pathway annotations indicate similar gene profiles for P. magdalenae and the closest species with annotated genomes: blind cave fish (Astyanax mexicanus) and red piranha (Pygocentrus nattereri). A comparative analysis showed similar genome traits to other characid species. The fully sequenced and annotated mitochondrial genome reproduces the taxonomic classification of P. magdalenae and confirms the low mitochondrial genetic divergence inside the Prochilodus genus. Phylogenomic analysis, using nuclear single-copy orthologous genes, also confirmed the evolutionary position of the species. This genome assembly provides a high-resolution genetic resource for sustainable P. magdalenae management in Colombia and, as the first genome assembly for the Prochilodontidae family, will contribute to fish genomics throughout South America.

In Silico Screening and Development of Microsatellite Markers for Genetic Analysis in Perca fluviatilis

Perca fluviatilis is an economically important species of freshwater fish. To understand the genetic structure of P. fluviatilis in China, 268 samples were collected from Wulungu Lake (WL), Jili Lake (JL), the Wulungu River (WR), and the Kalaeerqisi River (KR). These samples were then analyzed using microsatellite markers. A total of 98,425 microsatellite markers were developed based on the genomic data, and 29 polymorphic microsatellite markers were selected to analyze genetic diversity in this study. The number of alleles (Na) and observed heterozygosity (Ho) per population ranged from 4.621 (KR) to 11.172 (WL) and from 0.510 (KR) to 0.716 (JL), respectively. The results of the polymorphic information content (PIC) showed that the WL, JL, and WR populations were highly polymorphic (PIC≥ 0.5) and that the KR population was moderately polymorphic (0.25 ≤ PIC < 0.5). The genetic differentiation coefficient (Fst) among the four P. fluviatilis populations was 0.074, indicating moderate genetic differentiation among the populations in Xinjiang. The reason for the significant difference between the rivers and lakes could be the presence of a dam blocking the flow of P. fluviatilis. The development of microsatellite markers provides support for population genetics in the future. The evaluation of the genetic structure of P. fluviatilis in Xinjiang provides a reference for the reproduction and conservation of P. fluviatilis.

Whole-genome sequencing illuminates multifaceted targets of selection to humic substances in Eurasian perch

Extreme environments are inhospitable to the majority of species, but some organisms are able to survive in such hostile conditions due to evolutionary adaptations. For example, modern bony fishes have colonized various aquatic environments, including perpetually dark, hypoxic, hypersaline and toxic habitats. Eurasian perch (Perca fluviatilis) is among the few fish species of northern latitudes that is able to live in very acidic humic lakes. Such lakes represent almost "nocturnal" environments; they contain high levels of dissolved organic matter, which in addition to creating a challenging visual environment, also affects a large number of other habitat parameters and biotic interactions. To reveal the genomic targets of humic-associated selection, we performed whole-genome sequencing of perch originating from 16 humic and 16 clear-water lakes in northern Europe. We identified over 800,000 single nucleotide polymorphisms, of which >10,000 were identified as potential candidates under selection (associated with >3000 genes) using multiple outlier approaches. Our findings suggest that adaptation to the humic environment may involve hundreds of regions scattered across the genome. Putative signals of adaptation were detected in genes and gene families with diverse functions, including organism development and ion transportation. The observed excess of variants under selection in regulatory regions highlights the importance of adaptive evolution via regulatory elements, rather than via protein sequence modification. Our study demonstrates the power of whole-genome analysis to illuminate the multifaceted nature of humic adaptation and provides the foundation for further investigation of causal mutations underlying phenotypic traits of ecological and evolutionary importance.

Genomic analysis unveils mechanisms of northward invasion and signatures of plateau adaptation in the Asian house rat

The Asian house rat (AHR), Rattus tanezumi, has recently invaded the northern half of China. The AHR is a highly adaptive rat species that has also successfully conquered the Qinghai-Tibet Plateau (QTP) and replaced the brown rat (BR), R. norvegicus, at the edge of the QTP. Here, we assembled a draft genome of the AHR and explored the mechanisms of its northward invasion and the genetic basis underlying plateau adaptation in this species. Population genomic analyses revealed that the northwardly invasive AHRs consisted of two independent and genetically distinct populations which might result from multiple independent primary invasion events. One invasive population exhibited reduced genetic diversity and distinct population structure compared with its source population, while the other displayed preserved genetic polymorphisms and little genetic differentiation from its source population. Genes involved in G-protein coupled receptors and carbohydrate metabolism may contribute to the local adaptation of northern AHRs. In particular, RTN4 was identified as a key gene for AHRs in the QTP that favours adaptation to high-altitude hypoxia. Coincidently, the physiological performance and transcriptome profiles of hypoxia-exposed rats both showed better hypoxia adaptation in AHRs than in BRs that failed to colonize the heart of the QTP, which may have facilitated the replacement of the BR population by the invading AHRs at the edge of the QTP. This study provides profound insights into the multiple origins of the northwardly invasive AHR and the great tolerance to hypoxia in this species.

Effects of Garlic Allium sativum Powder on Nutrient Digestibility, Haematology, and Immune and Stress Responses in Eurasian Perch Perca fluviatilis Juveniles

The supplementation of fish diets with phytogenics can increase growth performance and can modulate immune system response. European perch Perca fluviatilis (initial weight 25.0 ± 0.4 g) were fed a diet including 0 (Control), 10 (G10), 20 (G20), and 30 (G30) g kg-1 garlic powder. No significant difference in the growth parameters and somatic indices were observed. Significantly higher fat digestibility was observed in G10 and G30 diets compared to in the control and G20 diets(p < 0.05). Significantly greater red blood cell and white blood cell counts were observed with the G10 diet (p < 0.05). Garlic significantly decreased serum cholesterol in all of the experimental groups. Serum albumin was significantly higher in the G10 and G20 diets (p < 0.05). Immediately after the overcrowding stress challenge, the garlic groups showed significantly higher cortisol levels than the control group, while no significant difference was observed in the glucose concentration among groups. At 1 h post-stress, all of the groups that had been fed a garlic-supplemented diet showed lower cortisol levels than the control group, and this trend was maintained at 6 and 24 h post stress (p < 0.05), and glucose level in all garlic groups was significantly lower than control (p < 0.05). Garlic at 10 g kg-1 in feed can improve apparent fat digestibility and selected blood parameters and can enhance resistance against high-density and net handling stress in Eurasian perch.

Chromosome-Level Assembly of the Atlantic Silverside Genome Reveals Extreme Levels of Sequence Diversity and Structural Genetic Variation

The levels and distribution of standing genetic variation in a genome can provide a wealth of insights about the adaptive potential, demographic history, and genome structure of a population or species. As structural variants are increasingly associated with traits important for adaptation and speciation, investigating both sequence and structural variation is essential for wholly tapping this potential. Using a combination of shotgun sequencing, 10x Genomics linked reads and proximity-ligation data (Chicago and Hi-C), we produced and annotated a chromosome-level genome assembly for the Atlantic silverside (Menidia menidia)-an established ecological model for studying the phenotypic effects of natural and artificial selection-and examined patterns of genomic variation across two individuals sampled from different populations with divergent local adaptations. Levels of diversity varied substantially across each chromosome, consistently being highly elevated near the ends (presumably near telomeric regions) and dipping to near zero around putative centromeres. Overall, our estimate of the genome-wide average heterozygosity in the Atlantic silverside is among the highest reported for a fish, or any vertebrate (1.32-1.76% depending on inference method and sample). Furthermore, we also found extreme levels of structural variation, affecting ∼23% of the total genome sequence, including multiple large inversions (> 1 Mb and up to 12.6 Mb) associated with previously identified haploblocks showing strong differentiation between locally adapted populations. These extreme levels of standing genetic variation are likely associated with large effective population sizes and may help explain the remarkable adaptive divergence among populations of the Atlantic silverside.

Genome assembly and annotation of the California harvester ant Pogonomyrmex californicus

The harvester ant genus Pogonomyrmex is endemic to arid and semiarid habitats and deserts of North and South America. The California harvester ant Pogonomyrmex californicus is the most widely distributed Pogonomyrmex species in North America. Pogonomyrmex californicus colonies are usually monogynous, i.e. a colony has one queen. However, in a few populations in California, primary polygyny evolved, i.e. several queens cooperate in colony founding after their mating flights and continue to coexist in mature colonies. Here, we present a genome assembly and annotation of P. californicus. The size of the assembly is 241 Mb, which is in agreement with the previously estimated genome size. We were able to annotate 17,889 genes in total, including 15,688 protein-coding ones with BUSCO (Benchmarking Universal Single-Copy Orthologs) completeness at a 95% level. The presented P. californicus genome assembly will pave the way for investigations of the genomic underpinnings of social polymorphism in the number of queens, regulation of aggression, and the evolution of adaptations to dry habitats.

Genomic signatures of the evolution of defence against its natural enemies in the poisonous and medicinal plant Datura stramonium (Solanaceae)

Tropane alkaloids and terpenoids are widely used in the medicine and pharmaceutic industry and evolved as chemical defenses against herbivores and pathogens in the annual herb Datura stramonium (Solanaceae). Here, we present the first draft genomes of two plants from contrasting environments of D. stramonium. Using these de novo assemblies, along with other previously published genomes from 11 Solanaceae species, we carried out comparative genomic analyses to provide insights on the genome evolution of D. stramonium within the Solanaceae family, and to elucidate adaptive genomic signatures to biotic and abiotic stresses in this plant. We also studied, in detail, the evolution of four genes of D. stramonium-Putrescine N-methyltransferase, Tropinone reductase I, Tropinone reductase II and Hyoscyamine-6S-dioxygenase-involved in the tropane alkaloid biosynthesis. Our analyses revealed that the genomes of D. stramonium show signatures of expansion, physicochemical divergence and/or positive selection on proteins related to the production of tropane alkaloids, terpenoids, and glycoalkaloids as well as on R defensive genes and other important proteins related with biotic and abiotic pressures such as defense against natural enemies and drought.

Draft Genome Assembly of the Freshwater Apex Predator Wels Catfish (Silurus glanis) Using Linked-Read Sequencing

The wels catfish (Silurus glanis) is one of the largest freshwater fish species in the world. This top predator plays a key role in ecosystem stability, and represents an iconic trophy-fish for recreational fishermen. S. glanis is also a highly valued species for its high-quality boneless flesh, and has been cultivated for over 100 years in Eastern and Central Europe. The interest in rearing S. glanis continues to grow; the aquaculture production of this species has almost doubled during the last decade. However, despite its high ecological, cultural and economic importance, the available genomic resources for S. glanis are very limited. To fulfill this gap we report a de novo assembly and annotation of the whole genome sequence of a female S. glanis. The linked-read based technology with 10X Genomics Chromium chemistry and Supernova assembler produced a highly continuous draft genome of S. glanis: ∼0.8Gb assembly (scaffold N₅₀ = 3.2 Mb; longest individual scaffold = 13.9 Mb; BUSCO completeness = 84.2%), which included 313.3 Mb of putative repeated sequences. In total, 21,316 protein-coding genes were predicted, of which 96% were annotated functionally from either sequence homology or protein signature searches. The highly continuous genome assembly will be an invaluable resource for aquaculture genomics, genetics, conservation, and breeding research of S. glanis.

Twelve quick steps for genome assembly and annotation in the classroom

Eukaryotic genome sequencing and de novo assembly, once the exclusive domain of well-funded international consortia, have become increasingly affordable, thus fitting the budgets of individual research groups. Third-generation long-read DNA sequencing technologies are increasingly used, providing extensive genomic toolkits that were once reserved for a few select model organisms. Generating high-quality genome assemblies and annotations for many aquatic species still presents significant challenges due to their large genome sizes, complexity, and high chromosome numbers. Indeed, selecting the most appropriate sequencing and software platforms and annotation pipelines for a new genome project can be daunting because tools often only work in limited contexts. In genomics, generating a high-quality genome assembly/annotation has become an indispensable tool for better understanding the biology of any species. Herein, we state 12 steps to help researchers get started in genome projects by presenting guidelines that are broadly applicable (to any species), sustainable over time, and cover all aspects of genome assembly and annotation projects from start to finish. We review some commonly used approaches, including practical methods to extract high-quality DNA and choices for the best sequencing platforms and library preparations. In addition, we discuss the range of potential bioinformatics pipelines, including structural and functional annotations (e.g., transposable elements and repetitive sequences). This paper also includes information on how to build a wide community for a genome project, the importance of data management, and how to make the data and results Findable, Accessible, Interoperable, and Reusable (FAIR) by submitting them to a public repository and sharing them with the research community.

Humic-acid-driven escape from eye parasites revealed by RNA-seq and target-specific metabarcoding

Background

Next generation sequencing (NGS) technologies are extensively used to dissect the molecular mechanisms of host-parasite interactions in human pathogens. However, ecological studies have yet to fully exploit the power of NGS as a rich source for formulating and testing new hypotheses.

Methods

We studied Eurasian perch (Perca fluviatilis) and its eye parasite (Trematoda, Diplostomidae) communities in 14 lakes that differed in humic content in order to explore host-parasite-environment interactions. We hypothesised that high humic content along with low pH would decrease the abundance of the intermediate hosts (gastropods), thus limiting the occurrence of diplostomid parasites in humic lakes. This hypothesis was initially invoked by whole eye RNA-seq data analysis and subsequently tested using PCR-based detection and a novel targeted metabarcoding approach.

Results

Whole eye transcriptome results revealed overexpression of immune-related genes and the presence of eye parasite sequences in RNA-seq data obtained from perch living in clear-water lakes. Both PCR-based and targeted-metabarcoding approach showed that perch from humic lakes were completely free from diplostomid parasites, while the prevalence of eye flukes in clear-water lakes that contain low amounts of humic substances was close to 100%, with the majority of NGS reads assigned to Tylodelphys clavata.

Conclusions

High intraspecific diversity of T. clavata indicates that massively parallel sequencing of naturally pooled samples represents an efficient and powerful strategy for shedding light on cryptic diversity of eye parasites. Our results demonstrate that perch populations in clear-water lakes experience contrasting eye parasite pressure compared to those from humic lakes, which is reflected by prevalent differences in the expression of immune-related genes in the eye. This study highlights the utility of NGS to discover novel host-parasite-environment interactions and provide unprecedented power to characterize the molecular diversity of cryptic parasites.

Transcriptional Changes in the Ovaries of Perch from Chernobyl

Fish have been highly exposed to radiation in freshwater systems after the Chernobyl Nuclear Power Plant (NPP) accident in 1986 and in freshwater and marine systems after the more recent Fukushima NPP accident in 2011. In the years after the accident, the radioactivity levels rapidly declined due to radioactive decay and environmental processes, but chronic lower dose exposures persisted. To gain insights into the long-term effects of environmental low dose radiation on fish ovaries development, a high-throughput transcriptomic approach including a de novo assembly was applied to different gonad phenotypes of female perch: developed gonads from reference lakes, developed/irradiated from medium contaminated lake, and both developed/irradiated and undeveloped from more highly contaminated lakes. This is the most comprehensive analysis to date of the gene responses in wildlife reproductive system to radiation. Some gene responses that were modulated in irradiated gonads were found to be involved in biological processes including cell differentiation and proliferation (ggnb2, mod5, rergl), cytoskeleton organization (k1C18, mtpn), gonad development (nell2, tcp4), lipid metabolism (ldah, at11b, nltp), reproduction (cyb5, cyp17A, ovos), DNA damage repair (wdhd1, rad51, hus1), and epigenetic mechanisms (dmap1). Identification of these genes provides a better understanding of the underlying molecular mechanisms underpinning the development of the gonad phenotypes of wild perch and how fish may respond to chronic exposure to radiation in their natural environment, though causal attribution of gene responses remains unclear in the undeveloped gonads.

Genome assembly of Chiococca alba uncovers key enzymes involved in the biosynthesis of unusual terpenoids

Chiococca alba (L.) Hitchc. (snowberry), a member of the Rubiaceae, has been used as a folk remedy for a range of health issues including inflammation and rheumatism and produces a wealth of specialized metabolites including terpenes, alkaloids, and flavonoids. We generated a 558 Mb draft genome assembly for snowberry which encodes 28,707 high-confidence genes. Comparative analyses with other angiosperm genomes revealed enrichment in snowberry of lineage-specific genes involved in specialized metabolism. Synteny between snowberry and Coffea canephora Pierre ex A. Froehner (coffee) was evident, including the chromosomal region encoding caffeine biosynthesis in coffee, albeit syntelogs of N-methyltransferase were absent in snowberry. A total of 27 putative terpene synthase genes were identified, including 10 that encode diterpene synthases. Functional validation of a subset of putative terpene synthases revealed that combinations of diterpene synthases yielded access to products of both general and specialized metabolism. Specifically, we identified plausible intermediates in the biosynthesis of merilactone and ribenone, structurally unique antimicrobial diterpene natural products. Access to the C. alba genome will enable additional characterization of biosynthetic pathways responsible for health-promoting compounds in this medicinal species.

Whole-genome sequencing of leopard coral grouper ( Plectropomus leopardus) and exploration of regulation mechanism of skin color and adaptive evolution

Leopard coral groupers belong to the Plectropomus genus of the Epinephelidae family and are important fish for coral reef ecosystems and the marine aquaculture industry. To promote future research of this species, a high-quality chromosome-level genome was assembled using PacBio sequencing and Hi-C technology. A 787.06 Mb genome was assembled, with 99.7% (784.57 Mb) of bases anchored to 24 chromosomes. The leopard coral grouper genome size was smaller than that of other groupers, which may be related to its ancient status among grouper species. A total of 22 317 protein-coding genes were predicted. This high-quality genome of the leopard coral grouper is the first genomic resource for Plectropomus and should provide a pivotal genetic foundation for further research. Phylogenetic analysis of the leopard coral grouper and 12 other fish species showed that this fish is closely related to the brown-marbled grouper. Expanded genes in the leopard coral grouper genome were mainly associated with immune response and movement ability, which may be related to the adaptive evolution of this species to its habitat. In addition, we also identified differentially expressed genes (DEGs) associated with carotenoid metabolism between red and brown-colored leopard coral groupers. These genes may play roles in skin color decision by regulating carotenoid content in these groupers.

CSA: A high-throughput chromosome-scale assembly pipeline for vertebrate genomes

BACKGROUND

Easy-to-use and fast bioinformatics pipelines for long-read assembly that go beyond the contig level to generate highly continuous chromosome-scale genomes from raw data remain scarce.

RESULT

Chromosome-Scale Assembler (CSA) is a novel computationally highly efficient bioinformatics pipeline that fills this gap. CSA integrates information from scaffolded assemblies (e.g., Hi-C or 10X Genomics) or even from diverged reference genomes into the assembly process. As CSA performs automated assembly of chromosome-sized scaffolds, we benchmark its performance against state-of-the-art reference genomes, i.e., conventionally built in a laborious fashion using multiple separate assembly tools and manual curation. CSA increases the contig lengths using scaffolding, local re-assembly, and gap closing. On certain datasets, initial contig N50 may be increased up to 4.5-fold. For smaller vertebrate genomes, chromosome-scale assemblies can be achieved within 12 h using low-cost, high-end desktop computers. Mammalian genomes can be processed within 16 h on compute-servers. Using diverged reference genomes for fish, birds, and mammals, we demonstrate that CSA calculates chromosome-scale assemblies from long-read data and genome comparisons alone. Even contig-level draft assemblies of diverged genomes are helpful for reconstructing chromosome-scale sequences. CSA is also capable of assembling ultra-long reads.

CONCLUSIONS

CSA can speed up and simplify chromosome-level assembly and significantly lower costs of large-scale family-level vertebrate genome projects.

Genomic Resources for Darters (Percidae: Etheostominae) Provide Insight into Postzygotic Barriers Implicated in Speciation

Comparative genomic approaches are increasingly being used to study the evolution of reproductive barriers in nonmodel species. Although numerous studies have examined prezygotic isolation in darters (Percidae), investigations into postzygotic barriers have remained rare due to long generation times and a lack of genomic resources. Orangethroat and rainbow darters naturally hybridize and provide a remarkable example of male-driven speciation via character displacement. Backcross hybrids suffer from high mortality, which appears to promote behavioral isolation in sympatry. To investigate the genomic architecture of postzygotic isolation, we used Illumina and PacBio sequencing to generate a chromosome-level, annotated assembly of the orangethroat darter genome and high-density linkage maps for orangethroat and rainbow darters. We also analyzed genome-wide RADseq data from wild-caught adults of both species and laboratory-generated backcrosses to identify genomic regions associated with hybrid incompatibles. Several putative chromosomal translocations and inversions were observed between orangethroat and rainbow darters, suggesting structural rearrangements may underlie postzygotic isolation. We also found evidence of selection against recombinant haplotypes and transmission ratio distortion in backcross hybrid genomes, providing further insight into the genomic architecture of genetic incompatibilities. Notably, regions with high levels of genetic divergence between species were enriched for genes associated with developmental and meiotic processes, providing strong candidates for postzygotic isolating barriers. These findings mark significant contributions to our understanding of the genetic basis of reproductive isolation between species undergoing character displacement. Furthermore, the genomic resources presented here will be instrumental for studying speciation in darters, the most diverse vertebrate group in North America.

Characterization of a Y-specific duplication/insertion of the anti-Mullerian hormone type II receptor gene based on a chromosome-scale genome assembly of yellow perch, Perca flavescens

Yellow perch, Perca flavescens, is an ecologically and economically important species native to a large portion of the northern United States and southern Canada and is also a promising candidate species for aquaculture. However, no yellow perch reference genome has been available to facilitate improvements in both fisheries and aquaculture management practices. By combining Oxford Nanopore Technologies long-reads, 10X Genomics Illumina short linked reads and a chromosome contact map produced with Hi-C, we generated a high-continuity chromosome-scale yellow perch genome assembly of 877.4 Mb. It contains, in agreement with the known diploid chromosome yellow perch count, 24 chromosome-size scaffolds covering 98.8% of the complete assembly (N50 = 37.4 Mb, L50 = 11). We also provide a first characterization of the yellow perch sex determination locus that contains a male-specific duplicate of the anti-Mullerian hormone type II receptor gene (amhr2by) inserted at the proximal end of the Y chromosome (chromosome 9). Using this sex-specific information, we developed a simple PCR genotyping assay which accurately differentiates XY genetic males (amhr2by+ ) from XX genetic females (amhr2by- ). Our high-quality genome assembly is an important genomic resource for future studies on yellow perch ecology, toxicology, fisheries and aquaculture research. In addition, characterization of the amhr2by gene as a candidate sex-determining gene in yellow perch provides a new example of the recurrent implication of the transforming growth factor beta pathway in fish sex determination, and highlights gene duplication as an important genomic mechanism for the emergence of new master sex determination genes.

Review of resistance to chronic ionizing radiation exposure under environmental conditions in multicellular organisms

Ionizing radiation resistance occurs among many phylogenetic groups and its mechanisms remain incompletely understood. Tolerances to acute and chronic irradiation do not always correlate because different mechanisms may be involved. The radioresistance phenomenon becomes even more complex in the field than in the laboratory because the effects of radioactive contamination on natural populations are intertwined with those of other factors, such as bioaccumulation of radionuclides, interspecific competition, seasonal variations in environmental conditions, and land use changes due to evacuation of humans from contaminated areas. Previous reviews of studies performed in radioactive sites like the Kyshtym, Chernobyl, and Fukushima accident regions, and of protracted irradiation experiments, often focused on detecting radiation effects at low doses in radiosensitive organisms. Here we review the literature with a different purpose: to identify organisms with high tolerance to chronic irradiation under environmental conditions, which maintained abundant populations and/or outcompeted more radiosensitive species at high dose rates. Taxa for which consistent evidence for radioresistance came from multiple studies conducted in different locations and at different times were found among plants (e.g. willow and birch trees, sedges), invertebrate and vertebrate animals (e.g. rotifers, some insects, crustaceans and freshwater fish). These organisms are not specialized "extremophiles", but tend to tolerate broad ranges of environmental conditions and stresses, have small genomes, reproduce quickly and/or disperse effectively over long distances. Based on these findings, resistance to radioactive contamination can be examined in a more broad context of chronic stress responses.

Domestication may affect the maternal mRNA profile in unfertilized eggs, potentially impacting the embryonic development of Eurasian perch (Perca fluviatilis)

Domestication is an evolutionary process during which we expect populations to progressively adapt to an environment controlled by humans. It is accompanied by genetic and presumably epigenetic changes potentially leading to modifications in the transcriptomic profile in various tissues. Reproduction is a key function often affected by this process in numerous species, regardless of the mechanism. The maternal mRNA in fish eggs is crucial for the proper embryogenesis. Our working hypothesis is that modifications of maternal mRNAs may reflect potential genetic and/or epigenetic modifications occurring during domestication and could have consequences during embryogenesis. Consequently, we investigated the trancriptomic profile of unfertilized eggs from two populations of Eurasian perch. These two populations differed by their domestication histories (F1 vs. F7+-at least seven generations of reproduction in captivity) and were genetically differentiated (FST = 0.1055, p<0.05). A broad follow up of the oogenesis progression failed to show significant differences during oogenesis between populations. However, the F1 population spawned earlier with embryos presenting an overall higher survivorship than those from the F7+ population. The transcriptomic profile of unfertilized eggs showed 358 differentially expressed genes between populations. In conclusion, our data suggests that the domestication process may influence the regulation of the maternal transcripts in fish eggs, which could in turn explain differences of developmental success.

The First Highly Contiguous Genome Assembly of Pikeperch (Sander lucioperca), an Emerging Aquaculture Species in Europe

The pikeperch (Sander lucioperca) is a fresh and brackish water Percid fish natively inhabiting the northern hemisphere. This species is emerging as a promising candidate for intensive aquaculture production in Europe. Specific traits like cannibalism, growth rate and meat quality require genomics based understanding, for an optimal husbandry and domestication process. Still, the aquaculture community is lacking an annotated genome sequence to facilitate genome-wide studies on pikeperch. Here, we report the first highly contiguous draft genome assembly of Sander lucioperca. In total, 413 and 66 giga base pairs of DNA sequencing raw data were generated with the Illumina platform and PacBio Sequel System, respectively. The PacBio data were assembled into a final assembly size of ~900 Mb covering 89% of the 1,014 Mb estimated genome size. The draft genome consisted of 1966 contigs ordered into 1,313 scaffolds. The contig and scaffold N50 lengths are 3.0 Mb and 4.9 Mb, respectively. The identified repetitive structures accounted for 39% of the genome. We utilized homologies to other ray-finned fishes, and ab initio gene prediction methods to predict 21,249 protein-coding genes in the Sander lucioperca genome, of which 88% were functionally annotated by either sequence homology or protein domains and signatures search. The assembled genome spans 97.6% and 96.3% of Vertebrate and Actinopterygii single-copy orthologs, respectively. The outstanding mapping rate (99.9%) of genomic PE-reads on the assembly suggests an accurate and nearly complete genome reconstruction. This draft genome sequence is the first genomic resource for this promising aquaculture species. It will provide an impetus for genomic-based breeding studies targeting phenotypic and performance traits of captive pikeperch.

A Multireference-Based Whole Genome Assembly for the Obligate Ant-Following Antbird, Rhegmatorhina melanosticta (Thamnophilidae)

Current generation high-throughput sequencing technology has facilitated the generation of more genomic-scale data than ever before, thus greatly improving our understanding of avian biology across a range of disciplines. Recent developments in linked-read sequencing (Chromium 10×) and reference-based whole-genome assembly offer an exciting prospect of more accessible chromosome-level genome sequencing in the near future. We sequenced and assembled a genome of the Hairy-crested Antbird (Rhegmatorhina melanosticta), which represents the first publicly available genome for any antbird (Thamnophilidae). Our objectives were to (1) assemble scaffolds to chromosome level based on multiple reference genomes, and report on differences relative to other genomes, (2) assess genome completeness and compare content to other related genomes, and (3) assess the suitability of linked-read sequencing technology for future studies in comparative phylogenomics and population genomics studies. Our R. melanosticta assembly was both highly contiguous (de novo scaffold N50 = 3.3 Mb, reference based N50 = 53.3 Mb) and relatively complete (contained close to 90% of evolutionarily conserved single-copy avian genes and known tetrapod ultraconserved elements). The high contiguity and completeness of this assembly enabled the genome to be successfully mapped to the chromosome level, which uncovered a consistent structural difference between R. melanosticta and other avian genomes. Our results are consistent with the observation that avian genomes are structurally conserved. Additionally, our results demonstrate the utility of linked-read sequencing for non-model genomics. Finally, we demonstrate the value of our R. melanosticta genome for future researchers by mapping reduced representation sequencing data, and by accurately reconstructing the phylogenetic relationships among a sample of thamnophilid species.