Seed and pollen gene dispersal in Taxus baccata, a dioecious conifer in the face of strong population fragmentation

Sexual interference revealed by joint study of male and female pollination success in chestnut

Most seed plants produce both pollen and ovules. In principle, pollen export could interfere with pollen import through self-pollination, resulting in ovule usurpation and reduced fruit set. Evidence for such interference exists under experimental settings but its importance under natural conditions is unknown. To test for sexual interference in nature, it is necessary to study together mating system, through paternity analyses, and fruit set, the proportion of flowers giving seeds or fruits. We developed a new model combining both processes, using chestnut (Castanea) as case study. We carried out a paternity analysis in an intensively studied plot of 273 trees belonging to three interfertile chestnut species and including a range of individuals with more or less functional stamens, resulting in a large data set of 1924 mating events. We then measured fruit set on 216 of these trees. Fruit set of male-fertile trees was much lower than that of male-sterile trees. Our process-based model shows that pollen is not limiting in the study site and hence cannot account for reduced fruit set. It also indicates that self-pollination is high (74%) but selfing rate is low (4%). Self-pollen is less competitive than cross-pollen, reducing sexual interference, but not sufficiently, as many ovules end up being self-fertilized, 95% of which abort before fruit formation, resulting in the loss of 46% of the fruit crop. These results suggest that the main cause of reduced reproductive potential in chestnut is sexual interference by self-pollen, raising questions on its evolutionary origins.

NMπ 2.0: Software update to minimize the risk of false positives among determinants of reproductive success

In plant populations, parentage analysis helps understand factors shaping individual reproductive success. However, estimating reproductive success determinants based on parentage counts requires decoupling the effects of individual fecundity and propagule dispersal. The neighbourhood model implemented in the NMπ software provides a standard solution for this problem based on the fixed-effects regression-like approach. Nonetheless, it has been recently shown that the method is prone to false discoveries when important fecundity determinants are omitted. To account for the unexplained variance in fecundity, the Bayesian approach was developed based on the new model (the hierarchical neighbourhood model; HNM). Here, I present the NMπ software update that allows the HNM approach to be used in the framework of a friendly interface. More importantly, the HNM approach is now made available for both dispersed (seedlings) and nondispersed (seeds with known mothers) progeny data. The Bayesian approach, among others, selects significant fecundity determinants, estimates the proportion of variance in reproductive potential explained by selected determinants (R² ), and provides individual female and male fecundity values. Although the software was designed to handle microsatellite marker data, a solution is proposed for large sets of single nucleotide polymorphisms. The program can be run on Windows (using either a terminal or a graphical interface) as well as (using a terminal) on Linux, or macOS platforms. In any case, NMπ can utilize multicore processors to speed up the analysis. The updated package containing the code, the executable file, the user manual, and example data is available at https://www.ukw.edu.pl/pracownicy/plik/igor_chybicki/3694/.

Impacts of Flowering Density on Pollen Dispersal and Gametic Diversity Are Scale Dependent

AbstractPollen dispersal is a key evolutionary and ecological process, but the degree to which variation in the density of concurrently flowering conspecific plants (i.e., coflowering density) shapes pollination patterns remains understudied. We monitored coflowering density and corresponding pollination patterns of the insect-pollinated palm Oenocarpus bataua in northwestern Ecuador and found that the influence of coflowering density on these patterns was scale dependent: high neighborhood densities were associated with reductions in pollen dispersal distance and gametic diversity of progeny arrays, whereas we observed the opposite pattern at the landscape scale. In addition, neighborhood coflowering density also impacted forward pollen dispersal kernel parameters, suggesting that low neighborhood densities encourage pollen movement and may promote gene flow and genetic diversity. Our work reveals how coflowering density at different spatial scales influences pollen movement, which in turn informs our broader understanding of the mechanisms underlying patterns of genetic diversity and gene flow within populations of plants.

Study of the Pollen Grain Metabolome under Deposition of Nitrogen and Phosphorus in Taxus baccata L. and Juniperus communis L

Nitrogen plays an important role in both quantitative and qualitative aspects of plant reproduction, including pollen grain compounds and seed production. Recent studies have pointed out that pollen grains produced by male plants of T. baccata and J. communis subjected to a long period of fertilizer supplementation have lower in vitro germination ability and higher nitrogen content. To gain molecular insights into these observations, we conducted GC-MS analysis of both species to characterize the metabolomes of dry, mature pollen grains, which allowed for the identification and quantification of more than 200 metabolites. The results demonstrated that fertilizer supplementation impacts the relative content of 14 metabolites in J. communis (9 downregulated and 5 upregulated) and 21 in T. baccata (6 downregulated and 15 upregulated). Although plants showed little similarity in patterns, in metabolite profiles, both up and down fold-changes were observed. This is the first report on the gymnosperm pollen grain metabolomic profile and changes induced by long-term nitrogen and phosphorus supplementation. Pollen grains produced by fertilizer-supplemented male individuals had significantly lower relative content of linolenic acid, 5,6-dihydrouracil, maltotriose, galactonic acid, D-xylulose, and glycerol-α-phosphate but higher content of sorbitol, glucosamine, and 1,5-anhydro-D-glucitol as well as n-acetyl-d-hexosamine, dimethyl phthalate, glycine, galactose-6-phosphate, D-fructose-6-phosphate, pyroglutamic acid, and 3-(3-hydroxyphenyl)-3-hydroxypropionic acid. Thus, in pollen grain samples earlier shown to have different germination abilities, the presence of different metabolites indicates a significant environmental impact on the quality of gymnosperm pollen grains.

Influence of Pollen Dispersal and Mating Pattern in Domestication of Intermediate Wheatgrass, a Novel Perennial Food Crop

Intermediate wheatgrass (IWG) is a perennial forage grass that is currently being domesticated as a grain crop. It is a primarily wind-pollinated outcrossing species and expresses severe inbreeding depression when self-pollinated. Characterization of pollen dispersal, mating parameters, and change in genetic diversity due to pollen movement is currently lacking in IWG. In this study, we examined pollen dispersal in an IWG selection nursery by evaluating 846 progeny from 15 mother plants and traced their parentage to 374 fathers. A set of 2,500 genomic loci was used to characterize the population. We assigned paternity to 769 (91%) progeny and the average number of fathers per mother plant was 37, from an average of 56 progeny examined per mother. An extensive number (80%) of pollination events occurred within 10 m of the mother plants. Pollination success was not correlated with trait attributes of the paternal genotypes. Mating system analysis confirmed that IWG is highly outcrossing and inbreeding was virtually absent. Neither genetic diversity nor the genome-estimated trait values of progeny were significantly affected by pollinator distance. The distance of pollinator in an IWG breeding nursery therefore was not found to be a major contributor in maintaining genetic diversity. These findings reveal the pollen dispersal model in IWG for the first time and its effect on genetic diversity, which will be valuable in designing future IWG breeding populations. Information generated and discussed in this study could be applied in understanding gene flow and genetic diversity of other open-pollinated species.

The Genetic Differences and Structure of Selected Important Populations of the Endangered Taxus baccata in the Czech Republic

Taxus baccata L. (common yew) is an endangered tree species in the Czech Republic. However, its natural occurrence has not been adequately protected in all areas of the country. The aim of this study is to determine whether the yew population in the newly established Mařeničky seed orchard (TS_L) enables mixing with other Czech yew populations. Using a set of nuclear microsatellites, the genetic diversity in the Lužické Mountains (TS_L) and in selected wild-provenance populations from the Czech Republic (Jílovské yews, TS_J; Březinské yews, TS_B, and yews from Moravský Karst, TS_M) was studied, as they could be donor sources for potential translocation activities. We observed that the level of genetic diversity within the four Czech yew units that were investigated was high. An analysis of the molecular variance (AMOVA) showed 7% variation among populations, and the genetic differentiation values were low to moderate (FST = 0.042–0.108). According to a STRUCTURE analysis, high genetic similarity was observed between the TS_L and TS_B units. Our results provide important genetic suggestions on how conservation management can be designed to maximize its success.

Estimating dispersal using close kin dyads: The kindisperse R package

Investigating dispersal in animal populations can be difficult, particularly for taxa that are hard to directly observe such as those that are small or rare. A promising solution may come from new approaches that use genome-wide sequence data to detect close kin dyads and estimate dispersal parameters from the distribution of these dyads. These methods have so far only been applied to mosquito populations. However, they should have broad applicability to a range of taxa, although no assessment has yet been made on their performance under different dispersal conditions and study designs. Here we develop an R package and shiny app, kindisperse, that can be used to estimate dispersal parameters from the spatial distribution of close kin. kindisperse can handle study designs that target different life stages and allows for a range of dispersal kernel shapes and organismal life histories; we provide implementation examples for a vertebrate (Antechinus) and an invertebrate (Aedes). We use simulations run in kindisperse to compare the performance of two published close kin methodologies, showing that one method produces unbiased estimates whereas the other produces downward-biased estimates. We also use kindisperse simulations to investigate how study design affects dispersal estimates, and we provide guidelines for the size and shape of sample sites as well as the number of close kin needed for accurate estimates. kindisperse is easily adaptable for application to a variety of research contexts ranging from invasive pests to threatened species where noninvasive DNA sampling can be used to detect close kin.

An Overview of the Practices and Management Methods for Enhancing Seed Production in Conifer Plantations for Commercial Use

Flowering, the beginning of the reproductive growth, is a significant stage in the growth and development of plants. Conifers are economically and ecologically important, characterized by straight trunks and a good wood quality and, thus, conifer plantations are widely distributed around the world. In addition, conifer species have a good tolerance to biotic and abiotic stress, and a stronger survival ability. Seeds of some conifer species, such as Pinus koraiensis, are rich in vitamins, amino acids, mineral elements and other nutrients, which are used for food and medicine. Although conifers are the largest (giant sequoia) and oldest living plants (bristlecone pine), their growth cycle is relatively long, and the seed yield is unstable. In the present work, we reviewed selected literature and provide a comprehensive overview on the most influential factors and on the methods and techniques that can be adopted in order to improve flowering and seed production in conifers species. The review revealed that flowering and seed yields in conifers are affected by a variety of factors, such as pollen, temperature, light, water availability, nutrients, etc., and a number of management techniques, including topping off, pruning, fertilization, hormone treatment, supplementary pollination, etc. has been developed for improving cone yields. Furthermore, several flowering-related genes (FT, Flowering locus T and MADS-box, MCMI, AGAMOUS, DEFICIENCES and SRF) that play a crucial role in flowering in coniferous trees were identified. The results of this study can be useful for forest managers and for enhancing seed yields in conifer plantations for commercial use.

Extensive pollen-mediated gene flow across intensively managed landscapes in an insect-pollinated shrub native to semiarid habitats

Our knowledge of the impact of landscape fragmentation on gene flow patterns is mainly drawn from tropical and temperate ecosystems, where landscape features, such as the distance of a tree to the forest edge, drive connectivity and mating patterns. Yet, the structure of arid and semiarid plant communities - with open canopies and a scattered distribution of trees - differs greatly from those that are well-characterized in the literature. As a result, we ignore whether the documented consequences of landscape fragmentation on plant mating and gene flow patterns also hold for native plant communities in arid and semiarid regions. We investigated the relative contribution of plant traits, pollinator activity, and individual neighbourhood in explaining variation in mating and gene flow patterns of an insect-pollinated semiarid arborescent shrub, Ziziphus lotus, at three sites embedded in highly altered agriculture landscapes. We used 14 SSRs, seed paternity analyses, and individual mixed effect mating models (MEMMi) to estimate the individual mating variables and the pollen dispersal kernel at each site. Individual spatial location, flower density, and floral visitation rate explained most of the variation of mating variables. Unexpectedly, individual correlated paternity was very low and shrubs surrounded by the most degraded matrix exhibited an increased fraction of pollen immigration and a high effective number of pollen donors per mother shrub. Overall, our results reveal that an active pollinator assemblage ensures highly efficient mating, and maintains pollen-mediated gene flow and notable connectivity levels, even in highly altered landscapes, potentially halting genetic isolation within and between distant sites.

Sex-specific spatial variation in fitness in the highly dimorphic Leucadendron rubrum

Sexual dimorphism in plants may emerge as a result of sex-specific selection on traits enhancing access to nutritive resources and/or to sexual partners. Here we investigated sex-specific differences in selection of sexually dimorphic traits and in the spatial distribution of effective fecundity (our fitness proxy) in a highly dimorphic dioecious wind-pollinated shrub, Leucadendron rubrum. In particular, we tested for the effect of density on male and female effective fecundity. We used spatial and genotypic data of parent and offspring cohorts to jointly estimate individual male and female effective fecundity on the one hand and pollen and seed dispersal kernels on the other hand. This methodology was adapted to the case of dioecious species. Explicitly modelling dispersal avoids the confounding effects of heterogeneous spatial distribution of mates and sampled seedlings on the estimation of effective fecundity. We also estimated selection gradients on plant traits while modelling sex-specific spatial autocorrelation in fecundity. Males exhibited spatial autocorrelation in effective fecundity at a smaller scale than females. A higher local density of plants was associated with lower effective fecundity in males but was not related to female effective fecundity. These results suggest sex-specific sensitivities to environmental heterogeneity in L. rubrum. Despite these sexual differences, we found directional selection for wider canopies and smaller leaves in both sexes, and no sexually antagonistic selection on strongly dimorphic traits in L. rubrum. Many empirical studies in animals similarly failed to detect sexually antagonistic selection in species expressing strong sexual dimorphism, and we discuss reasons explaining this common pattern.

Identification of determinants of pollen donor fecundity using the hierarchical neighborhood model

Individual differences in male reproductive success drive genetic drift and natural selection, altering genetic variation and phenotypic trait distributions in future generations. Therefore, identifying the determinants of reproductive success is important for understanding the ecology and evolution of plants. Here, based on the spatially explicit mating model (the neighborhood model), we develop a hierarchical probability model that links co-dominant genotypes of offspring and candidate parents with phenotypic determinants of male reproductive success. The model accounts for pollen dispersal, genotyping errors as well as individual variation in selfing, pollen immigration, and differentiation of immigrant pollen pools. Unlike the classic neighborhood model approach, our approach is specially designed to account for excessive variation (overdispersion) in male fecundity. We implemented a Bayesian estimation method (the Windows computer program available at: https://www.ukw.edu.pl/pracownicy/plik/igor_chybicki/1806/) that, among others, allows for selecting phenotypic variables important for male fecundity and assessing the fraction of variance in fecundity (R² ) explained by selected variables. Simulations showed that our method outperforms both the classic neighborhood model and the two-step approach, where fecundities and the effects of phenotypic variables are estimated separately. The analysis of two data examples showed that in wind-pollinated trees, male fecundity depends on both the amount of produced pollen and the ability to pollen spread. However, despite that the tree size was positively correlated with male fecundity, it explained only a fraction of the total variance in fecundity, indicating the presence of additional factors. Finally, case studies highlighted the importance of accounting for pollen dispersal in the estimation of fecundity determinants.

Landscape Genetics of Plants: Challenges and Opportunities

Dispersal is one of the most important but least understood processes in plant ecology and evolutionary biology. Dispersal of seeds maintains and establishes populations, and pollen and seed dispersal are responsible for gene flow within and among populations. Traditional views of dispersal and gene flow assume models that are governed solely by geographic distance and do not account for variation in dispersal vector behavior in response to heterogenous landscapes. Landscape genetics integrates population genetics with Geographic Information Systems (GIS) to evaluate the effects of landscape features on gene flow patterns (effective dispersal). Surprisingly, relatively few landscape genetic studies have been conducted on plants. Plants present advantages because their populations are stationary, allowing more reliable estimates of the effects of landscape features on effective dispersal rates. On the other hand, plant dispersal is intrinsically complex because it depends on the habitat preferences of the plant and its pollen and seed dispersal vectors. We discuss strategies to assess the separate contributions of pollen and seed movement to effective dispersal and to delineate the effects of plant habitat quality from those of landscape features that affect vector behavior. Preliminary analyses of seed dispersal for three species indicate that isolation by landscape resistance is a better predictor of the rates and patterns of dispersal than geographic distance. Rates of effective dispersal are lower in areas of high plant habitat quality, which may be due to the effects of the shape of the dispersal kernel or to movement behaviors of biotic vectors. Landscape genetic studies in plants have the potential to provide novel insights into the process of gene flow among populations and to improve our understanding of the behavior of biotic and abiotic dispersal vectors in response to heterogeneous landscapes.

Individual reproductive success in Norway spruce natural populations depends on growth rate, age and sensitivity to temperature

Quantifying the individual reproductive success and understanding its determinants is a central issue in evolutionary research for the major consequences that the transmission of genetic variation from parents to offspring has on the adaptive potential of populations. Here, we propose to distil the myriad of information embedded in tree-ring time series into a set of tree-ring-based phenotypic traits to be investigated as potential drivers of reproductive success in forest trees. By using a cross-disciplinary approach that combines parentage analysis and a thorough dendrophenotypic characterisation of putative parents, we assessed sex-specific relationships between such dendrophenotypic traits (i.e., age, growth rate and parameters describing sensitivity to climate and to extreme climatic events) and reproductive success in Norway spruce. We applied a full probability method for reconstructing parent-offspring relationships between 604 seedlings and 518 adult trees sampled within five populations from southern and central Europe. We found that individual female and male reproductive success was positively associated with tree growth rate and age. Female reproductive success was also positively influenced by the correlation between growth and the mean temperature of the previous vegetative season. Overall, our results showed that Norway spruce individuals with the highest fitness are those who are able to keep high-growth rates despite potential growth limitations caused by reproductive costs and climatic limiting conditions. Identifying such functional links between the individual ecophysiological behaviour and its evolutionary gain would increase our understanding on how natural selection shapes the genetic composition of forest tree populations over time.

Recent Fragmentation May Not Alter Genetic Patterns in Endangered Long-Lived Species: Evidence From Taxus cuspidata

Forestland fragmentation caused by overexploitation of forest resources can in principle reduce genetic diversity, limit gene flow and eventually lead to species developing strong genetic structure. However, the genetic consequences of recent anthropogenic fragmentation of tree species remain unclear. Taxus cuspidata, which has extremely small populations distributed mainly in Changbai Mt. in Northeast (NE) China, has recently endured severe habitat fragmentation. Here, we investigate the pattern of genetic diversity and structure, identify risk factors, predict the future distribution and finally provide guidelines for the conservation and management of this species. We used three chloroplast and two mitochondrial DNA fragments, which are both paternally inherited in yews but differ in mutation rates, to genotype a total of 265 individuals from 26 populations covering the distribution of the species in China. Both chloroplast and mitochondrial data showed high degrees of genetic diversity, extensive gene flow over the entire geographical range and historical stability of both effective population size and distribution of the species. However, ecological niche modeling suggests a decrease in suitable areas for this species by the years 2050 and 2070. The maintenance of high genetic diversity and the existence of sufficient gene flow suggest that recent fragmentation has not affected the genetic composition of the long-lived tree T. cuspidata. However, severe impacts of anthropogenic activities are already threatening the species. Conservation and management strategies should be implemented in order to protect the remnant populations.