<title>Background</title> <p>The jumonji C (JmjC)-domain containing proteins are an extending family of redox enzymes that catalyse a wide range of oxidation reactions. These proteins are found in all living organisms and are characterised by sharing the highly conserved JmjC domain. This domain was first described by Takeuchi and colleagues, who isolated a gene in the mouse which they called <italic>jumonji </italic>(meaning cruciform in Japanese), in regard to an abnormal phenotype visible in mutant embryos during neural plate development [##REF##7758946##1##]. Since then, the JmjC domain has found its entry into public domain databases and many other proteins containing this domain have been annotated or identified subsequently [##REF##10838566##2##, ####REF##11165500##3##, ##REF##16983801##4##, ##REF##16715513##5####16715513##5##]. Currently, more then 11,000 sequence entries of JmjC domain containing proteins are present in Uniprot, PFAM, Interpro and SMART, thus demonstrating the tremendous expansion of members in this protein superfamily. Structural studies showed that the JmjC domain has a very characteristic topology. It forms a double-stranded β-helix (DSBH) fold, also known as jelly-roll motif, that normally consists of eight antiparallel β-strands [##REF##16513174##6##]. This DSBH topology is characteristic for the cupin metalloenzymes [##REF##11738598##7##] and for the JmjC domain containing proteins [##REF##11165500##3##]. Cupin metalloenzymes constitute a large and diverse superfamily of proteins with enzymatically active and non-active members. They contain at least one DSBH or jelly-roll motive. Secondary structure elements that surround the DSBH are used further to define different subfamilies (for review see [##REF##16513174##6##]). The DSBH fold in JmjC domain containing proteins can form an enzymatically active pocket by coordinating di-iron (Fe(II)) and the co-substrate 2-oxoglutarate (2OG). In almost all cases, the two-electron oxidation of the "prime" substrate is coupled to the conversion of 2OG into succinate and CO<sub>2 </sub>[##REF##10607676##8##]. Therefore, this class of enzymes is also known as nonheme-Fe(II)-2-oxoglutarate-dependent dioxygenases.</p> <p>Oxidative reactions catalysed by 2OG-dependent dioxygenase are critical steps in the biosynthesis of metabolites, or post translational modification of interacting target proteins involving hydroxylation or demethylation and DNA/RNA repair of N-methylated nucleic acids [##REF##16513174##6##,##REF##15809658##9##]. As such, 2OG-dependent dioxygenases are of immanent biological importance. A number of JmjC domain containing proteins are involved in regulating key developmental processes during mammalian embryogenesis by acting as transcription factors. An example of this subgroup of JmjC domain containing proteins is jumonji itself (now called Jarid2), which can act as a transcriptional repressor. Jarid2 inhibits cyclin D1 transcription during heart and brain development and thereby controls cellular proliferation and differentiation in these organs during morphogenesis [##REF##10807864##10##, ####REF##10446263##11##, ##REF##10913339##12####10913339##12##]. Other examples of JmjC-domain proteins that function as transcriptional repressors are the retinoblastoma-binding protein 2 (Jarid1a, also known as RBP2), which is required for the differentiation of various cell types including muscle, bone and myeloid cells [##REF##15949438##13##] and hairless, a nuclear co-repressor of multiple nuclear hormone receptors such as the vitamin D receptor, the thyroid hormone receptor and the retinoic acid-receptor-related orphan receptor α [##REF##12847098##14##, ####REF##12215540##15##, ##REF##11641275##16####11641275##16##]. Mutations in hairless result in humans and mice in alopecia, a form of hereditary hair loss due to defects in the regeneration of hair follicle cells [##REF##9445480##17##, ####REF##9736769##18##, ##REF##15280217##19####15280217##19##]. More recently, JmjC-domain containing proteins have gained emerging importance in histone modification and chromatin regulation. The JmjC domain containing histone demethylases (JHDMs) can catalyse lysine demethylation of histones through an oxidative reaction that requires iron Fe(II) and 2OG. Histone methylation marks have important roles in regulating gene expression and are central to the control of epigenetic information that regulates cell fate and identity [##REF##16261189##20##]. Until now, JHDMs have been reported to be capable to reverse H3K36 (JHDM1/FBXL11, [##REF##16362057##21##], JMJD2A, [##REF##16677698##22##,##REF##16603238##23##], Rph1, [##REF##17371840##24##,##REF##17369256##25##], Jhd2, [##REF##17369256##25##], Gis1, [##REF##17369256##25##], Jhd1, [##REF##17369256##25##]), H3K9 (JHDM2A/JMJD1A, [##REF##17363312##26##]), and both H3K9 and H3K36 methylation (JHDM3A to JHDM3D/JMJD2A to JMJD2D, [##REF##16603238##23##,##REF##16732293##27##, ####REF##16738407##28##, ##REF##16732292##29##, ##REF##17277772##30####17277772##30##]). In addition, the JARID subfamily of JmjC-domain containing proteins have been shown to possess H3K4 demethylase activity in different species (JARID1A/RBP2, [##REF##17320161##31##,##REF##17320163##32##], JARID1C/SMCX, [##REF##17320161##31##,##REF##17320160##33##], JARID1D/SMCY, [##REF##17320162##34##], JARID1B/PLU-1/Yjr119Cp/Jhdp2, [##REF##17320161##31##,##REF##17310254##35##,##REF##17310255##36##], Lid, [##REF##17363312##26##,##REF##17351630##37##, ####REF##17351631##38##, ##REF##17311883##39####17311883##39##], and JMJ2 [##REF##17550896##40##]). Most recently, the JmjC-domain containing proteins UTX and JMJD3 have been shown to function as H3K27 specific demethylases and demonstrated to play an important role in animal body patterning and regulation of inflammatory responses [##REF##17713478##41##, ####REF##17825402##42##, ##REF##17851529##43##, ##REF##17761849##44####17761849##44##]. Through this activity the JHDMs have been postulated to regulate diverse cellular processes including control of cell proliferation, heterochromatin assembly and gene expression [##REF##17363312##26##,##REF##16738407##28##,##REF##17277772##30##,##REF##17320162##34##,##REF##16603237##45##,##REF##16762840##46##]. However, for most of the newly discovered JHDMs gene ablation studies <italic>in vivo </italic>have as yet not been carried out leaving the possibility open that in addition to histones other substrates might be targets of these proteins.</p> <p>Many of the different classes of JmjC domain containing proteins possess complex domain architectures by having at least one additional protein domain besides the JmjC domain. Here, the JHDMs and Jarid proteins are good examples. In addition to the JmjC domain they have PHD, Tudor, Bright/Arid and zinc finger domains, which have recently been shown to interact with methylated histone tails and are required for their catalytic function [##REF##17320160##33##,##REF##17351631##38##,##REF##16601153##47##]. In contrast, other members in this protein superfamily have no other recognisable protein domain apart from the JmjC domain. Important members of this subgroup of JmjC domain only proteins are the 2OG-dependent dioxygenase factor inhibiting HIF (Hif1an, also known as FIH) and the PHD hydroxylases, which are involved in the maintenance of oxygen homeostasis. Hif1an was the first structurally resolved mammalian JmjC domain containing protein [##REF##12432100##48##,##REF##12446723##49##] and was shown to function as an asparagine hydroxylase for the hypoxia inducible factor α (HIF1α), a transcription factor that is central to oxygen homeostasis in both physiological and pathophysiological processes (for review see [##REF##15122348##50##]). The PHD hydroxylases PHD1 to PHD3 also target HIF1α by hydroxylating proline residues. Under normal oxygen conditions the hydroxylation reactions carried out by PHDs induce the proteasomal degradation of the HIF1α subunit, whereas under hypoxia HIF1α is not hydroxylated and translocates from the cytoplasm to the nucleus where it activates transcription of genes needed for adaptation to low oxygen levels. AlkB is another member of the JmjC domain only proteins. As a 2OG-dependent oxygenase AlkB can repair DNA that is methylated at 1-methyladenine and 3-methylcytosine by oxidation of the methyl group [##REF##12226668##51##,##REF##12226667##52##]. AlkB and related proteins define a subfamily of small 2OG-dependent oxygenases that are involved in DNA and RNA repair [##REF##11276424##53##].</p> <p>In the past, we have analysed the function of the <italic>phosphatidylserine receptor </italic>(<italic>Ptdsr </italic>or also known as <italic>PSR </italic>or <italic>PtdSerR</italic>), a gene that also encodes a JmjC domain only protein [##REF##15345036##54##,##REF##15615595##55##]. First described as a cell surface receptor for the recognition and engulfment of apoptotic cells [##REF##10811223##56##], recent studies have raised serious doubts about its involvement in apoptotic cell clearance. Five groups independently reported a nuclear localisation of Ptdsr, which clearly contradicts the hypothesis that the protein might function as a transmembrane receptor [##REF##15193161##57##, ####REF##14729065##58##, ##REF##17522160##59##, ##REF##16317002##60##, ##REF##17534701##61####17534701##61##]. Furthermore, a comprehensive analysis carried out by us in <italic>Ptdsr</italic>-knockout mice demonstrated no impairment in the removal of apoptotic cells in <italic>Ptdsr</italic>-deficient embryos and macrophages ([##REF##15345036##54##] for <italic>Ptdsr</italic><sup><italic>tm</italic>1<italic>Gbf</italic></sup>) thereby contradicting two previous studies that reported defects in apoptotic cell clearance in two other <italic>Ptdsr</italic>-knockout mice ([##REF##14645847##62##] for <italic>Ptdsr</italic><sup><italic>tm</italic>1<italic>Flv</italic></sup>, and [##REF##14715629##63##] for <italic>Ptdsr</italic><sup><italic>tm</italic>1<italic>Ysfk</italic></sup>). All studies involving <italic>Ptdsr</italic>-knockout mice conclusively demonstrated that the gene is essential for normal development. <italic>Ptdsr</italic>-deficient mice are dying at birth and display severe differentiation defects in many tissues such as brain, eyes, lung, kidney, liver and intestine at different stages of embryogenesis [##REF##15345036##54##,##REF##14645847##62##,##REF##14715629##63##]. Furthermore, <italic>Ptdsr </italic>plays an important role during heart development and ablation of <italic>Ptdsr </italic>function is associated with complex cardio-pulmonary malformations that resemble the human congenital heart syndrome Tetralogy of Fallot [##REF##15615595##55##].</p> <p>Although <italic>Ptdsr </italic>is clearly an important differentiation-promoting gene during embryogenesis its molecular function is still unknown. Based on a bioinformatic analysis of the cloned <italic>Hydra </italic>homologue <italic>PSR</italic>, Cikala and colleagues were the first to suggest that the protein might function as a 2OG-dependent dioxygenase that is capable of modifying unknown nuclear targets [##REF##15193161##57##]. In recognition of the still unknown function of the "phosphatidylserine receptor gene" the International Committee for Standardized Genetic Nomenclature in Mice (ICSGNM) has most recently revised its name. Ptdsr has been renamed into <italic>jumonji domain containing gene 6 </italic>(<italic>Jmjd6</italic>). We will obey the guidelines of the ICSGNM and use this name and the gene symbol <italic>Jmjd6 </italic>throughout the remainder of this article.</p> <p>To gain more insights into possible Jmjd6 functions using sequence homology approaches, we developed a semiautomated process to extract all available Jmjd6 homologous protein and nucleic acid sequences from public databases. This comprehensive data set was used to identify proteins with sequence homology to Jmjd6 in 62 species across all major living phyla. In addition, this allowed us to identify new <italic>Jmjd6/JMJD6 </italic>splice variants in mice and humans and to analyse the genomic neighbourhood of the gene in vertebrates. Immediately downstream of <italic>Jmjd6 </italic>we found the gene <italic>1110005A03Rik</italic>, which is predicted to encode a putative methyltransferase, which is as yet not further biochemically characterised. Interestingly, we found <italic>1110005A03Rik </italic>to be affected by deletions in two previously reported <italic>Jmjd6 </italic>knockout mice. To predict possible structure-function relationships of Jmjd6, we employed a comparative modelling approach of the three-dimensional structure of the protein and annotated probable functional residues based on sequence conservation. To integrate Jmjd6 into the superfamily of JmjC-domain containing proteins, we assembled a non-redundant set of 313 DSBH-domain containing proteins by using our semiautomatic retrieval script and grouped these according to phylogenetic relationships. Taken together, our results further support a conserved enzymatic function of Jmjd6 most likely as a classical non-heme, Fe(II)- and 2OG-dependent dioxygenase. It also shows that Jmjd6 is closely related to Jmjd4 and JFP6, and constitutes together with these two proteins a subgroup of JmjC domain containing proteins. Moreover, we identified in protista, fungi and plants proteins with sequence homology to Jmjd6 that diverge from the classical "JmjC-only domain" structure. Some of these proteins have in comparison to the animalia Jmjd6 extended C-terminal regions that contain additional domains that might have been lost during evolution.</p> <p>With these approaches we establish here an overall framework of information that will facilitate not only further research regarding the functional characterisation of Jmjd6 but also of other members within the DSBH-domain containing superfamily of proteins.</p>
<title>Methods</title> <title>Bioinformatic approaches</title> <p>To investigate the splice variants of <italic>Jmjd6</italic>, the genomic sequence of <italic>Jmjd6 </italic>and adjacent genes on either side was extracted to perform an unrestricted e2g-based annotation of mouse ESTs with a seed-length of 18 and a minimum identity of 95% [##REF##15215398##67##]. The resulting list of ESTs that overlap with the known <italic>Jmjd6 </italic>coding-region was passed on to a self-developed script that automatically retrieves the sequences of all matches and performs an in-depth verification of each EST using Sim4 [##REF##9750195##68##]. The script sorts matches in different groups, according to the match quality. All matches with unidirectional splice sites, without internal gaps in the mRNA, and a &lt;95% identity score to the corresponding genomic regions leaving aside the poly-A tail were used for manual annotation.</p> <p>PipMaker sequence analysis [##REF##10779500##66##] was done by extracting genomic sequences of 50 to 100 kb encompassing the <italic>Jmjd6 </italic>locus and neighbouring loci of mouse, human, bovine, armadillo, opossum, zebrafish, chicken, western clawed frog and puffer fish. Sequences were analysed by RepeatMasker [##UREF##1##97##] and positions of all coding sequences/exons in <italic>Jmjd6</italic>, <italic>1110005A03Rik</italic>, and <italic>Sfrs2 </italic>were identified and annotated using the BLAST algorithm. The masked and annotated mouse sequence were then finally used for interspecies pairwise alignment with the sequences from the different species mentioned above using the program PipMaker [##REF##10779500##66##].</p> <p>To identify <italic>Jmjd6 </italic>orthologous sequences, a standard NCBI blastp-search against the nr database using amino acids 52–152 of the murine <italic>Jmjd6 </italic>protein was performed [##UREF##2##98##]. Identified sequences were validated using the "reciprocal best BLAST hit" criterion in a <italic>Mus musculus</italic>-restricted blastp-search against the nr database [##REF##9381173##72##].</p> <p>SMART-parser, a Perl-script for batch queries against SMART [##UREF##3##99##], data retrieval and processing was developed to identify proteins with similar domains as Jmjd6. SMART and SMART-parser were used to automatically identify, extract and annotate the sequences of all JmjC-, Pfam-JmjC, Blast-JmjC-, P4Hc-, Pfam-2OG-FeII-Oxy- and Pfam-PhyH-domains for later processing. These domain names have been used in the different public databases to indicate the presence of the JmjC domain containing DSBH fold proteins.</p> <p>Multiple sequence alignments were generated employing ClustalW with a default setup and five iterations of re-alignment. Visualisation, calculation of a background frequency corrected sequence logo, and calculation of amino acid conservation was done using LaTeX and TEXshade [##REF##10842735##100##]. A ClustalW alignment of all Jmjd6 orthologues (amino acids 44–329) identified in the prior approaches was used to generate a 3D-structure of Jmjd6 using the pipeline from the MPI Toolkit [##REF##16845021##101##]. To add 2OG and Fe(II) to the predicted structure and the predicted active site of Jmjd6, an iterative backbone annealing of the Hif1an structure (1h2k) on top of the Jmjd6 model was performed using the SwissPDB Viewer [##UREF##4##102##]. Afterwards, the coordinates of 2OG and Fe(II) were added to the Jmjd6 model and energies were minimised using DeepView. 3D representations of molecules were rendered using PyMOL [##UREF##5##103##].</p> <p>For a maximum likelihood phylogenetic analysis of Jmjd6 orthologues, we employed the PHYLIP toolbox [##UREF##6##104##]. A ClustalW alignment of all proteins with significant homology to Jmjd6 was used as an input for ProML, using the Jones-Taylor-Thornton algorithm with an estimated constant rate of mutations, global rearrangements and the more accurate iterated branch lengths option. The primary protein sequences of Jmjd6/JMJD6 were analysed for sumoylation sites using a SUMOplot prediction with a cutt-off of 90% probability [##UREF##7##105##]. Nuclear export signals (NES) were identified using the NetNES 1.1 server using a cut-off of 95% identity to the NES Hidden Markow Model in the database [##UREF##8##106##]. Secondary structure predictions were generated using the Jpred web service [##UREF##9##107##].</p> <title>Identification of <italic>Jmjd6 </italic>splice variants</title> <p>RNA was isolated from different cells (HEK 293-T, A549, mouse embryonic fibroblasts (MEFs) and reverse transcribed. Murine splice variants were amplified using the forward primer 5'-GGTGGTGCCTCTTCCCAACAAA-3' and the reverse primers 5'-TCCAGGTCAGGGTTGGGACAC-3' and 5'-AGCTAGAAGAGTCGCTGGAGCTGTC-3'. The amplicons of the primer pairs had a product size of 363 bp (Mm3), 407 bp (Mm2), and 432 bp (Mm1), respectively. All fragments were analysed by sequencing and the sequences were submitted to GenBank (Acc. No. <ext-link ext-link-type="gen" xlink:href="EF527404">EF527404</ext-link> and <ext-link ext-link-type="gen" xlink:href="EF527405">EF527405</ext-link>). Human cDNAs (C1244142, C1244152, C1244171-10, C1244035-10, C1244149, C1244264-10, C1244183, C1244246-10, C1244122-10, C1244260-10) were obtained from the BioChain Institute (Hayward, USA) and amplified according to the manufactures instructions using the forward primer 5'-GGGAACTCATCAAAGTGACCCGAGACG-3' and the reverse primers 5'-ACTGCCGCTGCCCGTGCTGTATC-3' and 5'-CGAGTCTGCGAGGACTGCCAACTCG-3'. The amplicons had a product size of 320 bp (Hs2), 344 bp (Hs1) and 438 bp (Hs3), respectively. All fragments were analysed by sequencing and the sequences were submitted to GenBank (Acc. No. <ext-link ext-link-type="gen" xlink:href="EF527406">EF527406</ext-link> and <ext-link ext-link-type="gen" xlink:href="EF527407">EF527407</ext-link>).</p> <title>Expression analysis of <italic>Jmjd6</italic>, <italic>1110005A03Rik </italic>and <italic>Sfrs2</italic></title> <p>RNA was isolated and reverse transcribed. Expression of <italic>Jmjd6, 1110005A03Rik, Sfrs2 </italic>in embryos or cells were analysed using the primers 5'-GTT CCA GCT CGT CAG ACT CG-3' and 5'-TGC CCC TAA GAC ATG ACC AC-3' for <italic>Jmjd6</italic>, 5'-GCC ACC ACA AGA CAT CAT TCT TG-3' and 5'-TCA GAT AAT TCA GCT TTA TGC CAG G-3' for <italic>1110005A03Rik</italic>, 5'-GCT CCA GAT CAA CCT CCA AG-3' and 5'-GCC ACC TGA GGC AGA TTA AA-3' for <italic>Sfrs2</italic>, respectively.</p> <title>Analysis of the bidirectional transcriptional unit</title> <p>An 841 bp and an 840 bp fragment was amplified from mouse C57BL/6J genomic DNA using the forward 5'-CTCGAGCTCTCGTAGTAGTTGTGCCGGGT-3' and reverse 5'-GCTAGCACTGTCCCTAAATGTGTCACTGGAGC-3' primer as well as the forward 3'-GCTAGCTCTCGTAGTAGTTGTGCCGGGT-5' and reverse 5'-CTCGAGACTGTCCCTAAATGTGTCACTGGAGC-3' primer, respectively. These fragments were integrated into a promoter-less luciferase expression vector (pGL3, Promega, Madison, USA) and these constructs were introduced into HEK 293-T cells, A549 cells or mouse embryonic fibroblasts (MEFs) using Nanofection (PAA, Pasching, Austria). After 24h of cultivation at 37°C and 5% CO<sub>2</sub>, cell extracts were generated to measure the luciferase-activity using the Dual-Luciferase-Reporter-Assay System (Promega) according to the manufactures instructions. The obtained results were normalized to pGL4.74 (Promega) luciferase-activity. The measurements were performed in triplicate and repeated at least three times.</p>
<title>Results</title> <title>Retrieval and analysis of a comprehensive set of Jmjd6 orthologs and other JmjC domain containing proteins</title> <p>We developed a semiautomatic process to extract all available Jmjd6 homologous protein and nucleic acid sequences from public databases. The computational pipeline consists of three parallel procedures based on "DNA", "protein", and "protein domain" searches (Figure ##FIG##0##1##). With this approach, we generated a candidate list of 11.376 sequence entries including all putative <italic>Jmjd6 </italic>orthologs and all JmjC/DSBH domain containing proteins present in Genbank, Ensembl, Uniprot, PFAM, Interpro, and SMART, independent of their annotation. The goal of our semiautomated computational pipeline was then to identify and systematically analyse all orthologous <italic>Jmjd6 </italic>loci, all transcripts encoded by the murine and human <italic>Jmjd6/JMJD6 </italic>loci, the degree and distribution of sequence conservations in Jmjd6 proteins, and their phylogenetic relationship to other JmjC/DSBH domain containing proteins. In the end, the combination of these analysis approaches enabled us to identify 62 proteins with significant sequence homology to Jmjd6 and to establish a set of 313 non-redundant DSBH domain containing proteins in human and eight model organisms.</p> <p>We started the "DNA" sequence based retrieval and analysis approach (depicted in yellow in Figure ##FIG##0##1##) by identification of all mouse <italic>Jmjd6 </italic>homologous sequences available in Ensembl (release 44). This candidate loci list contained homologous <italic>Jmjd6 </italic>loci in fully annotated genomes as well as <italic>Jmjd6 </italic>sequence fragments from ongoing sequencing projects. To identify conserved coding and non-coding sequences within these genomic regions we undertook a PipMaker analysis of retrieved sequences that showed a similar genomic organisation in comparison to the murine <italic>Jmjd6 </italic>locus. Next, we mapped all available murine and human ESTs in dbEST to the murine and human <italic>Jmjd6</italic>/<italic>JMJD6 </italic>loci, respectively, using the program e2g. Extensive validations of each EST employing Sim4 with our own perlscript for data retrieval, handling, examination and sorting of the results led to a comprehensive analysis of the murine and human <italic>Jmjd6</italic>/<italic>JMJD6 </italic>loci and a detailed annotation of all encoded transcripts.</p> <p>To identify all <italic>Jmjd6 </italic>orthologs, we employed a blastp search against Ensembl and nr databases using a unique N-terminal sequence of the murine Jmjd6 protein as a template ("protein" based search and analysis, depicted in Figure ##FIG##0##1## in green). The candidate matches were individually verified in a <italic>Mus musculus </italic>restricted blastp database search using the reciprocal best BLAST hit criterion and a perlscript for automatic data processing. These led to the identification of Jmjd6 proteins in other species that are closely related to and might be orthologs of the murine Jmjd6 protein. A full length ClustalW alignment of proteins with significant homology to Jmjd6 was processed using a TeXshade based conservation analysis and also as input for a 3D structure prediction of the catalytic core domain. By combining these results, we developed a novel way of predicting residues that are most likely essential for the catalytic activity of the protein based on their degree of conservation.</p> <p>To analyse the phylogenetic evolution of the Jmjd6 JmjC domain, we developed SMART-parser. SMART and SMART-parser were used to analyse the domain composition and extract the sequence information from a set of sequences that was derived from the Ensembl, Interpro, Uniprot and SMART databases ("DSBH-domain" routine, depicted in Figure ##FIG##0##1## in blue). This set comprises all JmjC-, Pfam-JmjC, Blast-JmjC-, P4Hc-, Pfam-2OG-FeII-Oxy- and Pfam-PhyH-domain containing sequences from human, mouse, zebrafish, pufferfish, nematode, fruit fly, and yeast. These different domain names have been used in the past to annotate the JmjC domain in DSBH fold containing proteins. The SMART-parser verified sequences were used as an input for a global ClustalW alignment. A PHYLIP processing of all ClustalW alignments allowed the analysis of the Jmjd6 phylogeny in relationship to all other identified DSBH domain containing proteins.</p> <title>Comparative analysis of <italic>Jmjd6 </italic>homologous loci in vertebrates</title> <p>In the mouse, <italic>Jmjd6 </italic>is localised in a gene dense region on distal chromosome 11. Using different gene targeting strategies three loss-of-function alleles of <italic>Jmjd6 </italic>have been generated that display conflicting phenotypes [##REF##15345036##54##,##REF##14645847##62##,##REF##14715629##63##]. To analyse the immediate genomic neighbourhood of <italic>Jmjd6 </italic>and to identify all transcripts encoded by the locus, we undertook a comparative sequence analysis of this gene region in different species. In Ensembl, the mouse <italic>Jmjd6 </italic>gene is annotated to consist of six exons compromising 6 kb (Ensembl release 44). Immediately downstream of <italic>Jmjd6 </italic>two genes are located, <italic>1110005A03Rik</italic>, an expressed sequence tag (EST) encoding a putative arginine protein methyltransferase [##REF##16257219##64##], and <italic>Sfrs2</italic>, encoding the splicing factor SC 35 [##REF##8530103##65##]. To analyse the genomic region encompassing <italic>Jmjd6 </italic>and <italic>1110005A03Rik </italic>in more detail, we performed a comprehensive database screen. First, we searched all fully sequenced genomes present in Ensembl build release 44 for an annotated <italic>Jmjd6 </italic>gene. We identified <italic>Jmjd6 </italic>homologous genes or genomic <italic>Jmjd6 </italic>fragments in 17 mammalian, 1 avian, 1 amphibian, 5 fish, 2 urochordata, and 4 invertebrate genomes (Table ##TAB##0##1##). The identified putative <italic>Jmjd6 </italic>orthologous loci were then searched for the presence of the <italic>1110005A03Rik </italic>gene. In 22 out of the 30 analysed genomes we found both genes located close together. Remarkably, in all analysed invertebrate species <italic>Jmjd6 </italic>and <italic>1110005A03Rik </italic>were not found to co-localise despite the presence of both genes at other chromosomal locations in the examined genomes. Next, we extracted 50 to 100 kb of genomic sequence encompassing the <italic>Jmjd6 </italic>locus and adjacent loci of these 22 species and used these for an initial PipMaker sequence analysis (data not shown). To carry this out, the mouse sequence was first analysed for repetitive sequence elements using RepeatMasker and secondly, the positions of all coding sequences/exons in <italic>Jmjd6</italic>, <italic>1110005A03Rik</italic>, and <italic>Sfrs2 </italic>were identified and annotated using the BLAST algorithm. The masked and annotated mouse sequence was then finally used for interspecies pairwise alignment with the human, bovine, armadillo, opossum, zebrafish, chicken, western clawed frog, and puffer fish sequences using the program PipMaker [##REF##10779500##66##]. These species were selected for the analysis either based on their large phylogenetic distance to the mouse or on the availability of high quality, full-coverage sequence information in the extracted regions. The resulting alignments of similar regions in the genomic sequences were summarised as "percent identity plot" (pip) displaying sequence identity in combination with the annotated repeats and transcriptional units (Figure ##FIG##1##2##). The pip shows a high degree of sequence conservation in the coding regions of the three genes and confirms the accuracy of the initial exon assignments annotated by Ensembl. However, a detailed inspection of the <italic>Jmjd6 </italic>locus identified further high-level conserved sequences (= 90%) in between the annotated exons 4 and 5 (red overlay in Figure ##FIG##1##2##). This sequence conservation suggested the presence of an up to now overseen additional exon, at least in the mammalian species. In addition, the pip confirms the close association of the <italic>Jmjd6 </italic>and <italic>1110005A03Rik </italic>genes in all analysed organisms. <italic>Jmjd6 </italic>and <italic>1110005A03Rik </italic>are arranged in a head-to-head transcriptional orientation with a short intergenic distance of 65 bp between both genes. This suggests that these two genes build a conserved transcriptional unit that possibly shares a bi-directional promoter region. This is further supported by the presence of a conserved short non-coding sequence in between the translational start sites of these two genes (orange overlay in Figure ##FIG##1##2##). Interestingly, we observed not only a conserved association of the <italic>Jmjd6 </italic>and <italic>1110005A03Rik </italic>genes in this genomic region, but also of <italic>Sfrs2</italic>, which is present in this region in all species with the exception of the puffer fish (Figure ##FIG##1##2##). In the mouse, <italic>Sfrs2 </italic>is located 6.3 kb distal of <italic>Jmjd6 </italic>and 170 bp distal to <italic>1110005A03Rik </italic>(Figure ##FIG##1##2##).</p> <p>To further test the hypothesis that an additional exon in <italic>Jmjd6 </italic>might generate alternative transcripts of this gene, we investigated all available mouse <italic>Jmjd6 </italic>EST data currently present in dbEST. Using e2g, a web-based server that allows an efficient mapping of large EST and cDNA datasets to genomic regions [##REF##15215398##67##], all <italic>Jmjd6 </italic>ESTs were mapped to the mouse locus and adjacent regions. The resulting EST hits represent the entire putative <italic>Jmjd6 </italic>coding region including the transcribed and untranslated regions (UTRs). These cDNA sequences were then extracted and individually aligned with the mouse <italic>Jmjd6 </italic>genomic sequence using Sim4, a computer program that allows an efficient assignment of cDNAs to large genomic regions by accounting for the position of introns [##REF##9750195##68##]. The Sim4 validated EST sequences were then manually curated and all predicted <italic>Jmjd6 </italic>transcripts were assembled and annotated to the genomic locus. As a result of this analysis, we identified a new variant of exon 4 in <italic>Jmjd6 </italic>and found an additional, before overseen exon in the conserved region between exons 4 and 5. This allowed us, solely based on the EST data, to predict the presence of two novel <italic>Jmjd6 </italic>splice variants in addition to the main transcript of 1.4 kb currently annotated in Ensembl (Figure ##FIG##2##3A##). To validate <italic>in vivo </italic>the existence of the predicted splice variants, we performed RT-PCR expression analysis. Various major organs from adult mice as well as whole embryos at different developmental stages (E) were used for the analysis. RT-PCR amplification yielded three DNA fragments indicative for the predicted transcripts in all analysed samples at varying levels (Figure ##FIG##2##3B##). Taken together, this confirmed the computational prediction of two new <italic>Jmjd6 </italic>splice variants encoded by the mouse locus and demonstrated that <italic>Jmjd6 </italic>consists of seven exons, with exon 4 occurring in two different variants 4a and 4b, respectively. At the protein level both new detected splice variants (transcripts 2 and 3 in Figure ##FIG##2##3A##) affect the C-terminal part of the Jmjd6 protein. In comparison to the main transcript annotated in Ensembl (transcript 1 in Figure ##FIG##2##3A##), either a stretch of 9 amino acids starting at position 302 is exchanged against 31 other residues (transcript Mm3 in Figure ##FIG##2##3C##) or a premature stop codon is introduced resulting in a truncated protein of 314 amino acids (transcript Mm2 in Figure ##FIG##2##3C##). Using the CBS-prediction service, we identified in this region of the protein a putative sumoylation site (probability score 92%) and a nuclear export signal (probability score 96%) using NetNES. Furthermore, this region contains the previously described putative AT-hook domain [##REF##15193161##57##] and the poly-serine stretch domain ([##REF##11536008##69##,##REF##16877497##70##] Figure ##FIG##2##3C##). All these potential interaction sites are not present in the protein encoded by transcript Mm2. Transcript Mm3 differs in comparison with the main transcript (Mm1) in a stretch of 31 amino replacing the putative AT-hook containing region (Figure ##FIG##2##3C##).</p> <p>Based on the identification of the highly conserved sequence coding for the new exon 5 in all mammalian species (Figure ##FIG##1##2## and data not shown), we performed a similar e2g and Sim4 analysis using the human <italic>JMJD6 </italic>locus and human EST database entries. This analysis confirmed the presence of exon 5 containing <italic>JMJD6 </italic>transcripts in humans and also the presence of two different variants of exon 4 in the human genome. Using again RT-PCR, we validated the expression of these transcripts <italic>in vivo</italic>. Similar to the mouse we found expression of the three transcripts in 10 human tissues [see Additional file ##SUPPL##0##1##]. These, as well as the mouse alternative transcripts, were further validated and confirmed by sequencing (Accession numbers EF527404 for Mm2 and EF527405 for Mm3 for mouse <italic>Jmjd6 </italic>and EF527406 for Hs2 and EF527407 for Hs3 for human <italic>JMJD6 </italic>transcript variants, respectively). In comparison to the new identified mouse transcript variants, both human transcript variants are predicted to generate truncated proteins. Neither of these predicted proteins appear to encode a nuclear export signal, sumoylation site, AT hook domain or poly-serine stretch at the C-terminal part of the protein (Fig. ##FIG##2##3C##). As an insertion of premature translation termination codons might result in nonsense-mediated decay (NMD) of the corresponding mRNA [##REF##16723977##71##], we analysed if the mouse transcript variant 2 and human splicing variants 2 and 3 might be subjected to degradation by the NMD pathway. Using RT-PCR assays for each of these variants we found that this is not the case. Primers located in the new introduced exons and at the 3'-end of each transcript yielded products of the expected size, indicating that these variants are indeed expressed and most likely functional [see Additional file ##SUPPL##1##2##]. The human and mouse splice variants seem to be expressed at different mRNA levels depending on the tissue analyzed (see Figure ##FIG##2##3## and [see Additional file ##SUPPL##1##2##]), with mouse splice variant Mm1 expressed in higher amounts in the developing embryo (Figure ##FIG##2##3##). Interestingly, an e2g and Sim4 analysis performed in non-mammalian species did not result in the recognition of mouse and human related <italic>Jmjd6 </italic>splice variants in these organisms (data not shown).</p> <p>To analyse the conserved, putative bi-directional promoter element in between <italic>Jmjd6 </italic>and <italic>1110005A03Rik</italic>, we tried to identify potential regulatory motifs in the genomic sequence encompassing the intergenic region, the 5'-UTRs and the first exons of both genes. We isolated a 840 bp fragment covering the genomic region between <italic>Jmjd6 </italic>and <italic>1110005A03Rik </italic>(Figure ##FIG##3##4A##) that contained several predicted transcription factor binding sites based on a MatInspector analysis (data not shown). This fragment was further analysed using a promoter-less luciferase reporter gene assay. In comparison to controls, the selected genomic fragment clearly possesses transcriptional promoter activity when transfected into primary mouse embryonic fibroblast cells (MEFs) and various cell lines (Figure ##FIG##3##4B## and data not shown). Moreover, this transcription activating activity was found to be orientation dependent. In the forward strand orientation (which controls <italic>in vivo 1110005A03Rik </italic>expression according to the genomic annotation in Ensembl) the fragment showed only half of the activity as compared to the reverse strand orientation (Figure ##FIG##3##4B##), which controls expression of <italic>Jmjd6 in vivo</italic>. This demonstrated that the short intergenic region between <italic>Jmjd6 </italic>and <italic>1110005A03Rik </italic>is capable to control basal expression of both genes, although this bi-directional transcriptional unit seems to activate both genes at different levels.</p> <p>We then went on to investigate the consequences on <italic>Jmjd6</italic>, <italic>1110005A03Rik </italic>and <italic>Sfrs2 </italic>expression of the different gene targeting strategies used to generate the three <italic>Jmjd6 </italic>mouse knockout alleles. In all three <italic>Jmjd6 </italic>mutant mice (<italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Flv</italic></sup>, generated by the Flavell laboratory; <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Ysfk</italic></sup>, generated by the Fukui laboratory; and <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Gbf</italic></sup>, generated in our laboratory) the first two or three exons of the gene have been deleted. However, an alignment of restriction enzyme sites used to construct the different gene targeting vectors reveals that in <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Flv </italic></sup>and <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Ysfk </italic></sup>mice, in contrast to the <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Gbf </italic></sup>mouse line, exon 1 of <italic>1110005A03Rik </italic>was deleted in addition [see Figure A in Additional file ##SUPPL##2##3##]. This was confirmed when we analysed the expression of all three genes in homozygous embryos harvested from the three mutant mouse lines. <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Flv</italic>-/- </sup>and <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Ysfk</italic>-/- </sup>animals show no expression of <italic>1110005A03Rik </italic>and <italic>Jmjd6</italic>, while the expression of <italic>Sfrs2 </italic>remains unaffected by the gene targeting. In <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Gbf </italic></sup>homozygous embryos <italic>1110005A03Rik</italic>, and <italic>Sfrs2 </italic>are expressed comparable to wild type levels [see Figure B in Additional file ##SUPPL##2##3##] and as previously shown [##REF##15345036##54##], only the expression of <italic>Jmjd6 </italic>is ablated in this mouse mutant line. This was further supported by semiquantitative expression analysis of <italic>1110005A03Rik </italic>in <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Gbf</italic>-/- </sup>embryos and primary fetal liver-derived macrophage cultures. RT-PCR analysis of serial diluted cDNA shows equal expression of <italic>1110005A03Rik </italic>in <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Gbf </italic></sup>animals in comparison to wildtype controls in all analysed samples [see Figure C in Additional file ##SUPPL##2##3##]. Thus, the <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Flv </italic></sup>and <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Ysfk</italic></sup>alleles represent double knockouts of the <italic>Jmjd6 </italic>and <italic>1110005A03Rik </italic>genes, while the <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Gbf </italic></sup>allele ablates only <italic>Jmjd6 </italic>gene expression.</p> <title>Phylogenetic analysis of the Jmjd6 protein</title> <p>Following the identification of Jmjd6 orthologues in 30 of the 34 species annotated in Ensembl (Table ##TAB##0##1##), which implied a broad distribution of <italic>Jmjd6 </italic>during evolution, we continued to search for additional homologous proteins in other species. To identify further candidate proteins, we performed a comprehensive database (nr) search using the blastp algorithm and the N-terminal part encompassing amino acids 53 to 152 of the mouse Jmjd6 protein. For the search strategy, we explicitly excluded the highly conserved JmjC domain to reduce the number of false positive hits representing Jmjd6-unrelated JmjC domain containing proteins. To ascertain that we select for the subsequent analysis in each species the sequence with the highest homology to Jmjd6, we used the "reciprocal best BLAST hit" criterion [##REF##9381173##72##]. The obtained full-length database entries were used for a <italic>Mus musculus</italic>-restricted blastp-search of the nr database and for each species only those sequences matching Jmjd6 as the best BLAST hit criterion were retained for the analysis. As some sequences are only available from early stage annotation projects (low-sequence coverage genomes, [##REF##17148474##73##]), we selected 54 out of 62 sequences due to quality measures as a working package for the calculation of a multiple sequence alignment by employing ClustalW. Following five iterations of re-alignment of the protein sequences the multiple sequence alignment was shaded with respect to amino acid identity/similarity and a protein sequence logo was generated and adjusted to the background frequency of each amino acid in the alignment. In addition, predicted secondary structure elements and SMART domains of the mouse Jmjd6 protein were added (Figure ##FIG##4##5##). The retrieved reciprocal best BLAST hits from the different species and the quality of the multiple sequence alignment allowed us to confirm the high homology of these 54 proteins to the murine Jmjd6 protein (Figure ##FIG##4##5##). Many of those are with no doubt true orthologous proteins of the mouse Jmjd6. These are found in the mammalia, aves, amphibia, actinopterygii, urchordata, echinodermata, insecta, nematoda, platyhelminthes, and cnideria and demonstrate the remarkable conservation of the <italic>Jmjd6 </italic>gene during evolution (Figure ##FIG##4##5## and [see Additional file ##SUPPL##3##4##]) as previously also noted by others ([##REF##10811223##56##,##REF##15193161##57##]). However, some proteins identified in plants, protista, fungi differ from the "JmjC only" domain structure characteristic for the Jmjd6 proteins of animalia. For example, some of them contain nearly five times more amino acids or have in addition to the JmjC-domain single F-box domains [see Additional files ##SUPPL##3##4## and ##SUPPL##4##9##]. Nevertheless, they are the closest relatives to the mouse Jmjd6 as defined by the best reciprocal BLAST hit criterion. Thus, it is difficult to decide if these proteins might be true orthologs of the animalia Jmjd6 proteins. If they are, it could mean that additional protein domains have been lost from an ancient Jmjd6 protein during evolution in the animalia lineage.</p> <p>To explore potentially important residues for Jmjd6 function we compared the degree of sequence conservation across the aligned sequences and included in the alignment predicted tertiary structure elements of the protein. First, our analysis showed the presence of a 2-His-1-carboxylate facial triad (HxDx<sub>n</sub>H), as a common structural motif among non-heme iron(II) enzymes, present in all identified Jmjd6 proteins (Figure ##FIG##4##5##), suggesting an evolutionary conserved enzymatic function of Jmjd6. We could identify within the predicted JmjC domain eight β-strands (I-VIII) that have been shown to form the characteristic jelly-roll topology in other structurally resolved JmjC domain containing proteins [##REF##16513174##6##]. The amino acids encompassing these β-strands are highly conserved in all proteins. 50% of all β-strands forming amino acids are identical in at least 85% of the analysed species suggesting the spatial and functional conservation of the JmjC-domain in all analysed proteins. The only exception we found is in between β-strands IV and V. At this position we found an insertion of sequences of variable amino acid conservation among species and of varying length as it has been previously described for other JmjC domain containing proteins [##REF##16513174##6##,##REF##16677698##22##,##REF##11276424##53##]. Second, within the C-terminal region of the Jmjd6 protein preceding the last β-strands of the recognised jelly-roll the degree of sequence conservation is decreased. In contrast to the metazoan species showing a high degree of sequence identity, all other species display a significantly reduced similarity in this region. Furthermore, only the metazoan species have a poly-serine stretch in the C-terminal region of the Jmjd6 protein (Figure ##FIG##4##5## and [##REF##16877497##70##]). We could not identify any comparable sequence motif in the other species but looking for chemical properties of amino acids present at this position in proteins of non-metazoan species we found an enrichment of polar amino acids in this part of the protein (data not shown). This suggests at least a conservation and fortification of physico-chemical properties in this region during evolution. Third, the part N-terminal of the JmjC domain shows in comparison to the C-terminus a higher degree of conservation. We found that 20% of the amino acids 59–140 are identical in at least 90% of the analysed proteins (Figure ##FIG##4##5##). In addition, we identified one highly conserved and three less conserved β-strands in this region. These findings suggest the presence of an important structural element of the Jmjd6 protein in close vicinity to the predicted JmjC domain (see below).</p> <title>Comparative modelling and analysis of the conserved putative catalytic core domain of Jmjd6</title> <p>Due to the fact that up to now no crystal structure of Jmjd6 is available, we decided to perform a homology/comparative modelling of the three-dimensional structure of the protein. The multiple sequence alignment of Jmjd6 homologous proteins (amino acids 44–329; Figure ##FIG##4##5##) was used for pairwise comparison of profile hidden Markov models (HMMs). This generated an alignment of the protein sequence that can be modelled with known related structures. This first step showed a significant probability (= 95%) of structural similarity of Jmjd6 to the three structurally resolved 2OG- and Fe(II)-dependent dioxygenases Hif1an (PDB entry <ext-link ext-link-type="pdb" xlink:href="1h2k">1h2k</ext-link>), putative asparaginyl hydroxylase – 2636534 from <italic>Bacillus subtilis </italic>(PDB entry <ext-link ext-link-type="pdb" xlink:href="1vrb">1vrb</ext-link>) and Jmjd2a (PDB entry <ext-link ext-link-type="pdb" xlink:href="2gp5">2gp5</ext-link>), respectively. The alignment generated by HHpred [##REF##15531603##74##] was finally prepared to create a structural model and MODELLER [##REF##8254673##75##] calculated this comparative model of the Jmjd6 protein based on the identified similarities containing all non-hydrogen atoms (Figure ##FIG##5##6##). The comparative modelling of Jmjd6 predicted that the Jmjd6 fold is a mixed alpha-beta structure composed of a major and a minor β-sheet surrounded by several α-helices. The structure showed the presence of the 2OG- and Fe(II)-dependent dioxygenase typical double-stranded β-helix (DSBH) motif as the protein core. It consists of two anti-parallel β-sheets composed of eight β-strands. The major sheet is defined by β-strands I, VIII, III, the minor sheet by β-strands VI and II, VII, IV, and V (Figure ##FIG##5##6## and [see Additional file ##SUPPL##5##7##]), which form a jelly-roll topology. In addition, we recognised in the model three β-strands extending the core DSBH, resulting in a seven-stranded major β-sheet (green arrows in Figure ##FIG##4##5## and orange arrows in Figure ##FIG##5##6A##). Interestingly, these additional strands were found to be paired anti-parallel to each other and to the major β-sheet.</p> <p>Following iterative backbone annealing based on the structure of Hif1an (PDB entry <ext-link ext-link-type="pdb" xlink:href="1h2k">1h2k</ext-link>), 2OG and Fe(II) were added to the 3D model of Jmjd6. After energy minimisation utilising the SwissPDB Viewer, the predicted active site of the Jmjd6-Fe(II)-2OG complex was analysed in detail. Here, we found the proper positioning of the 2-His-1-carboxylate facial triad chelating the iron atom in a manner similar to other 2OG- and Fe(II)-dependent dioxygenases (Figure ##FIG##5##6B##). The HxD motif (His 187 and Asp 189) is predicted to be located just after the sequence of DSBH β-strand II and the second Fe(II)-binding histidine (His 273) is located on DSBH β-strand VII. The observation of Fe(II) coordination by these three residues confirms their initially suggested assignment based on the primary sequence analysis of the multiple sequence alignment of Jmjd6 proteins (Figure ##FIG##4##5##). In addition, we could identify in the model several residues within the putative active core of Jmjd6 capable of positioning the co-substrate 2OG in the cavity formed by the interior of the DSBH. The proper positioning of 2OG relative to the Fe(II) seems to be achieved by an interaction of several residues of the active core that are close enough to the hydrogen-bond with the 2OG. Among these residues we predict Asn197 and Trp174 to interact with the 1-carboxylate oxygen of 2OG that is not ligated to the Fe(II). Asn197 and Trp174 are stabilised by Thr285. Interestingly, in Jmjd6, like in Hif1an and Jmjd2a, a lysine (Lys204, from DSBH β-strand IV) could be important for the binding of the C5'-carboxylate of 2OG. Most of the other structurally resolved 2OG- and Fe(II)-dependent dioxygenases have instead of the threonine 285 an arginine at this position, which is part of a RXS motif within the DSBH β-strand VIII [##REF##8223617##76##,##REF##10777523##77##]. Thus, Jmjd6 has in contrast to other 2OG-dependent dioxygenases no tripeptide RXS motif (Figure ##FIG##4##5##). In addition to these major 2OG interacting residues in Jmjd6, we could recognise additional residues (Phe133 and Val275) for putative hydrophobic interactions in the active core of the protein.</p> <p>Beside the high probability of the Jmjd6 structure deduced from the comparative modelling approach, we searched for an additional, independent validation of this prediction. To achieve this, we decided to combine phylogenetic sequence analysis data and our comparative modelling data to proof the reliability of the Jmjd6 protein structure prediction. Accordingly, the amino acid identity across all 54 analysed species was translated into an indicative colour code ranging from dark red, representing amino acid identity across all species, to dark blue for highly variable/non-conserved residues (Figure ##FIG##4##5##). This colour code was applied to the calculated Jmjd6 structure and the individual amino acids were coloured according to their conservation (Figure ##FIG##6##7##). This computational approach demonstrated a strong conservation of the predicted Jmjd6 catalytic triad during evolution with 100% sequence identity of the HxD motif in all analysed 54 species (Figure ##FIG##6##7A##). In addition, residues that are predicted to bind 2-OG (Asn197, Thr285, Trp174, Lys204, and Thr187) are also highly conserved (Figure ##FIG##6##7B##), strongly suggesting that Jmjd6 is indeed possessing a 2OG-dependent dioxygenase function. Furthermore, a surface model of the protein with annotated sequence conservation shows that the putative active centre is located in a deep and narrow pocket [see Additional file ##SUPPL##6##5##]. Interestingly, the surface protein model also shows conservation of an exposed N-terminal region between β-strands 2 and 3 of the major β-sheet that is extending the DSBH core indicating that these residues might be important for interaction with other proteins [see Additional file ##SUPPL##6##5##]. In contrast to the sequence conservation observed in the predicted catalytic core and the exposed N-terminal major β-sheet, other regions of the protein have been much less conserved during evolution. Thus, the major selective pressure seems to have been on the catalytic function of the protein, while other areas of the protein retained flexibility in their functional design.</p> <title>Integration of Jmjd6 into the JmjC domain containing protein family</title> <p>By using the computational pipeline described in Figure ##FIG##0##1##, we extracted a non-redundant set of 313 JmjC/DSBH domain-containing proteins from the Ensembl, Interpro, Uniprot and SMART databases. We selected these proteins from human and eight model organisms (mouse, zebrafish, pufferfish, tetraodon, nematode, fruitfly and yeast) based on the sequence homology of their JmjC domain to the mouse Jmjd6 and on the quality of the available sequence. Because genome-sequencing projects of many species are still ongoing, we had to exclude JmjC domain encoding genes from some species with low-sequence coverage genomes that did not match our quality criteria. Nevertheless, our comprehensive catalogue included 61 JmjC/DSBH-domain containing proteins from <italic>Homo sapiens</italic>, 57 from <italic>Mus musculus</italic>, 49 from <italic>Danio rerio</italic>, 53 from <italic>Tetraodon </italic>(<italic>T. nigroviridis </italic>and <italic>T. rubripes</italic>), 45 from <italic>Drosophila melanogaster</italic>, 27 from <italic>Caenorhabditis elegans</italic>, 6 from <italic>Saccharomyces cerevisae</italic>, and 10 from <italic>Schizosaccharomyces pombe </italic>[see Additional file ##SUPPL##8##8##]. The catalogue contains all classes of JmjC domain containing proteins that have been described so far, including "JmjC domain only" proteins. To classify these proteins into evolutionary subgroups and to investigate their evolutionary relationship we analysed the selected sequences. A multiple sequence alignment and subsequent bootstrapping analysis was carried out using ClustalW (1000 bootstrap trials, 111 seeds) and PHYLIP. Based on the similarity of the JmjC domain, the phylogenetic analysis revealed 25 evolutionary conserved subgroups [see Additional file ##SUPPL##7##6##]. These subgroups are in general accordance to previous studies, which reported the classification of 44 JmjC domain containing proteins into 4 related clusters [##REF##16715513##5##] and of 98 JmjC domain containing proteins into 7 subfamilies [##REF##16983801##4##]. Although, our analysis shows that especially the previously classified subgroup of "JmjC domain only family" [##REF##16983801##4##] splits into 8 further subgroups based on the sequence divergence of their JmjC domain (1110034B05Rik family, Hif1an, Hspbap1, Pla2g4b, Jmjd5, Jmjd4, Jmjd6, JFP6-like, and No66/Mina53 subgroups). As previously described [##REF##16983801##4##,##REF##16715513##5##], these proteins have beside their JmjC domain no other recognisable protein domain but they are clearly functionally distinct from each other regarding reported protein interaction partners, cellular localisation and physiological functions. For example, Hif1an is an asparagine protein hydroxylase for the HIF1α transcription factor, Hspbap1 is a cytosolic protein that associates with heat shock protein 27 kD, and Jmjd6, Mina53, No66 are distinct nuclear proteins with as yet unknown molecular functions. In addition to these groups, we included in our analysis the AlkB, PH4/Phy-3, Ogfod2, Ogfod1/Tpa1, Plod 1–3, Lepre1/Leprel 1–2, 1110031I02Rik, and Egln1–3/Egl-9/Hph subgroups that were not analysed before together with the recently identified JmjC domain containing histone demethylase proteins [##REF##16983801##4##]. Interestingly, we found that all subgroups could be classified independent of additional domain architectural features of the analysed proteins. As previously reported [##REF##16983801##4##,##REF##16715513##5##], the domain compositions of individual proteins classified into a given subgroup varied, indicating that neither their enzymatic specificity nor their functional relationships can be predicted solely based on the analysis of primary protein sequences. This is further supported by the identification of certain members within subgroups that obviously can not possess enzymatic functions as 2OG-dependent dioxygenases due to sequence changes within their predicted Fe(II) and/or 2OG-binding sites. For example, the Hairless protein within the Jmjd1/Hairless subgroup, the Epe1 protein in the Jhd1 subgroup and the whole Jarid2 protein subfamily contain substitutions in these binding sites indicating that these proteins are unlikely to be functional enzymes ([##REF##16983801##4##] and Figure ##FIG##7##8##).</p> <p>Interestingly, in some subgroups we observed a huge expansion in the number of family members. The P4H/Phy-3 family of prolyl-4-hydroxylases contained 22 proteins from <italic>Drosophila melanogaster</italic>. The P4H/Phy-3 subgroup compromises prolyl-4-hydroxylases involved in collagen biosynthesis or hypoxia sensing. Currently, it is unknown why especially the fly has so many members in this protein family. Another example is the AlkB subfamily of DNA and RNA demethylating enzymes. This group contains 45 proteins that split into the Alkbh1 to Alkbh8 protein family members. Apart from Alkbh1 to Alkbh3 that are known to be DNA/RNA repair enzymes [##REF##12226667##52##,##REF##16174769##78##,##REF##16482161##79##] the exact biochemical activity of most of the other subfamily members has still to be determined.</p> <p>Following phylogenetic analysis and tree generation a subtree encompassing the Jmjd6 protein was extracted for further analysis (Figure ##FIG##7##8##). By combining the information from the phylogenomic analysis with the annotation of protein subfamilies it is evident that Jmjd4 is closely related to Jmjd6 regarding sequence similarity. The JmjC domain of both proteins shows the highest sequence similarity among all 313 classified members (32% sequence identity and 50% sequence similarity). This indicates that <italic>Jmjd4 </italic>is the closest paralogue of <italic>Jmjd6</italic>.</p> <p>Moreover, a detailed sequence analysis of JmjC domain containing protein subfamilies that are related to the Jmjd6 and Jmjd4 subgroups revealed further insights into Fe (II) and 2OG coordinating mechanisms (Figure ##FIG##7##8##). The JFP6 subgroup has a HxHxH motif instead of the Hx(D/E)x<sub>n</sub>H motif that is commonly used by all other JmjC domain containing proteins for the coordination of Fe(II) (Figure ##FIG##7##8##, orange box). If this motif is compatible with an oxygenase activity, as one would expect, has to be investigated in the future. Furthermore, the Jmjd5 homologue of <italic>Tetraodon nigroviridis</italic>, all Phf2 proteins, and Epe1 have only a Hx(D/E) motif (red box in Figure ##FIG##7##8##) and therefore lack the second histidine, which is known to be important for iron binding. Moreover, the <italic>T</italic>. <italic>nigroviridis </italic>Phf2 homologue lacks in addition the lysine that stabilises the C5'-caboxylate of 2OG. Next to this motif, Tuskada et al. described residues in positions -3, +6–8 and +13–15, in respect to the HxD motif to be important for the catalytic function of histone demethylases [##REF##16362057##21##]. In the -3 position, the 11100034B05Rik, JFP6 and No66/Mina subgroups have a nonpolar side chain present in contrast to the polar side chains of serine/threonine residues in other subgroups. In the +6–8 position, members of the Pla2g4b, FBXL10/FBXL11, Phf2/Phf8, and individual members of the JFP6 and Jhd1 subgroup have a conserved tyrosine. In Epe1, this tyrosine has been shown to be important for heterochromatin destabilisation, which is independent of an Epe1 histone demethylation activity [##REF##16362057##21##,##REF##16762840##46##]. Another interesting feature is the predicted presence of an extended major β-sheet in most of the JmjC domain only protein families encompassing the subtree (Figure ##FIG##7##8##).</p>
<title>Discussion</title> <p>Here we report a systematic annotation and comparison of Jmjd6 homologous proteins across diverse phyla. With the recent description of the JmjC domain as a novel demethylase signature motif [##REF##16362057##21##] and the discovery that many proteins containing JmjC domains are involved in reversing histone methylation marks [##REF##16983801##4##,##REF##17342184##80##], this protein superfamily has attracted central interest in many different fields of biology. It has become quite clear that these enzymes have not only important functions in the epigenetic regulation of chromatin modification states but also in diverse other physiological contexts ranging from oxygen sensing to fatty acid metabolism and DNA repair [##REF##17301803##81##]. However, for most of these proteins it has remained difficult to unravel their molecular function or to predict their enzymatic targets solely based on sequence comparisons and analysis of their JmjC domains. In the past, we have been interested in the "phosphatidylserine receptor" (now renamed into jumonji domain containing 6), a JmjC domain containing protein, which was initially described to function as a type II transmembrane receptor for the engulfment of apoptotic cells [##REF##10811223##56##]. To compare the JmjC domain of Jmjd6 to the DSBH/JmjC domains of other members within this protein superfamily we designed a semiautomated sequence retrieval and analysis routine that allowed us to establish a catalogue of 313 DSBH/JmjC domain containing proteins in nine selected species. In addition, we analysed the <italic>Jmjd6 </italic>and neighbouring loci in vertebrates and investigated their modes of expression.</p> <title>Comparative analysis of the vertebrate <italic>Jmjd6 </italic>locus and immediate neighbouring loci</title> <p>With the ongoing genome sequencing projects in many different species and integrated efforts to annotate the functional content of fully sequenced genomes there is a need to complement existing automated genome annotations with a focused bioinformatics analysis of individual genes and gene families to make fully use of publicly available databases and resources. Such approaches might help to correct mistakes in automated annotation procedures and might add additional information to gene sets. This will finally lead to a more comprehensive and accurate annotation of gene functions.</p> <p>We used the PipMaker alignment tool to analyse the patterns of sequence conservation across the mouse <italic>Jmjd6 </italic>and the neighbouring loci <italic>1110005A03Rik </italic>and <italic>Sfrs2</italic>. To identify all protein coding exons as well as highly conserved non-coding sequence elements in this region we analysed genomic sequences of 22 different vertebrate species. This led to the identification of an up to now overseen additional <italic>Jmjd6 </italic>exon in mammals and to the detection of a conserved putative bi-directional promoter element in between the closely spaced <italic>Jmjd6 </italic>and <italic>1110005A03Rik </italic>genes. The re-annotated genomic region encompassing the <italic>Jmjd6</italic>, <italic>1110005A03Rik </italic>and <italic>Sfrs2 </italic>loci was then the basis for a comparison of gene-targeting strategies that have been used to generate three different loss-of-function alleles of <italic>Jmjd6 </italic>in the mouse [##REF##15345036##54##,##REF##14645847##62##,##REF##14715629##63##]. By aligning a map of restriction enzyme sites used for the construction of the different targeting vectors with the genomic region encompassing the three genes, we could show that in the <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Flv </italic></sup>allele and in the <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Ysfk </italic></sup>allele exon 1 of the <italic>111005A03Rik </italic>gene as well as the conserved intergenic region in between <italic>Jmjd6 </italic>and <italic>111005A03Rik </italic>was deleted upon homologous recombination in ES cells. An expression analysis of <italic>Jmjd6</italic>, <italic>111005A03Rik </italic>and <italic>Sfsr2 </italic>in the three different mouse alleles confirmed that <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Flv </italic></sup>and <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Ysfk </italic></sup>mutant mice are double knockouts of the <italic>Jmjd6 </italic>and <italic>111005A03Rik </italic>genes. In contrast, the <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Gbf </italic></sup>allele ablates only the expression of <italic>Jmjd6</italic>. Although all three independently generated knockout mice show to some degree overlapping phenotypes regarding defects during embryogenesis (e.g. disturbed eye development, brain malformations and defects in haematopoiesis) other controversial observations concerning the description of specific phenotypes have been reported [##REF##15345036##54##,##REF##14645847##62##,##REF##14715629##63##]. Foremost, the capabilities of the mutant mouse lines to engulf apoptotic cells have been described differently. While Kunisaki et al. [##REF##14715629##63##] and Li et al. [##REF##14645847##62##] reported that their knockout mice showed impairment in the clearance of apoptotic cells during embryogenesis and defects in apoptotic cell phagocytosis assays, no such differences in phagocytosis competence could be detected in <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Gbf </italic></sup>mutant mice [##REF##15345036##54##]. It has been suggested that differences in genetic backgrounds might account for these reported phenotypic variations [##UREF##0##82##, ####REF##15453906##83##, ##REF##17471263##84####17471263##84##]. However, in the light of the described results with <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Flv </italic></sup>and <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Ysfk </italic></sup>mutant mice being double knockouts of <italic>Jmjd6 </italic>and <italic>111005A03Rik </italic>the situation might be far more complex. It could well be that loss-of <italic>111005A03Rik </italic>function in <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Flv </italic></sup>and <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Ysfk </italic></sup>knockout mice might contribute in addition to the <italic>Jmjd6</italic>-deficiency to the phenotypic discrepancies observed. Future studies involving the generation and analysis of <italic>111005A03Rik </italic>knockout mice will address this possibility.</p> <p>Interestingly, the as yet functionally uncharacterised <italic>111005A03Rik </italic>gene contains a predicted arginine protein methyltransferase domain. We demonstrate here, that the conserved sequence element in the short intergenic region between <italic>Jmjd6 </italic>and <italic>111005A03Rik </italic>possess bi-directional promoter activities for both genes. Of note, the forward orientation of the fragment, which controls <italic>111005A03Rik </italic>expression <italic>in vivo </italic>showed only half of the transcriptional activity as compared to the reverse orientation, which controls <italic>Jmjd6 </italic>expression <italic>in vivo</italic>. This indicates that both genes might compose a bi-directional transcriptional unit with opposing transcriptional and regulatory activities. This gives room for speculations. If <italic>111005A03Rik </italic>has a function as a methyltransferase and the Jmjd6 protein possess demethylation activity as many other JmjC domain containing proteins are doing, both proteins might be involved in a common biochemical pathway targeting possibly a shared substrate. Interestingly, such counterbalancing principles of methylation and demethylation activity have been already suggested for other JmjC domain containing proteins [##REF##11276424##53##]. For example, the Alkbh8 protein contains an amino-terminal JmjC domain, which might possess DNA demethylase activity that is fused to a predicted C-terminal methylase domain. This domain architecture is conserved in the homologous Alkbh8 proteins in <italic>Drosophila </italic>(CG17807) and <italic>C. elegans </italic>(C14B1.10) suggesting the selection of a controlled methylation/demethylation regulatory mechanism during the evolution of these proteins.</p> <p>If <italic>Jmjd6 </italic>and <italic>111005A03Rik </italic>also establish such counteracting functions needs to be addressed in the future. To proof such a hypothesis, the molecular function of both genes need first to be characterised. Nevertheless, the additional deletion of the <italic>111005A03Rik </italic>gene in two of the three described <italic>Jmjd6 </italic>knockout mice might well explain differences reported for their apoptotic cell engulfment phenotypes. If <italic>111005A03Rik </italic>is functioning as a methyltransferase in the control of transcriptional regulation, it is possible that the expression of other genes targeted or regulated by 111005A03Rik and/or Jmjd6 might influence phagocytosis of apoptotic cells. In such a scenario the described engulfment defects might be related to secondary effects caused by deregulation of genes involved in processes of phagocytosis or cellular differentiation of phagocytes.</p> <p>Using PipMaker analysis based on sequence conservation we identified a before overseen exon in <italic>Jmjd6 </italic>loci of mammalian species. A subsequently undertaken e2g annotation of all available EST sequences to the <italic>Jmjd6/JMJD6 </italic>genomic regions lead to the discovery of two new alternative transcripts in mice and humans. At the protein level the most distal C-terminal part of the Jmjd6 protein (last 100 amino acids) is affected in all new identified transcripts. Either truncated polypeptides are generated that miss the poly-serine stretch, the putative AT-hook domain, and the predicted sumoylation and nuclear export signal (NES) or a stretch of eight amino acid is exchanged with 31 other residues. In this second splice variant identified in mice (Mm3) the predicted sumoylation site, the NES and the poly-serine stretch are retained but the putative AT-hook domain is missing (Figure ##FIG##2##3C##). Our 3D structural model of the Jmjd6 protein indicates that the sumoylation site, the NES and the AT-hook motif might be accessible for interacting proteins. The AT-hook domain is known to mediate binding to DNA in high-mobility group proteins [##REF##9742243##85##]. Jmjd6 is a nuclear protein, which has been shown to contain five independently active nuclear localisation signals [##REF##14729065##58##]. The NES might allow a shuttling of the protein from the nucleus to the cytoplasm. Likewise, the sumoylation site might be used for the regulation of interactions with other proteins.</p> <p>A comparison of the Jmjd6 sequence with the resolved protein structure of the Hif1an/HIF1α C-terminal transactivation domain co-complex indicates that interactions with other proteins might occur in the Jmjd6 C-terminal region [##REF##12446723##49##]. Alternatively, at this site homo-di- or homo-multimerization of Jmjd6 polypeptides might take place. For many JmjC/DSBH domain containing enzymes it is known that they can exist in various oligomeric forms. For example, Hif1an is a dimer in solution, germin forms a hexamer and HPPE a tetramer (for review see [##REF##16513174##6##]). Indeed, very recently it could be demonstrated that Jmjd6 forms covalent multimers through homomeric polymerisation of the 50 kD protein [##REF##17534701##61##]. Nuclear aggregated Jmjd6 products with masses of 100, 150, 200 and 250 kD were detected after immunoprecipitation indicating that assembly of homomultimers might be important for the protein function [##REF##17534701##61##]. Thus, the new alternative transcripts encoded by the <italic>Jmjd6 </italic>locus might possess a regulatory role by suppressing multimer formation or interactions with other proteins. Such possibilities need to be experimentally investigated in the future.</p> <title>Structural model of the Jmjd6 protein and evolutionary relationship to other JmjC/DSBH-domain containing proteins</title> <p>Recent crystallographic studies have helped to provide important insights in the structure and functions of JmjC domain containing proteins. More than ten structures of different members of this large oxygenase family have been solved demonstrating that the JmjC domain forms a characteristic DSBH structural fold composed of eight β-strands that build two β-sheets each containing four anti-parallel β-strands [##REF##16513174##6##]. The catalytic core of these metalloenzymes consists of a highly conserved Hx(E/D)x<sub>n</sub>H signature motif that supplies three chelating positions for the Fe(II) iron and two additional chelating positions are contributed by the substrate co-factor 2OG. Very often the binding sites for 2OG in the DSBH core are as well conserved among the members of this protein family. Due to high sequence similarity of Jmjd6 to other JmjC domain containing proteins [##REF##15193161##57##], we have embarked on the modelling of the Jmjd6 3D-structure to gain further insights into possible protein functions. Our analysis shows that Jmjd6 is indeed most likely a DSBH-fold containing 2OG-dependent dioxygenase enzyme, as previously first suggested by Cikala and colleagues [##REF##15193161##57##]. All residues needed for Fe(II) coordination, 2OG co-substrate binding and oxidative catalysis were found to be located within the predicted DSBH domain of the protein model and are evolutionary preserved in all analysed Jmjd6 homologous proteins probably including those of non-metazoan species. We identified in all orthologous Jmjd6 proteins from metazoan species the characteristic HxDx<sub>[n]</sub>H catalytical triad and our model predicts five conserved residues (Trp174, Thr184, Asn197, Lys204, and Thr285) as being important for interaction with a 2OG co-substrate. In addition, the extreme conservation of residues forming the eight β-strands of the DSBH domain suggests a spatial and functional conservation of the Jmjd6 JmjC domain during evolution. This makes it highly likely that all Jmjd6 proteins in the different animal kingdoms are indeed active 2OG-dependent dioxygenase enzymes. This might not be true for other members of the JmjC/DSBH fold superfamily of proteins. For example, the subfamily of JFP6-like proteins has a substitution in the Fe(II) coordinating residues changing the Fe(II) binding motif. Furthermore, hairless and Jarid2 proteins are most likely enzymatically inactive due to the missing HxD motifs in their DSBH protein cores [##REF##16983801##4##]. This suggests that JmjC domains might possess additional functions that are independent of any enzymatic activity.</p> <p>2OG-dependent dioxygenases are known to catalyse two types of posttranslational protein modifications. Described have been hydroxylation of the amino acid side chains of proline, asparagines, aspartic acid and tryptophan, and secondly demethylation reactions on lysine residues (for review see [##REF##16513174##6##,##REF##17301803##81##]). It has been proposed that also methylated arginines might be targets of JmjC domain containing proteins. In fact, most recently it has been suggested that Jmjd6 might possess an enzymatic activity as an arginine histone demethylase. It was reported that Jmjd6 is capable to catalyse the removal of di-methyl groups on histone H3 at arginine 2 (H3R2) and on histone 4 at arginine 3 (H4R3) [##REF##17947579##86##]. Asymmetric di-methylation of H3R2 has been associated with repression of transcriptional activity through its enrichment in heterochromatic loci and inactive euchromatic genes [##REF##17898714##87##,##REF##17898715##88##]. Thus, Jmjd6 might have a function in gene activation by erasing repressive chromatin marks. But it is also not unlikely that additional, non-histone substrates are targeted by Jmjd6.</p> <p>All types of posttranslational modifications catalysed by 2OG-dependent dioxygenases involve a common type of hydroxylation reaction. During the catalytic reaction, the splitting of molecular oxygen is coupled to the hydroxylation of the target substrate and the oxidative decarboxylation of 2OG to give succinate and CO<sub>2</sub>. In the case of demethylation, hydroxylation of a methyl group is followed immediately by deformylation producing formaldehyde as an additional reaction product. As noted before by others, it is impossible to predict the prime substrate of 2OG-dependent dioxygenase based on a primary sequence analysis of the JmjC domain [##REF##17342184##80##,##REF##17301803##81##]. Even JmjC domain containing enzymes of the same subfamily such as the JHDM1A and the JHDM2A histone demethylases that both recognise H3K36 methylation have very limited primary or secondary structure similarity [##REF##16677698##22##]. This suggests a requirement for an extended substrate recognition mechanism that probably involves protein-protein interactions in larger protein complexes or for a necessity of very defined interaction sites in the catalytic core that are important for determining substrate selectivity [##REF##17589501##89##]. This is further supported by a number of recent studies that implicated other methyl-dependent protein domains like Tudor, MBT, PHD fingers and chromo domains as reader domains in the recognition of the histone code [##REF##16415788##90##, ####REF##16728978##91##, ##REF##16728977##92##, ##REF##16728974##93##, ##REF##16728976##94####16728976##94##]. Jmjd6 contains a single JmjC domain and beside the predicted recognition sequences for nuclear in- and export, sumoylation and the AT-hook DNA-binding motif no other recognisable domain is present in the protein. Thus, the substrate specificity of Jmjd6 must be regulated by interaction with other proteins.</p> <p>What is quite clear is that the putative enzymatic function of Jmjd6 as a 2OG-dependent dioxygenase is linked to its important role in development and morphogenesis, at least in vertebrates. All mouse mutants conclusively point to an essential role of <italic>Jmjd6 </italic>during embryogenesis were the gene is involved in the differentiation of multiple tissues and organs. Also in zebrafish, the morpholino knockdown of the <italic>jmjd6 </italic>orthologue resulted in severe defects in somite, brain, heart and notochord development [##REF##15469976##95##]. This may or may not involve secondary effects on apoptosis induction and cell corpse clearance. Unlikely is a direct involvement of <italic>jmjd6 </italic>in apoptotic cell recognition or engulfment. Given that many 2OG-dependent dioxygenases have been implicated in the regulation of multiple pathways, global effects on gene expression and metabolism are not surprising when their enzymatic functions get perturbed or ablated.</p> <p>The function of Jmjd6 during development seems to be different in vertebrates and invertebrates. Inactivation of the orthologous genes in <italic>Drosophila </italic>and <italic>C. elegans </italic>were surprisingly not associated with embryonic lethality and no major developmental defects have been reported for the loss-of-function mutants [##REF##17522160##59##,##REF##14645848##96##]. Interestingly, mutant flies that lack <italic>dPSR </italic>show inappropriate apoptosis during eye development [##REF##17522160##59##]. Major phenotypes were observed in gain-of-function experiments in which the <italic>dPSR </italic>gene was overexpressed during development. Male transgenic flies that expressed <italic>dPSR </italic>under ubiquitous promoters developed a genital rotation phenotype. <italic>dPSR </italic>overexpression in the wing resulted in various wing blade fusion phenotypes and driving gene expression in the eye induced a small rough eye phenotype. Further genetic analysis pointed to an involvement of <italic>dPSR </italic>in the <italic>head involution defective </italic>(<italic>hid</italic>) and c-Jun-NH2 terminal kinase pathways. These suggested that loss of <italic>dPSR </italic>enhances apoptosis in the developing eye, while overexpression of <italic>dPSR </italic>protects from apoptosis. Of note, the analysis of mutant flies did not reveal a role for <italic>dPSR </italic>in the phagocytosis of apoptotic cells during development. Altogether, these observations support a function of <italic>dPSR </italic>in cell survival and differentiation that might have parallels in the vertebrate system [##REF##17522160##59##].</p> <p>When we compared the protein sequences of vertebrate and invertebrate orthologous Jmjd6 proteins we found divergence in the N- and C-terminal regions of the protein. Interestingly, the C-terminal poly-serine stretch is either shorter or missing in the invertebrate species. Moreover, we could not identify an additional exon in the invertebrate orthologous <italic>Jmjd6 </italic>genes that would generate alternative transcripts with effects on the C-terminal part of the protein. Therefore, it could well be that differences in this region of the protein are linked in vertebrates with different physiological functions of <italic>Jmjd6 </italic>protein, possibly through interactions with other proteins. This is also supported by identification of proteins in plants, fungi and protista that show partial sequence homology to the animalia Jmjd6 protein (Figure ##FIG##4##5## and [see Additional file ##SUPPL##3##4##]). These were identified through sequence search routines focussing on a conserved N-terminal sequence stretch outside of the highly conserved JmjC domain. Consequently, these proteins show in their N-terminal parts significant sequence homology to the mouse Jmjd6 protein, but are diverse to Jmjd6 in their C-terminal regions. Moreover, many of these proteins are longer then the animalia Jmjd6 proteins (up to 2.000 aa instead of around 400 aa for animalia Jmjd6 proteins) and contain often in addition to the JmjC domain a single F-box domain.</p> <p>Thus, it is very likely that they possess different or additional functions than the animalia Jmjd6 proteins and a possible homology of these proteins to Jmjd6 should be treated with caution. As the reciprocal high similarity suggests a relation to Jmjd6, it might be possible to postulate an ancestral "Jmjd6" gene that encoded a larger protein with additional domains (like F-box domains), which then have been lost during evolution of <italic>Jmd6 </italic>genes in animalia. Alternatively, the evolution of <italic>Jmjd6 </italic>genes in animalia could have been an independent event and these genes are not related in any way to the proteins we identified in plants, fungi and protista. In the future it might be possible to consolidate on one hypothesis, given that more functional data about these hypothetical proteins are being gathered and higher quality sequence information becomes available for phylogenetic analysis.</p>
<title>Conclusion</title> <p>Taken together, our integrated bioinformatic analysis provided new insights into possible structure-function relationships of Jmjd6 proteins and on the evolution of the whole superfamily of DSBH/JmjC domain containing proteins. Of note, the observations made on the protein model are not experimentally proven at the current state and await confirmation when the crystallographic structure of Jmjd6 has been obtained. In addition, we provide first evidence for a counteracting transcriptional regulative mechanism involving the co-localised <italic>Jmjd6 </italic>and <italic>1110005A03Rik </italic>genes.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>The <italic>jumonji C (JmjC) domain containing gene 6 </italic>(<italic>Jmjd6</italic>, previously known as phosphatidylserine receptor) has misleadingly been annotated to encode a transmembrane receptor for the engulfment of apoptotic cells. Given the importance of JmjC domain containing proteins in controlling a wide range of diverse biological functions, we undertook a comparative genomic analysis to gain further insights in <italic>Jmjd6 </italic>gene organisation, evolution, and protein function.</p> <title>Results</title> <p>We describe here a semiautomated computational pipeline to identify and annotate JmjC domain containing proteins. Using a sequence segment N-terminal of the Jmjd6 JmjC domain as query for a reciprocal BLAST search, we identified homologous sequences in 62 species across all major phyla. Retrieved <italic>Jmjd6 </italic>sequences were used to phylogenetically analyse corresponding loci and their genomic neighbourhood. This analysis let to the identification and characterisation of a bi-directional transcriptional unit compromising the <italic>Jmjd6 </italic>and <italic>1110005A03Rik </italic>genes and to the recognition of a new, before overseen <italic>Jmjd6 </italic>exon in mammals. Using expression studies, two novel <italic>Jmjd6 </italic>splice variants were identified and validated <italic>in vivo</italic>. Analysis of the <italic>Jmjd6 </italic>neighbouring gene <italic>1110005A03Rik </italic>revealed an incident deletion of this gene in two out of three earlier reported <italic>Jmjd6 </italic>knockout mice, which might affect previously described conflicting phenotypes. To determine potentially important residues for <italic>Jmjd6 </italic>function a structural model of the Jmjd6 protein was calculated based on sequence conservation. This approach identified a conserved double-stranded β<sup>-</sup>helix (DSBH) fold and a HxDx<sub>n</sub>H facial triad as structural motifs. Moreover, our systematic annotation in nine species identified 313 DSBH fold-containing proteins that split into 25 highly conserved subgroups.</p> <title>Conclusion</title> <p>We give further evidence that <italic>Jmjd6 </italic>most likely has a function as a nonheme-Fe(II)-2-oxoglutarate-dependent dioxygenase as previously suggested. Further, we provide novel insights into the evolution of Jmjd6 and other related members of the superfamily of JmjC domain containing proteins. Finally, we discuss possibilities of the involvement of <italic>Jmjd6 </italic>and <italic>1110005A03Rik </italic>in an antagonistic biochemical pathway.</p>
<title>Authors' contributions</title> <p>PH and JB designed the computational pipeline flow together. The domain based multiple sequence alignment analysis of all protein sequences, the phylogeny prediction there from and the PipMaker sequence analysis as well as the analysis of the different mouse mutant <italic>Jmjd6 </italic>alleles was carried out by JB. PH contributed the full-length Jmjd6 analysis, 3D structure prediction and developed the new method for conservation based residue shading. Furthermore, PH identified and validated new splice variants of <italic>Jmjd6</italic>, characterised the bi-directional <italic>Jmjd6 </italic>transcriptional unit and developed SMART-parser as well as additional Perl scripts for automated data processing. SE carried out parts of the <italic>Jmjd6 </italic>expression analysis. AL coordinated the study and wrote the paper. All authors read and approved the final manuscript.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>We would like to thank Ivonne Wegener for excellent technical assistance. Dr. Richard Flavell (Yale University, New Haven, USA) and Dr. Yoshinori Fukui (Kyushu University, Fukuoka, Japan) for making available <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Flv </italic></sup>knockout mice and RNA from <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Ysfk </italic></sup>mutant embryos, respectively. We are thankful to Dr. Bastian Pasche (HZI, Braunschweig) for giving critical comments on the manuscript. This work was supported by a grant from Deutsche Forschungsgemeinschaft (DFG LE 2130/2-1) to AL.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>Computational pipeline flow used for the retrieval and analysis of DSBH/JmjC-domain containing proteins</bold>. It contains three parallel procedures. Indicated in yellow colour are searches and analysis approaches based on DNA sequences, data processing approaches based on protein sequences are coloured green, and phylogenomic analysis methods based on DSBH protein domains are shown in blue. Initial queries are depicted as circles, databases as ovals, candidate gene and protein lists as cut-off boxes, and methods/analysis tools as normal boxes, respectively. Barrels indicate Perl script based analysis of retrieved results. Candidate lists of the individual workflows are shown and the final results of the analysis are depicted in bold. Arrows indicate the flow of information and the connections between different search routines and methods. Sources of queried databases and the analysis software used are given in the respective box elements.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>Percent identity plot of the chromosomal regions encompassing the <italic>Jmjd6 </italic>locus in vertebrate species</bold>. The genomic regions of eight vertebrate species identified as putative <italic>Jmjd6 </italic>orthologous loci were aligned with the mouse <italic>Jmjd6 </italic>locus and its neighbouring genes using PipMaker. Green, blue and yellow overlays highlight the conservation of the annotated exons/coding sequences of the individual genes found in the reference sequence. The orange overlay indicates conserved sequences of the first exons of the <italic>Jmjd6 </italic>and <italic>1110005A03Rik </italic>genes and sequences outside of the coding elements. The red overlay marks a high-scoring segment in all analysed mammalian species supporting the presence of at least one additional exon (No. 5) that is not included in the current Ensembl gene annotation (build 44). Transcriptional orientations of the genes, their exons, UTRs, repetitive elements, and CpG islands are in the symbol key at the figure bottom.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold>Identification and characterisation of new splice variants of <italic>Jmjd6</italic></bold>. <bold><italic>(A) </italic></bold>Two additional putative exons (green boxes) were identified in addition to the transcript annotated in Ensembl (variant 1) by screening all public EST databases. Alternative splicing using these two exons results in the generation of two additional transcripts of <italic>Jmjd6 </italic>(variants 2 and 3, respectively). Half arrows and orange lines in the schematic presentation of the transcripts highlight the combination of exons detected by RT-PCR. Numbering of <italic>Jmjd6 </italic>exons includes the additional identified exon variants 4b and 5. <bold><italic>(B) </italic></bold>Expression analysis of alternative <italic>Jmjd6 </italic>transcripts in adult mouse organs and embryonic stages (E) using RT-PCR confirmed the presence of all three transcripts at different mRNA expression levels in the analysed samples. Half arrows in (A) indicate the position of the PCR primers used for the amplification. Blue numbers and arrows in (B) point to PCR fragments amplified from splice variants 1 to 3. Amplification of the housekeeping <italic>β-actin </italic>gene was used as a loading control. <bold><italic>(C) </italic></bold>The effects of the different splicing events on the C-termini are shown in the gapped alignment from amino acid 300 onwards. The full length human JMJD6 and murine Jmjd6 proteins are represented in the upper part of the figure. Nearly all novel splice variants identified in the mouse (Mm 2) and in humans (Hs 2–3) are truncated at the C-terminus in comparison to splice variant 1 (Mm 1 and Hs 1) and are predicted to contain no poly-serine stretch ([Ser]<sub>n</sub>). In contrast, the JmjC domain is not affected (green box). Predicted AT-hook domain, sumoylation recognition site (Sumo), and nuclear export signals (NES) are annotated as depicted at the bottom.</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p><bold>Identification of the <italic>Jmjd6 </italic>– <italic>1110005A03Rik </italic>bi-directional transcriptional unit</bold>. <bold><italic>(A) </italic></bold>Schematic overview of the 65 bp non-coding and non-overlapping intergenic region between the murine <italic>Jmjd6 </italic>and <italic>1110005A03Rik </italic>(<italic>Rik</italic>) genes. The conserved sequence element identified using PipMaker in Figure 2 is shown in orange. 5'-untranslated regions (UTRs) are indicated in grey, and first exon coding sequences of the <italic>Jmjd6 </italic>gene in green and of the <italic>1110005A03Rik </italic>gene in blue, respectively. The region indicated by two half arrows was used for promoter activity assays presented in (B). It was amplified and cloned to yield an 840 bp test fragment <bold><italic>(B) </italic></bold>This putative bi-directional promoter fragment was used in a luciferase reporter gene assay to measure transcriptional activity. Basal promoter activity was found for both transcriptional orientations compared to a promoterless control (pGL3) vector. The 5'-3' orientation of the fragment (<italic>1110005A03Rik </italic>orientation, <italic>Rik</italic>) yielded only half of the transcriptional activity as it was measured for the 3'-5' orientation (<italic>Jmjd6 </italic>orientation, <italic>Jmjd6</italic>). The values reported for the transfection experiments are means ± standard deviation of three independent triplicate experiments. P value was determined using Student's t-test. *** P &lt; 0.001.</p></caption></fig> <fig position="float" id="F5"><label>Figure 5</label><caption><p><bold>Multiple sequence alignment of homologous Jmjd6 proteins</bold>. The sequences of 54 homologous Jmjd6 proteins were used for the calculation of a multiple sequence alignment (ClustalW). The numbering of the residues is according to the mouse protein sequence. Identical residues were shaded in blue and similar residues in magenta. Green shading indicates the predicted catalytic residues of the HxDx<sub>n</sub>H facial. Sequence insertions in individual species outside of the conserved core regions of the mouse Jmjd6 protein were masked out and the maximum numbers of amino acids in these insertions are indicated [n]. Based on the sequence conservation, a frequency corrected sequence logo was calculated (shown above the alignment) and the mouse sequence conservation was highlighted in colour dots according to their degree of conservation (shown as colour code in percent identity below the alignment). The JmjC domain (light blue bar) and additional identified secondary structure elements (red cylinders: α-helices; yellow and green arrows: β-strands) were annotated to the sequence as indicated at the bottom and the β-strands forming the jelly-roll of the JmjC domain were numbered (I-VIII). Green arrows represent the β-sheet that extends the major sheet outside of the jelly-roll.</p></caption></fig> <fig position="float" id="F6"><label>Figure 6</label><caption><p><bold>Predicted structure of Jmjd6 based on comparative modelling (stereoview)</bold>. Jmjd6 sequence alignments (Fig. 5) and HHpred were used to model Jmjd6 utilising the PDB structures of Hif1an (<ext-link ext-link-type="pdb" xlink:href="1h2k">1h2k</ext-link>), Putative asparaginyl hydroxylase – 2636534 (<ext-link ext-link-type="pdb" xlink:href="1vrb">1vrb</ext-link>) and Jmjd2a (<ext-link ext-link-type="pdb" xlink:href="2gp5">2gp5</ext-link>), respectively. <bold><italic>(A) </italic></bold>Stereoview of the predicted structure of Jmjd6 presented as a ribbon model. The conserved eight-stranded DSBH core found in all Fe(II) and 2-oxoglutarate (2OG)-dependent dioxygenases is coloured in blue. Additional β-strands attached to the major β-sheet are shown in orange. Helices are depicted in green. <bold><italic>(B) </italic></bold>View of the predicted active site of the Jmjd6-Fe(II)-2OG complex showing coordination of Fe(II) to 2OG, His187, Asp189, and His273. 2OG also ligates to Trp174, Asn197 and Lys204 with Thr285 stabilising Asn197 and Trp174. Additional important residues for putative interactions are shown in cyan and include hydrophobic interactions from Phe133 and Val275 as well as Thr184, which is involved in 2OG-binding in Hif1an. Interacting residues along with the 2OG co-substrate are shown as stick presentations, putative H-bond interactions as dotted lines. Fe(II) is depicted as an orange ball.</p></caption></fig> <fig position="float" id="F7"><label>Figure 7</label><caption><p><bold>Evolutional conservation of critical residues in Jmjd6 (stereoview)</bold>. <bold><italic>(A) </italic></bold>Ribbon presentation of the predicted Jmjd6 structure is coloured based on the sequence homology observed in 54 Jmjd6 homologous proteins (Fig. 5). <bold><italic>(B) </italic></bold>Annotation of the sequence identity onto the catalytic domain demonstrated the conservation of the critical residues within the active site of the Jmjd6-Fe(II)-2OG complex. Residues were shaded according to the heat map with dark red representing total conservation of residues (100% sequence identity across 54 analysed species), and blue colour representing 50% sequence identity.</p></caption></fig> <fig position="float" id="F8"><label>Figure 8</label><caption><p><bold>Phylogenetic relationship of Jmjd6 to other JmjC domain containing proteins</bold>. Detailed section of the phylogenetic tree of 313 JmjC proteins that is presented [in additional file ##SUPPL##7##6##]. Shown are 12 out of 25 identified DSBH fold containing protein subgroups as indicated with protein family names on the right side. Scale bar represents the relative phylogenetic distance as determined with PHYLIP. Bootstrap values are shown for values &lt;1000. The first column of the table shows the residues potentially involved in iron binding. The second, third and fourth columns show potentially important residues within the -3, +6–8 and +13–15 regions in respect to the Hx(D/E/H) motif, respectively. Residues discussed in the text are highlighted in red, orange, blue, green and yellow. The fifth column shows the predicted ability to build an extended major sheet with three antiparallel β-strands as determined by Jpred. A "+" indicates the presence of three beta strands, a "?" two beta-strands and a "-" no predicted beta-strands within a 50 aa window starting 80 aa before the first β-strand from the JmjC jelly-roll. H.s. = <italic>Homo sapiens</italic>, M.m. = <italic>Mus musculus</italic>, D.r. = <italic>Danio rerio</italic>, T.n. = <italic>Tetraodon nigroviridis</italic>, T.r. = <italic>Takifugu rubripes</italic>, D.m. = <italic>Drosophila melanogaster</italic>, C.e. = <italic>Caenorhabditis elegans</italic>, S.c. = <italic>Saccharomyces cerevisiae</italic>, S.p. = <italic>Schizosaccharomyces pombe</italic>.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Species distribution and chromosomal locations of <italic>Jmjd6 </italic>orthologs in mammalian and non-mammalian genomes</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>Species name</bold></td><td align="left"><bold>Common name</bold></td><td align="center"><bold>Genomic location</bold></td><td align="center"><bold>Cluster with <italic>1110005A03Rik</italic></bold></td><td align="center"><bold>Cluster with <italic>Sfrs2</italic></bold></td></tr></thead><tbody><tr><td align="left"><italic>Homo sapiens</italic></td><td align="left">Human</td><td align="center">17</td><td align="center"><bold>+</bold></td><td align="center"><bold>+</bold></td></tr><tr><td align="left"><italic>Pan troglodytes</italic></td><td align="left">Chimp</td><td align="center">17</td><td align="center"><bold>+</bold></td><td align="center"><bold>+</bold></td></tr><tr><td align="left"><italic>Macaca mulatta</italic></td><td align="left">Macaque</td><td align="center">16</td><td align="center"><bold>+</bold></td><td align="center"><bold>+</bold></td></tr><tr><td align="left"><italic>Mus musculus</italic></td><td align="left">Mouse</td><td align="center">11</td><td align="center"><bold>+</bold></td><td align="center"><bold>+</bold></td></tr><tr><td align="left"><italic>Rattus norvegicus</italic></td><td align="left">Rat</td><td align="center">10</td><td align="center"><bold>+</bold></td><td align="center"><bold>+</bold></td></tr><tr><td align="left"><italic>Cavia porcellus</italic></td><td align="left">Guniea pig</td><td align="center">Gene scaffold 591</td><td align="center"><bold>+</bold></td><td/></tr><tr><td align="left"><italic>Spermophilus tridecemlineatus</italic></td><td align="left">Ground squirrel</td><td align="center">Scaffold 5110</td><td align="center"><bold>+</bold></td><td/></tr><tr><td align="left"><italic>Tupaia belangeri</italic></td><td align="left">Tree shrew</td><td align="center">Scaffold 5346</td><td align="center"><bold>+</bold></td><td/></tr><tr><td align="left"><italic>Canis familiaris</italic></td><td align="left">Dog</td><td align="center">18</td><td align="center"><bold>+</bold></td><td align="center"><bold>+</bold></td></tr><tr><td align="left"><italic>Felis catus</italic></td><td align="left">Cat</td><td align="center">Scaffold 295</td><td align="center"><bold>+</bold></td><td/></tr><tr><td align="left"><italic>Bos taurus</italic></td><td align="left">Cow</td><td align="center">19</td><td align="center"><bold>+</bold></td><td align="center"><bold>+</bold></td></tr><tr><td align="left"><italic>Myotis lucifugus</italic></td><td align="left">Microbat</td><td align="center">Scaffold 462</td><td align="center"><bold>+</bold></td><td/></tr><tr><td align="left"><italic>Erinaceus europaeus</italic></td><td align="left">Hedgehoc</td><td align="center">Scaffold 670</td><td align="center"><bold>+</bold></td><td align="center"><bold>+</bold></td></tr><tr><td align="left"><italic>Dasypus novemcinctus</italic></td><td align="left">Armadillo</td><td align="center">Scaffold 584</td><td align="center"><bold>+</bold></td><td/></tr><tr><td align="left"><italic>Echinops telfairi</italic></td><td align="left">Tenrec</td><td align="center">Scaffold 620</td><td align="center"><bold>+</bold></td><td/></tr><tr><td align="left"><italic>Monodelphis domestica</italic></td><td align="left">Opossum</td><td align="center">2</td><td align="center"><bold>+</bold></td><td align="center"><bold>+</bold></td></tr><tr><td align="left"><italic>Ornithorhynchus anatinus</italic></td><td align="left">Platypus</td><td align="center">Ultra 430</td><td align="center"><bold>+</bold></td><td/></tr><tr><td align="left"><italic>Gallus gallus</italic></td><td align="left">Chicken</td><td align="center">18</td><td align="center"><bold>+</bold></td><td/></tr><tr><td align="left"><italic>Xenopus tropicalis</italic></td><td align="left">Western clawed frog</td><td align="center">Scaffold 178</td><td align="center"><bold>+</bold></td><td align="center"><bold>+</bold></td></tr><tr><td align="left"><italic>Gasterosteus aculeatus</italic></td><td align="left">Stickleback</td><td align="center">Group XI</td><td align="center"><bold>+</bold></td><td/></tr><tr><td align="left"><italic>Oryzias latipes</italic></td><td align="left">Medaka</td><td align="center">Scaffold 1014</td><td/><td/></tr><tr><td align="left"><italic>Tetraodon nigroviridis</italic></td><td align="left">Tetraodon</td><td align="center">SCAF 13844</td><td align="center"><bold>+</bold></td><td/></tr><tr><td align="left"><italic>Takifugu rubripes</italic></td><td align="left">Fugu</td><td align="center">Scaffold 519</td><td align="center"><bold>+</bold></td><td/></tr><tr><td align="left"><italic>Danio rerio</italic></td><td align="left">Zebrafish</td><td align="center">3</td><td align="center"><bold>+</bold></td><td align="center"><bold>+</bold></td></tr><tr><td align="left"><italic>Ciona intestinalis</italic></td><td align="left">Sea squirt</td><td align="center">3p</td><td/><td/></tr><tr><td align="left"><italic>Ciona savignyi</italic></td><td align="left">Sea squirt</td><td align="center">Reftig 1</td><td/><td/></tr><tr><td align="left"><italic>Drosophila melanogaster</italic></td><td align="left">Fruitfly</td><td align="center">3R</td><td/><td/></tr><tr><td align="left"><italic>Aedes aegypti</italic></td><td align="left">Aedes</td><td align="center">Super contig 1.887</td><td/><td/></tr><tr><td align="left"><italic>Anopheles gambiae</italic></td><td align="left">Mosquito</td><td align="center">X</td><td/><td/></tr><tr><td align="left"><italic>Caenorhabditis elegans</italic></td><td align="left">C. elegans</td><td align="center">IV</td><td/><td/></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p>Validation of <italic>JMJD6 </italic>splices variants in humans. RT-PCR expression analysis of <italic>JMJD6 </italic>splice variants in a human cDNA tissue panel obtained from the BioChain Institute (Hayward, USA). Arrows indicate amplicons from correctly predicted splice variants that were confirmed by sequencing. Multiplex PCR amplification results in one verified unspecific PCR fragment migrating at approximately 400 bp. Amplification of the housekeeping gene <italic>Ribosomal Protein 9 </italic>(<italic>RPS9</italic>) was used as loading control.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title><p>Absence of nonsense-mediated decay in <italic>Jmjd6 </italic>splice variants. Two additional putative exons (green boxes) were identified <bold>(A)</bold>. Alternative splicing using these two exons results in the generation of two additional transcripts of Jmjd6 (variants 2 and 3, respectively). Half arrows and orange lines in the schematic presentation of the transcripts highlight the combination of exons detected by RT-PCR in <bold>(B) </bold>and <bold>(C)</bold>. Primers were designed to bind to the respective exons shown. <bold>(B) </bold>Experimental validation of the two new predicted alternative Jmjd6 transcripts using RT-PCR and agarose gel electrophoresis. Expression analyses confirmed that both splice variants are not subject to nonsense-mediated decay in adult mouse organs and <bold>(C) </bold>human tissues. Amplification of the housekeeping β-actin gene was used as a RNA loading control.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional file 3</title><p>Comparison of strategies used to inactivate the <italic>Jmjd6 </italic>gene in mouse ES-cells. Schematic representation of the mouse locus according to the Ensembl annotation of the <italic>Jmjd6 </italic>gene and its neighbouring genes <bold><italic>(A)</italic></bold>. Arrows indicate transcriptional orientation of the genes along chromosome 11, from centromer to telomer. <italic>Jmjd6 </italic>exons are shown in green, exons of <italic>1110005A03Rik </italic>in blue, and of <italic>Sfrs2 </italic>in yellow. Individual exons are numbered. The strategies used to generate the different targeted knock out alleles are shown. Red boxes indicate gene regions deleted by homologous recombination. The <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Flv </italic></sup>allele was generated by replacing an <italic>Aat</italic>II/<italic>Spe</italic>I fragment by a neomycin resistance cassette. In the <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Ysfk </italic></sup>allele the <italic>Eco</italic>RI fragment and in the <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Gbf </italic></sup>allele the <italic>Aat</italic>II/<italic>Rsr</italic>II fragment were replaced respectively. Based on these chosen fragments, targeted inactivation of <italic>Jmjd6 </italic>seems to affect in two of the three alleles the neighbouring <italic>1110005A03Rik </italic>locus. <bold><italic>(B) </italic></bold>RT-PCR analysis of the expression of the <italic>Jmjd6</italic>, <italic>1110005A03Rik </italic>(<italic>Rik</italic>) and <italic>Sfrs2 </italic>genes in <italic>Jmjd6 </italic>wild type and <italic>Jmjd6 </italic>homozygous mutant embryos. The expression of <italic>1110005A03Rik </italic>is not detectable in the <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Flv </italic></sup>allele and in the <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Ysfk </italic></sup>allele, whereas expression of <italic>Sfrs2 </italic>is unaltered in all three mouse lines investigated <bold><italic>(C) </italic></bold>Semi-quantitative RT-PCR analysis of <italic>1110005A03Rik </italic>expression in <italic>Jmjd6</italic><sup><italic>tm</italic>1<italic>Gbf </italic></sup>wild type and homozygous mutant embryos (embryonic stage E13.5) and in fetal-liver derived macrophages shows that the level of expression is not altered due to the targeted inactivation of the <italic>Jmjd6 </italic>gene in this mouse line. RT-PCR expression analysis was performed in differential dilution steps of cDNA material (1, 1:10, 1:100). In all RT-PCR expression experiments amplification of the housekeeping gene <italic>β-actin </italic>severed as loading control.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S4"><caption><title>Additional file 4</title><p>Phylogenetic tree of live for identified putative Jmjd6 orthologs and their domain composition. Scientific classification of species with identified putative Jmjd6 orthologs was performed according to the NCBI taxonomy browser [##UREF##10##108##]. Schematic presentation of the grouping and categorisation of 62 species with identified Jmjd6 proteins. The individual species are shown on the left side, underlying colours highlight the kingdoms – animalia (red), plantae (green), protista (blue), fungi (yellow), and eubacteria (magenta), boxes on the right side represent the domain. The intermediate boxes show informative higher order ranks (e.g. genus, family, order, class, or phyla). On the left side, the domain composition for each protein is given. Length of polypeptides in amino acids (aa) are indicated on the top. Asterisks indicate partial or truncated protein sequences.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S9"><caption><title>Additional file 9</title><p>Jmjd6 homologous proteins. Catalogue of 62 identified Jmjd6 homologous proteins.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S7"><caption><title>Additional file 7</title><p>Features of structurally resolved 2OG-Fe(II)-dependent dioxygenases in comparison to Jmjd6. Table describing features of structurally resolved and experimentally characterised 2OG-Fe(II)-dependent dioxygenases in comparison to Jmjd6.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S5"><caption><title>Additional file 5</title><p>Evolutionary conservation of surface exposed residues in the Jmjd6 protein model (stereoview). The surface model was computed using PyMol and coloured according to the sequence conservation code below in percent sequence identity in 54 analysed species. The coordinating Fe(II) of the catalytic triad is depicted in black, the co-substrate 2OG is shown as stick representation.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S6"><caption><title>Additional file 6</title><p>Phylogenetic relationship of Jmjd6 to the superfamily of JmjC-domain containing proteins. In total 313 sequences of DSBH/JmjC domain containing proteins from humans and model organisms (mouse, zebrafish, pufferfish, nematode, fruit fly, and yeast) were analysed for their relationship to different subgroups of JmjC domain containing proteins. JmjC domain sequences of the proteins were analysed by multiple sequence alignment and subsequent bootstrap analysis using ClustalW (1000 bootstrap trials, 111 seeds) and PHYLIP (ProML). The resulting unrouted tree based on sequence similarity within the JmjC domain revealed 25 evolutionary conserved protein subgroups. These are highlighted in different colours according to the colour code on the left side. Five proteins that could not be grouped into subgroups are indicated in red letters. H.s. = <italic>Homo sapiens</italic>, M.m. = <italic>Mus musculus</italic>, D.r. = <italic>Danio rerio</italic>, T.n. = <italic>Tetraodon nigroviridis</italic>, T.r. = <italic>Takifugu rubripes</italic>, D.m. = <italic>Drosophila melanogaster</italic>, C.e. = <italic>Caenorhabditis elegans</italic>, S.c. = <italic>Saccharomyces cerevisiae</italic>, S.p. = <italic>Schizosaccharomyces pombe</italic>.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S8"><caption><title>Additional file 8</title><p>Catalogue of JmjC/DSBH domain containing proteins. Table of 313 JmjC/DSBH domain containing protein used for the phylogenetic analysis.</p></caption></supplementary-material>
<table-wrap-foot><p>Chromosomal locations have been taken from ENSEMBL, + indicates the presences of <italic>1110005A03Rik </italic>or <italic>Sfrs2 </italic>in the immediate genomic neighbourhood of <italic>Jmjd6</italic>.</p></table-wrap-foot>
<graphic xlink:href="1471-2164-9-293-1"/> <graphic xlink:href="1471-2164-9-293-2"/> <graphic xlink:href="1471-2164-9-293-3"/> <graphic xlink:href="1471-2164-9-293-4"/> <graphic xlink:href="1471-2164-9-293-5"/> <graphic xlink:href="1471-2164-9-293-6"/> <graphic xlink:href="1471-2164-9-293-7"/> <graphic xlink:href="1471-2164-9-293-8"/>
<media xlink:href="1471-2164-9-293-S1.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-293-S2.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-293-S3.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-293-S4.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-293-S9.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-293-S7.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-293-S5.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-293-S6.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-293-S8.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Weitzman"], "given-names": ["JB"], "article-title": ["The curious world of apoptotic cell clearance"], "source": ["J Biol"], "year": ["2004"], "volume": ["3"], "fpage": ["13"], "pub-id": ["10.1186/jbiol13"]}, {"article-title": ["RepeatMasker"]}, {"article-title": ["NCBI-BLAST"]}, {"article-title": ["SMART"]}, {"article-title": ["Swiss-Pdb-Viewer"]}, {"article-title": ["PyMOL"]}, {"article-title": ["PHYLIP"]}, {"article-title": ["SUMOplot"]}, {"article-title": ["NetNES"]}, {"article-title": ["Jpred-3"]}, {"article-title": ["NCBI-Taxonomy"]}]
BMC Genomics. 2008 Jun 18; 9:293