Elucidation of resistance signaling and identification of powdery mildew resistant mapping loci (ClaPMR2) during watermelon-Podosphaera xanthii interaction using RNA-Seq and whole-genome resequencing approach.

Watermelon is an important vegetable crop and is widely cultivated in USA with an approximate global production of > 100 million tons. Powdery mildew (PM) caused by Podosphaera xanthii is a major production-limiting factor on watermelon and other cucurbits. Numerous PM and multiple disease resistant (MDR) watermelon germplasm lines have been developed by the USDA in Charleston, SC. To gain a better understanding of the innate and activated molecular defense mechanisms involved during compatible and incompatible PM-watermelon interactions, we inoculated PM susceptible (USVL677-PMS) and resistant (USVL531-MDR) watermelon plants with 105 conidia ml-1 of P. xanthii. RNA-seq profiling was done on leaf samples collected at 0, 1, 3, and 8 days post inoculation (DPI). A total of 2,566 unique differentially expressed genes (DEGs) were identified between compatible and incompatible interactions with P. xanthii. The compatible interactions resulted in distinct plant gene activation (> twofold unique transcripts, 335:191:1762 :: 1:3:8 DPI) as compared to incompatible interaction (> twofold unique transcripts, 314:681:487 :: 1:3:8 DPI). Further, comparative whole-genome resequencing analysis of USVL531-PMR, USVL677-PMS and four introgressed PM resistant recombinant inbred lines (RIL, USVL531-PMR × USVL677-PMS) were performed to identify the region of PM resistance introgressed break points along with other traits inherent by USVL531-PMR by comparing the SNPs and InDels. Based on SNPs identification and CAPS markers, the resistance gene was identified as ClaPMR2, Citrullus lanatus PM Resistance gene 2 {Chr2 : 26750001 .. 26753327 (-)}, a NBS-LRR resistance protein (R) with homology to the Arabidopsis thaliana PM resistance protein, RPW8. The transcriptome data also revealed a complex regulatory network associated with the introgressed junctions mediated by PM resistance R proteins (R genes) that may involve multiple signal regulators and transducers, carbohydrate metabolism, cell wall modifications and the hormone-signaling pathway.

A Targeted Sequencing Assay for Serotyping Escherichia coli Using AgriSeq Technology.

The gold standard method for serotyping Escherichia coli has relied on antisera-based typing of the O- and H-antigens, which is labor intensive and often unreliable. In the post-genomic era, sequence-based assays are potentially faster to provide results, could combine O-serogrouping and H-typing in a single test, and could simultaneously screen for the presence of other genetic markers of interest such as virulence factors. Whole genome sequencing is one approach; however, this method has limited multiplexing capabilities, and only a small fraction of the sequence is informative for subtyping or identifying virulence potential. A targeted, sequence-based assay and accompanying software for data analysis would be a great improvement over the currently available methods for serotyping. The purpose of this study was to develop a high-throughput, molecular method for serotyping E. coli by sequencing the genes that are required for production of O- and H-antigens, as well as to develop software for data analysis and serotype identification. To expand the utility of the assay, targets for the virulence factors, Shiga toxins (stx 1, and stx 2) and intimin (eae) were included. To validate the assay, genomic DNA was extracted from O-serogroup and H-type standard strains and from Shiga toxin-producing E. coli, the targeted regions were amplified, and then sequencing libraries were prepared from the amplified products followed by sequencing of the libraries on the Ion S5™ sequencer. The resulting sequence files were analyzed via the SeroType Caller™ software for identification of O-serogroup, H-type, and presence of stx 1 , stx 2, and eae. We successfully identified 169 O-serogroups and 41 H-types. The assay also routinely detected the presence of stx 1a,c,d (3 of 3 strains), stx 2c-e,g (8 of 8 strains), stx 2f (1 strain), and eae (6 of 6 strains). Taken together, the high-throughput, sequence-based method presented here is a reliable alternative to antisera-based serotyping methods for E. coli.

Exome resequencing and GWAS for growth, ecophysiology, and chemical and metabolomic composition of wood of Populus trichocarpa.

BACKGROUND: Populus trichocarpa is an important forest tree species for the generation of lignocellulosic ethanol. Understanding the genomic basis of biomass production and chemical composition of wood is fundamental in supporting genetic improvement programs. Considerable variation has been observed in this species for complex traits related to growth, phenology, ecophysiology and wood chemistry. Those traits are influenced by both polygenic control and environmental effects, and their genome architecture and regulation are only partially understood. Genome wide association studies (GWAS) represent an approach to advance that aim using thousands of single nucleotide polymorphisms (SNPs). Genotyping using exome capture methodologies represent an efficient approach to identify specific functional regions of genomes underlying phenotypic variation. RESULTS: We identified 813 K SNPs, which were utilized for genotyping 461 P. trichocarpa clones, representing 101 provenances collected from Oregon and Washington, and established in California. A GWAS performed on 20 traits, considering single SNP-marker tests identified a variable number of significant SNPs (p-value < 6.1479E-8) in association with diameter, height, leaf carbon and nitrogen contents, and δ15N. The number of significant SNPs ranged from 2 to 220 per trait. Additionally, multiple-marker analyses by sliding-windows tests detected between 6 and 192 significant windows for the analyzed traits. The significant SNPs resided within genes that encode proteins belonging to different functional classes as such protein synthesis, energy/metabolism and DNA/RNA metabolism, among others. CONCLUSIONS: SNP-markers within genes associated with traits of importance for biomass production were detected. They contribute to characterize the genomic architecture of P. trichocarpa biomass required to support the development and application of marker breeding technologies.

Phenotypic and Genomic Local Adaptation across Latitude and Altitude in Populus trichocarpa.

Local adaptation to climate allows plants to cope with temporally and spatially heterogeneous environments, and parallel phenotypic clines provide a natural experiment to uncover the genomic architecture of adaptation. Though extensive effort has been made to investigate the genomic basis of local adaptation to climate across the latitudinal range of tree species, less is known for altitudinal clines. We used exome capture to genotype 451 Populus trichocarpa genotypes across altitudinal and latitudinal gradients spanning the natural species range, and phenotyped these trees for a variety of adaptive traits in two common gardens. We observed clinal variation in phenotypic traits across the two transects, which indicates climate-driven selection, and coupled gene-based genotype-phenotype and genotype-environment association scans to identify imprints of climatic adaptation on the genome. Although many of the phenotype- and climate-associated genes were unique to one transect, we found evidence of parallelism between latitude and altitude, as well as significant convergence when we compared our outlier genes with those putatively involved in climatic adaptation in two gymnosperm species. These results suggest that not only genomic constraint during adaptation to similar environmental gradients in poplar but also different environmental contexts, spatial scale, and perhaps redundant function among potentially adaptive genes and polymorphisms lead to divergent adaptive architectures.

Differential roles of melatonin in plant-host resistance and pathogen suppression in cucurbits.

Since the 1950s, research on the animal neurohormone, melatonin, has focused on its multiregulatory effect on patients suffering from insomnia, cancer, and Alzheimer's disease. In plants, melatonin plays major role in plant growth and development, and is inducible in response to diverse biotic and abiotic stresses. However, studies on the direct role of melatonin in disease suppression and as a signaling molecule in host-pathogen defense mechanism are lacking. This study provides insight on the predicted biosynthetic pathway of melatonin in watermelon (Citrullus lanatus), and how application of melatonin, an environmental-friendly immune inducer, can boost plant immunity and suppress pathogen growth where fungicide resistance and lack of genetic resistance are major problems. We evaluated the effect of spray-applied melatonin and also transformed watermelon plants with the melatonin biosynthetic gene SNAT (serotonin N-acetyltransferase) to determine the role of melatonin in plant defense. Increased melatonin levels in plants were found to boost resistance against the foliar pathogen Podosphaera xanthii (powdery mildew), and the soil-borne oomycete Phytophthora capsici in watermelon and other cucurbits. Further, transcriptomic data on melatonin-sprayed (1 mmol/L) watermelon leaves suggest that melatonin alters the expression of genes involved in both PAMP-mediated (pathogen-associated molecular pattern) and ETI-mediated (effector-triggered immunity) defenses. Twenty-seven upregulated genes were associated with constitutive defense as well as initial priming of the melatonin-induced plant resistance response. Our results indicate that developing strategies to increase melatonin levels in specialty crops such as watermelon can lead to resistance against diverse filamentous pathogens.

Convergent local adaptation to climate in distantly related conifers.

When confronted with an adaptive challenge, such as extreme temperature, closely related species frequently evolve similar phenotypes using the same genes. Although such repeated evolution is thought to be less likely in highly polygenic traits and distantly related species, this has not been tested at the genome scale. We performed a population genomic study of convergent local adaptation among two distantly related species, lodgepole pine and interior spruce. We identified a suite of 47 genes, enriched for duplicated genes, with variants associated with spatial variation in temperature or cold hardiness in both species, providing evidence of convergent local adaptation despite 140 million years of separate evolution. These results show that adaptation to climate can be genetically constrained, with certain key genes playing nonredundant roles.

Conservation and divergence of gene expression plasticity following c. 140 million years of evolution in lodgepole pine (Pinus contorta) and interior spruce (Picea glauca×Picea engelmannii).

Species respond to environmental stress through a combination of genetic adaptation and phenotypic plasticity, both of which may be important for survival in the face of climatic change. By characterizing the molecular basis of plastic responses and comparing patterns among species, it is possible to identify how such traits evolve. Here, we used de novo transcriptome assembly and RNAseq to explore how patterns of gene expression differ in response to temperature, moisture, and light regime treatments in lodgepole pine (Pinus contorta) and interior spruce (a natural hybrid population of Picea glauca and Picea engelmannii). We found wide evidence for an effect of treatment on expression within each species, with 6413 and 11,658 differentially expressed genes identified in spruce and pine, respectively. Comparing patterns of expression among these species, we found that 74% of all orthologs with differential expression had a pattern that was conserved in both species, despite 140 million yr of evolution. We also found that the specific treatments driving expression patterns differed between genes with conserved versus diverged patterns of expression. We conclude that natural selection has probably played a role in shaping plastic responses to environment in these species.

Evidence for extensive heterotrophic metabolism, antioxidant action, and associated regulatory events during winter hardening in Sitka spruce.

BACKGROUND: Cold acclimation in woody perennials is a metabolically intensive process, but coincides with environmental conditions that are not conducive to the generation of energy through photosynthesis. While the negative effects of low temperatures on the photosynthetic apparatus during winter have been well studied, less is known about how this is reflected at the level of gene and metabolite expression, nor how the plant generates primary metabolites needed for adaptive processes during autumn. RESULTS: The MapMan tool revealed enrichment of the expression of genes related to mitochondrial function, antioxidant and associated regulatory activity, while changes in metabolite levels over the time course were consistent with the gene expression patterns observed. Genes related to thylakoid function were down-regulated as expected, with the exception of plastid targeted specific antioxidant gene products such as thylakoid-bound ascorbate peroxidase, components of the reactive oxygen species scavenging cycle, and the plastid terminal oxidase. In contrast, the conventional and alternative mitochondrial electron transport chains, the tricarboxylic acid cycle, and redox-associated proteins providing reactive oxygen species scavenging generated by electron transport chains functioning at low temperatures were all active. CONCLUSIONS: A regulatory mechanism linking thylakoid-bound ascorbate peroxidase action with "chloroplast dormancy" is proposed. Most importantly, the energy and substrates required for the substantial metabolic remodeling that is a hallmark of freezing acclimation could be provided by heterotrophic metabolism.

Mining and visualization of microarray and metabolomic data reveal extensive cell wall remodeling during winter hardening in Sitka spruce (Picea sitchensis).

Microarray gene expression profiling is a powerful technique to understand complex developmental processes, but making biologically meaningful inferences from such studies has always been challenging. We previously reported a microarray study of the freezing acclimation period in Sitka spruce (Picea sitchensis) in which a large number of candidate genes for climatic adaptation were identified. In the current paper, we apply additional systems biology tools to these data to further probe changes in the levels of genes and metabolites and activities of associated pathways that regulate this complex developmental transition. One aspect of this adaptive process that is not well understood is the role of the cell wall. Our data suggest coordinated metabolic and signaling responses leading to cell wall remodeling. Co-expression of genes encoding proteins associated with biosynthesis of structural and non-structural cell wall carbohydrates was observed, which may be regulated by ethylene signaling components. At the same time, numerous genes, whose products are putatively localized to the endomembrane system and involved in both the synthesis and trafficking of cell wall carbohydrates, were up-regulated. Taken together, these results suggest a link between ethylene signaling and biosynthesis, and targeting of cell wall related gene products during the period of winter hardening. Automated Layout Pipeline for Inferred NEtworks (ALPINE), an in-house plugin for the Cytoscape visualization environment that utilizes the existing GeneMANIA and Mosaic plugins, together with the use of visualization tools, provided images of proposed signaling processes that became active over the time course of winter hardening, particularly at later time points in the process. The resulting visualizations have the potential to reveal novel, hypothesis-generating, gene association patterns in the context of targeted subcellular location.

Divergent selection and heterogeneous migration rates across the range of Sitka spruce (Picea sitchensis).

Gene flow and effective population size (N(e)) should depend on a population's position within its range: those near the edges are expected to have smaller N(e) and lower relative emigration rates, whereas those nearer the centre should have larger N(e) and higher relative emigration rates. In species with continuous ranges, this phenomenon may limit the ability of peripheral populations to respond to divergent selection. Here, we employ Sitka spruce as a model to test these predictions. We previously genotyped 339 single nucleotide polymorphisms (SNPs) in 410 individuals from 13 populations, and used these data to identify putative targets of divergent selection, as well as to explore the extent to which central-peripheral structure may impede adaptation. Fourteen SNPs had outlier F(ST) estimates suggestive of divergent selection, of which nine were previously associated with phenotypic variation in adaptive traits (timing of autumn budset and cold hardiness). Using coalescent simulations, we show that populations from near the centre of the range have higher effective populations sizes than those from the edges, and that central populations contribute more migrants to marginal populations than the reverse. Our results suggest that while divergent selection appears to have shaped allele frequencies among populations, asymmetrical movement of alleles from the centre to the edges of the species range may affect the adaptive capacity of peripheral populations. In southern peripheral populations, the movement of cold-adapted alleles from the north represents a significant impediment to adaptation under climate change, while in the north, movement of warm-adapted alleles from the south may enhance adaptation.