The Role of TET Proteins in the Epigenetic Regulation of Neural Gene Expression and Behavior

Understanding how genes affect behavior is critical to develop precise therapies for human behavioral disorders. The ability to investigate the relationship between genes and behavior has been greatly advanced over the last few decades due to progress in gene-targeting technology. Recently, the Tet gene family was discovered and implicated in epigenetic modification of DNA methylation by converting 5-methylcytosine to 5-hydroxymethylcytosine (5hmC). 5hmC and its catalysts, the TET proteins, are highly abundant in the postnatal brain but with unclear functions. To investigate their neural functions, we generated new lines of Tet1 and Tet3 mutant mice using a gene targeting approach. We designed both mutations to cause a frameshift by deleting the largest coding exon of Tet1 (Tet1Δe4) and the catalytic domain of Tet3 (Tet3Δe7-9). As Tet1 is also highly expressed in embryonic stem cells (ESCs), we generated Tet1 homozygous deleted ESCs through sequential targeting to compare the function of Tet1 in the brain to its role in ESCs. To test our hypothesis that TET proteins epigenetically regulate transcription of key neural genes important for normal brain function, we examined transcriptional and epigenetic differences in the Tet1Δe4 mouse brain. The oxytocin receptor (OXTR), a neural gene implicated in social behaviors, is suggested to be epigenetically regulated by an unknown mechanism. Interestingly, several human studies have found associations between OXTR DNA hypermethylation and a wide spectrum of behavioral traits and neuropsychiatric disorders including autism spectrum disorders. Here we report the first evidence for an epigenetic mechanism of Oxtr transcription as expression of Oxtr is reduced in the brains of Tet1Δe4-/- mice. Likewise, the CpG island overlapping the promoter of Oxtr is hypermethylated during early embryonic development and persists into adulthood. We also discovered altered histone modifications at the hypermethylated regions, indicating the loss of TET1 has broad effects on the chromatin structure at Oxtr. Unexpectedly, we discovered an array of novel mRNA isoforms of Oxtr that are selectively reduced in Tet1Δe4-/- mice. Additionally, Tet1Δe4-/- mice display increased agonistic behaviors and impaired maternal care and short-term memory. Our findings support a novel role for TET1 in regulating Oxtr expression by preventing DNA hypermethylation and implicate TET1 in social behaviors, offering novel insight into Oxtr epigenetic regulation and its role in neuropsychiatric disorders.

Phocid seal leptin: tertiary structure and hydrophobic receptor binding site preservation during distinct leptin gene evolution

The cytokine hormone leptin is a key signalling molecule in many pathways that control physiological functions. Although leptin demonstrates structural conservation in mammals, there is evidence of positive selection in primates, lagomorphs and chiropterans. We previously reported that the leptin genes of the grey and harbour seals (phocids) have significantly diverged from other mammals. Therefore we further investigated the diversification of leptin in phocids, other marine mammals and terrestrial taxa by sequencing the leptin genes of representative species. Phylogenetic reconstruction revealed that leptin diversification was pronounced within the phocid seals with a high dN/dS ratio of 2.8, indicating positive selection. We found significant evidence of positive selection along the branch leading to the phocids, within the phocid clade, but not over the dataset as a whole. Structural predictions indicate that the individual residues under selection are away from the leptin receptor (LEPR) binding site. Predictions of the surface electrostatic potential indicate that phocid seal leptin is notably different to other mammalian leptins, including the otariids. Cloning the grey seal leptin binding domain of LEPR confirmed that this was structurally conserved. These data, viewed in toto, support a hypothesis that phocid leptin divergence is unlikely to have arisen by random mutation. Based upon these phylogenetic and structural assessments, and considering the comparative physiology and varying life histories among species, we postulate that the unique phocid diving behaviour has produced this selection pressure. The Phocidae includes some of the deepest diving species, yet have the least modified lung structure to cope with pressure and volume changes experienced at depth. Therefore, greater surfactant production is required to facilitate rapid lung re-inflation upon surfacing, while maintaining patent airways. We suggest that this additional surfactant requirement is met by the leptin pulmonary surfactant production pathway which normally appears only to function in the mammalian foetus.

Functional analysis of the 5′ flanking region of the human G6PC3 gene: regulation of promoter activity by glucose, pyruvate, AMP kinase and the pentose phosphate pathway

G6PC3 is a widely expressed isoform of glucose-6-phosphatase, found in many foetal and adult tissues. Mutations in this gene cause developmental abnormalities and severe neutropenia due to abolition of glucose recycling between the cytoplasm and endoplasmic reticulum. Low G6PC3 expression as a result of promoter polymorphisms or dysregulation could produce similar outcomes. Here we investigated the regulation of human G6PC3 promoter activity. HeLa and H4IIE cells were transiently transfected with G6PC3 promoter coupled to the firefly luciferase gene, and promoter activity was measured by dual luciferase assay. Activity was highest in a 453 bp segment of the G6PC3 promoter, from − 455 to − 3 relative to the transcriptional start site. This promoter was unresponsive to glucostatic hormones. Its activity increased significantly between 1 and 5.5 mM glucose, and was not elevated further by glucose concentrations up to 25 mM. Pyruvate increased its activity, but β-hydroxybutyrate and sodium acetate did not. Promoter activity was reduced by inhibitors of hexokinase, glyceraldehyde phosphate dehydrogenase and the oxidative branch of the pentose phosphate pathway, but not by a transketolase inhibitor. Deletion of two adjacent Enhancer-boxes (− 274 to − 279 and − 299 to − 304) reduced promoter activity and abolished the glucose effect, suggesting they could function as a glucose response element. Deletion of an additional downstream 140 bp (− 140 to − 306) restored activity, but not the glucose response, suggesting the presence of repressor elements in this region. 5-Aminoimidazole-4-carboxamide 1-β-d-ribofuranoside (AICAR) reduced promoter activity, showing dependence on AMP-kinase. Regulation of the G6PC3 promoter is thus radically different to that of the hepatic isoform, G6PC. It is sensitive to carbohydrate, but not to fatty acid metabolites, and at much lower physiological concentrations. Based on these findings, we speculate that reduced G6PC3 expression could occur during hypoglycemic episodes in vivo, which are common in utero and in the postnatal period. If such episodes lower G6PC3 expression they could place the foetus or infant at risk of impaired immune function and development, and this possibility requires further examination both in vitro and in vivo.

2C-Cas9: a versatile tool for clonal analysis of gene function.

CRISPR/Cas9-mediated targeted mutagenesis allows efficient generation of loss-of-function alleles in zebrafish. To date, this technology has been primarily used to generate genetic knockout animals. Nevertheless, the study of the function of certain loci might require tight spatiotemporal control of gene inactivation. Here, we show that tissue-specific gene disruption can be achieved by driving Cas9 expression with the Gal4/UAS system. Furthermore, by combining the Gal4/UAS and Cre/loxP systems, we establish a versatile tool to genetically label mutant cell clones, enabling their phenotypic analysis. Our technique has the potential to be applied to diverse model organisms, enabling tissue-specific loss-of-function and phenotypic characterization of live and fixed tissues.

Synthetic regulation of eukaryotic gene expression by noncoding RNA

Synthetic biological systems promise to combine the spectacular diversity of biological functionality with engineering principles to design new life to address many pressing needs. As these engineered systems advance in sophistication, there is ever-greater need for customizable, situation-specific expression of desired genes. However, existing gene control platforms are generally not modular, or do not display performance requirements required for robust phenotypic responses to input signals. This work expands the capabilities of eukaryotic gene control in two important directions.

For development of greater modularity, we extend the use of synthetic self-cleaving ribozyme switches to detect changes in input protein levels and convey that information into programmed gene expression in eukaryotic cells. We demonstrate both up- and down-regulation of levels of an output transgene by more than 4-fold in response to rising input protein levels, with maximal output gene expression approaching the highest levels observed in yeast. In vitro experiments demonstrate protein-dependent ribozyme activity modulation. We further demonstrate the platform in mammalian cells. Our switch devices do not depend on special input protein activity, and can be tailored to respond to any input protein to which a suitable RNA aptamer can be developed. This platform can potentially be employed to regulate the expression of any transgene or any endogenous gene by 3’ UTR replacement, allowing for more complex cell state-specific reprogramming.

We also address an important concern with ribozyme switches, and riboswitch performance in general, their dynamic range. While riboswitches have generally allowed for versatile and modular regulation, so far their dynamic ranges of output gene modulation have been modest, generally at most 10-fold. We address this shortcoming by developing a modular genetic amplifier for near-digital control of eukaryotic gene expression. We combine ribozyme switch-mediated regulation of a synthetic TF with TF-mediated regulation of an output gene. The amplifier platform allows for as much as 20-fold regulation of output gene expression in response to input signal, with maximal expression approaching the highest levels observed in yeast, yet being tunable to intermediate and lower expression levels. EC50 values are more than 4 times lower than in previously best-performing non-amplifier ribozyme switches. The system design retains the modular-input architecture of the ribozyme switch platform, and the near-digital dynamic ranges of TF-based gene control.

Together, these developments suggest great potential for the wide applicability of these platforms for better-performing eukaryotic gene regulation, and more sophisticated, customizable reprogramming of cellular activity.

Evolution of developmental gene regulatory networks in echinoids

Developmental gene regulatory networks (dGRNs) are assemblages of regulatory genes that direct embryonic development of animal body plans and their morpho-logical structures. dGRNs exhibit recursively-wired circuitry that is encoded in the genome and executed during development. Alteration to the regulatory architecture of dGRNs causes variation in developmental programs both during the development of an individual organism and during the evolution of an individual lineage. The ex-planatory power of these networks is best exemplified by the global dGRN directing early development of the euechinoid sea urchin Strongylocentrotus purpuratus. This network consists of numerous regulatory genes engaging in hundreds of genomic regulatory transactions that collectively direct the delineation of early embryonic domains and the specification of cell lineages. Research on closely-related euechi-noid sea urchins, e.g. Lytechinus variegatus and Paracentrotus lividus, has revealed marked conservation of dGRN architecture in echinoid development, suggesting little appreciable alteration has occurred since their divergence in evolution at least 90 million years ago (mya).

We sought to test whether this observation extends to all sea urchins (echinoids) and undertook a systematic analysis of over 50 regulatory genes in the cidaroid sea urchin Eucidaris tribuloides, surveing their regulatory activity and function in a sea urchin that diverged from euechinoid sea urchins at least 268 mya. Our results revealed extensive alterations have occurred to all levels of echinoid dGRN archi-tecture since the cidaroid-euechinoid divergence. Alterations to mesodermal sub-circuits were particularly striking, including functional di˙erences in specification of non-skeletogenic mesenchyme (NSM), skeletogenic mesenchyme (SM), and en-domesodermal segregation. Specification of endomesodermal embryonic domains revealed that, while their underlying network circuitry had clearly diverged, regu-latory states established in pregastrular embryos of these two groups are strikingly similar. Analyses of E. tribuloides specification leading to the estab-lishment of dorsal-ventral (aboral-oral) larval polarity indicated that regulation of regulatory genes expressed in mesodermal embryonic domains had incurred significantly more alterations than those expressed in endodermal and ectodermal domains. Taken together, this study highlights the ability of dGRN architecture to buffer extensive alterations in the evolution and early development of echinoids and adds further support to the notion that alterations can occur at all levels of dGRN architecture and all stages of embryonic development.

Hormone-Mediated Gene-Specific Translation Regulation

The central role of translation regulation in the control of critical cellular processes has long been recognized. Yet the systematic exploration of quantitative changes in translation at a genome-wide scale in response to specific stimuli has only recently become technically feasible. Using a genetic approach, we have identified new Arabidopsis weak-ethylene insensitive mutants that also display defects in translation, which suggested the existence of a previously unknown molecular module involved in ethylene-mediated translation regulation of components of this signaling pathway. To explore this link in detail, we implemented for Arabidopsis the ribosome-footprinting technology, which enables the study of translation at a whole-genome level at single codon resolution[1]. Using ribosome-footprinting we examined the effects of short exposure to ethylene on the Arabidopsis translatome looking for ethylene-triggered changes in translation rates that could not be explained by changes in transcript levels. The results of this research, in combination with the characterization of a subset of the aforementioned weak-ethylene insensitive mutants that are defective in the UPF genes (core-components of the nonsense-mediated mRNA decay machinery), uncovered a translation-based branch of the ethylene signaling pathway[2]. In the presence of ethylene, translation of a negative regulator of ethylene signaling EBF2 is repressed, despite induced transcription of this gene. These translational effects of ethylene require the long 3´UTR of EBF2 (3´EBF2), which is recognized by the C-terminal end of the key ethylene-signaling protein EIN2 (EIN2C) in the cytoplasm once EIN2C is released from the ER-membrane by proteolytic cleavage. EIN2C binds the 3´EBF2, recruits the UPF proteins and moves to P-bodies, where the translation of EBF2 in inhibited despite its mRNA accumulation. Once the ethylene signal is withdrawn, the translation of the stored EBF2 mRNAs is resumed, thus rapidly dampening the ethylene response. These findings represent a mechanistic paradigm of gene-specific regulation of translation in response to a key growth regulator. Translation regulatory elements can be located in both 3′ and 5′ UTRs. We are now focusing on the ead1 and ead2 mutants, another set of ethylene-signaling mutants defective in translational regulation. Ribosome-footprinting on the ead1 mutant revealed an accumulation of translating ribosomes in the 5´UTRs of uORF-containing genes and reduction in the levels of ribosomes in the main ORF. The mutant is also impaired in the translation of GFP when this reporter is fused to WT 5´UTR of potential EAD1 targets but not when GFP is fused to the uORF-less versions of the same 5´UTRs. Our hypothesis is that EAD1/2 work as a complex that is required for the efficient translation of mRNAs that have common structural (complex 5´UTR with uORFs) and functional (regulation of key cellular processes) features. We are working towards the identification of the conditions where the EAD1 regulation of translation is required. [1] Ingolia, N. et al. (2009) Genome-Wide Analysis in Vivo of Translation with Nucleotide Resolution Using Ribosome Profiling. Science, 324; 218-222 [2] Merchante, C. et al. (2015) Gene-Specific Translation Regulation Mediated by the Hormone-Signaling Molecule EIN2. Cell, 163(3): 684-697

POPLAR GENE EXPRESSION DATA ANALYSIS PIPELINES

Analyzing large-scale gene expression data is a labor-intensive and time-consuming process. To make data analysis easier, we developed a set of pipelines for rapid processing and analysis poplar gene expression data for knowledge discovery. Of all pipelines developed, differentially expressed genes (DEGs) pipeline is the one designed to identify biologically important genes that are differentially expressed in one of multiple time points for conditions. Pathway analysis pipeline was designed to identify the differentially expression metabolic pathways. Protein domain enrichment pipeline can identify the enriched protein domains present in the DEGs. Finally, Gene Ontology (GO) enrichment analysis pipeline was developed to identify the enriched GO terms in the DEGs. Our pipeline tools can analyze both microarray gene data and high-throughput gene data. These two types of data are obtained by two different technologies. A microarray technology is to measure gene expression levels via microarray chips, a collection of microscopic DNA spots attached to a solid (glass) surface, whereas high throughput sequencing, also called as the next-generation sequencing, is a new technology to measure gene expression levels by directly sequencing mRNAs, and obtaining each mRNA’s copy numbers in cells or tissues. We also developed a web portal (http://sys.bio.mtu.edu/) to make all pipelines available to public to facilitate users to analyze their gene expression data. In addition to the analyses mentioned above, it can also perform GO hierarchy analysis, i.e. construct GO trees using a list of GO terms as an input.

Protein expression and gene polymorphism of CXCL10 in patients with colorectal cancer

Chemokines (chemotactic cytokines) promote leukocyte attraction to sites of inflammation and cancer. Certain chemokines promote and regulate neoplastic progression, including metastasis and angiogenesis. One such chemokine, CXCL10, was found to be expressed in colorectal cancer (CRC) tissue. To gain insight into the prognostic significance of CXCL10, we investigated whether the levels of this chemokine were altered in the colorectal tissue or plasma of CRC patients. Using Luminex technology for protein analyses, we observed a significantly higher CXCL10 protein level in cancer tissue compared to that in paired normal tissue. Moreover, significantly higher plasma levels of CXCL10 were detected in patients compared to those in control subjects and the plasma levels of CXCL10 in disseminated disease were found to be significantly higher compared to those in localized disease. The single‑nucleotide polymorphism rs8878, which has been described in exon 4 in the 3'‑untranslated region of the CXCL10 gene, was investigated using a TaqMan system. There were significant differences in genotype distribution and allelic frequencies between CRC patients and control subjects. In conclusion, altered CXCL10 protein concentrations in CRC tissues or plasma and the rs8878 genotype variant of CXCL10 may contribute to the prediction of clinical outcome.

Raman microscopy of formamide-intercalated kaolinites treated by controlled-rate thermal analysis technology

Raman spectroscopy of formamide-intercalated kaolinites treated using controlled-rate thermal analysis technology (CRTA), allowing the separation of adsorbed formamide from intercalated formamide in formamide-intercalated kaolinites, is reported. The Raman spectra of the CRTA-treated formamide-intercalated kaolinites are significantly different from those of the intercalated kaolinites, which display a combination of both intercalated and adsorbed formamide. An intense band is observed at 3629 cm-1, attributed to the inner surface hydroxyls hydrogen bonded to the formamide. Broad bands are observed at 3600 and 3639 cm-1, assigned to the inner surface hydroxyls, which are hydrogen bonded to the adsorbed water molecules. The hydroxyl-stretching band of the inner hydroxyl is observed at 3621 cm-1 in the Raman spectra of the CRTA-treated formamide-intercalated kaolinites. The results of thermal analysis show that the amount of intercalated formamide between the kaolinite layers is independent of the presence of water. Significant differences are observed in the CO stretching region between the adsorbed and intercalated formamide.

No evidence for association of ataxia-telangiectasia mutated gene T2119C and C3161G amino acid substitution variants with risk of breast cancer

Background There is evidence that certain mutations in the double-strand break repair pathway ataxia-telangiectasia mutated gene act in a dominant-negative manner to increase the risk of breast cancer. There are also some reports to suggest that the amino acid substitution variants T2119C Ser707Pro and C3161G Pro1054Arg may be associated with breast cancer risk. We investigate the breast cancer risk associated with these two nonconservative amino acid substitution variants using a large Australian population-based case–control study. Methods The polymorphisms were genotyped in more than 1300 cases and 600 controls using 5' exonuclease assays. Case–control analyses and genotype distributions were compared by logistic regression. Results The 2119C variant was rare, occurring at frequencies of 1.4 and 1.3% in cases and controls, respectively (P = 0.8). There was no difference in genotype distribution between cases and controls (P = 0.8), and the TC genotype was not associated with increased risk of breast cancer (adjusted odds ratio = 1.08, 95% confidence interval = 0.59–1.97, P = 0.8). Similarly, the 3161G variant was no more common in cases than in controls (2.9% versus 2.2%, P = 0.2), there was no difference in genotype distribution between cases and controls (P = 0.1), and the CG genotype was not associated with an increased risk of breast cancer (adjusted odds ratio = 1.30, 95% confidence interval = 0.85–1.98, P = 0.2). This lack of evidence for an association persisted within groups defined by the family history of breast cancer or by age. Conclusion The 2119C and 3161G amino acid substitution variants are not associated with moderate or high risks of breast cancer in Australian women.