Convergent copy number increase of genes associated with freshwater colonization in fishes

Draft Genome of Akame (Lates Japonicus) Reveals Possible Genetic Mechanisms for Long-Term Persistence and Adaptive Evolution with Low Genetic Diversity

It is known that some endangered species have persisted for thousands of years despite their very small effective population sizes and low levels of genetic polymorphisms. To understand the genetic mechanisms of long-term persistence in threatened species, we determined the whole genome sequences of akame (Lates japonicus), which has survived for a long time with extremely low genetic variations. Genome-wide heterozygosity in akame was estimated to be 3.3 to 3.4 × 10-4/bp, one of the smallest values in teleost fishes. Analysis of demographic history revealed that the effective population size in akame was around 1,000 from 30,000 years ago to the recent past. The relatively high ratio of nonsynonymous to synonymous heterozygosity in akame indicated an increased genetic load. However, a detailed analysis of genetic diversity in the akame genome revealed that multiple genomic regions, including genes involved in immunity, synaptic development, and olfactory sensory systems, have retained relatively high nucleotide polymorphisms. This implies that the akame genome has preserved the functional genetic variations by balancing selection, to avoid a reduction in viability and loss of adaptive potential. Analysis of synonymous and nonsynonymous nucleotide substitution rates has detected signs of positive selection in many akame genes, suggesting adaptive evolution to temperate waters after the speciation of akame and its close relative, barramundi (Lates calcarifer). Our results indicate that the functional genetic diversity likely contributed to the long-term persistence of this species by avoiding the harmful effects of the population size reduction.

Epigenetic diversity of genes with copy number variations among natural populations of the three-spined stickleback

Duplicated genes provide the opportunity for evolutionary novelty and adaptive divergence. In many cases, having more gene copies increases gene expression, which might facilitate adaptation to stressful or novel environments. Conversely, overexpression or misexpression of duplicated genes can be detrimental and subject to negative selection. In this scenario, newly duplicate genes may evade purifying selection if they are epigenetically silenced, at least temporarily, leading them to persist in populations as copy number variations (CNVs). In animals and plants, younger gene duplicates tend to have higher levels of DNA methylation and lower levels of gene expression, suggesting epigenetic regulation could promote the retention of gene duplications via expression repression or silencing. Here, we test the hypothesis that DNA methylation variation coincides with young duplicate genes that are segregating as CNVs in six populations of the three-spined stickleback that span a salinity gradient from 4 to 30 PSU. Using reduced-representation bisulfite sequencing, we found DNA methylation and CNV differentiation outliers rarely overlapped. Whereas lineage-specific genes and young duplicates were found to be highly methylated, just two gene CNVs showed a significant association between promoter methylation level and copy number, suggesting that DNA methylation might not interact with CNVs in our dataset. If most new duplications are regulated for dosage by epigenetic mechanisms, our results do not support a strong contribution from DNA methylation soon after duplication. Instead, our results are consistent with a preference to duplicate genes that are already highly methylated.

Evolutionary Rate Shifts in Coding and Regulatory Regions Underpin Repeated Adaptation to Sulfidic Streams in Poeciliid Fishes

Adaptation to extreme environments often involves the evolution of dramatic physiological changes. To better understand how organisms evolve these complex phenotypic changes, the repeatability and predictability of evolution, and possible constraints on adapting to an extreme environment, it is important to understand how adaptive variation has evolved. Poeciliid fishes represent a particularly fruitful study system for investigations of adaptation to extreme environments due to their repeated colonization of toxic hydrogen sulfide-rich springs across multiple species within the clade. Previous investigations have highlighted changes in the physiology and gene expression in specific species that are thought to facilitate adaptation to hydrogen sulfide-rich springs. However, the presence of adaptive nucleotide variation in coding and regulatory regions and the degree to which convergent evolution has shaped the genomic regions underpinning sulfide tolerance across taxa are unknown. By sampling across seven independent lineages in which nonsulfidic lineages have colonized and adapted to sulfide springs, we reveal signatures of shared evolutionary rate shifts across the genome. We found evidence of genes, promoters, and putative enhancer regions associated with both increased and decreased convergent evolutionary rate shifts in hydrogen sulfide-adapted lineages. Our analysis highlights convergent evolutionary rate shifts in sulfidic lineages associated with the modulation of endogenous hydrogen sulfide production and hydrogen sulfide detoxification. We also found that regions with shifted evolutionary rates in sulfide spring fishes more often exhibited convergent shifts in either the coding region or the regulatory sequence of a given gene, rather than both.

The prevalence of copy number increase at multiallelic copy number variants associated with cave colonization

Copy number variation is a common contributor to phenotypic diversity, yet its involvement in ecological adaptation is not easily discerned. Instances of parallelly evolving populations of the same species in a similar environment marked by strong selective pressures present opportunities to study the role of copy number variants (CNVs) in adaptation. By identifying CNVs that repeatedly occur in multiple populations of the derived ecotype and are not (or are rarely) present in the populations of the ancestral ecotype, the association of such CNVs with adaptation to the novel environment can be inferred. We used this paradigm to identify CNVs associated with recurrent adaptation of the Mexican tetra (Astyanax mexicanus) to cave environment. Using a read-depth approach, we detected CNVs from previously re-sequenced genomes of 44 individuals belonging to two ancestral surfaces and three derived cave populations. We identified 102 genes and 292 genomic regions that repeatedly diverge in copy number between the two ecotypes and occupy 0.8% of the reference genome. Functional analysis revealed their association with processes previously recognized to be relevant for adaptation, such as vision, immunity, oxygen consumption, metabolism, and neural function and we propose that these variants have been selected for in the cave or surface waters. The majority of the ecotype-divergent CNVs are multiallelic and display copy number increases in cavefish compared to surface fish. Our findings suggest that multiallelic CNVs - including gene duplications - and divergence in copy number provide a fast route to produce novel phenotypes associated with adaptation to subterranean life.

Genome-Wide Identification of Trachinotus ovatus Antimicrobial Peptides and Their Immune Response against Two Pathogen Challenges

Golden pompano, Trachinotus ovatus, as a highly nutritious commercially valuable marine fish, has become one of the preferred species for many fish farmers due to its rapid growth, wide adaptability, and ease of feeding and management. However, with the expansion of aquaculture scale, bacterial and parasitic diseases have also become major threats to the golden pompano industry. This study, based on comparative genomics, shows the possibility of preferential evolution of freshwater fish over marine fish by analyzing the phylogenetic relationships and divergence times of 14 marine fish and freshwater fish. Furthermore, we identified antimicrobial peptide genes from 14 species at the genomic level and found that the number of putative antimicrobial peptides may be related to species evolution. Subsequently, we classified the 341 identified AMPs from golden pompano into 38 categories based on the classification provided by the APD3. Among them, TCP represented the highest proportion, accounting for 23.2% of the total, followed by scolopendin, lectin, chemokine, BPTI, and histone-derived peptides. At the same time, the distribution of AMPs in chromosomes varied with type, and covariance analysis showed the frequency of its repeat events. Enrichment analysis and PPI indicated that AMP was mainly concentrated in pathways associated with disease immunity. In addition, our transcriptomic data measured the expression of putative AMPs of golden pompano in 12 normal tissues, as well as in the liver, spleen, and kidney infected with Streptococcus agalactiae and skin infected with Cryptocaryon irritans. As the infection with S. agalactiae and C. irritans progressed, we observed tissue specificity in the number and types of responsive AMPs. Positive selection of AMP genes may participate in the immune response through the MAPK signaling pathway. The genome-wide identification of antimicrobial peptides in the golden pompano provided a complete database of potential AMPs that can contribute to further understanding the immune mechanisms in pathogens. AMPs were expected to replace traditional antibiotics and be developed into targeted drugs against specific bacterial and parasitic pathogens for more precise and effective treatment to improve aquaculture production.

Functional characterization reveals a diverse array of metazoan fatty acid biosynthesis genes

Long-chain (≥C₂₀ ) polyunsaturated fatty acids (LC-PUFAs) are physiologically important fatty acids for most animals, including humans. Although most LC-PUFA production occurs in aquatic primary producers such as microalgae, recent research indicates the ability of certain groups of (mainly marine) invertebrates for endogenous LC-PUFA biosynthesis and/or bioconversion from dietary precursors. The genetic pathways for and mechanisms behind LC-PUFA biosynthesis remain unknown in many invertebrates to date, especially in non-model species. However, the numerous genomic and transcriptomic resources currently available can contribute to our knowledge of the LC-PUFA biosynthetic capabilities of metazoans. Within our previously generated transcriptome of the benthic harpacticoid copepod Platychelipus littoralis, we detected expression of one methyl-end desaturase, one front-end desaturase, and seven elongases, key enzymes responsible for LC-PUFA biosynthesis. To demonstrate their functionality, we characterized eight of them using heterologous expression in yeast. The P. littoralis methyl-end desaturase has Δ15/17/19 desaturation activity, enabling biosynthesis of α-linolenic acid, eicosapentaenoic acid and docosahexaenoic acid (DHA) from 18:2 n-6, 20:4 n-6 and 22:5 n-6, respectively. Its front-end desaturase has Δ4 desaturation activity from 22:5 n-3 to DHA, implying that P. littoralis has multiple pathways to produce this physiologically important fatty acid. All studied P. littoralis elongases possess varying degrees of elongation activity for saturated and unsaturated fatty acids, producing aliphatic hydrocarbon chains with lengths of up to 30 carbons. Our investigation revealed a functionally diverse range of fatty acid biosynthesis genes in copepods, which highlights the need to scrutinize the role that primary consumers could perform in providing essential nutrients to upper trophic levels.

Recent advances and current challenges in population genomics of structural variation in animals and plants

The field of population genomics has seen a surge of studies on genomic structural variation over the past two decades. These studies witnessed that structural variation is taxonomically ubiquitous and represent a dominant form of genetic variation within species. Recent advances in technology, especially the development of long-read sequencing platforms, have enabled the discovery of structural variants (SVs) in previously inaccessible genomic regions which unlocked additional structural variation for population studies and revealed that more SVs contribute to evolution than previously perceived. An increasing number of studies suggest that SVs of all types and sizes may have a large effect on phenotype and consequently major impact on rapid adaptation, population divergence, and speciation. However, the functional effect of the vast majority of SVs is unknown and the field generally lacks evidence on the phenotypic consequences of most SVs that are suggested to have adaptive potential. Non-human genomes are heavily under-represented in population-scale studies of SVs. We argue that more research on other species is needed to objectively estimate the contribution of SVs to evolution. We discuss technical challenges associated with SV detection and outline the most recent advances towards more representative reference genomes, which opens a new era in population-scale studies of structural variation.

Convergent copy number increase of genes associated with freshwater colonization in fishes

Copy number variation (CNV) can cause phenotypic changes. However, in contrast to amino acid substitutions and cis-regulatory changes, little is known about the functional categories of genes in which CNV is important for adaptation to novel environments. It is also unclear whether the same genes repeatedly change the copy numbers for adapting to similar environments. Here, we investigate CNV associated with freshwater colonization in fishes, which was observed multiple times across different lineages. Using 48 ray-finned fishes across diverse orders, we identified 23 genes whose copy number increases were associated with freshwater colonization. These genes showed enrichment for peptide receptor activity, hexosyltransferase activity and unsaturated fatty acid metabolism. We further revealed that three of the genes showed copy number increases in freshwater populations compared to marine ancestral populations of the stickleback genus Gasterosteus. These results indicate that copy number increases of genes involved in fatty acid metabolism (FADS2), immune function (PSMB8a) and thyroid hormone metabolism (UGT2) may be important for freshwater colonization at both the inter-order macroevolutionary scale and at the intra-genus microevolutionary scale. Further analysis across diverse taxa will help to understand the role of CNV in the adaptation to novel environments. This article is part of the theme issue 'Genetic basis of adaptation and speciation: from loci to causative mutations'.

Genetic basis of speciation and adaptation: from loci to causative mutations

Does evolution proceed in small steps or large leaps? How repeatable is evolution? How constrained is the evolutionary process? Answering these long-standing questions in evolutionary biology is indispensable for both understanding how extant biodiversity has evolved and predicting how organisms and ecosystems will respond to changing environments in the future. Understanding the genetic basis of phenotypic diversification and speciation in natural populations is key to properly answering these questions. The leap forward in genome sequencing technologies has made it increasingly easier to not only investigate the genetic architecture but also identify the variant sites underlying adaptation and speciation in natural populations. Furthermore, recent advances in genome editing technologies are making it possible to investigate the functions of each candidate gene in organisms from natural populations. In this article, we discuss how these recent technological advances enable the analysis of causative genes and mutations and how such analysis can help answer long-standing evolutionary biology questions.

This article is part of the theme issue ‘Genetic basis of adaptation and speciation: from loci to causative mutations’.