Do saline taxa evolve faster? Comparing relative rates of molecular evolution between freshwater and marine eukaryotes

Resolving Marine-Freshwater Transitions by Diatoms Through a Fog of Gene Tree Discordance

Despite the obstacles facing marine colonists, most lineages of aquatic organisms have colonized and diversified in freshwaters repeatedly. These transitions can trigger rapid morphological or physiological change and, on longer timescales, lead to increased rates of speciation and extinction. Diatoms are a lineage of ancestrally marine microalgae that have diversified throughout freshwater habitats worldwide. We generated a phylogenomic data set of genomes and transcriptomes for 59 diatom taxa to resolve freshwater transitions in one lineage, the Thalassiosirales. Although most parts of the species tree were consistently resolved with strong support, we had difficulties resolving a Paleocene radiation, which affected the placement of one freshwater lineage. This and other parts of the tree were characterized by high levels of gene tree discordance caused by incomplete lineage sorting and low phylogenetic signal. Despite differences in species trees inferred from concatenation versus summary methods and codons versus amino acids, traditional methods of ancestral state reconstruction supported six transitions into freshwaters, two of which led to subsequent species diversification. Evidence from gene trees, protein alignments, and diatom life history together suggest that habitat transitions were largely the product of homoplasy rather than hemiplasy, a condition where transitions occur on branches in gene trees not shared with the species tree. Nevertheless, we identified a set of putatively hemiplasious genes, many of which have been associated with shifts to low salinity, indicating that hemiplasy played a small but potentially important role in freshwater adaptation. Accounting for differences in evolutionary outcomes, in which some taxa became locked into freshwaters while others were able to return to the ocean or become salinity generalists, might help further distinguish different sources of adaptive mutation in freshwater diatoms.

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

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

Should we hail the Red King? Evolutionary consequences of a mutualistic lifestyle in genomes of lichenized ascomycetes

The Red Queen dynamic is often brought into play for antagonistic relationships. However, the coevolutionary effects of mutualistic interactions, which predict slower evolution for interacting organisms (Red King), have been investigated to a lesser extent. Lichens are a stable, mutualistic relationship of fungi and cyanobacteria and/or algae, which originated several times independently during the evolution of fungi. Therefore, they represent a suitable system to investigate the coevolutionary effect of mutualism on the fungal genome. We measured substitution rates and selective pressure of about 2000 protein-coding genes (plus the rDNA region) in two different classes of Ascomycota, each consisting of closely related lineages of lichenized and non-lichenized fungi. Our results show that independent lichenized clades are characterized by significantly slower rates for both synonymous and non-synonymous substitutions. We hypothesize that this evolutionary pattern is connected to the lichen life cycle (longer generation time of lichenized fungi) rather than a result of different selection strengths, which is described as the main driver for the Red Kind dynamic. This first empirical evidence of slower evolution in lichens provides an important insight on how biotic cooperative interactions are able to shape the evolution of symbiotic organisms.

Toxicity evaluation of butyl acrylate on the photosynthetic pigments, chlorophyll fluorescence parameters, and oxygen evolution activity of Phaeodactylum tricornutum and Platymonas subcordiformis

Butyl acrylate is a hazardous and noxious substance (HNS) listed in the top 50 chemicals that is most likely to be involved in HNS spilling incident. At present, information about toxicity effect of butyl acrylate on marine organisms was insufficient, especially on marine microalgae. Phaeodactylum tricornutum (P. tricornutum) and Platymonas subcordiformis (P. subcordiformis) were used as test organism to evaluate the toxic effect of butyl acrylate on their photosynthetic system. Results showed that chlorophyll a (Chl-a) content, the net photosynthetic oxygen evolution rate (NOR), and chlorophyll fluorescence parameters including maximal photochemical efficiency (Fv/Fm), electron transfer rate (ETR), photochemical quenching (qP), and non-photochemical quenching (NPQ) were all stimulated in the toxic dose of 5,10, and 25 mg/L while those were significantly inhibited in the highest concentration of 25 mg/L groups after 96 h. Meanwhile, it was also found that Fv/Fm would be a suitable indicator for evaluating the toxicity of butyl acrylate on the photosynthetic system of two marine microalgae according to the analysis of Pearson correlation coefficient and integrated biomarker response (IBR). Once butyl acrylate enters the marine ecosystem, the toxicity data obtained in this study could be used as a reference for evaluating the effect of butyl acrylate on the photosynthetic capacity of marine microalgae.

Effects of fluctuating temperature in open raceway ponds on the biomass accumulation and harvest efficiency of Spirulina in large-scale cultivation

It is of great significance to select strains with wide adaptability to temperature range for large-scale commercial cultivation of Spirulina. The aim of this study was to comprehend how the strain H-208 grew and whether this strain had any advantages in temperature adaptation compared with local production strain during the large-scale cultivation in Inner Mongolia. The results showed that the strain H-208 could adapt to the new environmental condition quickly, and the daily average biomass dry weight of strain H-208 was 49% and 52% more than that of production strain M-1 in first cycle (20.24 g/m2/day) and second cycle (16.90 g/m2/day) of acclimation experiment, respectively. The growth rate of strain H-208 was 0.055 and 0.066 g/L/day from July 22 to July 25 and from July 26 to July 29, respectively, while the growth rate of strain M-1 was only 0.036 and 0.032 g/L/day, respectively, during the same cultured days in 605-m2 raceway ponds before high temperature. The harvesting efficiency of H-208 and M-1 was 95.1% and 72.1% before high temperature, and that was 95.3% and 52.5% after being stressed by high temperature, respectively. Meanwhile, it was also observed that the filaments of the two strains contracted and their pitches were smaller than that before high temperature stress, especially the strain M-1. In 20-m2 raceway ponds of recovery experiment after high temperature, the percentage of daily average biomass dry weight of strain H-208 was 68% more than that of strain M-1, which demonstrated that strain H-208 could recover and grow rapidly, and its self-regulation ability was superior to that of strain M-1.

Analysis of Adaptive Evolution and Coevolution of rbcL Gene in the Genus Hildenbrandia (Rhodophyta)

The adaptive evolution and coevolution of the ribulose-1,5-bisphosphate carboxylase/oxygenase large subunit (rbcL) gene in the genus Hildenbrandia were studied based on phylogenetic tree construction and the physicochemical properties and the secondary structures of protein encoded by rbcL (Rubisco large subunit) were analyzed. The amino acids compositions and grand average of hydropathicity of freshwater H. rivularis and marine H. rubra were similar. Rubisco large subunit of Hildenbrandia was hydrophilic and the secondary structure was primarily composed of α-helixes and β-sheets, revealing the relatively stable structure of this protein. The predicted phosphorylation sites in H. rivularis and H. rubra were 33 and 36, respectively. No positive selection sites were detected in the genus Hildenbrandia, implying that rbcL gene evolved either neutrally or under purifying selection. A total of 41 coevolutionary groups were detected in the Rubisco large subunit of Hildenbrandia and the coevolving sites are in closer proximity in 3-dimensional structure of the protein. Despite the long evolutionary history, rbcL gene in genus Hildenbrandia under different environments is rather conservative.

On the Selection and Analysis of Clades in Comparative Evolutionary Studies

Researchers commonly present results of comparative studies of taxonomic groups. In this review, we criticize the focus on named clades, usually, comparably ranked groups such as families or orders, for comparative evolutionary analyses and question the general practice of using clades as units of analysis. The practice of analyzing sets of named groups persists despite widespread appreciation that the groups we have chosen to name are based on subjective human concerns rather than objective properties of nature. We demonstrate an effect of clade selection on results in one study and present some potential alternatives to selecting named clades for analysis that are relatively objective in clade choice. However, we note that these alternatives are only partial solutions for clade-based studies. The practice of analyzing named clades obviously is biased and problematic, but its issues portend broader problems with the general approach of employing clades as units of analysis. Most clade-based studies do not account for the nonindependence of clades, and the biological insight gained from demonstrating some pattern among a particular arbitrary sample of groups is arguable. [Clades; comparative biology; taxonomic groups.].

Toxicity assessment of p-choroaniline on Platymonas subcordiformis and its biodegradation

The use of p-chloroaniline (PCA) in various aspects leads to its existence and accumulation in the environment. Relevant researches showed that PCA was a prime toxic pollutant that had imposed a serious risk to public health and the environment. This paper investigated the toxicity effects of PCA on Platymonas subcordiformis (P. subcordiformis) and the biodegradation of PCA by the marine microalga. In the toxicity experiments, the EC₅₀ of PCA on P. subcordiformis at 24 h, 48 h, 72 h and 96 h was 41.42, 24.04, 17.15 and 13.05 mg L^-1, respectively. The pigment parameters including chlorophyll a, chlorophyll b, carotenoids, photosynthetic O₂ release rate, respiration O₂ consumption rate and the chlorophyll fluorescence parameters including Fv/Fm, ETR and qP decreased greatly while antioxidant enzyme activities (SOD, CAT) and the chlorophyll fluorescence parameter NPQ increased when P. subcordiformis exposed to PCA compared with the control group. Fv/Fm would be a suitable indicator for assessing the toxicity of PCA in marine environment based on the analysis of Pearson's correlation coefficient and Integrated Biomarker Response (IBR). The degradation assay in P. subcordiformis indicated that the green marine microalga had the ability to remove and degrade PCA, and the order of removal and degradation proportion of PCA was 2 mg L^-1 > 5 mg L^-1>10 mg L^-1. The maximum removal and biodegradation percentage was 54% and 34%, respectively.

The Molecular Basis of Freshwater Adaptation in Prawns: Insights from Comparative Transcriptomics of Three Macrobrachium Species

Elucidating the molecular basis of adaptation to different environmental conditions is important because adaptive ability of a species can shape its distribution, influence speciation, and also drive a variety of evolutionary processes. For crustaceans, colonization of freshwater habitats has significantly impacted diversity, but the molecular basis of this process is poorly understood. In the current study, we examined three prawn species from the genus Macrobrachium (M. australiense, M. tolmerum, and M. novaehollandiae) to better understand the molecular basis of freshwater adaptation using a comparative transcriptomics approach. Each of these species naturally inhabit environments with different salinity levels; here, we exposed them to the same experimental salinity conditions (0‰ and 15‰), to compare expression patterns of candidate genes that previously have been shown to influence phenotypic traits associated with freshwater adaptation (e.g., genes associated with osmoregulation). Differential gene expression analysis revealed 876, 861, and 925 differentially expressed transcripts under the two salinities for M. australiense, M. tolmerum, and M. novaehollandiae, respectively. Of these, 16 were found to be unannotated novel transcripts and may be taxonomically restricted or orphan genes. Functional enrichment and molecular pathway mapping revealed 13 functionally enriched categories and 11 enriched molecular pathways that were common to the three Macrobrachium species. Pattern of selection analysis revealed 26 genes with signatures of positive selection among pairwise species comparisons. Overall, our results indicate that the same key genes and similar molecular pathways are likely to be involved with freshwater adaptation widely across this decapod group; with nonoverlapping sets of genes showing differential expression (mainly osmoregulatory genes) and signatures of positive selection (genes involved with different life history traits).

Freshwater-to-marine transitions may explain the evolution of herbivory in the subgenus Mollienesia (genus Poecilia, mollies and guppies)

The ability of organisms to cross ecosystem boundaries is an important catalyst of evolutionary diversification. The genus Poecilia (mollies and guppies) is an excellent system for studying ecosystem transitions because species display a range of salinity and dietary preferences, with herbivory concentrated in the subgenus Mollienesia. We reconstructed ancestral habitats and diets across a phylogeny of the genus Poecilia, evaluated diversification rates and used phylogenetically independent contrasts to determine whether diet evolved in response to habitat transition in this group. The results suggest that ancestors of subgenus Mollienesia were exclusively herbivorous, whereas ancestral diets of other Poecilia included animals. We found that transitions across euryhaline boundaries occurred at least once in this group, probably after the divergence of the subgenus Mollienesia. Furthermore, increased salinity affiliation explained 24% of the decrease in animals in the gut, and jaw morphology was associated with the percentage of animals in the gut, but not with the percentage of species occupying saline habitats. These findings suggest that in the genus Poecilia, herbivory evolved in association with transitions from fresh to euryhaline habitats, and jaw morphology evolved in response to the appearance of herbivory. These results provide a rare example of increased diet diversification associated with the transition from freshwater to euryhaline habitats.

Toxicity effects of p-choroaniline on the growth, photosynthesis, respiration capacity and antioxidant enzyme activities of a diatom, Phaeodactylum tricornutu

The environmental risk issues of p-choroaniline have been concerned by the widespread application and transportation of this important chemical intermediate. The information about the toxicity of p-chloroaniline was mainly concentrated on freshwater organisms while the current knowledge on marine organisms was scarce yet. In this study, acute toxicity and toxic physiology characteristic of p-chloroaniline to Phaeodactylum tricornutum (P. tricornutum) were first determined. In the acute experiments, the effect of the p-choroaniline to P. tricornutum showed time- and dose-dependent response, which the half maximum effective concentration (EC₅₀) at 24 h, 48 h and 96 h was 35.35, 20.10 and 10.00 mgL^-1, respectively. Toxic physiology assays in P. tricornutum indicated that the p-choroaniline induced significant changes of photosynthetic pigments (Chl-a, Chl-b, Caro, Chl-a/b and Chl-(a+b)/Caro), Chlorophyll fluorescence parameters (Fv/Fm, ETR, qP and NPQ), rates of photosynthetic O₂ release and respiration O₂ consumption, and antioxidant enzyme activities (SOD, CAT). The obvious decrease of Fv/Fm, ETR and chl-a in low p-choroaniline treatments (≤ 5.00 mgL^-1) compared with the control could be observed, which implied that these parameters could be taken as sensitive indicators for the environmental assessment. Meanwhile, the activities of SOD and CAT significant increase in p-choroaniline stress after 24 h and the extent of the increase has fallen after 96 h. These toxicity data obtained here might provide available basic data for the ecological risk assessment of p-choroaniline pollution.

When the ribosomal DNA does not tell the truth: The case of the taxonomic position of Kurtia argillacea, an ericoid mycorrhizal fungus residing among Hymenochaetales

The nuclear ribosomal DNA (nuc-rDNA) is widely used for the identification and phylogenetic reconstruction of Agaricomycetes. However, nuc-rDNA-based phylogenies may sometimes be in conflict with phylogenetic relationships derived from protein coding genes. In this study, the taxonomic position of the basidiomycetous mycobiont that forms the recently discovered sheathed ericoid mycorrhiza was investigated, because its nuc-rDNA is highly dissimilar to any other available fungal sequences in terms of nucleotide composition and length, and its nuc-rDNA-based phylogeny is inconclusive and significantly disagrees with protein coding sequences and morphological data. In the present work, this mycobiont was identified as Kurtia argillacea (= Hyphoderma argillaceum) residing in the order Hymenochaetales (Basidiomycota). Bioinformatic screening of the Kurtia ribosomal DNA sequence indicates that it represents a gene with a non-standard substitution rate or nucleotide composition heterogeneity rather than a deep paralogue or a pseudogene. Such a phenomenon probably also occurs in other lineages of the Fungi and should be taken into consideration when nuc-rDNA (especially that with unusual nucleotide composition) is used as a sole marker for phylogenetic reconstructions. Kurtia argillacea so far represents the only confirmed non-sebacinoid ericoid mycorrhizal fungus in the Basidiomycota and its intriguing placement among mostly saprobic and parasitic Hymenochaetales begs further investigation of its eco-physiology.

De Novo Transcriptome Characterization of a Sterilizing Trematode Parasite (Microphallus sp.) from Two Species of New Zealand Snails

Snail-borne trematodes represent a large, diverse, and evolutionarily, ecologically, and medically important group of parasites, often imposing strong selection on their hosts and causing host morbidity and mortality. Even so, there are very few genomic and transcriptomic resources available for this important animal group. We help to fill this gap by providing transcriptome resources from trematode metacercariae infecting two congeneric snail species, Potamopyrgus antipodarum and P. estuarinus. This genus of New Zealand snails has gained prominence in large part through the development of P. antipodarum and its sterilizing trematode parasite Microphallus livelyi into a textbook model for host–parasite coevolutionary interactions in nature. By contrast, the interactions between Microphallus trematodes and P. estuarinus, an estuary-inhabiting species closely related to the freshwater P. antipodarum, are relatively unstudied. Here, we provide the first annotated transcriptome assemblies from Microphallus isolated from P. antipodarum and P. estuarinus. We also use these transcriptomes to produce genomic resources that will be broadly useful to those interested in host–parasite coevolution, local adaption, and molecular evolution and phylogenetics of this and other snail–trematode systems. Analyses of the two Microphallus transcriptomes revealed that the two trematode types are more genetically differentiated from one another than are the M. livelyi infecting different populations of P. antipodarum, suggesting that the Microphallus infecting P. estuarinus represent a distinct lineage. We also provide a promising set of candidate genes likely involved in parasitic infection and response to salinity stress.

Rates and patterns of molecular evolution in freshwater versus terrestrial insects

Insect lineages have crossed between terrestrial and aquatic habitats many times, for both immature and adult life stages. We explore patterns in molecular evolutionary rates between 42 sister pairs of related terrestrial and freshwater insect clades using publicly available protein-coding DNA sequence data from the orders Coleoptera, Diptera, Lepidoptera, Hemiptera, Mecoptera, Trichoptera, and Neuroptera. We furthermore test for habitat-associated convergent molecular evolution in the cytochrome c oxidase subunit I (COI) gene in general and at a particular amino acid site previously reported to exhibit habitat-linked convergence within an aquatic beetle group. While ratios of nonsynonymous-to-synonymous substitutions across available loci were higher in terrestrial than freshwater-associated taxa in 26 of 42 lineage pairs, a stronger trend was observed (20 of 31, p_binomial = 0.15, p_Wilcoxon = 0.017) when examining only terrestrial-aquatic pairs including fully aquatic taxa. We did not observe any widespread changes at particular amino acid sites in COI associated with habitat shifts, although there may be general differences in selection regime linked to habitat.