David T. Humphreys

Widespread occurrence of 5-methylcytosine in human coding and non-coding RNA

The modified base 5-methylcytosine (m5C) is well studied in DNA, but investigations of its prevalence in cellular RNA have been largely confined to tRNA and rRNA. In animals, the two m5C methyltransferases NSUN2 and TRDMT1 are known to modify specific tRNAs and have roles in the control of cell growth and differentiation. To map modified cytosine sites across a human transcriptome, we coupled bisulfite conversion of cellular RNA with next-generation sequencing. We confirmed 21 of the 28 previously known m5C sites in human tRNAs and identified 234 novel tRNA candidate sites, mostly in anticipated structural positions. Surprisingly, we discovered 10 275 sites in mRNAs and other non-coding RNAs. We observed that distribution of modified cytosines between RNA types was not random; within mRNAs they were enriched in the untranslated regions and near Argonaute binding regions. We also identified five new sites modified by NSUN2, broadening its known substrate range to another tRNA, the RPPH1 subunit of RNase P and two mRNAs. Our data demonstrates the widespread presence of modified cytosines throughout coding and non-coding sequences in a transcriptome, suggesting a broader role of this modification in the post-transcriptional control of cellular RNA function.

Neuronal MicroRNA Deregulation in Response to Alzheimer's Disease Amyloid-β

Normal brain development and function depends on microRNA (miRNA) networks to fine tune the balance between the transcriptome and proteome of the cell. These small non-coding RNA regulators are highly enriched in brain where they play key roles in neuronal development, plasticity and disease. In neurodegenerative disorders such as Alzheimers disease (AD), brain miRNA profiles are altered; thus miRNA dysfunction could be both a cause and a consequence of disease. Our study dissects the complexity of human AD pathology, and addresses the hypothesis that amyloid-β (Aβ) itself, a known causative factor of AD, causes neuronal miRNA deregulation, which could contribute to the pathomechanisms of AD. We used sensitive TaqMan low density miRNA arrays (TLDA) on murine primary hippocampal cultures to show that about half of all miRNAs tested were down-regulated in response to Aβ peptides. Time-course assays of neuronal Aβ treatments show that Aβ is in fact a powerful regulator of miRNA levels as the response of certain mature miRNAs is extremely rapid. Bioinformatic analysis predicts that the deregulated miRNAs are likely to affect target genes present in prominent neuronal pathways known to be disrupted in AD. Remarkably, we also found that the miRNA deregulation in hippocampal cultures was paralleled in vivo by a deregulation in the hippocampus of Aβ42-depositing APP23 mice, at the onset of Aβ plaque formation. In addition, the miRNA deregulation in hippocampal cultures and APP23 hippocampus overlaps with those obtained in human AD studies. Taken together, our findings suggest that neuronal miRNA deregulation in response to an insult by Aβ may be an important factor contributing to the cascade of events leading to AD.

Glioma microvesicles carry selectively packaged coding and non-coding RNAs which alter gene expression in recipient cells.

Interactions between glioma cells and their local environment are critical determinants of brain tumor growth, infiltration and neovascularisation. Communication with host cells and stroma via microvesicles represents one pathway by which tumors can modify their surroundings to achieve a tumor-permissive environment. Here we have taken an unbiased approach to identifying RNAs in glioma-derived microvesicles, and explored their potential to regulate gene expression in recipient cells. We find that glioma microvesicles are predominantly of exosomal origin and contain complex populations of coding and noncoding RNAs in proportions that are distinct from those in the cells from which they are derived. Microvesicles show a relative depletion in microRNA compared with their cells of origin, and are enriched in unusual or novel noncoding RNAs, most of which have no known function. Short-term exposure of brain microvascular endothelial cells to glioma microvesicles results in many gene expression changes in the endothelial cells, most of which cannot be explained by direct delivery of transcripts. Our data suggest that the scope of potential actions of tumor-derived microvesicles is much broader and more complex than previously supposed, and highlight a number of new classes of small RNA that remain to be characterized.

Target Gene Repression Mediated by miRNAs miR-181c and miR-9 Both of Which Are Down-regulated by Amyloid-β

MicroRNAs (miRNAs) are small non-coding RNA regulators of protein synthesis that are essential for normal brain development and function. Their profiles are significantly altered in neurodegenerative diseases such as Alzheimer’s disease (AD) that is characterized by amyloid-β (Aβ) and tau deposition in brain. How deregulated miRNAs contribute to AD is not understood, as their dysfunction could be both a cause and a consequence of disease. To address this question we had previously profiled miRNAs in models of AD. This identified miR-9 and -181c as being down-regulated by Aβ in hippocampal cultures. Interestingly, there was a remarkable overlap with those miRNAs that are deregulated in Aβ-depositing APP23 transgenic mice and in human AD tissue. While the Aβ precursor protein APP itself is a target of miRNA regulation, the challenge resides in identifying further targets. Here, we expand the repertoire of miRNA target genes by identifying the 3′ untranslated regions (3′ UTRs) of TGFBI, TRIM2, SIRT1 and BTBD3 as being repressed by miR-9 and -181c, either alone or in combination. Taken together, our study identifies putative target genes of miRNAs miR-9 and 181c, which may function in brain homeostasis and disease pathogenesis.

microRNA-Mediated Messenger RNA Deadenylation Contributes to Translational Repression in Mammalian Cells

Animal microRNAs (miRNAs) typically regulate gene expression by binding to partially complementary target sites in the 3′ untranslated region (UTR) of messenger RNA (mRNA) reducing its translation and stability. They also commonly induce shortening of the mRNA 3′ poly(A) tail, which contributes to their mRNA decay promoting function. The relationship between miRNA-mediated deadenylation and translational repression has been less clear. Using transfection of reporter constructs carrying three imperfectly matching let-7 target sites in the 3′ UTR into mammalian cells we observe rapid target mRNA deadenylation that precedes measureable translational repression by endogenous let-7 miRNA. Depleting cells of the argonaute co-factors RCK or TNRC6A can impair let-7-mediated repression despite ongoing mRNA deadenylation, indicating that deadenylation alone is not sufficient to effect full repression. Nevertheless, the magnitude of translational repression by let-7 is diminished when the target reporter lacks a poly(A) tail. Employing an antisense strategy to block deadenylation of target mRNA with poly(A) tail also partially impairs translational repression. On the one hand, these experiments confirm that tail removal by deadenylation is not strictly required for translational repression. On the other hand they show directly that deadenylation can augment miRNA-mediated translational repression in mammalian cells beyond stimulating mRNA decay. Taken together with published work, these results suggest a dual role of deadenylation in miRNA function: it contributes to translational repression as well as mRNA decay and is thus critically involved in establishing the quantitatively appropriate physiological response to miRNAs.

Complexity of Murine Cardiomyocyte miRNA Biogenesis, Sequence Variant Expression and Function

microRNAs (miRNAs) are critical to heart development and disease. Emerging research indicates that regulated precursor processing can give rise to an unexpected diversity of miRNA variants. We subjected small RNA from murine HL-1 cardiomyocyte cells to next generation sequencing to investigate the relevance of such diversity to cardiac biology. ∼40 million tags were mapped to known miRNA hairpin sequences as deposited in miRBase version 16, calling 403 generic miRNAs as appreciably expressed. Hairpin arm bias broadly agreed with miRBase annotation, although 44 miR* were unexpectedly abundant (>20% of tags); conversely, 33 -5p/-3p annotated hairpins were asymmetrically expressed. Overall, variability was infrequent at the 5′ start but common at the 3′ end of miRNAs (5.2% and 52.3% of tags, respectively). Nevertheless, 105 miRNAs showed marked 5′ isomiR expression (>20% of tags). Among these was miR-133a, a miRNA with important cardiac functions, and we demonstrated differential mRNA targeting by two of its prevalent 5′ isomiRs. Analyses of miRNA termini and base-pairing patterns around Drosha and Dicer cleavage regions confirmed the known bias towards uridine at the 5′ most position of miRNAs, as well as supporting the thermodynamic asymmetry rule for miRNA strand selection and a role for local structural distortions in fine tuning miRNA processing. We further recorded appreciable expression of 5 novel miR*, 38 extreme variants and 8 antisense miRNAs. Analysis of genome-mapped tags revealed 147 novel candidate miRNAs. In summary, we revealed pronounced sequence diversity among cardiomyocyte miRNAs, knowledge of which will underpin future research into the mechanisms involved in miRNA biogenesis and, importantly, cardiac function, disease and therapy.

The human Piwi protein Hiwi2 associates with tRNA-derived piRNAs in somatic cells

The Piwi-piRNA pathway is active in animal germ cells where its functions are required for germ cell maintenance and gamete differentiation. Piwi proteins and piRNAs have been detected outside germline tissue in multiple phyla, but activity of the pathway in mammalian somatic cells has been little explored. In particular, Piwi expression has been observed in cancer cells, but nothing is known about the piRNA partners or the function of the system in these cells. We have surveyed the expression of the three human Piwi genes, Hiwi, Hili and Hiwi2, in multiple normal tissues and cancer cell lines. We find that Hiwi2 is ubiquitously expressed; in cancer cells the protein is largely restricted to the cytoplasm and is associated with translating ribosomes. Immunoprecipitation of Hiwi2 from MDAMB231 cancer cells enriches for piRNAs that are predominantly derived from processed tRNAs and expressed genes, species which can also be found in adult human testis. Our studies indicate that a Piwi-piRNA pathway is present in human somatic cells, with an uncharacterised function linked to translation. Taking this evidence together with evidence from primitive organisms, we propose that this somatic function of the pathway predates the germline functions of the pathway in modern animals.

Methods to analyze microRNA-mediated control of mRNA translation.

MicroRNAs (miRs) are an important class of gene regulators that affect a wide range of biological processes. Despite the early recognition of miRs as translational regulators and intense interest in studying this phenomenon, it has so far not been possible to derive a consensus model for the underlying molecular mechanism(s). The potential of miRs to act in a combinatorial manner and to also promote mRNA decay creates conceptual and technical challenges in their study. Here, we discuss critical parameters in design and analysis of experiments used to study miR function including creation of synthetic miR and mRNA partners for assay of translational inhibition using luciferase reporters; measurement of mRNA stability after miR action; defining poly(A) tail length in miR target mRNA; determining the distribution of miRs and their target mRNAs in polysome profiles; and visualization of P-body components. We describe protocols for each of these procedures.

Male-lineage transmission of an acquired metabolic phenotype induced by grand-paternal obesity

Objective Parental obesity can induce metabolic phenotypes in offspring independent of the inherited DNA sequence. Here we asked whether such non-genetic acquired metabolic traits can be passed on to a second generation that has never been exposed to obesity, even as germ cells. Methods We examined the F1, F2, and F3 a/a offspring derived from F0 matings of obese prediabetic Avy/a sires and lean a/a dams. After F0, only lean a/a mice were used for breeding. Results We found that F1 sons of obese founder males exhibited defects in glucose and lipid metabolism, but only upon a post-weaning dietary challenge. F1 males transmitted these defects to their own male progeny (F2) in the absence of the dietary challenge, but the phenotype was largely attenuated by F3. The sperm of F1 males exhibited changes in the abundance of several small RNA species, including the recently reported diet-responsive tRNA-derived fragments. Conclusions These data indicate that induced metabolic phenotypes may be propagated for a generation beyond any direct exposure to an inducing factor. This non-genetic inheritance likely occurs via the actions of sperm noncoding RNA.

Genetic variation in the two-pore domain potassium channel, TASK-1, may contribute to an atrial substrate for arrhythmogenesis

The two-pore domain potassium channel, K2P3.1 (TASK-1) modulates background conductance in isolated human atrial cardiomyocytes and has been proposed as a potential drug target for atrial fibrillation (AF). TASK-1 knockout mice have a predominantly ventricular phenotype however, and effects of TASK-1 inactivation on atrial structure and function have yet to be demonstrated in vivo. The extent to which genetic variation in KCNK3, that encodes TASK-1, might be a determinant of susceptibility to AF is also unknown. To address these questions, we first evaluated the effects of transient knockdown of the zebrafish kcnk3a and kcnk3b genes and cardiac phenotypes were evaluated using videomicroscopy. Combined kcnk3a and kcnk3b knockdown in 72 hour post fertilization embryos resulted in lower heart rate (p<0.001), marked increase in atrial diameter (p<0.001), and mild increase in end-diastolic ventricular diameter (p=0.01) when compared with control-injected embryos. We next performed genetic screening of KCNK3 in two independent AF cohorts (373 subjects) and identified three novel KCNK3 variants. Two of these variants, present in one proband with familial AF, were located at adjacent nucleotides in the Kozak sequence and reduced expression of an engineered reporter. A third missense variant, V123L, in a patient with lone AF, reduced resting membrane potential and altered pH sensitivity in patch-clamp experiments, with structural modeling predicting instability in the vicinity of the TASK-1 pore. These in vitro data suggest that the double Kozak variants and V123L will have loss-of-function effects on ITASK. Cardiac action potential modeling predicted that reduced ITASK prolongs atrial action potential duration, and that this is potentiated by reciprocal changes in activity of other ion channel currents. Our findings demonstrate the functional importance of ITASK in the atrium and suggest that inactivation of TASK-1 may have diverse effects on atrial size and electrophysiological properties that can contribute to an arrhythmogenic substrate.