Rapid fine mapping of causative mutations from sets of unordered, contig-sized fragments of genome sequence

Background Traditional Map based Cloning approaches, used for the identification of desirable alleles, are extremely labour intensive and years can elapse between the mutagenesis and the detection of the polymorphism. High throughput sequencing based Mapping-by-sequencing approach requires an ordered genome assembly and cannot be used with fragmented, un-scaffolded draft genomes, limiting its application to model species and precluding many important organisms. Results We addressed this gap in resource and presented a computational method and software implementations called CHERIPIC (Computing Homozygosity Enriched Regions In genomes to Prioritise Identification of Candidate variants). We have successfully validated implementation of CHERIPIC using three different types of bulk segregant sequence data from Arabidopsis, maize and barley, respectively. Conclusions CHERIPIC allows users to rapidly analyse bulk segregant sequence data and we have made it available as a pre-packaged binary with all dependencies for Linux and MacOS and as Galaxy tool

SEPATH: benchmarking the search for pathogens in human tissue whole genome sequence data leads to template pipelines

Background : Human tissue is increasingly being whole genome sequenced as we transition into an era of genomic medicine. With this arises the potential to detect sequences originating from microorganisms, including pathogens amid the plethora of human sequencing reads. In cancer research, the tumorigenic ability of pathogens is being recognized, for example Helicobacter pylori and human papillomavirus in the cases of gastric non-cardia and cervical carcinomas respectively. As of yet, no benchmark has been carried out on the performance of computational approaches for bacterial and viral detection within host-dominated sequence data. Results : We present the results of benchmarking over 70 distinct combinations of tools and parameters on 100 simulated cancer datasets spiked with realistic proportions of bacteria. mOTUs2 and Kraken are the highest performing individual tools achieving median genus level F1-scores of 0.90 and 0.91 respectively. mOTUs2 demonstrates a high performance in estimating bacterial proportions. Employing Kraken on unassembled sequencing reads produces a good but variable performance depending on post-classification filtering parameters. These approaches are investigated on a selection of cervical and gastric cancer whole genome sequences where Alphapapillomavirus and Helicobacter are detected in addition to a variety of other interesting genera. Conclusions : We provide the top performing pipelines from this benchmark in a unifying tool called SEPATH, which is amenable to high throughput sequencing studies across a range of high-performance computing clusters. SEPATH provides a benchmarked and convenient approach to detect pathogens in tissue sequence data helping to determine the relationship between metagenomics and disease

A Downy Mildew Effector Attenuates Salicylic Acid–Triggered Immunity in Arabidopsis by Interacting with the Host Mediator Complex

Plants are continually exposed to pathogen attack but usually remain healthy because they can activate defences upon perception of microbes. However, pathogens have evolved to overcome plant immunity by delivering effectors into the plant cell to attenuate defence, resulting in disease. Recent studies suggest that some effectors may manipulate host transcription, but the specific mechanisms by which such effectors promote susceptibility remain unclear. We study the oomycete downy mildew pathogen of Arabidopsis, Hyaloperonospora arabidopsidis (Hpa), and show here that the nuclear-localized effector HaRxL44 interacts with Mediator subunit 19a (MED19a), resulting in the degradation of MED19a in a proteasome-dependent manner. The Mediator complex of ∼25 proteins is broadly conserved in eukaryotes and mediates the interaction between transcriptional regulators and RNA polymerase II. We found MED19a to be a positive regulator of immunity against Hpa. Expression profiling experiments reveal transcriptional changes resembling jasmonic acid/ethylene (JA/ET) signalling in the presence of HaRxL44, and also 3 d after infection with Hpa. Elevated JA/ET signalling is associated with a decrease in salicylic acid (SA)-triggered immunity (SATI) in Arabidopsis plants expressing HaRxL44 and in med19a loss-of-function mutants, whereas SATI is elevated in plants overexpressing MED19a. Using a PR1::GUS reporter, we discovered that Hpa suppresses PR1 expression specifically in cells containing haustoria, into which RxLR effectors are delivered, but not in nonhaustoriated adjacent cells, which show high PR1::GUS expression levels. Thus, HaRxL44 interferes with Mediator function by degrading MED19, shifting the balance of defence transcription from SA-responsive defence to JA/ET-signalling, and enhancing susceptibility to biotrophs by attenuating SA-dependent gene expression

Arabidopsis downy mildew effector HaRxL106 suppresses plant immunity by binding to RADICAL‐INDUCED CELL DEATH1

Summary - The oomycete pathogen Hyaloperonospora arabidopsidis (Hpa) causes downy mildew disease on Arabidopsis. To colonize its host, Hpa translocates effector proteins that suppress plant immunity into infected host cells. Here, we investigate the relevance of the interaction between one of these effectors, HaRxL106, and Arabidopsis RADICAL‐INDUCED CELL DEATH1 (RCD1). - We use pathogen infection assays as well as molecular and biochemical analyses to test the hypothesis that HaRxL106 manipulates RCD1 to attenuate transcriptional activation of defense genes. - We report that HaRxL106 suppresses transcriptional activation of salicylic acid (SA)‐induced defense genes and alters plant growth responses to light. HaRxL106‐mediated suppression of immunity is abolished in RCD1 loss‐of‐function mutants. We report that RCD1‐type proteins are phosphorylated, and we identified Mut9‐like kinases (MLKs), which function as phosphoregulatory nodes at the level of photoreceptors, as RCD1‐interacting proteins. An mlk1,3,4 triple mutant exhibits stronger SA‐induced defense marker gene expression compared with wild‐type plants, suggesting that MLKs also affect transcriptional regulation of SA signaling. - Based on the combined evidence, we hypothesize that nuclear RCD1/MLK complexes act as signaling nodes that integrate information from environmental cues an

Arabidopsis downy mildew effector HaRxL106 suppresses plant immunity by binding to RADICAL‐INDUCED CELL DEATH1

Summary - The oomycete pathogen Hyaloperonospora arabidopsidis (Hpa) causes downy mildew disease on Arabidopsis. To colonize its host, Hpa translocates effector proteins that suppress plant immunity into infected host cells. Here, we investigate the relevance of the interaction between one of these effectors, HaRxL106, and Arabidopsis RADICAL‐INDUCED CELL DEATH1 (RCD1). - We use pathogen infection assays as well as molecular and biochemical analyses to test the hypothesis that HaRxL106 manipulates RCD1 to attenuate transcriptional activation of defense genes. - We report that HaRxL106 suppresses transcriptional activation of salicylic acid (SA)‐induced defense genes and alters plant growth responses to light. HaRxL106‐mediated suppression of immunity is abolished in RCD1 loss‐of‐function mutants. We report that RCD1‐type proteins are phosphorylated, and we identified Mut9‐like kinases (MLKs), which function as phosphoregulatory nodes at the level of photoreceptors, as RCD1‐interacting proteins. An mlk1,3,4 triple mutant exhibits stronger SA‐induced defense marker gene expression compared with wild‐type plants, suggesting that MLKs also affect transcriptional regulation of SA signaling. - Based on the combined evidence, we hypothesize that nuclear RCD1/MLK complexes act as signaling nodes that integrate information from environmental cues an

EXPRSS: an Illumina based high-throughput expression-profiling method to reveal transcriptional dynamics

Background: Next Generation Sequencing technologies have facilitated differential gene expression analysis through RNA-seq and Tag-seq methods. RNA-seq has biases associated with transcript lengths, lacks uniform coverage of regions in mRNA and requires 10–20 times more reads than a typical Tag-seq. Most existing Tag-seq methods either have biases or not high throughput due to use of restriction enzymes or enzymatic manipulation of 5’ ends of mRNA or use of RNA ligations.  Results: We have developed EXpression Profiling through Randomly Sheared cDNA tag Sequencing (EXPRSS) that employs acoustic waves to randomly shear cDNA and generate sequence tags at a relatively defined position (~150-200 bp) from the 3′ end of each mRNA. Implementation of the method was verified through comparative analysis of expression data generated from EXPRSS, NlaIII-DGE and Affymetrix microarray and through qPCR quantification of selected genes. EXPRSS is a strand specific and restriction enzyme independent tag sequencing method that does not require cDNA length-based data transformations. EXPRSS is highly reproducible, is high-throughput and it also reveals alternative polyadenylation and polyadenylated antisense transcripts. It is cost-effective using barcoded multiplexing, avoids the biases of existing SAGE and derivative methods and can reveal polyadenylation position from paired-end sequencing.  Conclusions: EXPRSS Tag-seq provides sensitive and reliable gene expression data and enables high-throughput expression profiling with relatively simple downstream analysis

Probing formation of cargo/importin-α transport complexes in plant cells using a pathogen effector

Importin-αs are essential adapter proteins that recruit cytoplasmic proteins destined for active nuclear import to the nuclear transport machinery. Cargo proteins interact with the importin-α armadillo repeat domain via nuclear localization sequences (NLSs), short amino acids motifs enriched in Lys and Arg residues. Plant genomes typically encode several importin-α paralogs that can have both specific and partially redundant functions. Although some cargos are preferentially imported by a distinct importin-α it remains unknown how this specificity is generated and to what extent cargos compete for binding to nuclear transport receptors. Here we report that the effector protein HaRxL106 from the oomycete pathogen Hyaloperonospora arabidopsidis co-opts the host cell's nuclear import machinery. We use HaRxL106 as a probe to determine redundant and specific functions of importin-α paralogs from Arabidopsis thaliana. A crystal structure of the importin-α3/MOS6 armadillo repeat domain suggests that five of the six Arabidopsis importin-αs expressed in rosette leaves have an almost identical NLS-binding site. Comparison of the importin-α binding affinities of HaRxL106 and other cargos in vitro and in plant cells suggests that relatively small affinity differences in vitro affect the rate of transport complex formation in vivo. Our results suggest that cargo affinity for importin-α, sequence variation at the importin-α NLS-binding sites and tissue-specific expression levels of importin-αs determine formation of cargo/importin-α transport complexes in plant cells

Potential for diagnosis of infectious disease from the 100,000 Genomes Project Metagenomic Dataset: Recommendations for reporting results

The identification of microbiological infection is usually a diagnostic investigation, a complex process that is firstly initiated by clinical suspicion. With the emergence of high-throughput sequencing (HTS) technologies, metagenomic analysis has unveiled the power to identify microbial DNA/RNA from a diverse range of clinical samples (1). Metagenomic analysis of whole human genomes at the clinical/research interface bypasses the steps of clinical scrutiny and targeted testing and has the potential to generate unexpected findings relating to infectious and sometimes transmissible disease. There is no doubt that microbial findings that may have a significant impact on a patient’s treatment and their close contacts should be reported to those with clinical responsibility for the sample-donating patient. There are no clear recommendations on how such findings that are incidental, or outside the original investigation, should be handled. Here we aim to provide an informed protocol for the management of incidental microbial findings as part of the 100,000 Genomes Projectwhich may have broader application in this emerging field. As with any other clinical information, we aim to prioritise the reporting of data that are most likely to be of benefit to the patient and their close contacts. We also set out to minimize risks, costs and potential anxiety associated with the reporting of results that are unlikely to be of clinical significance. Our recommendations aim to support the practice of microbial metagenomics by providing a simplified pathway that can be applied to reporting the identification of potential pathogens from metagenomic datasets. Given that the ambition for UK sequenced human genomes over the next 5 years has been set to reach 5 million and the field of metagenomics is rapidly evolving, the guidance will be regularly reviewed and will likely adapt over time as experience develops

The nuclear immune receptor RPS4 is required for RRS1SLH1-dependent constitutive defense activation in Arabidopsis thaliana

Plant nucleotide-binding leucine-rich repeat (NB-LRR) disease resistance (R) proteins recognize specific ‘‘avirulent’’ pathogen effectors and activate immune responses. NB-LRR proteins structurally and functionally resemble mammalian Nod-like receptors (NLRs). How NB-LRR and NLR proteins activate defense is poorly understood. The divergently transcribed Arabidopsis R genes, RPS4 (resistance to Pseudomonas syringae 4) and RRS1 (resistance to Ralstonia solanacearum 1), function together to confer recognition of Pseudomonas AvrRps4 and Ralstonia PopP2. RRS1 is the only known recessive NBLRR R gene and encodes a WRKY DNA binding domain, prompting suggestions that it acts downstream of RPS4 for transcriptional activation of defense genes. We define here the early RRS1-dependent transcriptional changes upon delivery of PopP2 via Pseudomonas type III secretion. The Arabidopsis slh1 (sensitive to low humidity 1) mutant encodes an RRS1 allele (RRS1SLH1) with a single amino acid (leucine) insertion in the WRKY DNA-binding domain. Its poor growth due to constitutive defense activation is rescued at higher temperature. Transcription profiling data indicate that RRS1SLH1-mediated defense activation overlaps substantially with AvrRps4- and PopP2-regulated responses. To better understand the genetic basis of RPS4/RRS1-dependent immunity, we performed a genetic screen to identify suppressor of slh1 immunity (sushi) mutants. We show that many sushi mutants carry mutations in RPS4, suggesting that RPS4 acts downstream or in a complex with RRS1. Interestingly, several mutations were identified in a domain C-terminal to the RPS4 LRR domain. Using an Agrobacterium-mediated transient assay system, we demonstrate that the P-loop motif of RPS4 but not of RRS1SLH1 is required for RRS1SLH1 function. We also recapitulate the dominant suppression of RRS1SLH1 defense activation by wild type RRS1 and show this suppression requires an intact RRS1 P-loop. These analyses of RRS1SLH1 shed new light on mechanisms by which NB-LRR protein pairs activate defense signaling, or are held inactive in the absence of a pathogen effector

Albugo-imposed changes to tryptophan-derived antimicrobial metabolite biosynthesis may contribute to suppression of non-host resistance to Phytophthora infestans in Arabidopsis thaliana.

BACKGROUND: Plants are exposed to diverse pathogens and pests, yet most plants are resistant to most plant pathogens. Non-host resistance describes the ability of all members of a plant species to successfully prevent colonization by any given member of a pathogen species. White blister rust caused by Albugo species can overcome non-host resistance and enable secondary infection and reproduction of usually non-virulent pathogens, including the potato late blight pathogen Phytophthora infestans on Arabidopsis thaliana. However, the molecular basis of host defense suppression in this complex plant-microbe interaction is unclear. Here, we investigate specific defense mechanisms in Arabidopsis that are suppressed by Albugo infection. RESULTS: Gene expression profiling revealed that two species of Albugo upregulate genes associated with tryptophan-derived antimicrobial metabolites in Arabidopsis. Albugo laibachii-infected tissue has altered levels of these metabolites, with lower indol-3-yl methylglucosinolate and higher camalexin accumulation than uninfected tissue. We investigated the contribution of these Albugo-imposed phenotypes to suppression of non-host resistance to P. infestans. Absence of tryptophan-derived antimicrobial compounds enables P. infestans colonization of Arabidopsis, although to a lesser extent than Albugo-infected tissue. A. laibachii also suppresses a subset of genes regulated by salicylic acid; however, salicylic acid plays only a minor role in non-host resistance to P. infestans. CONCLUSIONS: Albugo sp. alter tryptophan-derived metabolites and suppress elements of the responses to salicylic acid in Arabidopsis. Albugo sp. imposed alterations in tryptophan-derived metabolites may play a role in Arabidopsis non-host resistance to P. infestans. Understanding the basis of non-host resistance to pathogens such as P. infestans could assist in development of strategies to elevate food security