Evolution of actin gene families of sea urchins

A single-cell RNA-seq analysis of Brachyury-expressing cell clusters suggests a morphogenesis-associated signal center of oral ectoderm in sea urchin embryos

Brachyury is a T-box family transcription factor and plays pivotal roles in morphogenesis. In sea urchin embryos, Brachyury, is expressed in the invaginating endoderm, and in the oral ectoderm of the invaginating mouth opening. The oral ectoderm is hypothesized to serve as a signaling center for oral (ventral)-aboral (dorsal) axis formation and to function as a ventral organizer. Our previous results of a single-cell RNA-seq (scRNA-seq) atlas of early Strongylocentrotus purpuratus embryos categorized the constituent cells into 22 clusters, in which the endoderm consists of three clusters and the oral ectoderm four clusters (Foster et al., 2020). Here we examined which clusters of cells expressed Brachyury in relation to the morphogenesis and the identity of the ventral organizer. Our results showed that cells of all three endoderm clusters expressed Brachyury in blastulae. Based on expression profiles of genes involved in the gene regulatory networks (GRNs) of sea urchin embryos, the three clusters are distinguishable, two likely derived from the Veg2 tier and one from the Veg1 tier. On the other hand, of the four oral-ectoderm clusters, cells of two clusters expressed Brachyury at the gastrula stage and genes that are responsible for the ventral organizer at the late blastula stage, but the other two clusters did not. At a single-cell level, most cells of the two oral-ectoderm clusters expressed organizer-related genes, nearly a half of which coincidently expressed Brachyury. This suggests that the ventral organizer contains Brachyury-positive cells which invaginate to form the stomodeum. This scRNA-seq study therefore highlights significant roles of Brachyury-expressing cells in body-plan formation of early sea urchin embryos, though cellular and molecular mechanisms for how Brachyury functions in these processes remain to be elucidated in future studies.

Characterization and phylogeny of two β-cytoskeletal actins fromHemibarbus mylodon(Cyprinidae, Cypriniformes), a threatened fish species in Korea

Complementary DNA and genomic sequences representing two different beta-actins were isolated from a threatened freshwater fish species Hemibarbus mylodon. The beta-actin 1 and 2 encoded an identical number of amino acids (375 aa), and shared 88.8 and 99.7% of identity at coding nucleotide and amino acid levels, respectively. Genomic open reading frame (ORF) sequences of both isoforms contained five translated exons interrupted by four introns with conserved GT/AG exon/intron boundary rule. Semi-quantitative RT-PCR showed that the two isoform mRNAs were ubiquitously detected in all tissues tested, but transcript levels were variable across tissues. Phylogenetic analysis showed that H. mylodon beta-actin 1 and 2 were clustered into two distinct major and minor branches of Cypriniformes, respectively. Comparisons of the 5'-upstream region and 3'-UTR of H. mylodon beta-actin 1 also showed a high degree of homology with those of the major teleost beta-actins and warmblooded vertebrate beta-cytoskeletal actins, suggesting their more recent common origin.

Cloning of black carp beta-actin gene and primarily detecting the function of its promoter region

A 3 338 bp DNA fragment including the open reading frame and 5'-flanking region of beta-actin gene for black carp genome was obtained through PCR amplification. Analysis of the sequencing results indicated the ORF of black carp beta-actin gene encoding a 375 amino acid protein that shares a high degree of conservation to other known actins. The black carp beta-actin sequence showed 100% identity to common carp, grass carp, and zebrafish, 99.2% identity to human and Norway rat beta-actin gene, 98.9% and 98.1% identity to chicken and Kenyan clawed frog beta-actin gene, respectively. The promoter region of black carp beta-actin gene was inserted into the promoterless pEGFP1 vector. The recombinant plasmid was microinjected into the fertilized eggs of mud loach before two-cell stage as well as transfected into HeLa cell line. GFP expression was found in 50% of mud loach embryos and 2/3 HeLa cells. The GFP expression could be observed in every part of the mud loach embryos, and in some embryos, the GFP was expressed in the whole body. Thus, the usefulness of black carp beta-actin promoter as a ubiquitous expression promoter was confirmed using the EGFP as a reporter gene.

Structure, expression, and transcriptional regulation of the Strongylocentrotus franciscanus spec gene family encoding intracellular calcium-binding proteins

The mechanisms by which gene expression patterns emerge during evolution are poorly understood. The sea urchin spec genes offer a useful means to investigate evolutionary mechanisms. Genes of the spec family from Strongylocentrotus purpuratus and Lytechinus pictus have identical patterns of aboral ectoderm-specific expression but exhibit species-specific differences in copy number, genomic structure, temporal expression, and cis-regulatory architecture. Here, we identify spec genes from a phylogenetic intermediate, Strongylocentrotus franciscanus, to gain insight into the evolution of the spec gene family and its transcriptional regulation. We identified two spec genes in the S. franciscanus genome, sfspec1a and sfspec1b, that were orthologous to spec1 from S. purpuratus. sfspec1b transcripts began to accumulate at the blastula stage and became progressively more abundant; this was reminiscent of spec expression in L. pictus but different from that in S. purpuratus. As expected, sfspec1b expression was restricted to aboral ectoderm cells. The six-exon structure of the sfspec1b genomic locus was identical to that of the S. purpuratus spec genes and was bounded by two repeat-spacer-repeat (RSR) repetitive sequence elements, which are conserved features of S. purpuratus spec genes and function as transcriptional enhancers. The enhancer activity of the sfspec1b RSRs was comparable to that of their S. purpuratus counterparts, although the placement and orientation of crucial cis-regulatory elements within the RSRs differed. We discovered a spec gene in S. franciscanus that was only distantly related to other spec genes but was highly conserved in S. purpuratus. Unexpectedly, this gene was expressed exclusively in endoderm lineages. Our results show that the evolution of spec cis-regulatory elements is highly dynamic and that substantial alterations can occur when maintaining or grossly modifying gene expression patterns.

Invertebrate Muscles: Muscle Specific Genes and Proteins

This is the first of a projected series of canonic reviews covering all invertebrate muscle literature prior to 2005 and covers muscle genes and proteins except those involved in excitation-contraction coupling (e.g., the ryanodine receptor) and those forming ligand- and voltage-dependent channels. Two themes are of primary importance. The first is the evolutionary antiquity of muscle proteins. Actin, myosin, and tropomyosin (at least, the presence of other muscle proteins in these organisms has not been examined) exist in muscle-like cells in Radiata, and almost all muscle proteins are present across Bilateria, implying that the first Bilaterian had a complete, or near-complete, complement of present-day muscle proteins. The second is the extraordinary diversity of protein isoforms and genetic mechanisms for producing them. This rich diversity suggests that studying invertebrate muscle proteins and genes can be usefully applied to resolve phylogenetic relationships and to understand protein assembly coevolution. Fully achieving these goals, however, will require examination of a much broader range of species than has been heretofore performed.

On combining protein sequences and nucleic acid sequences in phylogenetic analysis: the homeobox protein case

Amino acid encoding genes contain character state information that may be useful for phylogenetic analysis on at least two levels. The nucleotide sequence and the translated amino acid sequences have both been employed separately as character states for cladistic studies of various taxa, including studies of the genealogy of genes in multigene families. In essence, amino acid sequences and nucleic acid sequences are two different ways of character coding the information in a gene. Silent positions in the nucleotide sequence (first or third positions in codons that can accrue change without changing the identity of the amino acid that the triplet codes for) may accrue change relatively rapidly and become saturated, losing the pattern of historical divergence. On the other hand, non-silent nucleotide alterations and their accompanying amino acid changes may evolve too slowly to reveal relationships among closely related taxa. In general, the dynamics of sequence change in silent and non-silent positions in protein coding genes result in homoplasy and lack of resolution, respectively. We suggest that the combination of nucleic acid and the translated amino acid coded character states into the same data matrix for phylogenetic analysis addresses some of the problems caused by the rapid change of silent nucleotide positions and overall slow rate of change of non-silent nucleotide positions and slowly changing amino acid positions. One major theoretical problem with this approach is the apparent non-independence of the two sources of characters. However, there are at least three possible outcomes when comparing protein coding nucleic acid sequences with their translated amino acids in a phylogenetic context on a codon by codon basis. First, the two character sets for a codon may be entirely congruent with respect to the information they convey about the relationships of a certain set of taxa. Second, one character set may display no information concerning a phylogenetic hypothesis while the other character set may impact information to a hypothesis. These two possibilities are cases of non-independence, however, we argue that congruence in such cases can be thought of as increasing the weight of the particular phylogenetic hypothesis that is supported by those characters. In the third case, the two sources of character information for a particular codon may be entirely incongruent with respect to phylogenetic hypotheses concerning the taxa examined. In this last case the two character sets are independent in that information from neither can predict the character states of the other. Examples of these possibilities are discussed and the general applicability of combining these two sources of information for protein coding genes is presented using sequences from the homeobox region of 46 homeobox genes from Drosophila melanogaster to develop a hypothesis of genealogical relationship of these genes in this large multigene family.

Actin gene family evolution and the phylogeny of coleoid cephalopods (Mollusca: Cephalopoda)

Phylogenetic analysis conducted on a 784-bp fragment of 82 actin gene sequences of 44 coleoid cephalopod taxa, along with results obtained from genomic Southern blot analysis, confirmed the presence of at least three distinct actin loci in coleoids. Actin isoforms were characteri zed through phylogenetic analysis of representative cephalopod sequences from each of the three isoforms, along with translated actin cDNA sequences from a diverse array of metazoan taxa downloaded from GenBank. One of the three isoforms found in cephalopods was closely related to actin sequences expressed in the muscular tissues of other molluscs. A second isoform was most similar to cytoplasmic-specific actin amino acid sequences. The muscle type actins of molluscs were found to be distinct from those of arthropods, suggesting at least two independent derivations of muscle actins in the protostome lineage, although statistical support for this conclusion was lacking. Parsimony and maximum-likelihood analyses of two of the isoforms from which >30 orthologous coleoid sequences had been obtained (one of the cytoplasmic actins and the muscle actin) supported the monophyly of several higher-level coleoid taxa. These included the superorders Octopodiformes and Decapodiformes, the order Octopoda, the octopod suborder Incirrata, and the teuthoid suborder Myopsida. The monophyly of several taxonomic groups within the Decapodiformes was not supported, including the orders Teuthoidea and Sepioidea and the teuthoid suborder Oegopsida. Parametric bootstrap analysis conducted on the simulated cytoplasmic actin data set provided statistical support to reject the monophyly of the Sepioidea. Although parametric bootstrap analysis of the muscle actin isoform did not reject sepioid monophyly at the 5% level, the results (rejection at P: = 0.068) were certainly suggestive of sepioid nonmonophyly.

Identification of the clam species Ruditapes decussatus (Grooved carpet shell), Venerupis pullastra (Pullet carpet shell), and Ruditapes philippinarum (Japanese carpet shell)by PCR-RFLP

PCR-RFLP analysis has been applied to the identification of three clam species: Ruditapes decussatus (grooved carpet shell), Venerupis pullastra (pullet carpet shell), and Ruditapes philippinarum (Japanese carpet shell). PCR amplification was carried out using a set of primers designed from the DNA nucleotide sequences reported for alpha-actins from humans and various animals. Restriction endonuclease analysis based on sequence data of the PCR products of each clam species revealed the presence of species-specific polymorphic sites for MaeIII and RsaI endonucleases. Electrophoretic analysis of the amplicons digested with MaeIII and RsaI produced species-specific profiles that allowed the genetic identification of the three clam species.

Modifications in protein binding to upstream sequences of the sea urchin cytoplasmic actin gene CyIIa in comparison to its linked neighbors, CyI and CyIIb

The sequences corresponding to regions upstream of the ATG and transcription start site of the CyIIa cytoplasmic actin gene of the sea urchin Strongylocentrotus purpuratus were determined and compared to the genomically linked CyI and CyIIb actin genes. Sites of protein-DNA interaction in the CyIIa upstream sequences were identified by DNase I footprinting. The similarity between CyIIa and CyI (and CyIIb) upstream sequences was limited and included a consensus octamer sequence, serum response element (SRE) and some short sequences within the proximal promoter region. The octamer sequence was found to bind protein. A single DNase I hypersensitive site was detected within the SRE and to two flanking nucleotides, but otherwise, the SRE did not appear to be protected. This is in contrast to strong protein binding to the CyIIb SRE. A region in the CyIIa gene with limited identity to the functionally significant protein binding site D in CyI also did not bind protein. Four additional sites of protein-DNA interaction were identified in CyIIa upstream sequences. One of these is similar to a protein binding site previously located in the CyI upstream sequences, whereas the others appear to be unique. These data indicate that the CyIIa upstream sequences differ extensively from those of CyI. The pattern of CyIIa expression is likely a consequence of these alternations in DNA sequence and protein-DNA interactions.

Cloning and sequencing of a full-length sea bream (Sparus aurata) beta-actin cDNA

A full-length cDNA clone encoding beta-actin (beta-actin) was isolated from a sea bream (Sparus aurata) liver cDNA library. Sequencing of this clone reveals an open reading frame encoding a 375 amino acid protein that shares a high degree of conservation to other known actins. The sea bream beta-actin sequence showed 98% identity to carp and human beta-actin and 95% and 94% identity to sea squirt and Dictyostelium cytoplasmic actins, respectively.