Disease-causing 7.4 kb cis-regulatory deletion disrupting conserved non-coding sequences and their interaction with the FOXL2 promotor: implications for mutation screening.

To date, the contribution of disrupted potentially cis-regulatory conserved non-coding sequences (CNCs) to human disease is most likely underestimated, as no systematic screens for putative deleterious variations in CNCs have been conducted. As a model for monogenic disease we studied the involvement of genetic changes of CNCs in the cis-regulatory domain of FOXL2 in blepharophimosis syndrome (BPES). Fifty-seven molecularly unsolved BPES patients underwent high-resolution copy number screening and targeted sequencing of CNCs. Apart from three larger distant deletions, a de novo deletion as small as 7.4 kb was found at 283 kb 5' to FOXL2. The deletion appeared to be triggered by an H-DNA-induced double-stranded break (DSB). In addition, it disrupts a novel long non-coding RNA (ncRNA) PISRT1 and 8 CNCs. The regulatory potential of the deleted CNCs was substantiated by in vitro luciferase assays. Interestingly, Chromosome Conformation Capture (3C) of a 625 kb region surrounding FOXL2 in expressing cellular systems revealed physical interactions of three upstream fragments and the FOXL2 core promoter. Importantly, one of these contains the 7.4 kb deleted fragment. Overall, this study revealed the smallest distant deletion causing monogenic disease and impacts upon the concept of mutation screening in human disease and developmental disorders in particular.

Variation in novel exons (RACEfrags) of the MECP2 gene in Rett syndrome patients and controls.

The study of transcription using genomic tiling arrays has lead to the identification of numerous additional exons. One example is the MECP2 gene on the X chromosome; using 5'RACE and RT-PCR in human tissues and cell lines, we have found more than 70 novel exons (RACEfrags) connecting to at least one annotated exon.. We sequenced all MECP2-connected exons and flanking sequences in 3 groups: 46 patients with the Rett syndrome and without mutations in the currently annotated exons of the MECP2 and CDKL5 genes; 32 patients with the Rett syndrome and identified mutations in the MECP2 gene; 100 control individuals from the same geoethnic group. Approximately 13 kb were sequenced per sample, (2.4 Mb of DNA resequencing). A total of 75 individuals had novel rare variants (mostly private variants) but no statistically significant difference was found among the 3 groups. These results suggest that variants in the newly discovered exons may not contribute to Rett syndrome. Interestingly however, there are about twice more variants in the novel exons than in the flanking sequences (44 vs. 21 for approximately 1.3 Mb sequenced for each class of sequences, p=0.0025). Thus the evolutionary forces that shape these novel exons may be different than those of neighboring sequences.

Primary ciliary dyskinesia associated with normal axoneme ultrastructure is caused by DNAH11 mutations.

Primary ciliary dyskinesia (PCD) is an inherited disorder characterized by perturbed or absent beating of motile cilia, which is referred to as Kartagener syndrome (KS) when associated with situs inversus. We present a German family in which five individuals have PCD and one has KS. PCD was confirmed by analysis of native and cultured respiratory ciliated epithelia with high-speed video microscopy. Respiratory ciliated cells from the affected individuals showed an abnormal nonflexible beating pattern with a reduced cilium bending capacity and a hyperkinetic beat. Interestingly, the axonemal ultrastructure of these respiratory cilia was normal and outer dynein arms were intact, as shown by electron microscopy and immunohistochemistry. Microsatellite analysis indicated genetic linkage to the dynein heavy chain DNAH11 on chromosome 7p21. All affected individuals carried the compound heterozygous DNAH11 mutations c.12384C>G and c.13552_13608del. Both mutations are located in the C-terminal domain and predict a truncated DNAH11 protein (p.Y4128X, p.A4518_A4523delinsQ). The mutations described here were not present in a cohort of 96 PCD patients. In conclusion, our findings support the view that DNAH11 mutations indeed cause PCD and KS, and that the reported DNAH11 nonsense mutations are associated with a normal axonemal ultrastructure and are compatible with normal male fertility.

DNAI1 mutations explain only 2% of primary ciliary dykinesia.

BACKGROUND: Primary ciliary dyskinesia (PCD) is a rare recessive hereditary disorder characterized by dysmotility to immotility of ciliated and flagellated structures. Its main symptoms are respiratory, caused by defective ciliary beating in the epithelium of the upper airways (nose, bronchi and paranasal sinuses). Impairing the drainage of inhaled microorganisms and particles leads to recurrent infections and pulmonary complications. To date, 5 genes encoding 3 dynein protein arm subunits (DNAI1, DNAH5 and DNAH11), the kinase TXNDC3 and the X-linked RPGR have been found to be mutated in PCD. OBJECTIVES: We proposed to determine the impact of the DNAI1 gene on a cohort of unrelated PCD patients (n = 104) recruited without any phenotypic preselection. METHODS: We used denaturing high-performance liquid chromatography and sequencing to screen for mutations in the coding and splicing site sequences of the gene DNAI1. RESULTS: Three mutations were identified: a novel missense variant (p.Glu174Lys) was found in 1 patient and 2 previously reported variants were identified (p.Trp568Ser in 1 patient and IVS1+2_3insT in 3 patients). Overall, mutations on both alleles of gene DNAI1 were identified in only 2% of our clinically heterogeneous cohort of patients. CONCLUSION: We conclude that DNAI1 gene mutation is not a common cause of PCD, and that major or several additional disease gene(s) still remain to be identified before a sensitive molecular diagnostic test can be developed for PCD.

Mutations in the TMPRSS3 gene are a rare cause of childhood nonsyndromic deafness in Caucasian patients.

Two loci for nonsyndromic recessive deafness located on chromosome 21q22.3 have previously been reported, DFNB8 and DFNB10. Recently a gene which encodes a transmembrane serine protease, TMPRSS3 or ECHOS1, was found to be responsible for both the DFNB8 and DFNB10 phenotypes. To determine the contribution of TMPRSS3 mutations in the general congenital/childhood nonsyndromic deaf population we performed mutation analysis of the TMPRSS3 gene in 448 unrelated deaf patients from Spain, Italy, Greece, and Australia who did not have the common 35delG GJB2 mutation. From the 896 chromosomes studied we identified two novel pathogenic mutations accounting for four mutant alleles and at least 16 nonpathogenic sequence variants. The pathogenic mutations were a 1-bp deletion resulting in a frameshift and an amino acid substitution in the LDLRA domain of TMPRSS3. From this and another study we estimate the frequency of TMPRSS3 mutations in our sample as 0.45%, and approximately 0.38% in the general Caucasian childhood deaf population. However, TMPRSS3 is still an important contributor to genetic deafness in populations with large consanguineous families.

Molecular and clinical characteristics in 46 families affected with Peutz-Jeghers syndrome.

Germline mutations of the tumor suppressor gene LKB1/STK11 are responsible for the Peutz-Jeghers syndrome (PJS), an autosomal-dominant disorder characterized by mucocutaneous pigmentation, hamartomatous polyps, and an increased risk of associated malignancies. In this study, we assessed the presence of pathogenic mutations in the LKB1/STK11 gene in 46 unrelated PJS families, and also carried genotype-phenotype correlation in regard of the development of cancer in 170 PJS patients belonging to these families. All LKB1/STK11 variants detected with single-strand conformational polymorphism were confirmed by direct sequencing, and those without LKB1/STK11 mutation were further submitted to Southern blot analysis for detection of deletions/rearrangements. Statistical analysis for genotype-phenotype correlation was performed. In 59% (27/46) of unrelated PJS cases, pathogenic mutations in the LKB1/STK11 gene, including 9 novel mutations, were identified. The new mutations were 2 splice site deletion-insertions, 2 missenses, 1 nonsense, and 4 abnormal splice sites. Genotype-phenotype analysis did not yield any significant differences between patients carrying mutations in LKB1/STK11 versus those without mutations, even with respect to primary biliary adenocarcinoma. This study presents the molecular characterization and cancer occurrence of a large cohort of PJS patients, increases the mutational spectrum of LKB1/STK11 allelic variants worldwide, and provides a new insight useful for clinical diagnosis and genetic counseling of PJS families.

GENCODE: producing a reference annotation for ENCODE.

BACKGROUND: The GENCODE consortium was formed to identify and map all protein-coding genes within the ENCODE regions. This was achieved by a combination of initial manual annotation by the HAVANA team, experimental validation by the GENCODE consortium and a refinement of the annotation based on these experimental results. RESULTS: The GENCODE gene features are divided into eight different categories of which only the first two (known and novel coding sequence) are confidently predicted to be protein-coding genes. 5' rapid amplification of cDNA ends (RACE) and RT-PCR were used to experimentally verify the initial annotation. Of the 420 coding loci tested, 229 RACE products have been sequenced. They supported 5' extensions of 30 loci and new splice variants in 50 loci. In addition, 46 loci without evidence for a coding sequence were validated, consisting of 31 novel and 15 putative transcripts. We assessed the comprehensiveness of the GENCODE annotation by attempting to validate all the predicted exon boundaries outside the GENCODE annotation. Out of 1,215 tested in a subset of the ENCODE regions, 14 novel exon pairs were validated, only two of them in intergenic regions. CONCLUSION: In total, 487 loci, of which 434 are coding, have been annotated as part of the GENCODE reference set available from the UCSC browser. Comparison of GENCODE annotation with RefSeq and ENSEMBL show only 40% of GENCODE exons are contained within the two sets, which is a reflection of the high number of alternative splice forms with unique exons annotated. Over 50% of coding loci have been experimentally verified by 5' RACE for EGASP and the GENCODE collaboration is continuing to refine its annotation of 1% human genome with the aid of experimental validation.

Evolutionary discrimination of mammalian conserved non-genic sequences (CNGs).

Analysis of the human and mouse genomes identified an abundance of conserved non-genic sequences (CNGs). The significance and evolutionary depth of their conservation remain unanswered. We have quantified levels and patterns of conservation of 191 CNGs of human chromosome 21 in 14 mammalian species. We found that CNGs are significantly more conserved than protein-coding genes and noncoding RNAS (ncRNAs) within the mammalian class from primates to monotremes to marsupials. The pattern of substitutions in CNGs differed from that seen in protein-coding and ncRNA genes and resembled that of protein-binding regions. About 0.3% to 1% of the human genome corresponds to a previously unknown class of extremely constrained CNGs shared among mammals.

Nineteen additional unpredicted transcripts from human chromosome 21.

The identification of all human chromosome 21 (HC21) genes is a necessary step in understanding the molecular pathogenesis of trisomy 21 (Down syndrome). The first analysis of the sequence of 21q included 127 previously characterized genes and predicted an additional 98 novel anonymous genes. Recently we evaluated the quality of this annotation by characterizing a set of HC21 open reading frames (C21orfs) identified by mapping spliced expressed sequence tags (ESTs) and predicted genes (PREDs), identified only in silico. This study underscored the limitations of in silico-only gene prediction, as many PREDs were incorrectly predicted. To refine the HC21 annotation, we have developed a reliable algorithm to extract and stringently map sequences that contain bona fide 3' transcript ends to the genome. We then created a specific 21q graphical display allowing an integrated view of the data that incorporates new ESTs as well as features such as CpG islands, repeats, and gene predictions. Using these tools we identified 27 new putative genes. To validate these, we sequenced previously cloned cDNAs and carried out RT-PCR, 5'- and 3'-RACE procedures, and comparative mapping. These approaches substantiated 19 new transcripts, thus increasing the HC21 gene count by 9.5%. These transcripts were likely not previously identified because they are small and encode small proteins. We also identified four transcriptional units that are spliced but contain no obvious open reading frame. The HC21 data presented here further emphasize that current gene prediction algorithms miss a substantial number of transcripts that nevertheless can be identified using a combination of experimental approaches and multiple refined algorithms.

The transmembrane serine protease (TMPRSS3) mutated in deafness DFNB8/10 activates the epithelial sodium channel (ENaC) in vitro.

TMPRSS3 encodes a transmembrane serine protease that contains both LDLRA and SRCR domains and is mutated in non-syndromic autosomal recessive deafness (DFNB8/10). To study its function, we cloned the mouse ortholog which maps to Mmu17, which is structurally similar to the human gene and encodes a polypeptide with 88% identity to the human protein. RT-PCR and RNA in situ hybridization on rat and mouse cochlea revealed that Tmprss3 is expressed in the spiral ganglion, the cells supporting the organ of Corti and the stria vascularis. RT-PCR on mouse tissues showed expression in the thymus, stomach, testis and E19 embryos. Transient expression of wild-type or tagged TMPRSS3 protein showed a primary localization in the endoplasmic reticulum. The epithelial amiloride-sensitive sodium channel (ENaC), which is expressed in many sodium-reabsorbing tissues including the inner ear and is regulated by membrane-bound channel activating serine proteases (CAPs), is a potential substrate of TMPRSS3. In the Xenopus oocyte expression system, proteolytic processing of TMPRSS3 was associated with increased ENaC mediated currents. In contrast, 6 TMPRSS3 mutants (D103G, R109W, C194F, W251C, P404L, C407R) causing deafness and a mutant in the catalytic triad of TMPRSS3 (S401A), failed to undergo proteolytic cleavage and activate ENaC. These data indicate that important signaling pathways in the inner ear are controlled by proteolytic cleavage and suggest: (i) the existence of an auto-catalytic processing by which TMPRSS3 would become active, and (ii) that ENaC could be a substrate of TMPRSS3 in the inner ear.

Mutations in the DNAH11 (axonemal heavy chain dynein type 11) gene cause one form of situs inversus totalis and most likely primary ciliary dyskinesia.

Primary ciliary dyskinesia (PCD; MIM 242650) is an autosomal recessive disorder of ciliary dysfunction with extensive genetic heterogeneity. PCD is characterized by bronchiectasis and upper respiratory tract infections, and half of the patients with PCD have situs inversus (Kartagener syndrome). We characterized the transcript and the genomic organization of the axonemal heavy chain dynein type 11 (DNAH11) gene, the human homologue of murine Dnah11 or lrd, which is mutated in the iv/iv mouse model with situs inversus. To assess the role of DNAH11, which maps on chromosome 7p21, we searched for mutations in the 82 exons of this gene in a patient with situs inversus totalis, and probable Kartagener syndrome associated with paternal uniparental disomy of chromosome 7 (patUPD7). We identified a homozygous nonsense mutation (R2852X) in the DNAH11 gene. This patient is remarkable because he is also homozygous for the F508del allele of the cystic fibrosis transmembrane conductance regulator (CFTR) gene. Sequence analysis of the DNAH11 gene in an additional 6 selected PCD sibships that shared DNAH11 alleles revealed polymorphic variants and an R3004Q substitution in a conserved position that might be pathogenic. We conclude that mutations in the coding region of DNAH11 account for situs inversus totalis and probably a minority of cases of PCD.

In vitro analyses of known and novel RUNX1/AML1 mutations in dominant familial platelet disorder with predisposition to acute myelogenous leukemia: implications for mechanisms of pathogenesis.

Familial platelet disorder with predisposition to acute myelogenous leukemia (FPD/AML) is an autosomal dominant familial platelet disorder characterized by thrombocytopenia and a propensity to develop AML. Mutation analyses of RUNX1 in 3 families with FPD/AML showing linkage to chromosome 21q22.1 revealed 3 novel heterozygous point mutations (K83E, R135fsX177 (IVS4 + 3delA), and Y260X). Functional investigations of the 7 FPD/AML RUNX1 Runt domain point mutations described to date (2 frameshift, 2 nonsense, and 3 missense mutations) were performed. Consistent with the position of the mutations in the Runt domain at the RUNX1-DNA interface, DNA binding of all mutant RUNX1 proteins was absent or significantly decreased. In general, missense and nonsense RUNX1 proteins retained the ability to heterodimerize with PEBP2beta/CBFbeta and inhibited transactivation of a reporter gene by wild-type RUNX1. Colocalization of mutant RUNX1 and PEBP2beta/CBFbeta in the cytoplasm was observed. These results suggest that the sequestration of PEBP2beta/CBFbeta by mutant RUNX1 may cause the inhibitory effects. While haploinsufficiency of RUNX1 causes FPD/AML in some families (deletions and frameshifts), mutant RUNX1 proteins (missense and nonsense) may also inhibit wild-type RUNX1, possibly creating a higher propensity to develop leukemia. This is consistent with the hypothesis that a second mutation has to occur, either in RUNX1 or another gene, to cause leukemia among individuals harboring RUNX1 FPD/AML mutations and that the propensity to acquire these additional mutations is determined, at least partially, by the initial RUNX1 mutation.