Identification of microRNAs expressed in two mosquito vectors, Aedes albopictus and Culex quinquefasciatus.

BACKGROUND: MicroRNAs (miRNAs) are small non-coding RNAs that post-transcriptionally regulate gene expression in a variety of organisms, including insects, vertebrates, and plants. miRNAs play important roles in cell development and differentiation as well as in the cellular response to stress and infection. To date, there are limited reports of miRNA identification in mosquitoes, insects that act as essential vectors for the transmission of many human pathogens, including flaviviruses. West Nile virus (WNV) and dengue virus, members of the Flaviviridae family, are primarily transmitted by Aedes and Culex mosquitoes. Using high-throughput deep sequencing, we examined the miRNA repertoire in Ae. albopictus cells and Cx. quinquefasciatus mosquitoes. RESULTS: We identified a total of 65 miRNAs in the Ae. albopictus C7/10 cell line and 77 miRNAs in Cx. quinquefasciatus mosquitoes, the majority of which are conserved in other insects such as Drosophila melanogaster and Anopheles gambiae. The most highly expressed miRNA in both mosquito species was miR-184, a miRNA conserved from insects to vertebrates. Several previously reported Anopheles miRNAs, including miR-1890 and miR-1891, were also found in Culex and Aedes, and appear to be restricted to mosquitoes. We identified seven novel miRNAs, arising from nine different precursors, in C7/10 cells and Cx. quinquefasciatus mosquitoes, two of which have predicted orthologs in An. gambiae. Several of these novel miRNAs reside within a ~350 nt long cluster present in both Aedes and Culex. miRNA expression was confirmed by primer extension analysis. To determine whether flavivirus infection affects miRNA expression, we infected female Culex mosquitoes with WNV. Two miRNAs, miR-92 and miR-989, showed significant changes in expression levels following WNV infection. CONCLUSIONS: Aedes and Culex mosquitoes are important flavivirus vectors. Recent advances in both mosquito genomics and high-throughput sequencing technologies enabled us to interrogate the miRNA profile in these two species. Here, we provide evidence for over 60 conserved and seven novel mosquito miRNAs, expanding upon our current understanding of insect miRNAs. Undoubtedly, some of the miRNAs identified will have roles not only in mosquito development, but also in mediating viral infection in the mosquito host.

The genomic analysis of erythrocyte microRNA expression in sickle cell diseases.

BACKGROUND: Since mature erythrocytes are terminally differentiated cells without nuclei and organelles, it is commonly thought that they do not contain nucleic acids. In this study, we have re-examined this issue by analyzing the transcriptome of a purified population of human mature erythrocytes from individuals with normal hemoglobin (HbAA) and homozygous sickle cell disease (HbSS). METHODS AND FINDINGS: Using a combination of microarray analysis, real-time RT-PCR and Northern blots, we found that mature erythrocytes, while lacking ribosomal and large-sized RNAs, contain abundant and diverse microRNAs. MicroRNA expression of erythrocytes was different from that of reticulocytes and leukocytes, and contributed the majority of the microRNA expression in whole blood. When we used microRNA microarrays to analyze erythrocytes from HbAA and HbSS individuals, we noted a dramatic difference in their microRNA expression pattern. We found that miR-320 played an important role for the down-regulation of its target gene, CD71 during reticulocyte terminal differentiation. Further investigation revealed that poor expression of miR-320 in HbSS cells was associated with their defective downregulation CD71 during terminal differentiation. CONCLUSIONS: In summary, we have discovered significant microRNA expression in human mature erythrocytes, which is dramatically altered in HbSS erythrocytes and their defect in terminal differentiation. Thus, the global analysis of microRNA expression in circulating erythrocytes can provide mechanistic insights into the disease phenotypes of erythrocyte diseases.

Intergenic and genic sequence lengths have opposite relationships with respect to gene expression.

Eukaryotic genomes are mostly composed of noncoding DNA whose role is still poorly understood. Studies in several organisms have shown correlations between the length of the intergenic and genic sequences of a gene and the expression of its corresponding mRNA transcript. Some studies have found a positive relationship between intergenic sequence length and expression diversity between tissues, and concluded that genes under greater regulatory control require more regulatory information in their intergenic sequences. Other reports found a negative relationship between expression level and gene length and the interpretation was that there is selection pressure for highly expressed genes to remain small. However, a correlation between gene sequence length and expression diversity, opposite to that observed for intergenic sequences, has also been reported, and to date there is no testable explanation for this observation. To shed light on these varied and sometimes conflicting results, we performed a thorough study of the relationships between sequence length and gene expression using cell-type (tissue) specific microarray data in Arabidopsis thaliana. We measured median gene expression across tissues (expression level), expression variability between tissues (expression pattern uniformity), and expression variability between replicates (expression noise). We found that intergenic (upstream and downstream) and genic (coding and noncoding) sequences have generally opposite relationships with respect to expression, whether it is tissue variability, median, or expression noise. To explain these results we propose a model, in which the lengths of the intergenic and genic sequences have opposite effects on the ability of the transcribed region of the gene to be epigenetically regulated for differential expression. These findings could shed light on the role and influence of noncoding sequences on gene expression.

The prevalence and regulation of antisense transcripts in Schizosaccharomyces pombe.

A strand-specific transcriptome sequencing strategy, directional ligation sequencing or DeLi-seq, was employed to profile antisense transcriptome of Schizosaccharomyces pombe. Under both normal and heat shock conditions, we found that polyadenylated antisense transcripts are broadly expressed while distinct expression patterns were observed for protein-coding and non-coding loci. Dominant antisense expression is enriched in protein-coding genes involved in meiosis or stress response pathways. Detailed analyses further suggest that antisense transcripts are independently regulated with respect to their sense transcripts, and diverse mechanisms might be potentially involved in the biogenesis and degradation of antisense RNAs. Taken together, antisense transcription may have profound impacts on global gene regulation in S. pombe.

A long non-coding RNA targets microRNA miR-34a to regulate colon cancer stem cell asymmetric division.

The roles of long non-coding RNAs (lncRNAs) in regulating cancer and stem cells are being increasingly appreciated. Its diverse mechanisms provide the regulatory network with a bigger repertoire to increase complexity. Here we report a novel LncRNA, Lnc34a, that is enriched in colon cancer stem cells (CCSCs) and initiates asymmetric division by directly targeting the microRNA miR-34a to cause its spatial imbalance. Lnc34a recruits Dnmt3a via PHB2 and HDAC1 to methylate and deacetylate the miR-34a promoter simultaneously, hence epigenetically silencing miR-34a expression independent of its upstream regulator, p53. Lnc34a levels affect CCSC self-renewal and colorectal cancer (CRC) growth in xenograft models. Lnc34a is upregulated in late-stage CRCs, contributing to epigenetic miR-34a silencing and CRC proliferation. The fact that lncRNA targets microRNA highlights the regulatory complexity of non-coding RNAs (ncRNAs), which occupy the bulk of the genome.

Cell Fate Specification and the Regulation of RNA-dependent DNA Methylation in the Arabidopsis Root Meristem

The Arabidopsis root apical meristem (RAM) is a complex tissue capable of generating all the cell types that ultimately make up the root. The work presented in this thesis takes advantage of the versatility of high-throughput sequencing to address two independent questions about the root meristem. Although a lot of information is known regarding the cell fate decisions that occur at the RAM, cortex specification and differentiation remain poorly understood. In the first part of this thesis, I used an ethylmethanesulfonate (EMS) mutagenized marker line to perform a forward genetics screen. The goal of this screen was to identify novel genes involved in the specification and differentiation of the cortex tissue. Mapping analysis from the results obtained in this screen revealed a new allele of BRASSINOSTEROID4 with abnormal marker expression in the cortex tissue. Although this allele proved to be non-cortex specific, this project highlights new technology that allows mapping of EMS-generated mutations without the need to map-cross or back-cross. In the second part of this thesis, using fluorescence activated cell sorting (FACS) coupled with high throughput sequencing, my collaborators and I generated single-base resolution whole genome DNA methylomes, mRNA transcriptomes, and smallRNA transcriptomes for six different populations of cell types in the Arabidopsis root meristem. We were able to discover that the columella is hypermethylated in the CHH context within transposable elements. This hypermethylation is accompanied by upregulation of the RNA-dependent DNA methylation pathway (RdDM), including higher levels of 24-nt silencing RNAs (siRNAs). In summary, our studies demonstrate the versatility of high-throughput sequencing as a method for identifying single mutations or to perform complex comparative genomic analyses.

Behavioural Genetics in the Postgenomics Era

There is growing evidence that the complexity of higher organisms does not correlate with the ‘complexity’ of the genome (the human genome contains fewer protein coding genes than corn, and many genes are preserved across species). Rather, complexity is associated with the complexity of the pathways and processes whereby the cell utilises the deoxyribonucleic acid molecule, and much else, in the process of phenotype formation. These pro- cesses include the activity of the epigenome, noncoding ribonucleic acids, alternative splicing and post-transla- tional modifications. Not accidentally, all of these pro- cesses appear to be of particular importance for the human brain, the most complex organ in nature. Because these processes can be highly environmentally reactive, they are a key to understanding behavioural plasticity and highlight the importance of the developmental process in explaining behavioural outcomes.

microRNA Regulation of Cellular Immunity

Immunity is broadly defined as a mechanism of protection against non-self entities, a process which must be sufficiently robust to both eliminate the initial foreign body and then be maintained over the life of the host. Life-long immunity is impossible without the development of immunological memory, of which a central component is the cellular immune system, or T cells. Cellular immunity hinges upon a naïve T cell pool of sufficient size and breadth to enable Darwinian selection of clones responsive to foreign antigens during an initial encounter. Further, the generation and maintenance of memory T cells is required for rapid clearance responses against repeated insult, and so this small memory pool must be actively maintained by pro-survival cytokine signals over the life of the host.

T cell development, function, and maintenance are regulated on a number of molecular levels through complex regulatory networks. Recently, small non-coding RNAs, miRNAs, have been observed to have profound impacts on diverse aspects of T cell biology by impeding the translation of RNA transcripts to protein. While many miRNAs have been described that alter T cell development or functional differentiation, little is known regarding the role that miRNAs have in T cell maintenance in the periphery at homeostasis.

In Chapter 3 of this dissertation, tools to study miRNA biology and function were developed. First, to understand the effect that miRNA overexpression had on T cell responses, a novel overexpression system was developed to enhance the processing efficiency and ultimate expression of a given miRNA by placing it within an alternative miRNA backbone. Next, a conditional knockout mouse system was devised to specifically delete miR-191 in a cell population expressing recombinase. This strategy was expanded to permit the selective deletion of single miRNAs from within a cluster to discern the effects of specific miRNAs that were previously inaccessible in isolation. Last, to enable the identification of potentially therapeutically viable miRNA function and/or expression modulators, a high-throughput flow cytometry-based screening system utilizing miRNA activity reporters was tested and validated. Thus, several novel and useful tools were developed to assist in the studies described in Chapter 4 and in future miRNA studies.

In Chapter 4 of this dissertation, the role of miR-191 in T cell biology was evaluated. Using tools developed in Chapter 3, miR-191 was observed to be critical for T cell survival following activation-induced cell death, while proliferation was unaffected by alterations in miR-191 expression. Loss of miR-191 led to significant decreases in the numbers of CD4+ and CD8+ T cells in the periphery lymph nodes, but this loss had no impact on the homeostatic activation of either CD4+ or CD8+ cells. These peripheral changes were not caused by gross defects in thymic development, but rather impaired STAT5 phosphorylation downstream of pro-survival cytokine signals. miR-191 does not specifically inhibit STAT5, but rather directly targets the scaffolding protein, IRS1, which in turn alters cytokine-dependent signaling. The defect in peripheral T cell maintenance was exacerbated by the presence of a Bcl-2YFP transgene, which led to even greater peripheral T cell losses in addition to developmental defects. These studies collectively demonstrate that miR-191 controls peripheral T cell maintenance by modulating homeostatic cytokine signaling through the regulation of IRS1 expression and downstream STAT5 phosphorylation.

The studies described in this dissertation collectively demonstrate that miR-191 has a profound role in the maintenance of T cells at homeostasis in the periphery. Importantly, the manipulation of miR-191 altered immune homeostasis without leading to severe immunodeficiency or autoimmunity. As much data exists on the causative agents disrupting active immune responses and the formation of immunological memory, the basic processes underlying the continued maintenance of a functioning immune system must be fully characterized to facilitate the development of methods for promoting healthy immune function throughout the life of the individual. These findings also have powerful implications for the ability of patients with modest perturbations in T cell homeostasis to effectively fight disease and respond to vaccination and may provide valuable targets for therapeutic intervention.

Visualizing Rare Watson-Crick-Like Tautomeric and Anionic Mismatches in DNA and RNA

The central dogma of molecular biology relies on the correct Watson-Crick (WC) geometry of canonical deoxyribonucleic acid (DNA) dG•dC and dA•dT base pairs to replicate and transcribe genetic information with speed and an astonishing level of fidelity. In addition, the Watson-Crick geometry of canonical ribonucleic acid (RNA) rG•rC and rA•rU base pairs is highly conserved to ensure that proteins are translated with high fidelity. However, numerous other potential nucleobase tautomeric and ionic configurations are possible that can give rise to entirely new pairing modes between the nucleotide bases. Very early on, James Watson and Francis Crick recognized their importance and in 1953 postulated that if bases adopted one of their less energetically disfavored tautomeric forms (and later ionic forms) during replication it could lead to the formation of a mismatch with a Watson-Crick-like geometry and could give rise to “natural mutations.”

Since this time numerous studies have provided evidence in support of this hypothesis and have expanded upon it; computational studies have addressed the energetic feasibilities of different nucleobases’ tautomeric and ionic forms in siico; crystallographic studies have trapped different mismatches with WC-like geometries in polymerase or ribosome active sites. However, no direct evidence has been given for (i) the direct existence of these WC-like mismatches in canonical DNA duplex, RNA duplexes, or non-coding RNAs; (ii) which, if any, tautomeric or ionic form stabilizes the WC-like geometry. This thesis utilizes nuclear magnetic resonance (NMR) spectroscopy and rotating frame relaxation dispersion (R1ρ RD) in combination with density functional theory (DFT), biochemical assays, and targeted chemical perturbations to show that (i) dG•dT mismatches in DNA duplexes, as well as rG•rU mismatches RNA duplexes and non-coding RNAs, transiently adopt a WC-like geometry that is stabilized by (ii) an interconnected network of rapidly interconverting rare tautomers and anionic bases. These results support Watson and Crick’s tautomer hypothesis, but additionally support subsequent hypotheses invoking anionic mismatches and ultimately tie them together. This dissertation shows that a common mismatch can adopt a Watson-Crick-like geometry globally, in both DNA and RNA, and whose geometry is stabilized by a kinetically linked network of rare tautomeric and anionic bases. The studies herein also provide compelling evidence for their involvement in spontaneous replication and translation errors.