The genome and population genomics of allopolyploid Coffea arabica reveal the diversification history of modern coffee cultivars.

Coffea arabica, an allotetraploid hybrid of Coffea eugenioides and Coffea canephora, is the source of approximately 60% of coffee products worldwide, and its cultivated accessions have undergone several population bottlenecks. We present chromosome-level assemblies of a di-haploid C. arabica accession and modern representatives of its diploid progenitors, C. eugenioides and C. canephora. The three species exhibit largely conserved genome structures between diploid parents and descendant subgenomes, with no obvious global subgenome dominance. We find evidence for a founding polyploidy event 350,000-610,000 years ago, followed by several pre-domestication bottlenecks, resulting in narrow genetic variation. A split between wild accessions and cultivar progenitors occurred ~30.5 thousand years ago, followed by a period of migration between the two populations. Analysis of modern varieties, including lines historically introgressed with C. canephora, highlights their breeding histories and loci that may contribute to pathogen resistance, laying the groundwork for future genomics-based breeding of C. arabica.

Ecological and genomic vulnerability to climate change across native populations of Robusta coffee (Coffea canephora).

The assessment of population vulnerability under climate change is crucial for planning conservation as well as for ensuring food security. Coffea canephora is, in its native habitat, an understorey tree that is mainly distributed in the lowland rainforests of tropical Africa. Also known as Robusta, its commercial value constitutes a significant revenue for many human populations in tropical countries. Comparing ecological and genomic vulnerabilities within the species' native range can provide valuable insights about habitat loss and the species' adaptive potential, allowing to identify genotypes that may act as a resource for varietal improvement. By applying species distribution models, we assessed ecological vulnerability as the decrease in climatic suitability under future climatic conditions from 492 occurrences. We then quantified genomic vulnerability (or risk of maladaptation) as the allelic composition change required to keep pace with predicted climate change. Genomic vulnerability was estimated from genomic environmental correlations throughout the native range. Suitable habitat was predicted to diminish to half its size by 2050, with populations near coastlines and around the Congo River being the most vulnerable. Whole-genome sequencing revealed 165 candidate SNPs associated with climatic adaptation in C. canephora, which were located in genes involved in plant response to biotic and abiotic stressors. Genomic vulnerability was higher for populations in West Africa and in the region at the border between DRC and Uganda. Despite an overall low correlation between genomic and ecological vulnerability at broad scale, these two components of vulnerability overlap spatially in ways that may become damaging. Genomic vulnerability was estimated to be 23% higher in populations where habitat will be lost in 2050 compared to regions where habitat will remain suitable. These results highlight how ecological and genomic vulnerabilities are relevant when planning on how to cope with climate change regarding an economically important species.

DNA methylation dynamics during early plant life.

BACKGROUND: Cytosine methylation is crucial for gene regulation and silencing of transposable elements in mammals and plants. While this epigenetic mark is extensively reprogrammed in the germline and early embryos of mammals, the extent to which DNA methylation is reset between generations in plants remains largely unknown. RESULTS: Using Arabidopsis as a model, we uncovered distinct DNA methylation dynamics over transposable element sequences during the early stages of plant development. Specifically, transposable elements and their relics show invariably high methylation at CG sites but increasing methylation at CHG and CHH sites. This non-CG methylation culminates in mature embryos, where it reaches saturation for a large fraction of methylated CHH sites, compared to the typical 10-20% methylation level observed in seedlings or adult plants. Moreover, the increase in CHH methylation during embryogenesis matches the hypomethylated state in the early endosperm. Finally, we show that interfering with the embryo-to-seedling transition results in the persistence of high CHH methylation levels after germination, specifically over sequences that are targeted by the RNA-directed DNA methylation (RdDM) machinery. CONCLUSION: Our findings indicate the absence of extensive resetting of DNA methylation patterns during early plant life and point instead to an important role of RdDM in reinforcing DNA methylation of transposable element sequences in every cell of the mature embryo. Furthermore, we provide evidence that this elevated RdDM activity is a specific property of embryogenesis.

MitoRS, a method for high throughput, sensitive, and accurate detection of mitochondrial DNA heteroplasmy.

BACKGROUND: Mitochondrial dysfunction is linked to numerous pathological states, in particular related to metabolism, brain health and ageing. Nuclear encoded gene polymorphisms implicated in mitochondrial functions can be analyzed in the context of classical genome wide association studies. By contrast, mitochondrial DNA (mtDNA) variants are more challenging to identify and analyze for several reasons. First, contrary to the diploid nuclear genome, each cell carries several hundred copies of the circular mitochondrial genome. Mutations can therefore be present in only a subset of the mtDNA molecules, resulting in a heterogeneous pool of mtDNA, a situation referred to as heteroplasmy. Consequently, detection and quantification of variants requires extremely accurate tools, especially when this proportion is small. Additionally, the mitochondrial genome has pseudogenized into numerous copies within the nuclear genome over the course of evolution. These nuclear pseudogenes, named NUMTs, must be distinguished from genuine mtDNA sequences and excluded from the analysis. RESULTS: Here we describe a novel method, named MitoRS, in which the entire mitochondrial genome is amplified in a single reaction using rolling circle amplification. This approach is easier to setup and of higher throughput when compared to classical PCR amplification. Sequencing libraries are generated at high throughput exploiting a tagmentation-based method. Fine-tuned parameters are finally applied in the analysis to allow detection of variants even of low frequency heteroplasmy. The method was thoroughly benchmarked in a set of experiments designed to demonstrate its robustness, accuracy and sensitivity. The MitoRS method requires 5 ng total DNA as starting material. More than 96 samples can be processed in less than a day of laboratory work and sequenced in a single lane of an Illumina HiSeq flow cell. The lower limit for accurate quantification of single nucleotide variants has been measured at 1% frequency. CONCLUSIONS: The MitoRS method enables the robust, accurate, and sensitive analysis of a large number of samples. Because it is cost effective and simple to setup, we anticipate this method will promote the analysis of mtDNA variants in large cohorts, and may help assessing the impact of mtDNA heteroplasmy on metabolic health, brain function, cancer progression, or ageing.

Genome Sequence of Lactobacillus fermentum Strain NCC2970 (CNCM I-5068).

Lactobacillus fermentum NCC2970 (CNCM I-5068) is a lactic acid bacterium originating from the Nestle Culture Collection. Here, we disclose its full 1.9-Gb genome sequence comprising one chromosome with no plasmid.

DNA DAMAGE BINDING PROTEIN2 Shapes the DNA Methylation Landscape.

In eukaryotes, DNA repair pathways help to maintain genome integrity and epigenomic patterns. However, the factors at the nexus of DNA repair and chromatin modification/remodeling remain poorly characterized. Here, we uncover a previously unrecognized interplay between the DNA repair factor DNA DAMAGE BINDING PROTEIN2 (DDB2) and the DNA methylation machinery in Arabidopsis thaliana Loss-of-function mutation in DDB2 leads to genome-wide DNA methylation alterations. Genetic and biochemical evidence indicate that at many repeat loci, DDB2 influences de novo DNA methylation by interacting with ARGONAUTE4 and by controlling the local abundance of 24-nucleotide short interfering RNAs (siRNAs). We also show that DDB2 regulates active DNA demethylation mediated by REPRESSOR OF SILENCING1 and DEMETER LIKE3. Together, these findings reveal a role for the DNA repair factor DDB2 in shaping the Arabidopsis DNA methylation landscape in the absence of applied genotoxic stress.

Oral Phage Therapy of Acute Bacterial Diarrhea With Two Coliphage Preparations: A Randomized Trial in Children From Bangladesh.

BACKGROUND: Antibiotic resistance is rising in important bacterial pathogens. Phage therapy (PT), the use of bacterial viruses infecting the pathogen in a species-specific way, is a potential alternative. METHOD: T4-like coliphages or a commercial Russian coliphage product or placebo was orally given over 4 days to Bangladeshi children hospitalized with acute bacterial diarrhea. Safety of oral phage was assessed clinically and by functional tests; coliphage and Escherichia coli titers and enteropathogens were determined in stool and quantitative diarrhea parameters (stool output, stool frequency) were measured. Stool microbiota was studied by 16S rRNA gene sequencing; the genomes of four fecal Streptococcus isolates were sequenced. FINDINGS: No adverse events attributable to oral phage application were observed (primary safety outcome). Fecal coliphage was increased in treated over control children, but the titers did not show substantial intestinal phage replication (secondary microbiology outcome). 60% of the children suffered from a microbiologically proven E. coli diarrhea; the most frequent diagnosis was ETEC infections. Bacterial co-pathogens were also detected. Half of the patients contained phage-susceptible E. coli colonies in the stool. E. coli represented less than 5% of fecal bacteria. Stool ETEC titers showed only a short-lived peak and were otherwise close to the replication threshold determined for T4 phage in vitro. An interim analysis after the enrollment of 120 patients showed no amelioration in quantitative diarrhea parameter by PT over standard care (tertiary clinical outcome). Stool microbiota was characterized by an overgrowth with Streptococcus belonging to the Streptococcus gallolyticus and Streptococcus salivarius species groups, their abundance correlated with quantitative diarrhea outcome, but genome sequencing did not identify virulence genes. INTERPRETATION: Oral coliphages showed a safe gut transit in children, but failed to achieve intestinal amplification and to improve diarrhea outcome, possibly due to insufficient phage coverage and too low E. coli pathogen titers requiring higher oral phage doses. More knowledge is needed on in vivo phage-bacterium interaction and the role of E. coli in childhood diarrhea for successful PT. FUNDING: The study was supported by a grant from Nestlé Nutrition and Nestlé Health Science. The trial was registered with Identifier NCT00937274 at ClinicalTrials.gov.

Fully Closed Genome Sequences of Five Type Strains of the Genus Cronobacter and One Cronobacter sakazakii Strain.

Cronobacteris associated with infant infections and the consumption of reconstituted infant formula. Here we sequenced and closed six genomes ofC. condimenti(T),C. muytjensii(T),C. universalis(T),C. malonaticus(T),C. dublinensis(T), andC. sakazakiithat can be used as reference genomes in single nucleotide polymorphism (SNP)-based next-generation sequencing (NGS) analysis for source tracking investigations.

Genome-wide negative feedback drives transgenerational DNA methylation dynamics in Arabidopsis.

Epigenetic variations of phenotypes, especially those associated with DNA methylation, are often inherited over multiple generations in plants. The active and inactive chromatin states are heritable and can be maintained or even be amplified by positive feedback in a transgenerational manner. However, mechanisms controlling the transgenerational DNA methylation dynamics are largely unknown. As an approach to understand the transgenerational dynamics, we examined long-term effect of impaired DNA methylation in Arabidopsis mutants of the chromatin remodeler gene DDM1 (Decrease in DNA Methylation 1) through whole genome DNA methylation sequencing. The ddm1 mutation induces a drastic decrease in DNA methylation of transposable elements (TEs) and repeats in the initial generation, while also inducing ectopic DNA methylation at hundreds of loci. Unexpectedly, this ectopic methylation can only be seen after repeated self-pollination. The ectopic cytosine methylation is found primarily in the non-CG context and starts from 3' regions within transcription units and spreads upstream. Remarkably, when chromosomes with reduced DNA methylation were introduced from a ddm1 mutant into a DDM1 wild-type background, the ddm1-derived chromosomes also induced analogous de novo accumulation of DNA methylation in trans. These results lead us to propose a model to explain the transgenerational DNA methylation redistribution by genome-wide negative feedback. The global negative feedback, together with local positive feedback, would ensure robust and balanced differentiation of chromatin states within the genome.