Evolutionary expression differences of creatine synthesis-related genes: Implications for skeletal muscle metabolism in fish

Integrating metabolomics and transcriptomics to analyze the differences of breast muscle quality and flavor formation between Daweishan mini chicken and broiler

The quality and flavor of chicken are affected by muscle metabolites and related regulatory genes, and the molecular regulation mechanism of meat quality is different among different breeds of chicken. In this study, 40 one-day-old Daweishan mini chicken (DM) and Cobb broiler (CB) were selected from each group, with 4 replicates and 10 chickens in each replicate. The chickens were reared until 90 d of age under the same management conditions. Then, metabolomics and transcriptomics data of 90-day-old DM (n = 4) and CB (n = 4) were integrated to analyze metabolites affecting breast muscle quality and flavor, and to explore the important genes regulating meat quality and flavor related metabolites. The results showed that a total of 38 significantly different metabolites (SDMs) and 420 differentially expressed genes (DEGs) were detected in the breast muscle of the 2 breeds. Amino acid and lipid metabolism may be the cause of meat quality and flavor difference between DM and CB chickens, involving metabolites such as L-methionine, betaine, N6, N6, N6-Trimethyl-L-lysine, L-anserine, glutathione, glutathione disulfide, L-threonine, N-Acetyl-L-aspartic acid, succinate, choline, DOPC, SOPC, alpha-linolenic acid, L-palmitoylcarnitine, etc. Important regulatory genes with high correlation with flavor amino acids (GATM, GSTO1) and lipids (PPARG, LPL, PLIN1, SCD, ANGPTL4, FABP7, GK, B4GALT6, UGT8, PLPP4) were identified by correlation analysis, and the gene-metabolite interaction network of breast muscle mass and flavor formation in DM chicken was constructed. This study showed that there were significant differences in breast metabolites between DM and CB chickens, mainly in amino acid and lipid metabolites. These 2 kinds of substances may be the main reasons for the difference in breast muscle quality and flavor between the 2 breeds. In general, this study could provide a theoretical basis for further research on the molecular regulatory mechanism of the formation of breast muscle quality and flavor differences between DM and CB chickens, and provide a reference for the development, utilization and genetic breeding of high-quality meat chicken breeds.

Generation and characterization of cytochrome P450 3A74 CRISPR/Cas9 knockout bovine foetal hepatocyte cell line (BFH12)

In human, the cytochrome P450 3A (CYP3A) subfamily of drug-metabolizing enzymes (DMEs) is responsible for a significant number of phase I reactions, with the CYP3A4 isoform superintending the hepatic and intestinal metabolism of diverse endobiotic and xenobiotic compounds. The CYP3A4-dependent bioactivation of chemicals may result in hepatotoxicity and trigger carcinogenesis. In cattle, four CYP3A genes (CYP3A74, CYP3A76, CYP3A28 and CYP3A24) have been identified. Despite cattle being daily exposed to xenobiotics (e.g., mycotoxins, food additives, drugs and pesticides), the existing knowledge about the contribution of CYP3A in bovine hepatic metabolism is still incomplete. Nowadays, CRISPR/Cas9 mediated knockout (KO) is a valuable method to generate in vivo and in vitro models for studying the metabolism of xenobiotics. In the present study, we successfully performed CRISPR/Cas9-mediated KO of bovine CYP3A74, human CYP3A4-like, in a bovine foetal hepatocyte cell line (BFH12). After clonal expansion and selection, CYP3A74 ablation was confirmed at the DNA, mRNA, and protein level. The subsequent characterization of the CYP3A74 KO clone highlighted significant transcriptomic changes (RNA-sequencing) associated with the regulation of cell cycle and proliferation, immune and inflammatory response, as well as metabolic processes. Overall, this study successfully developed a new CYP3A74 KO in vitro model by using CRISPR/Cas9 technology, which represents a novel resource for xenobiotic metabolism studies in cattle. Furthermore, the transcriptomic analysis suggests a key role of CYP3A74 in bovine hepatocyte cell cycle regulation and metabolic homeostasis.

Metabolomic and Lipidomic Tools for Tracing Fish Escapes from Aquaculture Facilities

During adverse atmospheric events, enormous damage can occur at marine aquaculture facilities, as was the case during Storm Gloria in the southeastern Spanish Mediterranean in January 2020, with massive fish escapes. Fishes that escape were caught by professional fishermen. The objective of this study was to identify biomarkers in fish that enable differentiation among wild fish, escaped farm-raised fish, and farm-raised fish kept in aquaculture facilities until their slaughter. We focused on gilthead sea bream (Sparus aurata). We used nuclear magnetic resonance to search for possible biomarkers. We found that wild gilthead sea bream showed higher levels of taurine and trimethylamine-N-oxide (TMAO) in their muscle and higher levels of ω-3 fatty acids, whereas farm-escaped and farmed gilthead sea bream raised until slaughter exhibit higher levels of ω-6 fatty acids. From choline, carnitine, creatinine, betaine, or lecithin, trimethylamine (TMA) is synthesized in the intestine by the action of bacterial microflora. In the liver, TMA is oxidized to TMAO and transported to muscle cells. The identified biomarkers will improve the traceability of gilthead sea bream by distinguishing wild specimens from those raised in aquaculture.

First report on the toxicity of SARS-CoV-2, alone and in combination with polyethylene microplastics in neotropical fish

The COVID-19 pandemic has caused unprecedented negative impacts in the modern era, including economic, social, and public health losses. On the other hand, the potential effects that the input of SARS-CoV-2 in the aquatic environment from sewage may represent on non-target organisms are not well known. In addition, it is not yet known whether the association of SARS-CoV-2 with other pollutants, such as microplastics (MPs), may further impact the aquatic biota. Thus, we aimed to evaluate the possible ecotoxicological effects of exposure of male adults Poecilia reticulata, for 15 days, to inactivated SARS-CoV-2 (0.742 pg/L; isolated SARS.CoV2/SP02.2020.HIAE.Br) and polyethylene MP (PE MPs) (7.1 × 10⁴ particles/L), alone and in combination, from multiple biomarkers. Our data suggest that exposure to SARS-CoV-2 induced behavioral changes (in the open field test), nephrotoxic effect (inferred by the increase in creatinine), hepatotoxic effect (inferred by the increase in bilirubin production), imbalance in the homeostasis of Fe, Ca, and Mg, as well as an anticholinesterase effect in the animals [marked by the reduction of acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) activity]. On the other hand, exposure to PE MPs induced a genotoxic effect (assessed by the comet assay), as well as an increase in enzyme activity alpha-amylase, alkaline phosphatase, and carboxylesterases. However, we did not show synergistic, antagonistic, or additive effects caused by the combined exposure of P. reticulata to SARS-CoV-2 and PE MPs. Principal component analysis (PCA) and values from the "Integrated Biomarker Response" index indicate that exposure to SARS-CoV-2 was determinant for a more prominent effect in the evaluated animals. Therefore, our study sheds light on the ecotoxicity of the new coronavirus in non-target organisms and ratifies the need for more attention to the impacts of COVID-19 on aquatic biota.

Comparative Analysis of miRNA-mRNA Regulation in the Testes of Gobiocypris rarus following 17α-Methyltestosterone Exposure

17α-Methyltestosterone (17MT), a synthetic organic compound commonly found in sewage waters, can affect reproduction in aquatic animals, such as tilapia and yellow catfish. In the present study, male Gobiocypris rarus were exposed to 25, 50, and 100 ng/L of 17α-methyltestosterone (17MT) for 7 days. We first analyzed miRNA- and RNA-seq results to determine miRNA-target gene pairs and then developed miRNA-mRNA interactive networks after 17MT administration. Total weights, total lengths, and body lengths were not significantly different between the test groups and control groups. The paraffin slice method was applied to testes of G. rarus in the MT exposure and control groups. We found that there were more mature sperm (S) and fewer secondary spermatocytes (SSs) and spermatogonia (SGs) in the testes of control groups. As 17MT concentration increased, fewer and fewer mature sperm (S) were observed in the testes of male G. rarus. The results showed that FSH, 11-KT, and E2 were significantly higher in individuals exposed to 25 ng/L 17MT compared with the control groups. VTG, FSH, LH, 11-KT, and E2 were significantly lower in the 50 ng/L 17MT exposure groups compared to the control groups. VTG, FSH, LH, 11-KT, E2, and T were significantly lower in the groups exposed to 100 ng/L 17MT. High-throughput sequencing revealed 73,449 unigenes, 1205 known mature miRNAs, and 939 novel miRNAs in the gonads of G. rarus. With miRNA-seq, 49 (MT25-M vs. Con-M), 66 (MT50-M vs. Con-M), and 49 (MT100-M vs. Con-M) DEMs were identified in the treatment groups. Five mature miRNAs (miR-122-x, miR-574-x, miR-430-y, lin-4-x, and miR-7-y), as well as seven differentially expressed genes (soat2, inhbb, ihhb, gatm, faxdc2, ebp, and cyp1a1), which may be associated with testicular development, metabolism, apoptosis, and disease response, were assayed using qRT-PCR. Furthermore, miR-122-x (related to lipid metabolism), miR-430-y (embryonic development), lin-4-x (apoptosis), and miR-7-y (disease) were differentially expressed in the testes of 17MT-exposed G. rarus. This study highlights the role of miRNA-mRNA pairs in the regulation of testicular development and immune response to disease and will facilitate future studies on the miRNA-RNA-associated regulation of teleost reproduction.

Supply and demand of creatine and glycogen in broiler chicken embryos

Optimal embryonic development and growth of meat-type chickens (broilers) rely on incubation conditions (oxygen, heat, and humidity), on nutrients and on energy resources within the egg. Throughout incubation and according to the embryo's energy balance, the main energy storage molecules (creatine and glycogen) are continuously utilized and synthesized, mainly in the embryonic liver, breast muscle, and the extraembryonic yolk sac (YS) tissue. During the last phase of incubation, as the embryo nears hatching, dynamic changes in energy metabolism occur. These changes may affect embryonic survival, hatchlings' uniformity, quality and post hatch performance of broilers, hence, being of great importance to poultry production. Here, we followed the dynamics of creatine and glycogen from embryonic day (E) 11 until hatch and up to chick placement at the farm. We showed that creatine is stored mainly in the breast muscle while glycogen is stored mainly in the YS tissue. Analysis of creatine synthesis genes revealed their expression in the liver, kidney, YS tissue and in the breast muscle, suggesting a full synthesis capacity in these tissues. Expression analysis of genes involved in gluconeogenesis, glycogenesis, and glycogenolysis, revealed that glycogen metabolism is most active in the liver. Nevertheless, due to the relatively large size of the breast muscle and YS tissue, their contribution to glycogen metabolism in embryos is valuable. Towards hatch, post E19, creatine levels in all tissues increased while glycogen levels dramatically decreased and reached low levels at hatch and at chick placement. This proves the utmost importance of creatine in energy supply to late-term embryos and hatchlings.

Fish Allergenicity Modulation Using Tailored Enriched Diets—Where Are We?

Food allergy is an abnormal immune response to specific proteins in a certain food. The chronicity, prevalence, and the potential fatality of food allergy, make it a serious socio-economic problem. Fish is considered the third most allergenic food in the world, affecting part of the world population with a higher incidence in children and adolescents. The main allergen in fish, responsible for the large majority of fish-allergic reactions in sensitized patients, is a small and stable calcium-binding muscle protein named beta-parvalbumin. Targeting the expression or/and the 3D conformation of this protein by adding specific molecules to fish diets has been the innovative strategy of some researchers in the fields of fish allergies and nutrition. This has shown promising results, namely when the apo-form of β-parvalbumin is induced, leading in the case of gilthead seabream to a 50% reduction of IgE-reactivity in fish allergic patients.

Effect of creatine and EDTA supplemented diets on European seabass (Dicentrarchus labrax) allergenicity, fish muscle quality and omics fingerprint

The relatively easy access to fish worldwide, alongside the increase of aquaculture production contributes to increased fish consumption which result in higher prevalence of respective allergies. Allergies to fish constitute a significant concern worldwide. β-parvalbumin is the main elicitor for IgE-mediated reactions. Creatine, involved in the muscle energy metabolism, and ethylenediamine tetraacetic acid (EDTA), a calcium chelator, are potential molecules to modulate parvalbumin. The purpose of this study was to test creatine (2, 5 and 8%) and EDTA (1.5, 3 and 4.5%) supplementation in fish diets to modulate β-parvalbumin expression and structure and its allergenicity in farmed European seabass (Dicentrarchus labrax) while assessing its effects on the end-product quality. Fish welfare and muscle quality parameters were evaluated by plasma metabolites, rigor mortis, muscle pH and sensory and texture analysis. Proteomics was used to assess alterations in muscle proteome profile and metabolic fingerprinting by Fourier transform infrared spectroscopy was used to assess the liver metabolic profile. In addition, IgE-reactivity to parvalbumin was analysed using fish allergic patient sera. Metabolic fingerprinting of liver tissue revealed no major alterations in infrared spectra with creatine supplementation, while with EDTA, only absorption bands characteristic of lipids were altered. Comparative proteomics showed up regulation of (tropo) myosin and phosphoglycerate mutase 2 with Creatine supplementation. In the case of EDTA proteomics showed up regulation of proteins involved in cellular and ion homeostasis. Allergenicity seems not to be modulated with creatine or EDTA supplementation as no decreased expression levels were found and IgE-binding reactivity showed no quantitative differences.

Metabolic Changes During Growth and Reproductive Phases in the Liver of Female Goldfish (Carassius auratus)

Hormones of the brain-pituitary-peripheral axis regulate metabolism, gonadal maturation, and growth in vertebrates. In fish, reproduction requires a significant energy investment to metabolically support the production of hundreds of eggs and billions of sperms in females and males, respectively. This study used an LC-MS-based metabolomics approach to investigate seasonally-related changes in metabolic profile and energy allocation patterns in female goldfish liver. We measured basal metabolic profile in female goldfish at three phases of the reproductive cycle, including 1) Maximum growth period in postovulatory regressed phase, 2) mid recrudescence in fish with developing follicles, and 3) late recrudescence when the ovary contains mature ovulatory follicles. We also investigated changes in the liver metabolism following acute treatments with GnRH and GnIH, known to be involved in controlling reproduction and growth in goldfish. Chemometrics combined with pathway-driven bioinformatics revealed significant changes in the basal and GnRH/GnIH-induced hepatic metabolic profile, indicating that metabolic energy allocation is regulated to support gonadal development and growth at different reproductive cycles. Overall, the findings support the hypothesis that hormonal control of reproduction involves accompanying metabolic changes to energetically support gonadotropic and somatotropic activities in goldfish and other oviparous vertebrates.

GATM and GAMT synthesize creatine locally throughout the mammalian body and within oligodendrocytes of the brain

The enzymes glycine amidinotransferase, mitochondrial (GATM also known as AGAT) and guanidinoacetate N-methyltransferase (GAMT) function together to synthesize creatine from arginine, glycine, and S-Adenosyl methionine. Deficiency in either enzyme or the creatine transporter, CT1, results in a devastating neurological disorder, Cerebral Creatine Deficiency Syndrome (CCDS). To better understand the pathophysiology of CCDS, we mapped the distribution of GATM and GAMT at single cell resolution, leveraging RNA sequencing analysis combined with in vivo immunofluorescence (IF). Using the mouse as a model system, we find that GATM and GAMT are coexpressed in several tissues with distinct and overlapping cellular sources, implicating local synthesis as an important mechanism of creatine metabolism in numerous organs. Extending previous findings at the RNA level, our analysis demonstrates that oligodendrocytes express the highest level of Gatm and Gamt of any cell type in the body. We confirm this finding in the mouse brain by IF, where GATM localizes to the mitochondria of oligodendrocytes, whereas both oligodendrocytes and cerebral cortical neurons express GAMT. Interestingly, the latter is devoid of GATM. Single nucleus assay for transposase-accessible chromatin sequencing (snATAC-seq) analysis of 4 brain regions highlights a similar primacy of oligodendrocytes in the expression of GATM and GAMT in the human central nervous system. Importantly, an active putative regulatory element within intron 2 of human GATM is detected in oligodendrocytes but not neurons.

Arginine metabolism and its functions in growth, nutrient utilization, and immunonutrition of fish

Fish have limited ability in endogenous biosynthesis of arginine. Arginine is an indispensable amino acid for fish, and the arginine requirement varies with fish species and fish size. Recent studies on fish have demonstrated that arginine influences nutrient metabolism, stimulates insulin release, is involved in nonspecific immune responses and antioxidant responses, and elevates disease resistance. Specifically, arginine can regulate energy homeostasis via modulating the adenosine 5'-monophosphate (AMP)-activated protein kinase (AMPK) pathway, and also regulate protein synthesis via activating the target of rapamycin (TOR) signaling pathway. The present article reviews pertinent knowledge of arginine in fish, including dietary quantitative requirements, endogenous anabolism and catabolism, regulation of the endocrine and metabolic systems, and immune-regulatory functions under pathogenic challenge. Our findings showed that further data about the distribution of arginine after intake into specific cells, its sub-cellular sensor to initiate downstream signaling pathways, and its effects on fish mucosal immunity, especially the adaptive immune response against pathogenic infection in different species, are urgently needed.

Proteomic profile and morphological characteristics of skeletal muscle from the fast- and slow-growing yellow perch (Perca flavescens)

The objective of the present study was to compare skeletal muscle proteomic profiles, histochemical characteristics, and expression levels of myogenic regulatory factors (MRFs) between fast- versus slow-growing yellow perch Perca flavescens and identify the proteins/peptides that might play a crucial role in the muscle growth dynamic. Yellow perch were nursed in ponds for 6 weeks from larval stage and cultured in two meter diameter tanks thereafter. The fingerlings were graded to select the top 10% and bottom 10% fish which represented fast- and slow-growing groups (31 yellow perch per each group). Our statistical analyses showed 18 proteins that had different staining intensities between fast- and slow-growing yellow perch. From those proteins 10 showed higher expression in slow-growers, and 8 demonstrated higher expression in fast-growers. Fast-growing yellow perch with a greater body weight was influenced by both the muscle fiber hypertrophy and mosaic hyperplasia compared to slow-growing fish. These hyperplastic and hypertrophic growth in fast-grower were associated with not only metabolic enzymes, including creatine kinase, glycogen phosphorylase, and aldolase, but also myoD and myogenin as MRFs. Overall, the results of the present study contribute to the identification of different expression patterns of gene products in fast- and slow-growing fish associated with their muscle growth.

Transcriptome-wide study revealed m6A regulation of embryonic muscle development in Dingan goose (Anser cygnoides orientalis)

BACKGROUND

The number of myofiber is determined during the embryonic stage and does not increase during the postnatal period for birds, including goose. Thus, muscle production of adult goose is pre-determined during embryogenesis. Previous studies show N6-methyladenosine (m6A) is an important regulator for skeletal muscle development of birds and miRNAs play as a co-regulator for the skeletal muscle development in birds. Herein, we sequenced m6A and miRNA transcriptomes to investigate the profiles of m6A and their potential mechanism of regulating breast muscle development in Dingan Goose.

RESULTS

We selected embryonic 21th day (E21) and embryonic 30th day (E30) to investigate the roles of transcriptome-wide m6A modification combining with mRNAs and miRNAs in goose breast muscle development. In this study, m6A peaks were mainly enriched in coding sequence (CDS) and start codon and397 genes were identified as differentially methylated genes (DMGs). GO and KEGG analysis showed that DMGs were highly related to cellular and metabolic process and that most DMGs were enriched in muscle-related pathways including Wnt signaling pathway, mTOR signaling and FoxO signaling pathway. Interestingly, a negative correlation between m6A methylation level and mRNA abundance was found through the analysis of m6A-RNA and RNA-seq data. Besides, we found 26 muscle-related genes in 397 DMGs. We also detected 228 differentially expressed miRNAs (DEMs), and further found 329 genes shared by the target genes of DEMs and DMGs (m6A-miRNA-genes), suggesting a tightly relationship between DEMs and DMGs. Among the m6A-miRNA-genes, we found 10 genes are related to breast muscle development. We further picked out an m6A-miRNA-gene, PDK3, from the 10 genes to visualize it and the result showed differentially methylated peaks on the mRNA transcript consistent with our m6A-seq results.

CONCLUSION

GO and KEGG of DMGs between E21 and E30 showed most DMGs were muscle-related. In total, 228 DEMs were found, and the majority of DMGs were overlapped with the targets of DEGs. The differentially methylated peaks along with an m6A-miRNA-gene, PDK3, showed the similar results with m6A-seq results. Taken together, the results presented here provide a reference for further investigation of embryonic skeletal muscle development mechanism in goose.

Insights into early ontogenesis: characterization of stress and development key genes of pikeperch (Sander lucioperca) in vivo and in vitro

There are still numerous difficulties in the successful farming of pikeperch in the anthropogenic environment of various aquaculture systems, especially during early developmental steps in the hatchery. To investigate the physiological processes involved on the molecular level, we determined the basal expression patterns of 21 genes involved in stress and immune responses and early ontogenesis of pikeperch between 0 and 175 days post hatch (dph). Their transcription patterns most likely reflect the challenges of growth and feed conversion. The gene coding for apolipoprotein A (APOE) was strongly expressed at 0 dph, indicating its importance for yolk sac utilization. Genes encoding bone morphogenetic proteins 4 and 7 (BMP4, BMP7), creatine kinase M (CKM), and SRY-box transcription factor 9 (SOX9) were highly abundant during the peak phases of morphological changes and acclimatization processes at 4-18 dph. The high expression of genes coding for peroxisome proliferator-activated receptors alpha and delta (PPARA, PPARD) at 121 and 175 dph, respectively, suggests their importance during this strong growth phase of juvenile stages. As an alternative experimental model to replace further in vivo investigations of ontogenetically important processes, we initiated the first approach towards a long-lasting primary cell culture from whole pikeperch embryos. The present study provides a set of possible biomarkers to support the monitoring of pikeperch farming and provides a first basis for the establishment of a suitable cell model of this emerging aquaculture species.

Dietary arginine and citrulline supplementation modulates the immune condition and inflammatory response of European seabass

The present study was designed to determine the modulatory effects of arginine and citrulline dietary supplementation on the immune condition and inflammatory response of European seabass, Dicentrarchus labrax. Four diets were manufactured: a control diet (CTRL) was formulated to meet the indispensable amino acids profile established for seabass. Based on this formulation, three other diets were supplemented with l-arginine at two different levels (0.5% and 1%, ARG1 and ARG2, respectively) and l-citrulline at 0.5% (CIT). Fish were fed these diets for 2 or 4 weeks under controlled conditions. At the end of 4 weeks, fish from all dietary treatments were intraperitoneally-injected with Photobacterium damselae piscicida and sampled after 4, 24 our 48 h. Immune status was characterized by a lymphocyte time-dependent decrease regardless of dietary treatment, whereas peroxidase values dropped in time in fish fed ARG1 and ARG2 and was lower at 4 weeks in fish fed ARG1 than in fish fed CTRL. Up-regulation of several genes was more evident in ARG1-and CIT-fed fish, though pro-inflammatory cytokines were down-regulated by CIT dietary treatment. Following immune stimulation, seabass fed ARG1 showed a decrease in neutrophils and monocytes circulating numbers. On the other hand, expression of 17 selected immune and inflammatory responses genes was barely affected by dietary treatments. Based on the analyzed parameters, results suggest an active role of dietary arginine/citrulline supplementation in modulating immune defences that seem to translate into a suppressed immune repertoire, mostly at the cell response level. The observed changes due to citrulline dietary supplementation were in part similar to those caused by arginine, suggesting that citrulline might have been used by macrophages as an arginine precursor and then engaged in similar immune-impairment leading mechanisms.

Parallel selection on ecologically relevant gene functions in the transcriptomes of highly diversifying salmonids

Background

Salmonid fishes are characterised by a very high level of variation in trophic, ecological, physiological, and life history adaptations. Some salmonid taxa show exceptional potential for fast, within-lake diversification into morphologically and ecologically distinct variants, often in parallel; these are the lake-resident charr and whitefish (several species in the genera Salvelinus and Coregonus). To identify selection on genes and gene categories associated with such predictable diversifications, we analysed 2702 orthogroups (4.82 Mbp total; average 4.77 genes/orthogroup; average 1783 bp/orthogroup). We did so in two charr and two whitefish species and compared to five other salmonid lineages, which do not evolve in such ecologically predictable ways, and one non-salmonid outgroup.

Results

All selection analyses are based on Coregonus and Salvelinus compared to non-diversifying taxa. We found more orthogroups were affected by relaxed selection than intensified selection. Of those, 122 were under significant relaxed selection, with trends of an overrepresentation of serine family amino acid metabolism and transcriptional regulation, and significant enrichment of behaviour-associated gene functions. Seventy-eight orthogroups were under significant intensified selection and were enriched for signalling process and transcriptional regulation gene ontology terms and actin filament and lipid metabolism gene sets. Ninety-two orthogroups were under diversifying/positive selection. These were enriched for signal transduction, transmembrane transport, and pyruvate metabolism gene ontology terms and often contained genes involved in transcriptional regulation and development. Several orthogroups showed signs of multiple types of selection. For example, orthogroups under relaxed and diversifying selection contained genes such as ap1m2, involved in immunity and development, and slc6a8, playing an important role in muscle and brain creatine uptake. Orthogroups under intensified and diversifying selection were also found, such as genes syn3, with a role in neural processes, and ctsk, involved in bone remodelling.

Conclusions

Our approach pinpointed relevant genomic targets by distinguishing among different kinds of selection. We found that relaxed, intensified, and diversifying selection affect orthogroups and gene functions of ecological relevance in salmonids. Because they were found consistently and robustly across charr and whitefish and not other salmonid lineages, we propose these genes have a potential role in the replicated ecological diversifications.