Novel ITPA variants identified by whole genome sequencing and RNA sequencing.

Approximately 80% of rare diseases have a genetic cause, and an accurate genetic diagnosis is necessary for disease management, prognosis prediction, and genetic counseling. Whole-exome sequencing (WES) is a cost-effective approach for exploring the genetic cause, but several cases often remain undiagnosed. We combined whole genome sequencing (WGS) and RNA sequencing (RNA-seq) to identify the pathogenic variants in an unsolved case using WES. RNA-seq revealed aberrant exon 4 and exon 6 splicing of ITPA. WGS showed a previously unreported splicing donor variant, c.263+1G>A, and a novel heterozygous deletion, including exon 6. Detailed examination of the breakpoint indicated the deletion caused by recombination between Alu elements in different introns. The proband was found to have developmental and epileptic encephalopathies caused by variants in the ITPA gene. The combination of WGS and RNA-seq may be effective in diagnosing conditions in proband who could not be diagnosed using WES.

Strategic validation of variants of uncertain significance in ECHS1 genetic testing.

BACKGROUND: Enoyl-CoA hydratase short-chain 1 (ECHS1) is an enzyme involved in the metabolism of branched chain amino acids and fatty acids. Mutations in the ECHS1 gene lead to mitochondrial short-chain enoyl-CoA hydratase 1 deficiency, resulting in the accumulation of intermediates of valine. This is one of the most common causative genes in mitochondrial diseases. While genetic analysis studies have diagnosed numerous cases with ECHS1 variants, the increasing number of variants of uncertain significance (VUS) in genetic diagnosis is a major problem. METHODS: Here, we constructed an assay system to verify VUS function for ECHS1 gene. A high-throughput assay using ECHS1 knockout cells was performed to index these phenotypes by expressing cDNAs containing VUS. In parallel with the VUS validation system, a genetic analysis of samples from patients with mitochondrial disease was performed. The effect on gene expression in cases was verified by RNA-seq and proteome analysis. RESULTS: The functional validation of VUS identified novel variants causing loss of ECHS1 function. The VUS validation system also revealed the effect of the VUS in the compound heterozygous state and provided a new methodology for variant interpretation. Moreover, we performed multiomics analysis and identified a synonymous substitution p.P163= that results in splicing abnormality. The multiomics analysis complemented the diagnosis of some cases that could not be diagnosed by the VUS validation system. CONCLUSIONS: In summary, this study uncovered new ECHS1 cases based on VUS validation and omics analysis; these analyses are applicable to the functional evaluation of other genes associated with mitochondrial disease.

Involvement of GLCCI1 in mouse spermatogenesis.

Spermatid production is a complex regulatory process in which coordination between hormonal control and apoptosis plays a pivotal role in maintaining a balanced number of sperm cells. Apoptosis in spermatogenesis is controlled by pro-apoptotic and anti-apoptotic molecules. Hormones involved in the apoptotic process during spermatogenesis include gonadotrophins, sex hormones, and glucocorticoid (GC). GC acts broadly as an apoptosis inducer by binding to its receptor (glucocorticoid receptor: GR) during organ development processes, such as spermatogenesis. However, the downstream pathway induced in GC-GR signaling and the apoptotic process during spermatogenesis remains poorly understood. We reported previously that GC induces full-length glucocorticoid-induced transcript 1 (GLCCI1-long), which functions as an anti-apoptotic mediator in thymic T cell development. Here, we demonstrate that mature murine testis expresses a novel isoform of GLCCI1 protein (GLCCI1-short) in addition to GLCCI1-long. We demonstrate that GLCCI1-long is expressed in spermatocytes along with GR. In contrast, GLCCI1-short is primarily expressed in spermatids where GR is absent; instead, the estrogen receptor is expressed. GLCCI1-short also binds to LC8, which is a known mediator of the anti-apoptotic effect of GLCCI1-long. A luciferase reporter assay revealed that ß-estradiol treatment synergistically increased Glcci1-short promotor-driven luciferase activity in Erα-overexpressing cells. Together with the evidence that the conversion of testosterone to estrogen is preceded by aromatase expression in spermatids, we hypothesize that estrogen induces GLCCI1-short, which, in turn, may function as a novel anti-apoptotic mediator in mature murine testis.

Severe spinal cord hypoplasia due to a novel ATAD3A compound heterozygous deletion.

Biallelic deletions extending into the ATPase family AAA-domain containing protein 3A (ATAD3A) gene lead to infantile lethality with severe pontocerebellar hypoplasia (PCH). However, only 12 such cases have been reported worldwide to date, and the genotype-phenotype correlations are not well understood. We describe cases associated with the same novel biallelic deletions of the ATAD3A and ATAD3B/3A regions in Japanese siblings with severe spinal cord hypoplasia and multiple malformations, including PCH, leading to neonatal death. The ATAD3A protein is essential for normal interaction between mitochondria and endoplasmic reticulum and is important for mitochondrial biosynthesis. The cases were evaluated using whole-genome sequencing for genetic diagnosis of mitochondrial disease. Spinal cord lesions associated with biallelic compound heterozygous deletion extending into the ATAD3A gene have not been reported. In addition, the ATAD3A deletion was 19 base pairs long, which is short compared with those reported previously. This deletion introduced a frameshift, resulting in a premature termination codon, and was expected to be a null allele. The pathological findings of the atrophic spinal cord showed gliosis and tissue destruction of the gray and white matter. We describe spinal cord lesions as a new central nervous system phenotype associated with a biallelic compound heterozygous deletion extending into the ATAD3A gene. Biallelic ATAD3A deletions should be considered in cases of mitochondrial disease with spinal cord hypoplasia and PCH.

Impact of measuring heteroplasmy of a pathogenic mitochondrial DNA variant at the single-cell level in individuals with mitochondrial disease.

Pathogenic mitochondrial DNA heteroplasmy has mainly been assessed with bulk sequencing in individuals with mitochondrial disease. However, the distribution of heteroplasmy at the single-cell level in skin fibroblasts obtained from individuals, together with detailed clinical and biochemical information, remains to be investigated. We used the mitochondrial DNA single-cell assay for the transposase-accessible chromatin sequencing method. Skin fibroblasts were obtained from six individuals with mitochondrial disease and pathogenic m.3243A>G variants of differing severity. Different distributions of heteroplasmy at the single-cell level were identified in skin fibroblasts from all six individuals. Four individuals with different outcomes showed similar averaged heteroplasmy rates with normal mitochondrial respiratory chain enzyme activity, while the distribution of single-cell heteroplasmy patterns differed among the individuals. This study showed different heteroplasmy distribution patterns at the single-cell level in individuals with the m.3243A>G variant, who had a similar averaged heteroplasmy rates with normal mitochondrial respiratory chain enzyme activity. Whether such different heteroplasmy distribution patterns explain the different clinical outcomes should be assessed further in future studies. Measuring heteroplasmy of pathogenic mitochondrial DNA variants at the single-cell level could be important in individuals with mitochondrial disease.

Neonatal-onset mitochondrial disease: clinical features, molecular diagnosis and prognosis.

OBJECTIVE: Neonatal-onset mitochondrial disease has not been fully characterised owing to its heterogeneity. We analysed neonatal-onset mitochondrial disease in Japan to clarify its clinical features, molecular diagnosis and prognosis. DESIGN: Retrospective observational study from January 2004 to March 2020. SETTING: Population based. PATIENTS: Patients (281) with neonatal-onset mitochondrial disease diagnosed by biochemical and genetic approaches. INTERVENTIONS: None. MAIN OUTCOME MEASURES: Disease types, initial symptoms, biochemical findings, molecular diagnosis and prognosis. RESULTS: Of the 281 patients, multisystem mitochondrial disease was found in 194, Leigh syndrome in 26, cardiomyopathy in 38 and hepatopathy in 23 patients. Of the 321 initial symptoms, 236 occurred within 2 days of birth. Using biochemical approaches, 182 patients were diagnosed by mitochondrial respiratory chain enzyme activity rate and 89 by oxygen consumption rate. The remaining 10 patients were diagnosed using a genetic approach. Genetic analysis revealed 69 patients had nuclear DNA variants in 36 genes, 11 of 15 patients had mitochondrial DNA variants in five genes and four patients had single large deletion. The Cox proportional hazards regression analysis showed the effects of Leigh syndrome (HR=0.15, 95% CI 0.04 to 0.63, p=0.010) and molecular diagnosis (HR=1.87, 95% CI 1.18 to 2.96, p=0.008) on survival. CONCLUSIONS: Neonatal-onset mitochondrial disease has a heterogenous aetiology. The number of diagnoses can be increased, and clarity regarding prognosis can be achieved by comprehensive biochemical and molecular analyses using appropriate tissue samples.

Development of Leigh syndrome with a high probability of cardiac manifestations in infantile-onset patients with m.14453G > A.

The m.14453G > A mutation in MT-ND6 has been described in a few patients with mitochondrial encephalomyopathy, lactic acidosis, and stroke-like episodes or Leigh syndrome.However, the clinical spectrum and molecular characteristics are unclear.Here, we present four infantile-onset patients with m.14453G > A-associated Leigh syndrome. All four patients had brainstem lesions with basal ganglia lesions, and two patients had cardiac manifestations. Decreased ND6 protein expression and immunoreactivity were observed in patient-derived samples. There was no clear correlation between heteroplasmy levels and onset age or between heteroplasmy levels and phenotype; however, infantile onset was associated with Leigh syndrome.

Genome sequencing and RNA-seq analyses of mitochondrial complex I deficiency revealed Alu insertion-mediated deletion in NDUFV2.

Isolated complex I deficiency is the most common cause of pediatric mitochondrial disease. Exome sequencing (ES) has revealed many complex I causative genes. However, there are limitations associated with identifying causative genes by ES analysis. In this study, we performed multiomics analysis to reveal the causal variants. We here report two cases with mitochondrial complex I deficiency. In both cases, ES identified a novel c.580G>A (p.Glu194Lys) variant in NDUFV2. One case additionally harbored c.427C>T (p.Arg143*), but no other variants were observed in the other case. RNA sequencing showed aberrant exon splicing of NDUFV2 in the unsolved case. Genome sequencing revealed a novel heterozygous deletion in NDUFV2, which included one exon and resulted in exon skipping. Detailed examination of the breakpoint revealed that an Alu insertion-mediated rearrangement caused the deletion. Our report reveals that combined use of transcriptome sequencing and GS was effective for diagnosing cases that were unresolved by ES.

Long-term prognosis and genetic background of cardiomyopathy in 223 pediatric mitochondrial disease patients.

BACKGROUND: Cardiomyopathy is a risk factor for poor prognosis in pediatric patients with mitochondrial disease. However, other risk factors including genetic factors related to poor prognosis in mitochondrial disease has yet to be fully elucidated. METHODS AND RESULTS: Between January 2004 and September 2019, we enrolled 223 consecutive pediatric mitochondrial disease patients aged <18 years with a confirmed genetic diagnosis, including 114 with nuclear gene mutations, 89 patients with mitochondrial DNA (mtDNA) point mutations, 11 with mtDNA single large-scale deletions and 9 with chromosomal aberrations. Cardiomyopathy at baseline was observed in 46 patients (21%). Hazard ratios (HR) and 95% confidence intervals (CI) were calculated for all-cause mortality. Over a median follow-up of 36 months (12-77), there were 85 deaths (38%). The overall survival rate was significantly lower in patients with cardiomyopathy than in those without (p < 0.001, log-rank test). By multivariable analysis, left ventricular (LV) hypertrophy (HR = 4.6; 95% CI: 2.8-7.3), neonatal onset (HR = 2.9; 95% CI: 1.8-4.5) and chromosomal aberrations (HR = 2.9; 95% CI: 1.3-6.5) were independent predictors of all-cause mortality. Patients with LV hypertrophy with neonatal onset and/or chromosomal aberrations had higher mortality (100% in 21 patients) than those with LV hypertrophy alone (71% in 14 patients). CONCLUSION: In pediatric patients with mitochondrial disease, cardiomyopathy was common (21%) and was associated with increased mortality. LV hypertrophy, neonatal onset and chromosomal aberrations were independent predictors of all-cause mortality. Prognosis is particularly unfavorable if LV hypertrophy is combined with neonatal onset and/or chromosomal aberrations.

Clinical heterogeneity in patients with m.4412G > A MT-TM mutation and different heteroplasmy levels.

The identification of the m.4412G > A MT-TM (mt-tRNAMet) mutation was first reported in 2019. The affected individual presented with childhood-onset seizures and myopathy and bilateral basal ganglia changes, with heteroplasmy levels in muscle as high as 90%. Here, we describe another adult-onset patient with the same mutation and additional phenotypes, including hearing impairment, cerebellar ataxia, progressive dementia, and myopathy. The 10% heteroplasmy level observed in skin fibroblasts from this patient are lower than those in the previously reported patient. Our report suggests possible clinical heterogeneity in patients with mitochondrial tRNA mutations based on heteroplasmy levels.

Fatal perinatal mitochondrial cardiac failure caused by recurrent de novo duplications in the ATAD3 locus.

BACKGROUND: In about half of all patients with a suspected monogenic disease, genomic investigations fail to identify the diagnosis. A contributing factor is the difficulty with repetitive regions of the genome, such as those generated by segmental duplications. The ATAD3 locus is one such region, in which recessive deletions and dominant duplications have recently been reported to cause lethal perinatal mitochondrial diseases characterized by pontocerebellar hypoplasia or cardiomyopathy, respectively. METHODS: Whole exome, whole genome and long-read DNA sequencing techniques combined with studies of RNA and quantitative proteomics were used to investigate 17 subjects from 16 unrelated families with suspected mitochondrial disease. FINDINGS: We report six different de novo duplications in the ATAD3 gene locus causing a distinctive presentation including lethal perinatal cardiomyopathy, persistent hyperlactacidemia, and frequently corneal clouding or cataracts and encephalopathy. The recurrent 68 Kb ATAD3 duplications are identifiable from genome and exome sequencing but usually missed by microarrays. The ATAD3 duplications result in the formation of identical chimeric ATAD3A/ATAD3C proteins, altered ATAD3 complexes and a striking reduction in mitochondrial oxidative phosphorylation complex I and its activity in heart tissue. CONCLUSIONS: ATAD3 duplications appear to act in a dominant-negative manner and the de novo inheritance infers a low recurrence risk for families, unlike most pediatric mitochondrial diseases. More than 350 genes underlie mitochondrial diseases. In our experience the ATAD3 locus is now one of the five most common causes of nuclear-encoded pediatric mitochondrial disease but the repetitive nature of the locus means ATAD3 diagnoses may be frequently missed by current genomic strategies. FUNDING: Australian NHMRC, US Department of Defense, Japanese AMED and JSPS agencies, Australian Genomics Health Alliance and Australian Mito Foundation.

A homozygous variant in NDUFA8 is associated with developmental delay, microcephaly, and epilepsy due to mitochondrial complex I deficiency.

Mitochondrial complex I deficiency is caused by pathogenic variants in mitochondrial and nuclear genes associated with complex I structure and assembly. We report the case of a patient with　NDUFA8-related mitochondrial disease. The patient presented with developmental delay, microcephaly, and epilepsy. His fibroblasts showed apparent biochemical defects in mitochondrial complex I. Whole-exome sequencing revealed that the patient carried a homozygous variant in NDUFA8. His fibroblasts showed a reduction in the protein expression level of not only NDUFA8, but also the other complex I subunits, consistent with assembly defects. The enzyme activity of complex I and oxygen consumption rate were restored by reintroducing wild-typeNDUFA8 cDNA into patient fibroblasts. The functional properties of the variant in NDUFA8 were also investigated using NDUFA8 knockout cells expressing wild-type or mutated NDUFA8 cDNA. These experiments further supported the pathogenicity of the variant in complex I assembly. This is the first report describing that the loss of NDUFA8, which has not previously been associated with mitochondrial disease, causes severe defect in the assembly of mitochondrial complex I, leading to progressive neurological and developmental abnormalities.

ANISEED 2019: 4D exploration of genetic data for an extended range of tunicates.

ANISEED (https://www.aniseed.cnrs.fr) is the main model organism database for the worldwide community of scientists working on tunicates, the vertebrate sister-group. Information provided for each species includes functionally-annotated gene and transcript models with orthology relationships within tunicates, and with echinoderms, cephalochordates and vertebrates. Beyond genes the system describes other genetic elements, including repeated elements and cis-regulatory modules. Gene expression profiles for several thousand genes are formalized in both wild-type and experimentally-manipulated conditions, using formal anatomical ontologies. These data can be explored through three complementary types of browsers, each offering a different view-point. A developmental browser summarizes the information in a gene- or territory-centric manner. Advanced genomic browsers integrate the genetic features surrounding genes or gene sets within a species. A Genomicus synteny browser explores the conservation of local gene order across deuterostome. This new release covers an extended taxonomic range of 14 species, including for the first time a non-ascidian species, the appendicularian Oikopleura dioica. Functional annotations, provided for each species, were enhanced through a combination of manual curation of gene models and the development of an improved orthology detection pipeline. Finally, gene expression profiles and anatomical territories can be explored in 4D online through the newly developed Morphonet morphogenetic browser.

High-Throughput Protein Production Combined with High- Throughput SELEX Identifies an Extensive Atlas of Ciona robusta Transcription Factor DNA-Binding Specificities.

Transcription factors (TFs) control gene transcription, binding to specific DNA motifs located in cis-regulatory elements across the genome. The identification of TF-binding motifs is thus an important aspect to understand the role of TFs in gene regulation. SELEX, Systematic Evolution of Ligands by EXponential enrichment, is an efficient in vitro method, which can be used to determine the DNA-binding specificity of TFs. Thanks to the development of high-throughput (HT) DNA cloning system and protein production technology, the classical SELEX assay has be extended to high-throughput scale (HT-SELEX).We report here the detailed protocol for the cloning, production, and purification of 420 Ciona robusta DNA BD. 263 Ciona robusta TF DNA-binding domain proteins were purified in milligram quantities and analyzed by HT-SELEX. The identification of 139 recognition sequences generates an atlas of protein-DNA-binding specificities that is crucial for the understanding of the gene regulatory network (GRN) of Ciona robusta. Overall, our analysis suggests that the Ciona robusta repertoire of sequence-specific transcription factors comprises less than 500 genes. The protocols for high-throughput protein production and HT-SELEX described in this article for the study of Ciona robusta TF DNA-binding specificity are generic and have been successfully applied to a wide range of TFs from other species, including human, mouse, and Drosophila.

The interaction landscape between transcription factors and the nucleosome.

Nucleosomes cover most of the genome and are thought to be displaced by transcription factors in regions that direct gene expression. However, the modes of interaction between transcription factors and nucleosomal DNA remain largely unknown. Here we systematically explore interactions between the nucleosome and 220 transcription factors representing diverse structural families. Consistent with earlier observations, we find that the majority of the studied transcription factors have less access to nucleosomal DNA than to free DNA. The motifs recovered from transcription factors bound to nucleosomal and free DNA are generally similar. However, steric hindrance and scaffolding by the nucleosome result in specific positioning and orientation of the motifs. Many transcription factors preferentially bind close to the end of nucleosomal DNA, or to periodic positions on the solvent-exposed side of the DNA. In addition, several transcription factors usually bind to nucleosomal DNA in a particular orientation. Some transcription factors specifically interact with DNA located at the dyad position at which only one DNA gyre is wound, whereas other transcription factors prefer sites spanning two DNA gyres and bind specifically to each of them. Our work reveals notable differences in the binding of transcription factors to free and nucleosomal DNA, and uncovers a diverse interaction landscape between transcription factors and the nucleosome.

ANISEED 2017: extending the integrated ascidian database to the exploration and evolutionary comparison of genome-scale datasets.

ANISEED (www.aniseed.cnrs.fr) is the main model organism database for tunicates, the sister-group of vertebrates. This release gives access to annotated genomes, gene expression patterns, and anatomical descriptions for nine ascidian species. It provides increased integration with external molecular and taxonomy databases, better support for epigenomics datasets, in particular RNA-seq, ChIP-seq and SELEX-seq, and features novel interactive interfaces for existing and novel datatypes. In particular, the cross-species navigation and comparison is enhanced through a novel taxonomy section describing each represented species and through the implementation of interactive phylogenetic gene trees for 60% of tunicate genes. The gene expression section displays the results of RNA-seq experiments for the three major model species of solitary ascidians. Gene expression is controlled by the binding of transcription factors to cis-regulatory sequences. A high-resolution description of the DNA-binding specificity for 131 Ciona robusta (formerly C. intestinalis type A) transcription factors by SELEX-seq is provided and used to map candidate binding sites across the Ciona robusta and Phallusia mammillata genomes. Finally, use of a WashU Epigenome browser enhances genome navigation, while a Genomicus server was set up to explore microsynteny relationships within tunicates and with vertebrates, Amphioxus, echinoderms and hemichordates.

Impact of cytosine methylation on DNA binding specificities of human transcription factors.

The majority of CpG dinucleotides in the human genome are methylated at cytosine bases. However, active gene regulatory elements are generally hypomethylated relative to their flanking regions, and the binding of some transcription factors (TFs) is diminished by methylation of their target sequences. By analysis of 542 human TFs with methylation-sensitive SELEX (systematic evolution of ligands by exponential enrichment), we found that there are also many TFs that prefer CpG-methylated sequences. Most of these are in the extended homeodomain family. Structural analysis showed that homeodomain specificity for methylcytosine depends on direct hydrophobic interactions with the methylcytosine 5-methyl group. This study provides a systematic examination of the effect of an epigenetic DNA modification on human TF binding specificity and reveals that many developmentally important proteins display preference for mCpG-containing sequences.

Myt1l safeguards neuronal identity by actively repressing many non-neuronal fates.

Normal differentiation and induced reprogramming require the activation of target cell programs and silencing of donor cell programs. In reprogramming, the same factors are often used to reprogram many different donor cell types. As most developmental repressors, such as RE1-silencing transcription factor (REST) and Groucho (also known as TLE), are considered lineage-specific repressors, it remains unclear how identical combinations of transcription factors can silence so many different donor programs. Distinct lineage repressors would have to be induced in different donor cell types. Here, by studying the reprogramming of mouse fibroblasts to neurons, we found that the pan neuron-specific transcription factor Myt1-like (Myt1l) exerts its pro-neuronal function by direct repression of many different somatic lineage programs except the neuronal program. The repressive function of Myt1l is mediated via recruitment of a complex containing Sin3b by binding to a previously uncharacterized N-terminal domain. In agreement with its repressive function, the genomic binding sites of Myt1l are similar in neurons and fibroblasts and are preferentially in an open chromatin configuration. The Notch signalling pathway is repressed by Myt1l through silencing of several members, including Hes1. Acute knockdown of Myt1l in the developing mouse brain mimicked a Notch gain-of-function phenotype, suggesting that Myt1l allows newborn neurons to escape Notch activation during normal development. Depletion of Myt1l in primary postmitotic neurons de-repressed non-neuronal programs and impaired neuronal gene expression and function, indicating that many somatic lineage programs are actively and persistently repressed by Myt1l to maintain neuronal identity. It is now tempting to speculate that similar 'many-but-one' lineage repressors exist for other cell fates; such repressors, in combination with lineage-specific activators, would be prime candidates for use in reprogramming additional cell types.

DNA-dependent formation of transcription factor pairs alters their binding specificity.

Gene expression is regulated by transcription factors (TFs), proteins that recognize short DNA sequence motifs. Such sequences are very common in the human genome, and an important determinant of the specificity of gene expression is the cooperative binding of multiple TFs to closely located motifs. However, interactions between DNA-bound TFs have not been systematically characterized. To identify TF pairs that bind cooperatively to DNA, and to characterize their spacing and orientation preferences, we have performed consecutive affinity-purification systematic evolution of ligands by exponential enrichment (CAP-SELEX) analysis of 9,400 TF-TF-DNA interactions. This analysis revealed 315 TF-TF interactions recognizing 618 heterodimeric motifs, most of which have not been previously described. The observed cooperativity occurred promiscuously between TFs from diverse structural families. Structural analysis of the TF pairs, including a novel crystal structure of MEIS1 and DLX3 bound to their identified recognition site, revealed that the interactions between the TFs were predominantly mediated by DNA. Most TF pair sites identified involved a large overlap between individual TF recognition motifs, and resulted in recognition of composite sites that were markedly different from the individual TF's motifs. Together, our results indicate that the DNA molecule commonly plays an active role in cooperative interactions that define the gene regulatory lexicon.

Conservation of transcription factor binding specificities across 600 million years of bilateria evolution.

Divergent morphology of species has largely been ascribed to genetic differences in the tissue-specific expression of proteins, which could be achieved by divergence in cis-regulatory elements or by altering the binding specificity of transcription factors (TFs). The relative importance of the latter has been difficult to assess, as previous systematic analyses of TF binding specificity have been performed using different methods in different species. To address this, we determined the binding specificities of 242 Drosophila TFs, and compared them to human and mouse data. This analysis revealed that TF binding specificities are highly conserved between Drosophila and mammals, and that for orthologous TFs, the similarity extends even to the level of very subtle dinucleotide binding preferences. The few human TFs with divergent specificities function in cell types not found in fruit flies, suggesting that evolution of TF specificities contributes to emergence of novel types of differentiated cells.