Advancing human disease research with fish evolutionary mutant models

Identification of AGO2 and PLEC genes polymorphisms in Hu sheep and their relationship with body size traits

The body size traits are major traits in livestock, which intuitively displays the development of the animal's bones and muscles. This study used PCR amplification, Sanger sequencing, KASPar genotyping, and quantitative real-time reverse transcription PCR (qRT-PCR) to analyze the Single-nucleotide polymorphism and expression characteristics of Argonaute RISC catalytic component 2 (AGO2) and Plectin (PLEC) genes in Hu sheep. Two intron mutations were found in Hu sheep, which were AGO2 g.51700 A > C and PLEC g.23157 C > T, respectively. Through association analysis of two mutation sites and body size traits, it was found that AGO2 g.51700 A > C mainly affects the chest and cannon circumference of Hu sheep of while PLEC g.23157 C mainly affects body height and body length. The combined genotypes of AGO2 and PLEC genes with body size traits showed SNPs at the AGO2 g.51700 A > C and PLEC g.23157 C > T loci significantly improved the body size traits of Hu sheep. In addition, the AGO2 gene has the highest expression levels in the heart, rumen, and tail fat, and the PLEC gene is highly expressed in the heart. These two loci can provide new research ideas for improving the body size traits of Hu sheep.

Phylogenomic analyses of all species of swordtail fishes (genus Xiphophorus) show that hybridization preceded speciation

Hybridization has been recognized to play important roles in evolution, however studies of the genetic consequence are still lagging behind in vertebrates due to the lack of appropriate experimental systems. Fish of the genus Xiphophorus are proposed to have evolved with multiple ancient and ongoing hybridization events. They have served as an informative research model in evolutionary biology and in biomedical research on human disease for more than a century. Here, we provide the complete genomic resource including annotations for all described 26 Xiphophorus species and three undescribed taxa and resolve all uncertain phylogenetic relationships. We investigate the molecular evolution of genes related to cancers such as melanoma and for the genetic control of puberty timing, focusing on genes that are predicted to be involved in pre-and postzygotic isolation and thus affect hybridization. We discovered dramatic size-variation of some gene families. These persisted despite reticulate evolution, rapid speciation and short divergence time. Finally, we clarify the hybridization history in the entire genus settling disputed hybridization history of two Southern swordtails. Our comparative genomic analyses revealed hybridization ancestries that are manifested in the mosaic fused genomes and show that hybridization often preceded speciation.

Transposable Element Expression Profiles in Premalignant Pigment Cell Lesions and Melanoma of Xiphophorus

Transposable elements (TEs) are characterized by their ability to change their genomic position. Through insertion or recombination leading to deletions and other chromosomal aberrations, they can cause genetic instability. The extent to which they thereby exert regulatory influence on cellular functions is unclear. To better characterize TEs in processes such as carcinogenesis, we used the well-established Xiphophorus melanoma model. By transcriptome sequencing, we show that an increasing total number in transposons correlates with progression of malignancy in melanoma samples from Xiphophorus interspecific hybrids. Further, by comparing the presence of TEs in the parental genomes of Xiphophorus maculatus and Xiphophorus hellerii, we could show that even in closely related species, genomic location and spectrum of TEs are considerably different.

High-throughput and genome-scale targeted mutagenesis using CRISPR in a nonmodel multicellular organism, Bombyx mori

Large-scale genetic mutant libraries are powerful approaches to interrogating genotype-phenotype correlations and identifying genes responsible for certain environmental stimuli, both of which are the central goal of life science study. We produced the first large-scale CRISPR-Cas9-induced library in a nonmodel multicellular organism, Bombyx mori We developed a piggyBac-delivered binary genome editing strategy, which can simultaneously meet the requirements of mixed microinjection, efficient multipurpose genetic operation, and preservation of growth-defect lines. We constructed a single-guide RNA (sgRNA) plasmid library containing 92,917 sgRNAs targeting promoters and exons of 14,645 protein-coding genes, established 1726 transgenic sgRNA lines following microinjection of 66,650 embryos, and generated 300 mutant lines with diverse phenotypic changes. Phenomic characterization of mutant lines identified a large set of genes responsible for visual phenotypic or economically valuable trait changes. Next, we performed pooled context-specific positive screens for tolerance to environmental pollutant cadmium exposure, and identified KWMTBOMO12902 as a strong candidate gene for breeding applications in sericulture industry. Collectively, our results provide a novel and versatile approach for functional B. mori genomics, as well as a powerful resource for identifying the potential of key candidate genes for improving various economic traits. This study also shows the effectiveness, practicality, and convenience of large-scale mutant libraries in other nonmodel organisms.

Phylogenomics analyses of all species of Swordtails (Genus Xiphophorus ) highlights hybridization precedes speciation

Hybridization has been recognized as an important driving force for evolution, however studies of the genetic consequence and its cause are still lagging behind in vertebrates due to the lack of appropriate experimental systems. Fish of the central American genus Xiphophorus were proposed to have evolved with multiple ancient and ongoing hybridization events, and served as a valuable research model in evolutionary biology and in biomedical research on human disease for more than a century. Here, we provide the complete genome resource and its annotation of all 26 Xiphophorus species. On this dataset we resolved the so far conflicting phylogeny. Through comparative genomic analyses we investigated the molecular evolution of genes related to melanoma, for a main sexually selected trait and for the genetic control of puberty timing, which are predicted to be involved in pre-and postzygotic isolation and thus to influence the probability of interspecific hybridization in Xiphophorus . We demonstrate dramatic size-variation of some gene families across species, despite the reticulate evolution and short divergence time. Finally, we clarify the hybridization history in the genus Xiphophorus genus, settle the long dispute on the hybridization origin of two Southern swordtails, highlight hybridizations precedes speciation, and reveal the distribution of hybridization ancestry remaining in the fused genome.

Cold-Driven Hemoglobin Evolution in Antarctic Notothenioid Fishes Prior to Hemoglobin Gene Loss in White-Blooded Icefishes

Expression of multiple hemoglobin isoforms with differing physiochemical properties likely helps species adapt to different environmental and physiological conditions. Antarctic notothenioid fishes inhabit the icy Southern Ocean and display fewer hemoglobin isoforms, each with less affinity for oxygen than temperate relatives. Reduced hemoglobin multiplicity was proposed to result from relaxed selective pressure in the cold, thermally stable, and highly oxygenated Antarctic waters. These conditions also permitted the survival and diversification of white-blooded icefishes, the only vertebrates living without hemoglobin. To understand hemoglobin evolution during adaptation to freezing water, we analyzed hemoglobin genes from 36 notothenioid genome assemblies. Results showed that adaptation to frigid conditions shaped hemoglobin gene evolution by episodic diversifying selection concomitant with cold adaptation and by pervasive evolution in Antarctic notothenioids compared to temperate relatives, likely a continuing adaptation to Antarctic conditions. Analysis of hemoglobin gene expression in adult hematopoietic organs in various temperate and Antarctic species further revealed a switch in hemoglobin gene expression underlying hemoglobin multiplicity reduction in Antarctic fish, leading to a single hemoglobin isoform in adult plunderfishes and dragonfishes, the sister groups to icefishes. The predicted high hemoglobin multiplicity in Antarctic fish embryos based on transcriptomic data, however, raises questions about the molecular bases and physiological implications of diverse hemoglobin isoforms in embryos compared to adults. This analysis supports the hypothesis that the last common icefish ancestor was vulnerable to detrimental mutations affecting the single ancestral expressed alpha- and beta-globin gene pair, potentially predisposing their subsequent loss.

Novel mitochondrial genome rearrangements including duplications and extensive heteroplasmy could underlie temperature adaptations in Antarctic notothenioid fishes

Mitochondrial genomes are known for their compact size and conserved gene order, however, recent studies employing long-read sequencing technologies have revealed the presence of atypical mitogenomes in some species. In this study, we assembled and annotated the mitogenomes of five Antarctic notothenioids, including four icefishes (Champsocephalus gunnari, C. esox, Chaenocephalus aceratus, and Pseudochaenichthys georgianus) and the cold-specialized Trematomus borchgrevinki. Antarctic notothenioids are known to harbor some rearrangements in their mt genomes, however the extensive duplications in icefishes observed in our study have never been reported before. In the icefishes, we observed duplications of the protein coding gene ND6, two transfer RNAs, and the control region with different copy number variants present within the same individuals and with some ND6 duplications appearing to follow the canonical Duplication-Degeneration-Complementation (DDC) model in C. esox and C. gunnari. In addition, using long-read sequencing and k-mer analysis, we were able to detect extensive heteroplasmy in C. aceratus and C. esox. We also observed a large inversion in the mitogenome of T. borchgrevinki, along with the presence of tandem repeats in its control region. This study is the first in using long-read sequencing to assemble and identify structural variants and heteroplasmy in notothenioid mitogenomes and signifies the importance of long-reads in resolving complex mitochondrial architectures. Identification of such wide-ranging structural variants in the mitogenomes of these fishes could provide insight into the genetic basis of the atypical icefish mitochondrial physiology and more generally may provide insights about their potential role in cold adaptation.

High resolution genomes of multiple Xiphophorus species provide new insights into microevolution, hybrid incompatibility, and epistasis

Because of diverged adaptative phenotypes, fish species of the genus Xiphophorus have contributed to a wide range of research for a century. Existing Xiphophorus genome assemblies are not at the chromosomal level and are prone to sequence gaps, thus hindering advancement of the intra- and inter-species differences for evolutionary, comparative, and translational biomedical studies. Herein, we assembled high-quality chromosome-level genome assemblies for three distantly related Xiphophorus species, namely, X. maculatus, X. couchianus, and X. hellerii Our overall goal is to precisely assess microevolutionary processes in the clade to ascertain molecular events that led to the divergence of the Xiphophorus species and to progress understanding of genetic incompatibility to disease. In particular, we measured intra- and inter-species divergence and assessed gene expression dysregulation in reciprocal interspecies hybrids among the three species. We found expanded gene families and positively selected genes associated with live bearing, a special mode of reproduction. We also found positively selected gene families are significantly enriched in nonpolymorphic transposable elements, suggesting the dispersal of these nonpolymorphic transposable elements has accompanied the evolution of the genes, possibly by incorporating new regulatory elements in support of the Britten-Davidson hypothesis. We characterized inter-specific polymorphisms, structural variants, and polymorphic transposable element insertions and assessed their association to interspecies hybridization-induced gene expression dysregulation related to specific disease states in humans.

Establishment, Long-Term Maintenance, and Characterization of Primary Liver Cells from Astyanax mexicanus

The tetra fish species Astyanax mexicanus comprises two morphotypes: cavefish that live in caves and surface fish that inhabit rivers and lakes. Because cavefish have adapted to the nutrient-poor conditions in their habitat whereas the surface fish populations can be used as a proxy for the ancestral condition, this species has become a powerful model system for understanding genetic variation underlying metabolic adaptation. The liver plays a critical role in glucose and fat metabolism in the body and hence is an important tissue for studying altered metabolism in health and disease. Cavefish morphs of A. mexicanus have been shown to develop fatty livers and exhibit massive differences in gene expression and chromatin architecture. Primary cell lines from various tissues have become invaluable tools for biochemical, toxicology, and cell biology experiments, as well as genetic and genomic analyses. To enhance the utility of the model system by enabling an expanded set of biochemical and in vitro experiments, we developed protocols for the isolation and maintenance of primary liver cells from A. mexicanus surface fish and cavefish. We also describe methods that can be used for primary cell characterization, including cloning, characterization of cell growth pattern, and lentivirus transduction. © 2023 Wiley Periodicals LLC. Basic Protocol 1: Primary culture of liver cells Support Protocol 1: Maintenance of A. mexicanus primary liver cells Support Protocol 2: Banking of A. mexicanus primary liver cells Support Protocol 3: Recovery of A. mexicanus primary liver cells Support Protocol 4: Primary liver cell cloning Support Protocol 5: Characterization of A. mexicanus primary liver cell growth pattern Basic Protocol 2: Lentiviral transduction of A. mexicanus primary liver cells.

Pathways to polar adaptation in fishes revealed by long-read sequencing

Long-read sequencing is driving a new reality for genome science in which highly contiguous assemblies can be produced efficiently with modest resources. Genome assemblies from long-read sequences are particularly exciting for understanding the evolution of complex genomic regions that are often difficult to assemble. In this study, we utilized long-read sequencing data to generate a high-quality genome assembly for an Antarctic eelpout, Ophthalmolycus amberensis, the first for the globally distributed family Zoarcidae. We used this assembly to understand how O. amberensis has adapted to the harsh Southern Ocean and compared it to another group of Antarctic fishes: the notothenioids. We showed that selection has largely acted on different targets in eelpouts relative to notothenioids. However, we did find some overlap; in both groups, genes involved in membrane structure, thermal tolerance and vision have evidence of positive selection. We found evidence for historical shifts of transposable element activity in O. amberensis and other polar fishes, perhaps reflecting a response to environmental change. We were specifically interested in the evolution of two complex genomic loci known to underlie key adaptations to polar seas: haemoglobin and antifreeze proteins (AFPs). We observed unique evolution of the haemoglobin MN cluster in eelpouts and related fishes in the suborder Zoarcoidei relative to other Perciformes. For AFPs, we identified the first species in the suborder with no evidence of afpIII sequences (Cebidichthys violaceus) in the genomic region where they are found in all other Zoarcoidei, potentially reflecting a lineage-specific loss of this cluster. Beyond polar fishes, our results highlight the power of long-read sequencing to understand genome evolution.

A novel gnotobiotic experimental system for Atlantic salmon (Salmo salar L.) reveals a microbial influence on mucosal barrier function and adipose tissue accumulation during the yolk sac stage

Gnotobiotic models have had a crucial role in studying the effect that commensal microbiota has on the health of their animal hosts. Despite their physiological and ecological diversity, teleost fishes are still underrepresented in gnotobiotic research. Moreover, a better understanding of host-microbe interactions in farmed fish has the potential to contribute to sustainable global food supply. We have developed a novel gnotobiotic experimental system that includes the derivation of fertilized eggs of farmed and wild Atlantic salmon, and gnotobiotic husbandry of fry during the yolk sac stage. We used a microscopy-based approach to estimate the barrier function of the skin mucus layer and used this measurement to select the derivation procedure that minimized adverse effects on the skin mucosa. We also used this method to demonstrate that the mucus barrier was reduced in germ-free fry when compared to fry colonized with two different bacterial communities. This alteration in the mucus barrier was preceded by an increase in the number of cells containing neutral mucosubstances in the anterior segment of the body, but without changes in the number of cells containing acidic substances in any of the other segments studied along the body axis. In addition, we showed how the microbial status of the fry temporarily affected body size and the utilization of internal yolk stores during the yolk sac stage. Finally, we showed that the presence of bacterial communities associated with the fry, as well as their composition, affected the size of adipose tissue. Fry colonized with water from a lake had a larger visceral adipose tissue depot than both conventionally raised and germ-free fry. Together, our results show that this novel gnotobiotic experimental system is a useful tool for the study of host-microbe interactions in this species of aquacultural importance.

Extreme intraspecific divergence in mitochondrial haplotypes makes the threespine stickleback fish an emerging evolutionary mutant model for mito-nuclear interactions

Mitochondrial DNA is primarily maternally inherited in most animals and evolves about 10 times faster than biparentally inherited nuclear DNA. Mitochondrial dysfunction (mt-dys) arises when interactions between the co-evolving mitochondrial and nuclear genomes are perturbed in essential processes like oxidative phosphorylation (OXPHOS). Over time mt-dys can lead to mitochondrial diseases (mt-diseases), which are surprisingly prevalent and include common diseases such as Alzheimer's, Parkinson's, and diabetes. Unfortunately, the strong impact that intraspecific mitochondrial and nuclear genetic variation has on mt-disease complicates its study and the development of effective treatments. Animal models have advanced our understanding of mt-disease but their relevance to human conditions is often limited by their relatively low nuclear genetic diversity. Many traditional laboratory models also typically have a single mitochondrial haplotype (mitotype), in stark contrast to over 5,000 mitotypes in humans worldwide. The threespine stickleback fish has an evolutionary history that has made it a favorable evolutionary mutant model (EMM) for studying mito-nuclear interactions and possibly mt-diseases. EMMs are species with naturally evolved states that mimic maladaptive human diseases. In threespine stickleback, a period of isolation followed by introgression of the mitochondrial genome from a sister species resulted in the maintenance of two distinct mitochondrial haplotypes which continue to segregate within many populations of wild stickleback. The existence of two mitogenomes segregating in numerous genetically diverse populations provides a unique system for exploring complex mito-nuclear dynamics. Here we provide the first complete coding region analysis of the two threespine stickleback mitotypes, whose mitogenomic divergence exceeds that of other mammalian models for mitochondrial disease and even that between ancient and modern humans. We find that divergence is not uniform across the mitogenome, but primarily impacts protein coding genes, and significantly impacts proteins in Complex I of OXPHOS. The full characterization of these highly divergent intraspecific mitotypes provides a foundation for the development of threespine stickleback as an EMM for mito-nuclear interactions.

Promoting validation and cross-phylogenetic integration in model organism research

Model organism (MO) research provides a basic understanding of biology and disease due to the evolutionary conservation of the molecular and cellular language of life. MOs have been used to identify and understand the function of orthologous genes, proteins, cells and tissues involved in biological processes, to develop and evaluate techniques and methods, and to perform whole-organism-based chemical screens to test drug efficacy and toxicity. However, a growing richness of datasets and the rising power of computation raise an important question: How do we maximize the value of MOs? In-depth discussions in over 50 virtual presentations organized by the National Institutes of Health across more than 10 weeks yielded important suggestions for improving the rigor, validation, reproducibility and translatability of MO research. The effort clarified challenges and opportunities for developing and integrating tools and resources. Maintenance of critical existing infrastructure and the implementation of suggested improvements will play important roles in maintaining productivity and facilitating the validation of animal models of human biology and disease.

A parasite outbreak in notothenioid fish in an Antarctic fjord

Climate changes can promote disease outbreaks, but their nature and potential impacts in remote areas have received little attention. In a hot spot of biodiversity on the West Antarctic Peninsula, which faces among the fastest changing climates on Earth, we captured specimens of two notothenioid fish species affected by large skin tumors at an incidence never before observed in the Southern Ocean. Molecular and histopathological analyses revealed that X-cell parasitic alveolates, members of a genus we call Notoxcellia, are the etiological agent of these tumors. Parasite-specific molecular probes showed that xenomas remained within the skin but largely outgrew host cells in the dermis. We further observed that tumors induced neovascularization in underlying tissue and detrimentally affected host growth and condition. Although many knowledge gaps persist about X-cell disease, including its mode of transmission and life cycle, these findings reveal potentially active biotic threats to vulnerable Antarctic ecosystems.

Fish as Model Systems to Study Epigenetic Drivers in Human Self-Domestication and Neurodevelopmental Cognitive Disorders

Modern humans exhibit phenotypic traits and molecular events shared with other domesticates that are thought to be by-products of selection for reduced aggression. This is the human self-domestication hypothesis. As one of the first types of responses to a novel environment, epigenetic changes may have also facilitated early self-domestication in humans. Here, we argue that fish species, which have been recently domesticated, can provide model systems to study epigenetic drivers in human self-domestication. To test this, we used in silico approaches to compare genes with epigenetic changes in early domesticates of European sea bass with genes exhibiting methylation changes in anatomically modern humans (comparison 1), and neurodevelopmental cognitive disorders considered to exhibit abnormal self-domestication traits, i.e., schizophrenia, Williams syndrome, and autism spectrum disorders (comparison 2). Overlapping genes in comparison 1 were involved in processes like limb morphogenesis and phenotypes like abnormal jaw morphology and hypopigmentation. Overlapping genes in comparison 2 affected paralogue genes involved in processes such as neural crest differentiation and ectoderm differentiation. These findings pave the way for future studies using fish species as models to investigate epigenetic changes as drivers of human self-domestication and as triggers of cognitive disorders.

Single-cell Iso-Sequencing enables rapid genome annotation for scRNAseq analysis

Single-cell RNA sequencing is a powerful technique that continues to expand across various biological applications. However, incomplete 3'-UTR annotations can impede single-cell analysis resulting in genes that are partially or completely uncounted. Performing single-cell RNA sequencing with incomplete 3'-UTR annotations can hinder the identification of cell identities and gene expression patterns and lead to erroneous biological inferences. We demonstrate that performing single-cell isoform sequencing in tandem with single-cell RNA sequencing can rapidly improve 3'-UTR annotations. Using threespine stickleback fish (Gasterosteus aculeatus), we show that gene models resulting from a minimal embryonic single-cell isoform sequencing dataset retained 26.1% greater single-cell RNA sequencing reads than gene models from Ensembl alone. Furthermore, pooling our single-cell sequencing isoforms with a previously published adult bulk Iso-Seq dataset from stickleback, and merging the annotation with the Ensembl gene models, resulted in a marginal improvement (+0.8%) over the single-cell isoform sequencing only dataset. In addition, isoforms identified by single-cell isoform sequencing included thousands of new splicing variants. The improved gene models obtained using single-cell isoform sequencing led to successful identification of cell types and increased the reads identified of many genes in our single-cell RNA sequencing stickleback dataset. Our work illuminates single-cell isoform sequencing as a cost-effective and efficient mechanism to rapidly annotate genomes for single-cell RNA sequencing.

Bone microstructure and bone mineral density are not systemically different in Antarctic icefishes and related Antarctic notothenioids

Ancestors of the Antarctic icefishes (family Channichthyidae) were benthic and had no swim bladder, making it energetically expensive to rise from the ocean floor. To exploit the water column, benthopelagic icefishes were hypothesized to have evolved a skeleton with “reduced bone,” which gross anatomical data supported. Here, we tested the hypothesis that changes to icefish bones also occurred below the level of gross anatomy. Histology and micro‐CT imaging of representative craniofacial bones (i.e., ceratohyal, frontal, dentary, and articular) of extant Antarctic fish species specifically evaluated two features that might cause the appearance of “reduced bone”: bone microstructure (e.g., bone volume fraction and structure linear density) and bone mineral density (BMD, or mass of mineral per volume of bone). Measures of bone microstructure were not consistently different in bones from the icefishes Chaenocephalus aceratus and Champsocephalus gunnari, compared to the related benthic notothenioids Notothenia coriiceps and Gobionotothen gibberifrons. Some quantitative measures, such as bone volume fraction and structure linear density, were significantly increased in some icefish bones compared to homologous bones of non‐icefish. However, such differences were rare, and no microstructural measures were consistently different in icefishes across all bones and species analyzed. Furthermore, BMD was similar among homologous bones of icefish and non‐icefish Antarctic notothenioids. In summary, “reduced bone” in icefishes was not due to systemic changes in bone microstructure or BMD, raising the prospect that “reduced bone” in icefish occurs only at the gross anatomic level (i.e., smaller or fewer bones). Given that icefishes exhibit delayed skeletal development compared to non‐icefish Antarctic fishes, combining these phenotypic data with genomic data might clarify genetic changes driving skeletal heterochrony.