Common bean (Phaseolus vulgaris) is the most important food legume in the world. Although this crop is very important to both the developed and developing world as a means of dietary protein supply, resources available in common bean are limited. Global transcriptome analysis is important to better understand gene expression, genetic variation, and gene structure annotation in addition to other important features. However, the number and description of common bean sequences are very limited, which greatly inhibits genome and transcriptome research. Here we used 454 pyrosequencing to obtain a substantial transcriptome dataset for common bean.

Histone modifications such as methylation and acetylation play a significant role in controlling gene expression in unstressed and stressed plants. Genome-wide analysis of such stress-responsive modifications and genes in non-model crops is limited. We report the genome-wide profiling of histone methylation (H3K9_me2) and acetylation (H4K12_ac) in common bean (Phaseolus vulgaris L.) under rust (Uromyces appendiculatus) stress using two high-throughput approaches, chromatin immunoprecipitation sequencing (ChIP-Seq) and RNA sequencing (RNA-Seq). ChIP-Seq analysis revealed 1,235 and 556 histone methylation and acetylation responsive genes from common bean leaves treated with the rust pathogen at 0, 12 and 84 hour-after-inoculation (hai), while RNA-Seq analysis identified 145 and 1,763 genes differentially expressed between mock-inoculated and inoculated plants. The combined ChIP-Seq and RNA-Seq analyses identified some key defense responsive genes (calmodulin, cytochrome p450, chitinase, DNA Pol II, and LRR) and transcription factors (WRKY, bZIP, MYB, HSFB3, GRAS, NAC, and NMRA) in bean-rust interaction. Differential methylation and acetylation affected a large proportion of stress-responsive genes including resistant (R) proteins, detoxifying enzymes, and genes involved in ion flux and cell death. The genes identified were functionally classified using Gene Ontology (GO) and EuKaryotic Orthologous Groups (KOGs). The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis identified a putative pathway with ten key genes involved in plant-pathogen interactions. This first report of an integrated analysis of histone modifications and gene expression involved in the bean-rust interaction as reported here provides a comprehensive resource for other epigenomic regulation studies in non-model species under stress.

The complex process of formation of storage roots (SRs) from adventitious roots affects sweetpotato yield. Identifying the genes that are uniquely expressed in the SR forming cultivated species, Ipomoea batatas (Ib), and its immediate ancestral species, Ipomoea trifida (It), which does not form SRs, may provide insights into the molecular mechanisms underlying SR formation in sweetpotato. Illumina paired-end sequencing generated ~208 and ~200 million reads for Ib and It, respectively. Trinity assembly of the reads resulted in 98,317 transcripts for Ib and 275,044 for It, after post-assembly removal of trans-chimeras. From these sequences, we identified 4,865 orthologous genes in both Ib and It, 60 paralogous genes in Ib and 2,286 paralogous genes in It. Among paralogous gene sets, transcripts encoding the transcription factor RKD, which may have a role in nitrogen regulation and starch formation, and rhamnogalacturonate lyase (RGL) family proteins, which produce the precursors of cell wall polysaccharides, were found only in Ib. In addition, transcripts encoding a K+ efflux antiporter (KEA5) and the ERECTA protein kinase, which function in phytohormonal regulation and root proliferation, respectively, were also found only in Ib. qRT-PCR indicated that starch and sucrose metabolism genes, such as those encoding ADP-glucose pyrophosphorylase and beta-amylase, showed lower expression in It than Ib, whereas lignin genes such as caffeoyl-CoA O-methyltransferase (CoMT) and cinnamyl alcohol dehydrogenase (CAD) showed higher expression in It than Ib. A total of 7,067 and 9,650 unique microsatellite markers, 1,037,396 and 495,931 single nucleotide polymorphisms (SNPs) and 103,439 and 69,194 InDels in Ib and It, respectively, were also identified from this study. The detection of genes involved in the biosynthesis of RGL family proteins, the transcription factor RKD, and genes encoding a K+ efflux antiporter (KEA5) and the ERECTA protein kinase only in I. batatas indicate that these genes may have important functions in SR formation in sweetpotato. Potential molecular markers (SNPs, simple sequence repeats and InDels) and sequences identified in this study may represent a valuable resource for sweetpotato gene annotation and may serve as important tools for improving SR formation in sweetpotato through breeding.

A birth and death process is frequently used for modeling the size of a gene family that may vary along the branches of a phylogenetic tree. Under the birth and death model, maximum likelihood methods have been developed to estimate the birth and death rate and the sizes of ancient gene families (numbers of gene copies at the internodes of the phylogenetic tree). This paper aims to provide a Bayesian approach for estimating parameters in the birth and death model.

Common bean (Phaseolus vulgaris L.) is one of the most important legumes in the world. Several diseases severely reduce bean production and quality; therefore, it is very important to better understand disease resistance in common bean in order to prevent these losses. More than 70 resistance (R) genes which confer resistance against various pathogens have been cloned from diverse plant species. Most R genes share highly conserved domains which facilitates the identification of new candidate R genes from the same species or other species. The goals of this study were to isolate expressed R gene-like sequences (RGLs) from 454-derived transcriptomic sequences and expressed sequence tags (ESTs) of common bean, and to develop RGL-tagged molecular markers.

Key message Transcriptomes of two switchgrass genotypes representing the upland and lowland ecotypes will be key tools in switchgrass genome annotation and biotic and abiotic stress functional genomics. Abstract Switchgrass (Panicum virgatum L.) is an important bioenergy feedstock for cellulosic ethanol production. We report genome-wide transcriptome profiling of two contrasting tetraploid switchgrass genotypes, VS16 and AP13, representing the upland and lowland ecotypes, respectively. A total of 268 million Illumina short reads (50 nt) were generated, of which, 133 million were obtained in AP13 and the rest 135 million in VS16. More than 90% of these reads were mapped to the switchgrass reference genome (V1.1). We identified 6619 and 5369 differentially expressed genes in VS16 and AP13, respectively. Gene ontology and KEGG pathway analysis identified key genes that regulate important pathways including C4 photosynthesis, photorespiration and phenylpropanoid metabolism. A series of genes (33) involved in photosynthetic pathway were up-regulated in AP13 but only two genes showed higher expression in VS16. We identified three dicarboxylate transporter homologs that were highly expressed in AP13. Additionally, genes that mediate drought, heat, and salinity tolerance were also identified. Vesicular transport proteins, syntaxin and signal recognition particles were seen to be up-regulated in VS16. Analyses of selected genes involved in biosynthesis of secondary metabolites, plant–pathogen interaction, membrane transporters, heat, drought and salinity stress responses confirmed significant variation in the relative expression reflected in RNA-Seq data between VS16 and AP13 genotypes. The phenylpropanoid pathway genes identified here are potential targets for biofuel conversion.