Drivers of interlineage variability in mitogenomic evolutionary rates in Platyhelminthes.

Studies of forces driving interlineage variability in the evolutionary rates (both sequence and architecture) of mitochondrial genomes often produce contradictory results. Flatworms (Platyhelminthes) exhibit the fastest-evolving mitogenomic sequences among all bilaterian phyla. To test the effects of multiple factors previously associated with different aspects of mitogenomic evolution, we used mitogenomes of 223 flatworm species, phylogenetic multilevel regression models, and causal inference. Thermic host environment (endothermic vs. ectothermic) had nonsignificant impacts on both sequence evolution and mitogenomic size. Mitogenomic gene order rearrangements (GORR) were mostly positively correlated with mitogenomic size (R2 ≈ 20-30%). Longevity was not (negatively) correlated with sequence evolution in flatworms. The predominantly free-living "turbellaria" exhibited much shorter branches and faster-evolving mitogenomic architecture than parasitic Neodermata. As a result, "parasitism" had a strong explanatory power on the branch length variability (>90%), and there was a negative correlation between GORR and branch length. However, the stem branch of Neodermata comprised 63.6% of the total average branch length. This evolutionary period was also marked by a high rate of gene order rearrangements in the ancestral Neodermata. We discuss how this period of rapid evolution deep in the evolutionary history may have decoupled sequence evolution rates from longevity and GORR, and overestimated the explanatory power of "parasitism". This study shows that impacts of variables often vary across lineages, and stresses the importance accounting for the episodic nature of evolutionary patterns in studies of mitogenomic evolution.

The mitochondrial genome of Heterosentis pseudobagri (Wang & Zhang, 1987) Pichelin & Cribb, 1999 reveals novel aspects of tRNA genes evolution in Acanthocephala.

BACKGROUND: Acanthocephala is a clade of obligate endoparasites whose mitochondrial genomes (mitogenomes) and evolution remain relatively poorly understood. Previous studies reported that atp8 is lacking from acanthocephalan mitogenomes, and that tRNA genes often have nonstandard structures. Heterosentis pseudobagri (Arhythmacanthidae) is an acanthocephalan fish endoparasite for which no molecular data are currently available, and biological information is unavailable in the English language. Furthermore, there are currently no mitogenomes available for Arhythmacanthidae. METHODS: We sequenced its mitogenome and transcriptome, and conducted comparative mitogenomic analyses with almost all available acanthocephalan mitogenomes. RESULTS: The mitogenome had all genes encoded on the same strand and unique gene order in the dataset. Among the 12 protein-coding genes, several genes were highly divergent and annotated with difficulty. Moreover, several tRNA genes could not be identified automatically, so we had to identify them manually via a detailed comparison with orthologues. As common in acanthocephalans, some tRNAs lacked either the TWC arm or the DHU arm, but in several cases, we annotated tRNA genes only on the basis of the conserved narrow central segment comprising the anticodon, while the flanking 5' and 3' ends did not exhibit any resemblance to orthologues and they could not be folded into a tRNA secondary structure. We corroborated that these are not sequencing artefacts by assembling the mitogenome from transcriptomic data. Although this phenomenon was not observed in previous studies, our comparative analyses revealed the existence of highly divergent tRNAs in multiple acanthocephalan lineages. CONCLUSIONS: These findings indicate either that multiple tRNA genes are non-functional or that (some) tRNA genes in (some) acanthocephalans might undergo extensive posttranscriptional tRNA processing which restores them to more conventional structures. It is necessary to sequence mitogenomes from yet unrepresented lineages and further explore the unusual patterns of tRNA evolution in Acanthocephala.

Whole-genome sequencing identifies rare missense variants of WNT16 and ERVW-1 causing the systemic lupus erythematosus.

Systemic lupus erythematosus (SLE, OMIM 152700) is a rare autoimmune disease with high heritability that affects ~0.1% of the population. Previous studies have revealed several common variants with small effects in European and East Asian SLE patients. However, there is still no rare variant study on Chinese SLE patients using the whole-genome sequencing technology (WGS). Here, we designed a family based WGS study to identify novel rare variants with large effects. Based on large-scale allele frequency data from the gnomAD database, we identified rare protein-coding gene variants with disruptive and sequence-altering impacts in SLE patients. We found that the burden of rare variants was significantly higher than that of common variants in patients, suggesting a larger effect of rare variants on the SLE pathogenesis. We identified the pathogenic risk of rare missense variants with significant odds ratios (p < 0.05) in two genes, including WNT16 (NC_000007.14:g.121329757G > C, NP_057171.2:p.(Ala86Pro) and 7 g.121329760G > C, NP_057171.2:p.(Ala87Pro)), which explains five out of seven patients covering all three families but are absent from all controls, and ERVW-1 (NC_000007.14:g.92469882A > G, NP_001124397.1:p.(Leu167Pro), rs74545114; NC_000007.14:g.92469907G > A, NP_001124397.1:p.(Arg159Cys), rs201142302; NC_000007.14:g.92469919G > A, NP_001124397.1:p.(His155Tyr), rs199552228), which explains the other two patients. None of these variants were identified in any of the controls. These associations are supported by known gene expression studies in SLE patients based on literature review. We further tested the wild and mutant types using the luciferase assays and qPCR in cells. We found that WNT16 can activate the canonical Wnt/ß-catenin pathway while the mutant cannot. Additionally, the wild ERVW-1 expression can be significantly up-regulated by cAMP while the mutant cannot. Our study provides the first direct genetic and in vitro evidence for the pathogenic risk of mutant WNT16 and ERVW-1, which may facilitate the design of precision therapy for SLE.

Inverted base composition skews and discontinuous mitochondrial genome architecture evolution in the Enoplea (Nematoda).

BACKGROUND: Within the class Enoplea, the earliest-branching lineages in the phylum Nematoda, the relatively highly conserved ancestral mitochondrial architecture of Trichinellida is in stark contrast to the rapidly evolving architecture of Dorylaimida and Mermithida. To better understand the evolution of mitogenomic architecture in this lineage, we sequenced the mitogenome of a fish parasite Pseudocapillaria tomentosa (Trichinellida: Capillariidae) and compared it to all available enoplean mitogenomes. RESULTS: P. tomentosa exhibited highly reduced noncoding regions (the largest was 98 bp), and a unique base composition among the Enoplea. We attributed the latter to the inverted GC skew (0.08) in comparison to the ancestral skew in Trichinellidae (-0.43 to -0.37). Capillariidae, Trichuridae and Longidoridae (Dorylaimida) generally exhibited low negative or low positive skews (-0.1 to 0.1), whereas Mermithidae exhibited fully inverted low skews (0 to 0.05). This is indicative of inversions in the strand replication order or otherwise disrupted replication mechanism in the lineages with reduced/inverted skews. Among the Trichinellida, Trichinellidae and Trichuridae have almost perfectly conserved architecture, whereas Capillariidae exhibit multiple rearrangements of tRNA genes. In contrast, Mermithidae (Mermithida) and Longidoridae (Dorylaimida) exhibit almost no similarity to the ancestral architecture. CONCLUSIONS: Longidoridae exhibited more rearranged mitogenomic architecture than the hypervariable Mermithidae. Similar to the Chromadorea, the evolution of mitochondrial architecture in enoplean nematodes exhibits a strong discontinuity: lineages possessing a mostly conserved architecture over tens of millions of years are interspersed with lineages exhibiting architectural hypervariability. As Longidoridae also have some of the smallest metazoan mitochondrial genomes, they contradict the prediction that compact mitogenomes should be structurally stable. Lineages exhibiting inverted skews appear to represent the intermediate phase between the Trichinellidae (ancestral) and fully derived skews in Chromadorean mitogenomes (GC skews = 0.18 to 0.64). Multiple lines of evidence (CAT-GTR analysis in our study, a majority of previous mitogenomic results, and skew disruption scenarios) support the Dorylaimia split into two sister-clades: Dorylaimida + Mermithida and Trichinellida. However, skew inversions produce strong base composition biases, which can hamper phylogenetic and other evolutionary studies, so enoplean mitogenomes have to be used with utmost care in evolutionary studies.

A chromosome-level genome assembly of Cairina moschata and comparative genomic analyses.

BACKGROUND: The Muscovy duck (Cairina moschata) is an economically important duck species, with favourable growth and carcass composition parameters in comparison to other ducks. However, limited genomic resources for Muscovy duck hinder our understanding of its evolution and genetic diversity. RESULTS: We combined linked-reads sequencing technology and reference-guided methods for de novo genome assembly. The final draft assembly was 1.12 Gbp with 29 autosomes, one sex chromosome and 4,583 unlocalized scaffolds with an N50 size of 77.35 Mb. Based on universal single-copy orthologues (BUSCO), the draft genome assembly completeness was estimated to be 93.30 %. Genome annotation identified 15,580 genes, with 15,537 (99.72 %) genes annotated in public databases. We conducted comparative genomic analyses and found that species-specific and rapidly expanding gene families (compared to other birds) in Muscovy duck are mainly involved in Calcium signaling, Adrenergic signaling in cardiomyocytes, and GnRH signaling pathways. In comparison to the common domestic duck (Anas platyrhynchos), we identified 104 genes exhibiting strong signals of adaptive evolution (Ka/Ks > 1). Most of these genes were associated with immune defence pathways (e.g. IFNAR1 and TLR5). This is indicative of the existence of differences in the immune responses between the two species. Additionally, we combined divergence and polymorphism data to demonstrate the "faster-Z effect" of chromosome evolution. CONCLUSIONS: The chromosome-level genome assembly of Muscovy duck and comparative genomic analyses provide valuable resources for future molecular ecology studies, as well as the evolutionary arms race between the host and influenza viruses.

Slow crabs - fast genomes: Locomotory capacity predicts skew magnitude in crustacean mitogenomes.

Base composition skews (G-C/G+C) of mitochondrial genomes are believed to be primarily driven by mutational pressure, which is positively correlated with metabolic rate. In marine animals, metabolic rate is also positively correlated with locomotory capacity. Given the central role of mitochondria in energy metabolism, we hypothesised that selection for locomotory capacity should be positively correlated with the strength of purifying selection (dN/dS), and thus be negatively correlated with the skew magnitude. Therefore, these two models assume diametrically opposite associations between the metabolic rate and skew magnitude: positive correlation in the prevailing paradigm, and negative in our working hypothesis. We examined correlations between the skew magnitude, metabolic rate, locomotory capacity, and several other variables previously associated with mitochondrial evolution on 287 crustacean mitogenomes. Weakly locomotory taxa had higher skew magnitude and ω (dN/dS) values, but not the gene order rearrangement rate. Skew and ω magnitudes were correlated. Multilevel regression analyses indicated that three competing variables, body size, gene order rearrangement rate, and effective population size, had negligible impacts on the skew magnitude. In most crustacean lineages selection for locomotory capacity appears to be the primary factor determining the skew magnitude. Contrary to the prevailing paradigm, this implies that adaptive selection outweighs nonadaptive selection (mutation pressure) in crustaceans. However, we found indications that effective population size (nonadaptive factor) may outweigh the impact of locomotory capacity in sessile crustaceans (Thecostraca). In conclusion, skew magnitude is a product of the interplay between adaptive and nonadaptive factors, the balance of which varies among lineages.

Evolutionary history of inversions in directional mutational pressures in crustacean mitochondrial genomes: Implications for evolutionary studies.

Inversions of the origin of replication (ORI) in mitochondrial genomes produce asymmetrical mutational pressures that can cause strong base composition skews. Due to skews often being overlooked, the total number of crustacean lineages that underwent ORI events remains unknown. We analysed skews, cumulative skew plots, conserved sequence motifs, and mitochondrial architecture of all 965 available crustacean mitogenomes (699 unique species). We found indications of an ORI in 159 (22.7%) species, and mapped these to 23 ORI events: 16 identified with confidence and 7 putative (13 newly proposed, and for 5 we improved the resolution). Two ORIs occurred at or above the order level: Isopoda and Copepoda. Shifts in skew plots are not a precise tool for identifying the replication mechanism. We discuss how ORIs can produce mutational bursts in mitogenomes and show how these can interfere with various types of evolutionary studies. Phylogenetic analyses were plagued by artefactual clustering, and ORI lineages exhibited longer branches, a higher number of synonymous substitutions, higher mutational saturation, and higher compositional heterogeneity. ORI events also affected codon usage and protein properties. We discuss how this may have caused erroneous interpretation of data in previous studies that did not account for skew patterns.

Disrupted architecture and fast evolution of the mitochondrial genome of Argeia pugettensis (Isopoda): implications for speciation and fitness.

BACKGROUND: Argeia pugettensis is an isopod species that parasitizes other crustaceans. Its huge native geographic range spans the Pacific from China to California, but molecular data are available only for a handful of specimens from North-American populations. We sequenced and characterised the complete mitogenome of a specimen collected in the Yellow Sea. RESULTS: It exhibited a barcode (cox1) similarity level of only 87-89% with North-American populations, which is unusually low for conspecifics. Its mitogenome is among the largest in isopods (≈16.5 Kbp), mostly due to a large duplicated palindromic genomic segment (2 Kbp) comprising three genes. However, it lost a segment comprising three genes, nad4L-trnP-nad6, and many genes exhibited highly divergent sequences in comparison to isopod orthologues, including numerous mutations, deletions and insertions. Phylogenetic and selection analyses corroborated that this is one of the handful of most rapidly evolving available isopod mitogenomes, and that it evolves under highly relaxed selection constraints (as opposed to positive selection). However, its nuclear 18S gene is highly conserved, which suggests that rapid evolution is limited to its mitochondrial genome. The cox1 sequence analysis indicates that elevated mitogenomic evolutionary rates are not shared by North-American conspecifics, which suggests a breakdown of cox1 barcoding in this species. CONCLUSIONS: A highly architecturally disrupted mitogenome and decoupling of mitochondrial and nuclear rates would normally be expected to have strong negative impacts on the fitness of the organism, so the existence of this lineage is a puzzling evolutionary question. Additional studies are needed to assess the phylogenetic breadth of this disrupted mitochondrial architecture and its impact on fitness.

Gene expression patterns indicate that a high-fat-high-carbohydrate diet causes mitochondrial dysfunction in fish.

Expensive and unsustainable fishmeal is increasingly being replaced with cheaper lipids and carbohydrates as sources of energy in aquaculture. Although it is known that the excess of lipids and carbohydrates has negative effects on nutrient utilization, growth, metabolic homeostasis, and health of fish, our current understanding of mechanisms behind these effects is limited. To improve the understanding of diet-induced metabolic disorders (both in fish and other vertebrates), we conducted an eight-week high-fat-high-carbohydrate diet feeding trial on blunt snout bream (Megalobrama amblycephala), and studied gene expression changes (transcriptome and qPCR) in the liver. Disproportionately large numbers of differentially expressed genes were associated with mitochondrial metabolism, neurodegenerative diseases (Alzheimer's, Huntington's, and Parkinson's), and functional categories indicative of liver dysfunction. A high-fat-high-carbohydrate diet may have caused mitochondrial dysfunction, and possibly downregulated the mitochondrial biogenesis in the liver. While the relationship between diet and neurodegenerative disorders is well-established in mammals, this is the first report of this connection in fish. We propose that fishes should be further explored as a potentially promising model to study the mechanisms of diet-associated neurodegenerative disorders in humans.

Identification, characterization and expression in response to Aeromonas hydrophila challenge of five interferon regulatory factors in Megalobrama amblycephala.

Interferon regulatory factor (Irf) family represents one of the most important transcription factor families, with multiple biological roles. In this study, we characterized five Irf family members (irf4a, irf4b, irf6, irf8 and irf10) in Megalobrama amblycephala at the cDNA and (predicted) amino acid levels, analyzed them phylogenetically, and developed gene-specific primers for qPCR analysis. All five irfs were constitutively expressed in all examined tissues, but their transcription was significantly higher in lymphoid organs and tissues, such as kidney, spleen and intestine. Exceptions were irf8, which was expressed at a high level in heart and brain tissues, and irf6, expressed at low levels in most tissues. After a bacterial immune challenge with Aeromonas hydrophila, the expression of irfs in liver was up-regulated: mairf4a 8.12-fold, mairf4b 29.9-fold, mairf6 1.38-fold and mairf10 1.65-fold (mairf8 was an exception: 0.07-fold). In spleen, kidney, intestine and gills, transcript levels of studied irfs increased only at specific time-points. The results suggested that irfs are involved in the immune response to bacterial infection in M. amblycephala, which will help elucidate the biological functions of irfs in the immune system of teleost fish.

Mitochondrial Genomes of Two Thaparocleidus Species (Platyhelminthes: Monogenea) Reveal the First rRNA Gene Rearrangement among the Neodermata.

Phylogenetic framework for the closely related Ancylodiscoidinae and Ancyrocephalinae subfamilies remains contentious. As this issue was never studied using a large molecular marker, we sequenced the first two Ancylodiscoidinae mitogenomes: Thaparocleidus asoti and Thaparocleidus varicus. Both mitogenomes had two non-coding regions (NCRs) that contained a number of repetitive hairpin-forming elements (RHE). Due to these, the mitogenome of T. asoti (16,074 bp) is the longest among the Monogenea; especially large is its major NCR, with 3500 bp, approximately 1500 bp of which could not be sequenced (thus, the total mitogenome size is ≈ 17,600 bp). Although RHEs have been identified in other monopisthocotyleans, they appear to be independently derived in different taxa. The presence of RHEs may have contributed to the high gene order rearrangement rate observed in the two mitogenomes, including the first report of a transposition of rRNA genes within the Neodermata. Phylogenetic analyses using mitogenomic dataset produced Dactylogyrinae embedded within the Ancyrocephalinae (paraphyly), whereas Ancylodiscoidinae formed a sister-group with them. This was also supported by the gene order analysis. 28S rDNA dataset produced polyphyletic Dactylogyridae and Ancyrocephalinae. The phylogeny of the two subfamilies shall have to be further evaluated with more data.

Population size may shape the accumulation of functional mutations following domestication.

BACKGROUND: Population genetics theory predicts an important role of differences in the effective population size (N e ) among species on shaping the accumulation of functional mutations by regulating the selection efficiency. However, this correlation has never been tested in domesticated animals. RESULTS: Here, we synthesized 62 whole genome data in eight domesticated species (cat, dog, pig, goat, sheep, chicken, cattle and horse) and compared domesticates with their wild (or ancient) relatives. Genes with significantly different selection pressures (revealed by nonsynonymous/synonymous substitution rate ratios, Ka/Ks or ω) between domesticated (Dω) and wild animals (Wω) were determined by likelihood-ratio tests. Species-level effective population sizes (N e ) were evaluated by the pairwise sequentially Markovian coalescent (PSMC) model, and Dω/Wω were calculated for each species to evaluate the changes in accumulation of functional mutations after domestication relative to pre-domestication period. Correlation analysis revealed that the most recent (~ 10.000 years ago) N e (s) are positively correlated with Dω/Wω. This result is consistent with the corollary of the nearly neutral theory, that higher N e could boost the efficiency of positive selection, which might facilitate the overall accumulation of functional mutations. In addition, we also evaluated the accumulation of radical and conservative mutations during the domestication transition as: Dradical/Wradical and Dconservative/Wconservative, respectively. Surprisingly, only Dradical/Wradical ratio exhibited a positive correlation with N e (p < 0.05), suggesting that domestication process might magnify the accumulation of radical mutations in species with larger N e . CONCLUSIONS: Our results confirm the classical population genetics theory prediction and highlight the important role of species' N e in shaping the patterns of accumulation of functional mutations, especially radical mutations, in domesticated animals. The results aid our understanding of the mechanisms underlying the accumulation of functional mutations after domestication, which is critical for understanding the phenotypic diversification associated with this process.

Three new Diplozoidae mitogenomes expose unusual compositional biases within the Monogenea class: implications for phylogenetic studies.

BACKGROUND: As the topologies produced by previous molecular and morphological studies were contradictory and unstable (polytomy), evolutionary relationships within the Diplozoidae family and the Monogenea class (controversial relationships among the Discocotylinea, Microcotylinea and Gastrocotylinea suborders) remain unresolved. Complete mitogenomes carry a relatively large amount of information, sufficient to provide a much higher phylogenetic resolution than traditionally used morphological traits and/or single molecular markers. However, their implementation is hampered by the scarcity of available monogenean mitogenomes. Therefore, we sequenced and characterized mitogenomes belonging to three Diplozoidae family species, and conducted comparative genomic and phylogenomic analyses for the entire Monogenea class. RESULTS: Taxonomic identification was inconclusive, so two of the species were identified merely to the genus level. The complete mitogenomes of Sindiplozoon sp. and Eudiplozoon sp. are 14,334 bp and 15,239 bp in size, respectively. Paradiplozoon opsariichthydis (15,385 bp) is incomplete: an approximately 2000 bp-long gap within a non-coding region could not be sequenced. Each genome contains the standard 36 genes (atp8 is missing). G + T content and the degree of GC- and AT-skews of these three mitogenome (and their individual elements) were higher than in other monogeneans. nad2, atp6 and nad6 were the most variable PCGs, whereas cox1, nad1 and cytb were the most conserved. Mitochondrial phylogenomics analysis, conducted using concatenated amino acid sequences of all PCGs, indicates that evolutionary relationships of the three genera are: (Eudiplozoon, (Paradiplozoon, Sindiplozoon)); and of the three suborders: (Discocotylinea, (Microcotylinea, Gastrocotylinea)). These intergeneric relationships were also supported by the skewness and principal component analyses. CONCLUSIONS: Our results show that nad2, atp6 and nad6 (fast-evolving) would be better candidates than cox1 (slow-evolving) for species identification and population genetics studies in Diplozoidae. Nucleotide bias and codon and amino acid usage patterns of the three diplozoid mitogenomes are more similar to cestodes and trematodes than to other monogenean flatworms. This unusual mutational bias was reflected in disproportionately long branches in the phylogram. Our study offsets the scarcity of molecular data for the subclass Polyopisthocotylea to some extent, and might provide important new insights into the evolutionary history of the three genera and three suborders.

Metabolite and gene expression profiles suggest a putative mechanism through which high dietary carbohydrates reduce the content of hepatic betaine in Megalobrama amblycephala.

BACKGROUND: High-carbohydrate diets (HCD) are favoured by the aquaculture industry for economic reasons, but they can produce negative impacts on growth and induce hepatic steatosis. We hypothesised that the mechanism behind this is the reduction of hepatic betaine content. OBJECTIVE: We further explored this mechanism by supplementing betaine (1%) to the diet of a farmed fish Megalobrama amblycephala. METHODS: Four diet groups were designed: control (CD, 27.11% carbohydrates), high-carbohydrate (HCD, 36.75% carbohydrates), long-term betaine (LBD, 35.64% carbohydrates) and short-term betaine diet (SBD; 12 weeks HCD + 4 weeks LBD). We analysed growth performance, body composition, liver condition, and expression of genes and profiles of metabolites associated with betaine metabolism. RESULTS: HCD resulted in poorer growth and liver health (compared to CD), whereas LBD improved these parameters (compared to HCD). HCD induced the expression of genes associated with glucose, serine and cystathionine metabolisms, and (non-significantly, p = .20) a betaine-catabolizing enzyme betaine-homocysteine-methyltransferase; and decreased the content of betaine, methionine, S-adenosylhomocysteine and carnitine. Betaine supplementation (LBD) reversed these patterns, and elevated betaine-homocysteine-methyltransferase, S-adenosylmethionine and S-adenosylhomocysteine (all p ≤ .05). CONCLUSION: We hypothesise that HCD reduced the content of hepatic betaine by enhancing the activity of metabolic pathways from glucose to homocysteine, reflected in increased glycolysis, serine metabolism, cystathionine metabolism and homocysteine remethylation. Long-term dietary betaine supplementation improved the negative impacts of HCD, inculding growth parameters, body composition, liver condition, and betaine metabolism. However, betaine supplementation may have caused a temporary disruption in the metabolic homeostasis.

Chemotactic effect of ß-defensin 1 on macrophages in Megalobrama amblycephala.

Besides their function as a physical barrier against pathogens, ß-defensins possess the ability to induce direct or indirect chemotaxis in leukocytes of mammals. However little is known about the ability of defensins to guide the migration of macrophages in fish. The objective of our study was to investigate whether ß-defensin 1 (maBD1) can recruit leukocytes (specifically macrophages) in vivo and in vitro in a farmed cyprinid fish Megalobrama amblycephala. The M. amblycephala ß-defensin 1 (maBD1) gene was amplified from the head-kidney transcriptome. Synthetic maBD1 polypeptide (as well as its N-terminus half, but not the C-terminus half) was capable of inducing the migration of leukocytes (specifically macrophages) at concentrations from 26.0 µg/mL to 52.0 µg/mL in head kidney tissue in vitro. When injected intraperitoneally in vivo, the number of leukocytes in the peritoneal cavity was in positive correlation with the maBD1 concentration. maBD1 also induced the expression of two proinflammatory cytokines (IL-1beta and TNF-alpha) in spleen, head and body kidney, and hepatopancreas. These results strongly indicate that BD1 has a chemoattractant capacity for macrophages, as well as the ability to modulate the expression of proinflammatory cytokines in fish.

Transcriptomics, metabolomics and histology indicate that high-carbohydrate diet negatively affects the liver health of blunt snout bream (Megalobrama amblycephala).

BACKGROUND: Global trend of the introduction of high levels of relatively cheap carbohydrates to reduce the amount of costly protein in the aquatic animal feed production has affected the aquaculture of an economically important cyprinid fish, blunt snout bream (Megalobrama amblycephala). This dietary shift has resulted in increased prevalence of metabolic disorders, often causing economic losses. High dietary intake of carbohydrates, associated with obesity, is one of the major causes of non-alcoholic fatty liver disease (NAFLD) in humans. RESULTS: We have conducted an eight-week feeding trial to better understand how a high-carbohydrate diet (HCBD) affects the liver health in this fish. Hepatosomatic index and lipid content were significantly (P < 0.05) higher in the HCBD group. Histology results also suggested pathological changes in the livers of HCBD group, with excessive lipid accumulation and indication of liver damage. Metabolomics and serum biochemistry analyses showed that a number of metabolites indicative of liver damage were increased in the HCBD group. This group also exhibited low levels of betaine, which is a metabolite crucial for maintaining the healthy liver functions. Transcriptomic and qPCR analyses indicated that HCBD had a strong impact on the expression of a large number of genes associated with the NAFLD and insulin signalling pathways, which may lead to the development of insulin resistance in hepatocytes, pathological liver changes, and eventually the NAFLD. CONCLUSIONS: Transcriptomics, metabolomics and histology results all indicate early symptoms of liver damage. However whether these would actually lead to the development of NAFLD after a longer period of time, remains inconclusive. Additionally, a very high number of upregulated genes in the HCBD group associated with several neurodegenerative diseases is a strong indication of neurodegenerative changes caused by the high-carbohydrate diet in blunt snout bream. This suggests that fish might present a good model to study neurodegenerative changes associated with high-carbohydrate diet in humans.

The complete mitochondrial genome of parasitic nematode Camallanus cotti: extreme discontinuity in the rate of mitogenomic architecture evolution within the Chromadorea class.

BACKGROUND: Complete mitochondrial genomes are much better suited for the taxonomic identification and phylogenetic studies of nematodes than morphology or traditionally-used molecular markers, but they remain unavailable for the entire Camallanidae family (Chromadorea). As the only published mitogenome in the Camallanina suborder (Dracunculoidea superfamily) exhibited a unique gene order, the other objective of this research was to study the evolution of mitochondrial architecture in the Spirurida order. Thus, we sequenced the complete mitogenome of the Camallanus cotti fish parasite and conducted structural and phylogenomic comparative analyses with all available Spirurida mitogenomes. RESULTS: The mitogenome is exceptionally large (17,901 bp) among the Chromadorea and, with 46 (pseudo-) genes, exhibits a unique architecture among nematodes. Six protein-coding genes (PCGs) and six tRNAs are duplicated. An additional (seventh) tRNA (Trp) was probably duplicated by the remolding of tRNA-Ser2 (missing). Two pairs of these duplicated PCGs might be functional; three were incomplete and one contained stop codons. Apart from Ala and Asp, all other duplicated tRNAs are conserved and probably functional. Only 19 unique tRNAs were found. Phylogenomic analysis included Gnathostomatidae (Spirurina) in the Camallanina suborder. CONCLUSIONS: Within the Nematoda, comparable PCG duplications were observed only in the enoplean Mermithidae family, but those result from mitochondrial recombination, whereas characteristics of the studied mitogenome suggest that likely rearrangement mechanisms are either a series of duplications, transpositions and random loss events, or duplication, fragmentation and subsequent reassembly of the mitogenome. We put forward a hypothesis that the evolution of mitogenomic architecture is extremely discontinuous, and that once a long period of stasis in gene order and content has been punctuated by a rearrangement event, such a destabilised mitogenome is much more likely to undergo subsequent rearrangement events, resulting in an exponentially accelerated evolutionary rate of mitogenomic rearrangements. Implications of this model are particularly important for the application of gene order similarity as an additive source of phylogenetic information. Chromadorean nematodes, and particularly Camallanina clade (with C. cotti as an example of extremely accelerated rate of rearrangements), might be a good model to further study this discontinuity in the dynamics of mitogenomic evolution.

Molecular identification and functional characterisation of the interferon regulatory factor 1 in the blunt snout bream (Megalobrama amblycephala).

Interferon regulatory factors (IRFs) play a key role in mediating the host response against pathogen infection and other important biological processes. This is the first report of an IRF family member in blunt snout bream Megalobrama amblycephala. The complete cDNA of M. amblycephala (Ma) IRF1 gene has 1422 nucleotides (nt.), with an open reading frame of 858 nt, encoding a polypeptide of 285 amino acids. The putative MaIRF1 polypeptide shared significant structural homology with known IRF1 homologs: a conserved IRF domain was found at the N-terminal and an IRF association domain 2 at the C-terminal. Phylogenetic analysis showed that MaIRF1 amino acid sequence clustered with other teleost IRF1s, with a grass carp ortholog exhibiting the highest similarity. MaIRF1 mRNA expression patterns were studied using quantitative real-time PCR in healthy fish tissues and after a challenge with Aeromonas hydrophila bacterium. It was constitutively expressed in all examined tissues: the highest in blood, the lowest in muscle. The expression after A. hydrophila challenge was up-regulated in liver, spleen and kidney, but down-regulated in intestine and gills. At the protein level, similar expression patterns were observed in liver and gills. Patterns differed in intestine (up-regulation), spleen (down-regulation) and kidney (expression mostly unchanged). This study indicates that MaIRF1 gene plays an important role in the blunt snout bream immune system, hence providing an important base for further studies.

Blunt Snout Bream (Megalobrama amblycephala) MyD88 and TRAF6: characterisation, comparative homology modelling and expression.

MyD88 and TRAF6 play an essential role in the innate immune response in most animals. This study reports the full-length MaMyD88 and MaTRAF6 genes identified from the blunt snout bream (Megalobrama amblycephala) transcriptome profile. MaMyD88 is 2501 base pairs (bp) long, encoding a putative protein of 284 amino acids (aa), including the N-terminal DEATH domain of 78 aa and the C-terminal TIR domain of 138 aa. MaTRAF6 is 5474 bp long, encoding a putative protein of 542 aa, including the N-terminal low-complexity region, RING domain (40 aa), a coiled-coil region (64 aa) and C-terminal MATH domain (147 aa). Coding regions of MaMyD88 and MaTRAF6 genomic sequences consisted of five and six exons, respectively. Physicochemical and functional characteristics of the proteins were analysed. Alpha helices were dominant in the secondary structure of the proteins. Homology models of the MaMyD88 and MaTRAF6 domains were constructed applying the comparative modelling method. RT-qPCR was used to analyse the expression of MaMyD88 and MaTRAF6 mRNA transcripts in response to Aeromonas hydrophila challenge. Both genes were highly upregulated in the liver, spleen and kidney during the first 24 h after the challenge. While MyD88 and TRAF6 have been reported in various aquatic species, this is the first report and characterisation of these genes in blunt snout bream. This research also provides evidence of the important roles of these two genes in the blunt snout bream innate immune system.

Geography, phylogeny and host switch drive the coevolution of parasitic Gyrodactylus flatworms and their hosts.

BACKGROUND: Gyrodactylus is a lineage of monogenean flatworm ectoparasites exhibiting many features that make them a suitable model to study the host-parasite coevolutionary dynamics. Previous coevolutionary studies of this lineage mainly relied on low-power datasets (a small number of samples and a single molecular marker) and (now) outdated algorithms. METHODS: To investigate the coevolutionary relationship of gyrodactylids and their fish hosts in high resolution, we used complete mitogenomes (including two newly sequenced Gyrodactylus species), a large number of species in the single-gene dataset, and four different coevolutionary algorithms. RESULTS: The overall coevolutionary fit between the parasites and hosts was consistently significant. Multiple indicators confirmed that gyrodactylids are generally highly host-specific parasites, but several species could parasitize either multiple (more than 5) or phylogenetically distant fish hosts. The molecular dating results indicated that gyrodactylids tend to evolve towards high host specificity. Speciation by host switch was identified as a more important speciation mode than co-speciation. Assuming that the ancestral host belonged to Cypriniformes, we inferred four major host switch events to non-Cypriniformes hosts (mostly Salmoniformes), all of which occurred deep in the evolutionary history. Despite their relative rarity, these events had strong macroevolutionary consequences for gyrodactylid diversity. For example, in our dataset, 57.28% of all studied gyrodactylids parasitized only non-Cypriniformes hosts, which implies that the evolutionary history of more than half of all included lineages could be traced back to these major host switch events. The geographical co-occurrence of fishes and gyrodactylids determined the host use by these gyrodactylids, and geography accounted for most of the phylogenetic signal in host use. CONCLUSIONS: Our findings suggest that the coevolution of Gyrodactylus flatworms and their hosts is largely driven by geography, phylogeny, and host switches.