A homology guide for Pacific salmon genus Oncorhynchus resolves patterns of ohnolog retention, resolution and local adaptation following the salmonid-specific whole-genome duplication event

Long-read de novo genome assembly of Gulf toadfish (Opsanus beta)

BACKGROUND

The family Batrachoididae are a group of ecologically important teleost fishes with unique life histories, behavior, and physiology that has made them popular model organisms. Batrachoididae remain understudied in the realm of genomics, with only four reference genome assemblies available for the family, with three being highly fragmented and not up to current assembly standards. Among these is the Gulf toadfish, Opsanus beta, a model organism for serotonin physiology which has recently been bred in captivity.

RESULTS

Here we present a new, de novo genome and transcriptome assemblies for the Gulf toadfish using PacBio long read technology. The genome size of the final assembly is 2.1 gigabases, which is among the largest teleost genomes. This new assembly improves significantly upon the currently available reference for Opsanus beta with a final scaffold count of 62, of which 23 are chromosome scale, an N50 of 98,402,768, and a BUSCO completeness score of 97.3%. Annotation with ab initio and transcriptome-based methods generated 41,076 gene models. The genome is highly repetitive, with ~ 70% of the genome composed of simple repeats and transposable elements. Satellite DNA analysis identified potential telomeric and centromeric regions.

CONCLUSIONS

This improved assembly represents a valuable resource for future research using this important model organism and to teleost genomics more broadly.

Activin Signaling Pathway Specialization During Embryonic and Skeletal Muscle Development in Rainbow Trout (Oncorhynchus mykiss)

Activin signaling is essential for proper embryonic, skeletal muscle, and reproductive development. Duplication of the pathway in teleost fish has enabled diversification of gene function across the pathway but how gene duplication influences the function of activin signaling in non-mammalian species is poorly understood. Full characterization of activin receptor signaling pathway expression was performed across embryonic development and during early skeletal muscle growth in rainbow trout (RBT, Oncorhynchus mykiss). Rainbow trout are a model salmonid species that have undergone two additional rounds of whole genome duplication. A small number of genes were expressed early in development and most genes increased expression throughout development. There was limited expression of activin Ab in RBT embryos despite these genes exhibiting significantly elevated expression in post-hatch skeletal muscle. CRISPR editing of the activin Aa1 ohnolog and subsequent production of meiotic gynogenetic offspring revealed that biallelic disruption of activin Aa1 did not result in developmental defects, as occurs with knockout of activin A in mammals. The majority of gynogenetic offspring exhibited homozygous activin Aa1 genotypes (wild type, in-frame, or frameshift) derived from the mosaic founder female. The research identifies mechanisms of specialization among the duplicated activin ohnologs across embryonic development and during periods of high muscle growth in larval and juvenile fish. The knowledge gained provides insights into potential viable gene-targeting approaches for engineering the activin receptor signaling pathway and establishes the feasibility of employing meiotic gynogenesis as a tool for producing homozygous F1 genome-edited fish for species with long-generation times, such as salmonids.

Gene Function is a Driver of Activin Signaling Pathway Evolution Following Whole-Genome Duplication in Rainbow Trout (Oncorhynchus mykiss)

The genomes of plant and animal species are influenced by ancestral whole-genome duplication (WGD) events, which have profound impacts on the regulation and function of gene networks. To gain insight into the consequences of WGD events, we characterized the sequence conservation and expression patterns of ohnologs in the highly duplicated activin receptor signaling pathway in rainbow trout (RBT). The RBT activin receptor signaling pathway is defined by tissue-specific expression of inhibitors and ligands and broad expression of receptors and Co-Smad signaling molecules. Signaling pathway ligands exhibited shared expression, while inhibitors and Smad signaling molecules primarily express a single dominant ohnolog. Our findings suggest that gene function influences ohnolog evolution following duplication of the activin signaling pathway in RBT.

Polyploid plants take cytonuclear perturbations in stride

Hybridization in plants is often accompanied by nuclear genome doubling (allopolyploidy), which has been hypothesized to perturb interactions between nuclear and organellar (mitochondrial and plastid) genomes by creating imbalances in the relative copy number of these genomes and producing genetic incompatibilities between maternally derived organellar genomes and the half of the allopolyploid nuclear genome from the paternal progenitor. Several evolutionary responses have been predicted to ameliorate these effects, including selection for changes in protein sequences that restore cytonuclear interactions; biased gene retention/expression/conversion favoring maternal nuclear gene copies; and fine-tuning of relative cytonuclear genome copy numbers and expression levels. Numerous recent studies, however, have found that evolutionary responses are inconsistent and rarely scale to genome-wide generalities. The apparent robustness of plant cytonuclear interactions to allopolyploidy may reflect features that are general to allopolyploids such as the lack of F2 hybrid breakdown under disomic inheritance, and others that are more plant-specific, including slow sequence divergence in organellar genomes and preexisting regulatory responses to changes in cell size and endopolyploidy during development. Thus, cytonuclear interactions may only rarely act as the main barrier to establishment of allopolyploid lineages, perhaps helping to explain why allopolyploidy is so pervasive in plant evolution.

Plasma proteomic profiling of bacterial cold water disease-resistant and -susceptible rainbow trout lines and biomarker discovery

Genetic variation for disease resistance is present in salmonid fish; however, the molecular basis is poorly understood, and biomarkers of disease susceptibility/resistance are unavailable. Previously, we selected a line of rainbow trout for high survival following standardized challenge with Flavobacterium psychrophilum (Fp), the causative agent of bacterial cold water disease. The resistant line (ARS-Fp-R) exhibits over 60 percentage points higher survival compared to a reference susceptible line (ARS-Fp-S). To gain insight into the differential host response between genetic lines, we compared the plasma proteomes from day 6 after intramuscular challenge. Pooled plasma from unhandled, PBS-injected, and Fp-injected groups were simultaneously analyzed using a TMT 6-plex label, and the relative abundance of 513 proteins was determined. Data are available via ProteomeXchange, with identifier PXD041308, and the relative protein abundance values were compared to mRNA measured from a prior, whole-body RNA-seq dataset. Our results identified a subset of differentially abundant intracellular proteins was identified, including troponin and myosin, which were not transcriptionally regulated, suggesting that these proteins were released into plasma following pathogen-induced tissue damage. A separate subset of high-abundance, secreted proteins were transcriptionally regulated in infected fish. The highest differentially expressed protein was a C1q family member (designated complement C1q-like protein 3; C1q-LP3) that was upregulated over 20-fold in the infected susceptible line while only modestly upregulated, 1.8-fold, in the infected resistant line. Validation of biomarkers was performed using immunoassays and C1q-LP3, skeletal muscle troponin C, cathelcidin 2, haptoglobin, leptin, and growth and differentiation factor 15 exhibited elevated concentration in susceptible line plasma. Complement factor H-like 1 exhibited higher abundance in the resistant line compared to the susceptible line in both control and challenged fish and thus was a baseline differentiator between lines. C1q-LP3 and STNC were elevated in Atlantic salmon plasma following experimental challenge with Fp. In summary, these findings further the understanding of the differential host response to Fp and identifies salmonid biomarkers that may have use for genetic line evaluation and on-farm health monitoring.