Dissemination of methicillin-resistant Staphylococcus aureus USA300 sequence type 8 lineage in Latin America.

BACKGROUND: Methicillin-resistant Staphylococus aureus (MRSA) is an important nosocomial and community-associated (CA) pathogen. Recently, a variant of the MRSA USA300 clone emerged and disseminated in South America, causing important clinical problems. METHODS: S. aureus isolates were prospectively collected (2006-2008) from 32 tertiary hospitals in Colombia, Ecuador, Peru, and Venezuela. MRSA isolates were subjected to antimicrobial susceptibility testing and pulsed-field gel electrophoresis and were categorized as health care-associated (HA)-like or CA-like clones on the basis of genotypic characteristics and detection of genes encoding Panton-Valentine leukocidin and staphylococcal cassette chromosome (SCC) mec IV. In addition, multilocus sequence typing of representative isolates of each major CA-MRSA pulsotype was performed, and the presence of USA300-associated toxins and the arcA gene was investigated for all isolates categorized as CA-MRSA. RESULTS: A total of 1570 S. aureus were included; 651 were MRSA (41%)--with the highest rate of MRSA isolation in Peru (62%) and the lowest in Venezuela (26%)--and 71%, 27%, and 2% were classified as HA-like, CA-like, and non-CA/HA-like clones, respectively. Only 9 MRSA isolates were confirmed to have reduced susceptibility to glycopeptides (glycopeptide-intermediate S. aureus phenotype). The most common pulsotype (designated ComA) among the CA-like MRSA strains was found in 96% of isolates, with the majority (81%) having a < or =6-band difference with the USA300-0114 strain. Representative isolates of this clone were sequence type 8; however, unlike the USA300-0114 strain, they harbored a different SCCmec IV subtype and lacked arcA (an indicator of the arginine catabolic mobile element). CONCLUSION: A variant CA-MRSA USA300 clone has become established in South America and, in some countries, is endemic in hospital settings.

Analysis of clonality and antibiotic resistance among early clinical isolates of Enterococcus faecium in the United States.

BACKGROUND: The Enterococcus faecium genogroup, referred to as clonal complex 17 (CC17), seems to possess multiple determinants that increase its ability to survive and cause disease in nosocomial environments. METHODS: Using 53 clinical and geographically diverse US E. faecium isolates dating from 1971 to 1994, we determined the multilocus sequence type; the presence of 16 putative virulence genes (hyl(Efm), esp(Efm), and fms genes); resistance to ampicillin (AMP) and vancomycin (VAN); and high-level resistance to gentamicin and streptomycin. RESULTS: Overall, 16 different sequence types (STs), mostly CC17 isolates, were identified in 9 different regions of the United States. The earliest CC17 isolates were part of an outbreak that occurred in 1982 in Richmond, Virginia. The characteristics of CC17 isolates included increases in resistance to AMP, the presence of hyl(Efm) and esp(Efm), emergence of resistance to VAN, and the presence of at least 13 of 14 fms genes. Eight of 41 of the early isolates with resistance to AMP, however, were not in CC17. CONCLUSIONS: Although not all early US AMP isolates were clonally related, E. faecium CC17 isolates have been circulating in the United States since at least 1982 and appear to have progressively acquired additional virulence and antibiotic resistance determinants, perhaps explaining the recent success of this species in the hospital environment.

Clinical and microbiological aspects of linezolid resistance mediated by the cfr gene encoding a 23S rRNA methyltransferase.

The cfr (chloramphenicol-florfenicol resistance) gene encodes a 23S rRNA methyltransferase that confers resistance to linezolid. Detection of linezolid resistance was evaluated in the first cfr-carrying human hospital isolate of linezolid and methicillin-resistant Staphylococcus aureus (designated MRSA CM-05) by dilution and diffusion methods (including Etest). The presence of cfr was investigated in isolates of staphylococci colonizing the patient's household contacts and clinical isolates recovered from patients in the same unit where MRSA CM-05 was isolated. Additionally, 68 chloramphenicol-resistant Colombian MRSA isolates recovered from hospitals between 2001 and 2004 were screened for the presence of the cfr gene. In addition to erm(B), the erm(A) gene was also detected in CM-05. The isolate belonged to sequence type 5 and carried staphylococcal chromosomal cassette mec type I. We were unable to detect the cfr gene in any of the human staphylococci screened (either clinical or colonizing isolates). Agar and broth dilution methods detected linezolid resistance in CM-05. However, the Etest and disk diffusion methods failed to detect resistance after 24 h of incubation. Oxazolidinone resistance mediated by the cfr gene is rare, and acquisition by a human isolate appears to be a recent event in Colombia. The detection of cfr-mediated linezolid resistance might be compromised by the use of the disk diffusion or Etest method.

Identification and characterization of the herpes simplex virus type 2 gene encoding the essential capsid protein ICP32/VP19c.

We describe the characterization of the herpes simplex virus type 2 (HSV-2) gene encoding infected cell protein 32 (ICP32) and virion protein 19c (VP19c). We also demonstrate that the HSV-1 UL38/ORF.553 open reading frame (ORF), which has been shown to specify a viral protein essential for capsid formation (B. Pertuiset, M. Boccara, J. Cebrian, N. Berthelot, S. Chousterman, F. Puvian-Dutilleul, J. Sisman, and P. Sheldrick, J. Virol. 63: 2169-2179, 1989), must encode the cognate HSV type 1 (HSV-1) ICP32/VP19c protein. The region of the HSV-2 genome deduced to contain the gene specifying ICP32/VP19c was isolated and subcloned, and the nucleotide sequence of 2,158 base pairs of HSV-2 DNA mapping immediately upstream of the gene encoding the large subunit of the viral ribonucleotide reductase was determined. This region of the HSV-2 genome contains a large ORF capable of encoding two related 50,538- and 49,472-molecular-weight polypeptides. Direct evidence that this ORF encodes HSV-2 ICP32/VP19c was provided by immunoblotting experiments that utilized antisera directed against synthetic oligopeptides corresponding to internal portions of the predicted polypeptides encoded by the HSV-2 ORF or antisera directed against a TrpE/HSV-2 ORF fusion protein. The type-common immunoreactivity of the two antisera and comparison of the primary amino acid sequences of the predicted products of the HSV-2 ORF and the equivalent genomic region of HSV-1 provided evidence that the HSV-1 UL38 ORF encodes the HSV-1 ICP32/VP19c. Analysis of the expression of the HSV-1 and HSV-2 ICP32/VP19c cognate proteins indicated that there may be differences in their modes of synthesis. Comparison of the predicted structure of the HSV-2 ICP32/VP19c protein with the structures of related proteins encoded by other herpes viruses suggested that the internal capsid architecture of the herpes family of viruses varies substantially.

NMR structure of a complex between the VirB9/VirB7 interaction domains of the pKM101 type IV secretion system.

Type IV secretion (T4S) systems translocate DNA and protein effectors through the double membrane of Gram-negative bacteria. The paradigmatic T4S system in Agrobacterium tumefaciens is assembled from 11 VirB subunits and VirD4. Two subunits, VirB9 and VirB7, form an important stabilizing complex in the outer membrane. We describe here the NMR structure of a complex between the C-terminal domain of the VirB9 homolog TraO (TraO(CT)), bound to VirB7-like TraN from plasmid pKM101. TraO(CT) forms a beta-sandwich around which TraN winds. Structure-based mutations in VirB7 and VirB9 of A. tumefaciens show that the heterodimer interface is conserved. Opposite this interface, the TraO structure shows a protruding three-stranded beta-appendage, and here, we supply evidence that the corresponding region of VirB9 of A. tumefaciens inserts in the membrane and protrudes extracellularly. This complex structure elucidates the molecular basis for the interaction between two essential components of a T4S system.

Evidence for VirB4-mediated dislocation of membrane-integrated VirB2 pilin during biogenesis of the Agrobacterium VirB/VirD4 type IV secretion system.

Agrobacterium VirB2 pilin is required for assembly of the VirB/VirD4 type IV secretion system (T4SS). The propilin is processed by signal sequence cleavage and covalent linkage of the N and C termini, and the cyclized pilin integrates into the inner membrane (IM) as a pool for assembly of the secretion channel and T pilus. Here, by use of the substituted cysteine accessibility method (SCAM), we defined the VirB2 IM topology and then identified distinct contributions of the T4SS ATPase subunits to the pilin structural organization. Labeling patterns of Cys-substituted pilins exposed to the membrane-impermeative, thiol-reactive reagent 3-(N-maleimidopropionyl)biocytin (MPB) supported a topology model in which two hydrophobic stretches comprise transmembrane domains, an intervening hydrophilic loop (residues 90 to 94) is cytoplasmic, and the hydrophilic N and C termini joined at residues 48 and 121 form a periplasmic loop. Interestingly, the VirB4 ATPase, but not a Walker A nucleoside triphosphate (NTP) binding motif mutant, induced (i) MPB labeling of Cys94, a residue that in the absence of the ATPase is located in the cytoplasmic loop, and (ii) release of pilin from the IM upon osmotic shock. These findings, coupled with evidence for VirB2-VirB4 complex formation by coimmunoprecipitation, support a model in which VirB4 functions as a dislocation motor to extract pilins from the IM during T4SS biogenesis. The VirB11 ATPase functioned together with VirB4 to induce a structural change in the pilin that was detectable by MPB labeling, suggestive of a role for VirB11 as a modulator of VirB4 dislocase activity.

A trilocus sequence typing scheme for hospital epidemiology and subspecies differentiation of an important nosocomial pathogen, Enterococcus faecalis.

In this study, we present a trilocus sequence typing (TLST) scheme based on intragenic regions of two antigenic genes, ace and salA (encoding a collagen/laminin adhesin and a cell wall-associated antigen, respectively), and a gene associated with antibiotic resistance, lsa (encoding a putative ABC transporter), for subspecies differentiation of Enterococcus faecalis. Each of the alleles was analyzed using 50 E. faecalis isolates representing 42 diverse multilocus sequence types (ST(M); based on seven housekeeping genes) and four groups of clonally linked (by pulsed-field gel electrophoresis [PFGE]) isolates. The allelic profiles and/or concatenated sequences of the three genes agreed with multilocus sequence typing (MLST) results for typing of 49 of the 50 isolates; in addition to the one exception, two isolates were found to have identical TLST types but were single-locus variants (differing by a single nucleotide) by MLST and were therefore also classified as clonally related by MLST. TLST was also comparable to PFGE for establishing short-term epidemiological relationships, typing all isolates classified as clonally related by PFGE with the same type. TLST was then applied to representative isolates (of each PFGE subtype and isolation year) of a collection of 48 hospital isolates and demonstrated the same relationships between isolates of an outbreak strain as those found by MLST and PFGE. In conclusion, the TLST scheme described here was shown to be successful for investigating short-term epidemiology in a hospital setting and may provide an alternative to MLST for discriminating isolates.

Enterococcus faecalis PcfC, a spatially localized substrate receptor for type IV secretion of the pCF10 transfer intermediate.

Upon sensing of peptide pheromone, Enterococcus faecalis efficiently transfers plasmid pCF10 through a type IV secretion (T4S) system to recipient cells. The PcfF accessory factor and PcfG relaxase initiate transfer by catalyzing strand-specific nicking at the pCF10 origin of transfer sequence (oriT). Here, we present evidence that PcfF and PcfG spatially coordinate docking of the pCF10 transfer intermediate with PcfC, a membrane-bound putative ATPase related to the coupling proteins of gram-negative T4S machines. PcfC and PcfG fractionated with the membrane and PcfF with the cytoplasm, yet all three proteins formed several punctate foci at the peripheries of pheromone-induced cells as monitored by immunofluorescence microscopy. A PcfC Walker A nucleoside triphosphate (NTP) binding site mutant (K156T) fractionated with the E. faecalis membrane and also formed foci, whereas PcfC deleted of its N-terminal putative transmembrane domain (PcfCDelta N103) distributed uniformly throughout the cytoplasm. Native PcfC and mutant proteins PcfCK156T and PcfCDelta N103 bound pCF10 but not pcfG or Delta oriT mutant plasmids as shown by transfer DNA immunoprecipitation, indicating that PcfC binds only the processed form of pCF10 in vivo. Finally, purified PcfCDelta N103 bound DNA substrates and interacted with purified PcfF and PcfG in vitro. Our findings support a model in which (i) PcfF recruits PcfG to oriT to catalyze T-strand nicking, (ii) PcfF and PcfG spatially position the relaxosome at the cell membrane to stimulate substrate docking with PcfC, and (iii) PcfC initiates substrate transfer through the pCF10 T4S channel by an NTP-dependent mechanism.

Mutagenicity of herpes simplex virus type 1 in a shuttle vector

The shuttle vector plasmid pZ189 was used to find the kinds of mutations that are induced by herpes simplex virus type-1 (HSV-1). In cells infected by HSV-1 the frequency of mutation in supF gene, the mutagenesis marker, was increased over background by from two- to seven-fold, reaching 0.14-0.45%. No increase was induced by infection by vaccinia virus under the same conditions. Mutagenesis was an early event, showing a four-fold increase in mutation frequency at only two hours after infection, and peaking at a seven-fold increase at four hours after infection. DNA sequencing and gel electrophoresis analysis were performed on 105 HSV-1 induced mutants and 65 spontaneous mutants and provided the following information: (1) A change in plasmid size was seen in 54% of HSV-1 related mutants, compared with only 37% of spontaneous mutants. (2) Among point mutations, the predominant type was G:C to A:T transition, which accounted for 51% of point mutations in mutants isolated from cells infected with HSV-1, and 32% of point mutations in spontaneous mutants. (3) Deletions of DNA were seen in HSV-1 related mutants at a frequency of 40%, compared with 29% in spontaneous mutants. The HSV-1 related deletions were about half the length of spontaneous mutants and three contained short filler sequences. (4) Fifteen (15%) of HSV-1 induced mutants revealed the altered restriction patterns on agarose gel electrophoresis analysis and were due either to rearrangements of plasmid DNA, and/or to insertion of sequences derived from chromosomal DNA (seven plasmids). No insertions of DNA from HSV-1 were detected. Among spontaneous mutants, only 5 (7.7%) were rearrangements and none had inserted chromosomal DNA. (5) DNA sequence analysis of seven plasmids with inserted chromosomal DNA revealed that four cases had repetitive DNA sequences integrated and the other three were unidentified sequences from the GenBank database. Three repetitive DNA included $\alpha$ satellite, Alu and KpnI family sequences. The other sequence was identified as tRNA-like component. The observed mutations have implications for the mechanism of malignant transformation of cells by HSV-1. ^

{\it In vivo\/} induction of cytotoxic T lymphocytes against human immunodeficiency virus type 1 by synthetic viral peptides

Cytotoxic T lymphocytes (CTLs) play an important role in the suppression of initial viremia after acute infection with the human immunodeficiency virus (HIV), the causative agent of acquired immune deficiency syndrome (AIDS). Most HIV-infected individuals attain a high titer of anti-HIV antibodies within weeks of infection; however this antibody-mediated immune response appears not to be protective. In addition, anti-HIV antibodies can be detrimental to the immune response to HIV through enhancement of infection and participating in autoimmune reactions as a result of HIV protein mimicry of self antigens. Thus induction and maintenance of a strong HIV-specific CTL immune response in the absence of anti-HIV antibodies has been proposed to be the most effective means of controlling of HIV infection. Immunization with synthetic peptides representing HIV-specific CTL epitopes provides a way to induce specific CTL responses, while avoiding stimulation of anti-HIV antibody. This dissertation examines the capacity of synthetic peptides from the V3 loop region of the gp120 envelope protein from several different strain of HIV-1 to induce HIV-specific, MHC-restricted CD8$\sp+$ CTL response in vivo in a mouse model. Seven synthetic peptides representative of sequences found throughout North America, Europe, and Central Africa have been shown to prime CTLs in vivo. In the case of the MN strain of HIV-1, a 13 amino acid sequence defining the epitope is most efficient for optimal induction of specific CTL, whereas eight to nine amino acid sequences that could define the epitope were not immunogenic. In addition, synthesis of peptides with specific amino acid substitutions that are important for either MHC binding or T cell receptor recognition resulted in peptides that exhibited increased immunogenicity and induced CTLs that displayed altered specificity. V3 loop peptides from HIV-1 MN, SC, and Z321 induced a CTL population that was broadly cross-reactive against strains of HIV-1 found throughout the world. This research confirms the potential efficacy of using synthetic peptides for in vivo immunization to induce HIV-specific CTL-mediated responses and provides a basis for further research into development of synthetic peptide-based vaccines. ^

REPSA: A method for determining the sequence specificity of DNA binding ligands

We developed a novel combinatorial method termed restriction endonuclease protection selection and amplification (REPSA) to identify consensus binding sites of DNA-binding ligands. REPSA uses a unique enzymatic selection based on the inhibition of cleavage by a type IIS restriction endonuclease, an enzyme that cleaves DNA at a site distal from its recognition sequence. Sequences bound by a ligand are protected from cleavage while unprotected sequences are cleaved. This enzymatic selection occurs in solution under mild conditions and is dependant only on the DNA-binding ability of the ligand. Thus, REPSA is useful for a broad range of ligands including all classes of DNA-binding ligands, weakly binding ligands, mixed populations of ligands, and unknown ligands. Here I describe REPSA and the application of this method to select the consensus DNA-binding sequences of three representative DNA-binding ligands; a nucleic acid (triplex-forming single-stranded DNA), a protein (the TATA-binding protein), and a small molecule (Distamycin A). These studies generated new information regarding the specificity of these ligands in addition to establishing their DNA-binding sequences. ^

Non B-DNA structures and genomic instability associated with myotonic dystrophy type 2 (CCTG).(CAGG) repeats

There are many diseases associated with the expansion of DNA repeats in humans. Myotonic dystrophy type 2 is one of such diseases, characterized by expansions of a (CCTG)•(CAGG) repeat tract in intron 1 of zinc finger protein 9 (ZNF9) in chromosome 3q21.3. The DM2 repeat tract contains a flanking region 5' to the tract that consists of a polymorphic repetitive sequence (TG)14-25(TCTG)4-11(CCTG) n. The (CCTG)•(CAGG) repeat is typically 11-26 repeats in persons without the disease, but can expand up to 11,000 repeats in affected individuals, which is the largest expansion seen in DNA repeat diseases to date. This DNA tract remains one of the least characterized disease-associated DNA repeats, and mechanisms causing the repeat expansion in humans have yet to be elucidated. Alternative, non B-DNA structures formed by the expanded repeats are typical in DNA repeat expansion diseases. These sequences may promote instability of the repeat tracts. I determined that slipped strand structure formation occurs for (CCTG)•(CAGG) repeats at a length of 42 or more. In addition, Z-DNA structure forms in the flanking human sequence adjacent to the (CCTG)•(CAGG) repeat tract. I have also performed genetic assays in E. coli cells and results indicate that the (CCTG)•(CAGG) repeats are more similar to the highly unstable (CTG)•(CAG) repeat tracts seen in Huntington's disease and myotonic dystrophy type 1, than to those of the more stable (ATTCT)•(AGAAT) repeat tracts of spinocerebellar ataxia type 10. This instability, however, is RecA-independent in the (CCTG)•(CAGG) and (ATTCT)•(AGAAT) repeats, whereas the instability is RecA-dependent in the (CTG)•(CAG) repeats. Structural studies of the (CCTG)•(CAGG) repeat tract and the flanking sequence, as well as genetic selection assays may reveal the mechanisms responsible for the repeat instability in E. coli, and this may lead to a better understanding of the mechanisms contributing to the human disease state. ^

The chronobiology of the menstrual life sequence of American women

Seasonal variation in menarche, menstrual cycle length and menopause was investigated using Tremin Trust data. Too, self-reported hot flash data for women with natural and surgically-induced menopause were analyzed for rhythms.^ Menarche data from approximately 600 U.S. women born between 1940 and 1970 revealed a 6-month rhythm (first acrophase in January, double amplitude of 58%M). A notable shift from a December-January peak in menarche for those born in the 1940s and 1950s to an August-September peak for those born in the 1960s was observed. Groups of girls 8-14 and 15-17 yr old at menarche exhibited a seasonal difference in the pattern of menarche occurrence of about 6 months in relation to each other. Girls experiencing menarche during August-October were statistically significantly younger than those experiencing it at other times. Season of birth was not associated with season of menarche.^ The lengths of approximately 150,000 menstrual intervals of U.S. women were analyzed for seasonality. Menstrual intervals possibly disturbed by natural (e.g., childbirth) or other events (e.g., surgery, medication) were excluded. No 6- or 12-month rhythmicities were found for specific interval lengths (14-24, 25-31 and 32-56 days) or ages in relation to menstrual interval (9-11, 12-13, 15-19, 20-24, 25-39, 40-44 and 44 yr old and older).^ Hot flash data of 14 women experiencing natural menopause (NM) and 11 experiencing surgically-induced menopause (SIM) did not differ in frequency of hot flashes. Hot flashes in NM women exhibited 12- and 8-hr, but not 24-hr rhythmicities. Hot flashes in SIM women exhibited 24- and 12-hr, but not 8-hr, rhythmicities. Regardless of type of menopause, women with a peak frequency in hot flashes during the morning (0400 through 0950) were distinguishable from those with such in the evening (1600 through 2159).^ Data from approximately 200 U.S. women revealed a 6-month rhythm in menopause with first peak in May. No significant 12-month variation in menopause was detected by Cosinor analysis. Season of birth and age at menopause were not associated with season of menopause. Age at menopause declined significantly over the years for women born between 1907 and 1926, inclusive. ^

Transcriptional regulation of cell lineage-specific expression of the murine type I collagen genes and of osteoblast differentiation

The aim of my project is to examine the mechanisms of cell lineage-specific transcriptional regulation of the two type I collagen genes by characterizing critical cis-acting elements and trans-acting factors. I hypothesize that the transcription factors that are involved in the cell lineage-specific expression of these genes may have a larger essential role in cell lineage commitment and differentiation. I first examined the proximal promoters of the proα1(I) and the proα2(I) collagen genes for cell type-specific DNA-protein interactions, using in vitro DNaseI and in vivo DMS footprinting. These experiments demonstrated that the cis-acting elements in these promoters are accessible to ubiquitous DNA-binding proteins in fibroblasts that express these genes, but not in other cells that do not express these genes. I speculate that in type I collagen-expressing cells, cell type-specific enhancer elements facilitate binding of ubiquitous proteins to the proximal promoters of these genes. Subsequently, examination of the upstream promoter of the proα(I) collagen gene by transgenic mice experiments delineated a 117 bp sequence (-1656 to -1540 bp) as the minimum element required for osteoblast-specific expression. This 117 bp element contained two segments that appeared to have different functions: (1) the A-segment, which was necessary to obtain osteoblast-specific expression and (2) the C-segment, which was dispensable for osteoblast-specific expression, but was necessary to obtain high-level expression. In experiments to identify trans-acting factors that bind to the 117 bp element, I have demonstrated that the cell lineage-restricted homeodomain proteins, Dlx2, Dlx5 and mHOX, bound to the A-segment and that the ubiquitous transcription factor, Sp1, bound to the C-segment of this element. These results suggested a model where the binding of cell lineage-restricted proteins to the A-segment and of ubiquitous proteins to the C-segment of the 117 bp element of the proα1 (I) collagen gene activated this gene in osteoblasts. These results, combined with additional evidence that Dlx2, Dlx5 and mHOX are probably involved in osteoblast differentiation, support my hypothesis that the transcription factors involved in osteoblast-specific expression of type I collagen genes may have essential role in osteoblast lineage commitment and differentiation. ^

Evolution of cis-regulatory elements in a repetitive sequence adjacent to a sea urchin aboral ectoderm -specific gene

The creation, preservation, and degeneration of cis-regulatory elements controlling developmental gene expression are fundamental genome-level evolutionary processes about which little is known. In this study, critical differences in cis-regulatory elements controlling the expression of the sea urchin aboral ectoderm-specific spec genes were identified and explored. In genomes of species within the Strongylocentrotidae family, multiple copies of a repetitive sequence element termed RSR were present, but RSRs were not detected in genomes of species outside Strongylocentrotidae. RSRs are invariably associated with spec genes, and in Strongylocentrotus purpuratus, the spec2a RSR functioned as a transcriptional enhancer displaying greater activity than RSRs from the spec1 or spec2c paralogs. Single base-pair differences at two cis-regulatory elements within the spec2a RSR greatly increased the binding affinities of four transcription factors: SpCCAAT-binding factor at one element and SpOtx, SpGoosecoid, and SpGATA-E at another. The cis-regulatory elements to which SpCCAAT-binding factor, SpOtx, SpGoosecoid, and SpGATA-E bound were recent evolutionary acquisitions that could act either to activate or repress transcription, depending on the cell type. These elements were found in the spec2a RSR ortholog in Strongylocentrotus pallidus but not in the RSR orthologs of Strongylocentrotus droebachiensis or Hemicentrotus pulcherrimus. These results indicate that spec genes exhibit a dynamic pattern of cis-regulatory element evolution while stabilizing selection preserves their aboral ectoderm expression domain. ^