Fish lateral line innovation: insights into the evolutionary genomic dynamics of a unique mechanosensory organ

Fish Evo-Devo: Moving Toward Species-Specific and Knowledge-Based Interactome

A knowledge-based interactome maps interactions among proteins and molecules within a cell using experimental data, computational predictions, and literature mining. These interactomes are vital for understanding cellular functions, pathways, and the evolutionary conservation of protein interactions. They reveal how interactions regulate growth, differentiation, and development. Transitioning to functionally validated interactomes is crucial in evolutionary developmental biology (Evo-Devo), especially for non-model species, to uncover unique regulatory networks, evolutionary novelties, and reliable gene interaction models. This enhances our understanding of complex trait evolution across species. The European Evo-Devo 2024 conference in Helsinki hosted the first fish-specific Evo-Devo symposium, highlighting the growing interest in fish models. Advances in genome annotation, genome editing, imaging, and molecular screening are expanding fish Evo-Devo research. High-throughput molecular data have enabled the deduction of gene regulatory networks. The next steps involve creating species-specific interactomes, validating them functionally, and integrating additional molecular data to deepen the understanding of complex regulatory interactions in fish Evo-Devo. This short review aims to address the logical steps for this transition, as well as the necessities and limitations of this journey.

Neurotrophins and their receptors in the peripheral nervous system and non-nervous tissue of fish

Trophic factors, such as neurotrophins, are fundamental for cellular processes including differentiation, growth, survival, and regeneration. These molecules exhibit significant morphological and phylogenetic conservation throughout the animal kingdom, indicating conserved functions. In fish, the oldest and most diverse group of vertebrates, neurotrophins, and their receptors play pivotal roles not only within the central nervous system but also in various peripheral tissues. They are distributed in mechanosensory, muscle, skin, respiratory, circulatory, digestive, endocrine, urinary, reproductive, and immune systems, suggesting their involvement in the development and maintenance of all tissues/organs/systems. Despite this broad distribution, studies focusing on these molecules outside of the central nervous system have been limited to just 12 fish species. These investigations have revealed diverse expression patterns across different ages and tissues/organs/systems, expanding our comprehension of their functions beyond the central and peripheral nervous systems. Notably, BDNF and NT-3 are prominently expressed outside the central nervous system, particularly in mechanosensory and digestive tissues, whereas NGF is predominantly observed in mechanosensory and urinary systems. The expression and localization of neurotrophins and their receptors vary among organs, underscoring tissue-specific roles. Further research is imperative to decipher the precise functions and mechanisms of action of neurotrophins and their receptors in diverse fish tissues. Enhanced efforts are needed to include a broader range of fish species in these studies to advance our understanding of these agents in complex vertebrates, thereby shedding light on tissue development, regeneration, and maintenance, with potential implications for addressing organ-related issues.

Neural crest cells as a source of microevolutionary variation

Vertebrates have some of the most complex and diverse features in animals, from varied craniofacial morphologies to colorful pigmentation patterns and elaborate social behaviors. All of these traits have their developmental origins in a multipotent embryonic lineage of neural crest cells. This "fourth germ layer" is a vertebrate innovation and the source of a wide range of adult cell types. While others have discussed the role of neural crest cells in human disease and animal domestication, less is known about their role in contributing to adaptive changes in wild populations. Here, we review how variation in the development of neural crest cells and their derivatives generates considerable phenotypic diversity in nature. We focus on the broad span of traits under natural and sexual selection whose variation may originate in the neural crest, with emphasis on behavioral factors such as intraspecies communication that are often overlooked. In all, we encourage the integration of evolutionary ecology with developmental biology and molecular genetics to gain a more complete understanding of the role of this single cell type in trait covariation, evolutionary trajectories, and vertebrate diversity.

Positive selection as a key player for SARS-CoV-2 pathogenicity: Insights into ORF1ab, S and E genes

The human β-coronavirus SARS-CoV-2 epidemic started in late December 2019 in Wuhan, China. It causes Covid-19 disease which has become pandemic. Each of the five-known human β-coronaviruses has four major structural proteins (E, M, N and S) and 16 non-structural proteins encoded by ORF1a and ORF1b together (ORF1ab) that are involved in virus pathogenicity and infectivity. Here, we performed detailed positive selection analyses for those six genes among the four previously known human β-coronaviruses and within 38 SARS-CoV-2 genomes to assess signatures of adaptive evolution using maximum likelihood approaches. Our results suggest that three genes (E, S and ORF1ab genes) are under strong signatures of positive selection among human β-coronavirus, influencing codons that are located in functional important protein domains. The E protein-coding gene showed signatures of positive selection in two sites, Asp 66 and Ser 68, located inside a putative transmembrane α-helical domain C-terminal part, which is preferentially composed by hydrophilic residues. Such Asp and Ser sites substitutions (hydrophilic residues) increase the stability of the transmembrane domain in SARS-CoV-2. Moreover, substitutions in the spike (S) protein S1 N-terminal domain have been found, all of them were located on the S protein surface, suggesting their importance in viral transmissibility and survival. Furthermore, evidence of strong positive selection was detected in three of the SARS-CoV-2 nonstructural proteins (NSP1, NSP3, NSP16), which are encoded by ORF1ab and play vital roles in suppressing host translation machinery, viral replication and transcription and inhibiting the host immune response. These results are insightful to assess the role of positive selection in the SARS-CoV-2 encoded proteins, which will allow to better understand the virulent pathogenicity of the virus and potentially identifying targets for drug or vaccine strategy design.

Evolutionary genomics of mammalian lung cancer genes reveals signatures of positive selection in APC, RB1 and TP53

Lung cancer is the type of cancer causing most deaths in humans, with 234,030 new cases of lung cancer diagnosed in the United States in 2018. Recently, Tumor suppressor genes (TSGs) or the control of its pathway became promising drug targets for cancer therapy. A diverse group of TSGs is involved in progression and metastasis of lung cancer. Here, we surveyed nine highly significant mutated genes in 20 mammalian genomes to assess signatures of adaptive evolution using maximum likelihood approaches. We found that three genes (APC, RB1, and TP53) are under strong positive selection, influencing amino acids located in functionally important protein domains, such as three sites in APC found in the APC_N_CC domain, which is responsible for the binding to beta-catenin armadillo repeats that regulate beta-catenin level (beta-catenin is a transcription factor and its misregulation lead to malignant transformation of normal cells). Such sites substitutions mostly increase the stability of the domain. Moreover, substitution of some other sites found in important motifs, such as codon 47 (proline-directed kinase motif) in TP53, modify the phosphorylation activity of TP53 playing a key role in cancer risk. Our findings will open recommendation to drug targeting sites and will foster further research to understand better these proteins function.

Avian Binocularity and Adaptation to Nocturnal Environments: Genomic Insights from a Highly Derived Visual Phenotype

Typical avian eyes are phenotypically engineered for photopic vision (daylight). In contrast, the highly derived eyes of the barn owl (Tyto alba) are adapted for scotopic vision (dim light). The dramatic modifications distinguishing barn owl eyes from other birds include: 1) shifts in frontal orientation to improve binocularity, 2) rod-dominated retina, and 3) enlarged corneas and lenses. Some of these features parallel mammalian eye patterns, which are hypothesized to have initially evolved in nocturnal environments. Here, we used an integrative approach combining phylogenomics and functional phenotypes of 211 eye-development genes across 48 avian genomes representing most avian orders, including the stem lineage of the scotopic-adapted barn owl. Overall, we identified 25 eye-development genes that coevolved under intensified or relaxed selection in the retina, lens, cornea, and optic nerves of the barn owl. The agtpbp1 gene, which is associated with the survival of photoreceptor populations, was pseudogenized in the barn owl genome. Our results further revealed that barn owl retinal genes responsible for the maintenance, proliferation, and differentiation of photoreceptors experienced an evolutionary relaxation. Signatures of relaxed selection were also observed in the lens and cornea morphology-associated genes, suggesting that adaptive evolution in these structures was essentially structural. Four eye-development genes (ephb1, phactr4, prph2, and rs1) evolved in positive association with the orbit convergence in birds and under relaxed selection in the barn owl lineage, likely contributing to an increased reliance on binocular vision in the barn owl. Moreover, we found evidence of coevolutionary interactions among genes that are expressed in the retina, lens, and optic nerve, suggesting synergetic adaptive events. Our study disentangles the genomic changes governing the binocularity and low-light perception adaptations of barn owls to nocturnal environments while revealing the molecular mechanisms contributing to the shift from the typical avian photopic vision to the more-novel scotopic-adapted eye.

The Divergent Genomes of Teleosts

Boasting nearly 30,000 species, teleosts account for half of all extant vertebrates and approximately 98% of all ray-finned fish species (Actinopterygii). Teleosts are also the largest and most diverse group of vertebrates, exhibiting an astonishing level of morphological, physiological, and behavioral diversity. Previous studies had indicated that the teleost lineage has experienced an additional whole-genome duplication event. Recent comparative genomic analyses of teleosts and other bony vertebrates using spotted gar (a nonteleost ray-finned fish) and elephant shark (a cartilaginous fish) as outgroups have revealed several divergent features of teleost genomes. These include an accelerated evolutionary rate of protein-coding and nucleotide sequences, a higher rate of intron turnover, loss of many potential cis-regulatory elements and shorter conserved syntenic blocks. A combination of these divergent genomic features might have contributed to the evolution of the amazing phenotypic diversity and morphological innovations of teleosts.

Positive Selection Linked with Generation of Novel Mammalian Dentition Patterns

A diverse group of genes are involved in the tooth development of mammals. Several studies, focused mainly on mice and rats, have provided a detailed depiction of the processes coordinating tooth formation and shape. Here we surveyed 236 tooth-associated genes in 39 mammalian genomes and tested for signatures of selection to assess patterns of molecular adaptation in genes regulating mammalian dentition. Of the 236 genes, 31 (∼13.1%) showed strong signatures of positive selection that may be responsible for the phenotypic diversity observed in mammalian dentition. Mammalian-specific tooth-associated genes had accelerated mutation rates compared with older genes found across all vertebrates. More recently evolved genes had fewer interactions (either genetic or physical), were associated with fewer Gene Ontology terms and had faster evolutionary rates compared with older genes. The introns of these positively selected genes also exhibited accelerated evolutionary rates, which may reflect additional adaptive pressure in the intronic regions that are associated with regulatory processes that influence tooth-gene networks. The positively selected genes were mainly involved in processes like mineralization and structural organization of tooth specific tissues such as enamel and dentin. Of the 236 analyzed genes, 12 mammalian-specific genes (younger genes) provided insights on diversification of mammalian teeth as they have higher evolutionary rates and exhibit different expression profiles compared with older genes. Our results suggest that the evolution and development of mammalian dentition occurred in part through positive selection acting on genes that previously had other functions.

The Genome 10K Project: a way forward

The Genome 10K Project was established in 2009 by a consortium of biologists and genome scientists determined to facilitate the sequencing and analysis of the complete genomes of 10,000 vertebrate species. Since then the number of selected and initiated species has risen from ∼26 to 277 sequenced or ongoing with funding, an approximately tenfold increase in five years. Here we summarize the advances and commitments that have occurred by mid-2014 and outline the achievements and present challenges of reaching the 10,000-species goal. We summarize the status of known vertebrate genome projects, recommend standards for pronouncing a genome as sequenced or completed, and provide our present and future vision of the landscape of Genome 10K. The endeavor is ambitious, bold, expensive, and uncertain, but together the Genome 10K Consortium of Scientists and the worldwide genomics community are moving toward their goal of delivering to the coming generation the gift of genome empowerment for many vertebrate species.

Molecular evolution of multiple-level control of heme biosynthesis pathway in animal kingdom

Adaptation of enzymes in a metabolic pathway can occur not only through changes in amino acid sequences but also through variations in transcriptional activation, mRNA splicing and mRNA translation. The heme biosynthesis pathway, a linear pathway comprised of eight consecutive enzymes in animals, provides researchers with ample information for multiple types of evolutionary analyses performed with respect to the position of each enzyme in the pathway. Through bioinformatics analysis, we found that the protein-coding sequences of all enzymes in this pathway are under strong purifying selection, from cnidarians to mammals. However, loose evolutionary constraints are observed for enzymes in which self-catalysis occurs. Through comparative genomics, we found that in animals, the first intron of the enzyme-encoding genes has been co-opted for transcriptional activation of the genes in this pathway. Organisms sense the cellular content of iron, and through iron-responsive elements in the 5′ untranslated regions of mRNAs and the intron-exon boundary regions of pathway genes, translational inhibition and exon choice in enzymes may be enabled, respectively. Pathway product (heme)-mediated negative feedback control can affect the transport of pathway enzymes into the mitochondria as well as the ubiquitin-mediated stability of enzymes. Remarkably, the positions of these controls on pathway activity are not ubiquitous but are biased towards the enzymes in the upstream portion of the pathway. We revealed that multiple-level controls on the activity of the heme biosynthesis pathway depend on the linear depth of the enzymes in the pathway, indicating a new strategy for discovering the molecular constraints that shape the evolution of a metabolic pathway.

Mexican blind cavefish use mouth suction to detect obstacles

Fishes commonly use their lateral line system to detect moving bodies such as prey and predators. A remarkable case is the Mexican blind cavefish Astyanax fasciatus who evolved the ability to detect non-moving obstacles. The swimming body of A. fasciatus generates fluid disturbances, whose alteration by an obstacle can be sensed by the fish's lateral line system. It is generally accepted that these alterations can provide information on the distance to the obstacle. We observed that A. fasciatus swimming in an unfamiliar environment open and close their mouths at high frequency (0.7-4.5 Hz), in order to generate suction flows. We hypothesized that repeated mouth suction generate a hydrodynamic velocity field, whose alterations by an obstacle induce pressure gradients in the neuromasts of the lateral line, and corresponding strong lateral line stimuli. We observed that the frequency and rates of mouth opening events varied with the fish's distance to obstacles, a hallmark of pulse-based navigation mechanisms such as echolocation. We formulated a mathematical model of this hitherto unrecognized mechanism of obstacle detection and parameterized it experimentally. This model suggests that suction flows induce lateral line stimuli that are weakly dependent on the fish's speed, and may be an order of magnitude stronger than the correspondent stimuli induced by the fish's gliding body. We illustrate that A. fasciatus can navigate non-visually using a combination of two deeply ancestral and highly conserved mechanisms of ray-finned fishes: the mechanism of sensing water motion by the lateral line system and the mechanism of generating water motion by mouth suction.

EASER: Ensembl Easy Sequence Retriever

The rapid advances in genome sequencing technologies have increased the pace at which biological sequence databases are becoming available to the broad scientific community. Thus, obtaining and preparing an appropriate sequence dataset is a crucial first step for all types of genomic analyses. Here, we present a script that can widely facilitate the easy, fast, and effortless downloading and preparation of a proper biological sequence dataset for various genomics studies. This script retrieves Ensembl defined genomic features, associated with a given Ensembl identifier. Coding (CDS) and genomic sequences can be easily retrieved based on a selected relationship from a set of relationship types, either considering all available organisms or a user specified subset of organisms. The script is very user-friendly and by default starts with an interactive mode if no command-line options are specified.