Jason Macrander

Simultaneous delimitation of species and quantification of interspecific hybridization in Amazonian peacock cichlids (genus cichla) using multi-locus data.

BackgroundIntrogression likely plays a significant role in evolution, but understanding the extent and consequences of this process requires a clear identification of species boundaries in each focal group. The delimitation of species, however, is a contentious endeavor. This is true not only because of the inadequacy of current tools to identify species lineages, but also because of the inherent ambiguity between natural populations and species paradigms. The result has been a debate about the supremacy of various species concepts and criteria. Here, we utilized multiple separate sources of molecular data, mtDNA, nuclear sequences, and microsatellites, to delimit species under a polytypic species concept (PTSC) and estimate the frequency and genomic extent of introgression in a Neotropical genus of cichlid fishes (Cichla). We compared our inferences of species boundaries and introgression under this paradigm to those when species are identified under a diagnostic species concept (DSC).ResultsWe find that, based on extensive molecular data and an inclusive species concept, 8 separate biological entities should be recognized rather than the 15 described species of Cichla. Under the PTSC, fewer individuals are expected to exhibit hybrid ancestry than under the DSC (~2% vs. ~12%), but a similar number of the species exhibit introgression from at least one other species (75% vs. 60%). Under either species concept, the phylogenetic breadth of introgression in this group is notable, with both sister species and species from different major mtDNA clades exhibiting introgression.ConclusionsIntrogression was observed to be a widespread phenomenon for delimited species in this group. While several instances of introgressive hybridization were observed in anthropogenically altered habitats, most were found in undisturbed natural habitats, suggesting that introgression is a natural but ephemeral part of the evolution of many tropical species. Nevertheless, even transient introgression may facilitate an increase in genetic diversity or transfer of adaptive mutations that have important consequences in the evolution of tropical biodiversity.

A RNA-seq approach to identify putative toxins from acrorhagi in aggressive and non-aggressive Anthopleura elegantissima polyps

BackgroundThe use of venom in intraspecific aggression is uncommon and venom-transmitting structures specifically used for intraspecific competition are found in few lineages of venomous taxa. Next-generation transcriptome sequencing allows robust characterization of venom diversity and exploration of functionally unique tissues. Using a tissue-specific RNA-seq approach, we investigate the venom composition and gene ontology diversity of acrorhagi, specialized structures used in intraspecific competition, in aggressive and non-aggressive polyps of the aggregating sea anemone Anthopleura elegantissima (Cnidaria: Anthozoa: Hexacorallia: Actiniaria: Actiniidae).ResultsCollectively, we generated approximately 450,000 transcripts from acrorhagi of aggressive and non-aggressive polyps. For both transcriptomes we identified 65 candidate sea anemone toxin genes, representing phospholipase A2s, cytolysins, neurotoxins, and acrorhagins. When compared to previously characterized sea anemone toxin assemblages, each transcriptome revealed greater within-species sequence divergence across all toxin types. The transcriptome of the aggressive polyp had a higher abundance of type II voltage gated potassium channel toxins/Kunitz-type protease inhibitors and type II acrorhagins. Using toxin-like proteins from other venomous taxa, we also identified 612 candidate toxin-like transcripts with signaling regions, potentially unidentified secretory toxin-like proteins. Among these, metallopeptidases and cysteine rich (CRISP) candidate transcripts were in high abundance. Furthermore, our gene ontology analyses identified a high prevalence of genes associated with “blood coagulation” and “positive regulation of apoptosis”, as well as “nucleoside: sodium symporter activity” and “ion channel binding”. The resulting assemblage of expressed genes may represent synergistic proteins associated with toxins or proteins related to the morphology and behavior exhibited by the aggressive polyp.ConclusionWe implement a multifaceted approach to investigate the assemblage of expressed genes specifically within acrorhagi, specialized structures used only for intraspecific competition. By combining differential expression, phylogenetic, and gene ontology analyses, we identify several candidate toxins and other potentially important proteins in acrorhagi of A. elegantissima. Although not all of the toxins identified are used in intraspecific competition, our analysis highlights some candidates that may play a vital role in intraspecific competition. Our findings provide a framework for further investigation into components of venom used exclusively for intraspecific competition in acrorhagi-bearing sea anemones and potentially other venomous animals.

Tissue-Specific Venom Composition and Differential Gene Expression in Sea Anemones

Cnidarians represent one of the few groups of venomous animals that lack a centralized venom transmission system. Instead, they are equipped with stinging capsules collectively known as nematocysts. Nematocysts vary in abundance and type across different tissues; however, the venom composition in most species remains unknown. Depending on the tissue type, the venom composition in sea anemones may be vital for predation, defense, or digestion. Using a tissue-specific RNA-seq approach, we characterize the venom assemblage in the tentacles, mesenterial filaments, and column for three species of sea anemone (Anemonia sulcata, Heteractis crispa, and Megalactis griffithsi). These taxa vary with regard to inferred venom potency, symbiont abundance, and nematocyst diversity. We show that there is significant variation in abundance of toxin-like genes across tissues and species. Although the cumulative toxin abundance for the column was consistently the lowest, contributions to the overall toxin assemblage varied considerably among tissues for different toxin types. Our gene ontology (GO) analyses also show sharp contrasts between conserved GO groups emerging from whole transcriptome analysis and tissue-specific expression among GO groups in our differential expression analysis. This study provides a framework for future characterization of tissue-specific venom and other functionally important genes in this lineage of simple bodied animals.

Multi-copy venom genes hidden in de novo transcriptome assemblies, a cautionary tale with the snakelocks sea anemone Anemonia sulcata (Pennant, 1977)

Using a partial transcriptome of the snakelocks anemone (Anemonia sulcata) we identify toxin gene candidates that were incorrectly assembled into several Trinity components. Our approach recovers hidden diversity found within some toxin gene families that would otherwise go undetected when using Trinity, a widely used program for venom-focused transcriptome reconstructions. Unidentified hidden transcripts may significantly impact conclusions made regarding venom composition (or other multi-copy conserved genes) when using Trinity or other de novo assembly programs.

Permanent Genetic Resources added to Molecular Ecology Resources Database 1 April 2012 - 31 May 2012: PERMANENT GENETIC RESOURCES NOTE

This article documents the addition of 123 microsatellite marker loci to the Molecular Ecology Resources Database. Loci were developed for the following species: Brenthis ino, Cichla orinocensis, Cichla temensis, Epinephelus striatus, Gobio gobio, Liocarcinus depurator, Macrolophus pygmaeus, Monilinia vaccinii‐corymbosi, Pelochelys cantorii, Philotrypesis josephi, Romanogobio vladykovi, Takydromus luyeanus and Takydromus viridipunctatus. These loci were cross‐tested on the following species: Cichla intermedia, Cichla ocellaris, Cichla pinima, Epinephelus acanthistius, Gobio carpathicus, Gobio obtusirostris, Gobio sp. 1, Gobio volgensis, Macrolophus costalis, Macrolophus melanotoma, Macrolophus pygmaeus, Romanogobio albipinnatus, Romanogobio banaticus, Romanogobio belingi, Romanogobio kesslerii, Romanogobio parvus, Romanogobio pentatrichus, Romanogobio uranoscopus, Takydromus formosanus, Takydromus hsuehshanesis and Takydromus stejnegeri.

Evolution of the Cytolytic Pore-Forming Proteins (Actinoporins) in Sea Anemones

Sea anemones (Cnidaria, Anthozoa, and Actiniaria) use toxic peptides to incapacitate and immobilize prey and to deter potential predators. Their toxin arsenal is complex, targeting a variety of functionally important protein complexes and macromolecules involved in cellular homeostasis. Among these, actinoporins are one of the better characterized toxins; these venom proteins form a pore in cellular membranes containing sphingomyelin. We used a combined bioinformatic and phylogenetic approach to investigate how actinoporins have evolved across three superfamilies of sea anemones (Actinioidea, Metridioidea, and Actinostoloidea). Our analysis identified 90 candidate actinoporins across 20 species. We also found clusters of six actinoporin-like genes in five species of sea anemone (Nematostella vectensis, Stomphia coccinea, Epiactis japonica, Heteractis crispa, and Diadumene leucolena); these actinoporin-like sequences resembled actinoporins but have a higher sequence similarity with toxins from fungi, cone snails, and Hydra. Comparative analysis of the candidate actinoporins highlighted variable and conserved regions within actinoporins that may pertain to functional variation. Although multiple residues are involved in initiating sphingomyelin recognition and membrane binding, there is a high rate of replacement for a specific tryptophan with leucine (W112L) and other hydrophobic residues. Residues thought to be involved with oligomerization were variable, while those forming the phosphocholine (POC) binding site and the N-terminal region involved with cell membrane penetration were highly conserved.

Population genetics of the speckled peacock bass (Cichla temensis), South America’s most important inland sport fishery

Abstract The Neotropics harbor the world’s most diverse freshwater fish fauna, with many of these species supporting major commercial, subsistence, or sport fisheries. Knowledge of population genetic structure is available for very few Neotropical fishes, thereby restricting management. To address this need, we examined population genetic variation in mtDNA control region sequences and twelve microsatellite loci in the speckled or barred peacock bass, Cichla temensis. Moderate and statistically significant genetic divergence among localities indicates that migration is low in this species, implying that populations inhabiting tributaries or even smaller spatial units should constitute management units. Analysis of molecular variance of mtDNA sequences identified six areas with largely exclusive haplotype clades, and a seventh area of high admixture, but major drainage basins harbored non-monophyletic haplotype groups. On the other hand, molecular variation in the microsatellite data was best explained by drainage basin and, subsequently, by the seven areas. Populations in these seven areas could be considered evolutionarily significant units (ESUs), and, therefore, we tested hypotheses explaining the discordant signal of mtDNA and microsatellite data using approximate Bayesian computation. This analysis indicated that the divergence of mtDNA clades preceded the divergence of contemporary ESUs across basins, with subsequent lineage sorting among ESUs due to reduced gene flow. Available genetic and ecological information indicates that C. temensis populations of major tributary rivers should be managed as separate stocks that likely are adapted to local environmental conditions.

Structure, folding and stability of a minimal homologue from Anemonia sulcata of the sea anemone potassium channel blocker ShK

Graphical abstract AsK132958 is a 29‐residue peptide identified in a transcriptomic study of Anemonia sulcata. It has the same disulfide framework and a similar structure to ShK. AsK132958 is not active against KV1.3 channels, owing to the lack of a Lys‐Tyr dyad and other functionally important amino acid residues. AsK132958 is more resistant to proteolysis than ShK. Introducing a Lys‐Tyr functional dyad to the AsK132958 structural scaffold may be a useful way of developing a proteolytically stable KV1.3 blocker. Figure. No caption available. HighlightsAsK132958 is one of the shortest peptides with a ShK/BgK‐like cysteine framework.AsK132958 is a structural homologue of ShK.Despite having an ShK‐like scaffold, AsK132958 is not active against KV1.3 or related potassium channels.AsK132958 is more resistant to proteases than ShK and is a promising scaffold for engineering other activities.AsK132958 could be an evolutionary precursor of peptides with ShK‐like scaffold and activity. ABSTRACT Peptide toxins elaborated by sea anemones target various ion‐channel sub‐types. Recent transcriptomic studies of sea anemones have identified several novel candidate peptides, some of which have cysteine frameworks identical to those of previously reported sequences. One such peptide is AsK132958, which was identified in a transcriptomic study of Anemonia sulcata and has a cysteine framework similar to that of ShK from Stichodactyla helianthus, but is six amino acid residues shorter. We have determined the solution structure of this novel peptide using NMR spectroscopy. The disulfide connectivities and structural scaffold of AsK132958 are very similar to those of ShK but the structure is more constrained. Toxicity assays were performed using grass shrimp (Palaemonetes sp) and Artemia nauplii, and patch‐clamp electrophysiology assays were performed to assess the activity of AsK132958 against a range of voltage‐gated potassium (KV) channels. AsK132958 showed no activity against grass shrimp, Artemia nauplii, or any of the KV channels tested, owing partly to the absence of a functional Lys‐Tyr dyad. Three AsK132958 analogues, each containing a Tyr in the vicinity of Lys19, were therefore generated in an effort to restore binding, but none showed activity against any of KV channels tested. However, AsK132958 and its analogues are less susceptible to proteolysis than that of ShK. Our structure suggests that Lys19, which might be expected to occupy the pore of the channel, is not sufficiently accessible for binding, and therefore that AsK132958 must have a distinct functional role that does not involve KV channels.

Conservation of DNA and Ligand Binding Properties of Retinoid X Receptor from the Placozoan Trichoplax adhaerens to human

Nuclear receptors are a superfamily of transcription factors restricted to animals. These transcription factors regulate a wide variety of genes with diverse roles in cellular homeostasis, development, and physiology. The origin and specificity of ligand binding within lineages of nuclear receptors (e.g., subfamilies) continues to be a focus of investigation geared toward understanding how the functions of these proteins were shaped over evolutionary history. Among early-diverging animal lineages, the retinoid X receptor (RXR) is first detected in the placozoan, Trichoplax adhaerens. To gain insight into RXR evolution, we characterized ligand- and DNA-binding activity of the RXR from T. adhaerens (TaRXR). Like bilaterian RXRs, TaRXR specifically bound 9-cis-retinoic acid, which is consistent with a recently published result and supports a conclusion that the ancestral RXR bound ligand. DNA binding site specificity of TaRXR was determined through protein binding microarrays (PBMs) and compared with human RXRɑ. The binding sites for these two RXR proteins were broadly conserved (∼85% shared high-affinity sequences within a targeted array), suggesting evolutionary constraint for the regulation of downstream genes. We searched for predicted binding motifs of the T. adhaerens genome within 1000 bases of annotated genes to identify potential regulatory targets. We identified 648 unique protein coding regions with predicted TaRXR binding sites that had diverse predicted functions, with enriched processes related to intracellular signal transduction and protein transport. Together, our data support hypotheses that the original RXR protein in animals bound a ligand with structural similarity to 9-cis-retinoic acid; the DNA motif recognized by RXR has changed little in more than 1 billion years of evolution; and the suite of processes regulated by this transcription factor diversified early in animal evolution.

Corrigendum to “Structure, folding and stability of a minimal homologue from Anemonia sulcata of the sea anemone potassium channel blocker ShK” [Peptides 99 (2018) 169–178](S0196978117303054)(10.1016/j.peptides.2017.10.001)

Corrigendum to “Structure, folding and stability of a minimal homologue from Anemonia sulcata of the sea anemone potassium channel blocker ShK” [Peptides 99 (2018) 169–178] Bankala Krishnarjuna, Christopher A. MacRaild, Punnepalli Sunanda, Rodrigo A.V. Morales, Steve Peigneur, Jason Macrander, Heidi H. Yu, Marymegan Daly, Srinivasarao Raghothama, Vikas Dhawan, Satendra Chauhan, Jan Tytgat, Michael W. Pennington, Raymond S. Nortona,