Functional conservation of a forebrain enhancer from the elephant shark (Callorhinchus milii) in zebrafish and mice

Background: The phylogenetic position of the elephant shark (Callorhinchus milii ) is particularly relevant to study the evolution of genes and gene regulation in vertebrates. Here we examine the evolution of Dlx homeobox gene regulation during vertebrate embryonic development with a particular focus on the forebrain. We first identified the elephant shark sequence orthologous to the URE2 cis -regulatory element of the mouse Dlx1/Dlx2 locus (herein named CmURE2). We then conducted a comparative study of the sequence and enhancer activity of CmURE2 with that of orthologous regulatory sequences from zebrafish and mouse. Results: The CmURE2 sequence shows a high percentage of identity with its mouse and zebrafish counterparts but is overall more similar to mouse URE2 (MmURE2) than to zebrafish URE2 (DrURE2). In transgenic zebrafish and mouse embryos, CmURE2 displayed enhancer activity in the forebrain that overlapped with that of DrURE2 and MmURE2. However, we detected notable differences in the activity of the three sequences in the diencephalon. Outside of the forebrain, CmURE2 shows enhancer activity in areas such as the pharyngeal arches and dorsal root ganglia where its' counterparts are also active. Conclusions: Our transgenic assays show that part of the URE2 enhancer activity is conserved throughout jawed vertebrates but also that new characteristics have evolved in the different groups. Our study demonstrates that the elephant shark is a useful outgroup to study the evolution of regulatory mechanisms in vertebrates and to address how changes in the sequence of cis -regulatory elements translate into changes in their regulatory activity

AMOD: a morpholino oligonucleotide selection tool

AMOD is a web-based program that aids in the functional evaluation of nucleotide sequences through sequence characterization and antisense morpholino oligonucleotide (target site) selection. Submitted sequences are analyzed by translation initiation site prediction algorithms and sequence-to-sequence comparisons; results are used to characterize sequence features required for morpholino design. Within a defined subsequence, base composition and homodimerization values are computed for all putative morpholino oligonucleotides. Using these properties, morpholino candidates are selected and compared with genomic and transcriptome databases with the goal to identify target-specific enriched morpholinos. AMOD has been used at the University of Minnesota to design ∼200 morpholinos for a functional genomics screen in zebrafish. The AMOD web server and a tutorial are freely available to both academic and commercial users at

Analysis of four DLX homeobox genes in autistic probands

BACKGROUND: Linkage studies in autism have identified susceptibility loci on chromosomes 2q and 7q, regions containing the DLX1/2 and DLX5/6 bigene clusters. The DLX genes encode homeodomain transcription factors that control craniofacial patterning and differentiation and survival of forebrain inhibitory neurons. We investigated the role that sequence variants in DLX genes play in autism by in-depth resequencing of these genes in 161 autism probands from the AGRE collection. RESULTS: Sequencing of exons, exon/intron boundaries and known enhancers of DLX1, 2, 5 and 6 identified several nonsynonymous variants in DLX2 and DLX5 and a variant in a DLX5/6intragenic enhancer. The nonsynonymous variants were detected in 4 of 95 families from which samples were sequenced. Two of these four SNPs were not observed in 378 undiagnosed samples from North American populations, while the remaining 2 were seen in one sample each. CONCLUSION: Segregation of these variants in pedigrees did not generally support a contribution to autism susceptibility by these genes, although functional analyses may provide insight into the biological understanding of these important proteins

zfishbook: connecting you to a world of zebrafish revertible mutants

zfishbook is an internet-based openly accessible database of revertible protein trap gene-breaking transposon (GBT) insertional mutants in the zebrafish, Danio rerio. In these lines, a monomeric red fluorescent protein (mRFP) is encoded by an artificial 3′ exon, resulting in a translational fusion to endogenous loci. The natural transparency of the zebrafish embryo and larvae greatly facilitates the expression annotation of tagged loci using new capillary-based SCORE imaging methods. Molecular annotation of each line is facilitated by cloning methods such as 5′-Rapid Amplification of cDNA Ends (RACE) and inverse polymerase chain reaction (PCR). zfishbook (http://zfishbook.org) represents a central hub for molecular, expression and mutational information about GBT lines from the International Zebrafish Protein Trap Consortium (IZPTC) that includes researchers from around the globe. zfishbook is open to community-wide contributions including expression and functional annotation. zfishbook also represents a central location for information on how to obtain these lines from diverse members of the IZPTC and integration within other zebrafish community databases including Zebrafish Information Network (ZFIN), Ensembl and National Center for Biotechnology Information

Functional Analysis of Human Hematopoietic Stem Cell Gene Expression Using Zebrafish

Although several reports have characterized the hematopoietic stem cell (HSC) transcriptome, the roles of HSC-specific genes in hematopoiesis remain elusive. To identify candidate regulators of HSC fate decisions, we compared the transcriptome of human umbilical cord blood and bone marrow CD34(+)CD33(−)CD38(−)Rho(lo)c-kit(+) cells, enriched for hematopoietic stem/progenitor cells with CD34(+)CD33(−)CD38(−)Rho(hi) cells, enriched in committed progenitors. We identified 277 differentially expressed transcripts conserved in these ontogenically distinct cell sources. We next performed a morpholino antisense oligonucleotide (MO)-based functional screen in zebrafish to determine the hematopoietic function of 61 genes that had no previously known function in HSC biology and for which a likely zebrafish ortholog could be identified. MO knock down of 14/61 (23%) of the differentially expressed transcripts resulted in hematopoietic defects in developing zebrafish embryos, as demonstrated by altered levels of circulating blood cells at 30 and 48 h postfertilization and subsequently confirmed by quantitative RT-PCR for erythroid-specific hbae1 and myeloid-specific lcp1 transcripts. Recapitulating the knockdown phenotype using a second MO of independent sequence, absence of the phenotype using a mismatched MO sequence, and rescue of the phenotype by cDNA-based overexpression of the targeted transcript for zebrafish spry4 confirmed the specificity of MO targeting in this system. Further characterization of the spry4-deficient zebrafish embryos demonstrated that hematopoietic defects were not due to more widespread defects in the mesodermal development, and therefore represented primary defects in HSC specification, proliferation, and/or differentiation. Overall, this high-throughput screen for the functional validation of differentially expressed genes using a zebrafish model of hematopoiesis represents a major step toward obtaining meaningful information from global gene profiling of HSCs

Predicting the binding preference of transcription factors to individual DNA k-mers

Motivation: Recognition of specific DNA sequences is a central mechanism by which transcription factors (TFs) control gene expression. Many TF-binding preferences, however, are unknown or poorly characterized, in part due to the difficulty associated with determining their specificity experimentally, and an incomplete understanding of the mechanisms governing sequence specificity. New techniques that estimate the affinity of TFs to all possible k-mers provide a new opportunity to study DNA–protein interaction mechanisms, and may facilitate inference of binding preferences for members of a given TF family when such information is available for other family members. Results: We employed a new dataset consisting of the relative preferences of mouse homeodomains for all eight-base DNA sequences in order to ask how well we can predict the binding profiles of homeodomains when only their protein sequences are given. We evaluated a panel of standard statistical inference techniques, as well as variations of the protein features considered. Nearest neighbour among functionally important residues emerged among the most effective methods. Our results underscore the complexity of TF–DNA recognition, and suggest a rational approach for future analyses of TF families. Contact: [email protected] Supplementary information: Supplementary data are available at Bioinformatics online.Canadian Institutes of Health ResearchOntario Research FundNational Institutes of Health (U.S.)National Human Genome Research Institute (U.S.

Recurrent Modification of a Conserved Cis-Regulatory Element Underlies Fruit Fly Pigmentation Diversity

The development of morphological traits occurs through the collective action of networks of genes connected at the level of gene expression. As any node in a network may be a target of evolutionary change, the recurrent targeting of the same node would indicate that the path of evolution is biased for the relevant trait and network. Although examples of parallel evolution have implicated recurrent modification of the same gene and cis-regulatory element (CRE), little is known about the mutational and molecular paths of parallel CRE evolution. In Drosophila melanogaster fruit flies, the Bric-à-brac (Bab) transcription factors control the development of a suite of sexually dimorphic traits on the posterior abdomen. Female-specific Bab expression is regulated by the dimorphic element, a CRE that possesses direct inputs from body plan (ABD-B) and sex-determination (DSX) transcription factors. Here, we find that the recurrent evolutionary modification of this CRE underlies both intraspecific and interspecific variation in female pigmentation in the melanogaster species group. By reconstructing the sequence and regulatory activity of the ancestral Drosophila melanogaster dimorphic element, we demonstrate that a handful of mutations were sufficient to create independent CRE alleles with differing activities. Moreover, intraspecific and interspecific dimorphic element evolution proceeded with little to no alterations to the known body plan and sex-determination regulatory linkages. Collectively, our findings represent an example where the paths of evolution appear biased to a specific CRE, and drastic changes in function were accompanied by deep conservation of key regulatory linkages. © 2013 Rogers et al

Recognition models to predict DNA-binding specificities of homeodomain proteins

Motivation: Recognition models for protein-DNA interactions, which allow the prediction of specificity for a DNA-binding domain based only on its sequence or the alteration of specificity through rational design, have long been a goal of computational biology. There has been some progress in constructing useful models, especially for C2H2 zinc finger proteins, but it remains a challenging problem with ample room for improvement. For most families of transcription factors the best available methods utilize k-nearest neighbor (KNN) algorithms to make specificity predictions based on the average of the specificities of the k most similar proteins with defined specificities. Homeodomain (HD) proteins are the second most abundant family of transcription factors, after zinc fingers, in most metazoan genomes, and as a consequence an effective recognition model for this family would facilitate predictive models of many transcriptional regulatory networks within these genomes