Fitness dynamics within a poplar hybrid zone: II. Impact of exotic sex on native poplars in an urban jungle

Sex-determining region complements traditionally used in phylogenetic studies nuclear and chloroplast sequences in investigation of Aigeiros Duby and Tacamahaca Spach poplars (genus Populus L., Salicaceae)

Members of the genus Populus L. play an important role in the formation of forests in the northern hemisphere and are used in urban landscaping and timber production. Populus species of closely related sections show extensive hybridization. Therefore, the systematics of the genus is rather complicated, especially for poplars of hybrid origin. We aimed to assess the efficiency of application of the sex-determining region (SDR) in addition to the nuclear and chloroplast genome loci traditionally used in phylogenetic studies of poplars to investigate relationships in sections Aigeiros Duby and Tacamahaca Spach. Targeted deep sequencing of NTS 5S rDNA, ITS, DSH 2, DSH 5, DSH 8, DSH 12, DSH 29, 6, 15, 16, X18, trnG-psbK-psbI, rps2-rpoC2, rpoC2-rpoC1, as well as SDR and ARR17 gene was performed for 379 poplars. The SDR and ARR17 gene together with traditionally used multicopy and single-copy loci of nuclear and chloroplast DNA allowed us to obtain a clustering that is most consistent with poplar systematics based on morphological data and to shed light on several controversial hypotheses about the origin of the studied taxa (for example, the inexpediency of separating P. koreana, P. maximowiczii, and P. suaveolens into different species). We present a scheme of relationships between species and hybrids of sections Aigeiros and Tacamahaca based on molecular genetic, morphological, and geographical data. The geographical proximity of species and, therefore, the possibility of hybridization between them appear to be more important than the affiliation of species to the same section. We speculate that sections Aigeiros and Tacamahaca are distinguished primarily on an ecological principle (plain and mountain poplars) rather than on a genetic basis. Joint analysis of sequencing data for the SDR and chloroplast genome loci allowed us to determine the ancestors of P. × petrovskoe - P. laurifolia (female tree) × P. × canadensis (male tree), and P. × rasumovskoe - P. nigra (female tree) × P. suaveolens (male tree). Thus, the efficiency of using the SDR for the study of poplars of sections Aigeiros and Tacamahaca and the prospects of its use for the investigation of species of the genus Populus were shown.

Genome Assembly and Sex-Determining Region of Male and Female Populus × sibirica

The genus Populus is presented by dioecious species, and it became a promising object to study the genetics of sex in plants. In this work, genomes of male and female Populus × sibirica individuals were sequenced for the first time. To achieve high-quality genome assemblies, we used Oxford Nanopore Technologies and Illumina platforms. A protocol for the isolation of long and pure DNA from young poplar leaves was developed, which enabled us to obtain 31 Gb (N50 = 21 kb) for the male poplar and 23 Gb (N50 = 24 kb) for the female one using the MinION sequencer. Genome assembly was performed with different tools, and Canu provided the most complete and accurate assemblies with a length of 818 Mb (N50 = 1.5 Mb) for the male poplar and 816 Mb (N50 = 0.5 Mb) for the female one. After polishing with Racon and Medaka (Nanopore reads) and then with POLCA (Illumina reads), assembly completeness was 98.45% (87.48% duplicated) for the male and 98.20% (76.77% duplicated) for the female according to BUSCO (benchmarking universal single-copy orthologs). A high proportion of duplicated BUSCO and the increased genome size (about 300 Mb above the expected) pointed at the separation of haplotypes in a large part of male and female genomes of P. × sibirica. Due to this, we were able to identify two haplotypes of the sex-determining region (SDR) in both assemblies; and one of these four SDR haplotypes, in the male genome, contained partial repeats of the ARR17 gene (Y haplotype), while the rest three did not (X haplotypes). The analysis of the male P. × sibirica SDR suggested that the Y haplotype originated from P. nigra, while the X haplotype is close to P. trichocarpa and P. balsamifera species. Moreover, we revealed a Populus-specific repeat that could be involved in translocation of the ARR17 gene or its part to the SDR of P. × sibirica and other Populus species. The obtained results expand our knowledge on SDR features in the genus Populus and poplar phylogeny.

Genome and Transcriptome Sequencing of Populus × sibirica Identified Sex-Associated Allele-Specific Expression of the CLC Gene

Transcriptome sequencing of leaves, catkin axes, and flowers from male and female trees of Populus × sibirica and genome sequencing of the same plants were performed for the first time. The availability of both genome and transcriptome sequencing data enabled the identification of allele-specific expression. Such an analysis was performed for genes from the sex-determining region (SDR). P. × sibirica is an intersectional hybrid between species from sections Aigeiros (Populus nigra) and Tacamahaca (Populus laurifolia, Populus suaveolens, or Populus × moskoviensis); therefore, a significant number of heterozygous polymorphisms were identified in the SDR that allowed us to distinguish between alleles. In the SDR, both allelic variants of the TCP (T-complex protein 1 subunit gamma), CLC (Chloride channel protein CLC-c), and MET1 (DNA-methyltransferase 1) genes were expressed in females, while in males, two allelic variants were expressed for TCP and MET1 but only one allelic variant prevailed for CLC. Targeted sequencing of TCP, CLC, and MET1 regions on a representative set of trees confirmed the sex-associated allele-specific expression of the CLC gene in generative and vegetative tissues of P. × sibirica. Our study brings new knowledge on sex-associated differences in Populus species.

Sex-specific polymorphism of MET1 and ARR17 genes in Populus × sibirica

The genus Populus is an effective model in tree genetics. This genus includes dioecious species and, recently, whole genome resequencing of P. trichocarpa and P. balsamifera enabled the identification of sex-linked regions and sex-associated single nucleotide polymorphisms (SNPs). These results created new opportunities to study sex determination in poplars. In the present work, we performed deep sequencing of genes encoding METHYLTRANSFERASE1 (MET1) and homolog of ARABIDOPSIS RESPONSE REGULATOR 17 (ARR17), which are localized in a sex-linked region of Populus genome and contain a number of sex-associated SNPs. Amplicon libraries for 38 samples of P. × sibirica (19 males and 19 females) were sequenced on MiSeq Illumina (300 nt paired-end reads) and approximately 4000× coverage was obtained for each sample. In total, from 80 to 179 SNPs were detected in poplar individuals for MET1, and from 16 to 49 SNPs were detected for ARR17. We identified 17 sex-specific SNPs (11 in MET1 and 6 in ARR17) - they were present in all males but absent in all females. For identified sex-specific SNP sites, females were homozygous, while males were heterozygous. Moreover, colocation of sex-specific SNPs confirming the XY sex-determination system of poplars was revealed: in one allelic variant, males had the same nucleotides as females, while in the other, sex-specific SNPs were present. Based on the data obtained, we developed and successfully applied a high-resolution melting-based approach for sex identification in poplars. The developed molecular markers are useful for distinguishing between male and female poplars in scientific research and can also be applied to select male-only genotypes for use in city landscaping and production of paper, pulp, and biofuel.

Genetics polymorphism of poplars from Moscow region based on high-throughput sequencing of ITS

Poplars are widely used in landscaping of Moscow due to the ability to effectively purify the air from harmful impurities and to release a large amount of oxygen. The genusPopulusis characterized by a high level of intraspecies polymorphism, as well as the presence of natural interspecies hybrids. The aim of our work was to evaluate the genetic diversity of poplars, which are growing on the territory of Moscow city by high-throughput sequencing of internal transcribed spacers of 45S rRNA genes (ITS sequences). Sequencing of ITS of 40 poplar plants was performed on Illumina platform (MiSeq) and about 3 000 reads were obtained for each sample in average. Bioinformatics analysis was performed using CLC Genomics Workbench tool. The involved set of poplars had a high level of genetic diversity – the number of single nucleotide polymorphisms (SNPs) detected in each genotype relative to the reference ITS1 and ITS2 sequences ofP. trichocarpavarying from 4 to 44. We showed that even trees which had been planted on the same territory and, probably, at the same time had significant genetic differences. It can be speculated that highly polymorphic plant material was used for planting poplars in Moscow. For some sites with SNPs, several variants of nucleotides were found in the same individual and the ratio of SNPs was different. We assume that close to 50/50 ratio is observed in interspecific hybrids due to genetic differences in the ITS sequences between maternal and paternal genotypes. For SNPs with a predominance of one of the variants, the presence of paralogues among numerous genomic copies of ITS sequences is more likely. The results of our work can provide a framework for molecular genetic markers application with the purpose ofPopulusspecies and interspecific hybrids identification, determination the origin of a number of natural hybrids, and monitoring the diversity of genusPopulusin the Moscow city.

PHENETIC ANALYSIS OF POPULUS NIGRA, P. LAURIFOLIA AND P. × JRTYSCHENSIS IN NATURAL HYBRIDIZATION ZONE

The wide spread of hybridization in the genus Populus, including spontaneous hybridization, caused by cultivars, requires studying the variability and inheritance of morphological traits by hybrids for initial tracking of these processes. The analysis of endogenous, intra- and inter-population variability was performed on 533 individual poplar trees in seven populations of P. nigra, seven populations of P. laurifolia and four populations of P. × jrtyschensis in the Tom river basin. On each specimen, 15 leaves from short mid-crown branches were collected to determine the shape of the leaf blade, the shape of its tip and base, as well as the branch morphotype. Some biometric indicators were proposed for geometric assessment of the leaf blade shapes of poplar species. The analysis showed that of all the traits examined only the leaf blade shape did not meet the criterion for “phene”, since it is usually represented by several forms in the crown of one and the same tree. In all the species studied, the level of their intra-population diversity was found to be much higher than the inter-population one. According to the increase of intra-population variability of qualitative traits, the taxa could be ranked as P. nigra < P. laurifolia < P. × jrtyschensis. The share of inter-population diversity differed among the species studied, accounting for 21.5% in P. laurifolia, 13.8% in P. nigra and 8.0% in P. × jrtyschensis. The P. laurifolia populations showed the greatest inter-population differentiation, most likely because of a disjunct distribution due to narrow specialization to the montane river environment. The lower differentiation in P. nigra is probably due to the facts that this species dominates the poplar stands of the Tom River basin; its populations are not fragmented and are linked by vast gene flow. In P. nigra, an increase in the diversity of qualitative characteristics and phenotypes was observed in populations confined to hybridization centers. Natural selection is most likely the factor governing the inter-population differentiation in P. × jrtyschensis, leading to the predominance of F1 hybrids in populations and hence to a sharp decrease in inter-population variability.

An Empirical Assessment of Transgene Flow from a Bt Transgenic Poplar Plantation

To assess the possible impact of transgenic poplar plantations on the ecosystem, we analyzed the frequency and distance of gene flow from a mature male transgenic Populus nigra plantation carrying the Bacillus thuringiensis toxin gene (Bt poplar) and the survival of Bt poplar seeds. The resultant Bt poplar seeds occurred at a frequency of ~0.15% at 0 m to ~0.02% at 500 m from the Bt poplar plantation. The germination of Bt poplar seeds diminished within three weeks in the field (germination rate from 68% to 0%) compared to 48% after three weeks of storage at 4°C. The survival rate of seedlings in the field was 0% without any treatment but increased to 1.7% under the addition of four treatments (cleaning and trimming, watering, weeding, and covering with plastic film to maintain moisture) after being seeded in the field for eight weeks. The results of this study indicate that gene flow originating from the Bt poplar plantation occurred at an extremely low level through pollen or seeds under natural conditions. This study provides first-hand field data on the extent of transgene flow in poplar plantations and offers guidance for the risk assessment of transgenic poplar plantations.

Fitness dynamics within a poplar hybrid zone: I. Prezygotic and postzygotic barriers impacting a native poplar hybrid stand

Hybridization and introgression are pervasive evolutionary phenomena that provide insight into the selective forces that maintain species boundaries, permit gene flow, and control the direction of evolutionary change. Poplar trees (Populus L.) are well known for their ability to form viable hybrids and maintain their distinct species boundaries despite this interspecific gene flow. We sought to quantify the hybridization dynamics and postzygotic fitness within a hybrid stand of balsam poplar (Populus balsamifera L.), eastern cottonwood (P. deltoides Marsh.), and their natural hybrids to gain insight into the barriers maintaining this stable hybrid zone. We observed asymmetrical hybrid formation with P. deltoides acting as the seed parent, but with subsequent introgression biased toward P. balsamifera. Native hybrids expressed fitness traits intermediate to the parental species and were not universally unfit. That said, native hybrid seedlings were absent from the seedling population, which may indicate additional selective pressures controlling their recruitment. It is imperative that we understand the selective forces maintaining this native hybrid zone in order to quantify the impact of exotic poplar hybrids on this native system.