Krystyna A. Kelly

baySeq: Empirical Bayesian methods for identifying differential expression in sequence count data

BackgroundHigh throughput sequencing has become an important technology for studying expression levels in many types of genomic, and particularly transcriptomic, data. One key way of analysing such data is to look for elements of the data which display particular patterns of differential expression in order to take these forward for further analysis and validation.ResultsWe propose a framework for defining patterns of differential expression and develop a novel algorithm, baySeq, which uses an empirical Bayes approach to detect these patterns of differential expression within a set of sequencing samples. The method assumes a negative binomial distribution for the data and derives an empirically determined prior distribution from the entire dataset. We examine the performance of the method on real and simulated data.ConclusionsOur method performs at least as well, and often better, than existing methods for analyses of pairwise differential expression in both real and simulated data. When we compare methods for the analysis of data from experimental designs involving multiple sample groups, our method again shows substantial gains in performance. We believe that this approach thus represents an important step forward for the analysis of count data from sequencing experiments.

Uniparental expression of PolIV-dependent siRNAs in developing endosperm of Arabidopsis

Most eukaryotes produce small RNA (sRNA) mediators of gene silencing that bind to Argonaute proteins and guide them, by base pairing, to an RNA target. MicroRNAs (miRNAs) that normally target messenger RNAs for degradation or translational arrest are the best-understood class of sRNAs. However, in Arabidopsis thaliana flowers, miRNAs account for only 5% of the sRNA mass and less than 0.1% of the sequence complexity. The remaining sRNAs form a complex population of more than 100,000 different small interfering RNAs (siRNAs) transcribed from thousands of loci. The biogenesis of most of the siRNAs in Arabidopsis are dependent on RNA polymerase IV (PolIV), a homologue of DNA-dependent RNA polymerase II. A subset of these PolIV-dependent (p4)-siRNAs are involved in stress responses, and others are associated with epigenetic modifications to DNA or chromatin; however, the biological role is not known for most of them. Here we show that the predominant phase of p4-siRNA accumulation is initiated in the maternal gametophyte and continues during seed development. Expression of p4-siRNAs in developing endosperm is specifically from maternal chromosomes. Our results provide the first evidence for a link between genomic imprinting and RNA silencing in plants.

The Arabidopsis RNA-Directed DNA Methylation Argonautes Functionally Diverge Based on Their Expression and Interaction with Target Loci

This study examines the basis of functional divergence amongst the closely related Arabidopsis AGO4, AGO6, and AGO9 proteins. Using associated small RNAs, these proteins direct DNA epigenetic modifications in a tissue-specific manner. Argonaute (AGO) effectors of RNA silencing bind small RNA (sRNA) molecules and mediate mRNA cleavage, translational repression, or epigenetic DNA modification. In many organisms, these targeting mechanisms are devolved to different products of AGO multigene families. To investigate the basis of AGO functional diversification, we characterized three closely related Arabidopsis thaliana AGOs (AGO4, AGO6, and AGO9) implicated in RNA-directed DNA methylation. All three AGOs bound 5′ adenosine 24-nucleotide sRNAs, but each exhibited different preferences for sRNAs from different heterochromatin-associated loci. This difference was reduced when AGO6 and AGO9 were expressed from the AGO4 promoter, indicating that the functional diversification was partially due to differential expression of the corresponding genes. However, the AGO4-directed pattern of sRNA accumulation and DNA methylation was not fully recapitulated with AGO6 or AGO9 expressed from the AGO4 promoter. Here, we show that sRNA length and 5′ nucleotide do not account for the observed functional diversification of these AGOs. Instead, the selectivity of sRNA binding is determined by the coincident expression of the AGO and sRNA-generating loci, and epigenetic modification is influenced by interactions between the AGO protein and the different target loci. These findings highlight the importance of tissue specificity and AGO-associated proteins in influencing epigenetic modifications.

Arabidopsis meiotic crossover hot spots overlap with H2A.Z nucleosomes at gene promoters

PRDM9 directs human meiotic crossover hot spots to intergenic sequence motifs, whereas budding yeast hot spots overlap regions of low nucleosome density (LND) in gene promoters. To investigate hot spots in plants, which lack PRDM9, we used coalescent analysis of genetic variation in Arabidopsis thaliana. Crossovers increased toward gene promoters and terminators, and hot spots were associated with active chromatin modifications, including H2A.Z, histone H3 Lys4 trimethylation (H3K4me3), LND and low DNA methylation. Hot spot–enriched A-rich and CTT-repeat DNA motifs occurred upstream and downstream, respectively, of transcriptional start sites. Crossovers were asymmetric around promoters and were most frequent over CTT-repeat motifs and H2A.Z nucleosomes. Pollen typing, segregation and cytogenetic analysis showed decreased numbers of crossovers in the arp6 H2A.Z deposition mutant at multiple scales. During meiosis, H2A.Z forms overlapping chromosomal foci with the DMC1 and RAD51 recombinases. As arp6 reduced the number of DMC1 or RAD51 foci, H2A.Z may promote the formation or processing of meiotic DNA double-strand breaks. We propose that gene chromatin ancestrally designates hot spots within eukaryotes and PRDM9 is a derived state within vertebrates.

Neurological and behavioral abnormalities, ventricular dilatation, altered cellular functions, inflammation, and neuronal injury in brains of mice due to common, persistent, parasitic infection

BackgroundWorldwide, approximately two billion people are chronically infected with Toxoplasma gondii with largely unknown consequences.MethodsTo better understand long-term effects and pathogenesis of this common, persistent brain infection, mice were infected at a time in human years equivalent to early to mid adulthood and studied 5–12 months later. Appearance, behavior, neurologic function and brain MRIs were studied. Additional analyses of pathogenesis included: correlation of brain weight and neurologic findings; histopathology focusing on brain regions; full genome microarrays; immunohistochemistry characterizing inflammatory cells; determination of presence of tachyzoites and bradyzoites; electron microscopy; and study of markers of inflammation in serum. Histopathology in genetically resistant mice and cytokine and NRAMP knockout mice, effects of inoculation of isolated parasites, and treatment with sulfadiazine or αPD1 ligand were studied.ResultsTwelve months after infection, a time equivalent to middle to early elderly ages, mice had behavioral and neurological deficits, and brain MRIs showed mild to moderate ventricular dilatation. Lower brain weight correlated with greater magnitude of neurologic abnormalities and inflammation. Full genome microarrays of brains reflected inflammation causing neuronal damage (Gfap), effects on host cell protein processing (ubiquitin ligase), synapse remodeling (Complement 1q), and also increased expression of PD-1L (a ligand that allows persistent LCMV brain infection) and CD 36 (a fatty acid translocase and oxidized LDL receptor that mediates innate immune response to beta amyloid which is associated with pro-inflammation in Alzheimers disease). Immunostaining detected no inflammation around intra-neuronal cysts, practically no free tachyzoites, and only rare bradyzoites. Nonetheless, there were perivascular, leptomeningeal inflammatory cells, particularly contiguous to the aqueduct of Sylvius and hippocampus, CD4+ and CD8+ T cells, and activated microglia in perivascular areas and brain parenchyma. Genetically resistant, chronically infected mice had substantially less inflammation.ConclusionIn outbred mice, chronic, adult acquired T. gondii infection causes neurologic and behavioral abnormalities secondary to inflammation and loss of brain parenchyma. Perivascular inflammation is prominent particularly contiguous to the aqueduct of Sylvius and hippocampus. Even resistant mice have perivascular inflammation. This mouse model of chronic T. gondii infection raises questions of whether persistence of this parasite in brain can cause inflammation or neurodegeneration in genetically susceptible hosts.

Epigenetic Remodeling of Meiotic Crossover Frequency in Arabidopsis thaliana DNA Methyltransferase Mutants

Meiosis is a specialized eukaryotic cell division that generates haploid gametes required for sexual reproduction. During meiosis, homologous chromosomes pair and undergo reciprocal genetic exchange, termed crossover (CO). Meiotic CO frequency varies along the physical length of chromosomes and is determined by hierarchical mechanisms, including epigenetic organization, for example methylation of the DNA and histones. Here we investigate the role of DNA methylation in determining patterns of CO frequency along Arabidopsis thaliana chromosomes. In A. thaliana the pericentromeric regions are repetitive, densely DNA methylated, and suppressed for both RNA polymerase-II transcription and CO frequency. DNA hypomethylated methyltransferase1 (met1) mutants show transcriptional reactivation of repetitive sequences in the pericentromeres, which we demonstrate is coupled to extensive remodeling of CO frequency. We observe elevated centromere-proximal COs in met1, coincident with pericentromeric decreases and distal increases. Importantly, total numbers of CO events are similar between wild type and met1, suggesting a role for interference and homeostasis in CO remodeling. To understand recombination distributions at a finer scale we generated CO frequency maps close to the telomere of chromosome 3 in wild type and demonstrate an elevated recombination topology in met1. Using a pollen-typing strategy we have identified an intergenic nucleosome-free CO hotspot 3a, and we demonstrate that it undergoes increased recombination activity in met1. We hypothesize that modulation of 3a activity is caused by CO remodeling driven by elevated centromeric COs. These data demonstrate how regional epigenetic organization can pattern recombination frequency along eukaryotic chromosomes.

Putative Arabidopsis THO/TREX mRNA export complex is involved in transgene and endogenous siRNA biosynthesis

RNA silencing in plants and some animals has a non–cell-autonomous effect due to an RNA signal that moves between cells or organs. To identify unique factors involved in this process, we analyzed a group of Arabidopsis mutants with defective spread of RNA silencing from a transgene expressed specifically in the phloem. These mutants accumulated reduced amounts of small interfering (si)RNA from the transgene locus and from endogenous loci TAS1, TAS2, and an inverted repeat locus IR71. The defect in TAS1 and TAS2 siRNA biogenesis is in the processing of a long siRNA precursor. We mapped the mutations to a gene encoding the Arabidopsis homolog of a protein, TEX1, which is involved in intracellular transport of RNA in animals. TEX1 is a component of the THO/TREX complex, and we show that the Arabidopsis TEX1 interacts with other predicted components of a THO/TREX complex. Correspondingly, we found at least two other components of the Arabidopsis THO core complex that are involved in RNA silencing. To reconcile the effect of these mutations on transgene and endogenous gene siRNA, we propose a mechanism in which THO/TREX processes or transports a long RNA molecule so that it can be a template for secondary siRNA production.

The De Novo Cytosine Methyltransferase DRM2 Requires Intact UBA Domains and a Catalytically Mutated Paralog DRM3 during RNA-Directed DNA Methylation in Arabidopsis thaliana

Eukaryotic DNA cytosine methylation can be used to transcriptionally silence repetitive sequences, including transposons and retroviruses. This silencing is stable between cell generations as cytosine methylation is maintained epigenetically through DNA replication. The Arabidopsis thaliana Dnmt3 cytosine methyltransferase ortholog DOMAINS REARRANGED METHYLTRANSFERASE2 (DRM2) is required for establishment of small interfering RNA (siRNA) directed DNA methylation. In mammals PIWI proteins and piRNA act in a convergently evolved RNA–directed DNA methylation system that is required to repress transposon expression in the germ line. De novo methylation may also be independent of RNA interference and small RNAs, as in Neurospora crassa. Here we identify a clade of catalytically mutated DRM2 paralogs in flowering plant genomes, which in A.thaliana we term DOMAINS REARRANGED METHYLTRANSFERASE3 (DRM3). Despite being catalytically mutated, DRM3 is required for normal maintenance of non-CG DNA methylation, establishment of RNA–directed DNA methylation triggered by repeat sequences and accumulation of repeat-associated small RNAs. Although the mammalian catalytically inactive Dnmt3L paralogs act in an analogous manner, phylogenetic analysis indicates that the DRM and Dnmt3 protein families diverged independently in plants and animals. We also show by site-directed mutagenesis that both the DRM2 N-terminal UBA domains and C-terminal methyltransferase domain are required for normal RNA–directed DNA methylation, supporting an essential targeting function for the UBA domains. These results suggest that plant and mammalian RNA–directed DNA methylation systems consist of a combination of ancestral and convergent features.

Rice perception of symbiotic arbuscular mycorrhizal fungi requires the karrikin receptor complex

Early stages of a beneficial relationship Plants benefit from widespread symbiosis with arbuscular mycorrhizal fungi. This symbiosis between plant and fungus aids plants in capturing mineral and micronutrients from the soil. Gutjahr et al. have now identified a component of an intracellular receptor, the hydrolase DWARF 14 LIKE, required in rice roots for initiating the symbiosis. A similar receptor detects karrikins in smoke that signal opportunity for fireweed to grow after a forest fire. Science, this issue p. 1521 Widely beneficial symbiosis between plant and fungi shares signaling components with wildfire ephemerals. In terrestrial ecosystems, plants take up phosphate predominantly via association with arbuscular mycorrhizal fungi (AMF). We identified loss of responsiveness to AMF in the rice (Oryza sativa) mutant hebiba, reflected by the absence of physical contact and of characteristic transcriptional responses to fungal signals. Among the 26 genes deleted in hebiba, DWARF 14 LIKE is, the one responsible for loss of symbiosis . It encodes an alpha/beta-fold hydrolase, that is a component of an intracellular receptor complex involved in the detection of the smoke compound karrikin. Our finding reveals an unexpected plant recognition strategy for AMF and a previously unknown signaling link between symbiosis and plant development.

RNAcentral: A vision for an international database of RNA sequences

During the last decade there has been a great increase in the number of noncoding RNA genes identified, including new classes such as microRNAs and piRNAs. There is also a large growth in the amount of experimental characterization of these RNA components. Despite this growth in information, it is still difficult for researchers to access RNA data, because key data resources for noncoding RNAs have not yet been created. The most pressing omission is the lack of a comprehensive RNA sequence database, much like UniProt, which provides a comprehensive set of protein knowledge. In this article we propose the creation of a new open public resource that we term RNAcentral, which will contain a comprehensive collection of RNA sequences and fill an important gap in the provision of biomedical databases. We envision RNA researchers from all over the world joining a federated RNAcentral network, contributing specialized knowledge and databases. RNAcentral would centralize key data that are currently held across a variety of databases, allowing researchers instant access to a single, unified resource. This resource would facilitate the next generation of RNA research and help drive further discoveries, including those that improve food production and human and animal health. We encourage additional RNA database resources and research groups to join this effort. We aim to obtain international network funding to further this endeavor.