Stability and potential application of spliced nuclear pre -mRNA introns

A major portion of this thesis work was dedicated to study the nature and significance of spliced introns. The initial work was focused on studying the IVS1$\sb{\rm C\beta 1}$ intron from a T-cell receptor (TCR)-$\beta$ gene. Compared to an intron lariat control from adenovirus pre-mRNA that was spliced in vitro, IVS1$\sb{\rm C\beta 1}$ was debranched less efficiently by HeLa S100 extracts, although IVS1$\sb{\rm C\beta 1}$ also used the consensus branchpoint in vivo. Subcellular-fractionation analysis showed that most IVS1$\sb{\rm C\beta 1}$ lariats cofractionated with pre-mRNA in the nucleus, consistent with the possibility that intron degradation releases splicing factors which will be available for further rounds of splicing. The half-life of IVS1$\sb{\rm C\beta 1}$ from the endogenous TCR-$\beta$ gene was measured using the general transcription inhibitor actinomycin D to be about $\sim$15 min, which was similar to that of unstable mRNAs such as c-myc mRNA.^ The general transcription inhibitor DRB was also used for intron stability analysis. Unexpectedly, DRB decreased intron and pre-mRNA levels only initially, it later increased the levels of intron-containing RNAs. Inhibition of transcription initiation appeared to be the major early effect (the reduction phase); whereas enhanced premature transcription termination was dominant later (the induction phase).^ Having established the procedures for studying in vivo spliced introns, this approach was applied to study the mechanism of nonsense-mediated downregulation (NMD), a phenomena in which premature termination codons (PTCs) decrease the levels of mRNAs. In this study, the novel intron-oriented approach was applied to study the mechanism of NMD. The levels of spliced introns immediately upstream and downstream of a PTC-bearing exon in a TCR-$\beta$ gene were identified and analyzed along with their pre-mRNA. Although PTC reduced the mRNA levels by 4 to 9 fold, the steady-state levels of spliced introns and the pre-mRNA-to-intron ratios were not significantly altered, indicating that the PTC did not significantly inhibit TCR-$\beta$ RNA splicing. Consistent with this conclusion, the half-lives of the PTC$\sp+$ and PTC$\sp-$ pre-mRNA were similar. The protein synthesis inhibitor cyclohexmide (CHX) upregulated the levels of the PTC$\sp+$ mRNA over 10 fold without affecting the levels of the spliced introns, suggesting that the reversal effect of CHX was through stabilization, not production. These results indicated that inhibition of splicing could not be the major mechanism for the NMD pathway of the TCR-$\beta$ gene, instead, suggesting that mRNA destabilization may be more important. (Abstract shortened by UMI.) ^

Genome sequence of Fusobacterium nucleatum subspecies polymorphum - a genetically tractable fusobacterium.

Fusobacterium nucleatum is a prominent member of the oral microbiota and is a common cause of human infection. F. nucleatum includes five subspecies: polymorphum, nucleatum, vincentii, fusiforme, and animalis. F. nucleatum subsp. polymorphum ATCC 10953 has been well characterized phenotypically and, in contrast to previously sequenced strains, is amenable to gene transfer. We sequenced and annotated the 2,429,698 bp genome of F. nucleatum subsp. polymorphum ATCC 10953. Plasmid pFN3 from the strain was also sequenced and analyzed. When compared to the other two available fusobacterial genomes (F. nucleatum subsp. nucleatum, and F. nucleatum subsp. vincentii) 627 open reading frames unique to F. nucleatum subsp. polymorphum ATCC 10953 were identified. A large percentage of these mapped within one of 28 regions or islands containing five or more genes. Seventeen percent of the clustered proteins that demonstrated similarity were most similar to proteins from the clostridia, with others being most similar to proteins from other gram-positive organisms such as Bacillus and Streptococcus. A ten kilobase region homologous to the Salmonella typhimurium propanediol utilization locus was identified, as was a prophage and integrated conjugal plasmid. The genome contains five composite ribozyme/transposons, similar to the CdISt IStrons described in Clostridium difficile. IStrons are not present in the other fusobacterial genomes. These findings indicate that F. nucleatum subsp. polymorphum is proficient at horizontal gene transfer and that exchange with the Firmicutes, particularly the Clostridia, is common.

Estimation of male-to-female ratios of mutation rate in primates and rodents

Haldane (1935) developed a method for estimating the male-to-female ratio of mutation rate ($\alpha$) by using sex-linked recessive genetic disease, but in six different studies using hemophilia A data the estimates of $\alpha$ varied from 1.2 to 29.3. Direct genomic sequencing is a better approach, but it is laborious and not readily applicable to non-human organisms. To study the sex ratios of mutation rate in various mammals, I used an indirect method proposed by Miyata et al. (1987). This method takes advantage of the fact that different chromosomes segregate differently between males and females, and uses the ratios of mutation rate in sequences on different chromosomes to estimate the male-to-female ratio of mutation rate. I sequenced the last intron of ZFX and ZFY genes in 6 species of primates and 2 species of rodents; I also sequenced the partial genomic sequence of the Ube1x and Ube1y genes of mice and rats. The purposes of my study in addition to estimation of $\alpha$'s in different mammalian species, are to test the hypothesis that most mutations are replication dependent and to examine the generation-time effect on $\alpha$. The $\alpha$ value estimated from the ZFX and ZFY introns of the six primate specise is ${\sim}$6. This estimate is the same as an earlier estimate using only 4 species of primates, but the 95% confidence interval has been reduced from (2, 84) to (2, 33). The estimate of $\alpha$ in the rodents obtained from Zfx and Zfy introns is ${\sim}$1.9, and that deriving from Ube1x and Ube1y introns is ${\sim}$2. Both estimates have a 95% confidence interval from 1 to 3. These two estimates are very close to each other, but are only one-third of that of the primates, suggesting a generation-time effect on $\alpha$. An $\alpha$ of 6 in primates and 2 in rodents are close to the estimates of the male-to-female ratio of the number of germ-cell divisions per generation in humans and mice, which are 6 and 2, respectively, assuming the generation time in humans is 20 years and that in mice is 5 months. These findings suggest that errors during germ-cell DNA replication are the primary source of mutation and that $\alpha$ decreases with decreasing length of generation time. ^

Characterization of the human tissue transglutaminase gene promoter

Transglutaminases are a family of calcium-dependent enzymes, that catalyze the covalent cross-linking of proteins by forming $\varepsilon(\gamma$-glutamyl)lysine isopeptide bonds. In order to investigate the molecular mechanisms regulating the expression of the tissue transglutaminase gene and to determine its biological functions, the goal of this research has been to clone and characterize the human tissue transglutaminase promoter. Thirteen clones of the tissue transglutaminase gene were obtained from the screening of a human placental genomic DNA library. A 1.74 Kb fragment derived from DNA located immediately upstream of the translation start site was subcloned and sequenced. Sequence analysis of this DNA fragment revealed that it contains a TATA box (TATAA), a CAAT box (GGACAAT), and a series of potential transcription factor binding sites and hormone response elements. Four regions of significant homology, a GC-rich region, a TG-rich region, an AG-rich region, and HR1, were identified by aligning 1.8 Kb of DNA flanking the human, mouse, and guinea pig tissue transglutaminase genes.^ To measure promoter activity, we subcloned the 1.74 Kb fragment of the tissue transglutaminase gene into a luciferase reporter vector to generate transglutaminase promoter/luciferase reporter constructs. Transfection experiments showed that this DNA segment includes a functional promoter with high constitutive activity. Deletion analysis revealed that the SP1 sites or corresponding sequences contribute to this activity. We investigated the role of DNA methylation in regulating the activity of the promoter and found that in vitro methylation of tissue transglutaminase promoter/luciferase reporter constructs suppressed their basal activity. Methylation of the promoter is inversely correlated with the expression of the tissue transglutaminase gene in vivo. These results suggest that DNA methylation may be one of the mechanisms regulating the expression of the gene. The tumor suppressor gene product p53 was also shown to inhibit the activity of the promoter, suggesting that induction of the tissue transglutaminase gene is not involved in the p53-dependent programmed cell death pathway. Although retinoids regulate the expression of the tissue transglutaminase gene in vivo, retinoid-inducible activity can not be identified in 3.7 Kb of DNA 5$\sp\prime$ to the tissue transglutaminase gene.^ The structure of the 5$\sp\prime$ end of the tissue transglutaminase gene was mapped. Alignment analysis of the human tissue transglutaminase gene with other human transglutaminases showed that tissue transglutaminase is the simplest member of transglutaminase superfamily. Transglutaminase genes show a conserved core of exons and introns but diverse N-terminuses and promoters. These observations suggest that key regulatory sequences and promoter elements have been appended upstream of the core transglutaminase gene to generate the diversity of regulated expression and regulated activity characteristic of the transglutaminase gene family. ^

Human TNF receptor p75/80 gene structure and promoter characterization

Tumor necrosis factor receptor p75/80 ((TNF-R p75/80) is a 75 kDa type 1 transmembrane protein expressed predominately on cells of hematopoietic lineage. TNF-R p75/80 belongs to the TNF receptor superfamily characterized by cysteine-rich extracellular regions composed of three to six disulfide-linked domains. In the present report, we have characterized, for the first time, the complete gene structure for human TNF-R p75/80 which spans approximately 43 kbp. The gene consists of 10 exons (ranging from 34 bp to 2.5 kbp) and 9 introns (343 bp to 19 kbp). Consensus elements for transcription factors involved in T cell development and activation were noted in the 5$\sp\prime$ flanking region including TCF-1, Ikaros, AP-1, CK-2, IL-6RE, ISRE, GAS, NF-$\kappa$B and SP1, as well as an unusually high GC content and CpG frequency that appears characteristic of some TNF-R family members. The unusual (GATA)$\sb{\rm n}$ and (GAA)(GGA) repeats found within intron 1 may prove useful for further genome analysis within the 1p36 chromosomal locus. The human TNF-R p75/80 gene structure will permit further assessment of its involvement in normal hematopoietic cell development and function, autoimmune disease, and non-random translocations in hematopoietic malignancies. The region 1.8 kb 5$\sp\prime$ of the ATG was able to drive luciferase expression when transfected into cell lines expressing TNF-R p75/80. Further characterization of the 5$\sp\prime$-regulatory region will aid in determining factors and signal transduction pathways involved in regulating TNF-R p75/80 expression. ^

myogenin: Human candidate gene and knockout mouse effect

The myogenin gene encodes an evolutionarily conserved basic helix-loop-helix transcription factor that regulates the expression of skeletal muscle-specific genes and its homozygous deletion results in mice who die of respiratory failure at birth. The histology of skeletal muscle in the myogenin null mice is reminiscent of that found in some severe congenital myopathy patients, many of whom also die of respiratory complications and provides the rationale that an aberrant human myogenin (myf4) coding region could be associated with some congenital myopathy conditions.^ With PCR, we found similarly sized amplimers for the three exons of the myogenin gene in 37 patient and 40 control samples. In contrast to the GeneBank sequence for human myogenin, we report several differences in flanking and coding regions plus an additional 659 and 498 bps in the first and second introns, respectively, in all patients and controls. We also find a novel (CA)-dinucleotide repeat in the second intron. No causative mutations were detected in the myogenin coding regions of genomic DNA from patients with severe congenital myopathy.^ Severe congenital myopathies in humans are often associated with respiratory complications and pulmonary hypoplasia. We have employed the myogenin null mouse, which lacks normal development of skeletal muscle fibers as a genetically defined severe congenital myopathy mouse model to evaluate the effect of absent fetal breathing movement on pulmonary development.^ Significant differences are observed at embryonic days E14, E17 and E20 of lung:body weight, total DNA and histologically, suggesting that the myogenin null lungs are hypoplastic. RT-PCR, in-situ immunofluorescence and EM reveal pneumocyte type II differentiation in both null and wild lungs as early as E14. However, at E14, myogenin null lungs have decreased BrdU incorporation while E17 through term, augmented cell death is detected in the myogenin null lungs, not seen in wild littermates. Absent mechanical forces appear to impair normal growth, but not maturation, of the developing lungs in myogenin null mouse.^ These investigations provide the basis for delineating the DNA sequence of the myogenin gene and and highlight the importance of skeletal muscle development in utero for normal lung organogenesis. My observation of no mutations within the coding regions of the human myogenin gene in DNA from patients with severe congenital myopathy do not support any association with this condition. ^

Global patterns of DNA polymorphism in noncoding regions in human populations

DNA sequence variation is currently a major source of data for studying human origins, evolution, and demographic history, and for detecting linkage association of complex diseases. In this dissertation, I investigated DNA variation in worldwide populations from two ∼10 kb autosomal regions on 22q11.2 (noncoding) and 1q24 (introns). A total of 75 variant sites were found among 128 human sequences in the 22q11.2 region, yielding an estimate of 0.088% for nucleotide diversity (π), and a total of 52 variant sites were found among 122 human sequences in the 1q24 region with an estimated π value of 0.057%. The data from these two regions and a 10 kb noncoding region on Xq13.3 all show a strong excess of low-frequency variants in comparison to that expected from an equilibrium population, indicating a relatively recent population expansion. The effective population sizes estimated from the three regions were 11,000, 12,700, and 8,600, respectively, which are close to the commonly used value of 10,000. In each of the two autosomal regions, the age of the most recent common ancestor (MRCA) was estimated to be older than 1 million years among all the sequences and ∼600,000 years among non-African sequences, providing first evidence from autosomal noncoding or intronic regions for a genetic history of humans much more ancient than the emergence of modern humans. The ancient genetic history of humans indicates no severe bottleneck during the evolution of humans in the last half million years; otherwise, much of the ancient genetic history would have been lost during a severe bottleneck. This study strongly suggests that both the “out of Africa” and the multiregional models are too simple for explaining the evolution of modern humans. A compilation of genome-wide data revealed that nucleotide diversity is highest in autosomal regions, intermediate in X-linked regions, and lowest in Y-linked regions. The data suggest the existence of background selection or selective sweep on Y-linked loci. In general, the nucleotide diversity in humans is low compared to that in chimpanzee and Drosophila populations. ^

The roles of UPF3a and UPF3b in nonsense mediated decay

Nonsense-mediated decay (NMD) degrades aberrant transcripts containing premature termination codons (PTCs). The T-cell receptor (TCR) locus undergoes error-prone rearrangements that frequently acquire PTCs. Transcripts harboring PTCs from this locus are downregulated much more than transcripts from non-rearranging genes. Efficient splicing is essential for this robust downregulation. ^ Here I show that TCR NMD is unique in another respect: it is not impaired by RNAi-mediated depletion of the NMD factor UPF3b. This differentiates TCR transcripts from classical NMD (assayed using β-globin or triose phosphate isomerase transcripts), which does depend on UPF3b. Depletion of UPF3a, which encodes a gene related to UPF3b, also had no effect on TCR NMD. Mapping experiments identified TCR sequences that when deleted or mutated caused a switch to UPF3b dependence. Since UPF3b dependence was invariably accompanied by less efficient RNA splicing, this suggests that UPF3b-dependent NMD occurs when transcripts are generated by inefficient splicing. Microarray analysis revealed the existence of many NMD-targeted mRNAs from wild-type genes whose downregulation is impervious to UPF3b depletion. This suggests the existence of an alternative NMD pathway independent of UPF3b that is widely used to downregulate the level of both normal and mutant transcripts. ^ During the course of my studies, I also found that the function of UPF3a is fundamentally distinct from that of UPF3b in several aspects. First, classical NMD failed to be impaired by UPF3a depletion, whereas it was reversed by UPF3b depletion. Second, UPF3a depletion had no effect on NMD elicited by tethered UPF2, whereas UPF3b depletion blocked this response. Thus, UPF3a does not function in classical NMD. Third, UPF3b depletion upregulated the expression of UPF3a, whereas UPF3a depletion had no effect on UPF3b expression. This suggests that a UPF3b-mediated feedback network exists that regulates the UPF3a expression. Lastly, UPF3a depletion but not UPF3b depletion significantly upregulated TCR precursor RNAs. This suggests that UPF3a, not UPF3b, functions in the surveillance of precursor RNAs, which typically contain many PTCs in the introns. Collectively, my data suggests that UPF3a and UPF3b are not functionally redundant, as previously thought, but instead have separable functions. ^

Regulation of RNA splicing by nonsense mutations

Translation termination as a result of premature nonsense codon-incorporation in a RNA transcript can lead to the production of aberrant proteins with gain-of-function or dominant negative properties that could have deletrious effects on the cell. T-cell Receptor (TCR) genes acquire premature termination codons two-thirds of the time as a result of the error-prone programmed rearrangement events that normally occur during T-cell development. My studies have focused on the fate of TCR precursor mRNAs in response to in-frame nonsense mutations. ^ Previous published studies from our laboratory have shown that TCR precursor mRNAs are subject to nonsense mediated upregulation of pre-mRNA (NMUP). In this dissertation, I performed substitution and deletion analysis to characterize specific regions of TCR which are required to elicit NMUP. I performed frame- and factor-dependence studies to determine its relationship with other nonsense codon induced responses using several approaches including (i) translation dependence studies (ii) deletion and mutational analysis, as well as (iii) siRNA mediated knockdown of proteins involved. I also addressed the underlying molecular mechanism for this pre-mRNA upregulation by (i) RNA half-life studies using a c-fos inducible promoter, and (ii) a variety of assays to determine pre-mRNA splicing efficiency. ^ Using these approaches, I have identified a region of TCR that is both necessary and sufficient to elicit (NMUP). I have also found that neither cytoplasmic translation machinery nor the protein UPF1 are involved in eliciting this nuclear event. I have shown that the NMUP can be induced not only by nonsense and frameshift mutations, but also missense mutations that disrupt a cis splicing element in the exon that contains the mutation. However, the effect of nonsense mutations on pre-mRNA is unique and distinguishable from that of missense mutations in that nonsense mutations can upregulate pre-mRNA in a frame-dependent manner. Lastly, I provide evidence that NMUP occurs by a mechanism in which nonsense mutations inhibit the splicing of introns. In summary, I have found that TCR precursor mRNAs are subject to multiple forces involving both RNA splicing and translation that can either increase or decrease the levels of these precursor mRNAs. ^

Functional Analysis of Drosophila Integrator Complex in snRNA 3' End Processing

Uridine-rich small nuclear RNAs (U snRNAs) play essential roles in eukaryotic gene expression by facilitating the removal of introns from mRNA precursors and the processing of the replication-dependent histone pre-mRNAs. Formation of the 3’ end of these snRNAs is carried out by a poorly characterized, twelve-membered protein complex named Integrator Complex. In the effort to understand Integrator Complex function in the formation of the snRNA 3’ end, we performed a functional RNAi screen in Drosophila S2 cells to identify protein factors required for snRNA 3’ end formation. This screen was conducted by using a fluorescence-based reporter that elicits GFP expression in response to a deficiency in snRNA processing. Besides scoring the known Integrator subunits, we identified Asunder and CG4785 as additional core members of the Integrator Complex. Additionally, we also found a conserved requirement for Cyclin C and Cdk8 in both fly and human snRNA 3’ end processing. We have further demonstrated that the kinase activity of Cdk8 is critical for snRNA 3’ end processing and is likely to function independent of its well-documented function within the Mediator Cdk8 module. Taken together, this work functionally defines the Drosophila Integrator Complex and demonstrates a novel function for Cyclin C/Cdk8 in snRNA 3’ end formation. This thesis work has also characterized an important functional interaction mediated by a microdomain within Integrator subunit 12 (IntS12) and IntS1 that is required for the activity of the Integrator Complex in processing the snRNA 3’ end. Through the development of a reporter-based functional RNAi-rescue assay in Drosophila S2 cells, we analyzed domains within IntS12 required for snRNA 3’ end formation. This analysis unexpectedly revealed that an N-terminal 30 amino acid region and not the highly conserved central PHD finger domain, is required for snRNA 3’ end cleavage. The IntS12 microdomain (1-45) functions autonomously, and is sufficient to interact and stabilize the putative scaffold protein IntS1. Our findings provide more details of the Integrator Complex for understanding the molecular mechanism of snRNA 3’ end processing. Moreover, these results lay the foundation for future studies of the complex through the identification of a novel functional domain within one subunit and the identification of additional subunits.

Molecular and evolutionary genetics of the X -linked visual pigment genes in humans and New World monkeys

Normal humans have one red and at least one green visual pigment genes. These genes are tightly linked as tandem repeats on the X chromosome and each of them has six exons. There is only one X-linked visual pigment gene in New World monkeys (NWMs) but the locus has three polymorphic alleles encoding red, yellow and green visual pigments, respectively. The spectral properties of the squirrel monkey and the marmoset (both NWMs) have been studied and partial sequences of the three alleles are available. To study the evolutionary history of these X-linked opsin genes in humans and NWMs, coding and intron sequences of the three squirrel monkey alleles and the three marmoset alleles were amplified by PCR followed by subcloning and sequencing. Introns 2 and 4 of the human red and green pigment genes were also sequenced. The results obtained are as follows: (1) The sequences of introns 2 and 4 of the human red and green opsin genes are significantly more similar between the two genes than are coding sequences, contrary to the usual situation where coding regions are better conserved in evolution than are introns. The high similarities in the two introns are probably due to recent gene conversion events during evolution of the human lineage. (2) Phylogenetic analysis of both intron and exon sequences indicates that the phylogenetic tree of the available primate opsin genes is the same as the species tree. The two human genes were derived from a gene duplication event after the divergence of the human and NWM lineages. The three alleles in each of the two NWM species diverged after the split of the two NWMs but have persisted in the population for at least 5 million years. (3) Allelic gene conversion might have occurred between the three squirrel monkey alleles. (4) A model of additive effect of hydroxyl-bearing amino acids on spectral tuning is proposed by treating some unknown variables as groups. Under the assumption that some residues have no effect, it is found that at least five amino acid residues, at positions 178 (3 nm), 180 (5 nm), 230 ($-$4 nm), 277 (9 nm) and 285 (13 nm), have linear spectral tuning effects. (5) Adaptive evolution of the opsin genes to different spectral peaks was observed at four residues that are important for spectral tuning. ^