Bulked segregant transcriptome analysis in pea identifies key expression markers for resistance to Peyronellaea pinodes.

Peyronellaea pinodes is a devastating pathogen of pea crop. Quantitative trait loci (QTL) associated with resistance have been identified, as well as genes differentially expressed between resistant and susceptible pea lines. The key question is which of these many genes located into these QTLs, or differentially expressed, are the key genes that distinguish resistant from susceptible plants and could be used as markers. To identify these key genes, in the present study we applied MACE (Massive Analysis of cDNA Ends) -Seq to a whole Recombinant Inbred Line population segregating for resistance to this disease and their parental lines and identified those genes which expression was more correlated with the level of resistance. We also compared gene expression profiles between the most resistant and the most susceptible families of the RIL population. A total of 6780 transcripts were differentially expressed between the parental lines after inoculation. Of them, 803 showed the same expression pattern in the bulks formed by the most resistant and most susceptible RIL families. These genes, showing a consistent expression pattern, could be used as expression markers to distinguish resistant from susceptible plants. The analysis of these genes also discovered the crucial mechanisms acting against P. pinodes.

A combined transcriptomic approach to identify candidates for an anti-tick vaccine blocking B. afzelii transmission.

Ixodes ricinus is the vector for Borrelia afzelii, the predominant cause of Lyme borreliosis in Europe, whereas Ixodes scapularis is the vector for Borrelia burgdorferi in the USA. Transcription of several I. scapularis genes changes in the presence of B. burgdorferi and contributes to successful infection. To what extend B. afzelii influences gene expression in I. ricinus salivary glands is largely unknown. Therefore, we measured expression of uninfected vs. infected tick salivary gland genes during tick feeding using Massive Analysis of cDNA Ends (MACE) and RNAseq, quantifying 26.179 unique transcripts. While tick feeding was the main differentiator, B. afzelii infection significantly affected expression of hundreds of transcripts, including 465 transcripts after 24 h of tick feeding. Validation of the top-20 B. afzelii-upregulated transcripts at 24 h of tick feeding in ten biological genetic distinct replicates showed that expression varied extensively. Three transcripts could be validated, a basic tail protein, a lipocalin and an ixodegrin, and might be involved in B. afzelii transmission. However, vaccination with recombinant forms of these proteins only marginally altered B. afzelii infection in I. ricinus-challenged mice for one of the proteins. Collectively, our data show that identification of tick salivary genes upregulated in the presence of pathogens could serve to identify potential pathogen-blocking vaccine candidates.

Epigenetics in non-classical monocytes support their pro-inflammatory gene expression.

Non-classical human monocytes are characterized by high-level expression of cytokines like TNF, but the mechanisms involved are elusive. We have identified miRNAs and CpG-methylation sites that are unique to non-classical monocytes, defined via CD14 and CD16 expression levels. For down-regulated miRNAs that are linked to up-regulated mRNAs the dominant gene ontology term was intracellular signal transduction. This included down-regulated miRNA-20a-5p and miRNA-106b-5p, which both are linked to increased mRNA for the TRIM8 signaling molecule. Methylation analysis revealed 16 hypo-methylated CpG sites upstream of 14 differentially increased mRNAs including 2 sites upstream of TRIM8. Consistent with a positive role in signal transduction, high TRIM8 levels went along with high basal TNF mRNA levels in non-classical monocytes. Since cytokine expression levels in monocytes strongly increase after stimulation with toll-like-receptor ligands, we have analyzed non-classical monocytes (defined via slan expression) after stimulation with lipopolysaccharide (LPS). LPS-stimulated cells continued to have low miRNA-20a and miRNA-106b and high TRIM8 mRNA levels and they showed a 10-fold increase in TNF mRNA. These data suggest that decreased miRNAs and CpG hypo-methylation is linked to enhanced expression of TRIM8 and that this can contribute to the increased TNF levels in non-classical human monocytes.

Identification of Tick Ixodes ricinus Midgut Genes Differentially Expressed During the Transmission of Borrelia afzelii Spirochetes Using a Transcriptomic Approach.

Lyme borreliosis is an emerging tick-borne disease caused by spirochetes Borrelia burgdorferi sensu lato. In Europe, Lyme borreliosis is predominantly caused by Borrelia afzelii and transmitted by Ixodes ricinus. Although Borrelia behavior throughout tick development is quite well documented, specific molecular interactions between Borrelia and the tick have not been satisfactorily examined. Here, we present the first transcriptomic study focused on the expression of tick midgut genes regulated by Borrelia. By using massive analysis of cDNA ends (MACE), we searched for tick transcripts expressed differentially in the midgut of unfed, 24h-fed, and fully fed I. ricinus nymphs infected with B. afzelii. In total, we identified 553 upregulated and 530 downregulated tick genes and demonstrated that B. afzelii interacts intensively with the tick. Technical and biological validations confirmed the accuracy of the transcriptome. The expression of five validated tick genes was silenced by RNA interference. Silencing of the uncharacterized protein (GXP_Contig_30818) delayed the infection progress and decreased infection prevalence in the target mice tissues. Silencing of other genes did not significantly affect tick feeding nor the transmission of B. afzelii, suggesting a possible role of these genes rather in Borrelia acquisition or persistence in ticks. Identification of genes and proteins exploited by Borrelia during transmission and establishment in a tick could help the development of novel preventive strategies for Lyme borreliosis.

Biofilm Formation of Listeria monocytogenes Strains Under Food Processing Environments and Pan-Genome-Wide Association Study.

Concerns about food contamination by Listeria monocytogenes are on the rise with increasing consumption of ready-to-eat foods. Biofilm production of L. monocytogenes is presumed to be one of the ways that confer its increased resistance and persistence in the food chain. In this study, a collection of isolates from foods and food processing environments (FPEs) representing persistent, prevalent, and rarely detected genotypes was evaluated for biofilm forming capacities including adhesion and sessile biomass production under diverse environmental conditions. The quantity of sessile biomass varied according to growth conditions, lineage, serotype as well as genotype but association of clonal complex (CC) 26 genotype with biofilm production was evidenced under cold temperature. In general, relative biofilm productivity of each strain varied inconsistently across growth conditions. Under our experimental conditions, there were no clear associations between biofilm formation efficiency and persistent or prevalent genotypes. Distinct extrinsic factors affected specific steps of biofilm formation. Sudden nutrient deprivation enhanced cellular adhesion while a prolonged nutrient deficiency impeded biofilm maturation. Salt addition increased biofilm production, moreover, nutrient limitation supplemented by salt significantly stimulated biofilm formation. Pan-genome-wide association study (Pan-GWAS) assessed genetic composition with regard to biofilm phenotypes for the first time. The number of reported genes differed depending on the growth conditions and the number of common genes was low. However, a broad overview of the ontology contents revealed similar patterns regardless of the conditions. Functional analysis showed that functions related to transformation/competence and surface proteins including Internalins were highly enriched.

Gene Expression Profiling and Fine Mapping Identifies a Gibberellin 2-Oxidase Gene Co-segregating With the Dominant Dwarfing Gene Ddw1 in Rye (Secale cereale L.).

The gibberellin (GA)-sensitive dwarfing gene Ddw1 provides an opportunity to genetically reduce plant height in rye. Genetic analysis in a population of recombinant inbred lines confirmed a monogenetic dominant inheritance of Ddw1. Significant phenotypic differences in PH between homo- and heterozygotic genotypes indicate an incomplete dominance of Ddw1. De novo transcriptome sequencing of Ddw1 mutant as well as tall genotypes resulted in 113,547 contigs with an average length of 318 bp covering 36.18 Mbp rye DNA. A hierarchical cluster analysis based on individual groups of rye homologs of functionally characterized rice genes controlling morphological or physiological traits including plant height, flowering time, and source activity identified the gene expression profile of stems at the begin of heading to most comprehensively mirror effects of Ddw1. Genome-wide expression profiling identified 186 transcripts differentially expressed between semi-dwarf and tall genotypes in stems. In total, 29 novel markers have been established and mapped to a 27.2 cM segment in the distal part of the long arm of chromosome 5R. Ddw1 could be mapped within a 0.4 cM interval co-segregating with a marker representing the C20-GA2-oxidase gene ScGA2ox12, that is up-regulated in stems of Ddw1 genotypes. The increased expression of ScGA2ox12 observed in semi-dwarf rye as well as structural alterations in transcript sequences associated with the ScGA2ox12 gene implicate, that Ddw1 is a dominant gain-of-function mutant. Integration of the target interval in the wheat reference genome sequence indicated perfect micro-colinearity between the Ddw1 locus and a 831 kb segment on chromosome 5A, which resides inside of a 11.21 Mb interval carrying the GA-sensitive dwarfing gene Rht12 in wheat. The potential of Ddw1 as a breeder's option to improve lodging tolerance in rye is discussed.

Loss of the Chr16p11.2 ASD candidate gene QPRT leads to aberrant neuronal differentiation in the SH-SY5Y neuronal cell model.

Background: Altered neuronal development is discussed as the underlying pathogenic mechanism of autism spectrum disorders (ASD). Copy number variations of 16p11.2 have recurrently been identified in individuals with ASD. Of the 29 genes within this region, quinolinate phosphoribosyltransferase (QPRT) showed the strongest regulation during neuronal differentiation of SH-SY5Y neuroblastoma cells. We hypothesized a causal relation between this tryptophan metabolism-related enzyme and neuronal differentiation. We thus analyzed the effect of QPRT on the differentiation of SH-SY5Y and specifically focused on neuronal morphology, metabolites of the tryptophan pathway, and the neurodevelopmental transcriptome. Methods: The gene dosage-dependent change of QPRT expression following Chr16p11.2 deletion was investigated in a lymphoblastoid cell line (LCL) of a deletion carrier and compared to his non-carrier parents. Expression of QPRT was tested for correlation with neuromorphology in SH-SY5Y cells. QPRT function was inhibited in SH-SY5Y neuroblastoma cells using (i) siRNA knockdown (KD), (ii) chemical mimicking of loss of QPRT, and (iii) complete CRISPR/Cas9-mediated knock out (KO). QPRT-KD cells underwent morphological analysis. Chemically inhibited and QPRT-KO cells were characterized using viability assays. Additionally, QPRT-KO cells underwent metabolite and whole transcriptome analyses. Genes differentially expressed upon KO of QPRT were tested for enrichment in biological processes and co-regulated gene-networks of the human brain. Results: QPRT expression was reduced in the LCL of the deletion carrier and significantly correlated with the neuritic complexity of SH-SY5Y. The reduction of QPRT altered neuronal morphology of differentiated SH-SY5Y cells. Chemical inhibition as well as complete KO of the gene were lethal upon induction of neuronal differentiation, but not proliferation. The QPRT-associated tryptophan pathway was not affected by KO. At the transcriptome level, genes linked to neurodevelopmental processes and synaptic structures were affected. Differentially regulated genes were enriched for ASD candidates, and co-regulated gene networks were implicated in the development of the dorsolateral prefrontal cortex, the hippocampus, and the amygdala. Conclusions: In this study, QPRT was causally related to in vitro neuronal differentiation of SH-SY5Y cells and affected the regulation of genes and gene networks previously implicated in ASD. Thus, our data suggest that QPRT may play an important role in the pathogenesis of ASD in Chr16p11.2 deletion carriers.

First genetic linkage map of Lathyrus cicera based on RNA sequencing-derived markers: Key tool for genetic mapping of disease resistance.

The Lathyrus cicera transcriptome was analysed in response to rust (Uromyces pisi) infection to develop novel molecular breeding tools with potential for genetic mapping of resistance in this robust orphan legume species. One RNA-seq library each was generated from control and rust-inoculated leaves from two L. cicera genotypes with contrasting quantitative resistance, de novo assembled into contigs and sequence polymorphisms were identified. In toto, 19,224 SNPs differentiate the susceptible from the partially resistant genotype's transcriptome. In addition, we developed and tested 341 expressed E-SSR markers from the contigs, of which 60.7% varied between the two L. cicera genotypes. A first L. cicera linkage map was created using part of the developed markers in a RIL population from the cross of the two genotypes. This map contains 307 markers, covered 724.2 cM and is organised in 7 major and 2 minor linkage groups, with an average mapping interval of 2.4 cM. The genic markers also enabled us to compare their position in L. cicera map with the physical position of the same markers mapped on Medicago truncatula genome, highlighting a high macrosyntenic conservation between both species. This study provides a large new set of genic polymorphic molecular markers with potential for mapping rust resistances. It represents the first step towards genomics-assisted precision breeding in L. cicera.

Functional characterization of the Hyles euphorbiae hawkmoth transcriptome reveals strong expression of phorbol ester detoxification and seasonal cold hardiness genes.

BACKGROUND: The European spurge hawkmoth, Hyles euphorbiae (Lepidoptera, Sphingidae), has been intensively studied as a model organism for insect chemical ecology, cold hardiness and evolution of species delineation. To understand species isolation mechanisms at a molecular level, this study aims at determining genetic factors underlying two adaptive ecological trait candidates, phorbol ester (TPA) detoxification and seasonal cold acclimation. METHOD: A draft transcriptome of H. euphorbiae was generated using Illumina sequencing, providing the first genomic resource for the hawkmoth subfamily Macroglossinae. RNA expression levels in tissues of experimental TPA feeding larvae and cooled pupae was compared to levels in control larvae and pupae using 26 bp RNA sequence tag libraries (DeepSuperSAGE). Differential gene expression was assessed by homology searches of the tags in the transcriptome. RESULTS: In total, 389 and 605 differentially expressed transcripts for detoxification and cold hardiness, respectively, could be identified and annotated with proteins. The majority (22 of 28) of differentially expressed detox transcripts of the four 'drug metabolism' enzyme groups (cytochrome P450 (CYP), carboxylesterases (CES), glutathione S-transferases (GST) and lipases) are up-regulated. Triacylglycerol lipase was significantly over proportionally annotated among up-regulated detox transcripts. We record several up-regulated lipases, GSTe2, two CESs, CYP9A21, CYP6BD6 and CYP9A17 as candidate genes for further H. euphorbiae TPA detoxification analyses. Differential gene expression of the cold acclimation treatment is marked by metabolic depression with enriched Gene Ontology terms among down-regulated transcripts almost exclusively comprising metabolism, aerobic respiration and dissimilative functions. Down-regulated transcripts include energy expensive respiratory proteins like NADH dehydrogenase, cytochrome oxidase and ATP synthase. Gene expression patterns show shifts in carbohydrate metabolism towards cryoprotectant production. The Glycolysis enzymes, G1Pase, A1e, Gpi and an Akr isoform are up-regulated. Glycerol, an osmolyte which lowers the body liquid supercooling point, appears to be the predominant polyol cryoprotectant in H. euphorbiae diapause pupae. Several protein candidates involved in glucose, glycerol, myo-inositol and potentially sorbitol and trehalose synthesis were identified. CONCLUSIONS: A majority of differently expressed transcripts unique for either detoxification or cold hardiness indicates highly specialized functional adaptation which may have evolved from general cell metabolism and stress response.The transcriptome and extracted candidate biomarkers provide a basis for further gene expression studies of physiological processes and adaptive traits in H. euphorbiae.

Deep Super-SAGE transcriptomic analysis of cold acclimation in lentil (Lens culinaris Medik.).

BACKGROUND: Frost is one of the main abiotic stresses limiting plant distribution and crop production. To cope with the stress, plants evolved adaptations known as cold acclimation or chilling tolerance to maximize frost tolerance. Cold acclimation is a progressive acquisition of freezing tolerance by plants subjected to low non-freezing temperatures which subsequently allows them to survive exposure to frost. Lentil is a cool season grain legume that is challenged by winter frost in some areas of its cultivation. RESULTS: To better understand the genetic base of frost tolerance differential gene expression in response to cold acclimation was investigated. Recombinant inbred lines (RILs) from the cross Precoz x WA8649041 were first classified as cold tolerant or cold susceptible according to their response to temperatures between -3 to -15 °C. Then, RILs from both extremes of the response curve were cold acclimated and the leaf transcriptomes of two bulks each of eight frost tolerant and seven cold susceptible RILs were investigated by Deep Super-SAGE transcriptome profiling. Thus, four RNA bulks were analysed: the acclimated susceptible, the acclimated tolerant and the respective controls (non-acclimated susceptible and non-acclimated tolerant). Approximately 16.5 million 26 nucleotide long Super-SAGE tags were sequenced in the four sets (between ~3 and 5.4 millions). In total, 133,077 different unitags, each representing a particular transcript isoform, were identified in these four sets. Tags which showed a significantly different abundance in any of the bulks (fold change ≥4.0 and a significant p-value <0.001) were selected and used to identify the corresponding lentil gene sequence. Three hundred of such lentil sequences were identified. Most of their known homologs coded for glycine-rich, cold and drought-regulated proteins, dormancy-associated proteins, proline-rich proteins (PRPs) and other membrane proteins. These were generally but not exclusively over-expressed in the acclimated tolerant lines. CONCLUSIONS: This set of candidate genes implicated in the response to frost in lentil represents an useful base for deeper and more detailed investigations into this important agronomic trait in future.

Reprint of: MicroRNA profiling of human intermediate monocytes.

Among the three human monocyte subsets, intermediate CD14++CD16+ monocytes have been characterized as particularly proinflammatory cells in experimental studies and as potential biomarkers of cardiovascular risk in clinical cohorts. To further substantiate the distinct role of intermediate monocytes within human monocyte heterogeneity, we assessed subset-specific expression of miRNAs as central epigenetic regulators of gene expression. We hypothesized that intermediate monocytes have a distinct miRNA profile compared to classical and non-classical monocytes. By using small RNA-seq we analyzed 662 miRNAs in the three monocyte subsets. We identified 38 miRNAs that are differentially expressed in intermediate monocytes compared to both classical and non-classical monocytes with a p value of <10-10, of which two miRNAs - miR-6087 (upregulated) and miR-150-5p (downregulated) - differed in their expression more than ten-fold. Pathway analysis of the 38 differentially expressed miRNAs linked intermediate monocytes to distinct biological processes such as gene regulation, cell differentiation, toll-like receptor signaling as well as antigen processing and presentation. Moreover, differentially expressed miRNAs were connected to those genes that we previously identified as markers of intermediate monocytes. In aggregation, we provide first genome-wide miRNA data in the context of monocyte heterogeneity, which substantiate the concept of monocyte trichotomy in human immunity. The identification of miRNAs that are specific for intermediate monocytes may allow to develop strategies, which particularly target this cell population while sparing the other two subsets.

Precision medicine for hepatocelluar carcinoma using molecular pattern diagnostics: results from a preclinical pilot study.

The aim of this study was to design a road map for personalizing cancer therapy in hepatocellular carcinoma (HCC) by using molecular pattern diagnostics. As an exploratory study, we investigated molecular patterns of tissues of two tumors from individual HCC patients, which in previous experiments had shown contrasting reactions to the phase 2 transforming growth factor beta receptor 1 inhibitor galunisertib. Cancer-driving molecular patterns encompass - inter alias - altered transcription profiles and somatic mutations in coding regions differentiating tumors from their respective peritumoral tissues and from each other. Massive analysis of cDNA ends and all-exome sequencing demonstrate a highly divergent transcriptional and mutational landscape, respectively, for the two tumors, that offers potential explanations for the tumors contrasting responses to galunisertib. Molecular pattern diagnostics (MPDs) suggest alternative, individual-tumor-specific therapies, which in both cases deviate from the standard sorafenib treatment and from each other. Suggested personalized therapies use kinase inhibitors and immune-focused drugs as well as low-toxicity natural compounds identified using an advanced bioinformatics routine included in the MPD protocol. The MPD pipeline we describe here for the prediction of suitable drugs for treatment of two contrasting HCCs may serve as a blueprint for the design of therapies for various types of cancer.

A Combined Comparative Transcriptomic, Metabolomic, and Anatomical Analyses of Two Key Domestication Traits: Pod Dehiscence and Seed Dormancy in Pea (Pisum sp.).

The origin of the agriculture was one of the turning points in human history, and a central part of this was the evolution of new plant forms, domesticated crops. Seed dispersal and germination are two key traits which have been selected to facilitate cultivation and harvesting of crops. The objective of this study was to analyze anatomical structure of seed coat and pod, identify metabolic compounds associated with water-impermeable seed coat and differentially expressed genes involved in pea seed dormancy and pod dehiscence. Comparative anatomical, metabolomics, and transcriptomic analyses were carried out on wild dormant, dehiscent Pisum elatius (JI64, VIR320) and cultivated, indehiscent Pisum sativum non-dormant (JI92, Cameor) and recombinant inbred lines (RILs). Considerable differences were found in texture of testa surface, length of macrosclereids, and seed coat thickness. Histochemical and biochemical analyses indicated genotype related variation in composition and heterogeneity of seed coat cell walls within macrosclereids. Liquid chromatography-electrospray ionization/mass spectrometry and Laser desorption/ionization-mass spectrometry of separated seed coats revealed significantly higher contents of proanthocyanidins (dimer and trimer of gallocatechin), quercetin, and myricetin rhamnosides and hydroxylated fatty acids in dormant compared to non-dormant genotypes. Bulk Segregant Analysis coupled to high throughput RNA sequencing resulted in identification of 770 and 148 differentially expressed genes between dormant and non-dormant seeds or dehiscent and indehiscent pods, respectively. The expression of 14 selected dormancy-related genes was studied by qRT-PCR. Of these, expression pattern of four genes: porin (MACE-S082), peroxisomal membrane PEX14-like protein (MACE-S108), 4-coumarate CoA ligase (MACE-S131), and UDP-glucosyl transferase (MACE-S139) was in agreement in all four genotypes with Massive analysis of cDNA Ends (MACE) data. In case of pod dehiscence, the analysis of two candidate genes (SHATTERING and SHATTERPROOF) and three out of 20 MACE identified genes (MACE-P004, MACE-P013, MACE-P015) showed down-expression in dorsal and ventral pod suture of indehiscent genotypes. Moreover, MACE-P015, the homolog of peptidoglycan-binding domain or proline-rich extensin-like protein mapped correctly to predicted Dpo1 locus on PsLGIII. This integrated analysis of the seed coat in wild and cultivated pea provides new insight as well as raises new questions associated with domestication and seed dormancy and pod dehiscence.

MicroRNA profiling of human intermediate monocytes.

Among the three human monocyte subsets, intermediate CD14++CD16+ monocytes have been characterized as particularly proinflammatory cells in experimental studies and as potential biomarkers of cardiovascular risk in clinical cohorts. To further substantiate the distinct role of intermediate monocytes within human monocyte heterogeneity, we assessed subset-specific expression of miRNAs as central epigenetic regulators of gene expression. We hypothesized that intermediate monocytes have a distinct miRNA profile compared to classical and non-classical monocytes. By using small RNA-seq we analyzed 662 miRNAs in the three monocyte subsets. We identified 38 miRNAs that are differentially expressed in intermediate monocytes compared to both classical and non-classical monocytes with a p value of <10-10, of which two miRNAs - miR-6087 (upregulated) and miR-150-5p (downregulated) - differed in their expression more than ten-fold. Pathway analysis of the 38 differentially expressed miRNAs linked intermediate monocytes to distinct biological processes such as gene regulation, cell differentiation, toll-like receptor signaling as well as antigen processing and presentation. Moreover, differentially expressed miRNAs were connected to those genes that we previously identified as markers of intermediate monocytes. In aggregation, we provide first genome-wide miRNA data in the context of monocyte heterogeneity, which substantiate the concept of monocyte trichotomy in human immunity. The identification of miRNAs that are specific for intermediate monocytes may allow to develop strategies, which particularly target this cell population while sparing the other two subsets.

Massive analysis of cDNA ends (MACE) reveals a co-segregating candidate gene for LpPg1 stem rust resistance in perennial ryegrass (Lolium perenne).

KEY MESSAGE: Molecular markers including a potential resistance gene co-segregating with the LpPg1 stem rust resistance locus in perennial ryegrass were identified by massive analysis of cDNA ends (MACE) transcriptome profiling. Stem rust caused by Puccinia graminis subsp. graminicola is a severe fungal disease in the forage crop perennial ryegrass and other grasses. The previously identified LpPg1 locus confers efficient resistance against the pathogen. The aim of this study was to identify candidate genes involved in rust resistance and to use them as a resource for the development of molecular markers for LpPg1. To identify such candidates, bulked segregant analysis was combined with NGS-based massive analysis of cDNA ends (MACE) transcriptome profiling. Total RNA was isolated from bulks of infected and non-infected leaf segments from susceptible and resistant genotypes of a full-sibling mapping population and their respective parental lines and MACE was performed. Bioinformatic analysis detected 330 resistance-specific SNPs in 178 transcripts and 341 transcripts that were exclusively expressed in the resistant bulk. The sequences of many of these transcripts were homologous to genes in distinct regions of chromosomes one and four of the model grass Brachypodium distachyon. Of these, 30 were genetically mapped to a 50.8 cM spanning region surrounding the LpPg1 locus. One candidate NBS-LRR gene co-segregated with the resistance locus. Quantitative analysis of gene expression suggests that LpPg1 mediates an efficient resistance mechanism characterized by early recognition of the pathogen, fast defense signaling and rapid induction of antifungal proteins. We demonstrate here that MACE is a cost-efficient, fast and reliable tool that detects polymorphisms for genetic mapping of candidate resistance genes and simultaneously reveals deep insight into the molecular and genetic base of resistance.

DNA methylation profiling reveals differences in the 3 human monocyte subsets and identifies uremia to induce DNA methylation changes during differentiation.

Human monocytes are a heterogeneous cell population consisting of 3 subsets: classical CD14++CD16-, intermediate CD14++CD16+ and nonclassical CD14+CD16++ monocytes. Via poorly characterized mechanisms, intermediate monocyte counts rise in chronic inflammatory diseases, among which chronic kidney disease is of particular epidemiologic importance. DNA methylation is a central epigenetic feature that controls hematopoiesis. By applying next-generation Methyl-Sequencing we now tested how far the 3 monocyte subsets differ in their DNA methylome and whether uremia induces DNA methylation changes in differentiating monocytes. We found that each monocyte subset displays a unique phenotype with regards to DNA methylation. Genes with differentially methylated promoter regions in intermediate monocytes were linked to distinct immunological processes, which is in line with results from recent gene expression analyses. In vitro, uremia induced dysregulation of DNA methylation in differentiating monocytes, which affected several transcription regulators important for monocyte differentiation (e.g., FLT3, HDAC1, MNT) and led to enhanced generation of intermediate monocytes. As potential mediator, the uremic toxin and methylation inhibitor S-adenosylhomocysteine induced shifts in monocyte subsets in vitro, and associated with monocyte subset counts in vivo. Our data support the concept of monocyte trichotomy and the distinct role of intermediate monocytes in human immunity. The shift in monocyte subsets that occurs in chronic kidney disease, a proinflammatory condition of substantial epidemiological impact, may be induced by accumulation of uremic toxins that mediate epigenetic dysregulation.

In planta Identification of Putative Pathogenicity Factors from the Chickpea Pathogen Ascochyta rabiei by De novo Transcriptome Sequencing Using RNA-Seq and Massive Analysis of cDNA Ends.

The most important foliar diseases in legumes worldwide are ascochyta blights. Up to now, in the Ascochyta-legume pathosystem most studies focused on the identification of resistance genes in the host, while very little is known about the pathogenicity factors of the fungal pathogen. Moreover, available data were often obtained from fungi growing under artificial conditions. Therefore, in this study we aimed at the identification of the pathogenicity factors of Ascochyta rabiei, causing ascochyta blight in chickpea. To identify potential fungal pathogenicity factors, we employed RNA-seq and Massive Analysis of cDNA Ends (MACE) to produce comprehensive expression profiles of A. rabiei genes isolated either from the fungus growing in absence of its host or from fungi infecting chickpea leaves. We further provide a comprehensive de novo assembly of the A. rabiei transcriptome comprising 22,725 contigs with an average length of 1178 bp. Since pathogenicity factors are usually secreted, we predicted the A. rabiei secretome, yielding 550 putatively secreted proteins. MACE identified 596 transcripts that were up-regulated during infection. An analysis of these genes identified a collection of candidate pathogenicity factors and unraveled the pathogen's strategy for infecting its host.

Identification of novel small ncRNAs in pollen of tomato.

BACKGROUND: The unprecedented role of sncRNAs in the regulation of pollen biogenesis on both transcriptional and epigenetic levels has been experimentally proven. However, little is known about their global regulation, especially under stress conditions. We used tomato pollen in order to identify pollen stage-specific sncRNAs and their target mRNAs. We further deployed elevated temperatures to discern stress responsive sncRNAs. For this purpose high throughput sncRNA-sequencing as well as Massive Analysis of cDNA Ends (MACE) were performed for three-replicated sncRNAs libraries derived from tomato tetrad, post-meiotic, and mature pollen under control and heat stress conditions. RESULTS: Using the omiRas analysis pipeline we identified known and predicted novel miRNAs as well as sncRNAs from other classes, responsive or not to heat. Differential expression analysis revealed that post-meiotic and mature pollen react most strongly by regulation of the expression of coding and non-coding genomic regions in response to heat. To gain insight to the function of these miRNAs, we predicted targets and annotated them to Gene Ontology terms. This approach revealed that most of them belong to protein binding, transcription, and Serine/Threonine kinase activity GO categories. Beside miRNAs, we observed differential expression of both tRNAs and snoRNAs in tetrad, post-meiotic, and mature pollen when comparing normal and heat stress conditions. CONCLUSIONS: Thus, we describe a global spectrum of sncRNAs expressed in pollen as well as unveiled those which are regulated at specific time-points during pollen biogenesis. We integrated the small RNAs into the regulatory network of tomato heat stress response in pollen.

Characterization and mapping of LanrBo: a locus conferring anthracnose resistance in narrow-leafed lupin (Lupinus angustifolius L.).

KEY MESSAGE: A novel and highly effective source of anthracnose resistance in narrow-leafed lupin was identified. Resistance was shown to be governed by a single dominant locus. Molecular markers have been developed, which can be used for selecting resistant genotypes in lupin breeding. A screening for anthracnose resistance of a set of plant genetic resources of narrow-leafed lupin (Lupinus angustifolius L.) identified the breeding line Bo7212 as being highly resistant to anthracnose (Colletotrichum lupini). Segregation analysis indicated that the resistance of Bo7212 is inherited by a single dominant locus. The corresponding resistance gene was given the designation LanrBo. Previously published molecular anchor markers allowed us to locate LanrBo on linkage group NLL-11 of narrow-leafed lupin. Using information from RNAseq data obtained with inoculated resistant vs. susceptible lupin entries as well as EST-sequence information from the model genome Lotus japonicus, additional SNP and EST markers linked to LanrBo were derived. A bracket of two LanrBo-flanking markers allows for precise marker-assisted selection of the novel resistance gene in narrow-leafed lupin breeding programs.

Lathyrus sativus transcriptome resistance response to Ascochyta lathyri investigated by deepSuperSAGE analysis.

Lathyrus sativus (grass pea) is a temperate grain legume crop with a great potential for expansion in dry areas or zones that are becoming more drought-prone. It is also recognized as a potential source of resistance to several important diseases in legumes, such as ascochyta blight. Nevertheless, the lack of detailed genomic and/or transcriptomic information hampers further exploitation of grass pea resistance-related genes in precision breeding. To elucidate the pathways differentially regulated during ascochyta-grass pea interaction and to identify resistance candidate genes, we compared the early response of the leaf gene expression profile of a resistant L. sativus genotype to Ascochyta lathyri infection with a non-inoculated control sample from the same genotype employing deepSuperSAGE. This analysis generated 14.387 UniTags of which 95.7% mapped to a reference grass pea/rust interaction transcriptome. From the total mapped UniTags, 738 were significantly differentially expressed between control and inoculated leaves. The results indicate that several gene classes acting in different phases of the plant/pathogen interaction are involved in the L. sativus response to A. lathyri infection. Most notably a clear up-regulation of defense-related genes involved in and/or regulated by the ethylene pathway was observed. There was also evidence of alterations in cell wall metabolism indicated by overexpression of cellulose synthase and lignin biosynthesis genes. This first genome-wide overview of the gene expression profile of the L. sativus response to ascochyta infection delivered a valuable set of candidate resistance genes for future use in precision breeding.