Mapping QTL associated with resistance to Pseudomonas syringae pv. actinidiae in kiwifruit (Actinidia chinensis var. chinensis).

Pseudomonas syringae pv. actinidiae (Psa) is a bacterial pathogen of kiwifruit. This pathogen causes leaf-spotting, cane dieback, wilting, cankers (lesions), and in severe cases, plant death. Families of diploid A. chinensis seedlings grown in the field show a range of susceptibilities to the disease with up to 100% of seedlings in some families succumbing to Psa. But the effect of selection for field resistance to Psa on the alleles that remain in surviving seedlings has not been assessed. The objective of this work was to analyse, the effect of plant removal from Psa on the allele frequency of an incomplete-factorial-cross population. This population was founded using a range of genotypically distinct diploid A. chinensis var. chinensis parents to make 28 F1 families. However, because of the diversity of these families, low numbers of surviving individuals, and a lack of samples from dead individuals, standard QTL mapping approaches were unlikely to yield good results. Instead, a modified bulk segregant analysis (BSA) overcame these drawbacks while reducing the costs of sampling and sample processing, and the complexity of data analysis. Because the method was modified, part one of this work was used to determine the signal strength required for a QTL to be detected with BSA. Once QTL detection accuracy was known, part two of this work analysed the 28 families from the incomplete-factorial-cross population that had multiple individuals removed due to Psa infection. Each family was assigned to one of eight bulks based on a single parent that contributed to the families. DNA was extracted in bulk by grinding sampled leaf discs together before DNA extraction. Each sample bulk was compared against a bulk made up of WGS data from the parents contributing to the sample bulk. The deviation in allele frequency from the expected allele frequency within surviving populations using the modified BSA method was able to identify 11 QTLs for Psa that were present in at least two analyses. The identification of these Psa resistance QTL will enable marker development to selectively breed for resistance to Psa in future kiwifruit breeding programs.

Differential expression analyses reveal extensive transcriptional plasticity induced by temperature in New Zealand silver trevally (Pseudocaranx georgianus).

Ectotherm species, such as marine fishes, depend on environmental temperature to regulate their vital functions. In finfish aquaculture production, being able to predict physiological responses in growth and other economic traits to temperature is crucial to address challenges inherent in the selection of grow-out locations. This will become an even more significant issue under the various predicted future climate change scenarios. In this study, we used the marine teleost silver trevally (Pseudocaranx georgianus), a species currently being explored as a candidate for aquaculture in New Zealand, as a model to study plasticity in gene expression patterns and growth in response to different temperatures. Using a captive study population, temperature conditions were experimentally manipulated for 1 month to mimic seasonal extremes. Phenotypic differences in growth were measured in 400 individuals, and gene expression patterns of pituitary gland and liver were determined in a subset of 100 individuals. Results showed that growth increased 50% in the warmer compared with the colder condition, suggesting that temperature has a large impact on metabolic activities associated with growth. A total of 265,116,678 single-end RNA sequence reads were aligned to the trevally genome, and 28,416 transcript models were developed (27,887 of these had GenBank accessions, and 17,980 unique gene symbols). Further filtering reduced this set to 8597 gene models. 39 and 238 differentially expressed genes (DEGs) were found in the pituitary gland and the liver, respectively (|log2FC| > 0.26, p-value < 0.05). Of these, 6 DEGs showed a common expression pattern between both tissues, all involved in housekeeping functions. Temperature-modulated growth responses were linked to major pathways affecting metabolism, cell regulation and signalling, previously shown to be important for temperature tolerance in other fish species. An interesting finding of this study was that genes linked to the reproductive system were up-regulated in both tissues in the high treatment, indicating the onset of sexual maturation. Few studies have investigated the thermal plasticity of the gene expression in the main organs of the somatotropic axis simultaneously. Our findings indicate that trevally exhibit substantial growth differences and predictable plastic regulatory responses to different temperature conditions. We identified a set of genes that provide a list of candidates for further investigations for selective breeding objectives and how populations may adapt to increasing temperatures.

Genome-wide analysis reveals the genetic stock structure of hoki (Macruronus novaezelandiae).

The assessment of the genetic structuring of biodiversity is crucial for management and conservation. This is particularly critical for widely distributed and highly mobile deep-water teleosts, such as hoki (Macruronus novaezelandiae). This species is significant to Maori people and supports the largest commercial fishery in New Zealand, but uncertainty about its stock structure presents a challenge for management. Here, we apply a comprehensive genomic analysis to shed light on the demographic structure of this species by (1) assembling the genome, (2) generating a catalogue of genome-wide SNPs to infer the stock structure and (3) identifying regions of the genome under selection. The final genome assembly used short and long reads and is near complete, representing 93.8% of BUSCO genes, and consisting of 566 contigs totalling 501 Mb. Whole-genome re-sequencing of 510 hoki sampled from 14 locations around New Zealand and Australia, at a read depth greater than 10×, produced 227,490 filtered SNPs. Analyses of these SNPs were able to resolve the stock structure of hoki into two genetically and geographically distinct clusters, one including the Australian and the other one all New Zealand locations, indicating genetic exchange between these regions is limited. Location differences within New Zealand samples were much more subtle (global F ST  = 0.0006), and while small and significant differences could be detected, they did not conclusively identify additional substructures. Ten putative adaptive SNPs were detected within the New Zealand samples, but these also did not geographically partition the dataset further. Contemporary and historical N e estimation suggest the current New Zealand population of hoki is large yet declining. Overall, our study provides the first genomic resources for hoki and provides detailed insights into the fine-scale population genetic structure to inform the management of this species.

The genome of New Zealand trevally (Carangidae: Pseudocaranx georgianus) uncovers a XY sex determination locus.

BACKGROUND: The genetic control of sex determination in teleost species is poorly understood. This is partly because of the diversity of mechanisms that determine sex in this large group of vertebrates, including constitutive genes linked to sex chromosomes, polygenic constitutive mechanisms, environmental factors, hermaphroditism, and unisexuality. Here we use a de novo genome assembly of New Zealand silver trevally (Pseudocaranx georgianus) together with sex-specific whole genome sequencing data to detect sexually divergent genomic regions, identify candidate genes and develop molecular makers. RESULTS: The de novo assembly of an unsexed trevally (Trevally_v1) resulted in a final assembly of 579.4 Mb in length, with a N50 of 25.2 Mb. Of the assembled scaffolds, 24 were of chromosome scale, ranging from 11 to 31 Mb in length. A total of 28,416 genes were annotated after 12.8 % of the assembly was masked with repetitive elements. Whole genome re-sequencing of 13 wild sexed trevally (seven males and six females) identified two sexually divergent regions located on two scaffolds, including a 6 kb region at the proximal end of chromosome 21. Blast analyses revealed similarity between one region and the aromatase genes cyp19 (a1a/b) (E-value < 1.00E-25, identity > 78.8 %). Males contained higher numbers of heterozygous variants in both regions, while females showed regions of very low read-depth, indicative of male-specificity of this genomic region. Molecular markers were developed and subsequently tested on 96 histologically-sexed fish (42 males and 54 females). Three markers amplified in absolute correspondence with sex (positive in males, negative in females). CONCLUSIONS: The higher number of heterozygous variants in males combined with the absence of these regions in females support a XY sex-determination model, indicating that the trevally_v1 genome assembly was developed from a male specimen. This sex system contrasts with the ZW sex-determination model documented in closely related carangid species. Our results indicate a sex-determining function of a cyp19a1a-like gene, suggesting the molecular pathway of sex determination is somewhat conserved in this family. The genomic resources developed here will facilitate future comparative work, and enable improved insights into the varied sex determination pathways in teleosts. The sex marker developed in this study will be a valuable resource for aquaculture selective breeding programmes, and for determining sex ratios in wild populations.

A novel hairpin library-based approach to identify NBS-LRR genes required for effector-triggered hypersensitive response in Nicotiana benthamiana.

BACKGROUND: PTI and ETI are the two major defence mechanisms in plants. ETI is triggered by the detection of pathogen effectors, or their activity, in the plant cell and most of the time involves internal receptors known as resistance (R) genes. An increasing number of R genes responsible for recognition of specific effectors have been characterised over the years; however, methods to identify R genes are often challenging and cannot always be translated to crop plants. RESULTS: We present a novel method to identify R genes responsible for the recognition of specific effectors that trigger a hypersensitive response (HR) in Nicotiana benthamiana. This method is based on the genome-wide identification of most of the potential R genes of N. benthamiana and a systematic silencing of these potential R genes in a simple transient expression assay. A hairpin-RNAi library was constructed covering 345 R gene candidates of N. benthamiana. This library was then validated using several previously described R genes. Our approach indeed confirmed that Prf, NRC2a/b and NRC3 are required for the HR that is mediated in N. benthamiana by Pto/avrPto (prf, NRC2a/b and NRC3) and by Cf4/avr4 (NRC2a/b and NRC3). We also confirmed that NRG1, in association with N, is required for the Tobacco Mosaic Virus (TMV)-mediated HR in N. benthamiana. CONCLUSION: We present a novel approach combining bioinformatics, multiple-gene silencing and transient expression assay screening to rapidly identify one-to-one relationships between pathogen effectors and host R genes in N. benthamiana. This approach allowed the identification of previously described R genes responsible for detection of avirulence determinants from Pseudomonas, Cladosporium and TMV, demonstrating that the method could be applied to any effectors/proteins originating from a broad range of plant pathogens that trigger an HR in N. benthamiana. Moreover, with the increasing availability of genome sequences from model and crop plants and pathogens, this approach could be implemented in other plants, accelerating the process of identification and characterization of novel resistance genes.

Microarray analysis of kiwifruit (Actinidia chinensis) bark following challenge by the sucking insect Hemiberlesia lataniae (Hemiptera: Diaspididae).

Both commercial and experimental genotypes of kiwifruit (Actinidia spp.) exhibit large differences in response to insect pests. An understanding of the vine's physiological response to insect feeding and its genetic basis will be important in assisting the development of varieties with acceptable levels of pest resistance. This experiment describes transcriptome changes observed in the bark of kiwifruit 2 and 7 days after the commencement of feeding by the armored scale insect pest, Hemiberlesia lataniae. Using a cDNA microarray consisting of 17,512 unigenes, we measured transcriptome changes and analyzed these into functional ontology categories using MapMan. Results are available in the GEO database GSE73922 and are described fully in Ref. Hill et al. (2015) [1]. After 7 days, transcripts associated with photosynthesis were down-regulated and secondary metabolism was up-regulated. Differential expression of transcripts associated with stress response was consistent with a defense response involving both effector and herbivore-triggered immunities, with predominant involvement of the salicylic acid phytohormonal pathway. This hypothesis was supported by the results of two laboratory experiments. The methods described here could be further adapted and applied to the study of plant responses to a wide range of sessile sucking pests.

Random Tagging Genotyping by Sequencing (rtGBS), an Unbiased Approach to Locate Restriction Enzyme Sites across the Target Genome.

UNLABELLED: Genotyping by sequencing (GBS) is a restriction enzyme based targeted approach developed to reduce the genome complexity and discover genetic markers when a priori sequence information is unavailable. Sufficient coverage at each locus is essential to distinguish heterozygous from homozygous sites accurately. The number of GBS samples able to be pooled in one sequencing lane is limited by the number of restriction sites present in the genome and the read depth required at each site per sample for accurate calling of single-nucleotide polymorphisms. Loci bias was observed using a slight modification of the Elshire et al. METHOD: some restriction enzyme sites were represented in higher proportions while others were poorly represented or absent. This bias could be due to the quality of genomic DNA, the endonuclease and ligase reaction efficiency, the distance between restriction sites, the preferential amplification of small library restriction fragments, or bias towards cluster formation of small amplicons during the sequencing process. To overcome these issues, we have developed a GBS method based on randomly tagging genomic DNA (rtGBS). By randomly landing on the genome, we can, with less bias, find restriction sites that are far apart, and undetected by the standard GBS (stdGBS) method. The study comprises two types of biological replicates: six different kiwifruit plants and two independent DNA extractions per plant; and three types of technical replicates: four samples of each DNA extraction, stdGBS vs. rtGBS methods, and two independent library amplifications, each sequenced in separate lanes. A statistically significant unbiased distribution of restriction fragment size by rtGBS showed that this method targeted 49% (39,145) of BamH I sites shared with the reference genome, compared to only 14% (11,513) by stdGBS.

Transcriptome Analysis of Kiwifruit (Actinidia chinensis) Bark in Response to Armoured Scale Insect (Hemiberlesia lataniae) Feeding.

The kiwifruit cultivar Actinidia chinensis 'Hort16A' is resistant to the polyphagous armoured scale insect pest Hemiberlesia lataniae (Hemiptera: Diaspididae). A cDNA microarray consisting of 17,512 unigenes selected from over 132,000 expressed sequence tags (ESTs) was used to measure the transcriptomic profile of the A. chinensis 'Hort16A' canes in response to a controlled infestation of H. lataniae. After 2 days, 272 transcripts were differentially expressed. After 7 days, 5,284 (30%) transcripts were differentially expressed. The transcripts were grouped into 22 major functional categories using MapMan software. After 7 days, transcripts associated with photosynthesis (photosystem II) were significantly down-regulated, while those associated with secondary metabolism were significantly up-regulated. A total of 643 transcripts associated with response to stress were differentially expressed. This included biotic stress-related transcripts orthologous with pathogenesis related proteins, the phenylpropanoid pathway, NBS-LRR (R) genes, and receptor-like kinase-leucine rich repeat signalling proteins. While transcriptional studies are not conclusive in their own right, results were suggestive of a defence response involving both ETI and PTI, with predominance of the SA signalling pathway. Exogenous application of an SA-mimic decreased H. lataniae growth on A. chinensis 'Hort16A' plants in two laboratory experiments.

An upstream open reading frame is essential for feedback regulation of ascorbate biosynthesis in Arabidopsis.

Ascorbate (vitamin C) is an essential antioxidant and enzyme cofactor in both plants and animals. Ascorbate concentration is tightly regulated in plants, partly to respond to stress. Here, we demonstrate that ascorbate concentrations are determined via the posttranscriptional repression of GDP-l-galactose phosphorylase (GGP), a major control enzyme in the ascorbate biosynthesis pathway. This regulation requires a cis-acting upstream open reading frame (uORF) that represses the translation of the downstream GGP open reading frame under high ascorbate concentration. Disruption of this uORF stops the ascorbate feedback regulation of translation and results in increased ascorbate concentrations in leaves. The uORF is predicted to initiate at a noncanonical codon (ACG rather than AUG) and encode a 60- to 65-residue peptide. Analysis of ribosome protection data from Arabidopsis thaliana showed colocation of high levels of ribosomes with both the uORF and the main coding sequence of GGP. Together, our data indicate that the noncanonical uORF is translated and encodes a peptide that functions in the ascorbate inhibition of translation. This posttranslational regulation of ascorbate is likely an ancient mechanism of control as the uORF is conserved in GGP genes from mosses to angiosperms.

Engineering the anthocyanin regulatory complex of strawberry (Fragaria vesca).

The woodland strawberry, Fragaria vesca is a model fruit for a number of rosaceous crops. We have engineered altered concentrations of anthocyanin in F. vesca, to determine the impact on plant growth and fruit quality. Anthocyanin concentrations were significantly increased by over-expression or decreased by knock-down of the R2R3 MYB activator, MYB10. In contrast, a potential bHLH partner for MYB10 (bHLH33) did not affect the anthocyanin pathway when knocked down using RNAi constructs. Metabolic analysis of fruits revealed that, of all the polyphenolics surveyed, only cyanidin, and pelargonidin glucoside, and coumaryl hexose were significantly affected by over-expression and knock down of MYB10. Using the F. vesca genome sequence, members of the MYB, bHLH, and WD40 families were examined. Global analysis of gene expression and targeted qPCR analysis of biosynthetic genes and regulators confirmed the effects of altering MYB10 expression, as well as the knock-down of bHLH33. Other members of the MYB transcription factor family were affected by the transgenes. Transient expression of strawberry genes in Nicotiana benthamiana revealed that MYB10 can auto-regulate itself, and potential repressors of MYB10. In tobacco, MYB10's activation of biosynthetic steps is inhibited by the strawberry repressor MYB1.

The draft genome sequence of European pear (Pyrus communis L. 'Bartlett').

We present a draft assembly of the genome of European pear (Pyrus communis) 'Bartlett'. Our assembly was developed employing second generation sequencing technology (Roche 454), from single-end, 2 kb, and 7 kb insert paired-end reads using Newbler (version 2.7). It contains 142,083 scaffolds greater than 499 bases (maximum scaffold length of 1.2 Mb) and covers a total of 577.3 Mb, representing most of the expected 600 Mb Pyrus genome. A total of 829,823 putative single nucleotide polymorphisms (SNPs) were detected using re-sequencing of 'Louise Bonne de Jersey' and 'Old Home'. A total of 2,279 genetically mapped SNP markers anchor 171 Mb of the assembled genome. Ab initio gene prediction combined with prediction based on homology searching detected 43,419 putative gene models. Of these, 1219 proteins (556 clusters) are unique to European pear compared to 12 other sequenced plant genomes. Analysis of the expansin gene family provided an example of the quality of the gene prediction and an insight into the relationships among one class of cell wall related genes that control fruit softening in both European pear and apple (Malus × domestica). The 'Bartlett' genome assembly v1.0 (http://www.rosaceae.org/species/pyrus/pyrus_communis/genome_v1.0) is an invaluable tool for identifying the genetic control of key horticultural traits in pear and will enable the wide application of marker-assisted and genomic selection that will enhance the speed and efficiency of pear cultivar development.

An ancient duplication of apple MYB transcription factors is responsible for novel red fruit-flesh phenotypes.

Anthocyanin accumulation is coordinated in plants by a number of conserved transcription factors. In apple (Malus × domestica), an R2R3 MYB transcription factor has been shown to control fruit flesh and foliage anthocyanin pigmentation (MYB10) and fruit skin color (MYB1). However, the pattern of expression and allelic variation at these loci does not explain all anthocyanin-related apple phenotypes. One such example is an open-pollinated seedling of cv Sangrado that has green foliage and develops red flesh in the fruit cortex late in maturity. We used methods that combine plant breeding, molecular biology, and genomics to identify duplicated MYB transcription factors that could control this phenotype. We then demonstrated that the red-flesh cortex phenotype is associated with enhanced expression of MYB110a, a paralog of MYB10. Functional characterization of MYB110a showed that it was able to up-regulate anthocyanin biosynthesis in tobacco (Nicotiana tabacum). The chromosomal location of MYB110a is consistent with a whole-genome duplication event that occurred during the evolution of apple within the Maloideae family. Both MYB10 and MYB110a have conserved function in some cultivars, but they differ in their expression pattern and response to fruit maturity.

GENCODE: producing a reference annotation for ENCODE.

BACKGROUND: The GENCODE consortium was formed to identify and map all protein-coding genes within the ENCODE regions. This was achieved by a combination of initial manual annotation by the HAVANA team, experimental validation by the GENCODE consortium and a refinement of the annotation based on these experimental results. RESULTS: The GENCODE gene features are divided into eight different categories of which only the first two (known and novel coding sequence) are confidently predicted to be protein-coding genes. 5' rapid amplification of cDNA ends (RACE) and RT-PCR were used to experimentally verify the initial annotation. Of the 420 coding loci tested, 229 RACE products have been sequenced. They supported 5' extensions of 30 loci and new splice variants in 50 loci. In addition, 46 loci without evidence for a coding sequence were validated, consisting of 31 novel and 15 putative transcripts. We assessed the comprehensiveness of the GENCODE annotation by attempting to validate all the predicted exon boundaries outside the GENCODE annotation. Out of 1,215 tested in a subset of the ENCODE regions, 14 novel exon pairs were validated, only two of them in intergenic regions. CONCLUSION: In total, 487 loci, of which 434 are coding, have been annotated as part of the GENCODE reference set available from the UCSC browser. Comparison of GENCODE annotation with RefSeq and ENSEMBL show only 40% of GENCODE exons are contained within the two sets, which is a reflection of the high number of alternative splice forms with unique exons annotated. Over 50% of coding loci have been experimentally verified by 5' RACE for EGASP and the GENCODE collaboration is continuing to refine its annotation of 1% human genome with the aid of experimental validation.

The Ensembl analysis pipeline.

The Ensembl pipeline is an extension to the Ensembl system which allows automated annotation of genomic sequence. The software comprises two parts. First, there is a set of Perl modules ("Runnables" and "RunnableDBs") which are 'wrappers' for a variety of commonly used analysis tools. These retrieve sequence data from a relational database, run the analysis, and write the results back to the database. They inherit from a common interface, which simplifies the writing of new wrapper modules. On top of this sits a job submission system (the "RuleManager") which allows efficient and reliable submission of large numbers of jobs to a compute farm. Here we describe the fundamental software components of the pipeline, and we also highlight some features of the Sanger installation which were necessary to enable the pipeline to scale to whole-genome analysis.

An overview of Ensembl.

Ensembl (http://www.ensembl.org/) is a bioinformatics project to organize biological information around the sequences of large genomes. It is a comprehensive source of stable automatic annotation of individual genomes, and of the synteny and orthology relationships between them. It is also a framework for integration of any biological data that can be mapped onto features derived from the genomic sequence. Ensembl is available as an interactive Web site, a set of flat files, and as a complete, portable open source software system for handling genomes. All data are provided without restriction, and code is freely available. Ensembl's aims are to continue to "widen" this biological integration to include other model organisms relevant to understanding human biology as they become available; to "deepen" this integration to provide an ever more seamless linkage between equivalent components in different species; and to provide further classification of functional elements in the genome that have been previously elusive.