Malt1 self-cleavage is critical for regulatory T cell homeostasis and anti-tumor immunity in mice.

Mucosa-associated lymphoid tissue 1 (Malt1) regulates immune cell function by mediating the activation of nuclear factor κB (NF-κB) signaling through both its adaptor and proteolytic function. Malt1 is also a target of its own protease activity and this self-cleavage further contributes to NF-κB activity. Until now, the functional distinction between Malt1 self-cleavage and its general protease function in regulating NF-κB signaling and immune activation remained unclear. Here we demonstrate, using a new mouse model, the importance of Malt1 self-cleavage in regulating expression of NF-κB target genes and subsequent T cell activation. Significantly, we further establish that Treg homeostasis is critically linked to Malt1 function via a Treg intrinsic and extrinsic mechanism. TCR-mediated Malt1 proteolytic activity and self-cleavage was found to drive Il2 expression in conventional CD4+ T cells, thereby regulating Il2 availability for Treg homeostasis. Remarkably, the loss of Malt1-mediated self-cleavage alone was sufficient to cause a significant Treg deficit resulting in increased anti-tumor immune reactivity without associated autoimmunity complications. These results establish for the first time that inhibition of MALT1 proteolytic activity could be a viable therapeutic strategy to augment anti-tumor immunity.

Small chromosomal regions position themselves autonomously according to their chromatin class.

The spatial arrangement of chromatin is linked to the regulation of nuclear processes. One striking aspect of nuclear organization is the spatial segregation of heterochromatic and euchromatic domains. The mechanisms of this chromatin segregation are still poorly understood. In this work, we investigated the link between the primary genomic sequence and chromatin domains. We analyzed the spatial intranuclear arrangement of a human artificial chromosome (HAC) in a xenospecific mouse background in comparison to an orthologous region of native mouse chromosome. The two orthologous regions include segments that can be assigned to three major chromatin classes according to their gene abundance and repeat repertoire: (1) gene-rich and SINE-rich euchromatin; (2) gene-poor and LINE/LTR-rich heterochromatin; and (3) gene-depleted and satellite DNA-containing constitutive heterochromatin. We show, using fluorescence in situ hybridization (FISH) and 4C-seq technologies, that chromatin segments ranging from 0.6 to 3 Mb cluster with segments of the same chromatin class. As a consequence, the chromatin segments acquire corresponding positions in the nucleus irrespective of their chromosomal context, thereby strongly suggesting that this is their autonomous property. Interactions with the nuclear lamina, although largely retained in the HAC, reveal less autonomy. Taken together, our results suggest that building of a functional nucleus is largely a self-organizing process based on mutual recognition of chromosome segments belonging to the major chromatin classes.

Identification of Intellectual Disability Genes in Female Patients with a Skewed X-Inactivation Pattern.

Intellectual disability (ID) is a heterogeneous disorder with an unknown molecular etiology in many cases. Previously, X-linked ID (XLID) studies focused on males because of the hemizygous state of their X chromosome. Carrier females are generally unaffected because of the presence of a second normal allele, or inactivation of the mutant X chromosome in most of their cells (skewing). However, in female ID patients, we hypothesized that the presence of skewing of X-inactivation would be an indicator for an X chromosomal ID cause. We analyzed the X-inactivation patterns of 288 females with ID, and found that 22 (7.6%) had extreme skewing (>90%), which is significantly higher than observed in the general population (3.6%; P = 0.029). Whole-exome sequencing of 19 females with extreme skewing revealed causal variants in six females in the XLID genes DDX3X, NHS, WDR45, MECP2, and SMC1A. Interestingly, variants in genes escaping X-inactivation presumably cause both XLID and skewing of X-inactivation in three of these patients. Moreover, variants likely accounting for skewing only were detected in MED12, HDAC8, and TAF9B. All tested candidate causative variants were de novo events. Hence, extreme skewing is a good indicator for the presence of X-linked variants in female patients.

Nxf7 deficiency impairs social exploration and spatio-cognitive abilities as well as hippocampal synaptic plasticity in mice.

Nuclear RNA export factors (NXF) are conserved in all metazoans and are deemed essential for shuttling RNA across the nuclear envelope and other post-transcriptional processes (such as mRNA metabolism, storage and stability). Disruption of human NXF5 has been implicated in intellectual and psychosocial disabilities. In the present report, we use recently described Nxf7 knockout (KO) mice as an experimental model to analyze in detail the behavioral consequences of clinical NXF5 deficiency. We examined male Nxf7 KO mice using an extended cognitive and behavioral test battery, and recorded extracellular field potentials in the hippocampal CA1 region. We observed various cognitive and behavioral changes including alterations in social exploration, impaired spatial learning and spatio-cognitive abilities. We also defined a new experimental paradigm to discriminate search strategies in Morris water maze and showed significant differences between Nxf7 KO and control animals. Furthermore, while we observed no difference in a nose poke suppression in an conditioned emotional response (CER) protocol, Nxf7 KO mice were impaired in discriminating between differentially reinforced cues in an auditory fear conditioning protocol. This distinct neurocognitive phenotype was accompanied by impaired hippocampal Long-term potentiation (LTP), while long-term depression (LTD) was not affected by Nxf7 deficiency. Our data demonstrate that disruption of murine Nxf7 leads to behavioral phenotypes that may relate to the intellectual and social deficits in patients with NXF5 deficiency.

MALT1 auto-proteolysis is essential for NF-κB-dependent gene transcription in activated lymphocytes.

Mucosa-associated lymphoid tissue 1 (MALT1) controls antigen receptor-mediated signalling to nuclear factor κB (NF-κB) through both its adaptor and protease function. Upon antigen stimulation, MALT1 forms a complex with BCL10 and CARMA1, which is essential for initial IκBα phosphorylation and NF-κB nuclear translocation. Parallel induction of MALT1 protease activity serves to inactivate negative regulators of NF-κB signalling, such as A20 and RELB. Here we demonstrate a key role for auto-proteolytic MALT1 cleavage in B- and T-cell receptor signalling. MALT1 cleavage occurred after Arginine 149, between the N-terminal death domain and the first immunoglobulin-like region, and did not affect its proteolytic activity. Jurkat T cells expressing an un-cleavable MALT1-R149A mutant showed unaltered initial IκBα phosphorylation and normal nuclear accumulation of NF-κB subunits. Nevertheless, MALT1 cleavage was required for optimal activation of NF-κB reporter genes and expression of the NF-κB targets IL-2 and CSF2. Transcriptome analysis confirmed that MALT1 cleavage after R149 was required to induce NF-κB transcriptional activity in Jurkat T cells. Collectively, these data demonstrate that auto-proteolytic MALT1 cleavage controls antigen receptor-induced expression of NF-κB target genes downstream of nuclear NF-κB accumulation.

De novo MECP2 duplications in two females with intellectual disability and unfavorable complete skewed X-inactivation.

Xq28 microduplications of MECP2 are a prominent cause of a severe syndromic form of intellectual disability (ID) in males. Females are usually unaffected through near to complete X-inactivation of the aberrant X chromosome (skewing). In rare cases, affected females have been described due to random X-inactivation. Here, we report on two female patients carrying de novo MECP2 microduplications on their fully active X chromosomes. Both patients present with ID and additional clinical features. Mono-allelic expression confirmed complete skewing of X-inactivation. Consequently, significantly enhanced MECP2 mRNA levels were observed. We hypothesize that the cause for the complete skewing is due to a more harmful mutation on the other X chromosome, thereby forcing the MECP2 duplication to become active. However, we could not unequivocally identify such a second mutation by array-CGH or exome sequencing. Our data underline that, like in males, increased MECP2 dosage in females can contribute to ID too, which should be taken into account in diagnostics.

Evidence for increased SOX3 dosage as a risk factor for X-linked hypopituitarism and neural tube defects.

Genomic duplications of varying lengths at Xq26-q27 involving SOX3 have been described in families with X-linked hypopituitarism. Using array-CGH we detected a 1.1 Mb microduplication at Xq27 in a large family with three males suffering from X-linked hypopituitarism. The duplication was mapped from 138.7 to 139.8 Mb, harboring only two annotated genes, SOX3 and ATP11C, and was shown to be a direct tandem copy number gain. Unexpectedly, the microduplication did not fully segregate with the disease in this family suggesting that SOX3 duplications have variable penetrance for X-linked hypopituitarism. In the same family, a female fetus presenting with a neural tube defect was also shown to carry the SOX3 copy number gain. Since we also demonstrated increased SOX3 mRNA levels in amnion cells derived from an unrelated t(X;22)(q27;q11) female fetus with spina bifida, we propose that increased levels of SOX3 could be a risk factor for neural tube defects.

Ubiquitin ligase HUWE1 regulates axon branching through the Wnt/ß-catenin pathway in a Drosophila model for intellectual disability.

We recently reported that duplication of the E3 ubiquitin ligase HUWE1 results in intellectual disability (ID) in male patients. However, the underlying molecular mechanism remains unknown. We used Drosophila melanogaster as a model to investigate the effect of increased HUWE1 levels on the developing nervous system. Similar to the observed levels in patients we overexpressed the HUWE1 mRNA about 2-fold in the fly. The development of the mushroom body and neuromuscular junctions were not altered, and basal neurotransmission was unaffected. These data are in agreement with normal learning and memory in the courtship conditioning paradigm. However, a disturbed branching phenotype at the axon terminals of the dorsal cluster neurons (DCN) was detected. Interestingly, overexpression of HUWE1 was found to decrease the protein levels of dishevelled (dsh) by 50%. As dsh as well as Fz2 mutant flies showed the same disturbed DCN branching phenotype, and the constitutive active homolog of ß-catenin, armadillo, could partially rescue this phenotype, our data strongly suggest that increased dosage of HUWE1 compromises the Wnt/ß-catenin pathway possibly by enhancing the degradation of dsh.

The mitochondrial solute carrier SLC25A5 at Xq24 is a novel candidate gene for non-syndromic intellectual disability.

Loss-of-function mutations in several different neuronal pathways have been related to intellectual disability (ID). Such mutations often are found on the X chromosome in males since they result in functional null alleles. So far, microdeletions at Xq24 reported in males always have been associated with a syndromic form of ID due to the loss of UBE2A. Here, we report on overlapping microdeletions at Xq24 that do not include UBE2A or affect its expression, in patients with non-syndromic ID plus some additional features from three unrelated families. The smallest region of overlap, confirmed by junction sequencing, harbors two members of the mitochondrial solute carrier family 25, SLC25A5 and SLC25A43. However, identification of an intragenic microdeletion including SLC25A43 but not SLC25A5 in a healthy boy excluded a role for SLC25A43 in cognition. Therefore, our findings point to SLC25A5 as a novel gene for non-syndromic ID. This highly conserved gene is expressed ubiquitously with high levels in cortex and hippocampus, and a presumed role in mitochondrial exchange of ADP/ATP. Our data indicate that SLC25A5 is involved in memory formation or establishment, which could add mitochondrial processes to the wide array of pathways that regulate normal cognitive functions.

Generation and characterization of an Nxf7 knockout mouse to study NXF5 deficiency in a patient with intellectual disability.

Members of the Nuclear eXport Factor (NXF) family are involved in the export of mRNA from the nucleus to the cytoplasm, or hypothesized to play a role in transport of cytoplasmic mRNA. We previously reported on the loss of NXF5 in a male patient with a syndromic form of intellectual disability. To study the functional role of NXF5 we identified the mouse counterpart. Based on synteny, mouse Nxf2 is the ortholog of human NXF5. However, we provide several lines of evidence that mouse Nxf7 is the actual functional equivalent of NXF5. Both Nxf7 and NXF5 are predominantly expressed in the brain, show cytoplasmic localization, and present as granules in neuronal dendrites suggesting a role in cytoplasmic mRNA metabolism in neurons. Nxf7 was primarily detected in the pyramidal cells of the hippocampus and in layer V of the cortex. Similar to human NXF2, mouse Nxf2 is highly expressed in testis and shows a nuclear localization. Interestingly, these findings point to a different evolutionary path for both NXF genes in human and mouse. We thus generated and validated Nxf7 knockout mice, which were fertile and did not present any gross anatomical or morphological abnormalities. Expression profiling in the hippocampus and the cortex did not reveal significant changes between wild-type and Nxf7 knockout mice. However, impaired spatial memory was observed in these KO mice when evaluated in the Morris water maze test. In conclusion, our findings provide strong evidence that mouse Nxf7 is the functional counterpart of human NXF5, which might play a critical role in mRNA metabolism in the brain.

Telomere length homeostasis and telomere position effect on a linear human artificial chromosome are dictated by the genetic background.

Telomere position effect (TPE) is the influence of telomeres on subtelomeric epigenetic marks and gene expression. Previous studies suggested that TPE depends on genetic background. As these analyses were performed on different chromosomes, cell types and species, it remains unclear whether TPE represents a chromosome-rather than genetic background-specific regulation. We describe the development of a Linear Human Artificial Chromosome (L-HAC) as a new tool for telomere studies. The L-HAC was generated through the Cre-loxP-mediated addition of telomere ends to an existing circular HAC (C-HAC). As it can be transferred to genetically distinct cell lines and animal models the L-HAC enables the study of TPE in an unprecedented manner. The HAC was relocated to four telomerase-positive cell lines via microcell-mediated chromosome transfer and subsequently to mice via blastocyst injection of L-HAC(+)-ES-cells. We could show consistent genetic background-dependent adaptation of telomere length and telomere-associated de novo subtelomeric DNA methylation in mouse ES-R1 cells as well as in mice. Expression of the subtelomeric neomycin gene was inversely correlated with telomere length and subtelomeric methylation. We thus provide a new tool for functional telomere studies and provide strong evidence that telomere length, subtelomeric chromatin marks and expression of subtelomeric genes are genetic background dependent.

Copy-number gains of HUWE1 due to replication- and recombination-based rearrangements.

We previously reported on nonrecurrent overlapping duplications at Xp11.22 in individuals with nonsyndromic intellectual disability (ID) harboring HSD17B10, HUWE1, and the microRNAs miR-98 and let-7f-2 in the smallest region of overlap. Here, we describe six additional individuals with nonsyndromic ID and overlapping microduplications that segregate in the families. High-resolution mapping of the 12 copy-number gains reduced the minimal duplicated region to the HUWE1 locus only. Consequently, increased mRNA levels were detected for HUWE1, but not HSD17B10. Marker and SNP analysis, together with identification of two de novo events, suggested a paternally derived intrachromosomal duplication event. In four independent families, we report on a polymorphic 70 kb recurrent copy-number gain, which harbors part of HUWE1 (exon 28 to 3' untranslated region), including miR-98 and let-7f-2. Our findings thus demonstrate that HUWE1 is the only remaining dosage-sensitive gene associated with the ID phenotype. Junction and in silico analysis of breakpoint regions demonstrated simple microhomology-mediated rearrangements suggestive of replication-based duplication events. Intriguingly, in a single family, the duplication was generated through nonallelic homologous recombination (NAHR) with the use of HUWE1-flanking imperfect low-copy repeats, which drive this infrequent NAHR event. The recurrent partial HUWE1 copy-number gain was also generated through NAHR, but here, the homologous sequences used were identified as TcMAR-Tigger DNA elements, a template that has not yet been reported for NAHR. In summary, we showed that an increased dosage of HUWE1 causes nonsyndromic ID and demonstrated that the Xp11.22 region is prone to recombination- and replication-based rearrangements.

JAK2 rearrangements, including the novel SEC31A-JAK2 fusion, are recurrent in classical Hodgkin lymphoma.

The genetics of classical Hodgkin lymphoma (cHL) is poorly understood. The finding of a JAK2-involving t(4;9)(q21;p24) in 1 case of cHL prompted us to characterize this translocation on a molecular level and to determine the prevalence of JAK2 rearrangements in cHL. We showed that the t(4;9)(q21;p24) leads to a novel SEC31A-JAK2 fusion. Screening of 131 cHL cases identified 1 additional case with SEC31A-JAK2 and 2 additional cases with rearrangements involving JAK2. We demonstrated that SEC31A-JAK2 is oncogenic in vitro and acts as a constitutively activated tyrosine kinase that is sensitive to JAK inhibitors. In vivo, SEC31A-JAK2 was found to induce a T-lymphoblastic lymphoma or myeloid phenotype in a murine bone marrow transplantation model. Altogether, we identified SEC31A-JAK2 as a chromosomal aberration characteristic for cHL and provide evidence that JAK2 rearrangements occur in a minority of cHL cases. Given the proven oncogenic potential of this novel fusion, our studies provide new insights into the pathogenesis of cHL and indicate that in at least some cases, constitutive activation of the JAK/STAT pathway is caused by JAK2 rearrangements. The finding that SEC31A-JAK2 responds to JAK inhibitors indicates that patients with cHL and JAK2 rearrangements may benefit from targeted therapies.

Allele-specific copy number analysis of tumors.

We present an allele-specific copy number analysis of the in vivo breast cancer genome. We describe a unique bioinformatics approach, ASCAT (allele-specific copy number analysis of tumors), to accurately dissect the allele-specific copy number of solid tumors, simultaneously estimating and adjusting for both tumor ploidy and nonaberrant cell admixture. This allows calculation of "ASCAT profiles" (genome-wide allele-specific copy-number profiles) from which gains, losses, copy number-neutral events, and loss of heterozygosity (LOH) can accurately be determined. In an early-stage breast carcinoma series, we observe aneuploidy (>2.7n) in 45% of the cases and an average nonaberrant cell admixture of 49%. By aggregation of ASCAT profiles across our series, we obtain genomic frequency distributions of gains and losses, as well as genome-wide views of LOH and copy number-neutral events in breast cancer. In addition, the ASCAT profiles reveal differences in aberrant tumor cell fraction, ploidy, gains, losses, LOH, and copy number-neutral events between the five previously identified molecular breast cancer subtypes. Basal-like breast carcinomas have a significantly higher frequency of LOH compared with other subtypes, and their ASCAT profiles show large-scale loss of genomic material during tumor development, followed by a whole-genome duplication, resulting in near-triploid genomes. Finally, from the ASCAT profiles, we construct a genome-wide map of allelic skewness in breast cancer, indicating loci where one allele is preferentially lost, whereas the other allele is preferentially gained. We hypothesize that these alternative alleles have a different influence on breast carcinoma development.

Coactivated platelet-derived growth factor receptor {alpha} and epidermal growth factor receptor are potential therapeutic targets in intimal sarcoma.

Intimal sarcoma (IS) is a rare, malignant, and aggressive tumor that shows a relentless course with a concomitant low survival rate and for which no effective treatment is available. In this study, 21 cases of large arterial blood vessel IS were analyzed by immunohistochemistry and fluorescence in situ hybridization and selectively by karyotyping, array comparative genomic hybridization, sequencing, phospho-kinase antibody arrays, and Western immunoblotting in search for novel diagnostic markers and potential molecular therapeutic targets. Ex vivo immunoassays were applied to test the sensitivity of IS primary tumor cells to the receptor tyrosine kinase (RTK) inhibitors imatinib and dasatinib. We showed that amplification of platelet-derived growth factor receptor α (PDGFRA) is a common finding in IS, which should be considered as a molecular hallmark of this entity. This amplification is consistently associated with PDGFRA activation. Furthermore, the tumors reveal persistent activation of the epidermal growth factor receptor (EGFR), concurrent to PDGFRA activation. Activated PDGFRA and EGFR frequently coexist with amplification and overexpression of the MDM2 oncogene. Ex vivo immunoassays on primary IS cells from one case showed the potency of dasatinib to inhibit PDGFRA and downstream signaling pathways. Our findings provide a rationale for investigating therapies that target PDGFRA, EGFR, or MDM2 in IS. Given the clonal heterogeneity of this tumor type and the potential cross-talk between the PDGFRA and EGFR signaling pathways, targeting multiple RTKs and aberrant downstream effectors might be required to improve the therapeutic outcome for patients with this disease.

ALK-positive large B-cell lymphomas with cryptic SEC31A-ALK and NPM1-ALK fusions.

We report 2 ALK-positive large B-cell lymphoma cases showing granular cytoplasmic and cytoplasmic/nuclear ALK immunostaining in which cryptic ALK rearrangements were identified by fluorescent in situ hybridization and molecular analysis. In the first case, the ALK-involving t(2;3)(p23;q27) masked the cryptic SEC31A-ALK fusion generated by an insertion of the 5' end of SEC31A (4q21) upstream of the 3' end of ALK. This rearrangement was associated with loss of the 5' end of ALK and duplication of SEC31A-ALK on der(20). In the second case with complex rearrangements of both chromosomes 2, a submicroscopic NPM1-ALK fusion created by insertion of the 3' end of ALK into the NPM1 locus was evidenced. Further studies of SEC31A-ALK showed that this variant fusion transforms IL3-dependent Ba/F3 cells to growth factor independence, and that the ALK inhibitor TAE-684 reduces cell proliferation and kinase activity of SEC31A-ALK and its downstream effectors ERK1/2, AKT, STAT3 and STAT5.

Comparative expressed sequence hybridization studies of t(11;18)(q21;q21)-positive and -negative gastric MALT lymphomas reveal both unique and overlapping gene programs.

Among the genetic abnormalities reported to occur in MALT lymphomas, the translocation t(11;18)(q21;q21) is of particular interest because it is exclusively documented in MALT lymphomas, mainly with gastrointestinal location. It results in the creation of a fusion protein API2-MALT1 that activates the transcription factor NF-kappaB through enhanced IKK gamma polyubiquitination. Here, we apply the recently developed molecular technique termed comparative expressed sequence hybridization to identify differentially expressed chromosomal regions related to the pathogenesis of gastric MALT lymphomas. By comparing t(11;18)(q21;q21)-positive gastric MALT lymphomas to their t(11;18)(q21;q21)-negative counterparts, we found that the location of the MALT1 break point determines a difference in expression pattern within the t(11;18)(q21;q21)-positive group. Moreover, we could define a gastric MALT lymphoma signature, which most likely comprises the regions and genes with significance in the development of MALT lymphomas, by comparing both t(11;18)(q21;q21)-positive and -negative MALT lymphomas to normal lymphoid tissue. Finally, a significant imprint of the marginal zone signature, established by comparing microdissected, splenic B follicles with and without marginal zone, was evident in the expression profile of MALT lymphoma, further supporting a marginal zone origin for this type of B-cell non-Hodgkin's lymphoma.

T-cell/histiocyte-rich large B-cell lymphoma shows transcriptional features suggestive of a tolerogenic host immune response.

BACKGROUND: Gene expression profiling has successfully identified the prognostic significance of the host response in lymphomas. The aggressive T-cell/histiocyte-rich large B-cell lymphoma and the indolent nodular lymphocyte-predominant Hodgkin's lymphoma are both characterized by a paucity of tumor cells embedded in an overwhelming background. The tumor cells of both lymphomas share several characteristics, while the cellular composition of their microenvironment is clearly different. DESIGN AND METHODS: We collected 33 cases of T-cell/histiocyte-rich large B-cell lymphoma and 56 cases of nodular lymphocyte-predominant Hodgkin's lymphoma and performed microarray gene expression profiling on ten cases of each lymphoma, to obtain a better understanding of the lymphoma host response. By quantitative reverse transcriptase polymerase chain reaction we verified that these 20 selected cases were representative of the entire population of T-cell/histiocyte-rich large B-cell and nodular lymphocyte-predominant Hodgkin's lymphomas. RESULTS: We observed that the microenvironment in nodular lymphocyte-predominant Hodgkin's lymphoma is molecularly very similar to a lymph node characterized by follicular hyperplasia, while the microenvironment in T-cell/histiocyte-rich large B-cell lymphoma is clearly different. The T-cell/histiocyte-rich large B-cell lymphoma signature is hallmarked by up-regulation of CCL8, interferon-gamma, indoleamine 2,3 dioxygenase, VSIG4 and Toll-like receptors. These features may be responsible for the recruitment and activation of T cells, macrophages and dendritic cells, characterizing the stromal component of this lymphoma, and may point towards innate immunity and a tumor tolerogenic immune response in T-cell/histiocyte-rich large B-cell lymphoma. CONCLUSIONS: The gene expression profile of T-cell/histiocyte-rich large B-cell lymphoma, in comparison with that of nodular lymphocyte-predominant Hodgkin's lymphoma, shows features suggestive of a distinct tolerogenic host immune response that may play a key role in the aggressive behavior of this lymphoma, and that may serve as a potential target for future therapy.

Dosage-dependent severity of the phenotype in patients with mental retardation due to a recurrent copy-number gain at Xq28 mediated by an unusual recombination.

We report on the identification of a 0.3 Mb inherited recurrent but variable copy-number gain at Xq28 in affected males of four unrelated families with X-linked mental retardation (MR). All aberrations segregate with the disease in the families, and the carrier mothers show nonrandom X chromosome inactivation. Tiling Xq28-region-specific oligo array revealed that all aberrations start at the beginning of the low copy repeat LCR-K1, at position 153.20 Mb, and end just distal to LCR-L2, at 153.54 Mb. The copy-number gain always includes 18 annotated genes, of which RPL10, ATP6AP1 and GDI1 are highly expressed in brain. From these, GDI1 is the most likely candidate gene. Its copy number correlates with the severity of clinical features, because it is duplicated in one family with nonsyndromic moderate MR, is triplicated in males from two families with mild MR and additional features, and is present in five copies in a fourth family with a severe syndromic form of MR. Moreover, expression analysis revealed copy-number-dependent increased mRNA levels in affected patients compared to control individuals. Interestingly, analysis of the breakpoint regions suggests a recombination mechanism that involves two adjacent but different sets of low copy repeats. Taken together, our data strongly suggest that an increased expression of GDI1 results in impaired cognition in a dosage-dependent manner. Moreover, these data also imply that a copy-number gain of an individual gene present in the larger genomic aberration that leads to the severe MECP2 duplication syndrome can of itself result in a clinical phenotype as well.

Computational methods for the detection of cis-regulatory modules.

Metazoan transcription regulation occurs through the concerted action of multiple transcription factors that bind co-operatively to cis-regulatory modules (CRMs). The annotation of these key regulators of transcription is lagging far behind the annotation of the transcriptome itself. Here, we give an overview of existing computational methods to detect these CRMs in metazoan genomes. We subdivide these methods into three classes: CRM scanners screen sequences for CRMs based on predefined models that often consist of multiple position weight matrices (PWMs). CRM builders construct models of similar CRMs controlling a set of co-regulated or co-expressed genes. CRM genome screeners screen sequences or complete genomes for CRMs as homotypic or heterotypic clusters of binding sites for any combination of transcription factors. We believe that CRM scanners are currently the most advanced methods, although their applicability is limited. Finally, we argue that CRM builders that make use of PWM libraries will benefit greatly from future advances and will prove to be most instrumental for the annotation of regulatory regions in metazoan genomes.