An evolutionarily conserved nested gene pair - Mab21 and Lrba/Nbea in metazoan

Eukaryotic CD-NTase, STING, and viperin proteins evolved via domain shuffling, horizontal transfer, and ancient inheritance from prokaryotes

Animals use a variety of cell-autonomous innate immune proteins to detect viral infections and prevent replication. Recent studies have discovered that a subset of mammalian antiviral proteins have homology to antiphage defense proteins in bacteria, implying that there are aspects of innate immunity that are shared across the Tree of Life. While the majority of these studies have focused on characterizing the diversity and biochemical functions of the bacterial proteins, the evolutionary relationships between animal and bacterial proteins are less clear. This ambiguity is partly due to the long evolutionary distances separating animal and bacterial proteins, which obscures their relationships. Here, we tackle this problem for 3 innate immune families (CD-NTases [including cGAS], STINGs, and viperins) by deeply sampling protein diversity across eukaryotes. We find that viperins and OAS family CD-NTases are ancient immune proteins, likely inherited since the earliest eukaryotes first arose. In contrast, we find other immune proteins that were acquired via at least 4 independent events of horizontal gene transfer (HGT) from bacteria. Two of these events allowed algae to acquire new bacterial viperins, while 2 more HGT events gave rise to distinct superfamilies of eukaryotic CD-NTases: the cGLR superfamily (containing cGAS) that has since diversified via a series of animal-specific duplications and a previously undefined eSMODS superfamily, which more closely resembles bacterial CD-NTases. Finally, we found that cGAS and STING proteins have substantially different histories, with STING protein domains undergoing convergent domain shuffling in bacteria and eukaryotes. Overall, our findings paint a picture of eukaryotic innate immunity as highly dynamic, where eukaryotes build upon their ancient antiviral repertoires through the reuse of protein domains and by repeatedly sampling a rich reservoir of bacterial antiphage genes.

Eukaryotic antiviral immune proteins arose via convergence, horizontal transfer, and ancient inheritance

Animals use a variety of cell-autonomous innate immune proteins to detect viral infections and prevent replication. Recent studies have discovered that a subset of mammalian antiviral proteins have homology to anti-phage defense proteins in bacteria, implying that there are aspects of innate immunity that are shared across the Tree of Life. While the majority of these studies have focused on characterizing the diversity and biochemical functions of the bacterial proteins, the evolutionary relationships between animal and bacterial proteins are less clear. This ambiguity is partly due to the long evolutionary distances separating animal and bacterial proteins, which obscures their relationships. Here, we tackle this problem for three innate immune families (CD-NTases [including cGAS], STINGs, and Viperins) by deeply sampling protein diversity across eukaryotes. We find that Viperins and OAS family CD-NTases are truly ancient immune proteins, likely inherited since the last eukaryotic common ancestor and possibly longer. In contrast, we find other immune proteins that arose via at least four independent events of horizontal gene transfer (HGT) from bacteria. Two of these events allowed algae to acquire new bacterial viperins, while two more HGT events gave rise to distinct superfamilies of eukaryotic CD-NTases: the Mab21 superfamily (containing cGAS) which has diversified via a series of animal-specific duplications, and a previously undefined eSMODS superfamily, which more closely resembles bacterial CD-NTases. Finally, we found that cGAS and STING proteins have substantially different histories, with STINGs arising via convergent domain shuffling in bacteria and eukaryotes. Overall, our findings paint a picture of eukaryotic innate immunity as highly dynamic, where eukaryotes build upon their ancient antiviral repertoires through the reuse of protein domains and by repeatedly sampling a rich reservoir of bacterial anti-phage genes.

DNA-methylome-assisted classification of patients with poor prognostic subventricular zone associated IDH-wildtype glioblastoma

Glioblastoma (GBM) derived from the “stem cell” rich subventricular zone (SVZ) may constitute a therapy-refractory subgroup of tumors associated with poor prognosis. Risk stratification for these cases is necessary but is curtailed by error prone imaging-based evaluation. Therefore, we aimed to establish a robust DNA methylome-based classification of SVZ GBM and subsequently decipher underlying molecular characteristics. MRI assessment of SVZ association was performed in a retrospective training set of IDH-wildtype GBM patients (n = 54) uniformly treated with postoperative chemoradiotherapy. DNA isolated from FFPE samples was subject to methylome and copy number variation (CNV) analysis using Illumina Platform and cnAnalysis450k package. Deep next-generation sequencing (NGS) of a panel of 130 GBM-related genes was conducted (Agilent SureSelect/Illumina). Methylome, transcriptome, CNV, MRI, and mutational profiles of SVZ GBM were further evaluated in a confirmatory cohort of 132 patients (TCGA/TCIA). A 15 CpG SVZ methylation signature (SVZM) was discovered based on clustering and random forest analysis. One third of CpG in the SVZM were associated with MAB21L2/LRBA. There was a 14.8% (n = 8) discordance between SVZM vs. MRI classification. Re-analysis of these patients favored SVZM classification with a hazard ratio (HR) for OS of 2.48 [95% CI 1.35–4.58], p = 0.004 vs. 1.83 [1.0–3.35], p = 0.049 for MRI classification. In the validation cohort, consensus MRI based assignment was achieved in 62% of patients with an intraclass correlation (ICC) of 0.51 and non-significant HR for OS (2.03 [0.81–5.09], p = 0.133). In contrast, SVZM identified two prognostically distinct subgroups (HR 3.08 [1.24–7.66], p = 0.016). CNV alterations revealed loss of chromosome 10 in SVZM– and gains on chromosome 19 in SVZM– tumors. SVZM– tumors were also enriched for differentially mutated genes (p < 0.001). In summary, SVZM classification provides a novel means for stratifying GBM patients with poor prognosis and deciphering molecular mechanisms governing aggressive tumor phenotypes.

Functional analysis of a de novo variant in the neurodevelopment and generalized epilepsy disease gene NBEA

Neurobeachin (NBEA) was initially identified as a candidate gene for autism. Recently, variants in NBEA have been associated with neurodevelopmental delay and childhood epilepsy. Here, we report on a novel NBEA missense variant (c.5899G > A, p.Gly1967Arg) in the Domain of Unknown Function 1088 (DUF1088) identified in a child enrolled in the Undiagnosed Diseases Network (UDN), who presented with neurodevelopmental delay and seizures. Modeling of this variant in the Caenorhabditis elegans NBEA ortholog, sel-2, indicated that the variant was damaging to in vivo function as evidenced by altered cell fate determination and trafficking of potassium channels in neurons. The variant effect was indistinguishable from that of the reference null mutation suggesting that the variant is a strong hypomorph or a complete loss-of-function. Our experimental data provide strong support for the molecular diagnosis and pathogenicity of the NBEA p.Gly1967Arg variant and the importance of the DUF1088 for NBEA function.

Zebrafish mab21l2 mutants possess severe defects in optic cup morphogenesis, lens and cornea development

BACKGROUND

Mutations in MAB21L2 result in severe ocular defects including microphthalmia, anophthalmia, coloboma, microcornea, and cataracts. The molecular and cellular underpinnings of these defects are unknown, as is the normal cellular function of MAB21L2. Zebrafish mab21l2 ^au10 mutants possess ocular defects resembling those in humans with MAB21L2 mutations, providing an excellent model to characterize mab21l2 functions during eye development.

RESULTS

mab21l2 ^-/- mutants possessed a host of ocular defects including microphthalmia and colobomas as well as small, disorganized lenses and cornea dysgenesis. Decreased proliferation, increased cell death, and defects in marker gene expression were detected in the lens. Cell death in the optic stalk was elevated in mab21l2 ^-/- mutants and the basement membrane between the edges of the choroid fissure failed to break down. Neuronal differentiation in the retina was normal, however. mab21l2 ^-/- mutant corneas were disorganized, possessed an increased number of cells, some of which proliferated ectopically, and failed to differentiate the corneal stroma.

CONCLUSIONS

mab21l2 function is required for morphogenesis and cell survival in the lens and optic cup, and basement membrane breakdown in the choroid fissure. mab21l2 function also regulates proliferation in the lens and cornea; in its absence, the lens is small and mispatterned, and corneal morphogenesis and patterning are also disrupted.

Generation and characterization of pathogenic Mab21l2(R51C) mouse model

MAB21L2(R51C) is one of the five documented MAB21L2 mutations in human patients with bilateral eye malformations identified via whole exome sequencing. In addition to the eye abnormality, patients with MAB21L2 ^R51C/+ mutation also have skeletal dysplasia and intellectual disability. To evaluate the pathology of this mutant allele systematically in understanding the functional role of MAB21L2 in human development, we introduce the R51C mutation into the mouse genome by CRISPR/Cas9 system to generate a mouse model for detailed characterization. The Mab21l2 ^R51C/+ mice have eyeless phenotype and skeletal abnormalities. Micro-computed tomography (micro-CT) analysis showed the Mab21l2 ^R51C/+ mice have no eye balls but with two abnormal tissues underneath the brain. Histological analysis revealed that the early eye development in the mutant embryos is interrupted. In addition, Mab21l2 ^R51C/+ mice also have joint fusion phenotype; the humerus is fused with radius, whereas femur is fused with tibia. Limbs in the mutant animals are distinctly shorter than the wild type; and deltoid tuberosities in humeri are absent in these Mab21l2 ^R51C/+ mice. In summary, we showed that our Mab21l2 ^R51C/+ mutant mice have recapitulated the pathological features in eye and bone of human patients. Further analyses of the mutant phenotype with molecular markers will provide insight on how MAB21L2 guides the optic differentiation and skeletogenesis, revealing specific underlying pathogenic mechanism of the MAB21L2(R51C) mutation.

Identification of Variants in RET and IHH Pathway Members in a Large Family With History of Hirschsprung Disease

BACKGROUND & AIMS

Hirschsprung disease (HSCR) is an inherited congenital disorder characterized by absence of enteric ganglia in the distal part of the gut. Variants in ret proto-oncogene (RET) have been associated with up to 50% of familial and 35% of sporadic cases. We searched for variants that affect disease risk in a large, multigenerational family with history of HSCR in a linkage region previously associated with the disease (4q31.3-q32.3) and exome wide.

METHODS

We performed exome sequencing analyses of a family in the Netherlands with 5 members diagnosed with HSCR and 2 members diagnosed with functional constipation. We initially focused on variants in genes located in 4q31.3-q32.3; however, we also performed an exome-wide analysis in which known HSCR or HSCR-associated gene variants predicted to be deleterious were prioritized for further analysis. Candidate genes were expressed in HEK293, COS-7, and Neuro-2a cells and analyzed by luciferase and immunoblot assays. Morpholinos were designed to target exons of candidate genes and injected into 1-cell stage zebrafish embryos. Embryos were allowed to develop and stained for enteric neurons.

RESULTS

Within the linkage region, we identified 1 putative splice variant in the lipopolysaccharide responsive beige-like anchor protein gene (LRBA). Functional assays could not confirm its predicted effect on messenger RNA splicing or on expression of the mab-21 like 2 gene (MAB21L2), which is embedded in LRBA. Zebrafish that developed following injection of the lrba morpholino had a shortened body axis and subtle gut morphological defects, but no significant reduction in number of enteric neurons compared with controls. Outside the linkage region, members of 1 branch of the family carried a previously unidentified RET variant or an in-frame deletion in the glial cell line derived neurotrophic factor gene (GDNF), which encodes a ligand of RET. This deletion was located 6 base pairs before the last codon. We also found variants in the Indian hedgehog gene (IHH) and its mediator, the transcription factor GLI family zinc finger 3 (GLI3). When expressed in cells, the RET-P399L variant disrupted protein glycosylation and had altered phosphorylation following activation by GDNF. The deletion in GDNF prevented secretion of its gene product, reducing RET activation, and the IHH-Q51K variant reduced expression of the transcription factor GLI1. Injection of morpholinos that target ihh reduced the number of enteric neurons to 13% ± 1.4% of control zebrafish.

CONCLUSIONS

In a study of a large family with history of HSCR, we identified variants in LRBA, RET, the gene encoding the RET ligand (GDNF), IHH, and a gene encoding a mediator of IHH signaling (GLI3). These variants altered functions of the gene products when expressed in cells and knockout of ihh reduced the number of enteric neurons in the zebrafish gut.

The BEACH Protein LRBA Promotes the Localization of the Heterotrimeric G-protein Golf to Olfactory Cilia

BEACH domain proteins are involved in membrane protein traffic and human diseases, but their molecular mechanisms are not understood. The BEACH protein LRBA has been implicated in immune response and cell proliferation, and human LRBA mutations cause severe immune deficiency. Here, we report a first functional and molecular phenotype outside the immune system of LRBA-knockout mice: compromised olfaction, manifesting in reduced electro-olfactogram response amplitude, impaired food-finding efficiency, and smaller olfactory bulbs. LRBA is prominently expressed in olfactory and vomeronasal chemosensory neurons of wild-type mice. Olfactory impairment in the LRBA-KO is explained by markedly reduced concentrations (20–40% of wild-type levels) of all three subunits α_olf, β₁ and γ₁₃ of the olfactory heterotrimeric G-protein, G_olf, in the sensory cilia of olfactory neurons. In contrast, cilia morphology and the concentrations of many other proteins of olfactory cilia are not or only slightly affected. LRBA is also highly expressed in photoreceptor cells, another cell type with a specialized sensory cilium and heterotrimeric G-protein-based signalling; however, visual function appeared unimpaired by the LRBA-KO. To our knowledge, this is the first observation that a BEACH protein is required for the efficient subcellular localization of a lipid-anchored protein, and of a ciliary protein.

Genome wide association study identifies variants in NBEA associated with migraine in bipolar disorder

BACKGROUND

Migraine is a common comorbidity among individuals with bipolar disorder, but the underlying mechanisms for this co-occurrence are poorly understood. The aim of this study was to investigate the genetic background of bipolar patients with and without migraine.

METHODS

We performed a genome-wide association analysis contrasting 460 bipolar migraneurs with 914 bipolar patients without migraine from the Bipolar Genome Study (BiGS).

RESULTS

We identified one genome-wide significant association between migraine in bipolar disorder patients and rs1160720, an intronic single nucleotide polymorphism (SNP) in the NBEA gene (P=2.97 × 10(-8), OR: 1.82, 95% CI: 1.47-2.25), although this was not replicated in a smaller sample of 289 migraine cases.

LIMITATIONS

Our study is based on self-reported migraine.

CONCLUSIONS

NBEA encodes neurobeachin, a scaffolding protein primarily expressed in the brain and involved in trafficking of vesicles containing neurotransmitter receptors. This locus has not previously been implicated in migraine per se. We found no evidence of association in data from the GWAS migraine meta-analysis consortium (n=118,710 participants) suggesting that the association might be specific to migraine co-morbid with bipolar disorder.

2mit, an intronic gene of Drosophila melanogaster timeless2, is involved in behavioral plasticity

Background

Intronic genes represent ~6% of the total gene complement in Drosophila melanogaster and ~85% of them encode for proteins. We recently characterized the D. melanogastertimeless2 (tim2) gene, showing its active involvement in chromosomal stability and light synchronization of the adult circadian clock. The protein coding gene named 2mit maps on the 11^thtim2 intron in the opposite transcriptional orientation.

Methodology/Principal Findings

Here we report the molecular and functional characterization of 2mit. The 2mit gene is expressed throughout Drosophila development, localizing mainly in the nervous system during embryogenesis and mostly in the mushroom bodies and ellipsoid body of the central complex in the adult brain. Insilico analyses revealed that 2mit encodes a putative leucine-Rich Repeat transmembrane receptor with intrinsically disordered regions, harboring several fully conserved functional interaction motifs in the cytosolic side. Using insertional mutations, tissue-specific over-expression, and down-regulation approaches, it was found that 2mit is implicated in adult short-term memory, assessed by a courtship conditioning assay. In D. melanogaster, tim2 and 2mit do not seem to be functionally related. Bioinformatic analyses identified 2MIT orthologs in 21 Drosophilidae, 4 Lepidoptera and in Apis mellifera. In addition, the tim2-2mit host-nested gene organization was shown to be present in A. mellifera and maintained among Drosophila species. Within the Drosophilidae 2mit-hosting tim2 intron, insilico approaches detected a neuronal specific transcriptional binding site which might have contributed to preserve the specific host-nested gene association across Drosophila species.

Conclusions/Significance

Taken together, these results indicate that 2mit, a gene mainly expressed in the nervous system, has a role in the behavioral plasticity of the adult Drosophila. The presence of a putative 2mit regulatory enhancer within the 2mit-hosting tim2 intron could be considered an evolutionary constraint potentially involved in maintaining the tim2-2mit host-nested chromosomal architecture during the evolution of Drosophila species.

The BEACH is hot: a LYST of emerging roles for BEACH-domain containing proteins in human disease

BEACH (named after 'Beige and Chediak-Higashi') is a conserved ∼280 residue domain, present in nine human BEACH domain containing proteins (BDCPs). Most BDCPs are large, containing a PH-like domain for membrane association preceding their BEACH domain, and containing WD40 and other domains for ligand binding. Recent studies found that mutations in individual BDCPs cause several human diseases. BDCP alterations affect lysosome size (LYST and NSMAF), apoptosis (NSMAF), autophagy (LYST, WDFY3, LRBA), granule size (LYST, NBEAL2, NBEA) or synapse formation (NBEA). However, the roles of each BDCP in these membrane events remain controversial. After reviewing studies on individual BDCPs, we propose a unifying hypothesis that BDCPs act as scaffolding proteins that facilitate membrane events, including both fission and fusion, determined by their binding partners. BDCPs may also bind each other, enabling fusion or fission of vesicles that are not necessarily of the same type. Such mechanisms explain why different BDCPs may have roles in autophagy; each BDCP is specific for the cell type or the cargo, but not necessarily specific for attaching to the autophagosome. Further elucidation of these mechanisms, preferably carrying out the same experiment on multiple BDCPs, and possibly using patients' cells, may identify potential targets for therapy.

Drosophila rugose is a functional homolog of mammalian Neurobeachin and affects synaptic architecture, brain morphology, and associative learning

Neurobeachin (Nbea) is implicated in vesicle trafficking in the regulatory secretory pathway, but details on its molecular function are currently unknown. We have used Drosophila melanogaster mutants for rugose (rg), the Drosophila homolog of Nbea, to further elucidate the function of this multidomain protein. Rg is expressed in a granular pattern reminiscent of the Golgi network in neuronal cell bodies and colocalizes with transgenic Nbea, suggesting a function in secretory regulation. In contrast to Nbea(-/-) mice, rg null mutants are viable and fertile and exhibit aberrant associative odor learning, changes in gross brain morphology, and synaptic architecture as determined at the larval neuromuscular junction. At the same time, basal synaptic transmission is essentially unaffected, suggesting that structural and functional aspects are separable. Rg phenotypes can be rescued by a Drosophila rg+ transgene, whereas a mouse Nbea transgene rescues aversive odor learning and synaptic architecture; it fails to rescue brain morphology and appetitive odor learning. This dissociation between the functional redundancy of either the mouse or the fly transgene suggests that their complex composition of numerous functional and highly conserved domains support independent functions. We propose that the detailed compendium of phenotypes exhibited by the Drosophila rg null mutant provided here will serve as a test bed for dissecting the different functional domains of BEACH (for beige and human Chediak-Higashi syndrome) proteins, such as Rugose, mouse Nbea, or Nbea orthologs in other species, such as human.

Deleterious mutations in LRBA are associated with a syndrome of immune deficiency and autoimmunity

Most autosomal genetic causes of childhood-onset hypogammaglobulinemia are currently not well understood. Most affected individuals are simplex cases, but both autosomal-dominant and autosomal-recessive inheritance have been described. We performed genetic linkage analysis in consanguineous families affected by hypogammaglobulinemia. Four consanguineous families with childhood-onset humoral immune deficiency and features of autoimmunity shared genotype evidence for a linkage interval on chromosome 4q. Sequencing of positional candidate genes revealed that in each family, affected individuals had a distinct homozygous mutation in LRBA (lipopolysaccharide responsive beige-like anchor protein). All LRBA mutations segregated with the disease because homozygous individuals showed hypogammaglobulinemia and autoimmunity, whereas heterozygous individuals were healthy. These mutations were absent in healthy controls. Individuals with homozygous LRBA mutations had no LRBA, had disturbed B cell development, defective in vitro B cell activation, plasmablast formation, and immunoglobulin secretion, and had low proliferative responses. We conclude that mutations in LRBA cause an immune deficiency characterized by defects in B cell activation and autophagy and by susceptibility to apoptosis, all of which are associated with a clinical phenotype of hypogammaglobulinemia and autoimmunity.

Subgroup-specific alternative splicing in medulloblastoma

Medulloblastoma comprises four distinct molecular variants: WNT, SHH, Group 3, and Group 4. We analyzed alternative splicing usage in 14 normal cerebellar samples and 103 medulloblastomas of known subgroup. Medulloblastoma samples have a statistically significant increase in alternative splicing as compared to normal fetal cerebella (2.3-times; P < 6.47E-8). Splicing patterns are distinct and specific between molecular subgroups. Unsupervised hierarchical clustering of alternative splicing events accurately assigns medulloblastomas to their correct subgroup. Subgroup-specific splicing and alternative promoter usage was most prevalent in Group 3 (19.4%) and SHH (16.2%) medulloblastomas, while observed less frequently in WNT (3.2%), and Group 4 (9.3%) tumors. Functional annotation of alternatively spliced genes reveals overrepresentation of genes important for neuronal development. Alternative splicing events in medulloblastoma may be regulated in part by the correlative expression of antisense transcripts, suggesting a possible mechanism affecting subgroup-specific alternative splicing. Our results identify additional candidate markers for medulloblastoma subgroup affiliation, further support the existence of distinct subgroups of the disease, and demonstrate an additional level of transcriptional heterogeneity between medulloblastoma subgroups.

Family size and turnover rates among several classes of small non-protein-coding RNA genes in Caenorhabditis nematodes

It is important to understand the forces that shape the size and evolutionary histories of gene families. Here, we investigated the evolution of non–protein-coding RNA genes in the genomes of Caenorhabditis nematodes. We specifically focused on nested arrangements, that is, cases in which an RNA gene is entirely contained in an intron of another gene. Comparing these arrangements between species simplifies the inference of orthology and, therefore, of evolutionary fates of nested genes. Two distinct patterns are evident in the data. Genes encoding small nuclear RNAs (snRNAs) and transfer RNAs form large families, which have persisted since before the common ancestor of Metazoa. Yet, individual genes die relatively rapidly, with few orthologs having survived since the divergence of Caenorhabditis elegans and Caenorhabditis briggsae. In contrast, genes encoding small nucleolar RNAs (snoRNAs) are either single-copy or form small families. Individual snoRNAs turn over at a relatively slow rate—most C. elegans genes have clearly identifiable orthologs in C. briggsae. We also found that in Drosophila, genes from larger snRNA families die at a faster rate than their counterparts from single-gene families. These results suggest that a relationship between family size and the rate of gene turnover may be a general feature of genome evolution.

Dendritic spine formation and synaptic function require neurobeachin

A challenge in neuroscience is to understand the mechanisms underlying synapse formation. Most excitatory synapses in the brain are built on spines, which are actin-rich protrusions from dendrites. Spines are a major substrate of brain plasticity, and spine pathologies are observed in various mental illnesses. Here we investigate the role of neurobeachin (Nbea), a multidomain protein previously linked to cases of autism, in synaptogenesis. We show that deletion of Nbea leads to reduced numbers of spinous synapses in cultured neurons from complete knockouts and in cortical tissue from heterozygous mice, accompanied by altered miniature postsynaptic currents. In addition, excitatory synapses terminate mostly at dendritic shafts instead of spine heads in Nbea mutants, and actin becomes less enriched synaptically. As actin and synaptopodin, a spine-associated protein with actin-bundling activity, accumulate ectopically near the Golgi apparatus of mutant neurons, a role emerges for Nbea in trafficking important cargo to pre- and postsynaptic compartments.

Most excitatory synapses in the brain are found on dendritic spines, but the mechanisms underlying synapse formation are poorly understood. Niesmann et al. investigate the role of neurobeachin in synaptogenesis, and find that its deletion leads to fewer spinous synapses and altered postsynaptic currents.

mab21l2 transgenics reveal novel expression patterns of mab21l1 and mab21l2, and conserved promoter regulation without sequence conservation

mab21l1 and mab21l2 paralogs have widespread and dynamic expression patterns during vertebrate development. Both genes are expressed in the developing eye, midbrain, neural tube, and branchial arches. Our goal was to identify promoter regions with activity in mab21l2 expression domains. Assays of mab21l2 promoter-EGFP constructs in zebrafish embryos confirm that constructs containing 7.2 or 4.9 kb of mab21l2 promoter region are sufficient to drive expression in known (e.g., tectum, branchial arches) and unexpected domains (e.g., lens and retinal amacrine cells). A comparative analysis identifies complementary and novel expression domains of endogenous mab21l2 (e.g., lens and ventral iridocorneal canal) and mab21l1 (e.g., retinal amacrine and ganglion cells). Interestingly, therefore, despite the absence of conserved non-coding elements, a 4.9-kb mab21l2 promoter is sufficient to recapitulate expression in tissues unique to mab21l1 or mab21l2.

When needles look like hay: how to find tissue-specific enhancers in model organism genomes

A major prerequisite for the investigation of tissue-specific processes is the identification of cis-regulatory elements. No generally applicable technique is available to distinguish them from any other type of genomic non-coding sequence. Therefore, researchers often have to identify these elements by elaborate in vivo screens, testing individual regions until the right one is found. Here, based on many examples from the literature, we summarize how functional enhancers have been isolated from other elements in the genome and how they have been characterized in transgenic animals. Covering computational and experimental studies, we provide an overview of the global properties of cis-regulatory elements, like their specific interactions with promoters and target gene distances. We describe conserved non-coding elements (CNEs) and their internal structure, nucleotide composition, binding site clustering and overlap, with a special focus on developmental enhancers. Conflicting data and unresolved questions on the nature of these elements are highlighted. Our comprehensive overview of the experimental shortcuts that have been found in the different model organism communities and the new field of high-throughput assays should help during the preparation phase of a screen for enhancers. The review is accompanied by a list of general guidelines for such a project.