Cocaine induces long-lasting changes in mesolimbic ‘reward’ circuits of the brain after cessation of use. These lingering changes include the neuronal plasticity that is thought to underlie the chronic relapsing nature of substance use disorders. Genes involved in neuronal plasticity also encode circular RNAs (circRNAs), which are stable, non-coding RNAs formed through the back-splicing of pre-mRNA. The Homer1 gene family, which encodes proteins associated with cocaine-induced plasticity, also encodes circHomer1. Based on preliminary evidence from shows cocaine-regulated changes in the ratio of circHomer1 and Homer1b mRNA in the nucleus accumbens (NAc), this study examined the relationship between circHomer1 and incentive motivation for cocaine by using different lengths of abstinence to vary the degree of motivation. Male and female rats were trained to self-administer cocaine (0.75 mg/kg/infusion, IV) or received a yoked saline infusion. Rats proceeded on an increasingly more difficult variable ratio schedule of lever pressing until they reached a variable ratio 5 schedule, which requires an average of 5 lever presses, and light and tone cues were delivered with the drug infusions. Rats were then tested for cocaine-seeking behavior in response to cue presentations without drug delivery either 1 or 21 days after their last self-administration session. They were sacrificed immediately after and circHomer1 and Homer1b expression was then measured from homogenate and synaptosomal fractions of NAc shell using RT-qPCR. Lever pressing during the cue reactivity test increased from 1 to 21 days of abstinence as expected. Results showed no group differences in synaptic circHomer1 expression, however, total circHomer1 expression was downregulated in 21d rats compared to controls. Lack of change in synaptic circHomer1 was likely due to trends toward different temporal changes in males versus females. Total Homer1b expression was higher in females, although there was no effect of cocaine abstinence. Further research investigating the time course of circHomer1 and Homer1b expression is warranted based on the inverse relationship between total circHomer1and cocaine-seeking behavior observed in this study.

Advances in sequencing technology have generated an enormous amount of data over the past decade. Equally advanced computational methods are needed to conduct comparative and functional genomic studies on these datasets, in particular tools that appropriately interpret indels within an evolutionary framework. The evolutionary history of indels is complex and often involves repetitive genomic regions, which makes identification, alignment, and annotation difficult. While previous studies have found that indel lengths in both deoxyribonucleic acid and proteins obey a power law, probabilistic models for indel evolution have rarely been explored due to their computational complexity. In my research, I first explore an application of an expectation-maximization algorithm for maximum-likelihood training of a codon substitution model. I demonstrate the training accuracy of the expectation-maximization on my substitution model. Then I apply this algorithm on a published 90 pairwise species dataset and find a negative correlation between the branch length and non-synonymous selection coefficient. Second, I develop a post-alignment fixation method to profile each indel event into three different phases according to its codon position. Because current codon-aware models can only identify the indels by placing the gaps between codons and lead to the misalignment of the sequences. I find that the mouse-rat species pair is under purifying selection by looking at the proportion difference of the indel phases. I also demonstrate the power of my sliding-window method by comparing the post-aligned and original gap positions. Third, I create an indel-phase moore machine including the indel rates of three phases, length distributions, and codon substitution models. Then I design a gillespie simulation that is capable of generating true sequence alignments. Next I develop an importance sampling method within the expectation-maximization algorithm that can successfully train the indel-phase model and infer accurate parameter estimates from alignments. Finally, I extend the indel phase analysis to the 90 pairwise species dataset across three alignment methods, including Mafft+sw method developed in chapter 3, coati-sampling methods applied in chapter 4, and coati-max method. Also I explore a non-linear relationship between the dN/dS and Zn/(Zn+Zs) ratio across 90 species pairs.

Type 1 diabetes (T1D) is the result of an autoimmune attack against the insulin-producing β-cells of the pancreas causing hyperglycemia and requiring the individual to rely on life-long exogenous insulin. With the age of onset typically occurring in childhood, there is increased physical and emotional stress to the child as well as caregivers to maintain appropriate glucose levels. The majority of T1D patients have antibodies to one or more antigens: insulin, IA-2, GAD65, and ZnT8. Although antibodies are detectable years before symptoms occur, the initiating factors and mechanisms of progression towards β-cell destruction are still not known. The search for new autoantibodies to elucidate the autoimmune process in diabetes has been slow, with proteome level screenings on native proteins only finding a few minor antigens. Post-translational modifications (PTM)—chemical changes that occur to the protein after translation is complete—are an unexplored way a self-protein could become immunogenic. This dissertation presents the first large sale screening of autoantibodies in T1D to nitrated proteins. The Contra Capture Protein Array (CCPA) allowed for fresh expression of hundreds of proteins that were captured on a secondary slide by tag-specific ligand and subsequent modification with peroxynitrite. The IgG and IgM humoral response of 48 newly diagnosed T1D subjects and 48 age-matched controls were screened against 1632 proteins highly or specifically expressed in pancreatic cells. Top targets at 95% specificity were confirmed with the same serum samples using rapid antigenic protein in situ display enzyme-linked immunosorbent assay (RAPID ELISA) a modified sandwich ELISA employing the same cell-free expression as the CCPA. For validation, 8 IgG and 5 IgM targets were evaluated with an independent serum sample set of 94 T1D subjects and 94 controls. The two best candidates at 90% specificity were estrogen receptor 1 (ESR1) and phosphatidylinositol 4-kinase type 2 beta (PI4K2B) which had sensitivities of 22% (p=.014) and 25% (p=.045), respectively. Receiver operating characteristic (ROC) analyses found an area under curve (AUC) of 0.6 for ESR1 and 0.58 for PI4K2B. These studies demonstrate the ability and value for high-throughput autoantibody screening to modified antigens and the frequency of Type 1 diabetes.

Glioblastoma (GBM), the most common and aggressive primary brain tumor affecting adults, is characterized by an aberrant yet druggable epigenetic landscape. The Histone Deacetylases (HDACs), a major family of epigenetic regulators, favor transcriptional repression by mediating chromatin compaction and are frequently overexpressed in human cancers, including GBM. Hence, over the last decade there has been considerable interest in using HDAC inhibitors (HDACi) for the treatment of malignant primary brain tumors. However, to date most HDACi tested in clinical trials have failed to provide significant therapeutic benefit to patients with GBM. This is because current HDACi have poor or unknown pharmacokinetic profiles, lack selectivity towards the different HDAC isoforms, and have narrow therapeutic windows. Isoform selectivity for HDACi is important given that broad inhibition of all HDACs results in widespread toxicity across different organs. Moreover, the functional roles of individual HDAC isoforms in GBM are still not well understood. Here, I demonstrate that HDAC1 expression increases with brain tumor grade and is correlated with decreased survival in GBM. I find that HDAC1 is the essential HDAC isoform in glioma stem cells and its loss is not compensated for by its paralogue HDAC2 or other members of the HDAC family. Loss of HDAC1 alone has profound effects on the glioma stem cell phenotype in a p53-dependent manner and leads to significant suppression of tumor growth in vivo. While no HDAC isoform-selective inhibitors are currently available, the second-generation HDACi quisinostat harbors high specificity for HDAC1. I show that quisinostat exhibits potent growth inhibition in multiple patient-derived glioma stem cells. Using a pharmacokinetics- and pharmacodynamics-driven approach, I demonstrate that quisinostat is a brain-penetrant molecule that reduces tumor burden in flank and orthotopic models of GBM and significantly extends survival both alone and in combination with radiotherapy. The work presented in this thesis thereby unveils the non-redundant functions of HDAC1 in therapy- resistant glioma stem cells and identifies a brain-penetrant HDACi with higher selectivity towards HDAC1 as a potent radiosensitizer in preclinical models of GBM. Together, these results provide a rationale for developing quisinostat as a potential adjuvant therapy for the treatment of GBM.

The partitioning of photosynthates between their sites of production (source) and their sites of utilization (sink) is a major determinant of crop yield and the potential of regulating this translocation promises substantial opportunities for yield increases. Ubiquitous overexpression of the plant type I proton pyrophosphatase (H+-PPase) in crops improves several valuable traits including salt tolerance and drought resistance, nutrient and water use efficiencies, and increased root biomass and yield. Originally, type I H+-PPases were described as pyrophosphate (PPi)-dependent proton pumps localized exclusively in vacuoles of mesophyll and meristematic tissues. It has been proposed that in the meristematic tissues, the role of this enzyme would be hydrolyzing PPi originated in biosynthetic reactions and favoring sink strength. Interestingly, this enzyme has been also localized at the plasma membrane of companion cells in the phloem which load and transport photosynthates from source leaves to sinks. Of note, the plasma membrane-localized H+-PPase could only function as a PPi-synthase in these cells due to the steep proton gradient between the apoplast and cytosol. The generated PPi would favor active sucrose loading through the sucrose/proton symporter in the phloem by promoting sucrose hydrolysis through the Sucrose Synthase pathway and providing the ATP required to maintain the proton gradient. To better understand these two different roles of type I H+-PPases, a series of Arabidopsis thaliana transgenic plants were generated. By expressing soluble pyrophosphatases in companion cells of Col-0 ecotype and H+-PPase mutants, impaired photosynthates partitioning was observed, suggesting phloem-localized H+-PPase could generate the PPi required for sucrose loading. Col-0 plants expressed with either phloem- or meristem-specific AVP1 overexpression cassette and the cross between the two tissue specific lines (Cross) were generated. The results showed that the phloem-specific AVP1-overexpressing plants had increased root hair elongation under limited nutrient conditions and both phloem- and meristem-overexpression of AVP1 contributed to improved rhizosphere acidification and drought resistance. It was concluded that H+-PPases localized in both sink and source tissues regulate plant growth and performance under stress through its versatile enzymatic functions (PPi hydrolase and synthase).

The RNA editing enzyme adenosine deaminase acting on double stranded RNA 2 (ADAR2) converts adenosine into inosine in regions of double stranded RNA. Here, it was discovered that this critical function of ADAR2 was dysfunctional in amyotrophic lateral sclerosis (ALS) mediated by the C9orf72 hexanucleotide repeat expansion, the most common genetic abnormality associated with ALS. Typically a nuclear protein, ADAR2 was localized in cytoplasmic accumulations in postmortem tissue from C9orf72 ALS patients. The mislocalization of ADAR2 was confirmed using immunostaining in a C9orf72 mouse model and motor neurons differentiated from C9orf72 patient induced pluripotent stem cells. Notably, the cytoplasmic accumulation of ADAR2 coexisted in neurons with cytoplasmic accumulations of TAR DNA binding protein 43 (TDP-43). Interestingly, ADAR2 overexpression in mammalian cell lines induced nuclear depletion and cytoplasmic accumulation of TDP-43, reflective of the pathology observed in ALS patients. The mislocalization of TDP-43 was dependent on the catalytic activity of ADAR2 and the ability of TDP-43 to bind directly to inosine containing RNA. In addition, TDP-43 nuclear export was significantly elevated in cells with increased RNA editing. Together these results describe a novel cellular mechanism by which alterations in RNA editing drive the nuclear export of TDP-43 leading to its cytoplasmic mislocalization. Considering the contribution of cytoplasmic TDP-43 to the pathogenesis of ALS, these findings represent a novel understanding of how the formation of pathogenic cytoplasmic TDP-43 accumulations may be initiated. Further research exploring this mechanism will provide insights into opportunities for novel therapeutic interventions.

The splicing of precursor messenger RNAs (pre-mRNAs) plays an essential role in dictating the mature mRNA profiles of eukaryotic cells. Mis-regulation of splicing, due to mutations in pre-mRNAs or in components of the splicing machinery, is associated with many diseases. Therefore, knowledge of pre-mRNA splicing mechanisms is required to understand gene expression regulation during states of homeostasis and disease, and for the development of therapeutic interventions.Splicing is catalyzed by the spliceosome, a dynamic and protein-rich ribozyme composed of five small nuclear ribonucleoproteins (snRNPs) and ~170 auxiliary factors. Early interactions that occur in prespliceosomal complexes formed by the 5′- and 3′-splice-site bound U1 and U2 snRNPs are responsible for committing introns for removal. However, the mechanisms underlying these early interactions remain to be fully characterized for understanding the influence of alternative splicing factors and the impact of recurrent disease-associated mutations in prespliceosomal proteins. The goal of my dissertation research was to delineate the role of the U1 small nuclear RNA (snRNA) during prespliceosome assembly. By applying a cellular minigene reporter assay and a variety of in vitro techniques including cell-free protein expression, UV-crosslinking, electrophoretic mobility shift assays, surface plasmon resonance, and RNA affinity purification, my work establishes critical roles for the U1 snRNA stem-loops 3 (SL3) and 4 (SL4) in formation of intron definition interactions during prespliceosome assembly. Previously, the SL4 of the U1 snRNA was shown to form a molecular bridge across introns by contacting the U2-specific splicing factor 3A1 (SF3A1). I identified the Ubiquitin-like domain of SF3A1 as a non-canonical RNA binding domain responsible for U1-SL4 binding. I also determined a role for the SL3 region of the U1 snRNA in splicing and characterized the spliceosomal RNA helicase UAP56 as an SL3 interacting protein. By knocking-down the SL3- and SL4-interacting proteins, I confirmed that U1 splicing activity in vivo relies on UAP56 and SF3A1 and that their functions are interdependent. These findings, in addition to the observations made using in vitro splicing assays, support a model whereby UAP56, through its interaction with U1-SL3, enhances the cross-intron interaction between U1-SL4 and SF3A1 to promote prespliceosome formation.

Schizophrenia, a debilitating neuropsychiatric disorder, affects 1% of the population. This multifaceted disorder is comprised of positive (hallucinations/psychosis), negative (social withdrawal/anhedonia) and cognitive symptoms. While treatments for schizophrenia have advanced over the past few years, high economic burdens are still conferred to society, totaling more than $34 billion in direct annual costs to the United States of America. Thus, a critical need exists to identify the factors that contribute towards the etiology of schizophrenia. This research aimed to determine the interactions between environmental factors and genetics in the etiology of schizophrenia. Specifically, this research shows that the immediate early gene, early growth response 3 (EGR3), which is upregulated in response to neuronal activity, resides at the center of a biological pathway to confer risk for schizophrenia. While schizophrenia-risk proteins including neuregulin 1 (NRG1) and N-methyl-D-aspartate receptors (NMDAR’s) have been identified upstream of EGR3, the downstream targets of EGR3 remain relatively unknown. This research demonstrates that early growth response 3 regulates the expression of the serotonin 2A-receptor (5HT2AR) in the frontal cortex following the physiologic stimulus, sleep deprivation. This effect is translated to the level of protein as 8 hours of sleep-deprivation results in the upregulation of 5HT2ARs, a target of antipsychotic medications. Additional downstream targets were identified following maximal upregulation of EGR3 through electroconvulsive stimulation (ECS). Both brain-derived neurotrophic factor (BDNF) and its epigenetic regulator, growth arrest DNA-damage-inducible 45 beta (GADD45B) are upregulated one-hour following ECS in the hippocampus and require the presence of EGR3. These proteins play important roles in both cellular proliferation and dendritic structural changes. Next, the effects of ECS on downstream neurobiological processes, hippocampal cellular proliferation and dendritic structural changes were examined. Following ECS, hippocampal cellular proliferationwas increased, and dendritic structural changes were observed in both wild-type and early growth response 3 knock-out (Egr3-/-) mice. Effects in the number of dendritic spines and dendritic complexity following ECS were not found to require EGR3. Collectively, these results demonstrate that neuronal activity leads to the regulation of schizophrenia risk proteins by EGR3 and point to a possible molecular mechanism contributing risk for schizophrenia.

Duchenne muscular dystrophy (DMD) is a lethal, X-linked disease characterized by progressive muscle degeneration. The condition is driven by out-of-frame mutations in the dystrophin gene, and the absence of a functional dystrophin protein ultimately leads to instability of the sarcolemma, skeletal muscle necrosis, and atrophy. While the structural changes that occur in dystrophic muscle are well characterized, resulting changes in muscle-specific gene expression that take place in dystrophin’s absence remain largely uncharacterized, as they are potentially obscured by the characteristic chronic inflammation in dystrophin deficient muscle.

The conservation of the dystrophin gene across metazoans suggests that both vertebrate and invertebrate model systems can provide valuable contributions to the understanding of DMD initiation and progression. Specifically, the invertebrate C. elegans possesses a dystrophin protein ortholog, dys-1, and a mild inflammatory response that is inactive in the muscle, allowing for the characterization of transcriptome rearrangements affecting disease progression independently of inflammation. Furthermore, C. elegans do not possess a satellite cell equivalent, meaning muscle regeneration does not occur. This makes C. elegans unique in that they allow for the study of dystrophin deficiencies without muscle regeneration that may obscure detection of subtle but consequential changes in gene expression.

I hypothesize that gaining a comprehensive definition of both the structural and signaling roles of dystrophin in C. elegans will improve the community’s understanding of the progression of DMD as a whole. To address this hypothesis, I have performed a phylogenetic analysis on the conservation of each member of the dystrophin associated protein complex (DAPC) across 10 species, established an in vivo system to identify muscle-specific changes in gene expression in the dystrophin-deficient C. elegans, and performed a functional analysis to test the biological significance of changes in gene expression identified in my sequencing results. The results from this study indicate that in C. elegans, dystrophin may have a signaling role early in development, and its absence may activate compensatory mechanisms that counteract disease progression. Furthermore, these findings allow for the identification of transcriptome changes that potentially serve as both independent drivers of disease and potential therapeutic targets for the treatment of DMD.

The successful reduction of CO2 and protons by a light-induced cobalt porphyrin/cytb562 hybrid metalloenzyme in water is reported. Incorporation of the porphyrin into a protein scaffold results in increases in CO and H2 production over naked porphyrin. Rational point mutations to the CoPPIX binding site of cytb562 modulate production, indicating possible further improvements in catalytic activity.