Rapid genome reshaping by multiple-gene loss after whole-genome duplication in teleost fish suggested by mathematical modeling

Evolution of two-pore domain potassium channels and their gene expression in zebrafish embryos

BACKGROUND

The two-pore domain potassium (K2P) channels are a major type of potassium channels that maintain the cell membrane potential by conducting passive potassium leak currents independent of voltage change. They play prominent roles in multiple physiological processes, including neuromodulation, perception of pain, breathing and mood control, and response to volatile anesthetics. Mutations in K2P channels have been linked to many human diseases, such as neuronal and cardiovascular disorders and cancers. Significant progress has been made to understand their protein structures, physiological functions, and pharmacological modifiers. However, their expression and function during embryonic development remain largely unknown.

RESULTS

We employed the zebrafish model and identified 23 k2p genes using BLAST search and gene cloning. We first analyzed vertebrate K2P channel evolution by phylogenetic and syntenic analyses. Our data revealed that the six subtypes of the K2P genes have already evolved in invertebrates long before the emergence of vertebrates. Moreover, the vertebrate K2P gene number increased, most likely due to two whole-genome duplications. Furthermore, we examined zebrafish k2p gene expression during early embryogenesis by in situ hybridization. Each subgroup's genes showed similar but distinct gene expression domains with some exceptions. Most of them were expressed in neural tissues consistent with their known function of neural excitability regulation. However, a few k2p genes were expressed temporarily in specific tissues or organs, suggesting that these K2P channels may be needed for embryonic development.

CONCLUSIONS

Our phylogenetic and developmental analyses of K2P channels shed light on their evolutionary history and potential roles during embryogenesis related to their physiological functions and human channelopathies.

Molecular Dating of the Teleost Whole Genome Duplication (3R) Is Compatible With the Expectations of Delayed Rediploidization

Vertebrate evolution has been punctuated by three whole genome duplication events that have been implicated causally in phenotypic evolution, from the origin of phenotypic novelties to explosive diversification. Arguably, the most dramatic of these is the 3R whole genome duplication event associated with the origin of teleost fishes which comprise more than half of all living vertebrate species. However, tests of a causal relationship between whole genome duplication and teleost diversification have proven difficult due to the challenge of establishing the timing of these phenomena. Here we show, based on molecular clock dating of concatenated gene alignments, that the 3R whole genome duplication event occurred in the early-middle Permian (286.18 to 267.20 million years ago; Ma), 52.02 to 12.84 million years (Myr) before the divergence of crown-teleosts in the latest Permian-earliest Late Triassic (254.36 to 234.16 Ma) and long before the major pulses of teleost diversification in Ostariophysi and Percomorpha (56.37 to 100.17 Myr and at least 139.24 to 183.29 Myr later, respectively). The extent of this temporal gap between putative cause and effect precludes 3R as a deterministic driver of teleost diversification. However, these age constraints remain compatible with the expectations of a prolonged rediploidization process following whole genome duplication which, through the effects of chromosome rearrangement and gene loss, remains a viable mechanism to explain the evolution of teleost novelties and diversification.

Plant genome evolution in the genus Eucalyptus is driven by structural rearrangements that promote sequence divergence

Genomes have a highly organized architecture (nonrandom organization of functional and nonfunctional genetic elements within chromosomes) that is essential for many biological functions, particularly gene expression and reproduction. Despite the need to conserve genome architecture, a high level of structural variation has been observed within species. As species separate and diverge, genome architecture also diverges, becoming increasingly poorly conserved as divergence time increases. However, within plant genomes, the processes of genome architecture divergence are not well described. Here we use long-read sequencing and de novo assembly of 33 phylogenetically diverse, wild and naturally evolving Eucalyptus species, covering 1-50 million years of diverging genome evolution to measure genome architectural conservation and describe architectural divergence. The investigation of these genomes revealed that following lineage divergence, genome architecture is highly fragmented by rearrangements. As genomes continue to diverge, the accumulation of mutations and the subsequent divergence beyond recognition of rearrangements become the primary driver of genome divergence. The loss of syntenic regions also contribute to genome divergence but at a slower pace than that of rearrangements. We hypothesize that duplications and translocations are potentially the greatest contributors to Eucalyptus genome divergence.

Single-cell transcriptomic analysis reveals genome evolution in predatory litostomatean ciliates

Many ciliated protists prey on other large microbial organisms, including other protists and microscopic metazoans. The ciliate class Litostomatea unites both predatory and endosymbiotic species. The evolution of predation ability in ciliates remains poorly understood, in part, due to a lack of genomic data. To fill this gap, we acquired the transcriptome profiles of six predatory litostomateans using single-cell sequencing technology and investigated their transcriptomic features. Our results show that: (1) in contrast to non-predatory ciliates, the predatory litostomateans have expanded gene families associated with transmembrane activity and reactive oxidative stress response pathways, potentially as a result of cellular behaviors such as fast contraction and extension; (2) the expansion of the calcium-activated BK potassium channel gene family, which hypothetically regulates cell contractility, is an ancient evolutionary event for the class Litostomatea, suggesting a rewired metabolism associated with the hunting behavior of predatory ciliates; and (3) three whole genome duplication (WGD) events have been detected in litostomateans, with genes associated with biosynthetic processes, transmembrane activity, and calcium-activated potassium channel activity being retained during the WGD events. In addition, we explored the evolutionary relationships among 17 ciliate species, including eight litostomateans, and provided a rich foundational dataset for future in-depth phylogenomic studies of Litostomatea. Our comprehensive analyses suggest that the rewired cellular metabolism via expanded gene families and WGD events might be the potential genetic basis for the predation ability of raptorial ciliates.

Long-lasting redundant gnrh1/3 expression in GnRH neurons enabled apparent switching of paralog usage during evolution

Expressed subtype of paralogous genes in functionally homologous cells sometimes show differences across species, the reasons for which have not been explained. The present study examined hypophysiotropic gonadotropin-releasing hormone (GnRH) neurons in vertebrates to investigate this mechanism. These neurons express either gnrh1 or gnrh3 paralogs, depending on the species, and apparent switching of the expressed paralogs in them occurred at least four times in vertebrate evolution. First, we found redundant expression of gnrh1 and gnrh3 in a single neuron in piranha and hypothesized that it may represent an ancestral GnRH system. Moreover, the gnrh1/gnrh3 enhancer of piranha induced reporter RFP/GFP co-expression in a single hypophysiotropic GnRH neuron in both zebrafish and medaka, whose GnRH neurons only express either gnrh3 or gnrh1. Thus, we propose that redundant expression of gnrh1/3 of relatively recent common ancestors may be the key to apparent switching of the paralog usage among present-day species.

A comprehensive evolutionary scenario for the origin and neofunctionalization of the Drosophila speciation gene Odysseus (OdsH)

Odysseus (OdsH) was the first speciation gene described in Drosophila related to hybrid sterility in offspring of mating between Drosophila mauritiana and Drosophila simulans. Its origin is attributed to the duplication of the gene unc-4 in the subgenus Sophophora. By using a much larger sample of Drosophilidae species, we showed that contrary to what has been previously proposed, OdsH origin occurred 62 MYA. Evolutionary rates, expression, and transcription factor-binding sites of OdsH evidence that it may have rapidly experienced neofunctionalization in male sexual functions. Furthermore, the analysis of the OdsH peptide allowed the identification of mutations of D. mauritiana that could result in incompatibility in hybrids. In order to find if OdsH could be related to hybrid sterility, beyond Sophophora, we explored the expression of OdsH in Drosophila arizonae and Drosophila mojavensis, a pair of sister species with incomplete reproductive isolation. Our data indicated that OdsH expression is not atypical in their male-sterile hybrids. In conclusion, we have proposed that the origin of OdsH occurred earlier than previously proposed, followed by neofunctionalization. Our results also suggested that its role as a speciation gene might be restricted to D. mauritiana and D. simulans.

Exploring the role of polymorphic interspecies structural variants in reproductive isolation and adaptive divergence in Eucalyptus

Structural variations (SVs) play a significant role in speciation and adaptation in many species, yet few studies have explored the prevalence and impact of different categories of SVs. We conducted a comparative analysis of long-read assembled reference genomes of closely related Eucalyptus species to identify candidate SVs potentially influencing speciation and adaptation. Interspecies SVs can be either fixed differences or polymorphic in one or both species. To describe SV patterns, we employed short-read whole-genome sequencing on over 600 individuals of Eucalyptus melliodora and Eucalyptus sideroxylon, along with recent high-quality genome assemblies. We aligned reads and genotyped interspecies SVs predicted between species reference genomes. Our results revealed that 49,756 of 58,025 and 39,536 of 47,064 interspecies SVs could be typed with short reads in E. melliodora and E. sideroxylon, respectively. Focusing on inversions and translocations, symmetric SVs that are readily genotyped within both populations, 24 were found to be structural divergences, 2,623 structural polymorphisms, and 928 shared structural polymorphisms. We assessed the functional significance of fixed interspecies SVs by examining differences in estimated recombination rates and genetic differentiation between species, revealing a complex history of natural selection. Shared structural polymorphisms displayed enrichment of potentially adaptive genes. Understanding how different classes of genetic mutations contribute to genetic diversity and reproductive barriers is essential for understanding how organisms enhance fitness, adapt to changing environments, and diversify. Our findings reveal the prevalence of interspecies SVs and elucidate their role in genetic differentiation, adaptive evolution, and species divergence within and between populations.

Modeling familial and sporadic Parkinson's disease in small fishes

The establishment of animal models for Parkinson's disease (PD) has been challenging. Nevertheless, once established, they will serve as valuable tools for elucidating the causes and pathogenesis of PD, as well as for developing new strategies for its treatment. Following the recent discovery of a series of PD causative genes in familial cases, teleost fishes, including zebrafish and medaka, have often been used to establish genetic PD models because of their ease of breeding and gene manipulation, as well as the high conservation of gene orthologs. Some of the fish lines can recapitulate PD phenotypes, which are often more pronounced than those in rodent genetic models. In addition, a new experimental teleost fish, turquoise killifish, can be used as a sporadic PD model, because it spontaneously manifests age-dependent PD phenotypes. Several PD fish models have already made significant contributions to the discovery of novel PD pathological features, such as cytosolic leakage of mitochondrial DNA and pathogenic phosphorylation in α-synuclein. Therefore, utilizing various PD fish models with distinct degenerative phenotypes will be an effective strategy for identifying emerging facets of PD pathogenesis and therapeutic modalities.

Symmetric expression of ohnologs encoding conserved antiviral responses in tetraploid common carp suggest absence of subgenome dominance after whole genome duplication

Allopolyploids often experience subgenome dominance, with one subgenome showing higher levels of gene expression and greater gene retention. Here, we address the functionality of both subgenomes of allotetraploid common carp (Cyprinus carpio) by analysing a functional network of interferon-stimulated genes (ISGs) crucial in anti-viral immune defence. As an indicator of subgenome dominance we investigated retainment of a core set of ohnologous ISGs. To facilitate our functional genomic analysis a high quality genome was assembled (WagV4.0). Transcriptome data from an in vitro experiment mimicking a viral infection was used to infer ISG expression. Transcriptome analysis confirmed induction of 88 ISG ohnologs on both subgenomes. In both control and infected states, average expression of ISG ohnologs was comparable between the two subgenomes. Also, the highest expressing and most inducible gene copies of an ohnolog pair could be derived from either subgenome. We found no strong evidence of subgenome dominance for common carp.

Radiation impacts gene redundancy and biofilm regulation of cryoconite microbiomes in Northern Hemisphere glaciers

BACKGROUND

Glaciers harbor diverse microorganisms adapted to extreme conditions with high radiation, fluctuating temperature, and low nutrient availability. In glacial ecosystems, cryoconite granules are hotspots of microbial metabolic activity and could influences the biogeochemical cycle on glacier surface. Climate change could influence glacier dynamics by changing regional meteorological factors (e.g., radiation, precipitation, temperature, wind, and evaporation). Moreover, meteorological factors not only influence glacier dynamics but also directly or indirectly influence cryoconite microbiomes. However, the relationship of the meteorological factors and cryoconite microbiome are poorly understood.

RESULTS

Here, we collected 88 metagenomes from 26 glaciers distributed in the Northern Hemisphere with corresponding public meteorological data to reveal the relationship between meteorological factors and variation of cryoconite microbiome. Our results showed significant differences in taxonomic and genomic characteristics between cryoconite generalists and specialists. Additionally, we found that the biogeography of both generalists and specialists was influenced by solar radiation. Specialists with smaller genome size and lower gene redundancy were more abundant under high radiation stress, implying that streamlined genomes are more adapted to high radiation conditions. Network analysis revealed that biofilm regulation is a ubiquitous function in response to radiation stress, and hub genes were associated with the formation and dispersion of biofilms.

CONCLUSION

These findings enhance our understanding of glacier cryoconite microbiome variation on a hemispheric scale and indicate the response mechanisms to radiation stress, which will support forecasts of the ecological consequences of future climate change. Video Abstract.

Molecular and functional characterization of two granulocyte colony stimulating factors in goldfish (Carassius auratus L.)

Granulocyte colony-stimulating factor (GCSF) is a member of the hematopoietic growth factor family that acts primarily on neutrophils and neutrophilic precursors to promote cell proliferation and differentiation. Although multiple GCSF genes have been found in teleosts, knowledge of their functions during fish hematopoietic development is still limited. Here, we report for the first time the molecular and functional characterization of two goldfish GCSFs (gfGCSF-a and gfGCSF-b). The open reading frame (ORF) of the gfGCSF-a and gfGCSF-b cDNA transcript consisted respectively of 624 bp and 678 bp with its ORF encoding 207 and 225 amino acids (aa), with a 17 aa signal peptide for each gene and a conserved domain of the IL-6 superfamily. Treatment of goldfish head kidney leukocytes (HKLs) with LPS increased gfGCSF-a and gfGCSF-b mRNA expression levels, also exposure of HKLs to either heat-killed or live A. hydrophila, induced transcriptional upregulation of gfGCSF-a and gfGCSF-b levels. Recombinant gfGCSF-a and gfGCSF-b protein (rgGCSF-a and rgGCSF-b) induced a dose-dependent production of TNFα and IL-1β from goldfish neutrophils. In vitro experiments showed rgGCSF-a and rgGCSF-b differentially promoted the proliferation and differentiation of leukocytes in goldfish. Furthermore, treatment of HKLs with rgGCSF-a showed significant upregulation of mRNA levels of the hematopoietic transcription factor GATA2, Runx1, MafB, and cMyb, while gfGCSF-b induces not only all four transcriptional factors mentioned above but also CEBPα. Our results indicate that goldfish GCSF-a and GCSF-b are important regulators of neutrophil proliferation and differentiation, which could stimulate different stages and lineages of hematopoiesis.

Genomic and Transcriptional Profiles of Kelch-like (klhl) Gene Family in Polyploid Carassius Complex

Genome duplication supplies raw genetic materials and has been thought to be essential for evolutionary innovation and ecological adaptation. Here, we select Kelch-like (klhl) genes to study the evolution of the duplicated genes in the polyploid Carassius complex, including amphidiploid C. auratus and amphitriploid C. gibelio. Phylogenetic, chromosomal location and read coverage analyses indicate that most of Carassius klhl genes exhibit a 2:1 relationship with zebrafish orthologs and confirm two rounds of polyploidy, an allotetraploidy followed by an autotriploidy, occurred during Carassius evolution. The lineage-specific expansion and biased retention/loss of klhl genes are also found in Carassius. Transcriptome analyses across eight adult tissues and seven embryogenesis stages reveal varied expression dominance and divergence between the two species. The expression of klhls in response to Carassius herpesvirus 2 infection shows different expression changes corresponding to distinct herpesvirus resistances in three C. gibelio gynogenetic clones. Finally, we find that most C. gibelio klhl genes possess three alleles except eight genes that have lost one or two alleles due to genome rearrangement. The allele expression bias is prosperous for Cgklhl genes and varies during embryogenesis owning to the sequential expression manner of the alleles. The current study provides global insights into the genomic and transcriptional evolution of duplicated genes in a given superfamily resulting from multiple rounds of polyploidization.

Dynamics of Gene Loss following Ancient Whole-Genome Duplication in the Cryptic Paramecium Complex

Whole-genome duplications (WGDs) have shaped the gene repertoire of many eukaryotic lineages. The redundancy created by WGDs typically results in a phase of massive gene loss. However, some WGD-derived paralogs are maintained over long evolutionary periods, and the relative contributions of different selective pressures to their maintenance are still debated. Previous studies have revealed a history of three successive WGDs in the lineage of the ciliate Paramecium tetraurelia and two of its sister species from the Paramecium aurelia complex. Here, we report the genome sequence and analysis of 10 additional P. aurelia species and 1 additional out group, revealing aspects of post-WGD evolution in 13 species sharing a common ancestral WGD. Contrary to the morphological radiation of vertebrates that putatively followed two WGD events, members of the cryptic P. aurelia complex have remained morphologically indistinguishable after hundreds of millions of years. Biases in gene retention compatible with dosage constraints appear to play a major role opposing post-WGD gene loss across all 13 species. In addition, post-WGD gene loss has been slower in Paramecium than in other species having experienced genome duplication, suggesting that the selective pressures against post-WGD gene loss are especially strong in Paramecium. A near complete lack of recent single-gene duplications in Paramecium provides additional evidence for strong selective pressures against gene dosage changes. This exceptional data set of 13 species sharing an ancestral WGD and 2 closely related out group species will be a useful resource for future studies on Paramecium as a major model organism in the evolutionary cell biology.

Clustering pattern and evolution characteristic of microRNAs in grass carp (Ctenopharyngodon idella)

BACKGROUND

A considerable fraction of microRNAs (miRNAs) are highly conserved, and certain miRNAs correspond to genomic clusters. The clustering of miRNAs can be advantageous, possibly by allowing coordinated expression. However, little is known about the evolutionary forces responsible for the loss and acquisition of miRNA and miRNA clusters.

RESULTS

The results demonstrated that several novel miRNAs arose throughout grass carp evolution. Duplication and de novo production were critical strategies for miRNA cluster formation. Duplicates accounted for a smaller fraction of the expansion in the grass carp miRNA than de novo creation. Clustered miRNAs are more conserved and change slower, whereas unique miRNAs usually have high evolution rates and low expression levels. The expression level of miRNA expression in clusters is strongly correlated.

CONCLUSIONS

This study examines the genomic distribution, evolutionary background, and expression regulation of grass carp miRNAs. Our findings provide novel insights into the genesis and development of miRNA clusters in teleost.

A Chromosome-level assembly of the Japanese eel genome, insights into gene duplication and chromosomal reorganization

Japanese eels (Anguilla japonica) are commercially important species, harvested extensively for food. Currently, this and related species (American and European eels) are challenging to breed on a commercial basis. As a result, the wild stock is used for aquaculture. Moreover, climate change, habitat loss, water pollution, and altered ocean currents affect eel populations negatively. Accordingly, the International Union for Conservation of Nature lists Japanese eels as endangered and on its red list. Here we presented a high-quality genome assembly for Japanese eels and demonstrated that large chromosome reorganizations occurred in the events of third-round whole-genome duplications (3R-WRDs). Several chromosomal fusions and fissions have reduced the ancestral protochromosomal number of 25 to 19 in the Anguilla lineage. A phylogenetic analysis of the expanded gene families showed that the olfactory receptors (group δ and ζ genes) and voltage-gated Ca2+ channels expanded significantly. Both gene families are crucial for olfaction and neurophysiology. Additional tandem and proximal duplications occurred following 3R-WGD to acquire immune-related genes for an adaptive advantage against various pathogens. The Japanese eel assembly presented here can be used to study other Anguilla species relating to evolution and conservation.

Salmonidae Genome: Features, Evolutionary and Phylogenetic Characteristics

The salmon family is one of the most iconic and economically important fish families, primarily possessing meat of excellent taste as well as irreplaceable nutritional and biological value. One of the most common and, therefore, highly significant members of this family, the Atlantic salmon (Salmo salar L.), was not without reason one of the first fish species for which a high-quality reference genome assembly was produced and published. Genomic advancements are becoming increasingly essential in both the genetic enhancement of farmed salmon and the conservation of wild salmon stocks. The salmon genome has also played a significant role in influencing our comprehension of the evolutionary and functional ramifications of the ancestral whole-genome duplication event shared by all Salmonidae species. Here we provide an overview of the current state of research on the genomics and phylogeny of the various most studied subfamilies, genera, and individual salmonid species, focusing on those studies that aim to advance our understanding of salmonid ecology, physiology, and evolution, particularly for the purpose of improving aquaculture production. This review should make potential researchers pay attention to the current state of research on the salmonid genome, which should potentially attract interest in this important problem, and hence the application of new technologies (such as genome editing) in uncovering the genetic and evolutionary features of salmoniforms that underlie functional variation in traits of commercial and scientific importance.

An atlas of fish genome evolution reveals delayed rediploidization following the teleost whole-genome duplication

Teleost fishes are ancient tetraploids descended from an ancestral whole-genome duplication that may have contributed to the impressive diversification of this clade. Whole-genome duplications can occur via self-doubling (autopolyploidy) or via hybridization between different species (allopolyploidy). The mode of tetraploidization conditions evolutionary processes by which duplicated genomes return to diploid meiotic pairing, and subsequent genetic divergence of duplicated genes (cytological and genetic rediploidization). How teleosts became tetraploid remains unresolved, leaving a fundamental gap in the interpretation of their functional evolution. As a result of the whole-genome duplication, identifying orthologous and paralogous genomic regions across teleosts is challenging, hindering genome-wide investigations into their polyploid history. Here, we combine tailored gene phylogeny methodology together with a state-of-the-art ancestral karyotype reconstruction to establish the first high-resolution comparative atlas of paleopolyploid regions across 74 teleost genomes. We then leverage this atlas to investigate how rediploidization occurred in teleosts at the genome-wide level. We uncover that some duplicated regions maintained tetraploidy for more than 60 million years, with three chromosome pairs diverging genetically only after the separation of major teleost families. This evidence suggests that the teleost ancestor was an autopolyploid. Further, we find evidence for biased gene retention along several duplicated chromosomes, contradicting current paradigms that asymmetrical evolution is specific to allopolyploids. Altogether, our results offer novel insights into genome evolutionary dynamics following ancient polyploidizations in vertebrates.

auts2 Features and Expression Are Highly Conserved during Evolution Despite Different Evolutionary Fates Following Whole Genome Duplication

The AUTS2 gene plays major roles during brain development and is associated with various neuropathologies including autism. Data in non-mammalian species are scarce, and the aim of our study was to provide a comprehensive analysis of auts2 evolution in teleost fish, which are widely used for in vivo functional analysis and biomedical purposes. Comparative genomics in 78 species showed that auts2a and auts2b originate from the teleost-specific whole genome duplication (TGD). auts2a, which is highly similar to human AUTS2, was almost systematically retained following TGD. In contrast, auts2b, which encodes for a shorter protein similar to a short human AUTS2 isoform, was lost more frequently and independently during evolution. RNA-seq analysis in 10 species revealed a highly conserved profile with predominant expression of both genes in the embryo, brain, and gonads. Based on protein length, conserved domains, and expression profiles, we speculate that the long human isoform functions were retained by auts2a, while the short isoform functions were retained by auts2a and/or auts2b, depending on the lineage/species. auts2a showed a burst in expression during medaka brain formation, where it was expressed in areas of the brain associated with neurodevelopmental disorders. Together, our data suggest a strong conservation of auts2 functions in vertebrates despite different evolutionary scenarios in teleosts.

Comparative Study of the Molecular Characterization, Evolution, and Structure Modeling of Digestive Lipase Genes Reveals the Different Evolutionary Selection Between Mammals and Fishes

Vertebrates need suitable lipases to digest lipids for the requirement of energy and essential nutrients; however, the main digestive lipase genes of fishes have certain controversies. In this study, two types of digestive lipase genes (pancreatic lipase (pl) and bile salt-activated lipase (bsal)) were identified in mammals and fishes. The neighborhood genes and key active sites of the two lipase genes were conserved in mammals and fishes. Three copies of PL genes were found in mammals, but only one copy of the pl gene was found in most of the fish species, and the pl gene was even completely absent in some fish species (e.g., zebrafish, medaka, and common carp). Additionally, the hydrophobic amino acid residues (Ile and Leu) which are important to pancreatic lipase activity were also absent in most of the fish species. The PL was the main digestive lipase gene in mammals, but the pl gene seemed not to be the main digestive lipase gene in fish due to the absence of the pl gene sequence and the important amino acid residues. In contrast, the bsal gene existed in all fish species, even two to five copies of bsal genes were found in most of the fishes, but only one copy of the BSAL gene was found in mammals. The amino acid residues of bile salt-binding sites and the three-dimensional (3D) structure modeling of Bsal proteins were conserved in most of the fish species, so bsal might be the main digestive lipase gene in fish. The phylogenetic analysis also indicated that pl or bsal showed an independent evolution between mammals and fishes. Therefore, we inferred that the evolutionary selection of the main digestive lipase genes diverged into two types between mammals and fishes. These findings will provide valuable evidence for the study of lipid digestion in fish.

A chromosome-level genome assembly of chia provides insights into high omega-3 content and coat color variation of its seeds

Chia (Salvia hispanica) is a functional food crop for humans. Although its seeds contain high omega-3 fatty acids, the seed yield of chia is still low. Genomic resources available for this plant are limited. We report the first high-quality chromosome-level genome sequence of chia. The assembled genome size was 347.6 Mb and covered 98.1% of the estimated genome size. A total of 31 069 protein-coding genes were predicted. The absence of recent whole-genome duplication and the relatively low intensity of transposable element expansion in chia compared to its sister species contribute to its small genome size. Transcriptome sequencing and gene duplication analysis reveal that the expansion of the fab2 gene family is likely to be related to the high content of omega-3 in seeds. The white seed coat color is determined by a single locus on chromosome 4. This study provides novel insights into the evolution of Salvia species and high omega-3 content, as well as valuable genomic resources for genetic improvement of important commercial traits of chia and its related species.

Species Tree Estimation and the Impact of Gene Loss Following Whole-Genome Duplication

Whole-genome duplication (WGD) occurs broadly and repeatedly across the history of eukaryotes and is recognized as a prominent evolutionary force, especially in plants. Immediately following WGD, most genes are present in two copies as paralogs. Due to this redundancy, one copy of a paralog pair commonly undergoes pseudogenization and is eventually lost. When speciation occurs shortly after WGD; however, differential loss of paralogs may lead to spurious phylogenetic inference resulting from the inclusion of pseudoorthologs–paralogous genes mistakenly identified as orthologs because they are present in single copies within each sampled species. The influence and impact of including pseudoorthologs versus true orthologs as a result of gene extinction (or incomplete laboratory sampling) are only recently gaining empirical attention in the phylogenomics community. Moreover, few studies have yet to investigate this phenomenon in an explicit coalescent framework. Here, using mathematical models, numerous simulated data sets, and two newly assembled empirical data sets, we assess the effect of pseudoorthologs on species tree estimation under varying degrees of incomplete lineage sorting (ILS) and differential gene loss scenarios following WGD. When gene loss occurs along the terminal branches of the species tree, alignment-based (BPP) and gene-tree-based (ASTRAL, MP-EST, and STAR) coalescent methods are adversely affected as the degree of ILS increases. This can be greatly improved by sampling a sufficiently large number of genes. Under the same circumstances, however, concatenation methods consistently estimate incorrect species trees as the number of genes increases. Additionally, pseudoorthologs can greatly mislead species tree inference when gene loss occurs along the internal branches of the species tree. Here, both coalescent and concatenation methods yield inconsistent results. These results underscore the importance of understanding the influence of pseudoorthologs in the phylogenomics era. [Coalescent method; concatenation method; incomplete lineage sorting; pseudoorthologs; single-copy gene; whole-genome duplication.]

A Population-Genetic Lens into the Process of Gene Loss Following Whole-Genome Duplication

Whole-genome duplications (WGDs) have occurred in many eukaryotic lineages. However, the underlying evolutionary forces and molecular mechanisms responsible for the long-term retention of gene duplicates created by WGDs are not well understood. We employ a population-genomic approach to understand the selective forces acting on paralogs and investigate ongoing duplicate-gene loss in multiple species of Paramecium that share an ancient WGD. We show that mutations that abolish protein function are more likely to be segregating in retained WGD paralogs than in single-copy genes, most likely because of ongoing nonfunctionalization post-WGD. This relaxation of purifying selection occurs in only one WGD paralog, accompanied by the gradual fixation of nonsynonymous mutations and reduction in levels of expression, and occurs over a long period of evolutionary time, "marking" one locus for future loss. Concordantly, the fitness effects of new nonsynonymous mutations and frameshift-causing indels are significantly more deleterious in the highly expressed copy compared with their paralogs with lower expression. Our results provide a novel mechanistic model of gene duplicate loss following WGDs, wherein selection acts on the sum of functional activity of both duplicate genes, allowing the two to wander in expression and functional space, until one duplicate locus eventually degenerates enough in functional efficiency or expression that its contribution to total activity is too insignificant to be retained by purifying selection. Retention of duplicates by such mechanisms predicts long times to duplicate-gene loss, which should not be falsely attributed to retention due to gain/change in function.

Phylogeny of NF-YA trans-activation splicing isoforms in vertebrate evolution

NF-Y is a trimeric pioneer Transcription Factor (TF) whose target sequence -the CCAAT box- is present in ~25% of mammalian promoters. We reconstruct the phylogenetic history of the regulatory NF-YA subunit in vertebrates. We find that in addition to the remarkable conservation of the subunits-interaction and DNA-binding parts, the Transcriptional Activation Domain (TAD) is also conserved (>90% identity among bony vertebrates). We infer the phylogeny of the alternatively spliced exon-3 and partial splicing events of exon-7 -7N and 7C- revealing independent clade-specific losses of these regions. These isoforms shape the TAD. Absence of exon-3 in basal deuterostomes, cartilaginous fishes and hagfish, but not in lampreys, suggests that the "short" isoform is primordial, with emergence of exon-3 in chordates. Exon 7N was present in the vertebrate common ancestor, while 7C is a molecular innovation of teleost fishes. RNA-seq analysis in several species confirms expression of all these isoforms. We identify 3 blocks of amino acids in the TAD shared across deuterostomes, yet structural predictions and sequence analyses suggest an evolutionary drive for maintenance of an Intrinsically Disordered Region -IDR- within the TAD. Overall, these data help reconstruct the logic for alternative splicing of this essential eukaryotic TF.

Polyploidy as a Fundamental Phenomenon in Evolution, Development, Adaptation and Diseases

DNA replication during cell proliferation is 'vertical' copying, which reproduces an initial amount of genetic information. Polyploidy, which results from whole-genome duplication, is a fundamental complement to vertical copying. Both organismal and cell polyploidy can emerge via premature cell cycle exit or via cell-cell fusion, the latter giving rise to polyploid hybrid organisms and epigenetic hybrids of somatic cells. Polyploidy-related increase in biological plasticity, adaptation, and stress resistance manifests in evolution, development, regeneration, aging, oncogenesis, and cardiovascular diseases. Despite the prevalence in nature and importance for medicine, agri- and aquaculture, biological processes and epigenetic mechanisms underlying these fundamental features largely remain unknown. The evolutionarily conserved features of polyploidy include activation of transcription, response to stress, DNA damage and hypoxia, and induction of programs of morphogenesis, unicellularity, and longevity, suggesting that these common features confer adaptive plasticity, viability, and stress resistance to polyploid cells and organisms. By increasing cell viability, polyploidization can provide survival under stressful conditions where diploid cells cannot survive. However, in somatic cells it occurs at the expense of specific function, thus promoting developmental programming of adult cardiovascular diseases and increasing the risk of cancer. Notably, genes arising via evolutionary polyploidization are heavily involved in cancer and other diseases. Ploidy-related changes of gene expression presumably originate from chromatin modifications and the derepression of bivalent genes. The provided evidence elucidates the role of polyploidy in evolution, development, aging, and carcinogenesis, and may contribute to the development of new strategies for promoting regeneration and preventing cardiovascular diseases and cancer.

Tissue-Targeted CRISPR-Cas9-Mediated Genome Editing of Multiple Homeologs in F0-Generation Xenopus laevis Embryos

Xenopus laevis frogs are a powerful developmental model that enables studies combining classical embryology and molecular manipulation. Because of the large embryo size, ease of microinjection, and ability to target tissues through established fate maps, X. laevis has become the predominant amphibian research model. Given that their allotetraploid genome has complicated the generation of gene knockouts, strategies need to be established for efficient mutagenesis of multiple homeologs to evaluate gene function. Here we describe a protocol to use CRISPR-Cas9-mediated genome editing to target either single alleles or multiple alloalleles in F0 X. laevis embryos. A single-guide RNA (sgRNA) is designed to target a specific DNA sequence encoding a critical protein domain. To mutagenize a gene with two alloalleles, the sgRNA is designed against a sequence that is common to both homeologs. This sgRNA, along with the Cas9 protein, is microinjected into the zygote to disrupt the genomic sequences in the whole embryo or into a specific blastomere for tissue-targeted effects. Error-prone repair of CRISPR-Cas9-generated DNA double-strand breaks leads to insertions and deletions creating mosaic gene lesions within the embryos. The genomic DNA isolated from each mosaic F0 embryo is sequenced, and software is applied to assess the nature of the mutations generated and degree of mosaicism. This protocol enables the knockout of genes within the whole embryo or in specific tissues in F0 X. laevis embryos to facilitate the evaluation of resulting phenotypes.

Molecular evolution of hatching enzymes and their paralogous genes in vertebrates

BACKGROUND

Hatching is identified as one of the most important events in the reproduction of oviparous vertebrates. The genes for hatching enzymes, which are vital in the hatching process, are conserved among vertebrates. However, especially in teleost, it is difficult to trace their molecular evolution in detail due to the presence of other C6astacins, which are the subfamily to which the genes for hatching enzymes belong and are highly diverged. In particular, the hatching enzyme genes are diversified with frequent genome translocations due to retrocopy.

RESULTS

In this study, we took advantage of the rapid expansion of whole-genome data in recent years to examine the molecular evolutionary process of these genes in vertebrates. The phylogenetic analysis and the genomic synteny analysis revealed C6astacin genes other than the hatching enzyme genes, which was previously considered to be retained only in teleosts, was also retained in the genomes of basal ray-finned fishes, coelacanths, and cartilaginous fishes. These results suggest that the common ancestor of these genes can be traced back to at least the common ancestor of the Gnathostomata. Moreover, we also found that many of the C6astacin genes underwent multiple gene duplications during vertebrate evolution, and the results of gene expression analysis in frogs implied that genes derived from hatching enzyme genes underwent neo-functionalization.

CONCLUSIONS

In this study, we describe in detail the molecular evolution of the C6astacin gene in vertebrates, which has not been summarized previously. The results revealed the presence of the previously unknown C6astacin gene in the basal-lineage of jawed vertebrates and large-scale gene duplication of hatching enzyme genes in amphibians. The comprehensive investigation reported in this study will be an important basis for studying the molecular evolution of the vertebrate C6astacin genes, hatching enzyme, and its paralogous genes and for identifying these genes without the need for gene expression and functional analysis.

Studying phenotypic variation and DNA methylation across development, ecology and evolution in the clonal marbled crayfish: a paradigm for investigating epigenotype-phenotype relationships in macro-invertebrates

Animals can produce different phenotypes from the same genome during development, environmental adaptation and evolution, which is mediated by epigenetic mechanisms including DNA methylation. The obligatory parthenogenetic marbled crayfish, Procambarus virginalis, whose genome and methylome are fully established, proved very suitable to study this issue in detail. Comparison between developmental stages and DNA methylation revealed low expression of Dnmt methylation and Tet demethylation enzymes from the spawned oocyte to the 256 cell embryo and considerably increased expression thereafter. The global 5-methylcytosine level was 2.78% at mid-embryonic development and decreased slightly to 2.41% in 2-year-old adults. Genetically identical clutch-mates raised in the same uniform laboratory setting showed broad variation in morphological, behavioural and life history traits and differences in DNA methylation. The invasion of diverse habitats in tropical to cold-temperate biomes in the last 20 years by the marbled crayfish was associated with the expression of significantly different phenotypic traits and DNA methylation patterns, despite extremely low genetic variation on the whole genome scale, suggesting the establishment of epigenetic ecotypes. The evolution of marbled crayfish from its parent species Procambarus fallax by autotriploidy a few decades ago was accompanied by a significant increase in body size, fertility and life span, a 20% reduction of global DNA methylation and alteration of methylation in hundreds of genes, suggesting that epigenetic mechanisms were involved in speciation and fitness enhancement. The combined analysis of phenotypic traits and DNA methylation across multiple biological contexts in the laboratory and field in marbled crayfish may serve as a blueprint for uncovering the role of epigenetic mechanisms in shaping of phenotypes in macro-invertebrates.

Identification and functional characterization of interleukin-12 receptor beta 1 and 2 in grass carp (Ctenopharyngodon idella)

Interleukin 12 (IL-12) binds its receptor complex of IL-12 receptor beta 1 (IL-12Rβ1) and IL-12Rβ2 to transduce cellular signaling in mammals. In teleosts, the function of Il-12 is drawing increasing attention, but molecular and functional features of Il-12 receptors remain obscure. Especially, the existence of multiple Il-12 isoforms in some fish species elicits the requirement to clarify their receptors. In this study, we isolated three cDNA sequences as Il-12 receptor candidates from grass carp, entitled as grass carp Il-12rβ1 (gcIl-12rβ1), gcIl-12rβ2a and gcIl-12rβ2b. In silico analysis showed that gcIl-12rβ1 and gcIl-12rβ2a shared the conserved gene locus and similar structure characteristics with their orthologues of zebrafish, frog, chicken, mouse and human, respectively. However, the Il-12rβ2b of grass carp and zebrafish was similar to IL-27Ra in non-fish species. Further locally installed BLAST and gene synteny analysis uncovered three gcIl-12 receptors being single copied genes. Tissue distribution assay revealed that gcil12rβ1 and gcil12rβ2a transcripts were predominantly expressed in head kidney, differing from the even distribution of gcil12rβ2b transcripts in all detected tissues. Subsequently, the binding ability and antagonistic effects of recombinant extracellular region of gcIl-12rβ1 with recombinant grass carp Il-12 (rgcIl-12) isoforms were explored, providing functional evidence of the newly cloned gcIl-12rβ1 being genuine orthologues of mammalian IL-12Rβ1. Moreover, our data showed that gcIl-12rβ1 and gcIl-12rβ2a but not gcIl-12rβ1 and gcIl-12rβ2b mediated the effects of rgcIl-12 isoforms on ifn-γ promoter activity, thereby revealing Il-12 receptor signaling in fish. These results identified grass carp Il-12 receptors, thereby advancing our understanding of Il-12 isoform signaling in fish.

Genome-Wide Reconstruction of Rediploidization Following Autopolyploidization across One Hundred Million Years of Salmonid Evolution

The long-term evolutionary impacts of whole-genome duplication (WGD) are strongly influenced by the ensuing rediploidization process. Following autopolyploidization, rediploidization involves a transition from tetraploid to diploid meiotic pairing, allowing duplicated genes (ohnologs) to diverge genetically and functionally. Our understanding of autopolyploid rediploidization has been informed by a WGD event ancestral to salmonid fishes, where large genomic regions are characterized by temporally delayed rediploidization, allowing lineage-specific ohnolog sequence divergence in the major salmonid clades. Here, we investigate the long-term outcomes of autopolyploid rediploidization at genome-wide resolution, exploiting a recent "explosion" of salmonid genome assemblies, including a new genome sequence for the huchen (Hucho hucho). We developed a genome alignment approach to capture duplicated regions across multiple species, allowing us to create 121,864 phylogenetic trees describing genome-wide ohnolog divergence across salmonid evolution. Using molecular clock analysis, we show that 61% of the ancestral salmonid genome experienced an initial "wave" of rediploidization in the late Cretaceous (85-106 Ma). This was followed by a period of relative genomic stasis lasting 17-39 My, where much of the genome remained tetraploid. A second rediploidization wave began in the early Eocene and proceeded alongside species diversification, generating predictable patterns of lineage-specific ohnolog divergence, scaling in complexity with the number of speciation events. Using gene set enrichment, gene expression, and codon-based selection analyses, we provide insights into potential functional outcomes of delayed rediploidization. This study enhances our understanding of delayed autopolyploid rediploidization and has broad implications for future studies of WGD events.

On the relationship between extant innate immune receptors and the evolutionary origins of jawed vertebrate adaptive immunity

For over half a century, deciphering the origins of the genomic loci that form the jawed vertebrate adaptive immune response has been a major topic in comparative immunogenetics. Vertebrate adaptive immunity relies on an extensive and highly diverse repertoire of tandem arrays of variable (V), diversity (D), and joining (J) gene segments that recombine to produce different immunoglobulin (Ig) and T cell receptor (TCR) genes. The current consensus is that a recombination-activating gene (RAG)-like transposon invaded an exon of an ancient innate immune VJ-bearing receptor, giving rise to the extant diversity of Ig and TCR loci across jawed vertebrates. However, a model for the evolutionary relationships between extant non-recombining innate immune receptors and the V(D)J receptors of the jawed vertebrate adaptive immune system has only recently begun to come into focus. In this review, we provide an overview of non-recombining VJ genes, including CD8β, CD79b, natural cytotoxicity receptor 3 (NCR3/NKp30), putative remnants of an antigen receptor precursor (PRARPs), and the multigene family of signal-regulatory proteins (SIRPs), that play a wide range of roles in immune function. We then focus in detail on the VJ-containing novel immune-type receptors (NITRs) from ray-finned fishes, as recent work has indicated that these genes are at least 50 million years older than originally thought. We conclude by providing a conceptual model of the evolutionary origins and phylogenetic distribution of known VJ-containing innate immune receptors, highlighting opportunities for future comparative research that are empowered by this emerging evolutionary perspective.

Rhythmic Clock Gene Expression in Atlantic Salmon Parr Brain

To better understand the complexity of clock genes in salmonids, a taxon with an additional whole genome duplication, an analysis was performed to identify and classify gene family members (clock, arntl, period, cryptochrome, nr1d, ror, and csnk1). The majority of clock genes, in zebrafish and Northern pike, appeared to be duplicated. In comparison to the 29 clock genes described in zebrafish, 48 clock genes were discovered in salmonid species. There was also evidence of species-specific reciprocal gene losses conserved to the Oncorhynchus sister clade. From the six period genes identified three were highly significantly rhythmic, and circadian in their expression patterns (per1a.1, per1a.2, per1b) and two was significantly rhythmically expressed (per2a, per2b). The transcriptomic study of juvenile Atlantic salmon (parr) brain tissues confirmed gene identification and revealed that there were 2,864 rhythmically expressed genes (p < 0.001), including 1,215 genes with a circadian expression pattern, of which 11 were clock genes. The majority of circadian expressed genes peaked 2 h before and after daylight. These findings provide a foundation for further research into the function of clock genes circadian rhythmicity and the role of an enriched number of clock genes relating to seasonal driven life history in salmonids.

ORTHOSCOPE*: a phylogenetic pipeline to infer gene histories from genome-wide data

Comparative genome-scale analyses of protein-coding gene sequences are employed to examine evidence for whole-genome duplication and horizontal gene transfer. For this purpose, an orthogroup should be delineated to infer evolutionary history regarding each gene, and results of all orthogroup analyses need to be integrated to infer a genome-scale history. An orthogroup is a set of genes descended from a single gene in the last common ancestor of all species under consideration. However, such analyses confront several problems: (1) analytical pipelines to infer all gene histories with methods comparing species and gene trees are not fully developed, and (2) without detailed analyses within orthogroups, evolutionary events of paralogous genes in the same orthogroup cannot be distinguished for genome-wide integration of results derived from multiple orthogroup analyses. Here I present an analytical pipeline, ORTHOSCOPE* (star), to infer evolutionary histories of animal/plant genes from genome-scale data. ORTHOSCOPE* estimates a tree for a specified gene, detects speciation/gene duplication events that occurred at nodes belonging to only one lineage leading to a species of interest, and then integrates results derived from gene trees estimated for all query genes in genome-wide data. Thus, ORTHOSCOPE* can be used to detect species nodes just after whole genome duplications as a first step of comparative genomic analyses. Moreover, by examining the presence or absence of genes belonging to species lineages with dense taxon sampling available from the ORTHOSCOPE web version, ORTHOSCOPE* can detect genes lost in specific lineages and horizontal gene transfers. This pipeline is available at https://github.com/jun-inoue/ORTHOSCOPE_STAR.

Deciphering reticulate evolution of the largest group of polyploid vertebrates, the subfamily cyprininae (Teleostei: Cypriniformes)

Despite the rarity of polyploidy in animals, some groups with polyploid species exhibit complicated and interesting patterns of reticulate evolution. Here we focus on fishes in the subfamily Cyprininae, the largest polyploid group of vertebrates. The large number of polyploid taxa poses significant challenges for phylogenetic and evolutionary studies on this subfamily. In this study, we cloned and sequenced three single-copy nuclear loci to investigate the evolution of polyploidy in the Cyprininae. Topologies of nuclear gene trees were compared with a newly reconstructed mitochondrial tree. The data provided herein corroborate the hybrid origins of the tribes Torini, Cyprinini, Spinibarbini, Barbini, and also Probarbini. Based on results from this study and previous studies, we hypothesize that at least 13 independent polyploidization events have occurred during the evolution of the Cyprininae. We offer hypotheses on the origin of each polyploid group and show that a diploid group or the diploid ancestor of a polyploid group might have served as progenitor of one or two other polyploid groups. The evolutionary history of Cyprinine (since its first divergence) can be divided into three stages: the "Diploid stage" (69.2-43.4 Ma or million years ago), the "Tetraploidization stage" (43.4-18.9 Ma), and the "Hexaploidization stage" (18.9 Ma to present). The second stage is when all tetraploidization events happened, while the last stage is when all hexaploidization events and most genus- or species-specific polyploidization events occurred. The post-polyploidization dynamics of polyploid groups are complicated and warrant more genomic level studies. We showed that the subfamily Cyprininae may represent a more complicated polyploid system than most, if not all, other vertebrates and some plants, if one or more of the following factors are considered: number of polyploid species, number of different ploidy levels, and number and type of independent polyploidization events.

Holosteans contextualize the role of the teleost genome duplication in promoting the rise of evolutionary novelties in the ray-finned fish innate immune system

Over 99% of ray-finned fishes (Actinopterygii) are teleosts, a clade that comprises half of all living vertebrate species that have diversified across virtually all fresh and saltwater ecosystems. This ecological breadth raises the question of how the immunogenetic diversity required to persist under heterogeneous pathogen pressures evolved. The teleost genome duplication (TGD) has been hypothesized as the evolutionary event that provided the substrate for rapid genomic evolution and innovation. However, studies of putative teleost-specific innate immune receptors have been largely limited to comparisons either among teleosts or between teleosts and distantly related vertebrate clades such as tetrapods. Here we describe and characterize the receptor diversity of two clustered innate immune gene families in the teleost sister lineage: Holostei (bowfin and gars). Using genomic and transcriptomic data, we provide a detailed investigation of the phylogenetic history and conserved synteny of gene clusters encoding diverse immunoglobulin domain-containing proteins (DICPs) and novel immune-type receptors (NITRs). These data demonstrate an ancient linkage of DICPs to the major histocompatibility complex (MHC) and reveal an evolutionary origin of NITR variable-joining (VJ) exons that predate the TGD by at least 50 million years. Further characterizing the receptor diversity of Holostean DICPs and NITRs illuminates a sequence diversity that rivals the diversity of these innate immune receptor families in many teleosts. Taken together, our findings provide important historical context for the evolution of these gene families that challenge prevailing expectations concerning the consequences of the TGD during actinopterygiian evolution.

Mirage: estimation of ancestral gene-copy numbers by considering different evolutionary patterns among gene families

Motivation

Reconstruction of gene copy number evolution is an essential approach for understanding how complex biological systems have been organized. Although various models have been proposed for gene copy number evolution, existing evolutionary models have not appropriately addressed the fact that different gene families can have very different gene gain/loss rates.

Results

In this study, we developed Mirage (MIxtuRe model for Ancestral Genome Estimation), which allows different gene families to have flexible gene gain/loss rates. Mirage can use three models for formulating heterogeneous evolution among gene families: the discretized Γ model, probability distribution-free model and pattern mixture (PM) model. Simulation analysis showed that Mirage can accurately estimate heterogeneous gene gain/loss rates and reconstruct gene-content evolutionary history. Application to empirical datasets demonstrated that the PM model fits genome data from various taxonomic groups better than the other heterogeneous models. Using Mirage, we revealed that metabolic function-related gene families displayed frequent gene gains and losses in all taxa investigated.

Availability and implementation

The source code of Mirage is freely available at https://github.com/fukunagatsu/Mirage.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Fossilized cell structures identify an ancient origin for the teleost whole-genome duplication

Teleost fishes comprise one-half of all vertebrate species and possess a duplicated genome. This whole-genome duplication (WGD) occurred on the teleost stem lineage in an ancient common ancestor of all living teleosts and is hypothesized as a trigger of their exceptional evolutionary radiation. Genomic and phylogenetic data indicate that WGD occurred in the Mesozoic after the divergence of teleosts from their closest living relatives but before the origin of the extant teleost groups. However, these approaches cannot pinpoint WGD among the many extinct groups that populate this 50- to 100-million-y lineage, preventing tests of the evolutionary effects of WGD. We infer patterns of genome size evolution in fossil stem-group teleosts using high-resolution synchrotron X-ray tomography to measure the bone cell volumes, which correlate with genome size in living species. Our findings indicate that WGD occurred very early on the teleost stem lineage and that all extinct stem-group teleosts known so far possessed duplicated genomes. WGD therefore predates both the origin of proposed key innovations of the teleost skeleton and the onset of substantial morphological diversification in the clade. Moreover, the early occurrence of WGD allowed considerable time for postduplication reorganization prior to the origin of the teleost crown group. This suggests at most an indirect link between WGD and evolutionary success, with broad implications for the relationship between genomic architecture and large-scale evolutionary patterns in the vertebrate Tree of Life.

Structural Identification of the Pacemaker Cells and Expression of Hyperpolarization-Activated Cyclic Nucleotide-Gated (HCN) Channels in the Heart of the Wild Atlantic Cod, Gadus morhua (Linnaeus, 1758)

Hyperpolarization-activated cyclic nucleotide-gated (HCN) channels are proteins that contain highly conserved functional domains and sequence motifs that are correlated with their unique biophysical activities, to regulate cardiac pacemaker activity and synaptic transmission. These pacemaker proteins have been studied in mammalian species, but little is known now about their heart distribution in lower vertebrates and c-AMP modulation. Here, we characterized the pacemaker system in the heart of the wild Atlantic cod (Gadus morhua), with respect to primary pacemaker molecular markers. Special focus is given to the structural, ultrastructural and molecular characterization of the pacemaker domain, through the expression of HCN channel genes and the immunohistochemistry of HCN isoforms, including the location of intracardiac neurons that are adjacent to the sinoatrial region of the heart. Similarly to zebrafish and mammals, these neurons are immunoreactive to ChAT, VAChT and nNOS. It has been shown that cardiac pacemaking can be modulated by sympathetic and parasympathetic pathways, and the existence of intracardiac neurons projecting back to the central nervous system provide a plausible link between them.

The digestive tract as an essential organ for water acquisition in marine teleosts: lessons from euryhaline eels

Adaptation to a hypertonic marine environment is one of the major topics in animal physiology research. Marine teleosts lose water osmotically from the gills and compensate for this loss by drinking surrounding seawater and absorbing water from the intestine. This situation is in contrast to that in mammals, which experience a net osmotic loss of water after drinking seawater. Water absorption in fishes is made possible by (1) removal of monovalent ions (desalinization) by the esophagus, (2) removal of divalent ions as carbonate (Mg/CaCO₃) precipitates promoted by HCO₃^- secretion, and (3) facilitation of NaCl and water absorption from diluted seawater by the intestine using a suite of unique transporters. As a result, 70-85% of ingested seawater is absorbed during its passage through the digestive tract. Thus, the digestive tract is an essential organ for marine teleost survival in the hypertonic seawater environment. The eel is a species that has been frequently used for osmoregulation research in laboratories worldwide. The eel possesses many advantages as an experimental animal for osmoregulation studies, one of which is its outstanding euryhalinity, which enables researchers to examine changes in the structure and function of the digestive tract after direct transfer from freshwater to seawater. In recent years, the molecular mechanisms of ion and water transport across epithelial cells (the transcellular route) and through tight junctions (the paracellular route) have been elucidated for the esophagus and intestine. Thanks to the rapid progress in analytical methods for genome databases on teleosts, including the eel, the molecular identities of transporters, channels, pumps and junctional proteins have been clarified at the isoform level. As 10 y have passed since the previous reviews on this subject, it seems relevant and timely to summarize recent progress in research on the molecular mechanisms of water and ion transport in the digestive tract in eels and to compare the mechanisms with those of other teleosts and mammals from comparative and evolutionary viewpoints. We also propose future directions for this research field to achieve integrative understanding of the role of the digestive tract in adaptation to seawater with regard to pathways/mechanisms including the paracellular route, divalent ion absorption, metabolon formation and cellular trafficking of transporters. Notably, some of these have already attracted practical attention in laboratories.

Functional Domains and Evolutionary History of the PMEL and GPNMB Family Proteins

The ancient paralogs premelanosome protein (PMEL) and glycoprotein nonmetastatic melanoma protein B (GPNMB) have independently emerged as intriguing disease loci in recent years. Both proteins possess common functional domains and variants that cause a shared spectrum of overlapping phenotypes and disease associations: melanin-based pigmentation, cancer, neurodegenerative disease and glaucoma. Surprisingly, these proteins have yet to be shown to physically or genetically interact within the same cellular pathway. This juxtaposition inspired us to compare and contrast this family across a breadth of species to better understand the divergent evolutionary trajectories of two related, but distinct, genes. In this study, we investigated the evolutionary history of PMEL and GPNMB in clade-representative species and identified TMEM130 as the most ancient paralog of the family. By curating the functional domains in each paralog, we identified many commonalities dating back to the emergence of the gene family in basal metazoans. PMEL and GPNMB have gained functional domains since their divergence from TMEM130, including the core amyloid fragment (CAF) that is critical for the amyloid potential of PMEL. Additionally, the PMEL gene has acquired the enigmatic repeat domain (RPT), composed of a variable number of imperfect tandem repeats; this domain acts in an accessory role to control amyloid formation. Our analyses revealed the vast variability in sequence, length and repeat number in homologous RPT domains between craniates, even within the same taxonomic class. We hope that these analyses inspire further investigation into a gene family that is remarkable from the evolutionary, pathological and cell biology perspectives.

The divergence of alternative splicing between ohnologs in teleost fishes

BACKGROUND

Gene duplication and alternative splicing (AS) are two distinct mechanisms generating new materials for genetic innovations. The evolutionary link between gene duplication and AS is still controversial, due to utilizing duplicates from inconsistent ages of duplication events in earlier studies. With the aid of RNA-seq data, we explored evolutionary scenario of AS divergence between duplicates with ohnologs that resulted from the teleost genome duplication event in zebrafish, medaka, and stickleback.

RESULTS

Ohnologs in zebrafish have fewer AS forms compared to their singleton orthologs, supporting the function-sharing model of AS divergence between duplicates. Ohnologs in stickleback have more AS forms compared to their singleton orthologs, which supports the accelerated model of AS divergence between duplicates. The evolution of AS in ohnologs in medaka supports a combined scenario of the function-sharing and the accelerated model of AS divergence between duplicates. We also found a small number of ohnolog pairs in each of the three teleosts showed significantly asymmetric AS divergence. For example, the well-known ovary-factor gene cyp19a1a has no AS form but its ohnolog cyp19a1b has multiple AS forms in medaka, suggesting that functional divergence between duplicates might have result from AS divergence.

CONCLUSIONS

We found that a combined scenario of function-sharing and accelerated models for AS evolution in ohnologs in teleosts and rule out the independent model that assumes a lack of correlation between gene duplication and AS. Our study thus provided insights into the link between gene duplication and AS in general and ohnolog divergence in teleosts from AS perspective in particular.

Evolution after Whole-Genome Duplication: Teleost MicroRNAs

MicroRNAs (miRNAs) are important gene expression regulators implicated in many biological processes, but we lack a global understanding of how miRNA genes evolve and contribute to developmental canalization and phenotypic diversification. Whole-genome duplication events likely provide a substrate for species divergence and phenotypic change by increasing gene numbers and relaxing evolutionary pressures. To understand the consequences of genome duplication on miRNA evolution, we studied miRNA genes following the teleost genome duplication (TGD). Analysis of miRNA genes in four teleosts and in spotted gar, whose lineage diverged before the TGD, revealed that miRNA genes were retained in ohnologous pairs more frequently than protein-coding genes, and that gene losses occurred rapidly after the TGD. Genomic context influenced retention rates, with clustered miRNA genes retained more often than nonclustered miRNA genes and intergenic miRNA genes retained more frequently than intragenic miRNA genes, which often shared the evolutionary fate of their protein-coding host. Expression analyses revealed both conserved and divergent expression patterns across species in line with miRNA functions in phenotypic canalization and diversification, respectively. Finally, major strands of miRNA genes experienced stronger purifying selection, especially in their seeds and 3'-complementary regions, compared with minor strands, which nonetheless also displayed evolutionary features compatible with constrained function. This study provides the first genome-wide, multispecies analysis of the mechanisms influencing metazoan miRNA evolution after whole-genome duplication.

A New Zebrafish Model for Pseudoxanthoma Elasticum

Calcification of various tissues is a significant health issue associated with aging, cancer and autoimmune diseases. There are both environmental and genetic factors behind this phenomenon and understanding them is essential for the development of efficient therapeutic approaches. Pseudoxanthoma elasticum (PXE) is a rare genetic disease, a prototype for calcification disorders, resulting from the dysfunction of ABCC6, a transport protein found in the membranes of cells. It is identified by excess calcification in a variety of tissues (e.g., eyes, skin, arteries) and currently it has no cure, known treatments target the symptoms only. Preclinical studies of PXE have been successful in mice, proving the usefulness of animal models for the study of the disease. Here, we present a new zebrafish (Danio rerio) model for PXE. By resolving some ambiguous assemblies in the zebrafish genome, we show that there are two functional and one non-functional paralogs for ABCC6 in zebrafish (abcc6a, abcc6b.1, and abcc6b.2, respectively). We created single and double mutants for the functional paralogs and characterized their calcification defects with a combination of techniques. Zebrafish deficient in abcc6a show defects in their vertebral calcification and also display ectopic calcification foci in their soft tissues. Our results also suggest that the impairment of abcc6b.1 does not affect this biological process.

Genome-wide identification and characterization of olfactory receptor genes in common carp (Cyprinus carpio)

The environment contains a large extent of chemical information, which could be detected as olfactory sense. Olfactory in vertebrates plays important roles on many aspects during life time, including localizing prey or food, avoiding predators, mating behavior and social communication. Considering the essential role of olfactory receptors in the specific recognition of diverse stimuli, understanding the evolutionary dynamics of olfactory receptors in teleost means a lot, especially in the allotetraploid common carp, who has undergone the fourth whole-genome duplication event. Here, we identified the whole set of olfactory receptor genes in representative teleosts and found a significant contraction in common carp when compared with other teleosts. Odorant receptor genes (OR) occupy the most among four groups of olfactory receptors, including 33 functional genes and 16 pseudogenes. Furthermore, 6 trace amine-associated receptor (TAAR) genes (including 1 pseudogene), 7 odorant-related-A receptor genes, and 10 olfactory C family receptor genes (including 3 pseudogenes) were identified in common carp. Phylogenetic and motif analysis were performed to illustrate the phylogenetic relationship and structural conservation of teleost olfactory receptors. Selection pressure analysis suggested that olfactory receptor groups in common carp were all under relaxed purifying-selection. Additionally, gene expression divergences for olfactory receptor genes were investigated during embryonic development stages of common carp. We aim to determine the abundance of common carp olfactory receptor genes, explore the evolutionary fate and expression dynamics, and provide some genomic clues for the evolution of polyploid olfactory after whole-genome duplication and for future studies of teleost olfactory.

A reappraisal of the phylogenetic placement of the Aquilegia whole-genome duplication

The accurate placement of an ancient whole-genome duplication (WGD) in relation to the lineage divergence is important. Here, we re-investigated the Aquilegia coerulea WGD and found it is more likely lineage-specific rather than shared by all eudicots.

Draft Genome of the Common Snapping Turtle, Chelydra serpentina, a Model for Phenotypic Plasticity in Reptiles

Turtles are iconic reptiles that inhabit a range of ecosystems from oceans to deserts and climates from the tropics to northern temperate regions. Yet, we have little understanding of the genetic adaptations that allow turtles to survive and reproduce in such diverse environments. Common snapping turtles, Chelydra serpentina, are an ideal model species for studying adaptation to climate because they are widely distributed from tropical to northern temperate zones in North America. They are also easy to maintain and breed in captivity and produce large clutch sizes, which makes them amenable to quantitative genetic and molecular genetic studies of traits like temperature-dependent sex determination. We therefore established a captive breeding colony and sequenced DNA from one female using both short and long reads. After trimming and filtering, we had 209.51Gb of Illumina reads, 25.72Gb of PacBio reads, and 21.72 Gb of Nanopore reads. The assembled genome was 2.258 Gb in size and had 13,224 scaffolds with an N50 of 5.59Mb. The longest scaffold was 27.24Mb. BUSCO analysis revealed 97.4% of core vertebrate genes in the genome. We identified 3.27 million SNPs in the reference turtle, which indicates a relatively high level of individual heterozygosity. We assembled the transcriptome using RNA-Seq data and used gene prediction software to produce 22,812 models of protein coding genes. The quality and contiguity of the snapping turtle genome is similar to or better than most published reptile genomes. The genome and genetic variants identified here provide a foundation for future studies of adaptation to climate.

Exon probe sets and bioinformatics pipelines for all levels of fish phylogenomics

Exon markers have a long history of use in phylogenetics of ray-finned fishes, the most diverse clade of vertebrates with more than 35,000 species. As the number of published genomes increases, it has become easier to test exons and other genetic markers for signals of ancient duplication events and filter out paralogues that can mislead phylogenetic analysis. We present seven new probe sets for current target-capture phylogenomic protocols that capture 1,104 exons explicitly filtered for paralogues using gene trees. These seven probe sets span the diversity of teleost fishes, including four sets that target five hyperdiverse percomorph clades which together comprise ca. 17,000 species (Carangaria, Ovalentaria, Eupercaria, and Syngnatharia + Pelagiaria combined). We additionally included probes to capture legacy nuclear exons and mitochondrial markers that have been commonly used in fish phylogenetics (despite some exons being flagged for paralogues) to facilitate integration of old and new molecular phylogenetic matrices. We tested these probes experimentally for 56 fish species (eight species per probe set) and merged new exon-capture sequence data into an existing data matrix of 1,104 exons and 300 ray-finned fish species. We provide an optimized bioinformatics pipeline to assemble exon capture data from raw reads to alignments for downstream analysis. We show that legacy loci with known paralogues are at risk of assembling duplicated sequences with target-capture, but we also assembled many useful orthologous sequences that can be integrated with many PCR-generated matrices. These probe sets are a valuable resource for advancing fish phylogenomics because targeted exons can easily be extracted from increasingly available whole genome and transcriptome data sets, and also may be integrated with existing PCR-based exon and mitochondrial data.

Global assessment of organ specific basal gene expression over a diurnal cycle with analyses of gene copies exhibiting cyclic expression patterns

Background

Studying functional divergences between paralogs that originated from genome duplication is a significant topic in investigating molecular evolution. Genes that exhibit basal level cyclic expression patterns including circadian and light responsive genes are important physiological regulators. Temporal shifts in basal gene expression patterns are important factors to be considered when studying genetic functions. However, adequate efforts have not been applied to studying basal gene expression variation on a global scale to establish transcriptional activity baselines for each organ. Furthermore, the investigation of cyclic expression pattern comparisons between genome duplication created paralogs, and potential functional divergence between them has been neglected. To address these questions, we utilized a teleost fish species, Xiphophorus maculatus, and profiled gene expression within 9 organs at 3-h intervals throughout a 24-h diurnal period.

Results

Our results showed 1.3–21.9% of genes in different organs exhibited cyclic expression patterns, with eye showing the highest fraction of cycling genes while gonads yielded the lowest. A majority of the duplicated gene pairs exhibited divergences in their basal level expression patterns wherein only one paralog exhibited an oscillating expression pattern, or both paralogs exhibit oscillating expression patterns, but each gene duplicate showed a different peak expression time, and/or in different organs.

Conclusions

These observations suggest cyclic genes experienced significant sub-, neo-, or non-functionalization following the teleost genome duplication event. In addition, we developed a customized, web-accessible, gene expression browser to facilitate data mining and data visualization for the scientific community.

Directional divergence of Ep300 duplicates in teleosts and its implications

Background

EP300 is a conserved protein in vertebrates, which serves as a key mediator of cellular homeostasis. Mutations and dysregulation of EP300 give rise to severe human developmental disorders and malignancy. Danio rerio is a promising model organism to study EP300 related diseases and drugs; however, the effect of EP300 duplicates derived from teleost-specific whole genome duplication should not just be neglected.

Results

In this study, we obtained EP300 protein sequences of representative teleosts, mammals and sauropsids, with which we inferred a highly supported maximum likelihood tree. We observed that Ep300 duplicates (Ep300a and Ep300b) were widely retained in teleosts and universally expressed in a variety of tissues. Consensus sequences of Ep300a and Ep300b had exactly the same distribution of conserved domains, suggesting that their functions should still be largely overlapped. We analyzed the molecular evolution of Ep300 duplicates in teleosts, using branch-site models, clade models and site models. The results showed that both duplicates were subject to strong positive selection; however, for an extant species, generally at most one copy was under positive selection. At the clade level, there were evident positive correlations between evolutionary rates, the number of positively selected sites and gene expression levels. In Ostariophysi, Ep300a were under stronger positive selection than Ep300b; in Neoteleostei, another species-rich teleost clade, the contrary was the case. We also modeled 3D structures of zf-TAZ domain and its flanking regions of Ep300a and Ep300b of D. rerio and Oryzias latipes and found that in either species the faster evolving copy had more short helixes.

Conclusions

Collectively, the two copies of Ep300 have undoubtedly experienced directional divergence in main teleost clades. The divergence of EP300 between teleosts and mammals should be greater than the divergence between different teleost clades. Further studies are needed to clarify to what extent the EP300 involved regulatory network has diverged between teleosts and mammals, which would also help explain the huge success of teleosts.

Synteny-Guided Resolution of Gene Trees Clarifies the Functional Impact of Whole-Genome Duplications

Whole-genome duplications (WGDs) have major impacts on the evolution of species, as they produce new gene copies contributing substantially to adaptation, isolation, phenotypic robustness, and evolvability. They result in large, complex gene families with recurrent gene losses in descendant species that sequence-based phylogenetic methods fail to reconstruct accurately. As a result, orthologs and paralogs are difficult to identify reliably in WGD-descended species, which hinders the exploration of functional consequences of WGDs. Here, we present Synteny-guided CORrection of Paralogies and Orthologies (SCORPiOs), a novel method to reconstruct gene phylogenies in the context of a known WGD event. WGDs generate large duplicated syntenic regions, which SCORPiOs systematically leverages as a complement to sequence evolution to infer the evolutionary history of genes. We applied SCORPiOs to the 320-My-old WGD at the origin of teleost fish. We find that almost one in four teleost gene phylogenies in the Ensembl database (3,394) are inconsistent with their syntenic contexts. For 70% of these gene families (2,387), we were able to propose an improved phylogenetic tree consistent with both the molecular substitution distances and the local syntenic information. We show that these synteny-guided phylogenies are more congruent with the species tree, with sequence evolution and with expected expression conservation patterns than those produced by state-of-the-art methods. Finally, we show that synteny-guided gene trees emphasize contributions of WGD paralogs to evolutionary innovations in the teleost clade.