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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

RKIP downregulates B-Raf kinase activity in melanoma cancer cells.

The Raf-MEK-ERK protein kinase cascade is a highly conserved signaling pathway that is pivotal in relaying environmental cues from the cell surface to the nucleus. Three Raf isoforms, which share great sequence and structure similarities, have been identified in mammalian cells. We have previously identified Raf kinase inhibitor protein (RKIP) as a negative regulator of the Raf-MEK-ERK signaling pathway by specifically binding to the Raf-1 isoform. We show here that RKIP also antagonizes kinase activity of the B-Raf isoform. Yeast two-hybrid and coimmunoprecipitation experiments indicated that RKIP specifically interacted with B-Raf. Ectopic expression of RKIP antagonized the kinase activity of B-Raf. We showed that the effects of RKIP on B-Raf functions were independent of its known inhibitory action on Raf-1. The expression levels of RKIP in melanoma cancer cell lines are low relative to primary melanocytes. Forced expression of RKIP partially reverted the oncogenic B-Raf kinase-transformed melanoma cancer cell line SK-Mel-28. The low expression of RKIP and its antagonistic action on B-Raf suggests that RKIP may play an important role in melanoma turmorgenesis.

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.

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.

Inhibition of dipeptidyl peptidase IV by novel inhibitors with pyrazolidine scaffold.

Inhibition of dipeptidyl peptidase IV (DPP-IV) activity has been reported to improve nutrient-stimulated insulin secretion through the stabilization of glucagon-like peptide (GLP-1). In the present study, we identified novel DPP-IV inhibitors of pyrazolidine derivatives (Compounds 1 and 2) and characterized their biological effects in vitro and in vivo. Compound 1, an isoleucine pyrazolidide with a phenyl urea group, inhibited rat plasma DPP-IV, porcine kidney DPP-IV, as well as human Caco-2 DPP-IV with IC(50) values of 1.70, 2.26, and 2.02 microM, respectively. Because of the poor pharmacokinetic properties of Compound 1, further optimization was carried out, leading to the discovery of Compound 2, which had similar in vitro activities. Compound 2 acted as a selective and competitive inhibitor of DPP-IV. MALDI-TOF mass spectrometric analysis proved that the compound (20 microM) effectively blocked the degradation of active GLP-1 peptide by 61%. Although similar in in vitro potency, marked improvement of in vivo efficacy and pharmacokinetic properties was seen with Compound 2. Oral administration of Compound 2 resulted in potent and rapid inhibition of circulating DPP-IV in C57BL/6J mice, with ED(50) values of 26mg/kg (s.c.) and 42mg/kg (p.o.). In addition, this compound improved glucose tolerance in ob/ob mice, as determined by an oral glucose tolerance test (OGTT). These results indicate that Compound 2 is a potent and selective DPP-IV inhibitor with oral anti-hyperglycemic activity in vivo.

Two-component signaling in the AAA + ATPase DctD: binding Mg2+ and BeF3- selects between alternate dimeric states of the receiver domain.

A Crystallogral structure is described for the Mg2+-BeF3--bound receiver domain of Sinorhizobium meliloti DctD bearing amino acid substitution E121K. Differences between the apo- and ligand-bound active sites are similar to those reported for other receiver domains. However, the off and on states of the DctD receiver domain are characterized by dramatically different dimeric structures, which supports the following hypothesis of signal transduction. In the off state, the receiver domain and coiled-coil linker form a dimer that inhibits oligomerization of the AAA+ ATPase domain. In this conformation, the receiver domain cannot be phosphorylated or bind Mg2+ and BeF3-. Instead, these modifications stabilize an alternative dimeric conformation that repositions the subunits by approximately 20 A, thus replacing the a4-b5-a5 interface with an a4-b5 interface. Reoriented receiver domains permit the ATPase domain to oligomerize and stimulate open complex formation by the s54 form of RNA polymerase. NtrC, which shares 38% sequence identity with DctD, works differently. Its activated receiver domain must facilitate oligomerization of its ATPase domain. Significant differences exist in the signaling surfaces of the DctD and NtrC receiver domains that may help explain how triggering the common two-component switch can variously regulate assembly of a AAA+ ATPase domain.

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.

Engineering the serine/threonine protein kinase Raf-1 to utilise an orthogonal analogue of ATP substituted at the N6 position.

One key area of protein kinase research is the identification of cognate substrates. The search for substrates is hampered by problems in unambiguously assigning substrates to a particular kinase in vitro and in vivo. One solution to this impasse is to engineer the kinase of interest to accept an ATP analogue which is orthogonal (unable to fit into the ATP binding site) for the wild-type enzyme and the majority of other kinases. The acceptance of structurally modified, gamma-(32)P-labelled, nucleotide analogue by active site-modified kinase can provide a unique handle by which the direct substrates of any particular kinase can be displayed in crude mixtures or cell lysates. We have taken this approach with the serine/threonine kinase Raf-1, which plays an essential role in the transduction of stimuli through the Ras-->Raf-->MEK-->ERK/MAP kinase cascade. This cascade plays essential roles in proliferation, differentiation and apoptosis. Here we detail the mutagenesis strategy for the ATP binding pocket of Raf-1, such that it can utilise an N(6)-substituted ATP analogue. We show that these mutations do not alter the substrate specificity and signal transduction through Raf-1. We screen a library of analogues to identify which are orthogonal for Raf-1, and show that mutant Raf-1 can utilise the orthogonal analogue N(6)(2-phenethyl) ATP in vitro to phosphorylate its currently only accepted substrate MEK. Importantly we show that our approach can be used to tag putative direct substrates of Raf-1 kinase with (32)P-N(6)(2-phenethyl) ATP in cell lysates.

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.