Genomic insights into salt adaptation in a desert poplar

Abstract

Despite the high economic and ecological importance of forests, our knowledge of the genomic evolution of trees under salt stress remains very limited. Here we report the genome sequence of the desert poplar, Populus euphratica, which exhibits high tolerance to salt stress. Its genome is very similar and collinear to that of the closely related mesophytic congener, P. trichocarpa. However, we find that several gene families likely to be involved in tolerance to salt stress contain significantly more gene copies within the P. euphratica lineage. Furthermore, genes showing evidence of positive selection are significantly enriched in functional categories related to salt stress. Some of these genes, and others within the same categories, are significantly upregulated under salt stress relative to their expression in another salt-sensitive poplar. Our results provide an important background for understanding tree adaptation to salt stress and facilitating the genetic improvement of cultivated poplars for saline soils.

ARTICLE

Received 1 Sep 2013 | Accepted 21 Oct 2013 | Published 21 Nov 2013

Tao Ma1,*, Junyi Wang2,*, Gongke Zhou3,*, Zhen Yue2, Quanjun Hu1, Yan Chen2, Bingbing Liu1, Qiang Qiu1, Zhuo Wang2, Jian Zhang1, Kun Wang1, Dechun Jiang1, Caiyun Gou2, Lili Yu2, Dongliang Zhan2, Ran Zhou1, Wenchun Luo1, Hui Ma1, Yongzhi Yang1, Shengkai Pan2, Dongming Fang2, Yadan Luo2, Xia Wang1, Gaini Wang1, Juan Wang1, Qian Wang1, Xu Lu1, Zhe Chen2, Jinchao Liu2, Yao Lu2, Ye Yin2, Huanming Yang2, Richard J. Abbott4, Yuxia Wu1, Dongshi Wan1, Jia Li1, Tongming Yin5, Martin Lascoux6, Stephen P. DiFazio7, Gerald A. Tuskan8, Jun Wang2,9 & Liu Jianquan1

Despite the high economic and ecological importance of forests, our knowledge of the genomic evolution of trees under salt stress remains very limited. Here we report the genome sequence of the desert poplar, Populus euphratica, which exhibits high tolerance to salt stress. Its genome is very similar and collinear to that of the closely related mesophytic congener,P. trichocarpa. However, we nd that several gene families likely to be involved in tolerance to salt stress contain signicantly more gene copies within the P. euphratica lineage. Furthermore, genes showing evidence of positive selection are signicantly enriched in functional categories related to salt stress. Some of these genes, and others within the same categories, are signicantly upregulated under salt stress relative to their expression in another salt-sensitive poplar. Our results provide an important background for understanding tree adaptation to salt stress and facilitating the genetic improvement of cultivated poplars for saline soils.

1 State Key Laboratory of Grassland Agro-Ecosystem, School of Life Sciences, Lanzhou University, Lanzhou 730000, China. 2 BGI-Shenzhen, Shenzhen 518083, China. 3 Key Laboratory of Biofuels and Shandong Provincial Key Laboratory of Energy Genetics, Qingdao Institute of Bioenergy and Bioprocess Technology, Chinese Academy of Sciences, Qingdao 266101, China. 4 School of Biology, Mitchell Building, University of St Andrews, St Andrews, Fife KY16 9TH, UK. 5 The Key Lab of Forest Genetics and Gene Engineering, Nanjing Forestry University, Nanjing 210037, China. 6 Department of Ecology and Genetics, Evolutionary Biology Centre, Uppsala University, Norbyvagen, 18D 75326 Uppsala, Sweden. 7 Department of Biology, West Virginia University, Morgantown, West Virginia 26506-6057, USA. 8 BioSciences Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA. 9 Department of Biology, University of Copenhagen, Copenhagen 1017, Denmark. * These authors contributed equally to this work. Correspondence and requests for materials should be addressed to J.L. (email: mailto:[email protected]

Web End [email protected] ) or to J.W. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 4:2797 | DOI: 10.1038/ncomms3797 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

& 2013 Macmillan Publishers Limited. All rights reserved.

DOI: 10.1038/ncomms3797

Genomic insights into salt adaptation in a desert poplar

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms3797

Forests dominate much of the terrestrial landscape1. However, forest trees rarely occur on saline soils and little is known of the genetic basis of their tolerance to salt stress

despite strong demand for their cultivation on highly saline soils in many parts of the world2. Members of the genus Populus are used as a model forest species for diverse studies not only because of their amenability to experimental and genetic manipulation, but also because of their high economic and ecological importance as the most widely cultivated tree throughout the northern hemisphere3,4. More than 30 wild Populus species occur across diverse habitats over a wide geographical range, thereby providing an excellent system for unravelling the genetic bases of adaptive divergence4. Populus euphratica Oliv., which is native to desert regions ranging from western China to North Africa, is characterized by extraordinary adaptation to salt stress58. Notably, at high salinity it maintains higher growth and photosynthetic rates than other poplar species9,10 and can survive concentrations of NaCl in nutrient solution up to 450 mM11.

In this study, we examine genomic differences between a xeric desert poplar and its mesophytic congener, P. trichocarpa, for which a high-quality reference genome is available12. We further examine gene expression differences following salt stress treatment in a comparison with another salt-sensitive congener,P. tomentosa. Our comparisons highlight the genetic bases of salt tolerance in the desert poplar.

ResultsGenome sequencing and assembly. Because of the limitations of next-generation sequencing for complex genome assembly13 and the high levels of polymorphism found in this non-domesticated and open-pollinated species (Supplementary Fig. S1), we employed a newly developed fosmid-pooling strategy14 to sequence and assemble the P. euphratica genome (Table 1 and Supplementary Methods). Hierarchical assembly using 67.1 Gb (B112 ) whole-

genome shotgun reads (Supplementary Table S1), combined with more than 200 high-quality reads from 66,240 fosmid clones

(Supplementary Table S2), yielded a nal assembly with a total length of 496.5 Mb (Supplementary Table S3), representing 83.7% of the P. euphratica nucleotide space (Supplementary Tables S4 and S5). The contig N50 of the assembled sequence was 40.4 Kb (longest, 728.4 Kb) and scaffold N50 was 482 Kb (longest, 8.8 Mb; Table 1), which were comparable to those of other plant genome assemblies generated by next-generation sequencing technology (Supplementary Table S6). Sequencing depth distribution showed that over 92.5% of the assembly was covered by more than 20

(Supplementary Figs S2 and S3), ensuring a high single-base accuracy. The heterozygosity level in P. euphratica was B0.5%

(Supplementary Tables S7 and S8, and Supplementary Fig. S4), which is almost twice that in P. trichocarpa (0.26%)12. The assembly covered 97.3% of the 516,712 Populus expressed sequence tags (Supplementary Table S9) and 97.7% of the 7 complete fosmids sequenced by Sanger sequencing (Supplementary Table S10 and Supplementary Fig. S5), without any obvious misassembly occurring. The coverage of the core eukaryotic genes was estimated

to be 94.35% for the P. euphratica assembly (Supplementary Table S11), which is comparable to the estimate for P. trichocarpa(93.95%). All of these statistics supported that our draft genome sequence has high contiguity, coverage and accuracy, further demonstrating the feasibility of this hierarchical approach for de novo sequencing and assembly of a complex genome with high heterozygosity14.

Genome annotation. Using a combination of homology-based searches and de novo annotation, we found that B44% of the

P. euphratica genome is composed of repetitive elements (Supplementary Table S12), similar to that of the P. trichocarpa genome (47%; Fig. 1). Long-terminal repeats (LTRs) were the most abundant repeat class, representing 36.7% and 33.1% of theP. euphratica and P. trichocarpa genomes, respectively (Supplementary Table S13). The distribution of repeat divergence rates revealed a peak of gypsy LTR at 11% in P. euphratica (Supplementary Fig. S6), which is likely to reect a relatively recent expansion of this LTR family in the P. euphratica lineage.

A total of 34,279 protein-coding genes were predicted to be present in the P. euphratica genome (Supplementary Table S14 and Supplementary Fig. S7), 96.6% of which were supported by expressed sequence tags and/or homology-based searching with only 3.4% derived solely from de novo gene predictions (Supplementary Fig. S8). Functional annotation conrmed that94.3% of the predicted genes had known homologues in protein databases (Supplementary Table S15). Small RNA sequencing data supported the occurrence of 152 conserved and 114 candidate novel microRNAs predicted from the P. euphratica genome (Supplementary Table S16, Supplementary Data 1 and 2, and Supplementary Figs S9S11), most of which were extensively up/downregulated in response to salt stress (Supplementary Table S17 and Supplementary Fig. S12). In addition, we also identied 764 transfer RNAs, 706 ribosomal RNAs and 4,826 small nuclear RNAs (Supplementary Table S18).

Genome evolution. In accordance with previous research12,15, the distribution of the fourfold degenerate synonymous sites of the third codons (4DTv) value between duplicated genes showed similar peaks (B0.09 and B0.59) in both P. euphratica and

P. trichocarpa genomes, suggesting that two ancient whole-genome duplication (WGD) events had occurred in the Populus lineage (Supplementary Table S19, Supplementary Figs S13 and S14). These shared WGDs were also conrmed by the extensive collinearity between the genomes of both species (Fig. 1). A total of 1,214 collinear blocks 410 kb in length, corresponding to 323 Mb (65% of the assembly) and 332 Mb (76%) in theP. euphratica and P. trichocarpa genomes, respectively, were identied (Supplementary Table S20). Assuming that the recent WGD occurred around 65 million years ago (Mya)12, divergence between P. euphratica and P. trichocarpa can be placed to B14

Mya (4DTv, B0.02), which approximates to that estimated from phylogenetic analysis (B8 Mya; Supplementary Fig. S15).

We identied and designed a total of 18,938 universal pairs of simple sequence repeat primers in the collinear regions, which

Table 1 | Statistics for the assembly of the Populus euphratica genome.

Assembly Contig N50* (bp) Longest contig (bp) Scaffold N50* (bp) Longest scaffold (bp) Normal whole-genome shotgun reads only 5,209 97,386 40,276 744,683

Fosmid reads only 6,519 46,382 10,282 69,305 Combined 40,438 728,449 482,055 8,759,900

*N50 refers to the size above which 50% of the total length of the sequence assembly can be found.

2 NATURE COMMUNICATIONS | 4:2797 | DOI: 10.1038/ncomms3797 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2013 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms3797 ARTICLE

Chr. 01 Chr. 02 Chr. 03 Chr. 04 Chr. 05 Chr. 06 Chr. 07 Chr. 08 Chr. 09 Chr. 10

Chr. 11 Chr. 12 Chr. 13 Chr. 14 Chr. 15 Chr. 16 Chr. 17 Chr. 18 Chr. 19

DNA

LINE

LTR/Copia

LTR/Gypsy

5 Mb

Figure 1 | Collinearity between the P. euphratica and P. trichocarpa genomes. The P. euphratica scaffolds (blue) inferred to be collinear are linked (grey lines) to P. trichocarpa chromosomes (orange). The proportion of the repeat elements (left) across the chromosomes is indicated for 400 kb sliding windows at 50 kb steps. DNA transposable elements are shown in purple, long interspersed elements (LINE) are shown in light blue, Copia and Gypsy elements of LTR retrotransposons are shown in yellow and green, respectively.

can be converted into genetic markers across most poplar species (Supplementary Data 3 and Supplementary Table S21). These simple sequence repeat markers, as well as the intraspecic or interspecic nucleotide variation in collinear regions, will facilitate genetic dissection of agronomically important traits and accelerate the genetic improvement of cultivated poplars, particularly for growth on saline soils.

Adaptation to a saline environment. Copy number within gene families has been reported to vary greatly between closely related, divergent species16,17. Both gene family and InterProScan domain analysis revealed that several gene families related to salt stress were substantially expanded in P. euphratica compared with other plant species (Fig. 2a, Supplementary Tables S22S24, and Supplementary Figs S16 and S17). For example, the HKT1 (high-afnity K transporter 1) gene family, which encodes Na /K transporters that have important roles in affecting or determining salt tolerance in plants18, expanded from one member in theP. trichocarpa genome to four in the P. euphratica genome. Three of these genes occurred as tandem duplicates in P. euphratica, which together corresponded to a HKT1-like pseudogene inP. trichocarpa (Fig. 2b). HKT1 transporters have a key role in limiting Na transport from roots to shoots in Arabidopsis19, and may account for the lower rate of ion uptake and transport recorded in P. euphratica9,20,21. The gene family encoding P-type H-ATPases also had more copies in the P. euphratica genome than in the P. trichocarpa genome (Fig. 2c). These P-type

ATPases provide the basic energy for Na /H antiporters by sustaining an electrochemical H gradient across the plasma membrane, thus making an important contribution to maintenance of low Na concentrations in P. euphratica22,23.

Other expanded gene families include those encoding antioxidative enzymes24, such as CAT (catalase) and GR1 (glutathione reductase), and genes involved in abscisic acid (ABA) signalling regulation25, such as GCR2 (G-protein-coupled receptor 2) and PLD (phospholipase D). Heat-shock proteins usually protect cells against salinity by controlling the proper folding and conformation of proteins26. Several of these families (for example, HSP20, HSP70 and HSP90) were expanded in the P. euphratica genome. In addition, the P. euphratica genome has

more copies of BADH (betaine aldehyde dehydrogenase) and GolS4 (galactinol synthase 4), which encode key enzymes involved in biosynthesis of critical solutes that have roles in osmotic adjustment pathways under salt stress27,28.

Adaptive divergence at the molecular level may also be reected by an increased rate of non-synonymous changes within genes involved in adaptation29. In collinear regions, we identied 18,262 high-condence 1:1 orthologous genes betweenP. euphratica and P. trichocarpa, with a mean protein similarity close to 98.94% (Supplementary Fig. S18). The genes with elevated pairwise genetic differentiation were primarily enriched in photosynthetic electron transport chain, heat acclimation, oxidoreductase activity and cation channel activity (Supplementary Table S25), indicating rapid evolution and/or adaptive divergence in these functions between P. euphratica andP. trichocarpa. Of the 6,545 high-condence orthologues identied among 10 plant species (Supplementary Fig. S16), we detected 57 positively selected genes (PSGs) in the P. euphratica lineage (Supplementary Table S26), which is signicantly greater than the number (35 PSGs) in the P. trichocarpa lineage (P-value 0.014 by the Fishers exact test). Compared with

P. trichocarpa PSGs, P. euphratica PSGs were signicantly enriched (P-value r0.05 by the Fishers exact test) in response to stimulus, cation binding and oxidoreductase activity (Fig. 2d). They included ENH1 (enhancer of sos3-1), which encodes a chloroplast-localized rubredoxin-like protein and has an important role in mediation of both ion homeostasis and reactive oxygen species detoxication30; CIPK1 (CBL-interacting protein kinase 1), a protein kinase interacting strongly with the calcium sensors CBL1 and CBL9, and alternatively controlling ABA-dependent and ABA-independent stress responses in Arabidopsis31; and PSD1 (phosphatidylserine decarboxylase 1) encoding a crucial enzyme catalysing production of phosphatidylethanolamine and therefore raising stress tolerance by increasing the exibility of cell membranes32 (Fig. 2c). Several genes encoding transcription factors such as HB40, bHLH87 and AP2/ERF, and oxidoreductases such as peroxidase, 2-oxoglutarate and Fe(II)-dependent oxygenase, also showed signs of positive selection.

To examine the genome-wide responses to salt stress of this desert poplar, we performed a series of deep transcriptome

NATURE COMMUNICATIONS | 4:2797 | DOI: 10.1038/ncomms3797 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

& 2013 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms3797

5

199

P. euphraticaP. trichocarpa

CpaHKT1;1

Chr. 18

<>><<<<>>><

4

PpeHKT1;1

<

>

<

< >

166

FveHKT1;1

% Expanded genes

PpeHKT1;2

3

< >

<<

<

<

<

<

PtrHKT1-like

146

156

2

PeuHKT1;3

PeuHKT1;4

PeuHKT1;2

PeuHKT1;1

PtrHKT1;1

75

1

24

2

11 3

82

0.5 kb

Ion

transport

Cell wall organization

AthHKT1;1 TpaHKT1;1 TpaHKT1;2 OsaHKT1;1

0

20 kb

Carbohydrate

metabolism

ATPase

activity

Oxidoreductase

activity

Scaffold31.1

21 21 23 25

5

3

ABA

P. euphraticaP. trichocarpa

GCR2

BCH1

4

PSD1

PLD

ZEP

HSPs

ROS

3

% PSGs

MYB GR1

ERF

Nucleus

BADH

GolS4

CAT

bZIP

2

Peroxidase

bHLH

CIPK1

WRKY

8

1

57

35

7

Salt stress

13

[Ca2+]

Vacuole

H+

1

2

5

0

0

Na+

SOS2

NHX1

Genome background

Response to stimulus

Cation binding

Oxidoreductase

activity

Unknown

SOS3

ENH1

HB

Na+

SOS1

H+

H+

ATP ADP+Pi

HKT1

H-ATPase

Na+

Figure 2 | Adaptation of P. euphratica to salt stress. (a) Comparison of the proportions of expanded genes in the P. euphratica and P. trichocarpa lineages relative to their common ancestor. The number of expanded genes of each class is indicated above each bar. (b) Tandem duplications of HKT1 genes. Note that the PtrHKT1-like gene is pseudogenized in P. trichocarpa. (c) PSGs and expanded key genes in salt-stress response pathways of P. euphratica. Boxes with borders indicate PSGs (red) and expanded (black) P. euphratica genes, and the lled colors correspond to their degree of regulation in FPKMtreatment/

FPKMcontrol in response to salt stress. (d) Comparison of the proportions of PSGs in the P. euphratica and P. trichocarpa lineages. The number of PSGs of each class is indicated above each bar. FPKM, fragments per kilobase of exon per million fragments mapped.

sequencings (Supplementary Table S27) that identied 6,727, 3,954 and 3,733 genes that were differentially expressed in salt-stressed calluses, leaves and roots of seedlings, respectively (Supplementary Data 46 and Supplementary Fig. S19). These differentially expressed genes (DEGs), which included those comprising expanded gene families and also those bearing the signature of positive selection (Fig. 2c, Supplementary Table S28, and Supplementary Figs S20 and S21), were similarly enriched in functional categories, such as oxidoreductase activity, transcription factor activity and ion transport (Fig. 3a). Several expanded gene families in the P. euphratica genome comprised transcription factors (Supplementary Table S24), for example, Myb, ERF, bZIP and WRKY, having a role in the regulation of gene expression in response to abiotic stress33. Some of these were extensively upregulated in response to salt stress (Fig. 2c; Supplementary Tables S28 and S29). Furthermore, the key genes regulating Na /H antiporters and controlling ion homeostasis34 (Supplementary Table S30), for example, NHX1 (Na /H exchanger 1), SOS2 (salt overly sensitive 2) and SOS3,

and those involved in the biosynthesis of ABA35,36, for example, BCH1 (b-carotene hydroxylase 1) and ZEP (zeaxanthin epoxidase), were upregulated in salt-stressed samples, which is consistent with previous research37.

We further compared the expression proles of theP. euphratica calluses in response to salt stress with those of theP. tomentosa (a salt-sensitive poplar21) calluses (Supplementary Table S31, Supplementary Figs S22 and S23). The results showed that many of the DEGs (2,278) were specic to P. euphratica (Fig 3b), and that more genes involved in cation transporter, oxidoreductase activity and response to abiotic stimulus were induced in P. euphratica than in P. tomentosa (Fig. 3c). Clustering analysis suggested that many of the DEGs exhibited different regulatory patterns in response to salinity between these two species (Fig. 3d). For example, the K uptake transporter KUP3 was extensively upregulated after 24 h of salt stress in

P. euphratica, but was maintained at control levels in P. tomentosa (Fig. 3e), indicating a critical role of this gene in controlling K homeostasis in P. euphratica. Transcript levels of this gene are

4 NATURE COMMUNICATIONS | 4:2797 | DOI: 10.1038/ncomms3797 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2013 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms3797 ARTICLE

Time (h)

0 3 (log10 P )

L R 6 12 24 48

P. euphratica

6 h 12 h 24 h 48 h

P. tomentosa

6 h 12 h 24 h 48 h

Transmembrane transporter activity

Transcription factor activity

Response to stimulus

Response to oxidative stress

Regulation of hormone levels

Photosystem

Phenylpropanoid biosynthesis

Oxidoreductase activity

Ion transport

Cellulose synthase

Cell wall organization

Cell communication

Carbohydrate metabolic

Log2 fold change

P. euphratica P. tomentosa

2,278 667 839

2 1 0 1 2

KUP3

NCL

400

400

200

200

P. euphraticaP. tomentosa

0

0

40

% Differential expressed

genes

40

PeNhaD1

20

SOS5

30

42

261 132

97 12

22

20

10

20

36

56

156

21

0

0

10

40

4

PeNhaD2

SDIR1

Cation transporter

Oxidoreductase

activity

Response to abiotic stimulus

Plastid

organization

0

2 3

Protein kinase

activity

20

2

0

0 0 6 12 24 48 (h) 0 6 12 24 48 (h)

Figure 3 | Comparative transcriptomics of P. euphratica and P. tomentosa under salt stress. (a) Functional category enrichment (P-value r0.05 by the Fishers exact test) of DEGs in P. euphratica leaves (L), roots (R) and time-course proles. (b) Venn diagram of the number of DEGs in P. euphratica and P.

tomentosa under salt stress (Supplementary Methods). (c) DEGs proportions in P. euphratica and P. tomentosa. The number of DEGs of each class is indicated above each bar. (d) Expression of the DEGs identied in P. euphratica and/or P. tomentosa. The heatmap was generated from hierarchical cluster analysis of genes. (e) Transcript levels of the genes showing different expression patterns between P. euphratica (blue) and P. tomentosa (orange). The transcript levels were determined by fragments per kilobase of exon per million fragments mapped (FPKM).

strongly induced by K starvation in Arabidopsis38. Previous research suggested that PeNhaD1, encoding a NhaD-type Na /

H antiporter, has a role in mediating sodium tolerance inP. euphratica39. Consistent with this, we identied two gene members encoding NhaD-type antiporters, both of which maintained transcript levels in P. euphratica, but which signicantly reduced transcript levels under salt stress after 12 h in P. tomentosa before regaining control levels after 24 h (Fig. 3e). Another transporter encoding gene, NCL (Na/Ca2 exchanger-like protein), involved in the maintenance of Ca2 homeostasis under salt stress in Arabidopsis40, was strongly upregulated in P. euphratica relative to its expression in salt-sensitive P. tomentosa. In addition, SOS5, a gene encoding a putative cell surface adhesion protein for the maintenance of cell wall integrity and architecture under salt stress in Arabidopsis41, was downregulated in P. tomentosa after 12 h of salt stress, in

contrast to maintenance of transcript levels recorded inP. euphratica. We further found that the gene SDIR1 (salt-and drought-induced ring nger 1), whose overexpression improves drought tolerance in transgenic rice42, was specically upregulated after 6 h and maintained high transcript levels until 12 h after salt stress in P. euphratica. Finally, transcription factors, such as ERF3 and NAC042, and the oxidoreductases AOX1D (alternative oxidase 1D) and NDA2 (alternative NAD(P)H dehydrogenase2), exhibited different expression patterns under salt stress inP. euphratica relative to those recorded in P. tomentosa (Supplementary Fig. S24).

DiscussionAbiotic stress factors, especially salinity and drought, restrict plant biomass production and pose an increasing threat to

NATURE COMMUNICATIONS | 4:2797 | DOI: 10.1038/ncomms3797 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

& 2013 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms3797

sustainable agriculture and forestry worldwide. Numerous studies have been conducted on the genetic and molecular mechanisms underlying salt tolerance in plants18,19,30,31, and have included genomic analyses of the extremophiles Thellungiella parvula16 and T. salsuginea43. However, our current understanding of these aspects of salt tolerance remains limited, especially for woody plants22,27. Populus euphratica is an excellent candidate for the analysis of salt tolerance22, as it displays apoplastic sodium accumulation and develops leaf succulence after prolonged salt exposure8. Consequently, in the last decade it has become a model for elucidating both physiological and molecular mechanisms of salt tolerance in tree species511. Using a newly developed fosmid-pooling strategy14, we sequenced and assembled the complex genome of P. euphratica with high heterozygosity and compared it with the closely related salt-sensitive model plant, P. trichocarpa.

We found that P. euphratica diverged from P. trichocarpa within the last 8 to 14 million years (Supplementary Fig. S15). Although both species shared at least two WGDs and exhibited extensive collinearity across the gene space (Fig. 1 and Supplementary Fig. S13), species-specic genes involved in stress tolerance, such as ion transport, ATPase activity, transcript factor activity and oxidoreductase activity, were selectively expanded and/or positively selected in the P. euphratica genome (Fig. 2 and Supplementary Tables S23S26). In this regard, the Na /K transporter HKT1 is of particular interest, because it is similarly expanded as tandem duplicated copies in both

T. parvula16 and T. salsuginea43. Further functional analysis of this gene family is needed to understand its critical role in salt tolerance in plants. In addition, other genes involved in ion transport and homeostasis, such as NhaD1, KUP3 and NCL, were distinctly upregulated under salt stress when compared with another salt-sensitive poplar, P. tomentosa. Our analyses taken together suggest that P. euphratica may have increased its salt tolerance through duplication and/or upregulation of multiple genes involved in ion transport and homeostasis. These ndings are important for an improved understanding of tree adaptation to salt stress and for accelerating the genetic improvement of cultivated poplars for growth on saline soils.

Methods

Genome sequencing and assembly. Genomic DNA was extracted from callus induced from P. euphratica shoots. Paired-end and mate-pair Illumina libraries were constructed with multiple insert sizes (13840 kb) according to the manufacturers instructions. In addition, 66,240 fosmid clones with 40 kb in length were randomly selected and two small insert (250 and 500 bp) libraries were constructed for each clone. All libraries were sequenced on the Illumina Genome Analyzer and HiSeq 2000 sequencing system. The raw reads were processed by removing low-quality reads, adapter sequences and possible contaminated reads. Then a hierarchical strategy14 was used for genome assembly (Supplementary Methods).

Gene prediction. We used the homology-based and de novo methods, as well as RNA-seq data, to predict genes in the P. euphratica genome. For homology-based gene prediction, protein sequences from ve plants (P. trichocarpa, Ricinus communis, Prunus persica, Cucumis sativus and Glycine max) were initially mapped onto the P. euphratica genome using TBLASTN and the homologous genome sequences were aligned against the matching proteins using GeneWise44 for accurate spliced alignments. Next, we used the de novo gene prediction methods Augustus45 and GenScan46 to predict protein-coding genes, using parameters trained for P. euphratica and A. thaliana. We then integrated the homologues and those from de novo approaches using GLEAN47 to produce a consensus gene set. In addition, we aligned all the RNA-seq reads to the reference genome by TopHat48, assembled the transcripts using Cufinks49 and predicted the open reading frames from the resultant data. Finally, we combined the GLEAN set with the gene models produced from RNA-seq to generate a more condent gene set.

Collinear block and genome duplication identication. Pairwise whole-genome alignment between P. euphratica and P. trichocarpa was constructed using the BLAST algorithm, and the scaffolds of P. euphratica were anchored to theP. trichocarpa corresponding chromosomes based on the consensus order of

matched regions. To detect the signature of genome duplication, the programme MCSCAN50 was used to dene a duplicated block. At least ve genes are required to call synteny. For each duplicated block, the 4DTv values were calculated and distributions were plotted.

Gene family clusters. The protein-coding genes from nine plant species(P. trichocarpa, Ricinus communis, Arabidopsis thaliana, Thellungiella parvula, Carica papaya, Fragaria vesca, Prunus persica, Vitis vinifera and Oryza sativa) were downloaded. The longest translation form was chosen to represent each gene, and stretches of genes encoding fewer than 50 amino acids were ltered out. The OrthoMCL51 method was then used to cluster all the genes into paralogous and orthologous groups. The 1,776 single-copy gene families obtained from this analysis were used to reconstruct phylogenies and estimate divergence time using MrBayes52 and the MCMCtree programme implemented in the Phylogenetic Analysis by Maximum Likelihood53. Calibration times were obtained from the TimeTree database (http://www.timetree.org/

Web End =http://www.timetree.org/).

Identication of PSGs. Using the orthologues identied by OrthoMCL as a raw data set, we rst masked sites with low quality (phred-like quality scoreo20) and that were detected as single nucleotide variants in P. euphratica coding sequences. We then aligned them using the codon option in the Probabilistic Alignment Kit54 programme for the detection of positive selection. Alignments shorter than 150 bp after removing sites with ambiguous data were discarded. Finally, we obtained 6,545 high-condence orthologues within the two poplar species and at least three of the other eight species (R. communis, A. thaliana, T. parvula, C. papaya, F. vesca,P. persica, V. vinifera and O. sativa), averaging B7.7 species per gene. These alignments together with an unrooted phylogenetic tree (constructed as described above) were used for subsequent molecular evolutionary analysis. For the estimation of the lineage-specic evolutionary rate, the values of Ka, Ks and Ka/Ks were calculated for 10,000 concatenated alignments constructed from 150 randomly chosen genes using the Codeml programme with the free-ratio model in the Phylogenetic Analysis by Maximum Likelihood53 package. To detect PSGs in either P. euphratica or P. trichocarpa lineage, the lineage was specied in turn as the foreground branch. We then used the optimized branch-site model55 in which likelihood ratio test P-values were computed, assuming that the null distribution was a 50:50 mixture of a w2-distribution with one degree of freedom and a point mass at zero. To minimize the false discovery rate, we manually ltered all PSGs with potential errors in their alignments.

Transcriptome sequencing and analysis. Total RNAs were extracted and strand-specic RNA-seq libraries were generated from samples using a cetyl trimethylammonium bromide procedure56 for transcriptome sequencing. The analysis was conducted on pooled samples of roots, leaves, ower buds, owers, xylem and phloem from two mature male P. euphratica trees and one mature female tree from the Talim Basin desert, Xinjiang, on control samples and also on salt-stressed samples (200 mM NaCl for 6, 12, 24 and 48 h) generated from the same calluses used in genome sequencing. RNA-seq libraries were sequenced on an Illumina Genome Analyzer platform. In addition, salt-stressed leaves and roots were collected and RNA samples were isolated for Illumina short-read sequencing. Three independent biological replicate samples were examined. The resulting reads were aligned to the P. euphratica genome sequences using TopHat48. After alignment, the count of mapped reads from each sample was derived and normalized to fragments per kilobase of exon per million fragments mapped for each predicted transcript using the Cufinks49 package. DEGs were identied using the programme Cuffdiff in the Cufinks package. We had tried to induce the calluses from P. trichocarpa, but they grew too slowly. We therefore sequenced transcriptomes of P. tomentosa calluses that had been subjected to salt-stress treatment (200 mM NaCl for 0, 6, 12, 24 and 48 h) and identied DEGs in this salt-sensitive species. In addition to the analysis as described for P. euphratica, we further assembled and annotated all reads from P. tomentosa using Trinity57 package (Supplementary Methods). We used the InParanoid58 software to identify 1:1 orthologues between P. euphratica and P. tomentosa, and aligned the coding sequences of the orthologues using Threaded Blockset Aligner59 to extract perfectly aligned consensus blocks. Finally, we counted the reads aligned to the consensus blocks for each sample and performed edgeR60 in R package to identify DEGs.

References

1. Food and Agricultural Organization of the United Nations. State of the Worlds Forests (FAO, 2011).

2. Oh, D. H., Dassanayake, M., Bohnert, H. J. & Cheeseman, J. M. Life at the extreme: lessons from the genome. Genome Biol. 13, 241 (2012).

3. Jansson, S. & Douglas, C. J. Populus: a model system for plant biology. Annu. Rev. Plant Biol. 58, 435458 (2007).

4. Jansson, S., Bhalarao, R. P. & Groover, A. T. Genetics and Genomics of Populus (Springer, 2010).

5. Gries, D. et al. Growth and water relations of Tamarix ramosissima and Populus euphratica on Taklamakan desert dunes in relation to depth to a permanent water table. Plant Cell Environ. 26, 725736 (2003).

6 NATURE COMMUNICATIONS | 4:2797 | DOI: 10.1038/ncomms3797 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2013 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms3797 ARTICLE

6. Brinker, M. et al. Linking the salt transcriptome with physiological responses of a salt-resistant Populus species as a strategy to identify genes important for stress acclimation. Plant Physiol. 154, 16971709 (2010).

7. Brosche, M. et al. Gene expression and metabolite proling of Populus euphratica growing in the Negev desert. Genome Biol. 6, R101 (2005).

8. Ottow, E. A. et al. Populus euphratica displays apoplastic sodium accumulation, osmotic adjustment by decreases in calcium and soluble carbohydrates, and develops leaf succulence under salt stress. Plant Physiol. 139, 17621772 (2005).

9. Janz, D. et al. Salt stress induces the formation of a novel type of pressure wood in two Populus species. N. Phytol. 194, 129141 (2012).

10. Wang, R. et al. Leaf photosynthesis, uorescence response to salinity and the relevance to chloroplast salt compartmentation and anti-oxidative stress in two poplars. Trees 21, 581591 (2007).

11. Gu, R. et al. Transcript identication and proling during salt stress and recovery of Populus euphratica. Tree Physiol. 24, 265276 (2004).

12. Tuskan, G. A. et al. The genome of black cottonwood, Populus trichocarpa (Torr. & Gray). Science 313, 15961604 (2006).

13. Alkan, C., Sajjadian, S. & Eichler, E. E. Limitations of next-generation genome sequence assembly. Nat. Meth. 8, 6165 (2011).

14. Zhang, G. et al. The oyster genome reveals stress adaptation and complexity of shell formation. Nature 490, 4954 (2012).

15. Jiao, Y. et al. Ancestral polyploidy in seed plants and angiosperms. Nature 473, 97100 (2011).

16. Dassanayake, M. et al. The genome of the extremophile crucifer Thellungiella parvula. Nat. Genet. 43, 913918 (2011).

17. Dassanayake, M., Oh, D. H., Hong, H., Bohnert, H. J. & Cheeseman, J. M. Transcription strength and halophytic lifestyle. Trends Plant Sci. 16, 13 (2011).

18. Ren, Z. H. et al. A rice quantitative trait locus for salt tolerance encodes a sodium transporter. Nat. Genet. 37, 11411146 (2005).

19. Davenport, R. J. et al. The Na transporter AtHKT1;1 controls retrieval of Na from the xylem in Arabidopsis. Plant Cell Environ. 30, 497507 (2007).

20. Chen, S., Li, J., Fritz, E., Wang, S. & Httermann, A. Sodium and chloride distribution in roots and transport in three poplar genotypes under increasing NaCl stress. For. Ecol. Manage 168, 217230 (2002).

21. Chen, S. et al. Effects of NaCl on shoot growth, transpiration, ion compartmentation, and transport in regenerated plants of Populus euphratica and Populus tomentosa. Can. J. For. Res. 33, 967975 (2003).

22. Chen, S. & Polle, A. Salinity tolerance of Populus. Plant Biol. 12, 317333 (2010).

23. Yang, Y. et al. A novel method to quantify H -ATPase-dependent Na transport across plasma membrane vesicles. Biochim. Biophys. Acta 1768, 20782088 2007.

24. Noctor, G. & Foyer, C. H. Ascorbate and glutathione: keeping active oxygen under control. Annu. Rev. Plant Physiol. Plant Mol. Biol. 49, 249279 (1998).

25. Hirayama, T. & Shinozaki, K. Perception and transduction of abscisic acid signals: keys to the function of the versatile plant hormone ABA. Trends Plant Sci. 12, 343351 (2007).

26. Wang, W., Vinocur, B., Shoseyov, O. & Altman, A. Role of plant heat-shock proteins and molecular chaperones in the abiotic stress response. Trends Plant Sci. 9, 244252 (2004).

27. Bartels, D. & Sunkar, R. Drought and salt tolerance in plants. Crit. Rev. Plant Sci. 24, 2358 (2005).

28. Taji, T. et al. Important roles of drought- and cold-inducible genes for galactinol synthase in stress tolerance in Arabidopsis thaliana. Plant J. 29, 417426 (2002).

29. Qiu, Q. et al. The yak genome and adaptation to life at high altitude. Nat. Genet. 44, 946949 (2012).

30. Zhu, J. et al. An enhancer mutant of Arabidopsis salt overly sensitive 3 mediates both ion homeostasis and the oxidative stress response. Mol. Cell Biol. 27, 52145224 (2007).

31. DAngelo, C. et al. Alternative complex formation of the Ca-regulated protein kinase CIPK1 controls abscisic acid-dependent and independent stress responses in Arabidopsis. Plant J. 48, 857872 (2006).

32. Larsson, K. E., Nystrom, B. & Liljenberg, C. A phosphatidylserine decarboxylase activity in root cells of oat (Avena sativa) is involved in altering membrane phospholipid composition during drought stress acclimation. Plant Physiol. 44, 211219 (2006).

33. Yamaguchi-Shinozaki, K. & Shinozaki, K. Organization of cis-acting regulatory elements in osmotic- and cold-stress-responsive promoters. Trends Plant Sci. 10, 8894 (2005).

34. Zhu, J. K. Plant salt tolerance. Trends Plant Sci. 6, 6671 (2001).35. Seo, M. & Koshiba, T. Complex regulation of ABA biosynthesis in plants. Trends Plant Sci. 7, 4148 (2002).

36. Du, H. et al. Characterization of the beta-carotene hydroxylase gene DSM2 conferring drought and oxidative stress resistance by increasing xanthophylls and abscisic acid synthesis in rice. Plant Physiol. 154, 13041318 (2010).

37. Qiu, Q. et al. Genome-scale transcriptome analysis of the desert poplar, Populus euphratica. Tree Physiol. 31, 452461 (2011).

38. Kim, E. J., Kwak, J. M., Uozumi, N. & Schroeder, J. I. AtKUP1: an Arabidopsis gene encoding high-afnity potassium transport activity. Plant Cell 10, 5162 (1998).

39. Ottow, E. A. et al. Molecular characterization of PeNhaD1: the rst member of the NhaD Na /H antiporter family of plant origin. Plant Mol. Biol. 58, 7588 (2005).

40. Wang, P. et al. A Na /Ca2 exchanger-like protein (AtNCL) involved in salt stress in Arabidopsis. J. Biol. Chem. 287, 4406244070 (2012).

41. Shi, H., Kim, Y., Guo, Y., Stevenson, B. & Zhu, J. K. The Arabidopsis SOS5 locus encodes a putative cell surface adhesion protein and is required for normal cell expansion. Plant Cell 15, 1932 (2003).

42. Gao, T. et al. OsSDIR1 overexpression greatly improves drought tolerance in transgenic rice. Plant Mol. Biol. 76, 145156 (2011).

43. Wu, H. J. et al. Insights into salt tolerance from the genome of Thellungiella salsuginea. Proc. Natl Acad. Sci. USA 109, 1221912224 (2012).

44. Birney, E., Clamp, M. & Durbin, R. GeneWise and Genomewise. Genome Res. 14, 988995 (2004).

45. Stanke, M. et al. AUGUSTUS: ab initio prediction of alternative transcripts. Nucleic Acids Res. 34, W435W439 (2006).

46. Salamov, A. A. & Solovyev, V. V. Ab initio gene nding in Drosophila genomic DNA. Genome Res. 10, 516522 (2000).

47. Elsik, C. G. et al. Creating a honey bee consensus gene set. Genome Biol. 8, R13 (2007).

48. Trapnell, C., Pachter, L. & Salzberg, S. L. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics 25, 11051111 (2009).

49. Trapnell, C. et al. Transcript assembly and quantication by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 28, 511515 (2010).

50. Tang, H. et al. Unraveling ancient hexaploidy through multiply-aligned angiosperm gene maps. Genome Res. 18, 19441954 (2008).

51. Li, L., Stoeckert, Jr. C. J. & Roos, D. S. OrthoMCL: identication of ortholog groups for eukaryotic genomes. Genome Res. 13, 21782189 (2003).

52. Huelsenbeck, J. P. & Ronquist, F. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics 17, 754755 (2001).

53. Yang, Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol. Biol. Evol. 24, 15861591 (2007).

54. Loytynoja, A. & Goldman, N. An algorithm for progressive multiple alignment of sequences with insertions. Proc. Natl Acad. Sci. USA 102, 1055710562 (2005).

55. Zhang, J., Nielsen, R. & Yang, Z. Evaluation of an improved branch-site likelihood method for detecting positive selection at the molecular Level. Mol. Biol. Evol. 22, 24722479 (2005).

56. Chang, S., Puryear, J. & Cairney, J. A simple and efcient method for isolating RNA from pine trees. Plant Mol. Biol. Rep. 11, 113116 (1993).

57. Grabherr, M. G. et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 29, 644652 (2011).

58. Berglund, A. C., Sjolund, E., Ostlund, G. & Sonnhammer, E. L. InParanoid 6: eukaryotic ortholog clusters with inparalogs. Nucleic Acids Res. 36, D263D266 (2008).

59. Blanchette, M. et al. Aligning multiple genomic sequences with the threaded blockset aligner. Genome Res. 14, 708715 (2004).

60. Robinson, M. D., McCarthy, D. J. & Smyth, G. K. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139140 (2010).

Acknowledgements

Financial support was provided by the National Key Project for Basic Research (2012CB114504), the National High Technology Research and Development Program of China (863 Program, No. 2013AA100605), the National Science and Technology Support Program (2013BAD22B01), the Fundamental Research Funds for the Central Universities (lzujbky-2009-k05), the International Collaboration 111 Projects of China, the 985 and 211 Projects of Lanzhou University and the Shenzhen Municipal Government (ZYC200903240077A).

Author contributions

Jianquan Liu designed and managed the project. Jun Wang led the genome sequencing. Juan Wang and Q.W. prepared the P. euphratica nucleic acid samples. Junyi Wang, Z.Y., G.Z., D.J., S.P., D.F., Yadan Luo, X.L. and H.Y. performed the DNA sequencing. Z.Y., Y.C., Z.W., C.G., L.Y., D.Z. and Jinchao Liu performed the genome assembly. Junyi Wang, Y.C., Z.W., X.W., G.W., Z.C., Yao Lu and Ye Yin performed the genome annotation. T.M., Q.H. and Q.Q. designed evolutionary analyses. T.M., G.Z., Y.C., J.Z., Q.Q., K.W., D.J., R.Z. and W.L. performed evolutionary analyses. Q.H., T.M., J.Z., K.W. and W.L. performed the synteny analyses. K.W., T.M., B.L. and H.M. performed the transcriptome analyses. Q.H., B.L. and Yongzhi Yang carried out data submission and

NATURE COMMUNICATIONS | 4:2797 | DOI: 10.1038/ncomms3797 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

& 2013 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms3797

database construction. T.M. and Jianquan Liu wrote the paper. G.A.T., S.P.D., R.J.A., Jun Wang, M.L., T.Y., Jia Li, D.W. and Y.W. revised the paper.

Additional information

Accession codes: The whole-genome shotgun project has been deposited in DDBJ/EMBL/ GenBank nucleotide core database under the accession code AOFL00000000. The version described in this paper is the rst version, AOFL01000000. All short-read data have been deposited in the Sequence Read Archive (SRA) under accession SRA061340. Raw sequence data of the transcriptomes have been deposited in the SRA under accession codes SRP028829 and SRP028830.

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Ma, T. et al. Genomic insights into salt adaptation in a desert poplar. Nat. Commun. 4:2797 doi: 10.1038/ncomms3797 (2013).

8 NATURE COMMUNICATIONS | 4:2797 | DOI: 10.1038/ncomms3797 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2013 Macmillan Publishers Limited. All rights reserved.