Fine-mapping and mutation analysis of TRPM1: a candidate gene for leopard complex (LP) spotting and congenital stationary night blindness in horses.

Leopard Complex spotting occurs in several breeds of horses and is caused by an incompletely dominant allele (LP). Homozygosity for LP is also associated with congenital stationary night blindness (CSNB) in Appaloosa horses. Previously, LP was mapped to a 6 cm region on ECA1 containing the candidate gene TRPM1 (Transient Receptor Potential Cation Channel, Subfamily M, Member 1) and decreased expression of this gene, measured by qRT-PCR, was identified as the likely cause of both spotting and ocular phenotypes. This study describes investigations for a mutation causing or associated with the Leopard Complex and CSNB phenotype in horses. Re-sequencing of the gene and associated splice sites within the 105 624 bp genomic region of TRPM1 led to the discovery of 18 SNPs. Most of the SNPs did not have a predictive value for the presence of LP. However, one SNP (ECA1:108,249,293 C>T) found within intron 11 had a strong (P < 0.0005), but not complete, association with LP and CSNB and thus is a good marker but unlikely to be causative. To further localize the association, 70 SNPs spanning over two Mb including the TRPM1 gene were genotyped in 192 horses from three different breeds segregating for LP. A single 173 kb haplotype associated with LP and CSNB (ECA1: 108,197,355- 108,370,150) was identified. Illumina sequencing of 300 kb surrounding this haplotype revealed 57 SNP variants. Based on their localization within expressed sequences or regions of high sequence conservation across mammals, six of these SNPs were considered to be the most likely candidate mutations. While the precise function of TRPM1 remains to be elucidated, this work solidifies its functional role in both pigmentation and night vision. Further, this work has identified several potential regulatory elements of the TRPM1 gene that should be investigated further in this and other species.

Stable alteration of pre-mRNA splicing patterns by modified U7 small nuclear RNAs.

In several forms of beta-thalassemia, mutations in the second intron of the beta-globin gene create aberrant 5' splice sites and activate a common cryptic 3' splice site upstream. As a result, the thalassemic beta-globin pre-mRNAs are spliced almost exclusively via the aberrant splice sites leading to a deficiency of correctly spliced beta-globin mRNA and, consequently, beta-globin. We have designed a series of vectors that express modified U7 snRNAs containing sequences antisense to either the aberrant 5' or 3' splice sites in the IVS2-705 thalassemic pre-mRNA. Transient expression of modified U7 snRNAs in a HeLa cell line stably expressing the IVS2-705 beta-globin gene restored up to 65% of correct splicing in a sequence-specific and dose-dependent manner. Cell lines that stably coexpressed IVS2-705 pre-mRNA and appropriately modified U7 snRNA exhibited up to 55% of permanent restoration of correct splicing and expression of full-length beta-globin protein. This novel approach provides a potential alternative to gene replacement therapies.

The TREAT-NMD DMD Global Database: analysis of more than 7,000 Duchenne muscular dystrophy mutations.

Analyzing the type and frequency of patient-specific mutations that give rise to Duchenne muscular dystrophy (DMD) is an invaluable tool for diagnostics, basic scientific research, trial planning, and improved clinical care. Locus-specific databases allow for the collection, organization, storage, and analysis of genetic variants of disease. Here, we describe the development and analysis of the TREAT-NMD DMD Global database (http://umd.be/TREAT_DMD/). We analyzed genetic data for 7,149 DMD mutations held within the database. A total of 5,682 large mutations were observed (80% of total mutations), of which 4,894 (86%) were deletions (1 exon or larger) and 784 (14%) were duplications (1 exon or larger). There were 1,445 small mutations (smaller than 1 exon, 20% of all mutations), of which 358 (25%) were small deletions and 132 (9%) small insertions and 199 (14%) affected the splice sites. Point mutations totalled 756 (52% of small mutations) with 726 (50%) nonsense mutations and 30 (2%) missense mutations. Finally, 22 (0.3%) mid-intronic mutations were observed. In addition, mutations were identified within the database that would potentially benefit from novel genetic therapies for DMD including stop codon read-through therapies (10% of total mutations) and exon skipping therapy (80% of deletions and 55% of total mutations).

Regulation of mRNA splicing in Moloney murine sarcoma virus mutants

Cells infected with the conditionally defective MuSVts110 mutant of Moloney murine sarcoma virus are transformed at 33$\sp\circ$C but appear morphologically normal at 39$\sp\circ$C. The molecular basis for this phenotype is as follows: MuSVts110 contains a 1487 nucleotide central deletion that has truncated the 3$\sp\prime$ end to the gag gene and the 5$\sp\prime$ end of the mos gene. The resulting gag-mos junction is out-of-frame and the v-mos protein is not expressed. At 33$\sp\circ$C or lower, a splicing event is activated such that a 431 base intron is removed to realign the gag and mos gene in-frame, allowing the expression of a transforming protein P85$\sp{gag-mos}$. Temperature-dependent splicing appeared to be an intrinsic property of MuSVts110 transcripts and not a general feature of pre-mRNA splicing in 6m2 cells since splicing activity of a heterologous transcript in the same cells did not vary with temperature. The possibility that the splice event was not temperature-sensitive, but that the accumulation of spliced transcript at the lower growth temperatures was due to its selective thermolability was ruled out as stability studies revealed that the relative turnover rates of the unspliced and spliced MuSVts110 transcripts were not affected by temperature.^ The consensus sequences containing the splice sites activated in the MuSVts110 mutant (5$\sp\prime$ gag and 3$\sp\prime$ mos) are present, but not utilized, in wild-type MuSV-124. To test the hypothesis that it was the reduction of the 1919 base intervening sequence in MuSV-124 to 431 bases in MuSVts110 which activated splicing, the identical 1487 base deletion was introduced into cloned wild-type MuSV-124 DNA to create the MuSVts110 equivalent, ts32.^ To examine conditions permissive for splicing, we assayed splice site activation in a series of MuSV-124 "intron-modification" mutants. Data suggest that splicing in wild-type MuSV-124 may be blocked due to the lack of a proximal branchpoint sequence, but can be activated by those intron mutations which reposition a branch site closer to the 3$\sp\prime$ splice site. (Abstract shortened with permission of author.) ^

Repair of thalassemic human β-globin mRNA in mammalian cells by antisense oligonucleotides

In one form of β-thalassemia, a genetic blood disorder, a mutation in intron 2 of the β-globin gene (IVS2-654) causes aberrant splicing of β-globin pre-mRNA and, consequently, β-globin deficiency. Treatment of mammalian cells stably expressing the IVS2-654 human β-globin gene with antisense oligonucleotides targeted at the aberrant splice sites restored correct splicing in a dose-dependent fashion, generating correct human β-globin mRNA and polypeptide. Both products persisted for up to 72 hr posttreatment. The oligonucleotides modified splicing by a true antisense mechanism without overt unspecific effects on cell growth and splicing of other pre-mRNAs. This novel approach in which antisense oligonucleotides are used to restore rather than to down-regulate the activity of the target gene is applicable to other splicing mutants and is of potential clinical interest.

Development of an epithelium-specific expression cassette with human DNA regulatory elements for transgene expression in lung airways

The efficient expression of therapeutic genes in target cells or tissues is an important component of efficient and safe gene therapy. Utilizing regulatory elements from the human cytokeratin 18 (K18) gene, including 5′ genomic sequences and one of its introns, we have developed a novel expression cassette that can efficiently express reporter genes, as well as the human cystic fibrosis transmembrane conductance regulator (CFTR) gene, in cultured lung epithelial cells. CFTR transcripts expressed from the native K18 enhancer/promoter include two alternative splicing products, due to the activation of two cryptic splice sites in the CFTR coding region. Modification of the K18 intron and CFTR cDNA sequences eliminated the cryptic splice sites without changing the CFTR amino acid sequence, and led to enhanced CFTR mRNA and protein expression as well as biological function. Transgenic expression analysis in mice showed that the modified expression cassette can direct efficient and epithelium-specific expression of the Escherichia coli LacZ gene in the airways of fetal lungs, with no detectable expression in lung fibroblasts or endothelial cells. This is the first expression cassette which selectively directs lung transgene expression for CFTR gene therapy to airway epithelia.

Ultrastructural Analysis of Transcription and Splicing in the Cell Nucleus after Bromo-UTP Microinjection

In this study we demonstrate, at an ultrastructural level, the in situ distribution of heterogeneous nuclear RNA transcription sites after microinjection of 5-bromo-UTP (BrUTP) into the cytoplasm of living cells and subsequent postembedding immunoelectron microscopic visualization after different labeling periods. Moreover, immunocytochemical localization of several pre-mRNA transcription and processing factors has been carried out in the same cells. This high-resolution approach allowed us to reveal perichromatin regions as the most important sites of nucleoplasmic RNA transcription and the perichromatin fibrils (PFs) as in situ forms of nascent transcripts. Furthermore, we show that transcription takes place in a rather diffuse pattern, without notable local accumulation of transcription sites. RNA polymerase II, heterogeneous nuclear ribonucleoprotein (hnRNP) core proteins, general transcription factor TFIIH, poly(A) polymerase, splicing factor SC-35, and Sm complex of small nuclear ribonucleoproteins (snRNPs) are associated with PFs. This strongly supports the idea that PFs are also sites of major pre-mRNA processing events. The absence of nascent transcripts, RNA polymerase II, poly(A) polymerase, and hnRNPs within the clusters of interchromatin granules rules out the possibility that this domain plays a role in pre-mRNA transcription and polyadenylation; however, interchromatin granule-associated zones contain RNA polymerase II, TFIIH, and Sm complex of snRNPs and, after longer periods of BrUTP incubation, also Br-labeled RNA. Their role in nuclear functions still remains enigmatic. In the nucleolus, transcription sites occur in the dense fibrillar component. Our fine structural results show that PFs represent the major nucleoplasmic structural domain involved in active pre-mRNA transcriptional and processing events.

Architectural limits on split genes

Exon/intron architecture varies across the eukaryotic kingdom with large introns and small exons the rule in vertebrates and the opposite in lower eukaryotes. To investigate the relationship between exon and intron size in pre-mRNA processing, internally expanded exons were placed in vertebrate genes with small and large introns. Both exon and intron size influenced splicing phenotype. Intron size dictated if large exons were efficiently recognized. When introns were large, large exons were skipped; when introns were small, the same large exons were included. Thus, large exons were incompatible for splicing if and only if they were flanked by large introns. Both intron and exon size became problematic at ≈500 nt, although both exon and intron sequence influenced the size at which exons and introns failed to be recognized. These results indicate that present-day gene architecture reflects at least in part limitations on exon recognition. Furthermore, these results strengthen models that invoke pairing of splice sites during recognition of pre-mRNAs, and suggest that vertebrate consensus sequences support pairing across either introns or exons.

The SRm160/300 splicing coactivator is required for exon-enhancer function

Exonic splicing enhancer (ESE) sequences are important for the recognition of splice sites in pre-mRNA. These sequences are bound by specific serine-arginine (SR) repeat proteins that promote the assembly of splicing complexes at adjacent splice sites. We have recently identified a splicing “coactivator,” SRm160/300, which contains SRm160 (the SR nuclear matrix protein of 160 kDa) and a 300-kDa nuclear matrix antigen. In the present study, we show that SRm160/300 is required for a purine-rich ESE to promote the splicing of a pre-mRNA derived from the Drosophila doublesex gene. The association of SRm160/300 and U2 small nuclear ribonucleoprotein particle (snRNP) with this pre-mRNA requires both U1 snRNP and factors bound to the ESE. Independently of pre-mRNA, SRm160/300 specifically interacts with U2 snRNP and with a human homolog of the Drosophila alternative splicing regulator Transformer 2, which binds to purine-rich ESEs. The results suggest a model for ESE function in which the SRm160/300 splicing coactivator promotes critical interactions between ESE-bound “activators” and the snRNP machinery of the spliceosome.

Protein–RNA interactions: a structural analysis

A detailed computational analysis of 32 protein–RNA complexes is presented. A number of physical and chemical properties of the intermolecular interfaces are calculated and compared with those observed in protein–double-stranded DNA and protein–single-stranded DNA complexes. The interface properties of the protein–RNA complexes reveal the diverse nature of the binding sites. van der Waals contacts played a more prevalent role than hydrogen bond contacts, and preferential binding to guanine and uracil was observed. The positively charged residue, arginine, and the single aromatic residues, phenylalanine and tyrosine, all played key roles in the RNA binding sites. A comparison between protein–RNA and protein–DNA complexes showed that whilst base and backbone contacts (both hydrogen bonding and van der Waals) were observed with equal frequency in the protein–RNA complexes, backbone contacts were more dominant in the protein–DNA complexes. Although similar modes of secondary structure interactions have been observed in RNA and DNA binding proteins, the current analysis emphasises the differences that exist between the two types of nucleic acid binding protein at the atomic contact level.

The human XPG gene: gene architecture, alternative splicing and single nucleotide polymorphisms

Defects in the XPG DNA repair endonuclease gene can result in the cancer-prone disorders xeroderma pigmentosum (XP) or the XP–Cockayne syndrome complex. While the XPG cDNA sequence was known, determination of the genomic sequence was required to understand its different functions. In cells from normal donors, we found that the genomic sequence of the human XPG gene spans 30 kb, contains 15 exons that range from 61 to 1074 bp and 14 introns that range from 250 to 5763 bp. Analysis of the splice donor and acceptor sites using an information theory-based approach revealed three splice sites with low information content, which are components of the minor (U12) spliceosome. We identified six alternatively spliced XPG mRNA isoforms in cells from normal donors and from XPG patients: partial deletion of exon 8, partial retention of intron 8, two with alternative exons (in introns 1 and 6) and two that retained complete introns (introns 3 and 9). The amount of alternatively spliced XPG mRNA isoforms varied in different tissues. Most alternative splice donor and acceptor sites had a relatively high information content, but one has the U12 spliceosome sequence. A single nucleotide polymorphism has allele frequencies of 0.74 for 3507G and 0.26 for 3507C in 91 donors. The human XPG gene contains multiple splice sites with low information content in association with multiple alternatively spliced isoforms of XPG mRNA.

Isolation of a Pax-6 homolog from the ribbonworm Lineus sanguineus.

The Pax-6 genes of vertebrates and Drosophila encode transcription factors with highly conserved paired- and homeodomains. They are expressed in the nervous system and the developing eyes. Loss-of-function mutations in mammals and flies lead to a reduction or absence of the eyes. By ectopic expression of Pax-6 in Drosophila ectopic eyes can be induced, indicating a determinative role in eye morphogenesis. We have isolated a Pax-6 homolog of the ribbonworm Lineus sanguineus. This gene shares extensive sequence identity and several conserved splice sites with the mammalian and Drosophila genes. During head regeneration the L. sanguineus Pax-6 homolog is expressed in the central nervous system, in the cerebral organ, and in the eye region. These findings support the hypothesis that Pax-6 was present in primitive metazoa before the evolutionary separation of vertebrates and arthropods and suggest a fundamental role in eye and central nervous system development.

The 3'-terminal exon of the family of steroid and phenol sulfotransferase genes is spliced at the N-terminal glycine of the universally conserved GXXGXXK motif that forms the sulfonate donor binding site.

The guinea pig estrogen sulfotransferase gene has been cloned and compared to three other cloned steroid and phenol sulfotransferase genes (human estrogen sulfotransferase, human phenol sulfotransferase, and guinea pig 3 alpha-hydroxysteroid sulfotransferase). The four sulfotransferase genes demonstrate a common outstanding feature: the splice sites for their 3'-terminal exons are identically located. That is, the 3'-terminal exon splice sites involve a glycine that constitutes the N-terminal glycine of an invariably conserved GXXGXXK motif present in all steroid and phenol sulfotransferases for which primary structures are known. This consistency strongly suggests that all steroid and phenol sulfotransferase genes will be similarly spliced. The GXXGXXK motif forms the active binding site for the universal sulfonate donor 3'-phosphoadenosine 5'-phosphosulfate. Amino acid sequence alignment of 19 cloned steroid and phenol sulfotransferases starting with the GXXGXXK motif indicates that the 3'-terminal exon for each steroid and phenol sulfotransferase gene encodes a similarly sized C-terminal fragment of the protein. Interestingly, on further analysis of the alignment, three distinct amino acid sequence patterns emerge. The presence of the conserved functional GXXGXXK motif suggests that the protein domains encoded by steroid and phenol sulfotransferase 3'-terminal exons have evolved from a common ancestor. Furthermore, it is hypothesized that during the course of evolution, the 3'-terminal exon further diverged into at least three sulfotransferase subdivisions: a phenol or aryl group, an estrogen or phenolic steroid group, and a neutral steroid group.

Isolation and characterization of the U1 and U2 snRNAs in the silk gland of the 5th instar larval stage silk moth Bombyx mori

In isolation and characterization studies, expression level U1 and U2 snRNA isoforms were obtained from the 5th instar larval stage silk gland (SG). The DNA content of the SG cells is approximately 200,000-fold higher compared to the usual (2N) somatic cells of B. mori due to endoreduplication. In this study, the existence of U1 and U2 snRNA isoforms in the SG of the organism is investigated. Bombyx mori U1 and U2-specific RT-PCR libraries from the silk gland were generated. Five U1 and eight U2 isoforms were isolated and characterized. Nucleotide differences, structural alterations, as well as protein and RNA interaction sites were analyzed in these variants. For the U1 snRNA variants, they were compared to the previously reported BmN isoforms. In all these U-snRNA variants, polymorphic sites do not predominate at the core of known functional sequences, which were interspecifically conserved. Variant sites and inter-species differences are located in moderately conserved regions. Free energy (ΔG) values for the entire U1 and U2 snRNA secondary structures and for the individual stem/loops domains of the isoforms were generated and compared to determine their structural stability. This will be the first time that U1 and U2 variants are shown specific for a development stage (larval) other than embryonic or adult. ^ Using phylogenetic analysis, evolutionary trees were generated for the U1 and U2 snRNAs using animal, plant, protista and fungal species. The resulting trees were boostrapped for robustness and rooted with the self-splicing RNA group II intron sequence from the cyanobacterium Calothrix. Using phylogenetic analyses, possible structural and functional evolutionary interdependence between the U1 and U2 snRNAs was investigated. ^