Plant-RRBS, a bisulfite and next-generation sequencing-based methylome profiling method enriching for coverage of cytosine positions

Background: Cytosine methylation in plant genomes is important for the regulation of gene transcription and transposon activity. Genome-wide methylomes are studied upon mutation of the DNA methyltransferases, adaptation to environmental stresses or during development. However, from basic biology to breeding programs, there is a need to monitor multiple samples to determine transgenerational methylation inheritance or differential cytosine methylation. Methylome data obtained by sodium hydrogen sulfite (bisulfite)-conversion and next-generation sequencing (NGS) provide genome- wide information on cytosine methylation. However, a profiling method that detects cytosine methylation state dispersed over the genome would allow high-throughput analysis of multiple plant samples with distinct epigenetic signatures. We use specific restriction endonucleases to enrich for cytosine coverage in a bisulfite and NGS-based profiling method, which was compared to whole-genome bisulfite sequencing of the same plant material. Methods: We established an effective methylome profiling method in plants, termed plant-reduced representation bisulfite sequencing (plant-RRBS), using optimized double restriction endonuclease digestion, fragment end repair, adapter ligation, followed by bisulfite conversion, PCR amplification and NGS. We report a performant laboratory protocol and a straightforward bioinformatics data analysis pipeline for plant-RRBS, applicable for any reference-sequenced plant species. Results: As a proof of concept, methylome profiling was performed using an Oryza sativa ssp. indica pure breeding line and a derived epigenetically altered line (epiline). Plant-RRBS detects methylation levels at tens of millions of cytosine positions deduced from bisulfite conversion in multiple samples. To evaluate the method, the coverage of cytosine positions, the intra-line similarity and the differential cytosine methylation levels between the pure breeding line and the epiline were determined. Plant-RRBS reproducibly covers commonly up to one fourth of the cytosine positions in the rice genome when using MspI-DpnII within a group of five biological replicates of a line. The method predominantly detects cytosine methylation in putative promoter regions and not-annotated regions in rice. Conclusions: Plant-RRBS offers high-throughput and broad, genome- dispersed methylation detection by effective read number generation obtained from reproducibly covered genome fractions using optimized endonuclease combinations, facilitating comparative analyses of multi-sample studies for cytosine methylation and transgenerational stability in experimental material and plant breeding populations

Exceptional river gorge formation from unexceptional floods

An understanding of rates and mechanisms of incision and knickpoint retreat in bedrock rivers is fundamental to perceptions of landscape response to external drivers, yet only sparse field data are available. Here we present eye witness accounts and quantitative surveys of rapid, amphitheatre-headed gorge formation in unweathered granite from the overtopping of a rock-cut dam spillway by small-moderate floods (~100–1,500 m3 s−1). The amount of erosion demonstrates no relationship with flood magnitude or bedload availability. Instead, structural pattern of the bedrock through faults and joints appears to be the primary control on landscape change. These discontinuities facilitate rapid erosion (>270 m headward retreat; ~100 m incision; and ~160 m widening over 6 years) principally through fluvial plucking and block topple. The example demonstrates the potential for extremely rapid transient bedrock erosion even when rocks are mechanically strong and flood discharges are moderate. These observations are relevant to perceived models of gorge formation and knickpoint retreat

Study of the reaction e^{+}e^{-} -->J/psi\pi^{+}\pi^{-} via initial-state radiation at BaBar

We study the process $e^+e^-\to J/\psi\pi^{+}\pi^{-}$ with initial-state-radiation events produced at the PEP-II asymmetric-energy collider. The data were recorded with the BaBar detector at center-of-mass energies 10.58 and 10.54 GeV, and correspond to an integrated luminosity of 454 $\mathrm{fb^{-1}}$. We investigate the $J/\psi \pi^{+}\pi^{-}$ mass distribution in the region from 3.5 to 5.5 $\mathrm{GeV/c^{2}}$. Below 3.7 $\mathrm{GeV/c^{2}}$ the $\psi(2S)$ signal dominates, and above 4 $\mathrm{GeV/c^{2}}$ there is a significant peak due to the Y(4260). A fit to the data in the range 3.74 -- 5.50 $\mathrm{GeV/c^{2}}$ yields a mass value $4244 \pm 5$ (stat) $\pm 4$ (syst)$\mathrm{MeV/c^{2}}$ and a width value $114 ^{+16}_{-15}$ (stat)$\pm 7$(syst)$\mathrm{MeV}$ for this state. We do not confirm the report from the Belle collaboration of a broad structure at 4.01 $\mathrm{GeV/c^{2}}$. In addition, we investigate the $\pi^{+}\pi^{-}$ system which results from Y(4260) decay

Generation of a genomic tiling array of the human Major Histocompatibility Complex (MHC) and its application for DNA methylation analysis

Background: The major histocompatibility complex (MHC) is essential for human immunity and is highly associated with common diseases, including cancer. While the genetics of the MHC has been studied intensively for many decades, very little is known about the epigenetics of this most polymorphic and disease-associated region of the genome.Methods: To facilitate comprehensive epigenetic analyses of this region, we have generated a genomic tiling array of 2 Kb resolution covering the entire 4 Mb MHC region. The array has been designed to be compatible with chromatin immunoprecipitation (ChIP), methylated DNA immunoprecipitation (MeDIP), array comparative genomic hybridization (aCGH) and expression profiling, including of non-coding RNAs. The array comprises 7832 features, consisting of two replicates of both forward and reverse strands of MHC amplicons and appropriate controls.Results: Using MeDIP, we demonstrate the application of the MHC array for DNA methylation profiling and the identification of tissue-specific differentially methylated regions (tDMRs). Based on the analysis of two tissues and two cell types, we identified 90 tDMRs within the MHC and describe their characterisation.Conclusion: A tiling array covering the MHC region was developed and validated. Its successful application for DNA methylation profiling indicates that this array represents a useful tool for molecular analyses of the MHC in the context of medical genomics

Phospho-ERK and AKT status, but not KRAS mutation status, are associated with outcomes in rectal cancer treated with chemoradiotherapy

<p>Abstract</p> <p>Background</p> <p><it>KRAS </it>mutations may predict poor response to radiotherapy. Downstream events from <it>KRAS</it>, such as activation of <it>BRAF</it>, AKT and ERK, may also confer prognostic information but have not been tested in rectal cancer (RC). Our objective was to explore the relationships of <it>KRAS </it>and <it>BRAF </it>mutation status with p-AKT and p-ERK and outcomes in RC.</p> <p>Methods</p> <p>Pre-radiotherapy RC tumor biopsies were evaluated. <it>KRAS </it>and <it>BRAF </it>mutations were assessed by pyrosequencing; p-AKT and p-ERK expression by immunohistochemistry.</p> <p>Results</p> <p>Of 70 patients, mean age was 58; 36% stage II, 56% stage III, and 9% stage IV. Responses to neoadjuvant chemoradiotherapy: 64% limited, 19% major, and 17% pathologic complete response. 64% were <it>KRAS </it>WT, 95% were <it>BRAF </it>WT. High p-ERK levels were associated with improved OS but not for p-AKT. High levels of p-AKT and p-ERK expression were associated with better responses. <it>KRAS </it>WT correlated with lower p-AKT expression but not p-ERK expression. No differences in OS, residual disease, or tumor downstaging were detected by <it>KRAS </it>status.</p> <p>Conclusions</p> <p><it>KRAS </it>mutation was not associated with lesser response to chemoradiotherapy or worse OS. High p-ERK expression was associated with better OS and response. Higher p-AKT expression was correlated with better response but not OS.</p

Sensitivity of jarrah (Eucalyptus marginata) to phosphate, phosphite, and arsenate pulses as influenced by fungal symbiotic associations

Many plant species adapted to P-impoverished soils, including jarrah (Eucalyptus marginata), develop toxicity symptoms when exposed to high doses of phosphate (Pi) and its analogs such as phosphite (Phi) and arsenate (AsV). The present study was undertaken to investigate the effects of fungal symbionts Scutellospora calospora, Scleroderma sp., and Austroboletus occidentalis on the response of jarrah to highly toxic pulses (1.5 mmol kg−1 soil) of Pi, Phi, and AsV. S. calospora formed an arbuscular mycorrhizal (AM) symbiosis while both Scleroderma sp. and A. occidentalis established a non-colonizing symbiosis with jarrah plants. All these interactions significantly improved jarrah growth and Pi uptake under P-limiting conditions. The AM fungal colonization naturally declines in AM-eucalypt symbioses after 2–3 months; however, in the present study, the high Pi pulse inhibited the decline of AM fungal colonization in jarrah. Four weeks after exposure to the Pi pulse, plants inoculated with S. calospora had significantly lower toxicity symptoms compared to non-mycorrhizal (NM) plants, and all fungal treatments induced tolerance against Phi toxicity in jarrah. However, no tolerance was observed for AsV-treated plants even though all inoculated plants had significantly lower shoot As concentrations than the NM plants. The transcript profile of five jarrah high-affinity phosphate transporter (PHT1 family) genes in roots was not altered in response to any of the fungal species tested. Interestingly, plants exposed to high Pi supplies for 1 day did not have reduced transcript levels for any of the five PHT1 genes in roots, and transcript abundance of four PHT1 genes actually increased. It is therefore suggested that jarrah, and perhaps other P-sensitive perennial species, respond positively to Pi available in the soil solution through increasing rather than decreasing the expression of selected PHT1 genes. Furthermore, Scleroderma sp. can be considered as a fungus with dual functional capacity capable of forming both ectomycorrhizal and non-colonizing associations, where both pathways are always accompanied by evident growth and nutritional benefits

Transcriptomic landscape of breast cancers through mRNA sequencing

Breast cancer is a heterogeneous disease with a poorly defined genetic landscape, which poses a major challenge in diagnosis and treatment. By massively parallel mRNA sequencing, we obtained 1.2 billion reads from 17 individual human tissues belonging to TNBC, Non-TNBC, and HER2-positive breast cancers and defined their comprehensive digital transcriptome for the first time. Surprisingly, we identified a high number of novel and unannotated transcripts, revealing the global breast cancer transcriptomic adaptations. Comparative transcriptomic analyses elucidated differentially expressed transcripts between the three breast cancer groups, identifying several new modulators of breast cancer. Our study also identified common transcriptional regulatory elements, such as highly abundant primary transcripts, including osteonectin, RACK1, calnexin, calreticulin, FTL, and B2M, and “genomic hotspots” enriched in primary transcripts between the three groups. Thus, our study opens previously unexplored niches that could enable a better understanding of the disease and the development of potential intervention strategies

Evolutionary Dynamics of the Ty3/Gypsy LTR Retrotransposons in the Genome of Anopheles gambiae

Ty3/gypsy elements represent one of the most abundant and diverse LTR-retrotransposon (LTRr) groups in the Anopheles gambiae genome, but their evolutionary dynamics have not been explored in detail. Here, we conduct an in silico analysis of the distribution and abundance of the full complement of 1045 copies in the updated AgamP3 assembly. Chromosomal distribution of Ty3/gypsy elements is inversely related to arm length, with densities being greatest on the X, and greater on the short versus long arms of both autosomes. Taking into account the different heterochromatic and euchromatic compartments of the genome, our data suggest that the relative abundance of Ty3/gypsy LTRrs along each chromosome arm is determined mainly by the different proportions of heterochromatin, particularly pericentric heterochromatin, relative to total arm length. Additionally, the breakpoint regions of chromosomal inversion 2La appears to be a haven for LTRrs. These elements are underrepresented more than 7-fold in euchromatin, where 33% of the Ty3/gypsy copies are associated with genes. The euchromatin on chromosome 3R shows a faster turnover rate of Ty3/gypsy elements, characterized by a deficit of proviral sequences and the lowest average sequence divergence of any autosomal region analyzed in this study. This probably reflects a principal role of purifying selection against insertion for the preservation of longer conserved syntenyc blocks with adaptive importance located in 3R. Although some Ty3/gypsy LTRrs show evidence of recent activity, an important fraction are inactive remnants of relatively ancient insertions apparently subject to genetic drift. Consistent with these computational predictions, an analysis of the occupancy rate of putatively older insertions in natural populations suggested that the degenerate copies have been fixed across the species range in this mosquito, and also are shared with the sibling species Anopheles arabiensis

Parental Genome Dosage Imbalance Deregulates Imprinting in Arabidopsis

In mammals and in plants, parental genome dosage imbalance deregulates embryo growth and might be involved in reproductive isolation between emerging new species. Increased dosage of maternal genomes represses growth while an increased dosage of paternal genomes has the opposite effect. These observations led to the discovery of imprinted genes, which are expressed by a single parental allele. It was further proposed in the frame of the parental conflict theory that parental genome imbalances are directly mirrored by antagonistic regulations of imprinted genes encoding maternal growth inhibitors and paternal growth enhancers. However these hypotheses were never tested directly. Here, we investigated the effect of parental genome imbalance on the expression of Arabidopsis imprinted genes FERTILIZATION INDEPENDENT SEED2 (FIS2) and FLOWERING WAGENINGEN (FWA) controlled by DNA methylation, and MEDEA (MEA) and PHERES1 (PHE1) controlled by histone methylation. Genome dosage imbalance deregulated the expression of FIS2 and PHE1 in an antagonistic manner. In addition increased dosage of inactive alleles caused a loss of imprinting of FIS2 and MEA. Although FIS2 controls histone methylation, which represses MEA and PHE1 expression, the changes of PHE1 and MEA expression could not be fully accounted for by the corresponding fluctuations of FIS2 expression. Our results show that parental genome dosage imbalance deregulates imprinting using mechanisms, which are independent from known regulators of imprinting. The complexity of the network of regulations between expressed and silenced alleles of imprinted genes activated in response to parental dosage imbalance does not support simple models derived from the parental conflict hypothesis