Silencing of StRIK in potato suggests a role in periderm related to RNA processing and stress.

BACKGROUND: The periderm is a protective barrier crucial for land plant survival, but little is known about genetic factors involved in its development and regulation. Using a transcriptomic approach in the cork oak (Q. suber) periderm, we previously identified an RS2-INTERACTING KH PROTEIN (RIK) homologue of unknown function containing a K homology (KH)-domain RNA-binding protein, as a regulatory candidate gene in the periderm. RESULTS: To gain insight into the function of RIK in the periderm, potato (S. tuberosum) tuber periderm was used as a model: the full-length coding sequence of RIK, hereafter referred to as StRIK, was isolated, the transcript profile analyzed and gene silencing in potato performed to analyze the silencing effects on periderm anatomy and transcriptome. The StRIK transcript accumulated in all vegetative tissues studied, including periderm and other suberized tissues such as root and also in wounded tissues. Downregulation of StRIK in potato by RNA interference (StRIK-RNAi) did not show any obvious effects on tuber periderm anatomy but, unlike Wild type, transgenic plants flowered. Global transcript profiling of the StRIK-RNAi periderm did show altered expression of genes associated with RNA metabolism, stress and signaling, mirroring the biological processes found enriched within the in silico co-expression network of the Arabidopsis orthologue. CONCLUSIONS: The ubiquitous expression of StRIK transcript, the flower associated phenotype and the differential expression of StRIK-RNAi periderm point out to a general regulatory role of StRIK in diverse plant developmental processes. The transcriptome analysis suggests that StRIK might play roles in RNA maturation and stress response in the periderm.

Arabidopsis sirtuins and poly(ADP-ribose) polymerases regulate gene expression in the day but do not affect circadian rhythms.

Nicotinamide-adenine dinucleotide (NAD) is involved in redox homeostasis and acts as a substrate for NADases, including poly(ADP-ribose) polymerases (PARPs) that add poly(ADP-ribose) polymers to proteins and DNA, and sirtuins that deacetylate proteins. Nicotinamide, a by-product of NADases increases circadian period in both plants and animals. In mammals, the effect of nicotinamide on circadian period might be mediated by the PARPs and sirtuins because they directly bind to core circadian oscillator genes. We have investigated whether PARPs and sirtuins contribute to the regulation of the circadian oscillator in Arabidopsis. We found no evidence that PARPs and sirtuins regulate the circadian oscillator of Arabidopsis or are involved in the response to nicotinamide. RNA-seq analysis indicated that PARPs regulate the expression of only a few genes, including FLOWERING LOCUS C. However, we found profound effects of reduced sirtuin 1 expression on gene expression during the day but not at night, and an embryo lethal phenotype in knockouts. Our results demonstrate that PARPs and sirtuins are not associated with NAD regulation of the circadian oscillator and that sirtuin 1 is associated with daytime regulation of gene expression.

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

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.

Juxtaposition of heterozygous and homozygous regions causes reciprocal crossover remodelling via interference during Arabidopsis meiosis.

During meiosis homologous chromosomes undergo crossover recombination. Sequence differences between homologs can locally inhibit crossovers. Despite this, nucleotide diversity and population-scaled recombination are positively correlated in eukaryote genomes. To investigate interactions between heterozygosity and recombination we crossed Arabidopsis lines carrying fluorescent crossover reporters to 32 diverse accessions and observed hybrids with significantly higher and lower crossovers than homozygotes. Using recombinant populations derived from these crosses we observed that heterozygous regions increase crossovers when juxtaposed with homozygous regions, which reciprocally decrease. Total crossovers measured by chiasmata were unchanged when heterozygosity was varied, consistent with homeostatic control. We tested the effects of heterozygosity in mutants where the balance of interfering and non-interfering crossover repair is altered. Crossover remodeling at homozygosity-heterozygosity junctions requires interference, and non-interfering repair is inefficient in heterozygous regions. As a consequence, heterozygous regions show stronger crossover interference. Our findings reveal how varying homolog polymorphism patterns can shape meiotic recombination.

High-throughput analysis of meiotic crossover frequency and interference via flow cytometry of fluorescent pollen in Arabidopsis thaliana.

During meiosis, reciprocal exchange between homologous chromosomes occurs as a result of crossovers (COs). CO frequency varies within genomes and is subject to genetic, epigenetic and environmental control. As robust measurement of COs is limited by their low numbers, typically 1-2 per chromosome, we adapted flow cytometry for use with Arabidopsis transgenic fluorescent protein-tagged lines (FTLs) that express eCFP, dsRed or eYFP fluorescent proteins in pollen. Segregation of genetically linked transgenes encoding fluorescent proteins of distinct colors can be used to detect COs. The fluorescence of up to 80,000 pollen grains per individual plant can be measured in 10-15 min using this protocol. A key element of CO control is interference, which inhibits closely spaced COs. We describe a three-color assay for the measurement of CO frequency in adjacent intervals and calculation of CO interference. We show that this protocol can be used to detect changes in CO frequency and interference in the fancm zip4 double mutant. By enabling high-throughput measurement of CO frequency and interference, these methods will facilitate genetic dissection of meiotic recombination control.

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.

Empirical Bayesian analysis of paired high-throughput sequencing data with a beta-binomial distribution.

BACKGROUND: Pairing of samples arises naturally in many genomic experiments; for example, gene expression in tumour and normal tissue from the same patients. Methods for analysing high-throughput sequencing data from such experiments are required to identify differential expression, both within paired samples and between pairs under different experimental conditions. RESULTS: We develop an empirical Bayesian method based on the beta-binomial distribution to model paired data from high-throughput sequencing experiments. We examine the performance of this method on simulated and real data in a variety of scenarios. Our methods are implemented as part of the RbaySeq package (versions 1.11.6 and greater) available from Bioconductor (http://www.bioconductor.org). CONCLUSIONS: We compare our approach to alternatives based on generalised linear modelling approaches and show that our method offers significant gains in performance on simulated data. In testing on real data from oral squamous cell carcinoma patients, we discover greater enrichment of previously identified head and neck squamous cell carcinoma associated gene sets than has previously been achieved through a generalised linear modelling approach, suggesting that similar gains in performance may be found in real data. Our methods thus show real and substantial improvements in analyses of high-throughput sequencing data from paired samples.

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.

Identifying small interfering RNA loci from high-throughput sequencing data.

MOTIVATION: Small interfering RNAs (siRNAs) are produced from much longer sequences of double-stranded RNA precursors through cleavage by Dicer or a Dicer-like protein. These small RNAs play a key role in genetic and epigenetic regulation; however, a full understanding of the mechanisms by which they operate depends on the characterization of the precursors from which they are derived. RESULTS: High-throughput sequencing of small RNA populations allows the locations of the double-stranded RNA precursors to be inferred. We have developed methods to analyse small RNA sequencing data from multiple biological sources, taking into account replicate information, to identify robust sets of siRNA precursors. Our methods show good performance on both a set of small RNA sequencing data in Arabidopsis thaliana and simulated datasets. AVAILABILITY: Our methods are available as the Bioconductor (www.bioconductor.org) package segmentSeq (version 1.5.6 and above).

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.

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.

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

BACKGROUND: High 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. RESULTS: We 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. CONCLUSIONS: Our 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.

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 Arabidopsis RNA-directed DNA methylation argonautes functionally diverge based on their expression and interaction with target loci.

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.

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.

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

BACKGROUND: Worldwide, approximately two billion people are chronically infected with Toxoplasma gondii with largely unknown consequences. METHODS: To 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 alphaPD1 ligand were studied. RESULTS: Twelve 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 Alzheimer's 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. CONCLUSION: In 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.

Epigenomic modifications predict active promoters and gene structure in Toxoplasma gondii.

Mechanisms of gene regulation are poorly understood in Apicomplexa, a phylum that encompasses deadly human pathogens like Plasmodium and Toxoplasma. Initial studies suggest that epigenetic phenomena, including histone modifications and chromatin remodeling, have a profound effect upon gene expression and expression of virulence traits. Using the model organism Toxoplasma gondii, we characterized the epigenetic organization and transcription patterns of a contiguous 1% of the T. gondii genome using custom oligonucleotide microarrays. We show that methylation and acetylation of histones H3 and H4 are landmarks of active promoters in T. gondii that allow us to deduce the position and directionality of gene promoters with >95% accuracy. These histone methylation and acetylation "activation" marks are strongly associated with gene expression. We also demonstrate that the pattern of histone H3 arginine methylation distinguishes certain promoters, illustrating the complexity of the histone modification machinery in Toxoplasma. By integrating epigenetic data, gene prediction analysis, and gene expression data from the tachyzoite stage, we illustrate feasibility of creating an epigenomic map of T. gondii tachyzoite gene expression. Further, we illustrate the utility of the epigenomic map to empirically and biologically annotate the genome and show that this approach enables identification of previously unknown genes. Thus, our epigenomics approach provides novel insights into regulation of gene expression in the Apicomplexa. In addition, with its compact genome, genetic tractability, and discrete life cycle stages, T. gondii provides an important new model to study the evolutionarily conserved components of the histone code.

Common inheritance of chromosome Ia associated with clonal expansion of Toxoplasma gondii.

Toxoplasma gondii is a globally distributed protozoan parasite that can infect virtually all warm-blooded animals and humans. Despite the existence of a sexual phase in the life cycle, T. gondii has an unusual population structure dominated by three clonal lineages that predominate in North America and Europe, (Types I, II, and III). These lineages were founded by common ancestors approximately10,000 yr ago. The recent origin and widespread distribution of the clonal lineages is attributed to the circumvention of the sexual cycle by a new mode of transmission-asexual transmission between intermediate hosts. Asexual transmission appears to be multigenic and although the specific genes mediating this trait are unknown, it is predicted that all members of the clonal lineages should share the same alleles. Genetic mapping studies suggested that chromosome Ia was unusually monomorphic compared with the rest of the genome. To investigate this further, we sequenced chromosome Ia and chromosome Ib in the Type I strain, RH, and the Type II strain, ME49. Comparative genome analyses of the two chromosomal sequences revealed that the same copy of chromosome Ia was inherited in each lineage, whereas chromosome Ib maintained the same high frequency of between-strain polymorphism as the rest of the genome. Sampling of chromosome Ia sequence in seven additional representative strains from the three clonal lineages supports a monomorphic inheritance, which is unique within the genome. Taken together, our observations implicate a specific combination of alleles on chromosome Ia in the recent origin and widespread success of the clonal lineages of T. gondii.

Characterisation of conserved non-coding sequences in vertebrate genomes using bioinformatics, statistics and functional studies.

We recently identified approximately 1400 conserved non-coding elements (CNEs) shared by the genomes of fugu (Takifugu rubripes) and human that appear to be associated with developmental regulation in vertebrates [Woolfe, A., Goodson, M., Goode, D.K., Snell, P., McEwen, G.K., Vavouri, T., Smith, S.F., North, P., Callaway, H., Kelly, K., Walter, K., Abnizova, I., Gilks, W., Edwards, Y.J.K., Cooke, J.E., Elgar, G., 2005. Highly conserved non-coding sequences are associated with vertebrate development. PLoS Biol. 3 (1), e7]. This study encompassed a multi-disciplinary approach using bioinformatics, statistical methods and functional assays to identify and characterise the CNEs. Using an in vivo enhancer assay, over 90% of tested CNEs up-regulate tissue-specific GFP expression. Here we review our group's research in the field of characterising non-coding sequences conserved in vertebrates. We take this opportunity to discuss our research in progress and present some results of new and additional analyses. These include a phylogenomics analysis of CNEs, sequence conservation patterns in vertebrate CNEs and the distribution of human SNPs in the CNEs. We highlight the usefulness of the CNE dataset to help correlate genetic variation in health and disease. We also discuss the functional analysis using the enhancer assay and the enrichment of predicted transcription factor binding sites for two CNEs. Public access to the CNEs plus annotation is now possible and is described. The content of this review was presented by Dr. Y.J.K. Edwards at the TODAI International Symposium on Functional Genomics of the Pufferfish, Tokyo, Japan, 3-6 November 2004.