Search

Transcriptional regulation of gene expression plays a significant role in establishing the diversity of human cell types and biological functions from a common set of genes. The components of regulatory control in the human genome include cis-acting elements that act across immense genomic distances to influence the spatial and temporal distribution of gene expression. Here we review the established categories of distant-acting regulatory elements, discussing the classical and contemporary evidence of their regulatory potential and clinical importance. Current efforts to identify regulatory sequences throughout the genome and elucidate their biological significance depend heavily on advances in sequence conservation-based analyses and on increasingly large-scale efforts applying transgenic technologies in model organisms. We discuss the advantages and limitations of sequence conservation as a predictor of regulatory function and present complementary emerging technologies now being applied to annotate regulatory elements in vertebrate genomes. [ABSTRACT FROM AUTHOR]

Abstract: Zebrafish transgenesis is a powerful and increasingly common strategy to assay vertebrate transcriptional regulatory control. Several challenges remain, however, to the broader application of this technique; they include increasing the rate with which transgenes can be analyzed and maximizing the informational value of the data generated. Presently, many rely on the injection of individual constructs and the analysis of resulting reporter expression in mosaic G0 embryos. Here, we contrast these approaches, examining whether injecting pooled transgene constructs can increase the efficiency with which regulatory sequences can be assayed, restricting analysis to the offspring of germ line transmitting transgenic zebrafish in an effort to reduce potential subjectivity. We selected a 64kb interval encompassing the human ASCL1 locus as our model interval and report the analysis of 9 highly conserved putative enhancers therein. We identified 32 transgene-positive zebrafish, transmitting one or more independent constructs displaying ASCL1-like regulatory control. Through examination of embryos harboring one or more transgenes, we demonstrate that five of the nine sequences account for the observed control and describe their likely roles in ASCL1 regulation. These data demonstrate the utility of this approach and its potential for further adaptation and higher throughput application. [Copyright &y& Elsevier]

The identification of common variants that contribute to the genesis of human inherited disorders remains a significant challenge. Hirschsprung disease (HSCR) is a multifactorial, non-mendelian disorder in which rare high-penetrance coding sequence mutations in the receptor tyrosine kinase RET contribute to risk in combination with mutations at other genes. We have used family-based association studies to identify a disease interval, and integrated this with comparative and functional genomic analysis to prioritize conserved and functional elements within which mutations can be sought. We now show that a common non-coding RET variant within a conserved enhancer-like sequence in intron 1 is significantly associated with HSCR susceptibility and makes a 20-fold greater contribution to risk than rare alleles do. This mutation reduces in vitro enhancer activity markedly, has low penetrance, has different genetic effects in males and females, and explains several features of the complex inheritance pattern of HSCR. Thus, common low-penetrance variants, identified by association studies, can underlie both common and rare diseases. [ABSTRACT FROM AUTHOR]

Clinical expression of Hirschsprung disease (HSCR) requires the interaction of multiple susceptibility genes. Molecular genetic analyses have revealed that interactions between mutations in the genes encoding the RET receptor tyrosine kinase and the endothelin receptor type B (EDNRB) are central to the genesis of HSCR. We have established two locus noncomplementation assays in mice, using allelic series at Ednrb in the context of Ret kinase-null heterozygotes, to understand the clinical presentation, incomplete penetrance, variation in length of aganglionic segment, and sex bias observed in human HSCR patients. Titration of Ednrb in the presence of half the genetic dose of Ret determines the presentation of an enteric phenotype in these strains, revealing or abrogating a sex bias in disease expression depending on the genotype at Ednrb. RET and EDNRB signaling pathways are also critical for the normal development of other tissues, including the kidneys and neural crest-derived melanocytes. Our data demonstrate that interaction between these genes is restricted to the enteric nervous system and does not affect renal, coat color, and retinal choroid development. [ABSTRACT FROM AUTHOR]

Genetic studies of Hirschsprung disease, a common congenital malformation, have identified eight genes with mutations that can be associated with this condition. Mutations at individual loci are, however, neither necessary nor sufficient to cause clinical disease. We conducted a genome-wide association study in 43 Mennonite family trios using 2,083 microsatellites and single-nucleotide polymorphisms and a new multipoint linkage disequilibrium method that searches for association arising from common ancestry. We identified susceptibility loci at 10q11, 13q22 and 16q23; the gene at 13q22 is EDNRB, encoding a G proteincoupled receptor (GPCR) and the gene at 10qll is RET, encoding a receptor tyrosine kinase (RTK). Statistically significant joint transmission of RET and EDNRB alleles in affected individuals and non-complementation of aganglionosis in mouse intercrosses between Ret null and the Ednrb hypomorphic piebald allele are suggestive of epistasis between EDNRB and RET. Thus, genetic interaction between mutations in RET and EDNRB is an underlying mechanism for this complex disorder. [ABSTRACT FROM AUTHOR]

Determines the role of Endothelin type-B (Ednrb) receptor and Endothelin 3 in the normal development of melanocytes and crest-derived lineages, with emphasis on type II Hirschsprung (HSCR) disease. Genetics of HSCR type II disease; Relationships between Ednrb, melanocytes and the regulation of genes involved in pigment cell development; Endothelins and the development of neural crest.

Simple Summary: Ring chromosomes (RCs) are formed when chromosome ends fuse to form a circular structure. RCs are seen in <10% of patients with blood cancers (i.e., leukemias, lymphomas, and myelomas), and are associated with poor disease outcomes. We sought to evaluate how frequently RCs arise in patients diagnosed with these cancers and determine if there are any associated genetic variants, regions, or chromosomal abnormalities. Most patients with RCs have mutations in the TP53 gene and gross chromosomal rearrangements, particularly among chromosomes 5, 7, 11, and 17. Most of these patients also possess marker chromosomes with genetic material of unknown origin. The majority of RCs that resemble identifiable chromosomes, as well as chromosomes lacking gross rearrangements, were seen in patients without TP53 mutations. Our data suggest that mechanisms underlying RC generation may arise differently in patient groups with the presence or absence of functional TP53, a potentially important distinction for clinical decision making. Ring chromosomes (RC) are present in <10% of patients with hematological malignancies and are associated with poor prognosis. Until now, only small cohorts of patients with hematological neoplasms and concomitant RCs have been cytogenetically characterized. Here, we performed a conventional chromosome analysis on metaphase spreads from >13,000 patients diagnosed with hematological malignancies at the Johns Hopkins University Hospital and identified 98 patients with RCs—90 with myeloid malignancies and 8 with lymphoid malignancies. We also performed a targeted Next-Generation Sequencing (NGS) assay, using a panel of 642 cancer genes, to identify whether these patients harbor relevant pathogenic variants. Cytogenetic analyses revealed that RCs and marker chromosomes of unknown origin are concurrently present in most patients by karyotyping, and 93% of patients with NGS data have complex karyotypes. A total of 72% of these individuals have pathogenic mutations in TP53, most of whom also possess cytogenetic abnormalities resulting in the loss of 17p, including the loss of TP53. All patients with a detected RC and without complex karyotypes also lack TP53 mutations but have pathogenic mutations in TET2. Further, 70% of RCs that map to a known chromosome are detected in individuals without TP53 mutations. Our data suggest that RCs in hematological malignancies may arise through different mechanisms, but ultimately promote widespread chromosomal instability. [ABSTRACT FROM AUTHOR]

Haploid cell lines are a valuable research tool with broad applicability for genetic assays. As such the fully haploid human cell line, eHAP1, has been used in a wide array of studies. However, the absence of a corresponding reference genome sequence for this cell line has limited the potential for more widespread applications to experiments dependent on available sequence, like capture-clone methodologies. We generated ~15× coverage Nanopore long reads from ten GridION flowcells and utilized this data to assemble a de novo draft genome using minimap and miniasm and subsequently polished using Racon. This assembly was further polished using previously generated, low-coverage, Illumina short reads with Pilon and ntEdit. This resulted in a hybrid eHAP1 assembly with >90% complete BUSCO scores. We further assessed the eHAP1 long read data for structural variants using Sniffles and identify a variety of rearrangements, including a previously established Philadelphia translocation. Finally, we demonstrate how some of these variants overlap open chromatin regions, potentially impacting regulatory regions. By integrating both long and short reads, we generated a high-quality reference assembly for eHAP1 cells. The union of long and short reads demonstrates the utility in combining sequencing platforms to generate a high-quality reference genome de novo solely from low coverage data. We expect the resulting eHAP1 genome assembly to provide a useful resource to enable novel experimental applications in this important model cell line. • Generation of a reference genome for the human haploid cell line, eHAP1 • Long but noisy Nanopore sequencing combined with accurate but short Illumina sequencing yields a polished reference genome • Structural variants, including the Philadelphia translocation, can be identified using exclusively long reads [ABSTRACT FROM AUTHOR]

To overcome the ethical and technical limitations of in vivo human disease models, the broader scientific community frequently employs model organism-derived cell lines to investigate disease mechanisms, pathways, and therapeutic strategies. Despite the widespread use of certain in vitro models, many still lack contemporary genomic analysis supporting their use as a proxy for the affected human cells and tissues. Consequently, it is imperative to determine how accurately and effectively any proposed biological surrogate may reflect the biological processes it is assumed to model. One such cellular surrogate of human disease is the established mouse neural precursor cell line, SN4741, which has been used to elucidate mechanisms of neurotoxicity in Parkinson disease for over 25 years. Here, we are using a combination of classic and contemporary genomic techniques – karyotyping, RT-qPCR, single cell RNA-seq, bulk RNA-seq, and ATAC-seq – to characterize the transcriptional landscape, chromatin landscape, and genomic architecture of this cell line, and evaluate its suitability as a proxy for midbrain dopaminergic neurons in the study of Parkinson disease. We find that SN4741 cells possess an unstable triploidy and consistently exhibits low expression of dopaminergic neuron markers across assays, even when the cell line is shifted to the non-permissive temperature that drives differentiation. The transcriptional signatures of SN4741 cells suggest that they are maintained in an undifferentiated state at the permissive temperature and differentiate into immature neurons at the non-permissive temperature; however, they may not be dopaminergic neuron precursors, as previously suggested. Additionally, the chromatin landscapes of SN4741 cells, in both the differentiated and undifferentiated states, are not concordant with the open chromatin profiles of ex vivo, mouse E15.5 forebrain- or midbrain-derived dopaminergic neurons. Overall, our data suggest that SN4741 cells may reflect early aspects of neuronal differentiation but are likely not a suitable proxy for dopaminergic neurons as previously thought. The implications of this study extend broadly, illuminating the need for robust biological and genomic rationale underpinning the use of in vitro models of molecular processes. [ABSTRACT FROM AUTHOR]

Multifactorial diseases are characterized by inter-individual variation in etiology, age of onset, and penetrance. These diseases tend to be relatively common and arise from the combined action of genetic and environmental factors; however, parsing the convoluted mechanisms underlying these gene-by-environment interactions presents a significant challenge to their study and management. For neurodegenerative disorders, resolving this challenge is imperative, given the enormous health and societal burdens they impose. The mechanisms by which genetic and environmental effects may act in concert to destabilize homeostasis and elevate risk has become a major research focus in the study of common disease. Emphasis is further being placed on determining the extent to which a unifying biological principle may account for the progressively diminishing capacity of a system to buffer disease phenotypes, as risk for disease increases. Data emerging from studies of common, neurodegenerative diseases are providing insights to pragmatically connect mechanisms of genetic and environmental risk that previously seemed disparate. In this review, we discuss evidence positing inflammation as a unifying biological principle of homeostatic destabilization affecting the risk, onset, and progression of neurodegenerative diseases. Specifically, we discuss how genetic variation associated with Alzheimer disease and Parkinson disease may contribute to pro-inflammatory responses, how such underlying predisposition may be exacerbated by environmental insults, and how this common theme is being leveraged in the ongoing search for effective therapeutic interventions. [ABSTRACT FROM AUTHOR]

The microtubule-associated protein tau (MAPT) gene is considered a strong genetic risk factor for Parkinson's disease (PD) in Caucasians. MAPT is located within an inversion region of high linkage disequilibrium designated as H1 and H2 haplotype, and contains eight other genes which have been implicated in neurodegeneration. The aim of the current study was to identify common coding variants in strong linkage disequilibrium (LD) within the associated loci on chr17q21 harboring MAPT. Sanger sequencing of coding exons in 90 Caucasian late-onset PD (LOPD) patients was performed. Specific gene sequencing for LRRC37A , LRRC37A2 , ARL17A and ARL17B was not possible given the high homology, presence of pseudogenes and copy number variants that are in the region, and therefore four genes (NSF , KANSL1 , SPPL2C , and CRHR1) were included in the analysis. Coding variants from these four genes that did not perfectly tag (r2 = 1) the MAPT H1/H2 haplotype were genotyped in an independent replication series of Caucasian PD cases (N = 851) and controls (N = 730). In the 90 LOPD cases we identified 30 coding variants. Eleven non-synonymous variants tagged the MAPT H1/H2 haplotype, including two SPPL2C variants (rs12185233 and rs12373123) that had high pathogenic combined annotation dependent depletion (CADD) scores of >20. In the replication series, the non-synonymous KANSL1 rs17585974 variant was in very strong LD with MAPT H1/H2 and had a high CADD score of 24.7. We have identified several non-synonymous variants across neighboring genes of MAPT that may warrant further genetic and functional investigation within the biological etiology of PD. • MAPT H1 is a consistent Parkinson's disease association locus on chr17q21. • Chr17q21H1 extended haplotype is a complex genetic region of inversion. • Other genes/variants in complete linkage disequilibrium may be responsible. • Potentially damaging variants in KANSL1 , NSF and SPPL2C. [ABSTRACT FROM AUTHOR]

Introduction: The microtubule-associated protein tau (MAPT) gene is considered a strong genetic risk factor for Parkinson's disease (PD) in Caucasians. MAPT is located within an inversion region of high linkage disequilibrium designated as H1 and H2 haplotype, and contains eight other genes which have been implicated in neurodegeneration. The aim of the current study was to identify common coding variants in strong linkage disequilibrium (LD) within the associated loci on chr17q21 harboring MAPT.Methods: Sanger sequencing of coding exons in 90 Caucasian late-onset PD (LOPD) patients was performed. Specific gene sequencing for LRRC37A, LRRC37A2, ARL17A and ARL17B was not possible given the high homology, presence of pseudogenes and copy number variants that are in the region, and therefore four genes (NSF, KANSL1, SPPL2C, and CRHR1) were included in the analysis. Coding variants from these four genes that did not perfectly tag (r2 = 1) the MAPT H1/H2 haplotype were genotyped in an independent replication series of Caucasian PD cases (N = 851) and controls (N = 730).Results: In the 90 LOPD cases we identified 30 coding variants. Eleven non-synonymous variants tagged the MAPT H1/H2 haplotype, including two SPPL2C variants (rs12185233 and rs12373123) that had high pathogenic combined annotation dependent depletion (CADD) scores of >20. In the replication series, the non-synonymous KANSL1 rs17585974 variant was in very strong LD with MAPT H1/H2 and had a high CADD score of 24.7.Conclusion: We have identified several non-synonymous variants across neighboring genes of MAPT that may warrant further genetic and functional investigation within the biological etiology of PD. [ABSTRACT FROM AUTHOR]

Trisomy for human chromosome 21 (Hsa21) results in Down syndrome (DS), one of the most genetically complex conditions compatible with human survival. Assessment of the physiological consequences of dosage-driven overexpression of individual Hsa21 genes during early embryogenesis and the resulting contributions to DS pathology in mammals are not tractable in a systematic way. A recent study looked at loss-of-function of a subset of Caenorhabditis elegans orthologs of Hsa21 genes and identified ten candidates with behavioral phenotypes, but the equivalent over-expression experiment has not been done. We turned to zebrafish as a developmental model and, using a number of surrogate phenotypes, we screened Hsa21 genes for effects on early embyrogenesis. We prepared a library of 164 cDNAs of conserved protein coding genes, injected mRNA into early embryos and evaluated up to 5 days post-fertilization (dpf). Twenty-four genes produced a gross morphological phenotype, 11 of which could be reproduced reliably. Seven of these gave a phenotype consistent with down regulation of the sonic hedgehog (Shh) pathway; two showed defects indicative of defective neural crest migration; one resulted consistently in pericardial edema; and one was embryonic lethal. Combinatorial injections of multiple Hsa21 genes revealed both additive and compensatory effects, supporting the notion that complex genetic relationships underlie end phenotypes of trisomy that produce DS. Together, our data suggest that this system is useful in the genetic dissection of dosage-sensitive gene effects on early development and can inform the contribution of both individual loci and their combinatorial effects to phenotypes relevant to the etiopathology of DS. [ABSTRACT FROM AUTHOR]

Calcium channel subunits, including CACNA1C, have been associated with multiple psychiatric disorders. Specifically, genome wide association studies (GWAS) have repeatedly identified the single nucleotide polymorphism (SNP) rs1006737 in intron 3 of CACNA1C to be strongly associated with schizophrenia and bipolar disorder. Here, we show that rs1006737 marks a quantitative trait locus for CACNA1C transcript levels. We test 16 SNPs in high linkage disequilibrium with rs1007637 and find one, rs4765905, consistently showing allele-dependent regulatory function in reporter assays. We find allele-specific protein binding for 13 SNPs including rs4765905. Using protein microarrays, we identify several proteins binding ≥3 SNPs, but not control sequences, suggesting possible functional interactions and combinatorial haplotype effects. Finally, using circular chromatin conformation capture, we show interaction of the disease-associated region including the 16 SNPs with the CACNA1C promoter and other potential regulatory regions. Our results elucidate the pathogenic relevance of one of the best-supported risk loci for schizophrenia and bipolar disorder. [ABSTRACT FROM AUTHOR]

Most variants implicated in common human disease by genome-wide association studies (GWAS) lie in noncoding sequence intervals. Despite the suggestion that regulatory element disruption represents a common theme, identifying causal risk variants within implicated genomic regions remains a major challenge. Here we present a new sequence-based computational method to predict the effect of regulatory variation, using a classifier (gkm-SVM) that encodes cell type-specific regulatory sequence vocabularies. The induced change in the gkm-SVM score, deltaSVM, quantifies the effect of variants. We show that deltaSVM accurately predicts the impact of SNPs on DNase I sensitivity in their native genomic contexts and accurately predicts the results of dense mutagenesis of several enhancers in reporter assays. Previously validated GWAS SNPs yield large deltaSVM scores, and we predict new risk-conferring SNPs for several autoimmune diseases. Thus, deltaSVM provides a powerful computational approach to systematically identify functional regulatory variants. [ABSTRACT FROM AUTHOR]

Pyridostigmine bromide (PB) is an FDA-approved drug for the treatment of myasthenia gravis and a prophylactic pre-treatment for organophosphate nerve agent poisoning. Current methods for evaluating nerve agent treatments include enzymatic studies and mammalian models. Rapid whole animal screening tools for assessing the effects of nerve agent pre-treatment and post-exposure drugs represent an underdeveloped area of research. We used zebrafish as a model for acute and chronic developmental exposure to PB and two related carbamate acetylcholinesterase (AChE) inhibitors, neostigmine bromide (NB) and physostigmine (PS). Lethal doses and gross morphological phenotypes resulting from exposure to sub-lethal doses of these compounds were determined. Quantitative analyses of motility impairment and AChE enzyme inhibition were used to determine optimal dosing conditions for evaluation of the effects of carbamate exposures on neuronal development; ~ 50% impairment of response to startle stimuli and > 50% inhibition of AChE activity were observed at 80 mM PB, 20 mM NB and 0.1 mM PS. PB induced stunted somite length, but no other phenotypic effects were observed. In contrast, NB and PS induced more severe phenotypic morphological defects than PB as well as neurite outgrowth mislocalization. Additionally, NB induced mislocalization of nicotinic acetylcholine receptors, resulting in impaired synapse formation. Taken together, these data suggest that altered patterns of neuronal connectivity contribute to the developmental neurotoxicity of carbamates and demonstrate the utility of the zebrafish model for distinguishing subtle structure-based differential effects of AChE inhibitors, which include nerve agents, pesticides and drugs. [ABSTRACT FROM AUTHOR]

Glutaminase (GLS), which converts glutamine to glutamate, plays a key role in cancer cell metabolism, growth, and proliferation. GLS is being explored as a cancer therapeutic target, but whether GLS inhibitors affect cancer cell--autonomous growth or the host microenvironment or have off-target effects is unknown. Here, we report that loss of one copy of Gls blunted tumor progression in an immune-competent MYC-mediated mouse model of hepatocellular carcinoma. Compared with results in untreated animals with MYC-induced hepatocellular carcinoma, administration of the GLS-specific inhibitor bis-2-(5-phenylacetamido-1,3,4-thiadiazol-2-yl)ethyl sulfide (BPTES) prolonged survival without any apparent toxicities. BPTES also inhibited growth of a MYC-dependent human B cell lymphoma cell line (P493) by blocking DNA replication, leading to cell death and fragmentation. In mice harboring P493 tumor xenografts, BPTES treatment inhibited tumor cell growth; however, P493 xenografts expressing a BPTES-resistant GLS mutant (GLS-K325A) or overexpressing GLS were not affected by BPTES treatment. Moreover, a customized Vivo-Morpholino that targets human GLS mRNA markedly inhibited P493 xenograft growth without affecting mouse Gls expression. Conversely, a Vivo-Morpholino directed at mouse Gls had no antitumor activity in vivo. Collectively, our studies demonstrate that GLS is required for tumorigenesis and support small molecule and genetic inhibition of GLS as potential approaches for targeting the tumor cell--autonomous dependence on GLS for cancer therapy. [ABSTRACT FROM AUTHOR]

We recently demonstrated that the gene encoding the RNA binding motif protein 24 (RBM24) is expressed during mouse cardiogenesis, and determined the developmental requirement for its zebrafish homologs Rbm24a and Rbm24b during cardiac development. We demonstrate here that both Rbm24a and Rbm24b are also required for normal somite and craniofacial development. Diminution of rbm24a or rbm24b gene products by morpholino knockdown resulted in significant disruption of somite formation. Detailed in situ hybridization-based analyses of a spectrum of somitogenesis-associated transcripts revealed reduced expression of the cyclic muscle pattering genes dlc and dld encoding Notch ligands, as well as their respective target genes her7, her1. By contrast expression of the Notch receptors notch1a and notch3 appears unchanged. Some RBM-family members have been implicated in pre-mRNA processing. Analysis of affected Notch-pathway mRNAs in rbm24a and rbm24b morpholino-injected embryos revealed aberrant transcript fragments of dlc and dld, but not her1 or her7, suggesting the reduction in transcription levels of Notch pathway components may result from aberrant processing of its ligands. These data imply a previously unknown requirement for Rbm24a and Rbm24b in somite and craniofacial development. Although we anticipate the influence of disrupting RBM24 homologs likely extends beyond the Notch pathway, our results suggest their perturbation may directly, or indirectly, compromise post-transcriptional processing, exemplified by imprecise processing of dlc and dld. [ABSTRACT FROM AUTHOR]

Abstract: In this study, we sought to establish a novel method to prospectively and dynamically identify live human oligodendrocyte precursor cells (OPCs) and oligodendrocyte lineage cells from brain dissociates and pluripotent stem cell culture. We selected a highly conserved enhancer element of the Sox10 gene, known as MCS5, which directs reporter expression to oligodendrocyte lineage cells in mouse and zebrafish. We demonstrate that lentiviral Sox10-MCS5 induced expression of GFP at high levels in a subpopulation of human CD140a/PDGFαR-sorted OPCs as well as their immature oligodendrocyte progeny. Furthermore, we show that almost all Sox10-MCS5:GFPhigh cells expressed OPC antigen CD140a and human OPCs expressing SOX10, OLIG2, and PDGFRA mRNAs could be prospectively identified using GFP based fluorescence activated cells sorting alone. Additionally, we established a human induced pluripotent cell (iPSC) line transduced with the Sox10-MCS5:GFP reporter using a Rex-Neo cassette. Similar to human primary cells, GFP expression was restricted to embryoid bodies containing both oligodendrocyte progenitor and oligodendrocyte cells and co-localized with NG2 and O4-positive cells respectively. As such, we have developed a novel reporter system that can track oligodendrocyte commitment in human cells, establishing a valuable tool to improve our understanding and efficiency of human oligodendrocyte derivation. [Copyright &y& Elsevier]

We take a comprehensive approach to the study of regulatory control of gene expression in melanocytes that proceeds from large-scale enhancer discovery facilitated by ChIP-seq; to rigorous validation in silico, in vitro, and in vivo; and finally to the use of machine learning to elucidate a regulatory vocabulary with genome-wide predictive power. We identify 2489 putative melanocyte enhancer loci in the mouse genome by ChIP-seq for EP300 and H3K4me1. We demonstrate that these putative enhancers are evolutionarily constrained, enriched for sequence motifs predicted to bind key melanocyte transcription factors, located near genes relevant to melanocyte biology, and capable of driving reporter gene expression in melanocytes in culture (86%; 43/50) and in transgenic zebrafish (70%; 7/10). Next, using the sequences of these putative enhancers as a training set for a supervised machine learning algorithm, we develop a vocabulary of 6-mers predictive of melanocyte enhancer function. Lastly, we demonstrate that this vocabulary has genome-wide predictive power in both the mouse and human genomes. This study provides deep insight into the regulation of gene expression in melanocytes and demonstrates a powerful approach to the investigation of regulatory sequences that can be applied to other cell types. [ABSTRACT FROM AUTHOR]

Illuminating the primary sequence encryption of enhancers is central to understanding the regulatory architecture of genomes. We have developed a machine learning approach to decipher motif patterns of hindbrain enhancers and identify 40,000 sequences in the human genome that we predict display regulatory control that includes the hindbrain. Consistent with their roles in hindbrain patterning, MEIS1, NKX6-1, as well as HOX and POU family binding motifs contributed strongly to this enhancer model. Predicted hindbrain enhancers are overrepresented at genes expressed in hindbrain and associated with nervous system development, and primarily reside in the areas of open chromatin. In addition, 77 (0.2%) of these predictions are identified as hindbrain enhancers on the VISTA Enhancer Browser, and 26,000 (60%) overlap enhancer marks (H3K4me1 or H3K27ac). To validate these putative hindbrain enhancers, we selected 55 elements distributed throughout our predictions and six low scoring controls for evaluation in a zebrafish transgenic assay. When assayed in mosaic transgenic embryos, 51/55 elements directed expression in the central nervous system. Furthermore, 30/34 (88%) predicted enhancers analyzed in stable zebrafish transgenic lines directed expression in the larval zebrafish hindbrain. Subsequent analysis of sequence fragments selected based upon motif clustering further confirmed the critical role of the motifs contributing to the classifier. Our results demonstrate the existence of a primary sequence code characteristic to hindbrain enhancers. This code can be accurately extracted using machine-learning approaches and applied successfully for de novo identification of hindbrain enhancers. This study represents a critical step toward the dissection of regulatory control in specific neuronal subtypes. [ABSTRACT FROM AUTHOR]

Abstract: The transcription factor SOX10 has essential roles in neural crest-derived cell populations, including myelinating Schwann cells—specialized glial cells responsible for ensheathing axons in the peripheral nervous system. Importantly, SOX10 directly regulates the expression of genes essential for proper myelin function. To date, only a handful of SOX10 target loci have been characterized in Schwann cells. Addressing this lack of knowledge will provide a better understanding of Schwann cell biology and candidate loci for relevant diseases such as demyelinating peripheral neuropathies. We have identified a highly-conserved SOX10 binding site within an alternative promoter at the SH3-domain kinase binding protein 1 (Sh3kbp1) locus. The genomic segment identified at Sh3kbp1 binds to SOX10 and displays strong promoter activity in Schwann cells in vitro and in vivo. Mutation of the SOX10 binding site ablates promoter activity, and ectopic expression of SOX10 in SOX10-negative cells promotes the expression of endogenous Sh3kbp1. Combined, these data reveal Sh3kbp1 as a novel target of SOX10 and raise important questions regarding the function of SH3KBP1 isoforms in Schwann cells. [Copyright &y& Elsevier]

Abstract: The zebrafish is an ideal model for elucidating the cellular and molecular mechanisms that underlie development of the peripheral nervous system. A transgenic line that selectively labels all the sensory circuits would be a valuable tool for such investigations. In this study, we describe such a line: the enhancer trap zebrafish line Tg(SKIV2L2:gfp) j1775 which expresses green fluorescent protein (gfp) in the peripheral sensory ganglia. We show that this transgene marks all peripheral ganglia and sensory nerves, beginning at the time when the neurons are first extending their processes, but does not label the efferent nerves. The trapped reporter is inserted just upstream of a previously poorly described gene: lhfpl4 on LG6. The expression pattern of this gene by in situ hybridization reveals a different, but overlapping, pattern of expression compared to that of the transgene. This pattern also does not mimic that of the gene (skiv2l2), which provided the promoter element in the construct. These findings indicate that reporter expression is not dictated by an endogenous enhancer element, but instead arises through an unknown mechanism. Regardless, this reporter line should prove to be a valuable tool in the investigation of peripheral nervous system formation in the zebrafish. [Copyright &y& Elsevier]

Background: The ERBB3 gene is essential for the proper development of the neural crest (NC) and its derivative populations such as Schwann cells. As with all cell fate decisions, transcriptional regulatory control plays a significant role in the progressive restriction and specification of NC derived lineages during development. However, little is known about the sequences mediating transcriptional regulation of ERBB3 or the factors that bind them. Results: In this study we identified three transcriptional enhancers at the ERBB3 locus and evaluated their regulatory potential in vitro in NC-derived cell types and in vivo in transgenic zebrafish. One enhancer, termed ERBB3•MCS6, which lies within the first intron of ERBB3, directs the highest reporter expression in vitro and also demonstrates epigenetic marks consistent with enhancer activity. We identify a consensus SOX10 binding site within ERBB3•MCS6 and demonstrate, in vitro, its necessity and sufficiency for the activity of this enhancer. Additionally, we demonstrate that transcription from the endogenous Erbb3 locus is dependent on Sox10. Further we demonstrate in vitro that Sox10 physically interacts with that ERBB3•MCS6. Consistent with its in vitro activity, we also show that ERBB3•MCS6 drives reporter expression in NC cells and a subset of its derivative lineages in vivo in zebrafish in a manner consistent with erbb3b expression. We also demonstrate, using morpholino analysis, that Sox10 is necessary for ERBB3•MCS6 expression in vivo in zebrafish. Conclusions: Taken collectively, our data suggest that ERBB3 may be directly regulated by SOX10, and that this control may in part be facilitated by ERBB3•MCS6. [ABSTRACT FROM AUTHOR]

Background: Myelin protein zero (MPZ) is a critical structural component of myelin in the peripheral nervous system. The MPZ gene is regulated, in part, by the transcription factors SOX10 and EGR2. Mutations in MPZ, SOX10, and EGR2 have been implicated in demyelinating peripheral neuropathies, suggesting that components of this transcriptional network are candidates for harboring disease-causing mutations (or otherwise functional variants) that affect MPZ expression. Methodology: We utilized a combination of multi-species sequence comparisons, transcription factor-binding site predictions, targeted human DNA re-sequencing, and in vitro and in vivo enhancer assays to study human non-coding MPZ variants. Principal Findings: Our efforts revealed a variant within the first intron of MPZ that resides within a previously described SOX10 binding site is associated with decreased enhancer activity, and alters binding of nuclear proteins. Additionally, the genomic segment harboring this variant directs tissue-relevant reporter gene expression in zebrafish. Conclusions: This is the first reported MPZ variant within a cis-acting transcriptional regulatory element. While we were unable to implicate this variant in disease onset, our data suggests that similar non-coding sequences should be screened for mutations in patients with neurological disease. Furthermore, our multi-faceted approach for examining the functional significance of non-coding variants can be readily generalized to study other loci important for myelin structure and function. [ABSTRACT FROM AUTHOR]

Expression profile analysis clusters Gpnmb with known pigment genes, Tyrp1, Dct, and Si. During development, Gpnmb is expressed in a pattern similar to Mitf, Dct and Si with expression vastly reduced in Mitf mutant animals. Unlike Dct and Si, Gpnmb remains expressed in a discrete population of caudal melanoblasts in Sox10-deficient embryos. To understand the transcriptional regulation of Gpnmb we performed a whole genome annotation of 2,460,048 consensus MITF binding sites, and cross-referenced this with evolutionarily conserved genomic sequences at the GPNMB locus. One conserved element, GPNMB-MCS3, contained two MITF consensus sites, significantly increased luciferase activity in melanocytes and was sufficient to drive expression in melanoblasts in vivo. Deletion of the 5′-most MITF consensus site dramatically reduced enhancer activity indicating a significant role for this site in Gpnmb transcriptional regulation. Future analysis of the Gpnmb locus will provide insight into the transcriptional regulation of melanocytes, and Gpnmb expression can be used as a marker for analyzing melanocyte development and disease progression. [ABSTRACT FROM AUTHOR]

Sox10 is a dynamically regulated transcription factor gene that is essential for the development of neural crest-derived and oligodendroglial populations. Developmental genes often require multiple regulatory sequences that integrate discrete and overlapping functions to coordinate their expression. To identify Sox10 cis-regulatory elements, we integrated multiple model systems, including cell-based screens and transposon-mediated transgensis in zebrafish, to scrutinize mammalian conserved, noncoding genomic segments at the mouse Sox10 locus. We demonstrate that eight of 11 Sox10 genomic elements direct reporter gene expression in transgenic zebrafish similar to patterns observed in transgenic mice, despite an absence of observable sequence conservation between mice and zebrafish. Multiple segments direct expression in overlapping populations of neural crest derivatives and glial cells, ranging from pan-Sox10 and pan-neural crest regulatory control to the modulation of expression in subpopulations of Sox10-expressing cells, including developing melanocytes and Schwann cells. Several sequences demonstrate overlapping spatial control, yet direct expression in incompletely overlapping developmental intervals. We were able to partially explain neural crest expression patterns by the presence of head to head SoxE family binding sites within two of the elements. Moreover, we were able to use this transcription factor binding site signature to identify the corresponding zebrafish enhancers in the absence of overall sequence homology. We demonstrate the utility of zebrafish transgenesis as a high-fidelity surrogate in the dissection of mammalian gene regulation, especially those with dynamically controlled developmental expression. [ABSTRACT FROM AUTHOR]

Remarkably, although cardiac disease accounts for the largest proportion of adult mortality and morbidity in the industrialized world, the genetic programs controlling early cardiogenesis are largely incompletely understood. To better understand this process, we set out to identify genes whose expression is enriched within early cardiac fated populations, obtaining the transcriptional signatures of mouse embryonic stem cells (mESCs) at defined intervals during their differentiation along a cardiac path. We compared the RNA profiles of cardiac precursors cells (CPCs) with time-matched non-CPCs and undifferentiated mESCs, using a transgenic mESC line harboring an Nkx2-5 cardiac-specific regulatory sequence driving green fluorescent protein (GFP) to facilitate selection of CPCs. We identify 176 transcripts that are significantly elevated in their abundance within CPCs compared with other assayed populations, predicting that they will likely play a role in cardiogenesis. Of note, approximately 24% (43/176) of the cardiogenic candidate transcripts have known roles in cardiac function or development. Importantly, we evaluated the biological relevance of a significant subset 31/133 (23%) of the remaining candidate genes by in situ hybridization at multiple time points during development (embryonic day, E7.5-9.5) and report that all were expressed in key cardiac structures during cardiogenesis. Furthermore 9/31, of which many were previously uncharacterized, were detected as early as the formation of the cardiac crescent. These data demonstrate the potential power of integrating genomic approaches with mESC differentiation to illuminate developmental processes, and provides a valuable resource that may be mined to further elucidate the genetic programs underlying cardiogenesis. [ABSTRACT FROM AUTHOR]

Although the differentiation of ES cells to cardiomyocytes has been firmly established, the extent to which corresponding cardiac precursor cells can contribute to other cardiac populations remains unclear. To determine the molecular and cellular characteristics of cardiac-fated populations derived from mouse ES (mES) cells, we isolated cardiac progenitor cells (CPCs) from differentiating mES cell cultures by using a reporter cell line that expresses GFP under the control of a cardiac-specific enhancer element of Nkx2-5, a transcription factor expressed early in cardiac development. This ES cell-derived CPC population initially expressed genetic markers of both stem cells and mesoderm, while differentiated CPCs displayed markers of 3 distinct cell lineages (cardiomyocytes, vascular smooth muscle cells, and endothelial cells)--Flk1 (also known as Kdr), c-Kit, and Nkx2-5, but not Brachyury--and subsequently expressed Isl1. Clonally derived CPCs also demonstrated this multipotent phenotype. By transcription profiling of CPCs, we found that mES cell-derived CPCs displayed a transcriptional signature that paralleled in vivo cardiac development. Additionally, these studies suggested the involvement of genes that we believe were previously unknown to play a role in cardiac development. Taken together, our data demonstrate that ES cell-derived CPCs comprise a multipotent precursor population capable of populating multiple cardiac lineages and suggest that ES cell differentiation is a valid model for studying development of multiple cardiac-fated tissues. [ABSTRACT FROM AUTHOR]

RecQ DNA helicases resolve Rad-51-mediated recombination and suppress aberrant homologous recombination. RecQ gene loss is associated with cancer susceptibility and increased mitotic recombination. We have developed an in vivo assay based on a zebrafish pigment mutant for suppression of RecQ activity, and demonstrate that zebrafish RecQ genes have conserved function in suppressing mitotic recombination. [ABSTRACT FROM AUTHOR]

Evolutionary sequence conservation is an accepted criterion to identify noncoding regulatory sequences. We have used a transposon-based transgenic assay in zebrafish to evaluate noncoding sequences at the zebrafish ret locus, conserved among teleosts, and at the human RET locus, conserved among mammals. Most teleost sequences directed ref-specific reporter gene expression, with many displaying overlapping regulatory control. The majority of human RET noncoding sequences also directed ref-specific expression in zebrafish. Thus, vast amounts of functional sequence information may exist that would not be detected by sequence similarity approaches. [ABSTRACT FROM AUTHOR]

The ability to discriminate between deleterious and neutral amino acid substitutions in the genes of patients remains a significant challenge in human genetics. The increasing availability of genomic sequence data from multiple vertebrate species allows inclusion of sequence conservation and physicochemical properties of residues to be used for functional prediction. In this study, the RET receptor tyrosine kinase serves as a model disease gene in which a broad spectrum (≥116) of disease-associated mutations has been identified among patients with Hirschsprung disease and multiple endocrine neoplasia type 2. We report the alignment of the human RET protein sequence with the orthologous sequences of 12 non-human vertebrates (eight mammalian, one avian, and three teleost species), their comparative analysis, the evolutionary to- pology of the RET protein, and predicted tolerance for all published missense mutations. We show that, although evolutionary conservation alone provides significant information to predict the effect of a RET mutation, a model that combines comparative sequence data with analysis of physiochemical properties in a quantitative framework provides far greater accuracy. Although the ability to discern the impact of a mutation is imperfect, our analyses permit substantial discrimination between predicted functional classes of RET mutations and disease severity even for a multigenic disease such as Hirschsprung disease. [ABSTRACT FROM AUTHOR]

Rats harboring the agu mutation have altered behavior and brain pathology resembling human Parkinsonian syndromes; notably, they have a movement disorder and age-progressive dysfunction and death of neurons in the midbrain (substantia nigra pars compacta) that use dopamine as a neurotransmitter. We present evidence that this phenotype is due to a mutation in the rat protein kinase Cγ (in rat, Prkcg; in mouse, Prkcc; in human, PRKCG) gene, which generates a premature stop codon, drastically reducing the level of synthesis of the catalytic domain of the brain-specific protein kinase Cγ protein. [ABSTRACT FROM AUTHOR]

Despite its recognized utility, the extent to which evolutionary sequence conservation-based approaches may systematically overlook functional noncoding sequences remains unclear. We have tiled across sequence encompassing the zebrafish phox2b gene, ultimately evaluating 48 amplicons corresponding to all noncoding sequences therein for enhancer activity in zebrafish. Post hoc analyses of this interval utilizing five commonly used measures of evolutionary constraint (AVID, MLAGAN, SLAGAN, phastCons, WebMCS) demonstrate that each systematically overlooks regulatory sequences. These established algorithms detected only 29%-61% of our identified regulatory elements, consistent with the suggestion that many regulatory sequences may not be readily detected by metrics of sequence constraint. However, we were able to discriminate functional from nonfunctional sequences based upon GC composition and identified position weight matrices (PWM), demonstrating that, in at least one case, deleting sequences containing a subset of these PWMs from one identified regulatory element abrogated its regulatory function. Collectively, these data demonstrate that the noncoding functional component of vertebrate genomes may far exceed estimates predicated on evolutionary constraint. [ABSTRACT FROM AUTHOR]