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.

Exposure to 3,3',4,4',5-Pentachlorobiphenyl (PCB126) Causes Widespread DNA Hypomethylation in Adult Zebrafish Testis.

Exposure to environmental toxicants during preconception has been shown to affect offspring health and epigenetic mechanisms such as DNA methylation are hypothesized to be involved in adverse outcomes. However, studies addressing the effects of exposure to environmental toxicants during preconception on epigenetic changes in gametes are limited. The objective of this study is to determine the effect of preconceptional exposure to a dioxin-like polychlorinated biphenyl (3,3',4,4',5-pentachlorobiphenyl [PCB126]) on DNA methylation and gene expression in testis. Adult zebrafish were exposed to 3 and 10 nM PCB126 for 24 h and testis tissue was sampled at 7 days postexposure for histology, DNA methylation, and gene expression profiling. Reduced representation bisulfite sequencing revealed 37 and 92 differentially methylated regions (DMRs) in response to 3 and 10 nM PCB126 exposures, respectively. Among them, 19 DMRs were found to be common between both PCB126 treatment groups. Gene ontology (GO) analysis of DMRs revealed that enrichment of terms such as RNA processing, iron-sulfur cluster assembly, and gluconeogenesis. Gene expression profiling showed differential expression of 40 and 1621 genes in response to 3 and 10 nM PCB126 exposures, respectively. GO analysis of differentially expressed genes revealed enrichment of terms related to xenobiotic metabolism, oxidative stress, and immune function. There is no overlap in the GO terms or individual genes between DNA methylation and RNA sequencing results, but functionally many of the altered pathways have been shown to cause spermatogenic defects.

Evolution of DNA Methylation Across Ecdysozoa.

DNA methylation is a crucial, abundant mechanism of gene regulation in vertebrates. It is less prevalent in many other metazoan organisms and completely absent in some key model species, such as Drosophila melanogaster and Caenorhabditis elegans. We report here a comprehensive study of the presence and absence of DNA methyltransferases (DNMTs) in 138 Ecdysozoa, covering Arthropoda, Nematoda, Priapulida, Onychophora, and Tardigrada. Three of these phyla have not been investigated for the presence of DNA methylation before. We observe that the loss of individual DNMTs independently occurred multiple times across ecdysozoan phyla. We computationally predict the presence of DNA methylation based on CpG rates in coding sequences using an implementation of Gaussian Mixture Modeling, MethMod. Integrating both analysis we predict two previously unknown losses of DNA methylation in Ecdysozoa, one within Chelicerata (Mesostigmata) and one in Tardigrada. In the early-branching Ecdysozoa Priapulus caudatus, we predict the presence of a full set of DNMTs and the presence of DNA methylation. We are therefore showing a very diverse and independent evolution of DNA methylation in different ecdysozoan phyla spanning a phylogenetic range of more than 700 million years.

Resilience to height loss of articular cartilage of osteoarthritic stifle joints of old pigs, compared with healthy cartilage from young pigs in a tribological pin-on-plate exposure, revealing similar friction forces.

INTRODUCTION: We saw a lack of data on the biomechanical behavior of degenerated articular cartilage (OA) compared with that of healthy cartilage, even though the susceptibility to wear and tear of articular cartilage plays a key role in the progression of osteoarthritis (OA). Therefore, we performed a comparison between naturally occurring OA and healthy cartilage from pigs, before and after tribological stress. AIM: The aim of the study was to compare OA-cartilage with healthy cartilage and to analyze the resilience to tribological shear stress, which will be measured as height loss (HL), and to friction forces of the cartilage layers. The findings will be substantiated in macro- and microscopical evaluations before and after tribological exposure. METHODS: We assessed stifle joints of fifteen old and sixteen young pigs from the local abattoir radiologically, macroscopically and histologically to determine possible OA alterations. We put pins from the femoral part of the joints and plates from the corresponding tibial plateaus in a pin-on-plate tribometer under stress for about two hours with about 1108 reciprocating cycles under a pressure of approximately 1 MPa. As a surrogate criterion of wear and tear, the HL was recorded in the tribometer. The heights of the cartilage layers measured before and after the tribological exposure were compared histologically. The condition of the cartilage before and after the tribological exposure was analyzed both macroscopically with an adapted ICRS score and microscopically according to Little et al. (2010). We assessed the friction forces acting between the surfaces of the cartilage pair-specimens. RESULTS: Articular cartilage taken from old pigs showed significant degenerative changes compared to that taken from the young animals. The macroscopic and microscopic scores showed strong alterations of the cartilage after the tribological exposure. There was a noticeable HL of the cartilage specimens after the first 100 to 300 cycles. The HL after tribological exposure was lower in the group of the old animals with 0.52 mm ± 0.23 mm than in the group of the young animals with 0.86 mm ± 0.26 mm (p < 0.0001). The data for the HL was validated by the histological height measurements with 0.50 mm ± 0.82 mm for the old and 0.79 mm ±0.53 mm for the young animals (p = 0.133). The friction forces measured at the cartilage of the old animals were 2.25 N ± 1.15 N and 1.89 N ± 1.45 N of the young animals (p = 0.3225). CONCLUSION: Unlike articular cartilage from young pigs, articular cartilage from old pigs showed OA alterations. Tribological shear stress exposure revealed that OA cartilage showed less HL than healthy articular cartilage. Tribological stress exposure in a pin-on-plate tribometer seemed to be an appropriate way to analyze the mechanical stability of articular cartilage, and the applied protocol could reveal weaknesses of the assessed cartilage tissue. Friction and HL seemed to be independent parameters when degenerated and healthy articular cartilage were assessed under tribological exposure in a pin-on- plate tribometer.

Palladium Intercalated into the Walls of Mesoporous Silica as Robust and Regenerable Catalysts for Hydrodeoxygenation of Phenolic Compounds.

Nanostructured noble-metal catalysts traditionally suffer from sintering under high operating temperatures, leading to durability issues and process limitations. The encapsulation of nanostructured catalysts to prevent loss of activity through thermal sintering, while maintaining accessibility of active sites, remains a great challenge in the catalysis community. Here, we report a robust and regenerable palladium-based catalyst, wherein palladium particles are intercalated into the three-dimensional framework of SBA-15-type mesoporous silica. The encapsulated Pd active sites remain catalytically active as demonstrated in high-temperature/pressure phenol hydrodeoxygenation reactions. The confinement of Pd particles in the walls of SBA-15 prevents particle sintering at high temperatures. Moreover, a partially deactivated catalyst containing intercalated particles is regenerated almost completely even after several reaction cycles. In contrast, Pd particles, which are not encapsulated within the SBA-15 framework, sinter and do not recover prior activity after a regeneration procedure.

Ecological plant epigenetics: Evidence from model and non-model species, and the way forward.

Growing evidence shows that epigenetic mechanisms contribute to complex traits, with implications across many fields of biology. In plant ecology, recent studies have attempted to merge ecological experiments with epigenetic analyses to elucidate the contribution of epigenetics to plant phenotypes, stress responses, adaptation to habitat, and range distributions. While there has been some progress in revealing the role of epigenetics in ecological processes, studies with non-model species have so far been limited to describing broad patterns based on anonymous markers of DNA methylation. In contrast, studies with model species have benefited from powerful genomic resources, which contribute to a more mechanistic understanding but have limited ecological realism. Understanding the significance of epigenetics for plant ecology requires increased transfer of knowledge and methods from model species research to genomes of evolutionarily divergent species, and examination of responses to complex natural environments at a more mechanistic level. This requires transforming genomics tools specifically for studying non-model species, which is challenging given the large and often polyploid genomes of plants. Collaboration among molecular geneticists, ecologists and bioinformaticians promises to enhance our understanding of the mutual links between genome function and ecological processes.

The Helicase Aquarius/EMB-4 Is Required to Overcome Intronic Barriers to Allow Nuclear RNAi Pathways to Heritably Silence Transcription.

Small RNAs play a crucial role in genome defense against transposable elements and guide Argonaute proteins to nascent RNA transcripts to induce co-transcriptional gene silencing. However, the molecular basis of this process remains unknown. Here, we identify the conserved RNA helicase Aquarius/EMB-4 as a direct and essential link between small RNA pathways and the transcriptional machinery in Caenorhabditis elegans. Aquarius physically interacts with the germline Argonaute HRDE-1. Aquarius is required to initiate small-RNA-induced heritable gene silencing. HRDE-1 and Aquarius silence overlapping sets of genes and transposable elements. Surprisingly, removal of introns from a target gene abolishes the requirement for Aquarius, but not HRDE-1, for small RNA-dependent gene silencing. We conclude that Aquarius allows small RNA pathways to compete for access to nascent transcripts undergoing co-transcriptional splicing in order to detect and silence transposable elements. Thus, Aquarius and HRDE-1 act as gatekeepers coordinating gene expression and genome defense.

Recent advances in RNA folding.

In the realm of nucleic acid structures, secondary structure forms a conceptually important intermediate level of description and explains the dominating part of the free energy of structure formation. Secondary structures are well conserved over evolutionary time-scales and for many classes of RNAs evolve slower than the underlying primary sequences. Given the close link between structure and function, secondary structure is routinely used as a basis to explain experimental findings. Recent technological advances, finally, have made it possible to assay secondary structure directly using high throughput methods. From a computational biology point of view, secondary structures have a special role because they can be computed efficiently using exact dynamic programming algorithms. In this contribution we provide a short overview of RNA folding algorithms, recent additions and variations and address methods to align, compare, and cluster RNA structures, followed by a tabular summary of the most important software suites in the fields.

RNA-bioinformatics: Tools, services and databases for the analysis of RNA-based regulation.

The importance of RNA-based regulation is becoming more and more evident. Genome-wide sequencing efforts have shown that the majority of the DNA in eukaryotic genomes is transcribed. Advanced high-throughput techniques like CLIP for the genome-wide detection of RNA-protein interactions have shown that post-transcriptional regulation by RNA-binding proteins matches the complexity of transcriptional regulation. The need for a specialized and integrated analysis of RNA-based data has led to the foundation of the RNA Bioinformatics Center (RBC) within the German Network of Bioinformatics Infrastructure (de.NBI). This paper describes the tools, services and databases provided by the RBC, and shows example applications. Furthermore, we have setup an RNA workbench within the Galaxy framework. For an easy dissemination, we offer a virtualized version of Galaxy (via Galaxy Docker) enabling other groups to use our RNA workbench in a very simple way.

Catalysis Meets Nonthermal Separation for the Production of (Alkyl)phenols and Hydrocarbons from Pyrolysis Oil.

A simple and efficient hydrodeoxygenation strategy is described to selectively generate and separate high-value alkylphenols from pyrolysis bio-oil, produced directly from lignocellulosic biomass. The overall process is efficient and only requires low pressures of hydrogen gas (5 bar). Initially, an investigation using model compounds indicates that MoCx /C is a promising catalyst for targeted hydrodeoxygenation, enabling selective retention of the desired Ar-OH substituents. By applying this procedure to pyrolysis bio-oil, the primary products (phenol/4-alkylphenols and hydrocarbons) are easily separable from each other by short-path column chromatography, serving as potential valuable feedstocks for industry. The strategy requires no prior fractionation of the lignocellulosic biomass, no further synthetic steps, and no input of additional (e.g., petrochemical) platform molecules.

CYP3A5 mediates basal and acquired therapy resistance in different subtypes of pancreatic ductal adenocarcinoma.

Although subtypes of pancreatic ductal adenocarcinoma (PDAC) have been described, this malignancy is clinically still treated as a single disease. Here we present patient-derived models representing the full spectrum of previously identified quasi-mesenchymal (QM-PDA), classical and exocrine-like PDAC subtypes, and identify two markers--HNF1A and KRT81--that enable stratification of tumors into different subtypes by using immunohistochemistry. Individuals with tumors of these subtypes showed substantial differences in overall survival, and their tumors differed in drug sensitivity, with the exocrine-like subtype being resistant to tyrosine kinase inhibitors and paclitaxel. Cytochrome P450 3A5 (CYP3A5) metabolizes these compounds in tumors of the exocrine-like subtype, and pharmacological or short hairpin RNA (shRNA)-mediated CYP3A5 inhibition sensitizes tumor cells to these drugs. Whereas hepatocyte nuclear factor 4, alpha (HNF4A) controls basal expression of CYP3A5, drug-induced CYP3A5 upregulation is mediated by the nuclear receptor NR1I2. CYP3A5 also contributes to acquired drug resistance in QM-PDA and classical PDAC, and it is highly expressed in several additional malignancies. These findings designate CYP3A5 as a predictor of therapy response and as a tumor cell-autonomous detoxification mechanism that must be overcome to prevent drug resistance.

A novel method to analyze 5-hydroxymethylcytosine in CpG sequences using maintenance DNA methyltransferase, DNMT1.

Hydroxymethylcytosine has been shown to be involved in DNA demethylation and gene expression. Although methods to determine the position of hydroxymethylcytosine at single-base resolution have been reported, these methods involve some difficulties. Here, we report a simple method to analyze hydroxymethylcytosine in the CpG sequence utilizing the maintenance DNA methylation activity of DNMT1, which selectively methylates hemi-methylated but not hemi-hydroxymethylated CpG sequences. The method enables monitoring of the dynamics of the hydroxymethylation state of a specific genome site.

Evolution of the unspliced transcriptome.

BACKGROUND: Despite their abundance, unspliced EST data have received little attention as a source of information on non-coding RNAs. Very little is know, therefore, about the genomic distribution of unspliced non-coding transcripts and their relationship with the much better studied regularly spliced products. In particular, their evolution has remained virtually unstudied. RESULTS: We systematically study the evidence on unspliced transcripts available in EST annotation tracks for human and mouse, comprising 104,980 and 66,109 unspliced EST clusters, respectively. Roughly one third of these are located totally inside introns of known genes (TINs) and another third overlaps exonic regions (PINs). Eleven percent are "intergenic", far away from any annotated gene. Direct evidence for the independent transcription of many PINs and TINs is obtained from CAGE tag and chromatin data. We predict more than 2000 3'UTR-associated RNA candidates for each human and mouse. Fifteen to twenty percent of the unspliced EST cluster are conserved between human and mouse. With the exception of TINs, the sequences of unspliced EST clusters evolve significantly slower than genomic background. Furthermore, like spliced lincRNAs, they show highly tissue-specific expression patterns. CONCLUSIONS: Unspliced long non-coding RNAs are an important, rapidly evolving, component of mammalian transcriptomes. Their analysis is complicated by their preferential association with complex transcribed loci that usually also harbor a plethora of spliced transcripts. Unspliced EST data, although typically disregarded in transcriptome analysis, can be used to gain insights into this rarely investigated transcriptome component. The frequently postulated connection between lack of splicing and nuclear retention and the surprising overlap of chromatin-associated transcripts suggests that this class of transcripts might be involved in chromatin organization and possibly other mechanisms of epigenetic control.

LocARNAscan: Incorporating thermodynamic stability in sequence and structure-based RNA homology search.

BACKGROUND: The search for distant homologs has become an import issue in genome annotation. A particular difficulty is posed by divergent homologs that have lost recognizable sequence similarity. This same problem also arises in the recognition of novel members of large classes of RNAs such as snoRNAs or microRNAs that consist of families unrelated by common descent. Current homology search tools for structured RNAs are either based entirely on sequence similarity (such as blast or hmmer) or combine sequence and secondary structure. The most prominent example of the latter class of tools is Infernal. Alternatives are descriptor-based methods. In most practical applications published to-date, however, the information contained in covariance models or manually prescribed search patterns is dominated by sequence information. Here we ask two related questions: (1) Is secondary structure alone informative for homology search and the detection of novel members of RNA classes? (2) To what extent is the thermodynamic propensity of the target sequence to fold into the correct secondary structure helpful for this task? RESULTS: Sequence-structure alignment can be used as an alternative search strategy. In this scenario, the query consists of a base pairing probability matrix, which can be derived either from a single sequence or from a multiple alignment representing a set of known representatives. Sequence information can be optionally added to the query. The target sequence is pre-processed to obtain local base pairing probabilities. As a search engine we devised a semi-global scanning variant of LocARNA's algorithm for sequence-structure alignment. The LocARNAscan tool is optimized for speed and low memory consumption. In benchmarking experiments on artificial data we observe that the inclusion of thermodynamic stability is helpful, albeit only in a regime of extremely low sequence information in the query. We observe, furthermore, that the sensitivity is bounded in particular by the limited accuracy of the predicted local structures of the target sequence. CONCLUSIONS: Although we demonstrate that a purely structure-based homology search is feasible in principle, it is unlikely to outperform tools such as Infernal in most application scenarios, where a substantial amount of sequence information is typically available. The LocARNAscan approach will profit, however, from high throughput methods to determine RNA secondary structure. In transcriptome-wide applications, such methods will provide accurate structure annotations on the target side. AVAILABILITY: Source code of the free software LocARNAscan 1.0 and supplementary data are available at http://www.bioinf.uni-leipzig.de/Software/LocARNAscan.

Protein-coding structured RNAs: A computational survey of conserved RNA secondary structures overlapping coding regions in drosophilids.

Functional RNA elements can be embedded also within exonic sequences coding for functional proteins. While not uncommon in viruses, only a few examples of this type have been described in some detail for eukaryotic genomes. Here we use RNAz and RNAcode, two comparative genomics methods that measure signatures of stabilizing selection acting on RNA secondary structure and peptide sequence, resp., to survey the fruit fly genomes. We estimate that there might be on the order of 1000 loci that are subject to dual selection pressure. The used genome-wide screens also expose the limitations of the currently available methods.

Genomic organization of eukaryotic tRNAs.

BACKGROUND: Surprisingly little is known about the organization and distribution of tRNA genes and tRNA-related sequences on a genome-wide scale. While tRNA gene complements are usually reported in passing as part of genome annotation efforts, and peculiar features such as the tandem arrangements of tRNA gene in Entamoeba histolytica have been described in some detail, systematic comparative studies are rare and mostly restricted to bacteria. We therefore set out to survey the genomic arrangement of tRNA genes and pseudogenes in a wide range of eukaryotes to identify common patterns and taxon-specific peculiarities. RESULTS: In line with previous reports, we find that tRNA complements evolve rapidly and tRNA gene and pseudogene locations are subject to rapid turnover. At phylum level, the distributions of the number of tRNA genes and pseudogenes numbers are very broad, with standard deviations on the order of the mean. Even among closely related species we observe dramatic changes in local organization. For instance, 65% and 87% of the tRNA genes and pseudogenes are located in genomic clusters in zebrafish and stickleback, resp., while such arrangements are relatively rare in the other three sequenced teleost fish genomes. Among basal metazoa, Trichoplax adherens has hardly any duplicated tRNA gene, while the sea anemone Nematostella vectensis boasts more than 17000 tRNA genes and pseudogenes. Dramatic variations are observed even within the eutherian mammals. Higher primates, for instance, have 616 +/- 120 tRNA genes and pseudogenes of which 17% to 36% are arranged in clusters, while the genome of the bushbaby Otolemur garnetti has 45225 tRNA genes and pseudogenes of which only 5.6% appear in clusters. In contrast, the distribution is surprisingly uniform across plant genomes. Consistent with this variability, syntenic conservation of tRNA genes and pseudogenes is also poor in general, with turn-over rates comparable to those of unconstrained sequence elements. Despite this large variation in abundance in Eukarya we observe a significant correlation between the number of tRNA genes, tRNA pseudogenes, and genome size. CONCLUSIONS: The genomic organization of tRNA genes and pseudogenes shows complex lineage-specific patterns characterized by an extensive variability that is in striking contrast to the extreme levels of sequence-conservation of the tRNAs themselves. The comprehensive analysis of the genomic organization of tRNA genes and pseudogenes in Eukarya provides a basis for further studies into the interplay of tRNA gene arrangements and genome organization in general.