Nociceptor neurons affect cancer immunosurveillance.

Solid tumours are innervated by nerve fibres that arise from the autonomic and sensory peripheral nervous systems1-5. Whether the neo-innervation of tumours by pain-initiating sensory neurons affects cancer immunosurveillance remains unclear. Here we show that melanoma cells interact with nociceptor neurons, leading to increases in their neurite outgrowth, responsiveness to noxious ligands and neuropeptide release. Calcitonin gene-related peptide (CGRP)-one such nociceptor-produced neuropeptide-directly increases the exhaustion of cytotoxic CD8+ T cells, which limits their capacity to eliminate melanoma. Genetic ablation of the TRPV1 lineage, local pharmacological silencing of nociceptors and antagonism of the CGRP receptor RAMP1 all reduced the exhaustion of tumour-infiltrating leukocytes and decreased the growth of tumours, nearly tripling the survival rate of mice that were inoculated with B16F10 melanoma cells. Conversely, CD8+ T cell exhaustion was rescued in sensory-neuron-depleted mice that were treated with local recombinant CGRP. As compared with wild-type CD8+ T cells, Ramp1-/- CD8+ T cells were protected against exhaustion when co-transplanted into tumour-bearing Rag1-deficient mice. Single-cell RNA sequencing of biopsies from patients with melanoma revealed that intratumoral RAMP1-expressing CD8+ T cells were more exhausted than their RAMP1-negative counterparts, whereas overexpression of RAMP1 correlated with a poorer clinical prognosis. Overall, our results suggest that reducing the release of CGRP from tumour-innervating nociceptors could be a strategy to improve anti-tumour immunity by eliminating the immunomodulatory effects of CGRP on cytotoxic CD8+ T cells.

Genome-wide DNA methylation dynamics following recent polyploidy in the allotetraploid Tragopogon miscellus (Asteraceae).

Polyploidy is an important evolutionary force, yet epigenetic mechanisms, such as DNA methylation, that regulate genome-wide expression of duplicated genes remain largely unknown. Here, we use Tragopogon (Asteraceae) as a model system to discover patterns and temporal dynamics of DNA methylation in recently formed polyploids. The naturally occurring allotetraploid Tragopogon miscellus formed in the last 95-100 yr from parental diploids Tragopogon dubius and T. pratensis. We profiled the DNA methylomes of these three species using whole-genome bisulfite sequencing. Genome-wide methylation levels in T. miscellus were intermediate between its diploid parents. However, nonadditive CG and CHG methylation occurred in transposable elements (TEs), with variation among TE types. Most differentially methylated regions (DMRs) showed parental legacy, but some novel DMRs were detected in the polyploid. Differentially methylated genes (DMGs) were also identified and characterized. This study provides the first assessment of both overall and locus-specific patterns of DNA methylation in a recent natural allopolyploid and shows that novel methylation variants can be generated rapidly after polyploid formation. Together, these results demonstrate that mechanisms to regulate duplicate gene expression may arise soon after allopolyploid formation and that these mechanisms vary among genes.

Natural variation in DNA methylation homeostasis and the emergence of epialleles.

In plants and mammals, DNA methylation plays a critical role in transcriptional silencing by delineating heterochromatin from transcriptionally active euchromatin. A homeostatic balance between heterochromatin and euchromatin is essential to genomic stability. This is evident in many diseases and mutants for heterochromatin maintenance, which are characterized by global losses of DNA methylation coupled with localized ectopic gains of DNA methylation that alter transcription. Furthermore, we have shown that genome-wide methylation patterns in Arabidopsis thaliana are highly stable over generations, with the exception of rare epialleles. However, the extent to which natural variation in the robustness of targeting DNA methylation to heterochromatin exists, and the phenotypic consequences of such variation, remain to be fully explored. Here we describe the finding that heterochromatin and genic DNA methylation are highly variable among 725 A. thaliana accessions. We found that genic DNA methylation is inversely correlated with that in heterochromatin, suggesting that certain methylation pathway(s) may be redirected to genes upon the loss of heterochromatin. This redistribution likely involves a feedback loop involving the DNA methyltransferase, CHROMOMETHYLASE 3 (CMT3), H3K9me2, and histone turnover, as highly expressed, long genes with a high density of CMT3-preferred CWG sites are more likely to be methylated. Importantly, although the presence of CG methylation in genes alone may not affect transcription, genes containing CG methylation are more likely to become methylated at non-CG sites and silenced. These findings are consistent with the hypothesis that natural variation in DNA methylation homeostasis may underlie the evolution of epialleles that alter phenotypes.

pENCODE: a plant encyclopedia of DNA elements.

ENCODE projects exist for many eukaryotes, including humans, but as of yet no defined project exists for plants. A plant ENCODE would be invaluable to the research community and could be more readily produced than its metazoan equivalents by capitalizing on the preexisting infrastructure provided from similar projects. Collecting and normalizing plant epigenomic data for a range of species will facilitate hypothesis generation, cross-species comparisons, annotation of genomes, and an understanding of epigenomic functions throughout plant evolution. Here, we discuss the need for such a project, outline the challenges it faces, and suggest ways forward to build a plant ENCODE.

Genome-Wide Reinforcement of DNA Methylation Occurs during Somatic Embryogenesis in Soybean.

Somatic embryogenesis is an important tissue culture technique that sometimes leads to phenotypic variation via genetic and/or epigenetic changes. To understand the genomic and epigenomic impacts of somatic embryogenesis, we characterized soybean (Glycine max) epigenomes sampled from embryos at 10 different stages ranging from 6 weeks to 13 years of continuous culture. We identified genome-wide increases in DNA methylation from cultured samples, especially at CHH sites. The hypermethylation almost exclusively occurred in regions previously possessing non-CG methylation and was accompanied by increases in the expression of genes encoding the RNA-directed DNA methylation (RdDM) machinery. The epigenomic changes were similar between somatic and zygotic embryogenesis. Following the initial global wave of hypermethylation, rare decay events of maintenance methylation were observed, and the extent of the decay increased with time in culture. These losses in DNA methylation were accompanied by downregulation of genes encoding the RdDM machinery and transcriptome reprogramming reminiscent of transcriptomes during late-stage seed development. These results reveal a process for reinforcing already silenced regions to maintain genome integrity during somatic embryogenesis over the short term, which eventually decays at certain loci over longer time scales.

Monitoring the interplay between transposable element families and DNA methylation in maize.

DNA methylation and epigenetic silencing play important roles in the regulation of transposable elements (TEs) in many eukaryotic genomes. A majority of the maize genome is derived from TEs that can be classified into different orders and families based on their mechanism of transposition and sequence similarity, respectively. TEs themselves are highly methylated and it can be tempting to view them as a single uniform group. However, the analysis of DNA methylation profiles in flanking regions provides evidence for distinct groups of chromatin properties at different TE families. These differences among TE families are reproducible in different tissues and different inbred lines. TE families with varying levels of DNA methylation in flanking regions also show distinct patterns of chromatin accessibility and modifications within the TEs. The differences in the patterns of DNA methylation flanking TE families arise from a combination of non-random insertion preferences of TE families, changes in DNA methylation triggered by the insertion of the TE and subsequent selection pressure. A set of nearly 70,000 TE polymorphisms among four assembled maize genomes were used to monitor the level of DNA methylation at haplotypes with and without the TE insertions. In many cases, TE families with high levels of DNA methylation in flanking sequence are enriched for insertions into highly methylated regions. The majority of the >2,500 TE insertions into unmethylated regions result in changes in DNA methylation in haplotypes with the TE, suggesting the widespread potential for TE insertions to condition altered methylation in conserved regions of the genome. This study highlights the interplay between TEs and the methylome of a major crop species.

Subgenome Dominance in an Interspecific Hybrid, Synthetic Allopolyploid, and a 140-Year-Old Naturally Established Neo-Allopolyploid Monkeyflower.

Recent studies have shown that one of the parental subgenomes in ancient polyploids is generally more dominant, having retained more genes and being more highly expressed, a phenomenon termed subgenome dominance. The genomic features that determine how quickly and which subgenome dominates within a newly formed polyploid remain poorly understood. To investigate the rate of emergence of subgenome dominance, we examined gene expression, gene methylation, and transposable element (TE) methylation in a natural, <140-year-old allopolyploid (Mimulus peregrinus), a resynthesized interspecies triploid hybrid (M. robertsii), a resynthesized allopolyploid (M. peregrinus), and progenitor species (M. guttatus and M. luteus). We show that subgenome expression dominance occurs instantly following the hybridization of divergent genomes and significantly increases over generations. Additionally, CHH methylation levels are reduced in regions near genes and within TEs in the first-generation hybrid, intermediate in the resynthesized allopolyploid, and are repatterned differently between the dominant and recessive subgenomes in the natural allopolyploid. Subgenome differences in levels of TE methylation mirror the increase in expression bias observed over the generations following hybridization. These findings provide important insights into genomic and epigenomic shock that occurs following hybridization and polyploid events and may also contribute to uncovering the mechanistic basis of heterosis and subgenome dominance.

Changes of gene expression but not cytosine methylation are associated with male parental care reflecting behavioural state, social context and individual flexibility.

Behaviour is often a front line response to changing environments. Recent studies show behavioural changes are associated with changes of gene expression; however, these studies have primarily focused on discrete behavioural states. We build on these studies by addressing additional contexts that produce qualitatively similar behavioural changes. We measured levels of gene expression and cytosine methylation, which is hypothesized to regulate the transcriptional architecture of behavioural transitions, within the brain during male parental care of the burying beetle Nicrophorus vespilloides in a factorial design. Male parenting is a suitably plastic behaviour because although male N. vespilloides typically do not provide direct care (i.e. feed offspring) when females are present, levels of feeding by a male equivalent to the female can be induced by removing the female. We examined three different factors: behavioural state (caring versus non-caring), social context (with or without a female mate) and individual flexibility (if a male switched to direct care after his mate was removed). The greatest number of differentially expressed genes were associated with behavioural state, followed by social context and individual flexibility. Cytosine methylation was not associated with changes of gene expression in any of the factors. Our results suggest a hierarchical association between gene expression and the different factors, but that this process is not controlled by cytosine methylation. Our results further suggest that the extent a behaviour is transient plays an underappreciated role in determining its underpinning molecular mechanisms.

A novel requirement for DROSHA in maintenance of mammalian CG methylation.

In mammals, faithful inheritance of genomic methylation patterns ensures proper gene regulation and cell behaviour, impacting normal development and fertility. Following establishment, genomic methylation patterns are transmitted through S-phase by the maintenance methyltransferase Dnmt1. Using a protein interaction screen, we identify Microprocessor component DROSHA as a novel DNMT1-interactor. Drosha-deficient embryonic stem (ES) cells display genomic hypomethylation that is not accounted for by changes in the levels of DNMT proteins. DNMT1-mediated methyltransferase activity is also reduced in these cells. We identify two transcripts that are specifically upregulated in Drosha- but not Dicer-deficient ES cells. Regions within these transcripts predicted to form stem-loop structures are processed by Microprocessor and can inhibit DNMT1-mediated methylation in vitro. Our results highlight DROSHA as a novel regulator of mammalian DNA methylation and we propose that DROSHA-mediated processing of RNA is necessary to ensure full DNMT1 activity. This adds to the DROSHA repertoire of non-miRNA dependent functions as well as implicating RNA in regulating DNMT1 activity and correct levels of genomic methylation.

Rational development of transformation in Clostridium thermocellum ATCC 27405 via complete methylome analysis and evasion of native restriction-modification systems.

A major barrier to both metabolic engineering and fundamental biological studies is the lack of genetic tools in most microorganisms. One example is Clostridium thermocellum ATCC 27405T, where genetic tools are not available to help validate decades of hypotheses. A significant barrier to DNA transformation is restriction-modification systems, which defend against foreign DNA methylated differently than the host. To determine the active restriction-modification systems in this strain, we performed complete methylome analysis via single-molecule, real-time sequencing to detect 6-methyladenine and 4-methylcytosine and the rarely used whole-genome bisulfite sequencing to detect 5-methylcytosine. Multiple active systems were identified, and corresponding DNA methyltransferases were expressed from the Escherichia coli chromosome to mimic the C. thermocellum methylome. Plasmid methylation was experimentally validated and successfully electroporated into C. thermocellum ATCC 27405. This combined approach enabled genetic modification of the C. thermocellum-type strain and acts as a blueprint for transformation of other non-model microorganisms.

The LSH/DDM1 Homolog MUS-30 Is Required for Genome Stability, but Not for DNA Methylation in Neurospora crassa.

LSH/DDM1 enzymes are required for DNA methylation in higher eukaryotes and have poorly defined roles in genome maintenance in yeast, plants, and animals. The filamentous fungus Neurospora crassa is a tractable system that encodes a single LSH/DDM1 homolog (NCU06306). We report that the Neurospora LSH/DDM1 enzyme is encoded by mutagen sensitive-30 (mus-30), a locus identified in a genetic screen over 25 years ago. We show that MUS-30-deficient cells have normal DNA methylation, but are hypersensitive to DNA damaging agents. MUS-30 is a nuclear protein, consistent with its predicted role as a chromatin remodeling enzyme, and levels of MUS-30 are increased following DNA damage. MUS-30 co-purifies with Neurospora WDR76, a homolog of yeast Changed Mutation Rate-1 and mammalian WD40 repeat domain 76. Deletion of wdr76 rescued DNA damage-hypersensitivity of Δmus-30 strains, demonstrating that the MUS-30-WDR76 interaction is functionally important. DNA damage-sensitivity of Δmus-30 is partially suppressed by deletion of methyl adenine glycosylase-1, a component of the base excision repair machinery (BER); however, the rate of BER is not affected in Δmus-30 strains. We found that MUS-30-deficient cells are not defective for DSB repair, and we observed a negative genetic interaction between Δmus-30 and Δmei-3, the Neurospora RAD51 homolog required for homologous recombination. Together, our findings suggest that MUS-30, an LSH/DDM1 homolog, is required to prevent DNA damage arising from toxic base excision repair intermediates. Overall, our study provides important new information about the functions of the LSH/DDM1 family of enzymes.

On the origin and evolutionary consequences of gene body DNA methylation.

In plants, CG DNA methylation is prevalent in the transcribed regions of many constitutively expressed genes (gene body methylation; gbM), but the origin and function of gbM remain unknown. Here we report the discovery that Eutrema salsugineum has lost gbM from its genome, to our knowledge the first instance for an angiosperm. Of all known DNA methyltransferases, only CHROMOMETHYLASE 3 (CMT3) is missing from E. salsugineum Identification of an additional angiosperm, Conringia planisiliqua, which independently lost CMT3 and gbM, supports that CMT3 is required for the establishment of gbM. Detailed analyses of gene expression, the histone variant H2A.Z, and various histone modifications in E. salsugineum and in Arabidopsis thaliana epigenetic recombinant inbred lines found no evidence in support of any role for gbM in regulating transcription or affecting the composition and modification of chromatin over evolutionary timescales.

Genome-wide redistribution of H3K27me3 is linked to genotoxic stress and defective growth.

H3K9 methylation directs heterochromatin formation by recruiting multiple heterochromatin protein 1 (HP1)-containing complexes that deacetylate histones and methylate cytosine bases in DNA. In Neurospora crassa, a single H3K9 methyltransferase complex, called the DIM-5,-7,-9, CUL4, DDB1 Complex (DCDC), is required for normal growth and development. DCDC-deficient mutants are hypersensitive to the genotoxic agent methyl methanesulfonate (MMS), but the molecular basis of genotoxic stress is unclear. We found that both the MMS sensitivity and growth phenotypes of DCDC-deficient strains are suppressed by mutation of embryonic ectoderm development or Su-(var)3-9; E(z); Trithorax (set)-7, encoding components of the H3K27 methyltransferase Polycomb repressive complex-2 (PRC2). Trimethylated histone H3K27 (H3K27me3) undergoes genome-wide redistribution to constitutive heterochromatin in DCDC- or HP1-deficient mutants, and introduction of an H3K27 missense mutation is sufficient to rescue phenotypes of DCDC-deficient strains. Accumulation of H3K27me3 in heterochromatin does not compensate for silencing; rather, strains deficient for both DCDC and PRC2 exhibit synthetic sensitivity to the topoisomerase I inhibitor Camptothecin and accumulate γH2A at heterochromatin. Together, these data suggest that PRC2 modulates the response to genotoxic stress.

Maize death acids, 9-lipoxygenase-derived cyclopente(a)nones, display activity as cytotoxic phytoalexins and transcriptional mediators.

Plant damage promotes the interaction of lipoxygenases (LOXs) with fatty acids yielding 9-hydroperoxides, 13-hydroperoxides, and complex arrays of oxylipins. The action of 13-LOX on linolenic acid enables production of 12-oxo-phytodienoic acid (12-OPDA) and its downstream products, termed "jasmonates." As signals, jasmonates have related yet distinct roles in the regulation of plant resistance against insect and pathogen attack. A similar pathway involving 9-LOX activity on linolenic and linoleic acid leads to the 12-OPDA positional isomer, 10-oxo-11-phytodienoic acid (10-OPDA) and 10-oxo-11-phytoenoic acid (10-OPEA), respectively; however, physiological roles for 9-LOX cyclopentenones have remained unclear. In developing maize (Zea mays) leaves, southern leaf blight (Cochliobolus heterostrophus) infection results in dying necrotic tissue and the localized accumulation of 10-OPEA, 10-OPDA, and a series of related 14- and 12-carbon metabolites, collectively termed "death acids." 10-OPEA accumulation becomes wound inducible within fungal-infected tissues and at physiologically relevant concentrations acts as a phytoalexin by suppressing the growth of fungi and herbivores including Aspergillus flavus, Fusarium verticillioides, and Helicoverpa zea. Unlike previously established maize phytoalexins, 10-OPEA and 10-OPDA display significant phytotoxicity. Both 12-OPDA and 10-OPEA promote the transcription of defense genes encoding glutathione S transferases, cytochrome P450s, and pathogenesis-related proteins. In contrast, 10-OPEA only weakly promotes the accumulation of multiple protease inhibitor transcripts. Consistent with a role in dying tissue, 10-OPEA application promotes cysteine protease activation and cell death, which is inhibited by overexpression of the cysteine protease inhibitor maize cystatin-9. Unlike jasmonates, functions for 10-OPEA and associated death acids are consistent with specialized roles in local defense reactions.

Base-resolution detection of N4-methylcytosine in genomic DNA using 4mC-Tet-assisted-bisulfite- sequencing.

Restriction-modification (R-M) systems pose a major barrier to DNA transformation and genetic engineering of bacterial species. Systematic identification of DNA methylation in R-M systems, including N(6)-methyladenine (6mA), 5-methylcytosine (5mC) and N(4)-methylcytosine (4mC), will enable strategies to make these species genetically tractable. Although single-molecule, real time (SMRT) sequencing technology is capable of detecting 4mC directly for any bacterial species regardless of whether an assembled genome exists or not, it is not as scalable to profiling hundreds to thousands of samples compared with the commonly used next-generation sequencing technologies. Here, we present 4mC-Tet-assisted bisulfite-sequencing (4mC-TAB-seq), a next-generation sequencing method that rapidly and cost efficiently reveals the genome-wide locations of 4mC for bacterial species with an available assembled reference genome. In 4mC-TAB-seq, both cytosines and 5mCs are read out as thymines, whereas only 4mCs are read out as cytosines, revealing their specific positions throughout the genome. We applied 4mC-TAB-seq to study the methylation of a member of the hyperthermophilc genus, Caldicellulosiruptor, in which 4mC-related restriction is a major barrier to DNA transformation from other species. In combination with MethylC-seq, both 4mC- and 5mC-containing motifs are identified which can assist in rapid and efficient genetic engineering of these bacteria in the future.

Study of seed hair growth in Populus tomentosa, an important character of female floral bud development.

BACKGROUND: Poplar seed hair is an environmental annoyance in northern China due to its abundance and widespread airborne distribution after maturation. The morphogenesis and molecular mechanisms of its development are not well understood, and little attention has been focused on the dynamics of its development. To better understand the mechanism of poplar seed hair development, paraffin sections were used to examine the initiation and elongation of poplar seed hairs. RNA-seq technology was also employed to provide a comprehensive overview of transcriptional changes that occur during seed hair development. RESULTS: The placenta at the base of ovary, was identified as the origin of seed hair development, which is in sharp contrast to cotton fibers that originate from epidermal cells of the seed coat. An enlarged cell nucleus in seed hair cells was also observed, which was supported by our gene ontology enrichment analysis. The significant enriched GO term of "endoreduplication" indicated that cycles of endoreduplication, bypassing normal mitosis, is the underlying mechanisms for the maintenance of the uni-cellular structure of seed hairs. By analyzing global changes in the transcriptome, many genes regulating cell cycle, cell elongation, cell well modification were identified. Additionally, in an analysis of differential expression, cellulose synthesis and cell wall biosynthesis-related biological processes were enriched, indicating that this component of fiber structure in poplar seed hairs is consistent with what is found in cotton fibers. Differentially expressed transcription factors exhibited a stage-specific up-regulation. A dramatic down-regulation was also revealed during the mid-to-late stage of poplar seed hair development, which may point to novel mechanisms regulating cell fate determination and cell elongation. CONCLUSIONS: This study revealed the initiation site of poplar seed hairs and also provided a comprehensive overview of transcriptome dynamics during the process of seed hair development. The high level of resolution on dynamic changes in the transcriptome provided in this study may serve as a valuable resource for developing a more complete understanding of this important biological process.

Identification of glutathione S-transferase genes responding to pathogen infestation in Populus tomentosa.

Stem blister canker, caused by Botryosphaeria dothidea, is becoming the most serious disease of poplar in China. The molecular basis of the poplar in response to stem blister canker is not well understood. To reveal the global transcriptional changes of poplar to infection by B. dothidea, Solexa paired-end sequencing of complementary DNAs (cDNAs) from control (NB) and pathogen-treated samples (WB) was performed, resulting in a total of 339,283 transcripts and 183,881 unigenes. A total of 206,586 transcripts were differentially expressed in response to pathogen stress (false discovery rate ≤0.05 and an absolute value of log2Ratio (NB/WB) ≥1). In enrichment analysis, energy metabolism and redox reaction-related macromolecules were accumulated significantly in Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analyses, indicating components of dynamic defense against the fungus. A total of 852 transcripts (575 upregulated and 277 downregulated transcripts) potentially involved in plant-pathogen interaction were also differentially regulated, including genes encoding proteins linked to signal transduction (putative leucine-rich repeat (LRR) protein kinases and calcium-binding proteins), defense (pathogenesis-related protein 1), and cofactors (jasmonate-ZIM-domain-containing proteins and heat shock proteins). Moreover, transcripts encoding glutathione S-transferase (GST) were accumulated to high levels, revealing key genes and proteins potentially related to pathogen resistance. Poplar RNA sequence data were validated by quantitative real-time PCR (RT-qPCR), which revealed a highly reliability of the transcriptomic profiling data.

A Novel Moderate Constitutive Promoter Derived from Poplar (Populus tomentosa Carrière).

A novel sequence that functions as a promoter element for moderate constitutive expression of transgenes, designated as the PtMCP promoter, was isolated from the woody perennial Populus tomentosa. The PtMCP promoter was fused to the GUS reporter gene to characterize its expression pattern in different species. In stable Arabidopsis transformants, transcripts of the GUS reporter gene could be detected by RT-PCR in the root, stem, leaf, flower and silique. Further histochemical and fluorometric GUS activity assays demonstrated that the promoter could direct transgene expression in all tissues and organs, including roots, stems, rosette leaves, cauline leaves and flowers of seedlings and maturing plants. Its constitutive expression pattern was similar to that of the CaMV35S promoter, but the level of GUS activity was significantly lower than in CaMV35S promoter::GUS plants. We also characterized the promoter through transient expression in transgenic tobacco and observed similar expression patterns. Histochemical GUS staining and quantitative analysis detected GUS activity in all tissues and organs of tobacco, including roots, stems, leaves, flower buds and flowers, but GUS activity in PtMCP promoter::GUS plants was significantly lower than in CaMV35S promoter::GUS plants. Our results suggested that the PtMCP promoter from poplar is a constitutive promoter with moderate activity and that its function is presumably conserved in different species. Therefore, the PtMCP promoter may provide a practical choice to direct moderate level constitutive expression of transgenes and could be a valuable new tool in plant genetic engineering.

Identification and characterization of the Populus AREB/ABF subfamily.

Abscisic acid (ABA) is a major plant hormone that plays an important role in responses to abiotic stresses. The ABA-responsive element binding protein/ABRE-binding factor (AREB/ABF) gene subfamily contains crucial transcription factors in the ABA-mediated signaling pathway. In this study, a total of 14 putative AREB/ABF members were identified in the Populus trichocarpa Torr. & Gray. genome using five AREB/ABF amino acid sequences from Arabidopsis thaliana L. as probes. The 14 putative Populus subfamily members showed high protein similarities, especially in the basic leucine zipper (bZIP) domain region. A neighbor-joining analysis combined with gene structure data revealed homology among the 14 genes. The expression patterns of the Populus AREB/ABF subfamily suggested that the most abundant transcripts of 11 genes occurred in leaf tissues, while two genes were most transcribed in root tissues. Significantly, eight Populus AREB/ABF gene members were upregulated after treatment with 100 µM exogenous ABA, while the other six members were downregulated. We identified the expression profiles of the subfamily members in Populus tissues and elucidated different response patterns of Populus AREB/ABF members to ABA stress. This study provided insight into the roles of Populus AREB/ABF homologues in plant response to abiotic stresses.