A chromosome conformation capture ordered sequence of the barley genome.

Cereal grasses of the Triticeae tribe have been the major food source in temperate regions since the dawn of agriculture. Their large genomes are characterized by a high content of repetitive elements and large pericentromeric regions that are virtually devoid of meiotic recombination. Here we present a high-quality reference genome assembly for barley (Hordeum vulgare L.). We use chromosome conformation capture mapping to derive the linear order of sequences across the pericentromeric space and to investigate the spatial organization of chromatin in the nucleus at megabase resolution. The composition of genes and repetitive elements differs between distal and proximal regions. Gene family analyses reveal lineage-specific duplications of genes involved in the transport of nutrients to developing seeds and the mobilization of carbohydrates in grains. We demonstrate the importance of the barley reference sequence for breeding by inspecting the genomic partitioning of sequence variation in modern elite germplasm, highlighting regions vulnerable to genetic erosion.

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.

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.

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.

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).

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.

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/.

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.

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.

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.

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.

Dissecting the genome of the polyploid crop oilseed rape by transcriptome sequencing.

Polyploidy complicates genomics-based breeding of many crops, including wheat, potato, cotton, oat and sugarcane. To address this challenge, we sequenced leaf transcriptomes across a mapping population of the polyploid crop oilseed rape (Brassica napus) and representative ancestors of the parents of the population. Analysis of sequence variation and transcript abundance enabled us to construct twin single nucleotide polymorphism linkage maps of B. napus, comprising 23,037 markers. We used these to align the B. napus genome with that of a related species, Arabidopsis thaliana, and to genome sequence assemblies of its progenitor species, Brassica rapa and Brassica oleracea. We also developed methods to detect genome rearrangements and track inheritance of genomic segments, including the outcome of an interspecific cross. By revealing the genetic consequences of breeding, cost-effective, high-resolution dissection of crop genomes by transcriptome sequencing will increase the efficiency of predictive breeding even in the absence of a complete genome sequence.

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.