Chromosome-scale genome assembly of the brown anole (Anolis sagrei), an emerging model species.

Rapid technological improvements are democratizing access to high quality, chromosome-scale genome assemblies. No longer the domain of only the most highly studied model organisms, now non-traditional and emerging model species can be genome-enabled using a combination of sequencing technologies and assembly software. Consequently, old ideas built on sparse sampling across the tree of life have recently been amended in the face of genomic data drawn from a growing number of high-quality reference genomes. Arguably the most valuable are those long-studied species for which much is already known about their biology; what many term emerging model species. Here, we report a highly complete chromosome-scale genome assembly for the brown anole, Anolis sagrei - a lizard species widely studied across a variety of disciplines and for which a high-quality reference genome was long overdue. This assembly exceeds the vast majority of existing reptile and snake genomes in contiguity (N50 = 253.6 Mb) and annotation completeness. Through the analysis of this genome and population resequence data, we examine the history of repetitive element accumulation, identify the X chromosome, and propose a hypothesis for the evolutionary history of fusions between autosomes and the X that led to the sex chromosomes of A. sagrei.

A single locus regulates a female-limited color pattern polymorphism in a reptile.

Animal coloration is often expressed in periodic patterns that can arise from differential cell migration, yet how these processes are regulated remains elusive. We show that a female-limited polymorphism in dorsal patterning (diamond/chevron) in the brown anole is controlled by a single Mendelian locus. This locus contains the gene CCDC170 that is adjacent to, and coexpressed with, the Estrogen receptor-1 gene, explaining why the polymorphism is female limited. CCDC170 is an organizer of the Golgi-microtubule network underlying a cell's ability to migrate, and the two segregating alleles encode structurally different proteins. Our agent-based modeling of skin development demonstrates that, in principle, a change in cell migratory behaviors is sufficient to switch between the two morphs. These results suggest that CCDC170 might have been co-opted as a switch between color patterning morphs, likely by modulating cell migratory behaviors.

Changes in selection pressure can facilitate hybridization during biological invasion in a Cuban lizard.

Hybridization is among the evolutionary mechanisms most frequently hypothesized to drive the success of invasive species, in part because hybrids are common in invasive populations. One explanation for this pattern is that biological invasions coincide with a change in selection pressures that limit hybridization in the native range. To investigate this possibility, we studied the introduction of the brown anole (Anolis sagrei) in the southeastern United States. We find that native populations are highly genetically structured. In contrast, all invasive populations show evidence of hybridization among native-range lineages. Temporal sampling in the invasive range spanning 15 y showed that invasive genetic structure has stabilized, indicating that large-scale contemporary gene flow is limited among invasive populations and that hybrid ancestry is maintained. Additionally, our results are consistent with hybrid persistence in invasive populations resulting from changes in natural selection that occurred during invasion. Specifically, we identify a large-effect X chromosome locus associated with variation in limb length, a well-known adaptive trait in anoles, and show that this locus is often under selection in the native range, but rarely so in the invasive range. Moreover, we find that the effect size of alleles at this locus on limb length is much reduced in hybrids among divergent lineages, consistent with epistatic interactions. Thus, in the native range, epistasis manifested in hybrids can strengthen extrinsic postmating isolation. Together, our findings show how a change in natural selection can contribute to an increase in hybridization in invasive populations.

CRISPR-Cas9 Gene Editing in Lizards through Microinjection of Unfertilized Oocytes.

CRISPR-Cas9-mediated gene editing has enabled the direct manipulation of gene function in many species. However, the reproductive biology of reptiles presents unique barriers for the use of this technology, and there are no reptiles with effective methods for targeted mutagenesis. Here, we demonstrate that the microinjection of immature oocytes within the ovaries of Anolis sagrei females enables the production of CRISPR-Cas9-induced mutations. This method is capable of producing F0 embryos and hatchlings with monoallelic or biallelic mutations. We demonstrate that these mutations can be transmitted through the germline to establish genetically modified strains of lizards. Direct tests of gene function can now be performed in Anolis lizards, an important model for studies of reptile evolution and development.

Pigeon foot feathering reveals conserved limb identity networks.

The tetrapod limb is a stunning example of evolutionary diversity, with dramatic variation not only among distantly related species, but also between the serially homologous forelimbs (FLs) and hindlimbs (HLs) within species. Despite this variation, highly conserved genetic and developmental programs underlie limb development and identity in all tetrapods, raising the question of how limb diversification is generated from a conserved toolkit. In some breeds of domestic pigeon, shifts in the expression of two conserved limb identity transcription factors, PITX1 and TBX5, are associated with the formation of feathered HLs with partial FL identity. To determine how modulation of PITX1 and TBX5 expression affects downstream gene expression, we compared the transcriptomes of embryonic limb buds from pigeons with scaled and feathered HLs. We identified a set of differentially expressed genes enriched for genes encoding transcription factors, extracellular matrix proteins, and components of developmental signaling pathways with important roles in limb development. A subset of the genes that distinguish scaled and feathered HLs are also differentially expressed between FL and scaled HL buds in pigeons, pinpointing a set of gene expression changes downstream of PITX1 and TBX5 in the partial transformation from HL to FL identity. We extended our analyses by comparing pigeon limb bud transcriptomes to chicken, anole lizard, and mammalian datasets to identify deeply conserved PITX1- and TBX5-responsive components of the limb identity program. Our analyses reveal a suite of predominantly low-level gene expression changes that are conserved across amniotes to regulate the identity of morphologically distinct limbs.

Bypassing Glutamic Acid Decarboxylase 1 (Gad1) Induced Craniofacial Defects with a Photoactivatable Translation Blocker Morpholino.

γ-Amino butyric acid (GABA) mediated signaling is critical in the central and enteric nervous systems, pancreas, lungs, and other tissues. It is associated with many neurological disorders and craniofacial development. Glutamic acid decarboxylase (GAD) synthesizes GABA from glutamate, and knockdown of the gad1 gene results in craniofacial defects that are lethal in zebrafish. To bypass this and enable observation of the neurological defects resulting from knocking down gad1 expression, a photoactivatable morpholino oligonucleotide (MO) against gad1 was prepared by cyclization with a photocleavable linker rendering the MO inactive. The cyclized MO was stable in the dark and toward degradative enzymes and was completely linearized upon brief exposure to 405 nm light. In the course of investigating the function of the ccMOs in zebrafish, we discovered that zebrafish possess paralogous gad1 genes, gad1a and gad1b. A gad1b MO injected at the 1-4 cell stage caused severe morphological defects in head development, which could be bypassed, enabling the fish to develop normally, if the fish were injected with a photoactivatable, cyclized gad1b MO and grown in the dark. At 1 day post fertilization (dpf), light activation of the gad1b MO followed by observation at 3 and 7 dpf led to increased and abnormal electrophysiological brain activity compared to wild type animals. The photocleavable linker can be used to cyclize and inactivate any MO, and represents a general strategy to parse the function of developmentally important genes in a spatiotemporal manner.

A PAGE screening approach for identifying CRISPR-Cas9-induced mutations in zebrafish.

The introduction of CRISPR-Cas9 technology for targeted mutagenesis has revolutionized reverse genetics and made genome editing a realistic option in many model organisms. One of the difficulties with this technique is screening for mutations in large numbers of samples. Many screening approaches for identifying CRISPR-Cas9 mutants have been published; however, in practice these methods are time consuming, expensive, or often yield false positives. This report describes a PCR-based screening approach using non-denaturing PAGE. This approach does not depend on the formation of heteroduplexes and reliably detects changes as small as 1 base-pair (bp) in nucleic acid length at the target site. This approach can be used to identify novel mutations and is also useful as a routine genotyping method.

PITX1 promotes chondrogenesis and myogenesis in mouse hindlimbs through conserved regulatory targets.

The PITX1 transcription factor is expressed during hindlimb development, where it plays a critical role in directing hindlimb growth and the specification of hindlimb morphology. While it is known that PITX1 regulates hindlimb formation, in part, through activation of the Tbx4 gene, other transcriptional targets remain to be elucidated. We have used a combination of ChIP-seq and RNA-seq to investigate enhancer regions and target genes that are directly regulated by PITX1 in embryonic mouse hindlimbs. In addition, we have analyzed PITX1 binding sites in hindlimbs of Anolis lizards to identify ancient PITX1 regulatory targets. We find that PITX1-bound regions in both mouse and Anolis hindlimbs are strongly associated with genes implicated in limb and skeletal system development. Gene expression analyses reveal a large number of misexpressed genes in the hindlimbs of Pitx1-/- mouse embryos. By intersecting misexpressed genes with genes that have neighboring mouse PITX1 binding sites, we identified 440 candidate targets of PITX1. Of these candidates, 68 exhibit ultra-conserved PITX1 binding events that are shared between mouse and Anolis hindlimbs. Among the ancient targets of PITX1 are important regulators of cartilage and skeletal muscle development, including Sox9 and Six1. Our data suggest that PITX1 promotes chondrogenesis and myogenesis in the hindlimb by direct regulation of several key members of the cartilage and muscle transcriptional networks.

Isl1 mediates mesenchymal expansion in the developing external genitalia via regulation of Bmp4, Fgf10 and Wnt5a.

Genital malformations are among the most common human birth defects, and both genetic and environmental factors can contribute to these malformations. Development of the external genitalia in mammals relies on complex signaling networks, and disruption of these signaling pathways can lead to genital defects. Islet-1 (ISL1), a member of the LIM/Homeobox family of transcription factors, has been identified as a major susceptibility gene for classic bladder exstrophy in humans, a common form of the bladder exstrophy-epispadias complex (BEEC), and is implicated in a role in urinary tract development. We report that deletion of Isl1 from the genital mesenchyme in mice led to hypoplasia of the genital tubercle and prepuce, with an ectopic urethral opening and epispadias-like phenotype. These mice also developed hydroureter and hydronephrosis. Identification of ISL1 transcriptional targets via ChIP-Seq and expression analyses revealed that Isl1 regulates several important signaling pathways during embryonic genital development, including the BMP, WNT, and FGF cascades. An essential function of Isl1 during development of the external genitalia is to induce Bmp4-mediated apoptosis in the genital mesenchyme. Together, these studies demonstrate that Isl1 plays a critical role during development of the external genitalia and forms the basis for a greater understanding of the molecular mechanisms underlying the pathogenesis of BEEC and urinary tract defects in humans.

Molecular shifts in limb identity underlie development of feathered feet in two domestic avian species.

Birds display remarkable diversity in the distribution and morphology of scales and feathers on their feet, yet the genetic and developmental mechanisms governing this diversity remain unknown. Domestic pigeons have striking variation in foot feathering within a single species, providing a tractable model to investigate the molecular basis of skin appendage differences. We found that feathered feet in pigeons result from a partial transformation from hindlimb to forelimb identity mediated by cis-regulatory changes in the genes encoding the hindlimb-specific transcription factor Pitx1 and forelimb-specific transcription factor Tbx5. We also found that ectopic expression of Tbx5 is associated with foot feathers in chickens, suggesting similar molecular pathways underlie phenotypic convergence between these two species. These results show how changes in expression of regional patterning genes can generate localized changes in organ fate and morphology, and provide viable molecular mechanisms for diversity in hindlimb scale and feather distribution.

Shared Enhancer Activity in the Limbs and Phallus and Functional Divergence of a Limb-Genital cis-Regulatory Element in Snakes.

The amniote phallus and limbs differ dramatically in their morphologies but share patterns of signaling and gene expression in early development. Thus far, the extent to which genital and limb transcriptional networks also share cis-regulatory elements has remained unexplored. We show that many limb enhancers are retained in snake genomes, suggesting that these elements may function in non-limb tissues. Consistent with this, our analysis of cis-regulatory activity in mice and Anolis lizards reveals that patterns of enhancer activity in embryonic limbs and genitalia overlap heavily. In mice, deletion of HLEB, an enhancer of Tbx4, produces defects in hindlimbs and genitalia, establishing the importance of this limb-genital enhancer for development of these different appendages. Further analyses demonstrate that the HLEB of snakes has lost hindlimb enhancer function while retaining genital activity. Our findings identify roles for Tbx4 in genital development and highlight deep similarities in cis-regulatory activity between limbs and genitalia.

Conserved regulation of hoxc11 by pitx1 in Anolis lizards.

Anolis lizards are an emerging model system for the study of limb development and evolution, but very little is known concerning the regulatory interactions that control limb patterning differences among Anolis species or what regulatory interactions are deeply conserved between Anolis and other tetrapod groups. Here we report the establishment of an embryonic limb micromass culture system that enables functional studies of forelimb and hindlimb gene regulatory networks in Anolis. Characterization of this culture system demonstrated that embryonic forelimb and hindlimb micromasses from different Anolis species are easy to sustain in culture for weeks, and the expression of forelimb and hindlimb-specific gene expression patterns are maintained for at least 8 days in culture. We tested the ability of this system to explore regulatory linkages between transcription factors and their putative target genes through the ectopic expression of a hindlimb-specific transcription factor, pitx1, in forelimb micromasses. We found that pitx1 expression in forelimb cells is sufficient to strongly induce the expression of hoxc11, a gene that normally exhibits hindlimb-restricted expression.

Clinical implications for loss or diminution of expression of Raf-1 kinase inhibitory protein and its phosphorylated form in ductal breast cancer.

Raf Kinase inhibitory protein (RKIP) is a well-established metastasis suppressor that is frequently downregulated in aggressive cancers. The impact of RKIP and its phosphorylated form on disease-free survival (DFS) and other clinicopathological parameters in breast cancer is yet to be discovered. To this end, we examined RKIP expression in 3 independent breast cancer cohorts. At the Protein level, loss or reduced total RKIP expression was associated with large-sized tumors characterized by high proliferative index, high-grade and diminished estrogen (ER) and progesterone receptor expression. Loss or diminution of RKIP expression was significantly associated with shorter DFS in all cohorts. Moreover, the complete loss of p-RKIP was an independent prognostic factor using multivariate analysis in operable invasive ductal breast cancer. We show for the first time that ER, partly, drives RKIP expression through MTA3-Snail axis. Consistent with this finding, we found that, at the mRNA level, RKIP expression varied significantly across the different molecular subtypes of breast cancer with the Luminal (ER+) subtype expressing high levels of RKIP and the more aggressive Claudin-low (ER-) subtype, which depicted the highest epithelial to mesenchymal transition (EMT) registered the lowest RKIP expression levels. In conclusion, loss of expression/diminution of RKIP or its phosphorylated form is associated with poor diseases-free survival in breast cancer. Determining the expression of RKIP and p-RKIP adds significant prognostic value to the management and subtyping of this disease.

Pitx1 broadly associates with limb enhancers and is enriched on hindlimb cis-regulatory elements.

Extensive functional analyses have demonstrated that the pituitary homeodomain transcription factor Pitx1 plays a critical role in specifying hindlimb morphology in vertebrates. However, much less is known regarding the target genes and cis-regulatory elements through which Pitx1 acts. Earlier studies suggested that the hindlimb transcription factors Tbx4, HoxC10, and HoxC11 might be transcriptional targets of Pitx1, but definitive evidence for direct regulatory interactions has been lacking. Using ChIP-Seq on embryonic mouse hindlimbs, we have pinpointed the genome-wide location of Pitx1 binding sites during mouse hindlimb development and identified potential gene targets for Pitx1. We determined that Pitx1 binding is significantly enriched near genes involved in limb morphogenesis, including Tbx4, HoxC10, and HoxC11. Notably, Pitx1 is bound to the previously identified HLEA and HLEB hindlimb enhancers of the Tbx4 gene and to a newly identified Tbx2 hindlimb enhancer. Moreover, Pitx1 binding is significantly enriched on hindlimb relative to forelimb-specific cis-regulatory features that are differentially marked by H3K27ac. However, our analysis revealed that Pitx1 also strongly associates with many functionally verified limb enhancers that exhibit similar levels of activity in the embryonic mesenchyme of forelimbs and hindlimbs. We speculate that Pitx1 influences hindlimb morphology both through the activation of hindlimb-specific enhancers as well as through the hindlimb-specific modulation of enhancers that are active in both sets of limbs.

Polycomb protein EZH2 regulates tumor invasion via the transcriptional repression of the metastasis suppressor RKIP in breast and prostate cancer.

Epigenetic modifications such as histone methylation play an important role in human cancer metastasis. Enhancer of zeste homolog 2 (EZH2), which encodes the histone methyltransferase component of the polycomb repressive complex 2 (PRC2), is overexpressed widely in breast and prostate cancers and epigenetically silences tumor suppressor genes. Expression levels of the novel tumor and metastasis suppressor Raf-1 kinase inhibitor protein (RKIP) have been shown to correlate negatively with those of EZH2 in breast and prostate cell lines as well as in clinical cancer tissues. Here, we show that the RKIP/EZH2 ratio significantly decreases with the severity of disease and is negatively associated with relapse-free survival in breast cancer. Using a combination of loss- and gain-of-function approaches, we found that EZH2 negatively regulated RKIP transcription through repression-associated histone modifications. Direct recruitment of EZH2 and suppressor of zeste 12 (Suz12) to the proximal E-boxes of the RKIP promoter was accompanied by H3-K27-me3 and H3-K9-me3 modifications. The repressing activity of EZH2 on RKIP expression was dependent on histone deacetylase promoter recruitment and was negatively regulated upstream by miR-101. Together, our findings indicate that EZH2 accelerates cancer cell invasion, in part, via RKIP inhibition. These data also implicate EZH2 in the regulation of RKIP transcription, suggesting a potential mechanism by which EZH2 promotes tumor progression and metastasis.

A new model for raf kinase inhibitory protein induced chemotherapeutic resistance.

Therapeutic resistance remains the most challenging aspect of treating cancer. Raf kinase inhibitory protein (RKIP) emerged as a molecule capable of sensitizing cancerous cells to radio- and chemotherapy. Moreover, this small evolutionary conserved molecule, endows significant resistance to cancer therapy when its expression is reduced or lost. RKIP has been shown to inhibit the Raf-MEK-ERK, NFκB, GRK and activate the GSK3ß signaling pathways. Inhibition of Raf-MEK-ERK and NFκB remains the most prominent pathways implicated in the sensitization of cells to therapeutic drugs. Our purpose was to identify a possible link between RKIP-KEAP 1-NRF2 and drug resistance. To that end, RKIP-KEAP 1 association was tested in human colorectal cancer tissues using immunohistochemistry. RKIP miRNA silencing and its inducible overexpression were employed in HEK-293 immortalized cells, HT29 and HCT116 colon cancer cell lines to further investigate our aim. We show that RKIP enhanced Kelch-like ECH-associated protein1 (KEAP 1) stability in colorectal cancer tissues and HT29 CRC cell line. RKIP silencing in immortalized HEK-293 cells (termed HEK-499) correlated significantly with KEAP 1 protein degradation and subsequent NRF2 addiction in these cells. Moreover, RKIP depletion in HEK-499, compared to control cells, bestowed resistance to supra physiological levels of H(2)O(2) and Cisplatin possibly by upregulating NF-E2-related nuclear factor 2 (NRF2) responsive genes. Similarly, we observed a direct correlation between the extent of apoptosis, after treatment with Adriamycin, and the expression levels of RKIP/KEAP 1 in HT29 but not in HCT116 CRC cells. Our data illuminate, for the first time, the NRF2-KEAP 1 pathway as a possible target for personalized therapeutic intervention in RKIP depleted cancers.

RKIP inhibits NF-kappaB in cancer cells by regulating upstream signaling components of the IkappaB kinase complex.

RKIP was first identified as an inhibitor of the Raf-MEK-ERK signaling pathway. RKIP was also found to play an important role in the NF-kappaB pathway. Genetic and biochemical studies demonstrated that RKIP functioned as a scaffold protein facilitating the phosphorylation of IkappaB by upstream kinases. However, contrary to what one would expect of a scaffold protein, our results show that RKIP has an overall inhibitory effect on the NF-kappaB transcriptional activities. Since NF-kappaB target gene expression is subject to negative regulation involving the optimal induction of negative regulators, our data support a hypothesis that RKIP inhibits NF-kappaB activity via the auto-regulatory feedback loop by rapidly inducing the expression and synthesis of inhibitors of NF-kappaB activation.

The RKIP (Raf-1 Kinase Inhibitor Protein) conserved pocket binds to the phosphorylated N-region of Raf-1 and inhibits the Raf-1-mediated activated phosphorylation of MEK.

The Raf-MEK-ERK pathway regulates many fundamental biological processes, and its activity is finely tuned at multiple levels. The Raf kinase inhibitory protein (RKIP) is a widely expressed negative modulator of the Raf-MEK-ERK signaling pathway. We have previously shown that RKIP inhibits the phosphorylation of MEK by Raf-1 through interfering with the formation of a kinase-substrate complex by direct binding to both Raf-1 and MEK. Here, we show that the evolutionarily conserved ligand-binding pocket of RKIP is required for its inhibitory activity towards the Raf-1 kinase mediated activation of MEK. Single amino acid substitutions of two of the conserved residues form the base and the wall of the pocket confers a loss-of-function phenotype on RKIP. Loss-of-function RKIP mutants still appear to bind to Raf-1. However the stability of the complexes formed between mutants and the N-region Raf-1 phosphopeptide were drastically reduced. Our results therefore suggest that the RKIP conserved pocket may constitute a novel phosphoamino-acid binding motif and is absolutely required for RKIP function.

Regulation of RKIP binding to the N-region of the Raf-1 kinase.

The Raf kinase inhibitory protein (RKIP) binds to Raf-1 interfering with binding of the MEK substrate and potentially also Raf-1 activation. In response to mitogen stimulation RKIP dissociates from Raf-1 and later re-associates. Here, using a combination of mutational approaches, biochemical studies, peptide arrays and plasmon surface resonance (BIAcore), we fine map and characterize a minimal 24 amino acid long RKIP binding domain in the Raf-1 N-region, which consists of constitutive elements at both flanks and a center element that is regulated by phosphorylation and enhances the re-binding of RKIP to Raf-1 in the later phase of mitogen stimulation.

Pyrazolidine derivatives with heteroaryl urea as dipeptidyl peptidase IV inhibitors.

In the continuation of efforts to modify the structure of our novel DP-IV inhibitors, a series of pyrazolidine derivatives with heteroaryl urea was synthesized and evaluated for their ability to inhibit dipeptidyl peptidase IV (DP-IV).