Genome-wide association study of 23 flowering phenology traits and 4 floral agronomic traits in tree peony (Paeonia section Moutan DC.) reveals five genes known to regulate flowering time.

Tree peony is a unique traditional flower in China, with large, fragrant, and colorful flowers. However, a relatively short and concentrated flowering period limits the applications and production of tree peony. A genome-wide association study (GWAS) was conducted to accelerate molecular breeding for the improvement of flowering phenology traits and ornamental phenotypes in tree peony. A diverse panel of 451 tree peony accessions was phenotyped for 23 flowering phenology traits and 4 floral agronomic traits over 3 years. Genotyping by sequencing (GBS) was used to obtain a large number of genome-wide single-nucleotide polymorphisms (SNPs) (107 050) for the panel genotypes, and 1047 candidate genes were identified by association mapping. Eighty-two related genes were observed during at least 2 years for flowering, and seven SNPs repeatedly identified for multiple flowering phenology traits over multiple years were highly significantly associated with five genes known to regulate flowering time. We validated the temporal expression profiles of these candidate genes and highlighted their possible roles in the regulation of flower bud differentiation and flowering time in tree peony. This study shows that GWAS based on GBS can be used to identify the genetic determinants of complex traits in tree peony. The results expand our understanding of flowering time control in perennial woody plants. Identification of markers closely related to these flowering phenology traits can be used in tree peony breeding programs for important agronomic traits.

Host genomic influence on bacterial composition in the switchgrass rhizosphere.

Host genetic variation can shape the diversity and composition of associated microbiomes, which may reciprocally influence host traits and performance. While the genetic basis of phenotypic diversity of plant populations in nature has been studied, comparatively little research has investigated the genetics of host effects on their associated microbiomes. Switchgrass (Panicum virgatum) is a highly outcrossing, perennial, grass species with substantial locally adaptive diversity across its native North American range. Here, we compared 383 switchgrass accessions in a common garden to determine the host genotypic influence on rhizosphere bacterial composition. We hypothesized that the composition and diversity of rhizosphere bacterial assemblages would differentiate due to genotypic differences between hosts (potentially due to root phenotypes and associated life history variation). We observed higher alpha diversity of bacteria associated with upland ecotypes and tetraploids, compared to lowland ecotypes and octoploids, respectively. Alpha diversity correlated negatively with flowering time and plant height, indicating that bacterial composition varies along switchgrass life history axes. Narrow-sense heritability (h2 ) of the relative abundance of 21 core bacterial families was observed. Overall compositional differences among tetraploids, due to genetic variation, supports widespread genotypic influence on the rhizosphere microbiome. Tetraploids were only considered due to complexities associated with the octoploid genomes. Lastly, a genome-wide association study identified 1861 single-nucleotide polymorphisms associated with 110 families and genes containing them related to potential regulatory functions. Our findings suggest that switchgrass genomic and life-history variation influences bacterial composition in the rhizosphere, potentially due to host adaptation to local environments.

Integrated analysis of miRNAome transcriptome and degradome reveals miRNA-target modules governing floral florescence development and senescence across early- and late-flowering genotypes in tree peony.

As a candidate national flower of China, tree peony has extremely high ornamental, medicinal and oil value. However, the short florescence and rarity of early-flowering and late-flowering varieties restrict further improvement of the economic value of tree peony. Specific miRNAs and their target genes engaged in tree peony floral florescence, development and senescence remain unknown. This report presents the integrated analysis of the miRNAome, transcriptome and degradome of tree peony petals collected from blooming, initial flowering, full blooming and decay stages in early-flowering variety Paeonia ostii 'Fengdan', an early-flowering mutant line of Paeonia ostii 'Fengdan' and late-flowering variety Paeonia suffruticosa 'Lianhe'. Transcriptome analysis revealed a transcript ('psu.G.00014095') which was annotated as a xyloglucan endotransglycosylase/hydrolase precursor XTH-25 and found to be differentially expressed across flower developmental stages in Paeonia ostii 'Fengdan' and Paeonia suffruticosa 'Lianhe'. The miRNA-mRNA modules were presented significant enrichment in various pathways such as plant hormone signal transduction, indole alkaloid biosynthesis, arachidonic acid metabolism, folate biosynthesis, fatty acid elongation, and the MAPK signaling pathway. Multiple miRNA-mRNA-TF modules demonstrated the potential functions of MYB-related, bHLH, Trihelix, NAC, GRAS and HD-ZIP TF families in floral florescence, development, and senescence of tree peony. Comparative spatio-temporal expression investigation of eight floral-favored miRNA-target modules suggested that transcript 'psu.T.00024044' and microRNA mtr-miR166g-5p are involved in the floral florescence, development and senescence associated agronomic traits of tree peony. The results might accelerate the understanding of the potential regulation mechanism in regards to floral florescence, development and abscission, and supply guidance for tree peony breeding of varieties with later and longer florescence characteristics.

A high-quality reference genome for Fraxinus pennsylvanica for ash species restoration and research.

Green ash (Fraxinus pennsylvanica) is the most widely distributed ash tree in North America. Once common, it has experienced high mortality from the non-native invasive emerald ash borer (EAB; Agrilus planipennis). A small percentage of native green ash trees that remain healthy in long-infested areas, termed "lingering ash," display partial resistance to the insect, indicating that breeding and propagating populations with higher resistance to EAB may be possible. To assist in ash breeding, ecology and evolution studies, we report the first chromosome-level assembly from the genus Fraxinus for F. pennsylvanica with over 99% of bases anchored to 23 haploid chromosomes, spanning 757 Mb in total, composed of 49.43% repetitive DNA, and containing 35,470 high-confidence gene models assigned to 22,976 Asterid orthogroups. We also present results of range-wide genetic variation studies, the identification of candidate genes for important traits including potential EAB-resistance genes, and an investigation of comparative genome organization among Asterids based on this reference genome platform. Residual duplicated regions within the genome probably resulting from a recent whole genome duplication event in Oleaceae were visualized in relation to wild olive (Olea europaea var. sylvestris). We used our F. pennsylvanica chromosome assembly to construct reference-guided assemblies of 27 previously sequenced Fraxinus taxa, including F. excelsior. Thus, we present a significant step forward in genomic resources for research and protection of Fraxinus species.

Transcriptome analysis of contrasting resistance to herbivory by Empoasca fabae in two shrub willow species and their hybrid progeny.

Short rotation woody biomass cultivars developed from fast-growing shrub species of willow (Salix spp.) have superior properties as perennial energy crops for the Northeast and Midwest US. However, the insect pest potato leafhopper (PLH) Empoasca fabae (Harris) can cause serious damage and reduce yield of susceptible genotypes. Currently, the willow cultivars in use display varying levels of susceptibility under PLH infestation. However, genes and markers for resistance to PLH are not yet available for marker-assisted selection in breeding. In this study, transcriptome differences between a resistant genotype 94006 (S. purpurea) and a susceptible cultivar 'Jorr' (S. viminalis), and their hybrid progeny were determined. Over 600 million RNA-Seq reads were generated and mapped to the Salix purpurea reference transcriptome. Gene expression analyses revealed the unique defense mechanism in resistant genotype 94006 that involves PLH-induced secondary cell wall modification. In the susceptible genotypes, genes involved in programed cell death were highly expressed, explaining the necrosis symptoms after PLH feeding. Overall, the discovery of resistance genes and defense mechanisms provides new resources for shrub willow breeding and research in the future.

Factoring Ecological, Societal, and Economic Considerations into Inoculant Development.

There are many paths toward effective microbial inoculants for agriculture. Considering what is practical for the present day technological and farming landscape should not limit our creativity in developing innovative technologies. However, factors including production costs, practicality of implementation, and technology adoption by farmers will drive the success of new management approaches.

NaCl-induced CsRCI2E and CsRCI2F interact with aquaporin CsPIP2;1 to reduce water transport in Camelina sativa L.

Rare cold-inducible 2 (RCI2) proteins are small hydrophobic proteins that are known to be localized in cellular membranes. The function of RCI2 proteins has been reported to be associated with low-temperature, salt, and drought stress tolerances as a membrane potential regulator; however, the specific functions are still unknown. The PIP2 (plasma membrane intrinsic protein 2) aquaporins are proteins that transport water and small solutes into the cell. The expression and activity of PIP2 proteins, like RCI2, are also related to salt- and drought-stress tolerance. In this study, we identified novel protein interactions between RCI2 and PIP2; 1, including protein accumulation changes in the bioenergy crop Camelina sativa L. under various NaCl stress conditions. Accumulation of both CsRCI2E and CsRCI2F proteins increased with NaCl stress; however, to differing levels depending on the NaCl stress intensity. A co-immunoprecipitation test revealed interaction between CsRCI2E-CsPIP2 and CsRCI2F-CsPIP2. Moreover, co-expression of the four CsRCI2 proteins with CsPIP2; 1 in Xenopus laevis oocytes reduced water transport activity. Furthermore, the abundance of CsPIP2; 1 protein was decreased under CsRCI2E and CsRCI2F co-expression. These results suggest that NaCl-induced expression of CsRCI2E and CsRCI2F contributes to the regulation of CsPIP2; 1.

Host-plant induced changes in microbial community structure and midgut gene expression in an invasive polyphage (Anoplophora glabripennis).

Polyphagous insect herbivores possess diverse mechanisms to overcome challenges of feeding in multiple plant species including, but not limited to, transcriptional plasticity and associations with obligate or facultative symbionts. The Asian longhorned beetle (Anoplophora glabripennis) is a polyphagous wood-feeder capable of developing on over 100 tree species and, like other polyphages, its genome contains amplifications of digestive and detoxification genes. This insect also possesses a diverse gut microbial community, which has the metabolic potential to augment digestive physiology. While the genomic repertoires of A. glabripennis and its microbial community have been studied previously, comparatively less is known about how the gut transcriptome and community change in response to feeding in different hosts. In this study, we show that feeding in two suitable hosts (Acer spp. and Populus nigra) altered the expression levels of multicopy genes linked to digestion and detoxification. However, feeding in a host with documented resistance (Populus tomentosa) induced changes in the transcriptome and community beyond what was observed in insects reared in P. nigra, including the downregulation of numerous ß-glucosidases, odorant binding proteins, and juvenile hormone binding proteins, the upregulation of several cuticular genes, and the loss of one major bacterial family from the gut community.

Hardwood Tree Genomics: Unlocking Woody Plant Biology.

Woody perennial angiosperms (i.e., hardwood trees) are polyphyletic in origin and occur in most angiosperm orders. Despite their independent origins, hardwoods have shared physiological, anatomical, and life history traits distinct from their herbaceous relatives. New high-throughput DNA sequencing platforms have provided access to numerous woody plant genomes beyond the early reference genomes of Populus and Eucalyptus, references that now include willow and oak, with pecan and chestnut soon to follow. Genomic studies within these diverse and undomesticated species have successfully linked genes to ecological, physiological, and developmental traits directly. Moreover, comparative genomic approaches are providing insights into speciation events while large-scale DNA resequencing of native collections is identifying population-level genetic diversity responsible for variation in key woody plant biology across and within species. Current research is focused on developing genomic prediction models for breeding, defining speciation and local adaptation, detecting and characterizing somatic mutations, revealing the mechanisms of gender determination and flowering, and application of systems biology approaches to model complex regulatory networks underlying quantitative traits. Emerging technologies such as single-molecule, long-read sequencing is being employed as additional woody plant species, and genotypes within species, are sequenced, thus enabling a comparative ("evo-devo") approach to understanding the unique biology of large woody plants. Resource availability, current genomic and genetic applications, new discoveries and predicted future developments are illustrated and discussed for poplar, eucalyptus, willow, oak, chestnut, and pecan.

Development of novel genic microsatellite markers from transcriptome sequencing in sugar maple (Acer saccharum Marsh.).

BACKGROUND: Sugar maple (Acer saccharum Marsh.) is a hardwood tree species native to northeastern North America and economically valued for its wood and sap. Yet, few molecular genetic resources have been developed for this species to date. Microsatellite markers have been a useful tool in population genetics, e.g., to monitor genetic variation and to analyze gene flow patterns. The objective of this study is to develop a reference transcriptome and microsatellite markers in sugar maple. FINDINGS: A set of 117,861 putative unique transcripts were assembled using 29.2 Gb of RNA sequencing data derived from different tissues and stress treatments. From this set of sequences a total of 1068 microsatellite motifs were identified. Out of 58 genic microsatellite markers tested on a population of 47 sugar maple trees in upper Michigan, 22 amplified well, of which 16 were polymorphic and 6 were monomorphic. Values for expected heterozygosity varied from 0.224 to 0.726 for individual loci. Of the 16 polymorphic markers, 15 exhibited transferability to other Acer L. species. CONCLUSIONS: Genic microsatellite markers can be applied to analyze genetic variation in potentially adaptive genes relative to genomic reference markers as a basis for the management of sugar maple genetic resources in the face of climate change.

Genome-wide association study identifies a major gene for beech bark disease resistance in American beech (Fagus grandifolia Ehrh.).

BACKGROUND: The American Beech tree (Fagus grandifolia Ehrh.), native to eastern North America, is ecologically important and provides high quality wood products. This species is susceptible to beech bark disease (BBD) and is facing high rates of mortality in North America. The disease occurs from an interaction between the woolly beech scale insect (Cryptococcus fagisuga), one of two species of the fungus Neonectria (N. faginata or N. ditissima), and American Beech trees. METHODS: In this case-control genome-wide association study (GWAS), we tested 16 K high quality SNPs using the Affymetrix Axiom 1.5 K - 50 K assay to genotype an association population of 514 individuals. We also conducted linkage analysis in a full-sib family of 115 individuals. Fisher's exact test and logistic regression tests were performed to test associations between SNPs and phenotypes. RESULTS: Association tests revealed four highly significant SNPs on chromosome (Chr) 5 for a single gene (Mt), which encodes a mRNA for metallothionein-like protein (metal ion binding) in Fagus sylvatica. Metallothioneins represent Cys-rich metal chelators able to coordinate metal atoms and may play an important role in the resistance mechanisms against beech scale insect. CONCLUSION: The GWAS study has identified a single locus of major effect contributing to beech bark disease resistance. Knowledge of this genetic locus contributing to resistance might be used in applied breeding, conservation and restoration programs.

High-quality genetic mapping with ddRADseq in the non-model tree Quercus rubra.

BACKGROUND: Restriction site associated DNA sequencing (RADseq) has the potential to be a broadly applicable, low-cost approach for high-quality genetic linkage mapping in forest trees lacking a reference genome. The statistical inference of linear order must be as accurate as possible for the correct ordering of sequence scaffolds and contigs to chromosomal locations. Accurate maps also facilitate the discovery of chromosome segments containing allelic variants conferring resistance to the biotic and abiotic stresses that threaten forest trees worldwide. We used ddRADseq for genetic mapping in the tree Quercus rubra, with an approach optimized to produce a high-quality map. Our study design also enabled us to model the results we would have obtained with less depth of coverage. RESULTS: Our sequencing design produced a high sequencing depth in the parents (248×) and a moderate sequencing depth (15×) in the progeny. The digital normalization method of generating a de novo reference and the SAMtools SNP variant caller yielded the most SNP calls (78,725). The major drivers of map inflation were multiple SNPs located within the same sequence (77% of SNPs called). The highest quality map was generated with a low level of missing data (5%) and a genome-wide threshold of 0.025 for deviation from Mendelian expectation. The final map included 849 SNP markers (1.8% of the 78,725 SNPs called). Downsampling the individual FASTQ files to model lower depth of coverage revealed that sequencing the progeny using 96 samples per lane would have yielded too few SNP markers to generate a map, even if we had sequenced the parents at depth 248×. CONCLUSIONS: The ddRADseq technology produced enough high-quality SNP markers to make a moderately dense, high-quality map. The success of this project was due to high depth of coverage of the parents, moderate depth of coverage of the progeny, a good framework map, an optimized bioinformatics pipeline, and rigorous premapping filters. The ddRADseq approach is useful for the construction of high-quality genetic maps in organisms lacking a reference genome if the parents and progeny are sequenced at sufficient depth. Technical improvements in reduced representation sequencing (RRS) approaches are needed to reduce the amount of missing data.

The green ash transcriptome and identification of genes responding to abiotic and biotic stresses.

BACKGROUND: To develop a set of transcriptome sequences to support research on environmental stress responses in green ash (Fraxinus pennsylvanica), we undertook deep RNA sequencing of green ash tissues under various stress treatments. The treatments, including emerald ash borer (EAB) feeding, heat, drought, cold and ozone, were selected to mimic the increasing threats of climate change and invasive pests faced by green ash across its native habitat. RESULTS: We report the generation and assembly of RNA sequences from 55 green ash samples into 107,611 putative unique transcripts (PUTs). 52,899 open reading frames were identified. Functional annotation of the PUTs by comparison to the Uniprot protein database identified matches for 63 % of transcripts and for 98 % of transcripts with ORFs. Further functional annotation identified conserved protein domains and assigned gene ontology terms to the PUTs. Examination of transcript expression across different RNA libraries revealed that expression patterns clustered based on tissues regardless of stress treatment. The transcripts from stress treatments were further examined to identify differential expression. Tens to hundreds of differentially expressed PUTs were identified for each stress treatment. A set of 109 PUTs were found to be consistently up or down regulated across three or more different stress treatments, representing basal stress response candidate genes in green ash. In addition, 1956 simple sequence repeats were identified in the PUTs, of which we identified 465 high quality DNA markers and designed flanking PCR primers. CONCLUSIONS: North American native ash trees have suffered extensive mortality due to EAB infestation, creating a need to breed or select for resistant green ash genotypes. Stress from climate change is an additional concern for longevity of native ash populations. The use of genomics could accelerate management efforts. The green ash transcriptome we have developed provides important sequence information, genetic markers and stress-response candidate genes.

Genome Sequence of Fusarium Isolate MYA-4552 from the Midgut of Anoplophora glabripennis, an Invasive, Wood-Boring Beetle.

The Fusarium solani species complex (FSSC) is a clade of environmentally ubiquitous fungi that includes plant, animal, and insect associates. Here, we report the draft genome sequence of the undescribed species FSSC 6 (isolate MYA-4552), housed in the gut of the wood-boring cerambycid beetle Anoplophora glabripennis.

Preliminary Genomic Characterization of Ten Hardwood Tree Species from Multiplexed Low Coverage Whole Genome Sequencing.

Forest health issues are on the rise in the United States, resulting from introduction of alien pests and diseases, coupled with abiotic stresses related to climate change. Increasingly, forest scientists are finding genetic/genomic resources valuable in addressing forest health issues. For a set of ten ecologically and economically important native hardwood tree species representing a broad phylogenetic spectrum, we used low coverage whole genome sequencing from multiplex Illumina paired ends to economically profile their genomic content. For six species, the genome content was further analyzed by flow cytometry in order to determine the nuclear genome size. Sequencing yielded a depth of 0.8X to 7.5X, from which in silico analysis yielded preliminary estimates of gene and repetitive sequence content in the genome for each species. Thousands of genomic SSRs were identified, with a clear predisposition toward dinucleotide repeats and AT-rich repeat motifs. Flanking primers were designed for SSR loci for all ten species, ranging from 891 loci in sugar maple to 18,167 in redbay. In summary, we have demonstrated that useful preliminary genome information including repeat content, gene content and useful SSR markers can be obtained at low cost and time input from a single lane of Illumina multiplex sequence.

Substantial genome synteny preservation among woody angiosperm species: comparative genomics of Chinese chestnut (Castanea mollissima) and plant reference genomes.

BACKGROUND: Chinese chestnut (Castanea mollissima) has emerged as a model species for the Fagaceae family with extensive genomic resources including a physical map, a dense genetic map and quantitative trait loci (QTLs) for chestnut blight resistance. These resources enable comparative genomics analyses relative to model plants. We assessed the degree of conservation between the chestnut genome and other well annotated and assembled plant genomic sequences, focusing on the QTL regions of most interest to the chestnut breeding community. RESULTS: The integrated physical and genetic map of Chinese chestnut has been improved to now include 858 shared sequence-based markers. The utility of the integrated map has also been improved through the addition of 42,970 BAC (bacterial artificial chromosome) end sequences spanning over 26 million bases of the estimated 800 Mb chestnut genome. Synteny between chestnut and ten model plant species was conducted on a macro-syntenic scale using sequences from both individual probes and BAC end sequences across the chestnut physical map. Blocks of synteny with chestnut were found in all ten reference species, with the percent of the chestnut physical map that could be aligned ranging from 10 to 39 %. The integrated genetic and physical map was utilized to identify BACs that spanned the three previously identified QTL regions conferring blight resistance. The clones were pooled and sequenced, yielding 396 sequence scaffolds covering 13.9 Mbp. Comparative genomic analysis on a microsytenic scale, using the QTL-associated genomic sequence, identified synteny from chestnut to other plant genomes ranging from 5.4 to 12.9 % of the genome sequences aligning. CONCLUSIONS: On both the macro- and micro-synteny levels, the peach, grape and poplar genomes were found to be the most structurally conserved with chestnut. Interestingly, these results did not strictly follow the expectation that decreased phylogenetic distance would correspond to increased levels of genome preservation, but rather suggest the additional influence of life-history traits on preservation of synteny. The regions of synteny that were detected provide an important tool for defining and cataloging genes in the QTL regions for advancing chestnut blight resistance research.

Characterization of the Transcriptome of the Xerophyte Ammopiptanthus mongolicus Leaves under Drought Stress by 454 Pyrosequencing.

BACKGROUND: Ammopiptanthus mongolicus (Maxim. Ex Kom.) Cheng f., an endangered ancient legume species, endemic to the Gobi desert in north-western China. As the only evergreen broadleaf shrub in this area, A. mongolicus plays an important role in the region's ecological-environmental stability. Despite the strong potential of A. mongolicus in providing new insights on drought tolerance, sequence information on the species in public databases remains scarce. To both learn about the role of gene expression in drought stress tolerance in A. mongolicus and to expand genomic resources for the species, transcriptome sequencing of stress-treated A. mongolicus plants was performed. RESULTS: Using 454 pyrosequencing technology, 8,480 and 7,474 contigs were generated after de novo assembly of RNA sequences from leaves of untreated and drought-treated plants, respectively. After clustering using TGICL and CAP3 programs, a combined assembly of all reads produced a total of 11,357 putative unique transcripts (PUTs). Functional annotation and classification of the transcripts were conducted by aligning the 11,357 PUTs against the public protein databases and nucleotide database (Nt). Between control and drought-treated plants, 1,620 differentially expressed genes (DEGs) were identified, of which 1,106 were up-regulated and 514 were down-regulated. The differential expression of twenty candidate genes in metabolic pathways and transcription factors families related to stress-response were confirmed by quantitative real-time PCR. Representatives of several large gene families, such as WRKY and P5CS, were identified and verified in A. mongolicus for the first time. CONCLUSIONS: The additional transcriptome resources, gene expression profiles, functional annotations, and candidate genes provide a more comprehensive understanding of the stress response pathways in xeric-adapted plant species such as A. mongolicus.

Functional genomics and microbiome profiling of the Asian longhorned beetle (Anoplophora glabripennis) reveal insights into the digestive physiology and nutritional ecology of wood feeding beetles.

BACKGROUND: Wood-feeding beetles harbor an ecologically rich and taxonomically diverse assemblage of gut microbes that appear to promote survival in woody tissue, which is devoid of nitrogen and essential nutrients. Nevertheless, the contributions of these apparent symbionts to digestive physiology and nutritional ecology remain uncharacterized in most beetle lineages. RESULTS: Through parallel transcriptome profiling of beetle- and microbial- derived mRNAs, we demonstrate that the midgut microbiome of the Asian longhorned beetle (Anoplophora glabripennis), a member of the beetle family Cerambycidae, is enriched in biosynthetic pathways for the synthesis of essential amino acids, vitamins, and sterols. Consequently, the midgut microbiome of A. glabripennis can provide essential nutrients that the beetle cannot obtain from its woody diet or synthesize itself. The beetle gut microbiota also produce their own suite of transcripts that can enhance lignin degradation, degrade hemicellulose, and ferment xylose and wood sugars. An abundance of cellulases from several glycoside hydrolase families are expressed endogenously by A. glabripennis, as well as transcripts that allow the beetle to convert microbe-synthesized essential amino acids into non-essential amino acids. A. glabripennis and its gut microbes likely collaborate to digest carbohydrates and convert released sugars and amino acid intermediates into essential nutrients otherwise lacking from their woody host plants. CONCLUSIONS: The nutritional provisioning capabilities of the A. glabripennis gut microbiome may contribute to the beetles' unusually broad host range. The presence of some of the same microbes in the guts of other Cerambycidae and other wood-feeding beetles suggests that partnerships with microbes may be a facilitator of evolutionary radiations in beetles, as in certain other groups of insects, allowing access to novel food sources through enhanced nutritional provisioning.

The Camelina aquaporin CsPIP2;1 is regulated by phosphorylation at Ser273, but not at Ser277, of the C-terminus and is involved in salt- and drought-stress responses.

Aquaporin (AQP) proteins are involved in water homeostasis in cells at all taxonomic levels of life. Phosphorylation of some AQPs has been proposed to regulate water permeability via gating of the channel itself. We analyzed plasma membrane intrinsic proteins (PIP) from Camelina and characterized their biological functions under both stressful and favorable conditions. A three-dimensional theoretical model of the Camelina AQP proteins was built by homology modeling which could prove useful in further functional characterization of AQPs. CsPIP2;1 was strongly and constitutively expressed in roots and leaves of Camelina, suggesting that this gene is related to maintenance of homeostasis during salt and drought stresses. CsPIP2s exhibited water channel activity in Xenopus oocytes. We then examined the roles of CsPIP2;1 phosphorylation at Ser273 and Ser277 in the regulation of water permeability using phosphorylation mutants. A single deletion strain of CsPIP2;1 was generated to serve as the primary host for testing AQP expression constructs. A Ser277 to alanine mutation (to prevent phosphorylation) did not change CsPIP2;1 water permeability while a Ser273 mutation to alanine did affect water permeability. Furthermore, a CsPIP2;1 point mutation when ectopically expressed in yeast resulted in lower growth in salt and drought conditions compared with controls, and confirmation of Ser273 as the phosphorylation site. Our results support the idea that post-translational modifications in the Ser273 regulatory domains of the C-terminus fine tune water flux through CsPIP2;1.

Genome and secretome analysis of the hemibiotrophic fungal pathogen, Moniliophthora roreri, which causes frosty pod rot disease of cacao: mechanisms of the biotrophic and necrotrophic phases.

BACKGROUND: The basidiomycete Moniliophthora roreri is the causal agent of Frosty pod rot (FPR) disease of cacao (Theobroma cacao), the source of chocolate, and FPR is one of the most destructive diseases of this important perennial crop in the Americas. This hemibiotroph infects only cacao pods and has an extended biotrophic phase lasting up to sixty days, culminating in plant necrosis and sporulation of the fungus without the formation of a basidiocarp. RESULTS: We sequenced and assembled 52.3 Mb into 3,298 contigs that represent the M. roreri genome. Of the 17,920 predicted open reading frames (OFRs), 13,760 were validated by RNA-Seq. Using read count data from RNA sequencing of cacao pods at 30 and 60 days post infection, differential gene expression was estimated for the biotrophic and necrotrophic phases of this plant-pathogen interaction. The sequencing data were used to develop a genome based secretome for the infected pods. Of the 1,535 genes encoding putative secreted proteins, 1,355 were expressed in the biotrophic and necrotrophic phases. Analysis of the data revealed secretome gene expression that correlated with infection and intercellular growth in the biotrophic phase and invasive growth and plant cellular death in the necrotrophic phase. CONCLUSIONS: Genome sequencing and RNA-Seq was used to determine and validate the Moniliophthora roreri genome and secretome. High sequence identity between Moniliophthora roreri genes and Moniliophthora perniciosa genes supports the taxonomic relationship with Moniliophthora perniciosa and the relatedness of this fungus to other basidiomycetes. Analysis of RNA-Seq data from infected plant tissues revealed differentially expressed genes in the biotrophic and necrotrophic phases. The secreted protein genes that were upregulated in the biotrophic phase are primarily associated with breakdown of the intercellular matrix and modification of the fungal mycelia, possibly to mask the fungus from plant defenses. Based on the transcriptome data, the upregulated secreted proteins in the necrotrophic phase are hypothesized to be actively attacking the plant cell walls and plant cellular components resulting in necrosis. These genes are being used to develop a new understanding of how this disease interaction progresses and to identify potential targets to reduce the impact of this devastating disease.