A new species of Populus and the extensive hybrid speciation arising from it on the Qinghai-Tibet Plateau.

While the diversity of species formation is broadly acknowledged, significant debate exists regarding the universal nature of hybrid species formation. Through an 18-year comprehensive study of all Populus species on the Qinghai-Tibet Plateau, 23 previously recorded species and 8 new species were identified. Based on morphological characteristics, these can be classified into three groups: species in section Leucoides, species with large leaves, and species with small leaves in section Tacamahaca. By conducting whole-genome re-sequencing of 150 genotypes from these 31 species, 2.28 million single nucleotide polymorphisms (SNPs) were identified. Phylogenetic analysis utilizing these SNPs not only revealed a highly intricate evolutionary network within the large-leaf species of section Tacamahaca but also confirmed that a new species, P. curviserrata, naturally hybridized with P. cathayana, P. szechuanica, and P. ciliata, resulting in 11 hybrid species. These findings indicate the widespread occurrence of hybrid species formation within this genus, with hybridization serving as a key evolutionary mechanism for Populus on the plateau. A novel hypothesis, "Hybrid Species Exterminating Their Ancestral Species (HSEAS)," is introduced to explain the mechanisms of hybrid species formation at three different scales: the entire plateau, the southeastern mountain region, and individual river valleys.

Drivers of genomic landscapes of differentiation across a Populus divergence gradient.

Speciation, the continuous process by which new species form, is often investigated by looking at the variation of nucleotide diversity and differentiation across the genome (hereafter genomic landscapes). A key challenge lies in how to determine the main evolutionary forces at play shaping these patterns. One promising strategy, albeit little used to date, is to comparatively investigate these genomic landscapes as progression through time by using a series of species pairs along a divergence gradient. Here, we resequenced 201 whole-genomes from eight closely related Populus species, with pairs of species at different stages along the divergence gradient to learn more about speciation processes. Using population structure and ancestry analyses, we document extensive introgression between some species pairs, especially those with parapatric distributions. We further investigate genomic landscapes, focusing on within-species (i.e. nucleotide diversity and recombination rate) and among-species (i.e. relative and absolute divergence) summary statistics of diversity and divergence. We observe relatively conserved patterns of genomic divergence across species pairs. Independent of the stage across the divergence gradient, we find support for signatures of linked selection (i.e. the interaction between natural selection and genetic linkage) in shaping these genomic landscapes, along with gene flow and standing genetic variation. We highlight the importance of investigating genomic patterns on multiple species across a divergence gradient and discuss prospects to better understand the evolutionary forces shaping the genomic landscapes of diversity and differentiation.

Combined transcriptome and proteome analysis provides insights into anthocyanin accumulation in the leaves of red-leaved poplars.

KEY MESSAGE: Anthocyanin was highly accumulated in the leaves of red-leaved poplars; Many structural genes involved in anthocyanin synthesis were significantly up-regulated in 'Quanhong' and 'Xuanhong'; TTG2, HYH, and HY5 may be directly involved in the regulation of anthocyanin synthesis in both red-leaved poplars. The red-leaved poplar cultivars 'Quanhong' and 'Xuanhong' are bud mutations of Populus deltoides cv. 'Zhonglin 2025'. These cultivars are valued for their beautiful shape, lack of flying catkins, and ornamental leaf colors. However, the understanding of the molecular mechanism of anthocyanin accumulation in the leaves of red-leaved poplars is still unclear. Here, we profiled the changes of pigment content, transcriptome and proteome expression in the leaves of three poplar cultivars and the results showed that the ratios of anthocyanin to total chlorophyll in both red-leaved poplars were higher than that in 'Zhonglin 2025', indicating that the anthocyanin was highly accumulated in the leaves of red-leaved poplars. Based on the results of combined transcriptome and proteome analysis, 15 and 11 differentially expressed genes/proteins involved in anthocyanin synthesis were screened in 'Quanhong' and 'Xuanhong', respectively, indicating that the two red-leaved poplar cultivars have slightly different patterns of regulating anthocyanin biosynthesis. Among the 120 transcription factors, 3 (HY5, HYH, and TTG2), may be directly involved in the regulation of anthocyanin synthesis in both red-leaved poplars. This study screens the candidate genes involved in anthocyanin accumulation in the leaves of red-leaved poplars and lays a foundation for further exploring the molecular mechanism of leaf red coloration in red-leaved poplars.

Molecular Karyotyping on Populus simonii × P. nigra and the Derived Doubled Haploid.

The molecular karyotype could represent the basic genetic make-up in a cell nucleus of an organism or species. A doubled haploid (DH) is a genotype formed from the chromosome doubling of haploid cells. In the present study, molecular karyotype analysis of the poplar hybrid Populus simonii × P. nigra (P. xiaohei) and the derived doubled haploids was carried out with labeled telomeres, rDNA, and two newly repetitive sequences as probes by fluorescence in situ hybridization (FISH). The tandem repeats, pPC349_XHY and pPD284_XHY, with high-sequence homology were used, and the results showed that they presented the colocalized distribution signal in chromosomes. For P. xiaohei, pPD284_XHY produced hybridizations in chromosomes 1, 5, 8, and 9 in the hybrid. The combination of pPD284_XHY, 45S rDNA, and 5S rDNA distinctly distinguished six pairs of chromosomes, and the three pairs of chromosomes showed a significant difference in the hybridization between homologous chromosomes. The repeat probes used produced similar FISH hybridizations in the DH; nevertheless, pPD284_XHY generated an additional hybridization site in the telomere region of chromosome 14. Moreover, two pairs of chromosomes showed differential hybridization distributions between homologous chromosomes. Comparisons of the distinguished chromosomes between hybrid and DH poplar showed that three pairs of chromosomes in the DH presented hybridization patterns that varied from those of the hybrid. The No. 8 chromosome in DH and one of the homologous chromosomes in P. xiaohei shared highly similar FISH patterns, which suggested the possibility of intact or mostly partial transfer of the chromosome between the hybrid and DH. Our study will contribute to understanding the genetic mechanism of chromosomal variation in P. xiaohei and derived DH plants.

Genome-Wide Comprehensive Analysis of the GASA Gene Family in Populus.

Gibberellic acid-stimulated Arabidopsis (GASA) proteins, as cysteine-rich peptides (CRPs), play roles in development and reproduction and biotic and abiotic stresses. Although the GASA gene family has been identified in plants, the knowledge about GASAs in Populus euphratica, the woody model plant for studying abiotic stress, remains limited. Here, we referenced the well-sequenced Populus trichocarpa genome, and identified the GASAs in the whole genome of P. euphratica and P. trichocarpa. 21 candidate genes in P. trichocarpa and 19 candidate genes in P. euphratica were identified and categorized into three subfamilies by phylogenetic analysis. Most GASAs with signal peptides were located extracellularly. The GASA genes in Populus have experienced multiple gene duplication events, especially in the subfamily A. The evolution of the subfamily A, with the largest number of members, can be attributed to whole-genome duplication (WGD) and tandem duplication (TD). Collinearity analysis showed that WGD genes played a leading role in the evolution of GASA genes subfamily B. The expression patterns of P. trichocarpa and P. euphratica were investigated using the PlantGenIE database and the real-time quantitative PCR (qRT-PCR), respectively. GASA genes in P. trichocarpa and P. euphratica were mainly expressed in young tissues and organs, and almost rarely expressed in mature leaves. GASA genes in P. euphratica leaves were also widely involved in hormone responses and drought stress responses. GUS activity assay showed that PeuGASA15 was widely present in various organs of the plant, especially in vascular bundles, and was induced by auxin and inhibited by mannitol dramatically. In summary, this present study provides a theoretical foundation for further research on the function of GASA genes in P. euphratica.

Chloroplast (Cp) Transcriptome of P. davidiana Dode×P. bolleana Lauch provides insight into the Cp drought response and Populus Cp phylogeny.

BACKGROUND: Raw second-generation (2G) lignocellulosic biomass materials have the potential for development into a sustainable and renewable source of energy. Poplar is regarded as a promising 2G material (P. davidiana Dode×P. bolleana Lauch, P. bolleana, P. davidiana, P. euphratica, et al). However, their large-scale commercialization still faces many obstacles. For example, drought prevents sufficient irrigation or rainfall, which can reduce soil moisture and eventually destroy the chloroplast, the plant photosynthetic organelle. Heterosis is widely used in the production of drought-tolerant materials, such as the superior clone "Shanxinyang" selected from the offspring of Populus davidiana Dode×Populus bolleana Lauch. Because it produces good wood and is easily genetically transformed, "Shanxinyang" has become a promising material for use in tree genetics. It is also one of the most abundant biofuel plants in northern China. Understanding the genetic features of chloroplasts, the cp transcriptome and physiology is crucial to elucidating the chloroplast drought-response model. RESULTS: In this study, the whole genome of "Shanxinyang" was sequenced. The chloroplast genome was assembled, and chloroplast structure was analysed and compared with that of other popular plants. Chloroplast transcriptome analysis was performed under drought conditions. The total length of the "Shanxinyang" chloroplast genome was 156,190 bp, the GC content was 36.75%, and the genome was composed of four typical areas (LSC, IRa, IRb, and SSC). A total of 114 simple repeats were detected in the chloroplast genome of "Shanxinyang". In cp transcriptome analysis, we found 161 up-regulated and 157 down-regulated genes under drought, and 9 cpDEGs was randomly selected to conduct reverse transcription (RT)-qPCR., in which the Log2 (fold change) was significantly consistent with the qPCR results. The analysis of chloroplast transcription under drought provided clues for understanding chloroplast function under drought. The phylogenetic position of "Shanxinyang" within Populus was analysed by using the chloroplast genome sequences of 23 Populus plants, showing that "Shanxinyang" belongs to Sect. Populus and is sister to Populus davidiana. Further, mVISTA analysis showed that the variation in non-coding (regulatory) regions was greater than that in coding regions, which suggests that further attention should be paid to the chloroplast in order to obtain new evolutionary or functional insights related to aspects of plant biology. CONCLUSIONS: Our findings indicate that complex prokaryotic genome regulation occurs when processing transcripts under drought stress. The results not only offer clues for understanding the chloroplast genome and transcription features in woody plants but also serve as a basis for future molecular studies on poplar species.

Evidence for widespread selection in shaping the genomic landscape during speciation of Populus.

Increasing our understanding of how evolutionary processes drive the genomic landscape of variation is fundamental to a better understanding of the genomic consequences of speciation. However, genome-wide patterns of within- and between- species variation have not been fully investigated in most forest tree species despite their global ecological and economic importance. Here, we use whole-genome resequencing data from four Populus species spanning the speciation continuum to reconstruct their demographic histories and investigate patterns of diversity and divergence within and between species. Using Populus trichocarpa as an outgroup species, we further infer the genealogical relationships and estimate the extent of ancient introgression among the three aspen species (Populus tremula, Populus davidiana and Populus tremuloides) throughout the genome. Our results show substantial variation in these patterns along the genomes with this variation being strongly predicted by local recombination rates and the density of functional elements. This implies that the interaction between recurrent selection and intrinsic genomic features has dramatically sculpted the genomic landscape over long periods of time. In addition, our findings provide evidence that, apart from background selection, recent positive selection and long-term balancing selection have also been crucial components in shaping patterns of genome-wide variation during the speciation process.

Element content and expression of genes of interest in guard cells are connected to spatiotemporal variations in stomatal conductance.

Element content and expression of genes of interest on single cell types, such as stomata, provide valuable insights into their specific physiology, improving our understanding of leaf gas exchange regulation. We investigated how far differences in stomatal conductance (gs ) can be ascribed to changes in guard cells functioning in amphistomateous leaves. gs was measured during the day on both leaf sides, on well-watered and drought-stressed trees (two Populus euramericana Moench and two Populus nigra L. genotypes). In parallel, guard cells were dissected for element content and gene expressions analyses. Both were strongly arranged according to genotype, and drought had the lowest impact overall. Normalizing the data by genotype highlighted a structure on the basis of leaf sides and time of day both for element content and gene expression. Guard cells magnesium, phosphorus, and chlorine were the most abundant on the abaxial side in the morning, where gs was at the highest. In contrast, genes encoding H+ -ATPase and aquaporins were usually more abundant in the afternoon, whereas genes encoding Ca2+ -vacuolar antiporters, K+ channels, and ABA-related genes were in general more abundant on the adaxial side. Our work highlights the unique physiology of each leaf side and their analogous rhythmicity through the day.

Genetic diversity and population structure of black cottonwood (Populus deltoides) revealed using simple sequence repeat markers.

BACKGROUND: Black cottonwood (Populus deltoides) is one of the keystone forest tree species, and has become the main breeding parents in poplar hybrid breeding. However, the genetic diversity and population structure of the introduced resources are not fully understood. RESULTS: In the present study, five loci containing null alleles were excluded and 15 pairs of SSR (simple sequence repeat) primers were used to analyze the genetic diversity and population structure of 384 individuals from six provenances (Missouri, Iowa, Washington, Louisiana, and Tennessee (USA), and Quebec in Canada) of P. deltoides. Ultimately, 108 alleles (Na) were detected; the expected heterozygosity (He) per locus ranged from 0.070 to 0.905, and the average polymorphic information content (PIC) was 0.535. The provenance 'Was' had a relatively low genetic diversity, while 'Que', 'Lou', and 'Ten' provenances had high genetic diversity, with Shannon's information index (I) above 1.0. The mean coefficient of genetic differentiation (Fst) and gene flow (Nm) were 0.129 and 1.931, respectively. Analysis of molecular variance (AMOVA) showed that 84.88% of the genetic variation originated from individuals. Based on principal coordinate analysis (PCoA) and STRUCTURE cluster analysis, individuals distributed in the Mississippi River Basin were roughly classified as one group, while those distributed in the St. Lawrence River Basin and Columbia River Basin were classified as another group. The cluster analysis based on the population level showed that provenance 'Iow' had a small gene flow and high degree of genetic differentiation compared with the other provenances, and was classified into one group. There was a significant relationship between genetic distance and geographical distance. CONCLUSIONS: P. deltoides resources have high genetic diversity and there is a moderate level of genetic differentiation among provenances. Geographical isolation and natural conditions may be the main factors causing genetic differences among individuals. Individuals reflecting population genetic information can be selected to build a core germplasm bank. Meanwhile, the results could provide theoretical support for the scientific management and efficient utilization of P. deltoides genetic resources, and promote the development of molecular marker-assisted breeding of poplar.

Comparative Proteomic Analysis of Plant-Pathogen Interactions in Resistant and Susceptible Poplar Ecotypes Infected with Botryosphaeria dothidea.

Poplar are important forestry species in China, but the Botryosphaeria dothidea pathogen causes serious economic losses worldwide. To identify candidate B. dothidea resistance proteins and explore the molecular mechanisms involved in poplar-pathogen interactions, proteomic responses of stem samples from resistant and susceptible poplar ecotypes to B. dothidea were investigated using nanoflow liquid chromatography-tandem mass spectrometry with label-free quantitative analysis. We identified 588 proteins, divided into 21 biological process categories including 48 oxidoreductases, 72 hydrolytic enzymes, 80 metabolic enzymes, and 29 proteins of unknown function. Differential proteome analysis revealed large differences between resistant Populus tomentosa Carr and susceptible Populus beijingensis Hsu ecotypes before and after inoculation. Among 102 identified proteins, 22 were highly upregulated in the resistant genotype but downregulated in the susceptible genotype. Proteins induced in P. tomentosa Carr in response to B. dothidea are associated with plant defenses including oxidoreductase activity (catalase, isocitrate dehydrogenase, and superoxide dismutase), phenylpropanoid biosynthesis and phenylalanine metabolism (alcohol dehydrogenase), photosynthesis (ATP synthase subunit alpha, ATP synthase gamma chain, photosystem I P700 chlorophyll a apoprotein A2, photosystem II CP47 chlorophyll apoprotein), carbon fixation (pyruvate kinase, triosephosphate isomerase, malic enzyme, phosphoglycerate kinase, ribulose-1,5-bisphosphate carboxylase, and ribulose bisphosphate carboxylase small chain), and glycolysis/gluconeogenesis (fructose-bisphosphate aldolase). Kyoto Encyclopedia of Genes and Genomes pathway analysis identified 168 proteins related to metabolic pathways, 41 proteins related to the biosynthesis of phenylpropanoids, and 36 proteins related to the biosynthesis of plant hormones, the biosynthesis of alkaloids derived from ornithine, lysine, and nicotinic acid, and photosynthesis in response to B. dothidea. Our findings provide insight into plant-pathogen interactions in resistant and susceptible poplar ecotypes infected with B. dothidea and could assist the development of novel strategies for fighting poplar canker disease.

PtrARF2.1 Is Involved in Regulation of Leaf Development and Lignin Biosynthesis in Poplar Trees.

Auxin response factors (ARFs) are important regulators modulating the expression of auxin-responsive genes in various biological processes in plants. In the Populus genome, a total of 39 ARF members have been identified, but their detailed functions are still unclear. In this study, six poplar auxin response factor 2 (PtrARF2) members were isolated from P. trichocarpa. Expression pattern analysis showed that PtrARF2.1 is highly expressed in leaf tissues compared with other PtrARF2 genes and significantly repressed by exogenous auxin treatment. PtrARF2.1 is a nuclear-localized protein without transcriptional activation activity. Knockdown of PtrARF2.1 by RNA interference (RNAi) in poplars led to the dwarf plant, altered leaf shape, and reduced size of the leaf blade, while overexpression of PtrARF2.1 resulted in a slight reduction in plant height and the similar leaf phenotype in contrast to the wildtype. Furthermore, histological staining analysis revealed an ectopic deposition of lignin in leaf veins and petioles of PtrARF2.1-RNAi lines. RNA-Seq analysis showed that 74 differential expression genes (DEGs) belonging to 12 transcription factor families, such as NAM, ATAF and CUC (NAC), v-myb avian myeloblastosis viral oncogene homolog (MYB), ethylene response factors (ERF) and basic helix-loop-helix (bHLH), were identified in PtrARF2.1-RNAi leaves and other 24 DEGs were associated with the lignin biosynthetic pathway. Altogether, the data indicate that PtrARF2.1 plays an important role in regulating leaf development and influences the lignin biosynthesis in poplars.

A New Propolis Type from Changbai Mountains in North-east China: Chemical Composition, Botanical Origin and Biological Activity.

Propolis is a bee product with a wide range of biological activities and its chemical compounds depend highly on the type of plant accessible to the bees. The Changbai Mountains are a major mountain range in Northeast China and are one of the major bee product-producing areas in China. In this study, we evaluated the total phenolic acids and flavonoid contents as well as the antioxidant activity of propolis sampled from the Changbai Mountains area (CBM). We identified the major compounds and qualified their contents by HPLC-ESI/MS and HPLC-UV, and found that the content of p-coumaric acid and an unknown peak (CBE) in CBM propolis was higher than in propolis from other parts of China. The unknown compound CBE was isolated, purified, and identified as benzyl p-coumarate by MS and NMR. Possible plant sources of CBM propolis are Populus davidiana dode and Populus simonii Carr, which widely distributed in the Changbai Mountains area. CBM propolis is a new propolis type, that could be an excellent raw material for health foods and pharmaceuticals.

Phylogeographic patterns of the desert poplar in Northwest China shaped by both geology and climatic oscillations.

BACKGROUND: The effects of historical geology and climatic events on the evolution of plants around the Qinghai-Tibetan Plateau region have been at the center of debate for years. To identify the influence of the uplift of the Tianshan Mountains and/or climatic oscillations on the evolution of plants in arid northwest China, we investigated the phylogeography of the Euphrates poplar (Populus euphratica) using chloroplast DNA (cpDNA) sequences and nuclear microsatellites, and estimated its historical distribution using Ecological Niche Modeling (ENM). RESULTS: We found that the Euphrates poplar differed from another desert poplar, P. pruinosa, in both nuclear and chloroplast DNA. The low clonal diversity in both populations reflected the low regeneration rate by seed/seedlings in many locations. Both cpDNA and nuclear markers demonstrated a clear divergence between the Euphrates poplar populations from northern and southern Xinjiang regions. The divergence time was estimated to be early Pleistocene based on cpDNA, and late Pleistocene using an Approximate Bayesian Computation analysis based on microsatellites. Estimated gene flow was low between these two regions, and the limited gene flow occurred mainly via dispersal from eastern regions. ENM analysis supported a wider distribution of the Euphrates poplar at 3 Ma, but a more constricted distribution during both the glacial period and the interglacial period. CONCLUSIONS: These results indicate that the deformation of the Tianshan Mountains has impeded gene flow of the Euphrates poplar populations from northern and southern Xinjiang, and the distribution constriction due to climatic oscillations further accelerated the divergence of populations from these regions. To protect the desert poplars, more effort is needed to encourage seed germination and seedling establishment, and to conserve endemic gene resources in the northern Xinjiang region.

Adaptive introgression and maintenance of a trispecies hybrid complex in range-edge populations of Populus.

In hybrid zones occurring in marginal environments, adaptive introgression from one species into the genomic background of another may constitute a mechanism facilitating adaptation at range limits. Although recent studies have improved our understanding of adaptive introgression in widely distributed tree species, little is known about the dynamics of this process in populations at the margins of species ranges. We investigated the extent of introgression between three species of the genus Populus sect. Tacamahaca (P. balsamifera, P. angustifolia and P. trichocarpa) at the margins of their distributions in the Rocky Mountain region of the United States and Canada. Using genotyping by sequencing (GBS), we analysed ~ 83,000 single nucleotide polymorphisms genotyped in 296 individuals from 29 allopatric and sympatric populations of the three species. We found a trispecies hybrid complex present throughout the zone of range overlap, including early as well as advanced generation backcross hybrids, indicating recurrent gene flow in this hybrid complex. Using genomic cline analysis, we found evidence of non-neutral patterns of introgression at 23% of loci in hybrids, of which 47% and 8% represented excess ancestry from P. angustifolia and P. balsamifera, respectively. Gene ontology analysis suggested these genomic regions were enriched for genes associated with photoperiodic regulation, metal ion transport, maintenance of redox homeostasis and cell wall metabolites involved in regulation of seasonal dormancy. Our study demonstrates the role of adaptive introgression in a multispecies hybrid complex in range-edge populations and has implications for understanding the evolutionary dynamics of adaptation in hybrid zones, especially at the margins of species distributions.

Adaptive evolution and segregating load contribute to the genomic landscape of divergence in two tree species connected by episodic gene flow.

Speciation often involves repeated episodes of genetic contact between divergent populations before reproductive isolation (RI) is complete. Whole-genome sequencing (WGS) holds great promise for unravelling the genomic bases of speciation. We have studied two ecologically divergent, hybridizing species of the 'model tree' genus Populus (poplars, aspens, cottonwoods), Populus alba and P. tremula, using >8.6 million single nucleotide polymorphisms (SNPs) from WGS of population pools. We used the genomic data to (i) scan these species' genomes for regions of elevated and reduced divergence, (ii) assess key aspects of their joint demographic history based on genomewide site frequency spectra (SFS) and (iii) infer the potential roles of adaptive and deleterious coding mutations in shaping the genomic landscape of divergence. We identified numerous small, unevenly distributed genome regions without fixed polymorphisms despite high overall genomic differentiation. The joint SFS was best explained by ancient and repeated gene flow and allowed pinpointing candidate interspecific migrant tracts. The direction of selection (DoS) differed between genes in putative migrant tracts and the remainder of the genome, thus indicating the potential roles of adaptive divergence and segregating deleterious mutations on the evolution and breakdown of RI. Genes affected by positive selection during divergence were enriched for several functionally interesting groups, including well-known candidate 'speciation genes' involved in plant innate immunity. Our results suggest that adaptive divergence affects RI in these hybridizing species mainly through intrinsic and demographic processes. Integrating genomic with molecular data holds great promise for revealing the effects of particular genetic pathways on speciation.

Analysis of cellulose synthase gene expression strategies in higher plants using RNA-sequencing data.

The transcriptomes from different organs and tissues of western poplar, eucalyptus, soybean and common bean were studied. The expression level of cellulose synthase genes was notably different in different types of tissues and organs in studied plants. For common bean and eucalyptus transcriptome the domination of certain cellulose synthase genes was typical. These prevailing genes made up more than 50 % of the total expression pull of cellulose synthases. On the contrary, cellulose synthase expression pulls of wes-tern poplar and soybean were distributed between multiple genes. The different expression strategies of CesA-genes may reflect a phylogenetic processes that occurred in genomes studied.

Genome-wide analysis of lectin receptor-like kinases in Populus.

BACKGROUND: Receptor-like kinases (RLKs) belong to a large protein family with over 600 members in Arabidopsis and over 1000 in rice. Among RLKs, the lectin receptor-like kinases (LecRLKs) possess a characteristic extracellular carbohydrate-binding lectin domain and play important roles in plant development and innate immunity. There are 75 and 173 LecRLKs in Arabidopsis and rice, respectively. However, little is known about LecRLKs in perennial woody plants. RESULTS: Here we report the genome-wide analysis of classification, domain architecture and expression of LecRLKs in the perennial woody model plant Populus. We found that the LecRLK family has expanded in Populus to a total of 231, including 180 G-type, 50 L-type and 1 C-type LecRLKs. Expansion of the Populus LecRLKs (PtLecRLKs) occurred partially through tandem duplication. Based on domain architecture and orientation features, we classified PtLecRLKs into eight different classes. RNA-seq-based transcriptomics analysis revealed diverse expression patterns of PtLecRLK genes among leaves, stems, roots, buds and reproductive tissues and organs. CONCLUSIONS: This study offers a comprehensive view of LecRLKs in the perennial woody model plant Populus and provides a foundation for functional characterization of this important family of receptor-like kinases.

Genetic origin and composition of a natural hybrid poplar Populus × jrtyschensis from two distantly related species.

BACKGROUND: The factors that contribute to and maintain hybrid zones between distinct species are highly variable, depending on hybrid origins, frequencies and fitness. In this study, we aimed to examine genetic origins, compositions and possible maintenance of Populus × jrtyschensis, an assumed natural hybrid between two distantly related species. This hybrid poplar occurs mainly on the floodplains along the river valleys between the overlapping distributions of the two putative parents. RESULTS: We collected 566 individuals from 45 typical populations of P. × jrtyschensis, P. nigra and P. laurifolia. We genotyped them based on the sequence variations of one maternally inherited chloroplast DNA (cpDNA) fragment and genetic polymorphisms at 20 SSR loci. We further sequenced eight nuclear genes for 168 individuals from 31 populations. Two groups of cpDNA haplotypes characteristic of P. nigra and P. laurifolia respectively were both recovered for P. × jrtyschensis. Genetic structures and coalescent tests of two sets of nuclear population genetic data suggested that P. × jrtyschensis originated from hybridizations between the two assumed parental species. All examined populations of P. × jrtyschensis comprise mainly F1 hybrids from interspecific hybridizations between P. nigra and P. laurifolia. In the habitats of P. × jrtyschensis, there are lower concentrations of soil nitrogen than in the habitats occupied by the other two species. CONCLUSIONS: Our extensive examination of the genetic composition of P. × jrtyschensis suggested that it is typical of F1-dominated hybrid zones. This finding plus the low concentration of soil nitrogen in the floodplain soils support the F1-dominated bounded hybrid superiority hypothesis of hybrid zone maintenance for this particular hybrid poplar.

Genomic and functional approaches reveal a case of adaptive introgression from Populus balsamifera (balsam poplar) in P. trichocarpa (black cottonwood).

Natural hybrid zones in forest trees provide systems to study the transfer of adaptive genetic variation by introgression. Previous landscape genomic studies in Populus trichocarpa, a keystone tree species, indicated genomic footprints of admixture with its sister species Populus balsamifera and identified candidate genes for local adaptation. Here, we explored the patterns of introgression and signals of local adaptation in P. trichocarpa and P. balsamifera, employing genome resequencing data from three chromosomes in pure species and admixed individuals from wild populations. Local ancestry analysis in admixed P. trichocarpa revealed a telomeric region in chromosome 15 with P. balsamifera ancestry, containing several candidate genes for local adaptation. Genomic analyses revealed signals of selection in certain genes in this region (e.g. PRR5, COMT1), and functional analyses based on gene expression variation and correlations with adaptive phenotypes suggest distinct functions of the introgressed alleles. In contrast, a block of genes in chromosome 12 paralogous to the introgressed region showed no signs of introgression or signatures of selection. We hypothesize that the introgressed region in chromosome 15 has introduced modular or cassette-like variation into P. trichocarpa. These linked adaptive mutations are associated with a block of genes in chromosome 15 that appear to have undergone neo- or subfunctionalization relative to paralogs in a duplicated region on chromosome 12 that show no signatures of adaptive variation. The association between P. balsamifera introgressed alleles with the expression of adaptive traits in P. trichocarpa supports the hypothesis that this is a case of adaptive introgression in an ecologically important foundation species.

Multilocus analysis of nucleotide variation and speciation in three closely related Populus (Salicaceae) species.

Historical tectonism and climate oscillations can isolate and contract the geographical distributions of many plant species, and they are even known to trigger species divergence and ultimately speciation. Here, we estimated the nucleotide variation and speciation in three closely related Populus species, Populus tremuloides, P. tremula and P. davidiana, distributed in North America and Eurasia. We analysed the sequence variation in six single-copy nuclear loci and three chloroplast (cpDNA) fragments in 497 individuals sampled from 33 populations of these three species across their geographic distributions. These three Populus species harboured relatively high levels of nucleotide diversity and showed high levels of nucleotide differentiation. Phylogenetic analysis revealed that P. tremuloides diverged earlier than the other two species. The cpDNA haplotype network result clearly illustrated the dispersal route from North America to eastern Asia and then into Europe. Molecular dating results confirmed that the divergence of these three species coincided with the sundering of the Bering land bridge in the late Miocene and a rapid uplift of the Qinghai-Tibetan Plateau around the Miocene/Pliocene boundary. Vicariance-driven successful allopatric speciation resulting from historical tectonism and climate oscillations most likely played roles in the formation of the disjunct distributions and divergence of these three Populus species.