Genetic diversity of equine herpesvirus 1 isolated from neurological, abortigenic and respiratory disease outbreaks

Equine herpesvirus 1 (EHV‐1) causes respiratory disease, abortion, neonatal death and neurological disease in equines and is endemic in most countries. The viral factors that influence EHV‐1 disease severity are poorly understood, and this has hampered vaccine development. However, the N752D substitution in the viral DNA polymerase catalytic subunit has been shown statistically to be associated with neurological disease. This has given rise to the term “neuropathic strain,” even though strains lacking the polymorphism have been recovered from cases of neurological disease. To broaden understanding of EHV‐1 diversity in the field, 78 EHV‐1 strains isolated over a period of 35 years were sequenced. The great majority of isolates originated from the United Kingdom and included in the collection were low passage isolates from respiratory, abortigenic and neurological outbreaks. Phylogenetic analysis of regions spanning 80% of the genome showed that up to 13 viral clades have been circulating in the United Kingdom and that most of these are continuing to circulate. Abortion isolates grouped into nine clades, and neurological isolates grouped into five. Most neurological isolates had the N752D substitution, whereas most abortion isolates did not, although three of the neurological isolates from linked outbreaks had a different polymorphism. Finally, bioinformatic analysis suggested that recombination has occurred between EHV‐1 clades, between EHV‐1 and equine herpesvirus 4, and between EHV‐1 and equine herpesvirus 8.

Uncovering hidden variation in polyploid wheat

Comprehensive reverse genetic resources, which have been key to understanding gene function in diploid model organisms, are missing in many polyploid crops. Young polyploid species such as wheat, which was domesticated less than 10,000 y ago, have high levels of sequence identity among subgenomes that mask the effects of recessive alleles. Such redundancy reduces the probability of selection of favorable mutations during natural or human selection, but also allows wheat to tolerate high densities of induced mutations. Here we exploited this property to sequence and catalog more than 10 million mutations in the protein-coding regions of 2,735 mutant lines of tetraploid and hexaploid wheat. We detected, on average, 2,705 and 5,351 mutations per tetraploid and hexaploid line, respectively, which resulted in 35-40 mutations per kb in each population. With these mutation densities, we identified an average of 23-24 missense and truncation alleles per gene, with at least one truncation or deleterious missense mutation in more than 90% of the captured wheat genes per population. This public collection of mutant seed stocks and sequence data enables rapid identification of mutations in the different copies of the wheat genes, which can be combined to uncover previously hidden variation. Polyploidy is a central phenomenon in plant evolution, and many crop species have undergone recent genome duplication events. Therefore, the general strategy and methods developed herein can benefit other polyploid crops.

Genome distribution of differential homoeologue contributions to leaf gene expression in bread wheat

Using a combination of de novo transcriptome assembly, a newly developed 9495-marker transcriptome SNP genetic linkage map and comparative genomics approaches, we developed an ordered set of nonredundant transcripts for each of the subgenomes of hexaploid wheat: A (47 160 unigenes), B (59 663 unigenes) and D (40 588 unigenes). We used these as reference sequences against which to map Illumina mRNA-Seq reads derived from young leaf tissue. Transcript abundance was quantified for each unigene. Using a three-way reciprocal BLAST approach, 15 527 triplet sets of homoeologues (one from each genome) were identified. Differential expression (P < 0.05) was identified for 5248 unigenes, with 2906 represented at greater abundance than their two homoeologues and 2342 represented at lower abundance than their two homoeologues. Analysis of gene ontology terms revealed no biases between homoeologues. There was no evidence of genomewide dominance effects, rather the more highly transcribed individual genes were distributed throughout all three genomes. Transcriptome display tile plot, a visualization approach based on CMYK colour space, was developed and used to assess the genome for regions of skewed homoeologue transcript abundance. Extensive striation was revealed, indicative of many small regions of genome dominance (transcripts of homoeologues from one genome more abundant than the others) and many larger regions of genome repression (transcripts of homoeologues from one genome less abundant than the others).

Simian Immunodeficiency Virus Infection of Chimpanzees (Pan troglodytes) Shares Features of Both Pathogenic and Non-pathogenic Lentiviral Infections

The virus-host relationship in simian immunodeficiency virus (SIV) infected chimpanzees is thought to be different from that found in other SIV infected African primates. However, studies of captive SIVcpz infected chimpanzees are limited. Previously, the natural SIVcpz infection of one chimpanzee, and the experimental infection of six chimpanzees was reported, with limited follow-up. Here, we present a long-term study of these seven animals, with a retrospective re-examination of the early stages of infection. The only clinical signs consistent with AIDS or AIDS associated disease was thrombocytopenia in two cases, associated with the development of anti-platelet antibodies. However, compared to uninfected and HIV-1 infected animals, SIVcpz infected animals had significantly lower levels of peripheral blood CD4+ T-cells. Despite this, levels of T-cell activation in chronic infection were not significantly elevated. In addition, while plasma levels of β2 microglobulin, neopterin and soluble TNF-related apoptosis inducing ligand (sTRAIL) were elevated in acute infection, these markers returned to near-normal levels in chronic infection, reminiscent of immune activation patterns in ‘natural host’ species. Furthermore, plasma soluble CD14 was not elevated in chronic infection. However, examination of the secondary lymphoid environment revealed persistent changes to the lymphoid structure, including follicular hyperplasia in SIVcpz infected animals. In addition, both SIV and HIV-1 infected chimpanzees showed increased levels of deposition of collagen and increased levels of Mx1 expression in the T-cell zones of the lymph node. The outcome of SIVcpz infection of captive chimpanzees therefore shares features of both non-pathogenic and pathogenic lentivirus infections.

The HIV-1/AIDS pandemic is the result of cross-species transmission of simian immunodeficiency virus (SIVcpz) from chimpanzees to humans. Many African primates are infected with SIV, but those studied in captivity generally do not develop disease. However, wild chimpanzees infected with SIVcpz are at increased risk of death and may develop an AIDS-like disease. It has therefore been suggested that the viral features which SIVcpz and HIV-1 share, that differentiate them from other species’ SIV, may be critical in the development of disease in both humans and chimpanzees. Here, we present a long-term follow-up of 7 SIVcpz infected chimpanzees, housed in primate centres in the US and Europe, under similar conditions to other studied models. These animals did not develop an AIDS-like disease, after up to 25 years of infection, and showed features similar to other species where disease rarely develops, such as limited immune activation in the blood. However, they also had significantly reduced CD4+ T-cells and disruption to the secondary lymphoid tissues, normally associated with pathogenic primate lentiviral infections. Thus, while SIVcpz infection of chimpanzees shares features of both pathogenic and non-pathogenic infections, disease has not developed in captivity.

PerR controls oxidative stress defence and aerotolerance but not motility-associated phenotypes of Campylobacter jejuni

The foodborne bacterial pathogen Campylobacter jejuni is an obligate microaerophile that is exposed to atmospheric oxygen during transmission through the food chain. Survival under aerobic conditions requires the concerted control of oxidative stress systems, which in C. jejuni are intimately connected with iron metabolism via the PerR and Fur regulatory proteins. Here, we have characterized the roles of C. jejuni PerR in oxidative stress and motility phenotypes, and its regulon at the level of transcription, protein expression and promoter interactions. Insertional inactivation of perR in the C. jejuni reference strains NCTC 11168, 81-176 and 81116 did not result in any growth deficiencies, but strongly increased survival in atmospheric oxygen conditions, and allowed growth around filter discs infused with up to 30 % H2O2 (8.8 M). Expression of catalase, alkyl hydroperoxide reductase, thioredoxin reductase and the Rrc desulforubrerythrin was increased in the perR mutant, and this was mediated at the transcriptional level as shown by electrophoretic mobility shift assays of the katA, ahpC and trxB promoters using purified PerR. Differential RNA-sequencing analysis of a fur perR mutant allowed the identification of eight previously unknown transcription start sites of genes controlled by Fur and/or PerR. Finally, inactivation of perR in C. jejuni did not result in reduced motility, and did not reduce killing of Galleria melonella wax moth larvae. In conclusion, PerR plays an important role in controlling oxidative stress resistance and aerobic survival of C. jejuni, but this role does not extend into control of motility and associated phenotypes.

Genome-guided investigation of plant natural product biosynthesis

The medicinal plant Madagascar periwinkle, Catharanthus roseus (L.) G. Don, produces hundreds of biologically active monoterpene-derived indole alkaloid (MIA) metabolites and is the sole source of the potent, expensive anti-cancer compounds vinblastine and vincristine. Access to a genome sequence would enable insights into the biochemistry, control, and evolution of genes responsible for MIA biosynthesis. However, generation of a near-complete, scaffolded genome is prohibitive to small research communities due to the expense, time, and expertise required. In this study, we generated a genome assembly for C. roseus that provides a near-comprehensive representation of the genic space that revealed the genomic context of key points within the MIA biosynthetic pathway including physically clustered genes, tandem gene duplication, expression sub-functionalization, and putative neo-functionalization. The genome sequence also facilitated high resolution co-expression analyses that revealed three distinct clusters of co-expression within the components of the MIA pathway. Coordinated biosynthesis of precursors and intermediates throughout the pathway appear to be a feature of vinblastine/vincristine biosynthesis. The C. roseus genome also revealed localization of enzyme-rich genic regions and transporters near known biosynthetic enzymes, highlighting how even a draft genome sequence can empower the study of high-value specialized metabolites.

The control of seed oil polyunsaturate content in the polyploid crop species Brassica napus

Many important plant species have polyploidy in their recent ancestry, complicating inferences about the genetic basis of trait variation. Although the principal locus controlling the proportion of polyunsaturated fatty acids (PUFAs) in seeds of Arabidopsis thaliana is known (fatty acid desaturase 2; FAD2), commercial cultivars of a related crop, oilseed rape (Brassica napus), with very low PUFA content have yet to be developed. We showed that a cultivar of oilseed rape with lower than usual PUFA content has non-functional alleles at three of the four orthologous FAD2 loci. To explore the genetic basis further, we developed an ethyl methanesulphonate mutagenised population, JBnaCAB_E, and used it to identify lines that also carried mutations in the remaining functional copy. This confirmed the hypothesised basis of variation, resulting in an allelic series of mutant lines showing a spectrum of PUFA contents of seed oil. Several lines had PUFA content of ~6 % and oleic acid content of ~84 %, achieving a long-standing industry objective: very high oleic, very low PUFA rapeseed without the use of genetic modification technology. The population contains a high rate of mutations and represents an important resource for research in B. napus.

NextClip: an analysis and read preparation tool for Nextera Long Mate Pair libraries

Summary: Illumina’s recently released Nextera Long Mate Pair (LMP) kit enables production of jumping libraries of up to 12 kb. The LMP libraries are an invaluable resource for carrying out complex assemblies and other downstream bioinformatics analyses such as the characterization of structural variants. However, LMP libraries are intrinsically noisy and to maximize their value, post-sequencing data analysis is required. Standardizing laboratory protocols and the selection of sequenced reads for downstream analysis are non-trivial tasks. NextClip is a tool for analyzing reads from LMP libraries, generating a comprehensive quality report and extracting good quality trimmed and deduplicated reads.

Availability and implementation: Source code, user guide and example data are available from https://github.com/richardmleggett/nextclip/.

Contact:Richard.Leggett@tgac.ac.uk

Supplementary information:Supplementary data are available at Bioinformatics online.

Barley whole exome capture: a tool for genomic research in the genus Hordeum and beyond

Advanced resources for genome-assisted research in barley (Hordeum vulgare) including a whole-genome shotgun assembly and an integrated physical map have recently become available. These have made possible studies that aim to assess genetic diversity or to isolate single genes by whole-genome resequencing and in silico variant detection. However such an approach remains expensive given the 5 Gb size of the barley genome. Targeted sequencing of the mRNA-coding exome reduces barley genomic complexity more than 50-fold, thus dramatically reducing this heavy sequencing and analysis load. We have developed and employed an in-solution hybridization-based sequence capture platform to selectively enrich for a 61.6 megabase coding sequence target that includes predicted genes from the genome assembly of the cultivar Morex as well as publicly available full-length cDNAs and de novo assembled RNA-Seq consensus sequence contigs. The platform provides a highly specific capture with substantial and reproducible enrichment of targeted exons, both for cultivated barley and related species. We show that this exome capture platform provides a clear path towards a broader and deeper understanding of the natural variation residing in the mRNA-coding part of the barley genome and will thus constitute a valuable resource for applications such as mapping-by-sequencing and genetic diversity analyzes.

Sequencing-based variant detection in the polyploid crop oilseed rape

BACKGROUND

The detection and exploitation of genetic variation underpins crop improvement. However, the polyploid nature of the genomes of many of our most important crops represents a barrier, particularly for the analysis of variation within genes. To overcome this, we aimed to develop methodologies based on amplicon sequencing that involve the incorporation of barcoded amplification tags (BATs) into PCR products.

RESULTS

A protocol was developed to tag PCR products with 5' 6-base oligonucleotide barcode extensions before pooling for sequencing library production using standard Illumina adapters. A computational method was developed for the de-convolution of products and the robust detection and scoring of sequence variants. Using this methodology, amplicons targeted to gene sequences were screened across a B. napus mapping population and the resulting allele scoring strings for 24 markers linkage mapped to the expected regions of the genome. Furthermore, using one-dimensional 8-fold pooling, 4608 lines of a B. napus mutation population were screened for induced mutations in a locus-specific amplicon (an orthologue of GL2.b) and mixed product of three co-amplified loci (orthologues of FAD2), identifying 10 and 41 mutants respectively.

CONCLUSIONS

The utilisation of barcode tags to de-convolute pooled PCR products in multiplexed, variation screening via Illumina sequencing provides a cost effective method for SNP genotyping and mutation detection and, potentially, markers for causative changes, even in polyploid species. Combining this approach with existing Illumina multiplexing workflows allows the analysis of thousands of lines cheaply and efficiently in a single sequencing run with minimal library production costs.

High-resolution transcriptional analysis of the regulatory influence of cell-to-cell signalling reveals novel genes that contribute to Xanthomonas phytopathogenesis

The bacterium Xanthomonas campestris is an economically important pathogen of many crop species and a model for the study of bacterial phytopathogenesis. In X. campestris, a regulatory system mediated by the signal molecule DSF controls virulence to plants. The synthesis and recognition of the DSF signal depends upon different Rpf proteins. DSF signal generation requires RpfF whereas signal perception and transduction depends upon a system comprising the sensor RpfC and regulator RpfG. Here we have addressed the action and role of Rpf/DSF signalling in phytopathogenesis by high-resolution transcriptional analysis coupled to functional genomics. We detected transcripts for many genes that were unidentified by previous computational analysis of the genome sequence. Novel transcribed regions included intergenic transcripts predicted as coding or non-coding as well as those that were antisense to coding sequences. In total, mutation of rpfF, rpfG and rpfC led to alteration in transcript levels (more than fourfold) of approximately 480 genes. The regulatory influence of RpfF and RpfC demonstrated considerable overlap. Contrary to expectation, the regulatory influence of RpfC and RpfG had limited overlap, indicating complexities of the Rpf signalling system. Importantly, functional analysis revealed over 160 new virulence factors within the group of Rpf-regulated genes.

A resource of mapped dissociation launch pads for targeted insertional mutagenesis in the Arabidopsis genome

We describe a new resource for targeted insertional mutagenesis in Arabidopsis using a maize (Zea mays) Activator/Dissociation (Ds) two-element system. The two components of the system, T-DNA vectors carrying a Ds launch pad and a stable Activator transposase source, were designed to simplify selection of transposition events and maximize their usefulness. Because Ds elements preferentially transpose to nearby genomic sites, they can be used in targeted mutagenesis of linked genes. To efficiently target all genes throughout the genome, we generated a large population of transgenic Arabidopsis plants containing the Ds launch pad construct, identified lines containing single Ds launch pad inserts, and mapped the positions of Ds launch pads in 89 lines. The integration sites of the Ds launch pads were relatively evenly distributed on all five chromosomes, except for a region of chromosomes 2 and 4 and the centromeric regions. This resource therefore provides access to the majority of the Arabidopsis genome for targeted tagging.