A balance between activating and repressive histone modifications regulates cystic fibrosis transmembrane conductance regulator (CFTR) expression in vivo

The genetic mechanisms that regulate CFTR, the gene responsible for cystic fibrosis, have been widely investigated in cultured cells. However, mechanisms responsible for tissue-specific and time-specific expression are not completely elucidated in vivo. Through the survey of public databases, we found that the promoter of CFTR was associated with bivalent chromatin in human embryonic stem (ES) cells. In this work, we analyzed fetal (at different stages of pregnancy) and adult tissues and showed that, in digestive and lung tissues, which expressed CFTR, H3K4me3 was maintained in the promoter. Histone acetylation was high in the promoter and in two intronic enhancers, especially in fetal tissues. In contrast, in blood cells, which did not express CFTR, the bivalent chromatin was resolved (the promoter was labeled by the silencing mark H3K27me3). Cis-regulatory sequences were associated with lowly acetylated histones. We also provide evidence that the tissue-specific expression of CFTR is not regulated by dynamic changes of DNA methylation in the promoter. Overall, this work shows that a balance between activating and repressive histone modifications in the promoter and intronic enhancers results in the fine regulation of CFTR expression during development, thereby ensuring tissue specificity.     

Amplification of CFTR Exon 9 Sequences to Multiple Locations in the Human Genome

Cloning and characterization of the cystic fibrosis transmembrane conductance regulator (CFTR) gene led to the identification and isolation of cDNA and genomic sequences that cross-hybridized to the first nucleotide binding fold of CFTR. DNA sequence analysis of these clones showed that the cross-hybridizing sequences corresponded to CFTR exon 9 and its flanking introns, juxtapositioned with two segments of LINE1 sequences. The CFTR sequence appeared to have been transcribed from the opposite direction of the gene, reversely transcribed, and co-integrated with the L1 sequences into a chromosome location distinct from that of the CFTR locus. Based on hybridization intensity and complexity of the restriction fragments, it was estimated that there were at least 10 copies of the “amplified” CFTR exon 9 sequences in the human genome. Furthermore, when DNA segments adjacent to the insertion site were used in genomic DNA blot hybridization analysis, multiple copies were also detected. The overall similarity between these CFTR exon 9-related sequences suggested that they were derived from a single retrotransposition event and subsequent sequence amplification. The amplification unit appeared to be greater than 30 kb. Physical mapping studies includingin situhybridization to human metaphase chromosomes showed that multiple copies of these amplified sequences (with and without the CFTR exon 9 insertion) were dispersed throughout the genome. These findings provide insight into the structure and evolution of the human genome.     

A balance between activating and repressive histone modifications regulates cystic fibrosis transmembrane conductance regulator (CFTR) expression in vivo

The genetic mechanisms that regulate CFTR, the gene responsible for cystic fibrosis, have been widely investigated in cultured cells. However, mechanisms responsible for tissue-specific and time-specific expression are not completely elucidated in vivo. Through the survey of public databases, we found that the promoter of CFTR was associated with bivalent chromatin in human embryonic stem (ES) cells. In this work, we analyzed fetal (at different stages of pregnancy) and adult tissues and showed that, in digestive and lung tissues, which expressed CFTR, H3K4me3 was maintained in the promoter. Histone acetylation was high in the promoter and in two intronic enhancers, especially in fetal tissues. In contrast, in blood cells, which did not express CFTR, the bivalent chromatin was resolved (the promoter was labeled by the silencing mark H3K27me3). Cis-regulatory sequences were associated with lowly acetylated histones. We also provide evidence that the tissue-specific expression of CFTR is not regulated by dynamic changes of DNA methylation in the promoter. Overall, this work shows that a balance between activating and repressive histone modifications in the promoter and intronic enhancers results in the fine regulation of CFTR expression during development, thereby ensuring tissue specificity.     

Intronic alternative splicing regulators identified by comparative genomics in nematodes

Many alternative splicing events are regulated by pentameric and hexameric intronic sequences that serve as binding sites for splicing regulatory factors. We hypothesized that intronic elements that regulate alternative splicing are under selective pressure for evolutionary conservation. Using a Wobble Aware Bulk Aligner genomic alignment of Caenorhabditis elegans and Caenorhabditis briggsae, we identified 147 alternatively spliced cassette exons that exhibit short regions of high nucleotide conservation in the introns flanking the alternative exon. In vivo experiments on the alternatively spliced let-2 gene confirm that these conserved regions can be important for alternative splicing regulation. Conserved intronic element sequences were collected into a dataset and the occurrence of each pentamer and hexamer motif was counted. We compared the frequency of pentamers and hexamers in the conserved intronic elements to a dataset of all C. elegans intron sequences in order to identify short intronic motifs that are more likely to be associated with alternative splicing. High-scoring motifs were examined for upstream or downstream preferences in introns surrounding alternative exons. Many of the high- scoring nematode pentamer and hexamer motifs correspond to known mammalian splicing regulatory sequences, such as (T)GCATG, indicating that the mechanism of alternative splicing regulation is well conserved in metazoans. A comparison of the analysis of the conserved intronic elements, and analysis of the entire introns flanking these same exons, reveals that focusing on intronic conservation can increase the sensitivity of detecting putative splicing regulatory motifs. This approach also identified novel sequences whose role in splicing is under investigation and has allowed us to take a step forward in defining a catalog of splicing regulatory elements for an organism. In vivo experiments confirm that one novel high-scoring sequence from our analysis, (T)CTATC, is important for alternative splicing regulation of the unc-52 gene.     

RNA structure is a key regulatory element in pathological ATM and CFTR pseudoexon inclusion events

Genomic variations deep in the intronic regions of pre-mRNA molecules are increasingly reported to affect splicing events. However, there is no general explanation why apparently similar variations may have either no effect on splicing or cause significant splicing alterations. In this work we have examined the structural architecture of pseudoexons previously described in ATM and CFTR patients. The ATM case derives from the deletion of a repressor element and is characterized by an aberrant 5′ss selection despite the presence of better alternatives. The CFTR pseudoexon instead derives from the creation of a new 5′ss that is used while a nearby pre-existing donor-like sequence is never selected. Our results indicate that RNA structure is a major splicing regulatory factor in both cases. Furthermore, manipulation of the original RNA structures can lead to pseudoexon inclusion following the exposure of unused 5′ss already present in their wild-type intronic sequences and prevented to be recognized because of their location in RNA stem structures. Our data show that intrinsic structural features of introns must be taken into account to understand the mechanism of pseudoexon activation in genetic diseases. Our observations may help to improve diagnostics prediction programmes and eventual therapeutic targeting.     

Comparative and genetic analysis of the porcine glucocerebrosidase (GBA) gene

The genomic sequence of the porcine (Sus scrofa) glucocerebrosidase (GBA) gene (5.7 kb), encoding glucocerebrosidase (glucosylceramidase; acid beta-glucosidase; EC 3.2.1.45), was determined and compared with human (Homo sapiens) GBA and GBAP (pseudogene). The porcine gene harbours 11 exons and 10 introns, and the genomic organization is identical with human GBA. The exon sequences, coding for signal peptide and mature protein, show 81% and 90% sequence identity, respectively, with the corresponding human GBA sequences. Short interspersed elements, SINEs (PREs), are present in introns 2, 4 and 7. There is no evidence of a pseudogene in pig. The deduced protein sequence of GBA consists of 39 amino acids of signal peptide (long form) and 497 amino acids of the mature protein; the latter shows 90% sequence identity with the human protein. Four polymorphisms were observed within the porcine gene: insertion/deletion of one of the two SINEs (PREs) in intron 2 (locus PREA); deletion of a 37- to 39-bp stretch in intron 4 (one direct repeat and 5′ end of PRE); deletion of a 47-bp stretch in the middle part of PRE in intron 4 (locus PREB); and single-base transition (C–T) in intron 6 (locus HaeIII–RFLP). GBA was assigned to chromosome 4q21 by FISH and was localized to the same region by linkage analysis and RH mapping, i.e., to the chromosome 4 segment where quantitative trait loci for growth and some carcass traits are located.     

A complete sequence of the rice sucrose synthase-1 (RSs1) gene

Using a fragment of the maize sucrose synthase gene Sh-1 as probe, the rice genome was shown to contain at least three genes encoding sucrose synthase. One of these genes was isolated from a genomic library, and its full sequence, including 1.7 kb of 5 flanking sequence and 0.9 kb of 3 flanking sequence, is reported. The new rice gene, designated RSs1, is highly homologous to maize Sh-1 (approx. 94% identity in derived amino acid sequence), and contains an identical intron-exon structure (16 exons and 15 introns). Both RSs1 and maize Sh-1 show similar sequence homologies to a second rice sucrose synthase gene described recently (designated RSs2, Yu et al. (1992) Plant Mol Biol 18: 139–142), although both the rice genes predict an extra 6 amino acids at the C-terminus of the protein when compared to the maize gene. The RSs1 5 flanking sequence contains a number of promoter-like sequences, including putative protein-binding regions similar to maize zein genes.     

Germline repertoire of the immunoglobulin V H 3 family in rhesus monkeys

 To facilitate molecular studies of antibody responses in rhesus monkeys (Macaca mulatta), we cloned and sequenced germline segments from its largest and most diverse immunoglobulin heavy-chain gene family, V _H 3. Using a PCR-based approach, we characterized 29 sequences, 20 with open reading frames (ORFs) and 9 pseudogenes. The leader sequences, introns, exons, and recombination signal sequences of M. mulatta V _H 3 gene segments are not strictly identical to those of humans, but the mature coding regions demonstrate, on average, greater than 90% sequence similarity. Although the framework regions are more highly conserved, the complementarity-determining regions (CDRs) also show strong similarities, and their predicted three-dimensional structures resemble those of their human homologues. In one instance, homologous macaque and human CDR1 sequences were 100% identical at the nucleotide level, and some CDR2s shared nucleotide identity as high as 96.5%. However, some rhesus V _H 3 ORFs have unusual structural features, including atypical CDR lengths and uncommon amino acids at structurally crucial positions. The similarity of rhesus and human V _H 3 homologues reinforces the notion that humoral immunity in this nonhuman primate species is an appropriate system for modeling human antibody responses. Received: 10 August 1999 / Revised: 30 December 1999     

Nucleotide sequence and genomic organization of a human T lymphocyte serine protease gene

We have determined the nucleotide sequence of 4508 base pairs of human genomic DNA which contain the human serine esterase gene from cytotoxic T lymphocytes (SECT) (equivalent to the 1-3E cDNA clone) and include 879 bp of 5' flanking DNA and 393 bp of 3' flanking DNA. The gene consists of five exons of 88, 148, 136, 261, and 257 nucleotides separated by four introns of 1043, 455, 205, and 643 nucleotides. The location of introns with respect to protein coding sequences in the SECT gene is identical to that of the human cathepsin G and murine granzyme B genes. Comparison of SECT gene exonic sequences to murine granzyme B-F cDNA sequences indicates similarities of 75 and 72% for granzymes B and C and 61, 59, and 61% for granzymes D, E, and F, respectively. The 5' flanking sequence of the SECT gene showed similarity only to the 5' flanking sequence of the murine granzyme B gene, indicating that these genes are homologous. Comparison of the SECT gene sequence to the human cathepsin G sequence indicated no similarity in the 5' flanking DNA although the exonic sequences show 64% sequence similarity overall and 45% sequence similarity in the respective 3' untranslated regions. These similarities suggest that the SECT and cathepsin G genes are members of the same family of serine protease genes. Evidence from high and low stringency Southern transfer analysis of human genomic DNA indicates the presence of another gene of at least 85% sequence similarity to the SECT gene.     

Discovery of intron polymorphisms in cultivated tomato using both tomato and Arabidopsis genomic information

A low level of genetic variation has limited the application of molecular markers for characterizing important traits in cultivated tomato. To detect polymorphisms in tomato conserved ortholog sets (COS), expressed sequence tags (ESTs) were searched against tomato and Arabidopsis genomic sequences to define the positions of introns. Introns were amplified from 12 different accessions of tomato by polymerase chain reaction and nucleotide sequences were determined by sequencing. Results indicated that there was a possibility of 71% to amplify introns from tomato genomic DNA through this approach. A total of 201 introns were sequenced from 86 COS unigenes. The intron positions and numbers were conserved between tomato and Arabidopsis, but average intron length was three times longer in tomato than in Arabidopsis. A total of 307 single nucleotide polymorphisms (SNPs) and 75 indels were detected in introns of 57 COS unigenes among 12 tomato lines. Within cultivated tomato germplasm 172 SNPs and 47 indels were detected in introns of 33 COS unigenes. In addition, 41 SNPs were identified in the exons of 27 COS unigenes. The frequency of SNPs was 2.4 times higher in introns than in exons in the 22 COS unigenes having both intronic and exonic polymorphisms. These results indicate that intronic regions may contain sufficient variation to develop sufficient marker resources for genome-wide analysis in cultivated tomato.     

Nucleotide sequence of the gene for human factor IX (antihemophilic factor B)

Two different human genomic DNA libraries were screened for the gene for blood coagulation factor IX by employing a cDNA for the human protein as a hybridization probe. Five overlapping lambda phages were identified that contained the gene for factor IX. The complete DNA sequence of about 38 kilobases for the gene and the adjacent 5' and 3' flanking regions was established by the dideoxy chain termination and chemical degradation methods. The gene contained about 33.5 kilobases of DNA, including seven introns and eight exons within the coding and 3' noncoding regions of the gene. The eight exons code for a prepro leader sequence and 415 amino acids that make up the mature protein circulating in plasma. The intervening sequences range in size from 188 to 9473 nucleotides and contain four Alu repetitive sequences, including one in intron A and three in intron F. A fifth Alu repetitive sequence was found immediately flanking the 3' end of the gene. A 50 base pair insert in intron A was found in a clone from one of the genomic libraries but was absent in clones from the other library. Intron A as well as the 3' noncoding region of the gene also contained alternating purine-pyrimidine sequences that provide potential left-handed helical DNA or Z-DNA structures for the gene. KpnI repetitive sequences were identified in intron D and the region flanking the 5' end of the gene. The 5' flanking region also contained a 1.9-kb HindIII subfamily repeat. The seven introns in the gene for factor IX were located in essentially the same position as the seven introns in the gene for human protein C, while the first three were found in positions identical with those in the gene for human prothrombin.     

Carbonic Anhydrase-Related Protein XI: Structure of the Gene in the Greater False Vampire Bat (Megaderma lyra) Compared with Human and Domestic Pig

Carbonic anhydrase-related protein XI (CA-RP XI) is a member of the α-carbonic anhydrase family (encoded by the gene CA-11), which has lost features of the active site required for enzymatic activity. Using PCR, we amplified CA-11 from genomic DNA of the bat Megaderma lyra. To elucidate the gene structure, we sequenced PCR products and compared their sequences with genomic and mRNA sequences known from human and domestic pig. We identified and sequenced eight introns in the bat CA-11. Five introns (introns 3–7) are located in identical or similar positions in other members of the vertebrate α-carbonic anhydrase gene family. Two 5′ introns and one 3′ intron are located in the regions of little or no sequence similarity with other members of the gene family. The low sequence similarity and additional introns suggest a separate evolutionary origin for the 5′ and 3′ portions of the CA-RP XI gene.     

Structure and Expression of the Human Lysyl Hydroxylase Gene (PLOD): Introns 9 and 16 Contain Alu Sequences at the Sites of Recombination in Ehlers-Danlos Syndrome Type VI Patients

Lysyl hydroxylase (EC 1.14.11.4) catalyzes the formation of hydroxylysine in collagens by the hydroxylation of lysine residues in peptide linkages. This enzyme activity is known to be reduced in patients with the type VI variant of the Ehlers-Danlos syndrome, and the first mutations in the lysyl hydroxylase gene (PLOD) have recently been identified. We have now isolated genomic clones for human lysyl hydroxylase and determined the complete structure of the gene, which contains 19 exons and a 5′ flanking region with characteristics shared by housekeeping genes. The constitutive expression of the gene in different tissues further suggests that lysyl hydroxylase has an essential function. We have sequenced the introns of the gene in the region where many mutations and rearrangements analyzed to date are concentrated. Intron 9 and intron 16 show extensive homology resulting from the many Alu sequences found in these introns. Intron 9 contains five and intron 16 eight Alu sequences. The high homology and many short identical or complementary sequences in these introns generate many potential recombination sites with the gene. The delineation of the structure of the lysyl hydroxylase gene contributes significantly to our understanding of the rearrangements in the genome of Ehlers-Danlos type VI patients.     

A comparative genomic analysis of the cow,pig, and human CFTR genes identifies potential intronic regulatory elements

The identification of sequences within noncoding regions of genes that are conserved between several species may indicate potential regulatory elements. This is important for genes with complex control mechanisms such as the cystic fibrosis transmembrane conductance regulator (CFTR). CFTR demonstrates similar patterns of temporal and spatial expression in human and sheep, but these differ significantly in mouse cftr. The complete sheep CFTR sequence is unavailable so we annotated BAC clones encompassing the CFTR gene from two other artiodactyl species (cow and pig) for comparative sequence analysis. Regions of introns 2, 3, 10, 17a, 18, and 21 and 3' flanking sequence corresponding to human CFTR DNase I hypersensitive sites (DHS) showed high homology in the cow and pig. Cross-species sequence conservation also enabled finer mapping of other human DHS, including those in introns 1, 16, and 20. Additional potential regulatory elements not associated with human DHS were also identified.     

Intron requirement for expression of the human purine nucleoside phosphorylase gene.

Abbreviated purine nucleoside phosphorylase (PNP) genes were engineered to determine the effect of introns on human PNP gene expression. PNP minigenes containing the first intron (complete or shortened from 2.9 kb down to 855 bp), the first two introns or all five PNP introns resulted in substantial human PNP isozyme expression after transient transfection of murine NIH 3T3 cells. Low level human PNP activity was observed after transfection with a PNP minigene containing the last three introns. An intronless PNP minigene construct containing the PNP cDNA fused to genomic flanking sequences resulted in undetectable human PNP activity. Heterogeneous, stable NIH 3T3 transfectants of intron-containing PNP minigenes (verified by Southern analysis), expressed high levels of PNP activity and contained appropriately processed 1.7 kb message visualized by northern analysis. Stable transfectants of the intronless PNP minigene (40-45 copies per haploid genome) contained no detectable human PNP isozyme or mRNA. Insertion of the 855 bp shortened intron 1 sequence in either orientation upstream or downstream of a chimeric PNP promoter-bacterial chloramphenicol acetyltransferase (CAT) gene resulted in a several-fold increase in CAT expression in comparison with the parental PNP-CAT construct. We conclude that human PNP gene expression at the mRNA and protein level is dependent on the presence of intronic sequences and that the level of PNP expression varies directly with the number of introns included. The disproportionately greatest effect of intron 1 can be explained by the presence of an enhancer-like element retained in the shortened 855 bp intron 1 sequence.     

Complete mutational screening of the CFTR gene in 120 patients with pulmonary disease

In order to determine the possible role of the cystic fibrosis transmembrane regulator (CFTR) gene in pulmonary diseases not due to cystic fibrosis, a complete screening of the CFTR gene was performed in 120 Italian patients with disseminated bronchiectasis of unknown cause (DBE), chronic bronchitis (CB), pulmonary emphysema (E), lung cancer (LC), sarcoidosis (S) and other forms of pulmonary disease. The 27 exons of the CFTR gene and their intronic flanking regions were analyzed by denaturing gradient gel electrophoresis and automatic sequencing. Mutations were detected in 11/23 DBE (P=0.009), 7/25 E, 5/27 CB, 5/26 LC, 5/8 S (P=0.013), 1/4 tuberculosis, and 1/5 pneumonia patients, and in 5/33 controls. Moreover, the IVS8–5T allele was detected in 6/25 E patients (P=0.038). Four new mutations were identified: D651N, 2377C/T, E826K, and P1072L. These results confirm the involvement of the CFTR gene in disseminated bronchiectasis of unknown origin, and suggest a possible role for CFTR gene mutations in sarcoidosis, and for the 5T allele in pulmonary emphysema. Received: 9 September 1998 / Accepted: 21 October 1998     

A CELLULOSE SYNTHASE (CESA) GENE FROM THE RED ALGA PORPHYRA YEZOENSIS (RHODOPHYTA)1

The cell walls of Porphyra species, like those of land plants, contain cellulose microfibrils that are synthesized by clusters of cellulose synthase enzymes (“terminal complexes”), which move in the plasma membrane. However, the morphologies of the Porphyra terminal complexes and the cellulose microfibrils they produce differ from those of land plants. To characterize the genetic basis for these differences, we have identified, cloned, and sequenced a cellulose synthase (CESA) gene from Porphyra yezoensis Ueda strain TU‐1. A partial cDNA sequence was identified in the P. yezoensis expressed sequence tag (EST) index using a land plant CESA sequence as a query. High‐efficiency thermal asymmetric interlaced PCR was used to amplify sequences upstream of the cDNA sequence from P. yezoensis genomic DNA. Using the resulting genomic sequences as queries, we identified additional EST sequences and a full‐length cDNA clone, which we named PyCESA1. The conceptual translation of PyCESA1 includes the four catalytic domains and the N‐ and C‐terminal transmembrane domains that characterize CESA proteins. Genomic PCR demonstrated that PyCESA1 contains no introns. Southern blot analysis indicated that P. yezoensis has at least three genomic sequences with high similarity to the cloned gene; two of these are pseudogenes based on analysis of amplified genomic sequences. The P. yezoensis CESA peptide sequence is most similar to cellulose synthase sequences from the oomycete Phytophthora infestans and from cyanobacteria. Comparing the CESA genes of P. yezoensis and land plants may facilitate identification of sequences that control terminal complex and cellulose microfibril morphology.     

SLC26A9 stimulates CFTR expression and function in human bronchial cell lines

We investigated the possible functional‐ and physical protein‐interactions between two airway Cl^? channels, SLC26A9 and CFTR. Bronchial CFBE41o‐ cell lines expressing CFTR_WT or CFTR_ΔF508 were transduced with SLC26A9. Immunoblots identified a migrating band corresponding to SLC26A9 present in whole‐cell lysates as on apical membrane of cells grown on polarized filters. CFTR levels were increased by the presence of SLC26A9 in both CFTR_WT and CFTR_ΔF508 cell lines. In CFBE41o‐ cells and CFBE41o‐/CFTR_WT cells transduced with SLC26A9, currents associated to the protein expression were not detected. However, the forskolin (FK)‐stimulated currents were enhanced in SLC26A9‐transduced cells compared to control cells. Therefore, the presence of SLC26A9 resulted in an increase in CFTR activity (same % of CFTR_(inh)‐172 or GlyH‐101 inhibition in both groups). In CFBE41o‐/CFTR_ΔF508 cells transduced with SLC26A9 (at 27°C), a current associated to the protein expression was also lacking. FK‐stimulated currents and level of CFTR_(inh)‐172 inhibition were not different in both groups. The presence of SLC26A9 in Xenopus oocytes expressing CFTR also enhanced the FK‐stimulated currents as compared to oocytes expressing CFTR alone. This stimulation was mostly linked to CFTR. An enhancement of FK‐stimulated currents was not found in oocytes co‐expressing SLC26A9 and CFTR_ΔF508. In conclusion, in both protein expression systems used, SLC26A9 stimulates CFTR activity but not that of CFTR_ΔF508. Our co‐immunoprecipitation studies demonstrate a physical interaction between both anion channels. We propose as an alternative hypothesis (not exclusive) to the known SLC26A9‐STAS domain/CFTR interaction, that SLC26A9 favors the biogenesis and/or stabilization of CFTR, leading to stimulated currents. J. Cell. Physiol. 226: 212–223, 2010. © 2010 Wiley‐Liss, Inc.