Leah Clissold

Associative transcriptomics of traits in the polyploid crop species Brassica napus

Association genetics can quickly and efficiently delineate regions of the genome that control traits and provide markers to accelerate breeding by marker-assisted selection. But most crops are polyploid, making it difficult to identify the required markers and to assemble a genome sequence to order those markers. To circumvent this difficulty, we developed associative transcriptomics, which uses transcriptome sequencing to identify and score molecular markers representing variation in both gene sequences and gene expression, and correlate this with trait variation. Applying the method in the recently formed tetraploid crop Brassica napus, we identified genomic deletions that underlie two quantitative trait loci for glucosinolate content of seeds. The deleted regions contained orthologs of the transcription factor HAG1 (At5g61420), which controls aliphatic glucosinolate biosynthesis in Arabidopsis thaliana. This approach facilitates the application of association genetics in a broad range of crops, even those with complex genomes.

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.

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.

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.

NextClip: an analysis and read preparation tool for Nextera long mate pair libraries

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

Uncovering hidden variation in polyploid wheat

Significance Pasta and bread wheat are polyploid species that carry multiple copies of each gene. Therefore, loss-of-function mutations in one gene copy are frequently masked by functional copies on other genomes. We sequenced the protein coding regions of 2,735 mutant lines and developed a public database including more than 10 million mutations. Researchers and breeders can search this database online, identify mutations in the different copies of their target gene, and request seeds to study gene function or improve wheat varieties. Mutations are being used to improve the nutritional value of wheat, increase the size of the wheat grains, and generate additional variability in flowering genes to improve wheat adaptation to new and changing environments. 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.

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.

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.

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.