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.

PRDM9 directs human meiotic crossover hotspots to intergenic sequence motifs, whereas budding yeast hotspots overlap low nucleosome density regions in gene promoters. To investigate hotspots in plants, which lack PRDM9, we used coalescent analysis of Arabidopsis genetic variation. Crossovers increase towards gene promoters and terminators, and hotspots are associated with active chromatin modifications, including H2A.Z, histone H3K4me3, low nucleosome density and low DNA methylation. Hotspot-enriched A-rich and CTT-repeat DNA motifs occur upstream and downstream of transcriptional start respectively. Crossovers are asymmetric around promoters and highest over CTT-motifs and H2A.Z-nucleosomes. Pollen-typing, segregation and cytogenetic analysis show decreased crossovers in the arp6 H2A.Z deposition mutant, at multiple scales. During meiosis H2A.Z and DMC1/RAD51 recombinases form overlapping chromosomal foci. As arp6 reduces DMC1/RAD51 foci, H2A.Z may promote formation or processing of meiotic DNA double-strand breaks. We propose that gene chromatin ancestrally designates hotspots within eukaryotes and PRDM9 is a derived state within vertebrates.

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 59 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 59 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 AGOassociated proteins in influencing epigenetic modifications.

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.

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.

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 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.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.

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.

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.DOI: http://dx.doi.org/10.7554/eLife.03708.001