Reduced genetic diversity and the success of the invasive peacock bass (Cichliformes: Cichlidae) / Diversidade genética reduzida e o sucesso da espécie invasora tucunaré (Cichliformes: Cichlidae)

Abstract Several species of Cichla successfully colonized lakes and reservoirs of Brazil, since the 1960's, causing serious damage to local wildlife. In this study, 135 peacock bass were collected in a reservoir complex in order to identify if they represented a single dominant species or multiple ones, as several Cichla species have been reported in the basin. Specimens were identified by color pattern, morphometric and meristic data, and using mitochondrial markers COI, 16S rDNA and Control Region (CR). Overlapping morphological data and similar coloration patterns prevented their identification using the taxonomic keys to species identification available in the literature. However, Bayesian and maximum likelihood from sequencing data demonstrated the occurrence of a single species, Cichla kelberi. A single haplotype was observed for the 16S and CR, while three were detected for COI, with a dominant haplotype present in 98.5% of the samples. The extreme low diversity of the transplanted C. kelberi evidenced a limited number of founding maternal lineages. The success of this colonization seems to rely mainly on abiotic factors, such as increased water transparency of lentic environments that favor visual predators that along with the absence of predators, have made C. kelberi a successful invader of these reservoirs.

Resumo Muitas espécies de Cichla colonizaram com sucesso lagos e reservatórios do Brasil desde os anos 1960, causando graves prejuízos à vida selvagem nesses locais. Neste estudo, 135 tucunarés foram coletados em um complexo de reservatórios a fim de identificar se representavam uma espécie dominante ou múltiplas espécies, uma vez que diversas espécies de Cichla foram registradas na bacia. Os espécimes foram identificados com base na coloração, dados morfométricos e merísticos, e por marcadores mitocondriais COI, 16S rDNA e Região Controle (RC). A sobreposição dos dados morfométricos e o padrão similar de coloração impediram a identificação utilizando as chaves de identificação disponíveis na literatura. Entretanto, as análises bayesiana e de máxima verossimilhança de dados moleculares demonstraram a ocorrência de uma única espécie, Cichla kelberi. Um único haplótipo foi observado para o 16S e RC, enquanto três foram detectados para o COI, com um haplótipo dominante presente em 98,5% das amostras. A baixa diversidade nos exemplares introduzidos de C. kelberi evidenciou um número limitado de linhagens maternas fundadoras. O sucesso da invasão parece depender de fatores abióticos, como a maior transparência da água de ambientes lênticos que favorece predadores visuais que, atrelado à ausência de predadores, fez do C. kelberi um invasor bem-sucedido nesses reservatórios.

Physiological and immune profiling of tilapia Oreochromis niloticus gills by high-throughput single-cell transcriptome sequencing.

The physiological and immune functions of fish gills are largely recognized, but their following functional heterogeneity at the single cell scale has been rarely reported. Here, we performed single cell RNA sequencing (scRNA-seq) on the gills of tilapia fish Oreochromis niloticus. We identified a total of 12 cell populations and analyzed their functional heterogeneity. To investigate the physiological function of O. niloticus gills, expression patterns of genes encoding ion transporters were selected from the identified H+-ATPase-rich cells (HR cells), Na+/K+-ATPase-rich cells (NaR cells), and pavement cells. Specific enrichment of ca4a, slc9a1a, and LOC100692482 in the HR cells of O. niloticus gills explained their functions in acid-base regulation. Genes encoding Ca2+ transporters, including atp2b1, LOC100696627, and LOC 100706765, were specifically expressed in the NaR cells. Pavement cells were presumably the main sites responsible for ammonia and urea transports in O. niloticus gills with specific enrichment of Rhbg and LOC100693008, respectively. The expression patterns of the four immune cell subtypes varied greatly, with B cells being enriched with the most immune-related GO terms. KEGG enrichment analysis showed that MAPK signaling pathway was the most enriched pathway among the four types of immune cells in O. niloticus gills. Our results are important in understanding the physiological and immune responses of fish gills at the cellular resolution.

Transcriptome analysis reveals the neuroactive receptor genes response to Streptococcus agalactiae infection in tilapia, Oreochromis niloticus.

The detailed crosstalk between the neuroendocrine and immune systems in Oreochromis niloticus, an economically important fish, in response to pathogenic infections, remains unclear. This study revealed the head kidney transcriptional profiles of O. niloticus upon infections with Streptococcus agalactiae, a prevalent pathogen known to cause severe meningitis. Twelve cDNA libraries of O. niloticus head kidney, representing four treatment time points (0, 6, 24, and 48 h), were constructed and a total of 2,528 differentially expressed genes were identified based on pairwise comparisons. KEGG pathway analysis revealed a significant enrichment of the 'neuroactive ligand-receptor interaction' pathway (ko04080), with 13 genes exhibiting differential expression during S. agalactiae infection. Among these, six neuroactive receptor genes (lepr, nr3c1, ptger4, thrb, tspo, and ß2-ar) were selected, cloned, and characterized. Although these genes are ubiquitously expressed, and in head kidney leukocytes, their expression was mainly observed in T cells, Mo/Mφ, and NCCs, which are characterized by antimicrobial responses. Furthermore, we examined the response patterns of these six neuroactive receptor genes to gram-positive (S. agalactiae) and gram-negative (Aeromonas hydrophila) bacteria in four different tissues. Notably, lepr, ptger4, tspo, and ß2-ar were upregulated in all selected tissues in response to S. agalactiae and A. hydrophila infections. However, nr3c1 and thrb were downregulated in response to S. agalactiae infection in the head kidney and spleen, whereas nr3c1 was upregulated, and thrb was unresponsive to A. hydrophila infection. Our findings provide a theoretical foundation for understanding new links between the neuroendocrine and immune systems during bacterial infection in teleost fish.

Visual opsin gene expression evolution in the adaptive radiation of cichlid fishes of Lake Tanganyika.

Tuning the visual sensory system to the ambient light is essential for survival in many animal species. This is often achieved through duplication, functional diversification, and/or differential expression of visual opsin genes. Here, we examined 753 new retinal transcriptomes from 112 species of cichlid fishes from Lake Tanganyika to unravel adaptive changes in gene expression at the macro-evolutionary and ecosystem level of one of the largest vertebrate adaptive radiations. We found that, across the radiation, all seven cone opsins-but not the rhodopsin-rank among the most differentially expressed genes in the retina, together with other vision-, circadian rhythm-, and hemoglobin-related genes. We propose two visual palettes characteristic of very shallow- and deep-water living species, respectively, and show that visual system adaptations along two major ecological axes, macro-habitat and diet, occur primarily via gene expression variation in a subset of cone opsin genes.

Two granzyme A/K homologs in Zebra mbuna have different specificities, one classical tryptase and one with chymase activity.

Granzymes A and K are two highly homologous serine proteases expressed by mammalian cytotoxic T cells (CTLs) and natural killer (NK) cells. The locus encoding these two proteases is the first of the hematopoietic serine protease loci to appear during vertebrate evolution. This locus is found in all jawed vertebrates including the cartilaginous fishes. Granzyme A is the most abundant of the different granzymes expressed by CTLs and NK cells and its potential function has been studied extensively for many years. However, no clear conclusions concerning its primary role in the immune defense has been obtained. In all mammals, there are only one copy each of granzyme A and K, whereas additional copies are found in both cartilaginous and ray finned fishes. In cichlids two of these copies seem to encode new members of the granzyme A/K family. These two new members appear to have changed primary specificity and to be pure chymases based on the amino acids in their active site substrate binding pockets. Interestingly, one of these gene copies is located in the middle of the granzyme A/K locus, while the other copy is present in another locus, the met-ase locus. We here present a detailed characterization of the extended cleavage specificity of one of these non-classical granzymes, a Zebra mbuna granzyme positioned in the granzyme A/K locus. This enzyme, named granzyme A2, showed a high preference for tyrosine in the P1 position of substrates, thereby being a strict chymase. We have also characterized one of the classical granzyme A/Ks of the Zebra mbuna, granzyme A1, which is a tryptase with preference for arginine in the P1 position of substrates. Based on their extended specificities, the two granzymes showed major similarities, but also some differences in preferred amino acids in positions surrounding the cleavable amino acid. Fish lack one of the hematopoietic serine protease loci of mammals, the chymase locus, where one of the major mast cell enzymes is located. An interesting question is now if cichlids have by compensatory mechanisms generated a mast cell chymase from another locus, and if similar chymotryptic enzymes have appeared also in other fish species.

Mitogenomic Characterization of Cameroonian Endemic Coptodon camerunensis (Cichliformes: Cichlidae) and Matrilineal Phylogeny of Old-World Cichlids.

The mitogenomic evolution of old-world cichlids is still largely incomplete in Western Africa. In this present study, the complete mitogenome of the Cameroon endemic cichlid, Coptodon camerunensis, was determined by next-generation sequencing. The mitogenome was 16,557 bp long and encoded with 37 genes (13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, and a control region). The C. camerunensis mitogenome is AT-biased (52.63%), as exhibited in its congener, Coptodon zillii (52.76% and 53.04%). The majority of PCGs start with an ATG initiation codon, except COI, which starts with a GTG codon and five PCGs and ends with the TAA termination codon and except seven PCGs with an incomplete termination codon. In C. camerunensis mitogenome, most tRNAs showed classical cloverleaf secondary structures, except tRNA-serine with a lack of DHU stem. Comparative analyses of the conserved blocks of two Coptodonini species control regions revealed that the CSB-II block was longer than other blocks and contained highly variable sites. Using 13 concatenated PCGs, the mitogenome-based Bayesian phylogeny easily distinguished all the examined old-world cichlids. Except for Oreochromini and Coptodinini tribe members, the majority of the taxa exhibited monophyletic clustering within their respective lineages. C. camerunensis clustered closely with Heterotilapia buttikoferi (tribe Heterotilapiini) and had paraphyletic clustering with its congener, C. zillii. The Oreochromini species also displayed paraphyletic grouping, and the genus Oreochromis showed a close relationship with Coptodinini and Heterotilapiini species. In addition, illustrating the known distribution patterns of old-world cichlids, the present study is congruent with the previous hypothesis and proclaims that prehistoric geological evolution plays a key role in the hydroclimate of the African continent during Mesozoic, which simultaneously disperses and/or colonizes cichlids in different ichthyological provinces and Rift Lake systems in Africa. The present study suggests that further mitogenomes of cichlid species are required, especially from western Africa, to understand their unique evolution and adaptation.

The role of antifreeze genes in the tolerance of cold stress in the Nile tilapia (Oreochromis niloticus).

BACKGROUND: Tilapia is one of the most essential farmed fishes in the world. It is a tropical and subtropical freshwater fish well adapted to warm water but sensitive to cold weather. Extreme cold weather could cause severe stress and mass mortalities in tilapia. The present study was carried out to investigate the effects of cold stress on the up-regulation of antifreeze protein (AFP) genes in Nile tilapia (Oreochromis niloticus). Two treatment groups of fish were investigated (5 replicates of 15 fish for each group in fibreglass tanks/70 L each): 1) a control group; the fish were acclimated to lab conditions for two weeks and the water temperature was maintained at 25 °C during the whole experimental period with feeding on a commercial diet (30% crude protein). 2) Cold stress group; the same conditions as the control group except for the temperature. Initially, the temperature was decreased by one degree every 12 h. The fish started showing death symptoms when the water temperature reached 6-8 °C. In this stage the tissue (muscle) samples were taken from both groups. The immune response of fish exposed to cold stress was detected and characterized using Differential Display-PCR (DD-PCR). RESULTS: The results indicated that nine different up-regulation genes were detected in the cold-stressed fish compared to the control group. These genes are Integrin-alpha-2 (ITGA-2), Gap junction gamma-1 protein-like (GJC1), WD repeat-containing protein 59 isoform X2 (WDRP59), NUAK family SNF1-like kinase, G-protein coupled receptor-176 (GPR-176), Actin cytoskeleton-regulatory complex protein pan1-like (PAN-1), Whirlin protein (WHRN), Suppressor of tumorigenicity 7 protein isoform X2 (ST7P) and ATP-binding cassette sub-family A member 1-like isoform X2 (ABCA1). The antifreeze gene type-II amplification using a specific PCR product of 600 bp, followed by cloning and sequencing analysis revealed that the identified gene is antifreeze type-II, with similarity ranging from 70 to 95%. The in-vitro transcribed gene induced an antifreeze protein with a molecular size of 22 kDa. The antifreeze gene, ITGA-2 and the WD repeat protein belong to the lectin family (sugar-protein). CONCLUSIONS: In conclusion, under cold stress, Nile tilapia express many defence genes, an antifreeze gene consisting of one open reading frame of approximately 0.6 kbp.

The effects of 100 nm pure boron nanoparticles on the primary antioxidant system of Nile tilapia (Oreochromisniloticus).

In this study, the effects of 100 nm boron nanoparticles (B-NPs) on the primary antioxidant status of Nile tilapia were researched via analysis of enzyme activities and related gene expressions. This is a new study which focuses on the relationship between B-NPs and oxidative stress that contribute to the literature in terms of its scope. Fish (n = 15) for each group were exposed to three different concentrations as 5, 25 (n2) and 125 (n3) mg/L during 96 h to see the response of the primary antioxidant system. According to the results, SOD expressions differed in all treatment groups compared to the control group (P < 0.05). CAT expressions were different in 5 and 125 mg/L groups compared to control and 25 mg/L groups (P < 0.05). GPX expressions were only different in 125 mg/L group (P < 0.05). The changes in enzyme activities of SOD and CAT were significantly different in 25 mg/L groups. GPX enzyme activities were not significant (P > 0.05). TBARS concentrations in 25 mg/L group were significantly different from those in the control and 125 mg/L groups (P < 0.05).

Cellular profiling of a recently-evolved social behavior in cichlid fishes.

Social behaviors are diverse in nature, but it is unclear how conserved genes, brain regions, and cell populations generate this diversity. Here we investigate bower-building, a recently-evolved social behavior in cichlid fishes. We use single nucleus RNA-sequencing in 38 individuals to show signatures of recent behavior in specific neuronal populations, and building-associated rebalancing of neuronal proportions in the putative homolog of the hippocampal formation. Using comparative genomics across 27 species, we trace bower-associated genome evolution to a subpopulation of glia lining the dorsal telencephalon. We show evidence that building-associated neural activity and a departure from quiescence in this glial subpopulation together regulate hippocampal-like neuronal rebalancing. Our work links behavior-associated genomic variation to specific brain cell types and their functions, and suggests a social behavior has evolved through changes in glia.

Interleukin-12 induces IFN-γ secretion and STAT signaling implying its potential regulation of Th1 cell response in Nile tilapia.

As a pleiotropic cytokine consisting of IL-12p35 and IL-12p40, Interleukin-12 (IL-12) features in inflammation regulation and anti-bacterial immunity. While IL-12 homologs have been identified in non-mammalian species, the precise mechanisms by which IL-12 contributes to early adaptive immune responses in vertebrates remain incompletely understood. Herein, an evolutionary conserved Oreochromis niloticus IL-12 (defined as OnIL-12) was identified by synteny characterization, structural comparisons and phylogenetic pattern of IL-12p35b and IL-12p40a. IL-12p35b and IL-12p40a exhibited widespread expression in lymphoid-related tissues of tilapia, while their mRNA expression in head-kidney demonstrated a significant increase after Edwardsiella piscicida infection. Compared with other lymphocytes, recombinant OnIL-12 (rOnIL-12) displayed stronger affinity binding to T cells. Although stimulation of lymphocytes with the p35b or p40a subunit resulted in a significant induction of IFN-γ expression, rOnIL-12 showed stronger potential to promote IFN-γ expression than these subunits. rOnIL-12 not only elevated the mRNA expression level Th1 cell-associated transcription factor T-bet in lymphocytes, but also increased the proportion of CD4-1+IFN-γ+ lymphocytes. Moreover, the mRNA and phosphorylation levels of STAT1, STAT3, STAT4 and STAT5 were enhanced by rOnIL-12. These findings will offer previous evidence for further exploration into the regulatory mechanisms of Th1 cellular immunity in early vertebrates.

Ancient standing genetic variation facilitated the adaptive radiation of Lake Victoria cichlids.

Cichlid fishes are textbook examples of explosive speciation and adaptive radiation, providing a great opportunity to understand how the genomic substrate yields extraordinary species diversity. Recently, we performed comparative genomic analyses of three Lake Victoria cichlids to reveal the genomic substrates underlying their rapid speciation and adaptation. We found that long divergent haplotypes derived from large-scale standing genetic variation, which originated before the adaptive radiation of Lake Victoria cichlids, may have contributed to their rapid diversification. In addition, the present study on genomic data from other East African cichlids suggested the reuse of alleles that may have originated in the ancestral lineages of Lake Tanganyika cichlids during cichlid evolution. Therefore, our results highlight that the primary factor that could drive repeated adaptive radiation across East African cichlids was allelic reuse from standing genetic variation to adapt to their own specific environment. In this report, we summarize the main results and discuss the evolutionary mechanisms of cichlids, based on our latest findings.

Population structure and genetic diversity of Nile tilapia (Oreochromis niloticus) using microsatellite markers from selected water bodies in southwest Ethiopia.

BACKGROUND: The Ethiopian southwest drainage has been designated a huge potential for fishery and aquaculture development. However, the genetic diversity of Nile tilapia along the main water bodies has yet to be uncovered. Therefore, this study used microsatellite markers to assess the genetic structure and diversity of Nile cichlid (Oreochromis niloticus) populations. METHODS: One hundred and sixty-seven fish were collected from Alwero Dam, Baro River, Gilgel Gibe I Dam and Gojeb River southwest of Ethiopia. Genomic DNA was extracted from 40 tissue samples using a Qiagen DNA tissue extraction kit and genetic polymorphism was estimated using GenAlEx6.502 software. RESULTS: The result showed that four loci were polymorphic and produced 3-6 alleles per locus, while the number of effective alleles (Ne) and the number of distinct alleles (Na) were 1.370 and 2.063, respectively. Observed heterozygosity (Ho) and expected heterozygosity (He) were expanded between 0.100.50 and 0.115-0.265, respectively. Baro River population revealed the highest mean expected heterozygosity of 0.265, whereas the lowest expected mean heterozygosity of 0.115 was from the Alwero Dam. The genetic differentiation between populations (FST) resulted in a moderate to a very high degree of divergence between the populations. The Alwero Dam and Gojeb River populations were the most divergent populations (FST = 0.524), unlike those from the Alwero Dam and Baro River. The genetic variation among individuals within populations was 69%, whereas 31% appeared among populations and no variation within individuals to the overall variance. CONCLUSION: This study confirmed the presence of significant genetic variations and differentiation among populations except between Gilgel Gibe I Dam and Gojeb River. This study highlights the genetic structure and differentiation amongst wild Nile tilapia (O. niloticus) populations in south-western Ethiopia and paves the way for selective aquaculture breeding programs in the study area.

Comparative analysis of quantitative phosphoproteomics between two tilapias (Oreochromis niloticus and Oreochromis aureus) under low-temperature stress.

As an important farmed fish, tilapia has poor tolerance to low-temperatures. At the same time, different tilapia strains have apparent differences in low-temperature tolerance. In this study, using the iTRAQ method, the phosphorylated proteomics of two tilapia strains (Oreochromis niloticus and Oreochromis aureus) with different tolerances to low-temperature stress were quantitatively and comparatively analyzed, to clarify the physiological mechanism of tilapia's response to low-temperature stress. Through the GO and IPR analyses of differentially phosphorylated proteins, a number of similarities in physiological activities and regulatory effects were found between the two tilapias in response to low-temperature stress. Many differentially phosphorylated proteins are mainly involved in lipid metabolism, cell proliferation and apoptosis. However, the difference in endurance of low temperature of these two tilapias might be related to the differences in categories, expression and modification level of genetic products which were involved in the aforementioned physiological processes. And meanwhile, the enrichment results of KEGG showed the changes of multiple immune-related and growth-related phosphorylated proteins in the cytokine-cytokine receptor interaction pathway in O. aureus are more prominent. Furthermore, the significantly enriched pathway of carbohydrate digestion and absorption in O. niloticus may indicate that low-temperature stress exerts a more severe impact on energy metabolism. The relative results would help elucidating the molecular mechanism by which tilapia responds to low-temperature stress, and developing culture of tilapia species.

Evolution of MHC IIB Diversity Across Cichlid Fish Radiations.

The genes of the major histocompatibility complex (MHC) are among the most polymorphic genes in vertebrates and crucial for their adaptive immune response. These genes frequently show inconsistencies between allelic genealogies and species phylogenies. This phenomenon is thought to be the result of parasite-mediated balancing selection maintaining ancient alleles through speciation events (trans-species polymorphism [TSP]). However, allele similarities may also arise from postspeciation mechanisms, such as convergence or introgression. Here, we investigated the evolution of MHC class IIB diversity in the cichlid fish radiations across Africa and the Neotropics by a comprehensive review of available MHC IIB DNA sequence information. We explored what mechanism explains the MHC allele similarities found among cichlid radiations. Our results showed extensive allele similarity among cichlid fish across continents, likely due to TSP. Functionality at MHC was also shared among species of the different continents. The maintenance of MHC alleles for long evolutionary times and their shared functionality may imply that certain MHC variants are essential in immune adaptation, even in species that diverged millions of years ago and occupy different environments.

The cichlid pharyngeal jaw novelty enhances evolutionary integration in the feeding apparatus.

The modified pharyngeal jaw system of cichlid fishes is widely viewed as a key innovation that substantially facilitated the evolutionary exuberance of this iconic evolutionary radiation. We conduct comparative phylogenetic analyses of integration, disparity, and rate of evolution among feeding-related, skeletal structures in Neotropical cichlids and North American centrarchids, which lack the specialized pharyngeal jaw. Contrasting evolutionary patterns in these two continental radiations, we test a classic decoupling hypothesis. Specifically, we ask whether the modified pharyngeal jaw in cichlids resulted in enhanced evolutionary independence of the oral and pharyngeal jaws, leading to increased diversity of trophic structures. Contrary to this prediction, we find significantly stronger evolutionary integration between the oral and pharyngeal jaws in cichlids compared to centrarchids, although the two groups do not differ in patterns of integration within each jaw system. Further, though we find no significant differences in disparity, centrarchids show faster rates of morphological evolution. Our results suggest that the modified pharyngeal jaw resulted in less evolutionary independence and slower rates of evolution within the feeding system. Thus, we raise the possibility that the cichlid novelty enhances feeding performance, but does not prompt increased morphological diversification within the feeding apparatus, as has long been thought.

DHA alleviates high glucose-induced mitochondrial dysfunction in Oreochromis niloticus by inhibiting DRP1-mediated mitochondrial fission.

Dynamin-related protein 1 (DRP1) is a key regulator in the maintenance of mammalian glucose homeostasis, but the relevant information remains poorly understood on aquatic animals. In the study, DRP1 is formally described for the first time in Oreochromis niloticus. DRP1 encodes a peptide of 673 amino acid residues that contained three conserved domains: a GTPase domain, a dynamin middle domain and a dynamin GTPase effector domain. DRP1 transcripts are widely distributed in all of the detected seven organs/tissues, and the highest mRNA levels in brain. High-carbohydrate (45 %) fed fish showed a significant upregulation of liver DRP1 expression than that of control (30 %) group. Glucose administration upregulated liver DRP1 expression, with peak values observed at 1 h; then its expression returned to the basal value at 12 h. In the in vitro study, DRP1 over-expression significantly decreased mitochondrial abundance in hepatocytes. DHA significantly increased mitochondrial abundance, transcriptions of mitochondrial transcription factor A (TFAM) and mitofusin 1 and 2 (MFN1 and MFN2) and complex II and III activities of high glucose-treated hepatocyte, whereas the opposite was true for DRP1, mitochondrial fission factor (MFF) and fission (FIS) expression. Together, these findings illustrated that O. niloticus DRP1 is highly conserved, and it participated in glucose control of fish. DHA could alleviate high glucose-induced mitochondrial dysfunction of fish by inhibiting DRP1-mediated mitochondrial fission.

Comparison of Toxicity of Nano and Bulk Titanium Dioxide on Nile Tilapia (Oreochromis niloticus): Acetylcholinesterase Activity Modulation and DNA Damage.

This study compared effects of low concentrations (0.05 and 0.1 mg/L) of nano-TiO2 and bulk-TiO2 on brain, gill and liver acetylcholinesterase (AChE) and erythrocytic DNA of Nile tilapia over 7 and 14 days exposure. Both TiO2 forms did not affect brain AChE activities. Bulk-TiO2 induced elevation of gill AChE activities only after 7 days while nano-TiO2 had no effect. Liver AChE activities were increased by 0.1 mg/L bulk- and nano-TiO2 to similar extents. At 7 days, erythrocytic DNA damage was induced only by 0.1 mg/L nano- and bulk-TiO2 to similar extents, but damage was not repaired to control levels at 7 days recovery period. At 14 days continuous exposure, DNA damage was induced by 0.05, 0.1 mg/L nano-TiO2 and 0.1 mg/L of bulk-TiO2 to similar extents. Results show that both forms of TiO2 can pose genotoxic hazards to fish populations under sub-chronic exposure. However, their neurotoxic potential was not evident.

The transcriptional state and chromatin landscape of cichlid jaw shape variation across species and environments.

Adaptive phenotypes are shaped by a combination of genetic and environmental forces, but how they interact remains poorly understood. Here, we utilize the cichlid oral jaw apparatus to better understand these gene-by-environment effects. First, we employed RNA-seq in bony and ligamentous tissues important for jaw opening to identify differentially expressed genes between species and across foraging environments. We used two Lake Malawi species adapted to different foraging habitats along the pelagic-benthic ecomorphological axis. Our foraging treatments were designed to force animals to employ either suction or biting/scraping, which broadly mimic pelagic or benthic modes of feeding. We found a large number of differentially expressed genes between species, and while we identified relatively few differences between environments, species differences were far more pronounced when they were challenged with a pelagic versus benthic foraging mode. Expression data carried the signature of genetic assimilation, and implicated cell cycle regulation in shaping the jaw across species and environments. Next, we repeated the foraging experiment and performed ATAC-seq procedures on nuclei harvested from the same tissues. Cross-referencing results from both analyses revealed subsets of genes that were both differentially expressed and differentially accessible. This reduced dataset implicated notable candidate genes including the Hedgehog effector, KIAA0586 and the ETS transcription factor, etv4, which connects environmental stress and craniofacial morphogenesis. Taken together, these data provide novel insights into the epigenetic, genetic and cellular bases of species- and environment-specific bone shapes.

Impact of genetic relatedness on reproductive behavior in Pelvicachromis pulcher, a biparental cichlid fish with mutual mate choice and ornamentation.

Inbreeding can result in inbreeding depression. Therefore, many species seek to avoid inbreeding. However, theory predicts that inbreeding can be beneficial. Accordingly, some species tolerate inbreeding or even prefer mating with close relatives. Evidence for active inbreeding, i.e., kin-mating preference was reported in the biparental African cichlid fish Pelvicachromis taeniatus. Related mating partners revealed better parental cooperation due to kin selection, a potential benefit of inbreeding. In this study, we investigated kin-mating preference in a genetically diverse, outbred F2-lab population of Pelvicachromis pulcher, a closely related species to P. taeniatus. Like P. taeniatus, this species shows mutual ornamentation and mate choice as well as intense biparental brood care. The F1 P. pulcher generation had revealed signs of inbreeding depression but no inbreeding avoidance. We studied mating behavior and aggression in trios consisting of a male P. pulcher, an unfamiliar sister, and an unfamiliar, unrelated female. Because the study focused on kin-mating patterns, female pairs were matched for body size and coloration. The results provide no evidence for inbreeding avoidance but rather suggest inbreeding preference. We also found no significant impact of inbreeding on offspring survival. The results suggest no inbreeding avoidance in P. pulcher; however, the strength of inbreeding preference and inbreeding depression seems to be variable. We discuss possible causes for this variation like context-dependent inbreeding depression. The number of eggs positively correlated with female body size and coloration. Furthermore, the female aggressiveness was positively correlated with female coloration indicating that coloration signal female dominance and quality.

Variation in parasite infection between replicates of speciation in Lake Victoria cichlid fish.

Because of potentially strong eco-evolutionary interactions with their hosts, parasites may initiate or enhance host diversification. The adaptive radiation of cichlid fish in Lake Victoria provides a good system to study the role of parasites at different stages of host speciation. We analyzed the macroparasite infection of four replicates of sympatric blue and red Pundamilia species pairs that vary in their age and extent of differentiation. Sympatric host species differed in parasite community composition and in the infection levels of some of these parasite taxa. Most infection differences were consistent between sampling years, indicating temporal consistency in parasite-mediated divergent selection between species. Infection differentiation increased linearly with genetic differentiation. However, significant infection differences between sympatric species were only found in the oldest, most strongly differentiated Pundamilia species pair. This is inconsistent with parasite-driven speciation. Next, we identified five distinct species of Cichlidogyrus, a genus of highly specific gill parasites that has radiated elsewhere in Africa. Infection profiles of species of Cichlidogyrus differed between sympatric cichlid species only in the oldest and most differentiated pair, again inconsistent with parasite-mediated speciation. To conclude, parasites may contribute to host differentiation after speciation, but do not initiate host speciation.