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

Collagen fleeces do not improve colonic anastomotic strength but increase bowel obstructions in an experimental rat model

To investigate whether a collagen fleece kept in place by fibrin glue might seal off a colorectal anastomosis, provide reinforcement, and subsequently improve anastomotic healing. Wistar rats underwent a 1-cm left-sided colonic resection followed by a 4-suture end-to-end anastomosis. They were then randomly assigned to one of three treatment groups: no additional intervention (control, n = 20), the anastomosis covered with fibrin glue (fibrin glue, n = 20), the anastomosis covered with a collagen fleece, kept in place with fibrin glue (collagen fleece, n = 21). At either 3 or 7 days follow-up, anastomotic bursting pressure was measured and tissue was obtained for histology and collagen content assessment after which animals were sacrificed. Three rats in the control (15%), three in the fibrin glue (15%), and one in the collagen group (4.8%) died due to anastomotic complications (P = 0.497). Anastomotic bursting pressures were not significantly different between groups at 3 and 7 days follow-up (P = 0.659 and P = 0.427, respectively). However, bowel obstructions occurred significantly more often in the collagen group compared to the control group (14/21 vs. 3/20, P = 0.003). Collagen contents were not different between groups, but histology showed a more severe inflammation in the collagen group compared to the other groups at both 3 and 7 days follow-up. A collagen fleece kept in place by fibrin glue does not improve healing of colonic anastomoses in rats. Moreover, this technique induces significantly more bowel obstructions in rats, warranting further study before being translated to a clinical settin

A somatic-mutational process recurrently duplicates germline susceptibility loci and tissue-specific super-enhancers in breast cancers

Somatic rearrangements contribute to the mutagenized landscape of cancer genomes. Here, we systematically interrogated rearrangements in 560 breast cancers by using a piecewise constant fitting approach. We identified 33 hotspots of large (>100 kb) tandem duplications, a mutational signature associated with homologous-recombination-repair deficiency. Notably, these tandem-duplication hotspots were enriched in breast cancer germline susceptibility loci (odds ratio (OR) = 4.28) and breast-specific 'super-enhancer' regulatory elements (OR = 3.54). These hotspots may b

A somatic-mutational process recurrently duplicates germline susceptibility loci and tissue-specific super-enhancers in breast cancers

Somatic rearrangements contribute to the mutagenized landscape of cancer genomes. Here, we systematically interrogated rearrangements in 560 breast cancers by using a piecewise constant fitting approach. We identified 33 hotspots of large (>100 kb) tandem duplications, a mutational signature associated with homologous-recombination-repair deficiency. Notably, these tandem-duplication hotspots were enriched in breast cancer germline susceptibility loci (odds ratio (OR) = 4.28) and breast-specific 'super-enhancer' regulatory elements (OR = 3.54). These hotspots may be sites of selective susceptibility to double-strand-break damage due to high transcriptional activity or, through incrementally increasing copy number, may be sites of secondary selective pressure. The transcriptomic consequences ranged from strong individual oncogene effects to weak but quantifiable multigene expression effects. We thus present a somatic-rearrangement mutational process affecting coding sequences and noncoding regulatory elements and contributing a continuum of driver consequences, from modest to strong effects, thereby supporting a polygenic model of cancer development.DG is supported by the EU-FP7-SUPPRESSTEM project. SN-Z is funded by a Wellcome Trust Intermediate Fellowship (WT100183MA) and is a Wellcome Beit Fellow. For more information, please visit the publisher's website

Defining the importance of landscape metrics for large branchiopod biodiversity and conservation: the case of the Iberian Peninsula and Balearic Islands

The deficiency in the distributional data of invertebrate taxa is one of the major impediments acting on the bias towards the low awareness of its conservation status. The present study sets a basic framework to understand the large branchiopods distribution in the Iberian Peninsula and Balearic Islands. Since the extensive surveys performed in the late 1980s, no more studies existed updating the information for the whole studied area. The present study fills the gap, gathering together all available information on large branchiopods distribution since 1995, and analysing the effect of human population density and several landscape characteristics on their distribution, taking into consideration different spatial scales (100 m, 1 km and 10 km). In overall, 28 large branchiopod taxa (17 anostracans, 7 notostracans and 4 spinicaudatans) are known to occur in the area. Approximately 30% of the sites hosted multiple species, with a maximum of 6 species. Significant positive co-occurring species pairs were found clustered together, forming 4 different associations of large branchiopod species. In general, species clustered in the same group showed similar responses to analysed landscape characteristics, usually showing a better fit at higher spatial scales.Brazilian Conselho Nacional de Desenvolvimento Cientifico e Tecnologico-CNPq [401045/2014-5]Spanish Ministry of Education, Culture and Sport [FPU014/06783]info:eu-repo/semantics/publishedVersio

Silent chromatin at the middle and ends: lessons from yeasts

Eukaryotic centromeres and telomeres are specialized chromosomal regions that share one common characteristic: their underlying DNA sequences are assembled into heritably repressed chromatin. Silent chromatin in budding and fission yeast is composed of fundamentally divergent proteins tat assemble very different chromatin structures. However, the ultimate behaviour of silent chromatin and the pathways that assemble it seem strikingly similar among Saccharomyces cerevisiae (S. cerevisiae), Schizosaccharomyces pombe (S. pombe) and other eukaryotes. Thus, studies in both yeasts have been instrumental in dissecting the mechanisms that establish and maintain silent chromatin in eukaryotes, contributing substantially to our understanding of epigenetic processes. In this review, we discuss current models for the generation of heterochromatic domains at centromeres and telomeres in the two yeast species

Genetically defined elevated homocysteine levels do not result in widespread changes of DNA methylation in leukocytes

BACKGROUND:DNA methylation is affected by the activities of the key enzymes and intermediate metabolites of the one-carbon pathway, one of which involves homocysteine. We investigated the effect of the well-known genetic variant associated with mildly elevated homocysteine: MTHFR 677C>T independently and in combination with other homocysteine-associated variants, on genome-wide leukocyte DNA-methylation. METHODS:Methylation levels were assessed using Illumina 450k arrays on 9,894 individuals of European ancestry from 12 cohort studies. Linear-mixed-models were used to study the association of additive MTHFR 677C>T and genetic-risk score (GRS) based on 18 homocysteine-associated SNPs, with genome-wide methylation. RESULTS:Meta-analysis revealed that the MTHFR 677C>T variant was associated with 35 CpG sites in cis, and the GRS showed association with 113 CpG sites near the homocysteine-associated variants. Genome-wide analysis revealed that the MTHFR 677C>T variant was associated with 1 trans-CpG (nearest gene ZNF184), while the GRS model showed association with 5 significant trans-CpGs annotated to nearest genes PTF1A, MRPL55, CTDSP2, CRYM and FKBP5. CONCLUSIONS:Our results do not show widespread changes in DNA-methylation across the genome, and therefore do not support the hypothesis that mildly elevated homocysteine is associated with widespread methylation changes in leukocytes

HRDetect is a predictor of BRCA1 and BRCA2 deficiency based on mutational signatures.

Approximately 1-5% of breast cancers are attributed to inherited mutations in BRCA1 or BRCA2 and are selectively sensitive to poly(ADP-ribose) polymerase (PARP) inhibitors. In other cancer types, germline and/or somatic mutations in BRCA1 and/or BRCA2 (BRCA1/BRCA2) also confer selective sensitivity to PARP inhibitors. Thus, assays to detect BRCA1/BRCA2-deficient tumors have been sought. Recently, somatic substitution, insertion/deletion and rearrangement patterns, or 'mutational signatures', were associated with BRCA1/BRCA2 dysfunction. Herein we used a lasso logistic regression model to identify six distinguishing mutational signatures predictive of BRCA1/BRCA2 deficiency. A weighted model called HRDetect was developed to accurately detect BRCA1/BRCA2-deficient samples. HRDetect identifies BRCA1/BRCA2-deficient tumors with 98.7% sensitivity (area under the curve (AUC) = 0.98). Application of this model in a cohort of 560 individuals with breast cancer, of whom 22 were known to carry a germline BRCA1 or BRCA2 mutation, allowed us to identify an additional 22 tumors with somatic loss of BRCA1 or BRCA2 and 47 tumors with functional BRCA1/BRCA2 deficiency where no mutation was detected. We validated HRDetect on independent cohorts of breast, ovarian and pancreatic cancers and demonstrated its efficacy in alternative sequencing strategies. Integrating all of the classes of mutational signatures thus reveals a larger proportion of individuals with breast cancer harboring BRCA1/BRCA2 deficiency (up to 22%) than hitherto appreciated (∼1-5%) who could have selective therapeutic sensitivity to PARP inhibition

Processed pseudogenes acquired somatically during cancer development

Cancer evolves by mutation, with somatic reactivation of retrotransposons being one such mutational process. Germline retrotransposition can cause processed pseudogenes, but whether this occurs somatically has not been evaluated. Here we screen sequencing data from 660 cancer samples for somatically acquired pseudogenes. We find 42 events in 17 samples, especially non-small cell lung cancer (5/27) and colorectal cancer (2/11). Genomic features mirror those of germline LINE element retrotranspositions, with frequent target-site duplications (67%), consensus TTTTAA sites at insertion points, inverted rearrangements (21%), 5′ truncation (74%) and polyA tails (88%). Transcriptional consequences include expression of pseudogenes from UTRs or introns of target genes. In addition, a somatic pseudogene that integrated into the promoter and first exon of the tumour suppressor gene, MGA, abrogated expression from that allele. Thus, formation of processed pseudogenes represents a new class of mutation occurring during cancer development, with potentially diverse functional consequences depending on genomic context

Landscape of somatic mutations in 560 breast cancer whole-genome sequences.

We analysed whole-genome sequences of 560 breast cancers to advance understanding of the driver mutations conferring clonal advantage and the mutational processes generating somatic mutations. We found that 93 protein-coding cancer genes carried probable driver mutations. Some non-coding regions exhibited high mutation frequencies, but most have distinctive structural features probably causing elevated mutation rates and do not contain driver mutations. Mutational signature analysis was extended to genome rearrangements and revealed twelve base substitution and six rearrangement signatures. Three rearrangement signatures, characterized by tandem duplications or deletions, appear associated with defective homologous-recombination-based DNA repair: one with deficient BRCA1 function, another with deficient BRCA1 or BRCA2 function, the cause of the third is unknown. This analysis of all classes of somatic mutation across exons, introns and intergenic regions highlights the repertoire of cancer genes and mutational processes operating, and progresses towards a comprehensive account of the somatic genetic basis of breast cancer