Sall4 regulates downstream patterning genes during limb regeneration.

Many salamanders can completely regenerate a fully functional limb. Limb regeneration is a carefully coordinated process involving several defined stages. One key event during the regeneration process is the patterning of the blastema to inform cells of what they must differentiate into. Although it is known that many genes involved in the initial development of the limb are re-used during regeneration, the exact molecular circuitry involved in this process is not fully understood. Several large-scale transcriptional profiling studies of axolotl limb regeneration have identified many transcription factors that are up-regulated after limb amputation. Sall4 is a transcription factor that has been identified to play essential roles in maintaining cells in an undifferentiated state during development and also plays a unique role in limb development. Inactivation of Sall4 during limb bud development results in defects in anterior-posterior patterning of the limb. Sall4 has been found to be up-regulated during limb regeneration in both Xenopus and salamanders, but to date it function has been untested. We confirmed that Sall4 is up-regulated during limb regeneration in the axolotl using qRT-PCR and identified that it is present in the skin cells and also in cells within the blastema. Using CRISPR technology we microinjected gRNAs specific for Sall4 complexed with cas9 protein into the blastema to specifically knockout Sall4 in blastema cells only. This resulted in limb regenerate defects, including missing digits, fusion of digit elements, and defects in the radius and ulna. This suggests that during regeneration Sall4 may play a similar role in regulating the specification of anterior-proximal skeletal elements.

An esculentin-1 homolog from a dark-spotted frog (Pelophylax nigromaculatus) possesses antibacterial and immunoregulatory properties.

BACKGROUND: Esculentin-1, initially discovered in the skin secretions of pool frogs (Pelophylax lessonae), has demonstrated broad-spectrum antimicrobial activity; however, its immunomodulatory properties have received little attention. RESULTS: In the present study, esculentin-1 cDNA was identified by analysing the skin transcriptome of the dark-spotted frog (Pelophylax nigromaculatus). Esculentin-1 from this species (esculentin-1PN) encompasses a signal peptide, an acidic spacer peptide, and a mature peptide. Sequence alignments with other amphibian esculentins-1 demonstrated conservation of the peptide, and phylogenetic tree analysis revealed its closest genetic affinity to esculentin-1P, derived from the Fukien gold-striped pond frog (Pelophylax fukienensis). Esculentin-1PN transcripts were observed in various tissues, with the skin exhibiting the highest mRNA levels. Synthetic esculentin-1PN demonstrated antibacterial activity against various pathogens, and esculentin-1PN exhibited bactericidal activity by disrupting cell membrane integrity and hydrolyzing genomic DNA. Esculentin-1PN did not stimulate chemotaxis in RAW264.7, a murine leukemic monocyte/macrophage cell line. However, it amplified the respiratory burst and augmented the pro-inflammatory cytokine gene (TNF-α and IL-1ß) expression in RAW264.7 cells. CONCLUSIONS: This novel finding highlights the immunomodulatory activity of esculentin-1PN on immune cells.

VCP/p97 Extracts Sterically Trapped Ku70/80 Rings from DNA in Double-Strand Break Repair.

During DNA double-strand break (DSB) repair, the ring-shaped Ku70/80 complex becomes trapped on DNA and needs to be actively extracted, but it has remained unclear what provides the required energy. By means of reconstitution of DSB repair on beads, we demonstrate here that DNA-locked Ku rings are released by the AAA-ATPase p97. To achieve this, p97 requires ATP hydrolysis, cooperates with the Ufd1-Npl4 ubiquitin-adaptor complex, and specifically targets Ku80 that is modified by K48-linked ubiquitin chains. In U2OS cells, chemical inhibition of p97 or siRNA-mediated depletion of p97 or its adapters impairs Ku80 removal after non-homologous end joining of DSBs. Moreover, this inhibition attenuates early steps in homologous recombination, consistent with p97-driven Ku release also affecting repair pathway choice. Thus, our data answer a central question regarding regulation of Ku in DSB repair and illustrate the ability of p97 to segregate even tightly bound protein complexes for release from DNA.

Drivers and suppressors of triple-negative breast cancer.

To identify regulators of triple-negative breast cancer (TNBC), gene expression profiles of malignant parts of TNBC (mTNBC) and normal adjacent (nadj) parts of the same breasts have been compared. We are interested in the roles of estrogen receptor ß (ERß) and the cytochrome P450 family (CYPs) as drivers of TNBC. We examined by RNA sequencing the mTNBC and nadj parts of five women. We found more than a fivefold elevation in mTNBC of genes already known to be expressed in TNBC: BIRC5/survivin, Wnt-10A and -7B, matrix metalloproteinases (MMPs), chemokines, anterior gradient proteins, and lysophosphatidic acid receptor and the known basal characteristics of TNBC, sox10, ROPN1B, and Col9a3. There were two unexpected findings: 1) a strong induction of CYPs involved in activation of fatty acids (CYP4), and in inactivation of calcitriol (CYP24A1) and retinoic acid (CYP26A1); and 2) a marked down-regulation of FOS, FRA1, and JUN, known tethering partners of ERß. ERß is expressed in 20 to 30% of TNBCs and is being evaluated as a target for treating TNBC. We used ERß+ TNBC patient-derived xenografts in mice and found that the ERß agonist LY500703 had no effect on growth or proliferation. Expression of CYPs was confirmed by immunohistochemistry in formalin-fixed and paraffin-embedded (FFPE) TNBC. In TNBC cell lines, the CYP4Z1-catalyzed fatty acid metabolite 20-hydroxyeicosatetraenoic acid (20-HETE) increased proliferation, while calcitriol decreased proliferation but only after inhibition of CYP24A1. We conclude that CYP-mediated pathways can be drivers of TNBC but that ERß is unlikely to be a tumor suppressor because the absence of its main tethering partners renders ERß functionless on genes involved in proliferation and inflammation.

Mutations of the Transcriptional Corepressor ZMYM2 Cause Syndromic Urinary Tract Malformations.

Congenital anomalies of the kidney and urinary tract (CAKUT) constitute one of the most frequent birth defects and represent the most common cause of chronic kidney disease in the first three decades of life. Despite the discovery of dozens of monogenic causes of CAKUT, most pathogenic pathways remain elusive. We performed whole-exome sequencing (WES) in 551 individuals with CAKUT and identified a heterozygous de novo stop-gain variant in ZMYM2 in two different families with CAKUT. Through collaboration, we identified in total 14 different heterozygous loss-of-function mutations in ZMYM2 in 15 unrelated families. Most mutations occurred de novo, indicating possible interference with reproductive function. Human disease features are replicated in X. tropicalis larvae with morpholino knockdowns, in which expression of truncated ZMYM2 proteins, based on individual mutations, failed to rescue renal and craniofacial defects. Moreover, heterozygous Zmym2-deficient mice recapitulated features of CAKUT with high penetrance. The ZMYM2 protein is a component of a transcriptional corepressor complex recently linked to the silencing of developmentally regulated endogenous retrovirus elements. Using protein-protein interaction assays, we show that ZMYM2 interacts with additional epigenetic silencing complexes, as well as confirming that it binds to FOXP1, a transcription factor that has also been linked to CAKUT. In summary, our findings establish that loss-of-function mutations of ZMYM2, and potentially that of other proteins in its interactome, as causes of human CAKUT, offering new routes for studying the pathogenesis of the disorder.

BMP signaling is essential for sustaining proximo-distal progression in regenerating axolotl limbs.

Amputation of a salamander limb triggers a regeneration process that is perfect. A limited number of genes have been studied in this context and even fewer have been analyzed functionally. In this work, we use the BMP signaling inhibitor LDN193189 on Ambystoma mexicanum to explore the role of BMPs in regeneration. We find that BMP signaling is required for proper expression of various patterning genes and that its inhibition causes major defects in the regenerated limbs. Fgf8 is downregulated when BMP signaling is blocked, but ectopic injection of either human or axolotl protein did not rescue the defects. By administering LDN193189 treatments at different time points during regeneration, we show clearly that limb regeneration progresses in a proximal to distal fashion. This demonstrates that BMPs play a major role in patterning of regenerated limbs and that regeneration is a progressive process like development.

An Aerolysin-like Pore-Forming Protein Complex Targets Viral Envelope to Inactivate Herpes Simplex Virus Type 1.

Because most of animal viruses are enveloped, cytoplasmic entry of these viruses via fusion with cellular membrane initiates their invasion. However, the strategies in which host cells counteract cytoplasmic entry of such viruses are incompletely understood. Pore-forming toxin aerolysin-like proteins (ALPs) exist throughout the animal kingdom, but their functions are mostly unknown. In this study, we report that ßγ-crystallin fused aerolysin-like protein and trefoil factor complex (ßγ-CAT), an ALP and trefoil factor complex from the frog Bombina maxima, directly blocks enveloped virus invasion by interfering with cytoplasmic entry. ßγ-CAT targeted acidic glycosphingolipids on the HSV type 1 (HSV-1) envelope to induce pore formation, as indicated by the oligomer formation of protein and potassium and calcium ion efflux. Meanwhile, ßγ-CAT formed ring-like oligomers of ∼10 nm in diameter on the liposomes and induced dye release from liposomes that mimic viral envelope. Unexpectedly, transmission electron microscopy analysis showed that the ßγ-CAT-treated HSV-1 was visibly as intact as the vehicle-treated HSV-1, indicating that ßγ-CAT did not lyse the viral envelope. However, the cytoplasmic entry of the ßγ-CAT-treated HSV-1 into HeLa cells was totally hindered. In vivo, topical application of ßγ-CAT attenuated the HSV-1 corneal infection in mice. Collectively, these results uncovered that ßγ-CAT possesses the capacity to counteract enveloped virus invasion with its featured antiviral-acting manner. Our findings will also largely help to illustrate the putative antiviral activity of animal ALPs.

Conserved keratin gene clusters in ancient fish: An evolutionary seed for terrestrial adaptation.

Type I and type II keratins are subgroups of intermediate filament proteins that provide toughness to the epidermis and protect it from water loss. In terrestrial vertebrates, the keratin genes form two major clusters, clusters 1 and 2, each of which is dominated by type I and II keratin genes. By contrast, such clusters are not observed in teleost fish. Although the diversification of keratins is believed to have made a substantial contribution to terrestrial adaptation, its evolutionary process has not been clarified. Here, we performed a comprehensive genomic survey of the keratin genes of a broad range of vertebrates. As a result, we found that ancient fish lineages such as elephant shark, reedfish, spotted gar, and coelacanth share both keratin gene clusters. We also discovered an expansion of keratin genes that form a novel subcluster in reedfish. Syntenic and phylogenetic analyses revealed that two pairs of krt18/krt8 keratin genes were shared among all vertebrates, thus implying that they encode ancestral type I and II keratin protein sets. We further revealed that distinct keratin gene subclusters, which show specific expressions in the epidermis of adult amphibians, stemmed from canonical keratin genes in non-terrestrial ancestors. Molecular evolutionary analyses suggested that the selective constraints were relaxed in the adult epidermal subclusters of amphibians as well as the novel subcluster of reedfish. The results of the present study represent the process of diversification of keratins through a series of gene duplications that could have facilitated the terrestrial adaptation of vertebrates.

Codon usage pattern and evolutionary forces of mitochondrial ND genes among orders of class Amphibia.

In this study, we used a bioinformatics approach to analyze the nucleotide composition and pattern of synonymous codon usage in mitochondrial ND genes in three amphibian groups, that is, orders Anura, Caudata, and Gymnophiona to identify the commonality and the differences of codon usage as no research work was reported yet. The high value of the effective number of codons revealed that the codon usage bias (CUB) was low in mitochondrial ND genes among the orders. Nucleotide composition analysis suggested that for each gene, the compositional features differed among Anura, Caudata, and Gymnophiona and the GC content was lower than AT content. Furthermore, a highly significant difference (p < .05) for GC content was found in each gene among the orders. The heat map showed contrasting patterns of codon usage among different ND genes. The regression of GC12 on GC3 suggested a narrow range of GC3 distribution and some points were located in the diagonal, indicating both mutation pressure and natural selection might influence the CUB. Moreover, the slope of the regression line was less than 0.5 in all ND genes among orders, indicating natural selection might have played the dominant role whereas mutation pressure had played a minor role in shaping CUB of ND genes across orders.

A novel strategy for engineering of a smart generation of immune ribonucleases against EGFR+ cells.

The overexpression of epidermal growth factor receptor (EGFR) could result in the development of solid tumors of prostate, breast, gastric, colorectal, ovarian, and head and neck, leading to carcinoma. Antibody therapies are ideal methods to overcome malignant diseases. However, immunoribonucleases are a new generation of antibodies in which an RNase binds to a specific antibody and shows a stronger ability to terminate cancer cells. In this study, we engineered Rana pipiens RNase to bind to the scFv of human antiepidermal growth factor receptor antibody. The molecular dynamic simulations confirmed protein stability and the ability of scFv-ranpirnase (rantoxin) to bind to epidermal growth factor receptor protein. Then, the rantoxin construct was synthesized in a pCDNA 3.1 Neo vector. CHO-K1 cells were used as expression hosts and the construct was transfected. Cells were selected by antibiotic therapies using neomycin, 120 mg/ml, and the high-yield colony was screened by real-time polymerase chain reaction (PCR) methods. Then, the recombinant protein production was confirmed using the sodium dodecyl sulfate polyacrylamide gel electrophoresis and western blot analyses. The molecular dynamic simulation (MDS) confirmed that the I467, S468, Q408, and H409 amino acids of EGFR bonded well to rantoxin. As revealed by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) and western blot analyses, the rantoxin production and PCR analysis showed that the T3 colony can produce rantoxin messenger RNA fourfold higher than the GAPDH gene. The immunotoxin function was assessed in A431 cancer cells and EGFR-negative HEK293 cells, and IC50  values were estimated to be 22.4 ± 3 and >620.4 ± 5 nM, respectively. The results indicated that the immunotoxins produced in this study had the potential for use as anticancer drugs.

Toward an understanding of tree frog (Hyla japonica) for predator deterrence.

Gene-encoded peptides with distinct potent bioactivities enable several animals to take advantage of fierce interspecific interaction, as seen in the skin secretion of amphibians. Unlike, most amphibian species that frequently switches terrestrial-aquatic habitats and hides easily from terrestrial predators, tree frogs of small body size are considered as the vulnerable prey in the arboreal habitat. Here, we show the structural and functional diversity of peptide families based on the skin transcriptome of Hyla japonica, which has evolved to be wrapped as an efficient chemical toolkit for defensive use in arboreal habitat. Generally, the presence of antimicrobial peptide and proteinase inhibitor families reveals the functional consistency of Hyla japonica skin compared to other amphibian species. Furthermore, we found that Anntoxin-like neurotoxins with high expression levels are species-specific in tree frogs. Interestingly, derivatives in the Anntoxin-like family exhibit multiple evolutionary traits in modifying the copy number, folding type, and three-dimensional architecture, which are considered essential for targeting the ion channels of terrestrial predators. Together, our study not only reveals the peptide diversity in the skin secretion of H. japonica, but also draws insights into the predator-deterring strategy for coping with arboreal habitat.

Expression and purification of the native C-amidated antimicrobial peptide maculatin 1.1.

Maculatin 1.1 (Mac1) is an antimicrobial peptide (AMP) from an Australian tree frog and exhibits low micromolar activity against Gram-positive bacteria. The antimicrobial properties of Mac1 are linked to its disruption of bacterial lipid membranes, which has been studied extensively by in vitro nuclear magnetic resonance (NMR) spectroscopy and biophysical approaches. Although in vivo NMR has recently proven effective in probing peptide-lipid interplay in live bacterial cells, direct structural characterisation of AMPs has been prohibited by low sensitivity and overwhelming background noise. To overcome this issue, we report a recombinant expression protocol to produce isotopically enriched Mac1. We utilized a double-fusion construct to alleviate toxicity against the Escherichia coli host and generate the native N-free and C-amidated termini Mac1 peptide. The SUMO and intein tags allowed native N-terminus and C-terminal amidation, respectively, to be achieved in a one-pot reaction. The protocol yielded 0.1 mg/L of native, uniformly 15 N-labelled, Mac1, which possessed identical structure and activity to peptide obtained by solid-phase peptide synthesis.

Amino acid residues important for D-galactose recognition by the F-type lectin, Ranaspumin-4.

F-type lectins are typically L-fucose binding proteins with characteristic L-fucose-binding and calcium-binding sequence motifs, and an F-type lectin fold. An exception is Ranaspumin-4, an F-type lectin of the Tungra frog, Engystomops pustulosus. Ranaspumin-4 is D-galactose specific and does not bind to L-fucose although it has the conserved L-fucose binding sequence motif and shares overall sequence similarity with other F-type lectins. Here, we report the detailed glycan-binding profile of wild-type Ranaspumin-4 using hemagglutination inhibition assays, flow cytometry assays and enzyme-linked lectin assays, and identify residues important for D-galactose recognition using rational site-directed mutagenesis. We demonstrate that Ranaspumin-4 binds to terminal D-galactose in α or ß linkage with preference for α1-3, α1-4, ß1-3, and ß1-4 linkages. Further, we find that a methionine residue (M31) in Ranaspumin-4 that occurs in place of a conserved Gln residue (in other F-type lectins), supports D-galactose recognition. Resides Q42 and F156 also likely aid in D-galactose recognition.

Neutralization of Bombina variegata peptide 8 suppresses retinal neovascularization in two different murine models: The oxygen-induced retinopathy model and the rhodopsin promoter/VEGF transgenic mouse model.

Bombina variegata 8 (Bv8), also known as prokineticin-2 (PK-2), is a potent pro-angiogenic factor. However, its role in retinal neovascularization (RNV) remains unknown. In this study, we explored the role of Bv8 in the pathogenesis of RNV. We found that the expression of Bv8 was significantly increased in two different models of retinal neovascularization: the oxygen-induced retinopathy (OIR) mouse model and the rhodopsin promoter (rho)/VEGF transgenic mouse model. Neutralization of Bv8 by intravitreal injections of its antibody, not only inhibited retinal and subretinal neovascularization but also decreased the mRNA and protein levels of several pro-angiogenic factors. Our in vitro assay showed that recombinant human Bv8 (RhBv8) protein promoted human retinal microvascular endothelial cells (HRECs) tube-formation, cell proliferation and vascular endothelial growth factor receptor 1 (VEGFR1) and receptor 2 (VEGFR2) expression. Our findings suggest that Bv8 could be used as a novel target for the treatment of RNV-related ocular diseases.

Cloning of Hynobius lichenatus (Tohoku hynobiid salamander) p53 and analysis of its expression in response to radiation.

BACKGROUND: Caudata species such as salamanders are easily affected by environmental changes, which can drastically reduce their population. The effects of acute X-rays and chronic γ-irradiation on Hynobius lichenatus, the Japanese Tohoku hynobiid salamander, are known. However, the expression of radiation-inducible genes, such as the DNA-damage checkpoint response gene p53, has not been analyzed in H. lichenatus. This has not occurred because there is no established method for mRNA quantification in H. lichenatus due to a lack of information on available nucleotide sequences corresponding to both radiation-inducible genes and endogenous control genes such as ACTB (ß-actin). RESULTS: In this study, we aimed to evaluate the effects of radiation on gene expression in H. lichenatus. Using RNA extracted from irradiated salamanders, we performed rapid amplification of cDNA ends (RACE) and cloned H. lichenatus ß-actin, glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and p53. We confirmed that the cloned cDNAs were able to synthesize salamander proteins by western blotting after transfection into cultured HEK293 cells. Proliferation assays using HEK293 cells stably expressing H. lichenatus p53 protein showed that this protein has antiproliferative effects, similar to that of mammalian p53. Furthermore, RT-qPCR analysis using gene-specific primers revealed that p53 mRNA expression in H. lichenatus was upregulated upon exposure to radiation. CONCLUSION: Our results suggest that H. lichenatus p53 protein take an important role in regulating the cellular responses to various stimuli as mammalian p53 does. Furthermore, our study provides novel data to select appropriate primers to analyze internal control mRNA expression in H. lichenatus and to evaluate p53 expression as a marker of radiation and environmental stimuli.

Alternative polyadenylation coordinates embryonic development, sexual dimorphism and longitudinal growth in Xenopus tropicalis.

RNA alternative polyadenylation contributes to the complexity of information transfer from genome to phenome, thus amplifying gene function. Here, we report the first X. tropicalis resource with 127,914 alternative polyadenylation (APA) sites derived from embryos and adults. Overall, APA networks play central roles in coordinating the maternal-zygotic transition (MZT) in embryos, sexual dimorphism in adults and longitudinal growth from embryos to adults. APA sites coordinate reprogramming in embryos before the MZT, but developmental events after the MZT due to zygotic genome activation. The APA transcriptomes of young adults are more variable than growing adults and male frog APA transcriptomes are more divergent than females. The APA profiles of young females were similar to embryos before the MZT. Enriched pathways in developing embryos were distinct across the MZT and noticeably segregated from adults. Briefly, our results suggest that the minimal functional units in genomes are alternative transcripts as opposed to genes.

The Xenopus animal cap transcriptome: building a mucociliary epithelium.

With the availability of deep RNA sequencing, model organisms such as Xenopus offer an outstanding opportunity to investigate the genetic basis of vertebrate organ formation from its embryonic beginnings. Here we investigate dynamics of the RNA landscape during formation of the Xenopus tropicalis larval epidermis. Differentiation of non-neural ectoderm starts at gastrulation and takes about one day to produce a functional mucociliary epithelium, highly related to the one in human airways. To obtain RNA expression data, uncontaminated by non-epidermal tissues of the embryo, we use prospective ectodermal explants called Animal Caps (ACs), which differentiate autonomously into a ciliated epidermis. Their global transcriptome is investigated at three key timepoints, with a cumulative sequencing depth of ∼108 reads per developmental stage. This database is provided as online Web Tool to the scientific community. In this paper, we report on global changes in gene expression, an unanticipated diversity of mRNA splicing isoforms, expression patterns of repetitive DNA Elements, and the complexity of circular RNAs during this process. Computationally we derive transcription factor hubs from this data set, which may help in the future to define novel genetic drivers of epidermal differentiation in vertebrates.

Phylogenomics reveals rapid, simultaneous diversification of three major clades of Gondwanan frogs at the Cretaceous-Paleogene boundary.

Frogs (Anura) are one of the most diverse groups of vertebrates and comprise nearly 90% of living amphibian species. Their worldwide distribution and diverse biology make them well-suited for assessing fundamental questions in evolution, ecology, and conservation. However, despite their scientific importance, the evolutionary history and tempo of frog diversification remain poorly understood. By using a molecular dataset of unprecedented size, including 88-kb characters from 95 nuclear genes of 156 frog species, in conjunction with 20 fossil-based calibrations, our analyses result in the most strongly supported phylogeny of all major frog lineages and provide a timescale of frog evolution that suggests much younger divergence times than suggested by earlier studies. Unexpectedly, our divergence-time analyses show that three species-rich clades (Hyloidea, Microhylidae, and Natatanura), which together comprise ∼88% of extant anuran species, simultaneously underwent rapid diversification at the Cretaceous-Paleogene (K-Pg) boundary (KPB). Moreover, anuran families and subfamilies containing arboreal species originated near or after the KPB. These results suggest that the K-Pg mass extinction may have triggered explosive radiations of frogs by creating new ecological opportunities. This phylogeny also reveals relationships such as Microhylidae being sister to all other ranoid frogs and African continental lineages of Natatanura forming a clade that is sister to a clade of Eurasian, Indian, Melanesian, and Malagasy lineages. Biogeographical analyses suggest that the ancestral area of modern frogs was Africa, and their current distribution is largely associated with the breakup of Pangaea and subsequent Gondwanan fragmentation.

Characterization and function of the T-box 1 gene in Chinese giant salamander Andrias davidianus.

We characterized the Andrias davidianus T-box 1 (Tbx1) gene. Tbx1 expression was high in testis and low in other examined tissues. Immunohistochemistry detected tbx1 expression in somatic and germ cells 62 days post-hatching (dph), prior to gonad differentiation. At 210 dph, after gonad differentiation, tbx1 was expressed in spermatogonia and testis somatic cells and in granulosa cells in ovary. Tbx1 expression was up-regulated in ovary after high temperature treatment. In the neomale, tbx1 expression showed a similar profile to normal males, and vice-versa for genetic male. Over-expression of tbx1 in females after injection of TBX1 protein down-regulated the female-biased genes cyp19a and foxl2 and up-regulated the male-biased amh gene. When tbx1 was knocked down by tbx1/siRNA, cyp19a and foxl2 expression was up-regulated, and expression of amh, cyp26a, dmrt1, and wt1 was down-regulated. Results suggest that tbx1 influenced sex-related gene expression and participates in regulation of A. davidianus testis development.

An annual cycle of gene regulation in the red-legged salamander mental gland: from hypertrophy to expression of rapidly evolving pheromones.

BACKGROUND: Cell differentiation is mediated by synchronized waves of coordinated expression for hundreds to thousands of genes, and must be regulated to produce complex tissues and phenotypes. For many animal species, sexual selection has driven the development of elaborate male ornaments, requiring sex-specific differentiation pathways. One such male ornament is the pheromone-producing mental gland of the red-legged salamander (Plethodon shermani). Mental gland development follows an annual cycle of extreme hypertrophy, production of pheromones for the ~ 2 month mating season, and then complete resorption before repeating the process in the following year. At the peak of the mating season, the transcriptional and translational machinery of the mental gland are almost exclusively redirected to the synthesis of rapidly evolving pheromones. Of these pheromones, Plethodontid Modulating Factor (PMF) has experienced an unusual history: following gene duplication, the protein coding sequence diversified from positive sexual selection while the untranslated regions have been conserved by purifying selection. The molecular underpinnings that bridge the processes of gland hypertrophy, pheromone synthesis, and conservation of the untranslated regions remain to be determined. RESULTS: Using Illumina sequencing, we prepared a de novo transcriptome of the mental gland at six stages of development. Differential expression analysis and immunohistochemistry revealed that the mental gland initially adopts a highly proliferative, almost tumor-like phenotype, followed by a rapid increase in pheromone mRNA and protein. One likely player in this transition is Cold Inducible RNA Binding Protein (CIRBP), which selectively and cooperatively binds the highly conserved PMF 3' UTR. CIRBP, along with other proteins associated with stress response, have seemingly been co-opted to aid in mental gland development by helping to regulate pheromone synthesis. CONCLUSIONS: The P. shermani mental gland utilizes a complex system of transcriptional and post-transcriptional gene regulation to facilitate its hypertrophication and pheromone synthesis. The data support the evolutionary interplay of coding and noncoding segments in rapid gene evolution, and necessitate the study of co-evolution between pheromone gene products and their transcriptional/translational regulators. Additionally, the mental gland could be a powerful emerging model of regulated tissue proliferation and subsequent resorption within the dermis and share molecular links to skin cancer biology.