Directional and stabilizing selection shaped morphological, reproductive, and physiological traits of the invader Solidago canadensis

Trait evolution in invasive plant species is important because it can impact demographic parameters key to invasion success. Invasive plant species often show phenotypic clines along geographic and climatic gradients. However, the relative contributions of natural selection and neutral evolutionary processes to phenotypic trait variation among populations of invasive plants remains unclear. A common method to assess whether a trait has been shaped by natural selection or neutral evolutionary processes is to compare the geographical pattern for the trait of interest to the divergence in neutral genetic loci (i.e., QST-FST comparisons). A redundancy analysis (RDA) can facilitate identification of putative agents of natural selection on the trait. Here, we employed both a QST-FST comparisons approach and RDA to infer whether natural selection shaped traits of invasive populations of S. canadensis in China and identify the potential environmental drivers of natural selection. We addressed two questions: (1) Does natural selection drive phenotypic trait variation among S. canadensis populations? (2) Do climatic, latitudinal, longitudinal, and altitudinal gradients drive patterns of genetic variation among S. canadensis populations? We found significant directional selection for several morphological and reproductive traits (i.e., QST >FST) and stabilizing selection for physiological traits (i.e., QST < FST). The RDA showed that stem biomass of S. canadensis was strongly positively correlated with longitude, while leaf width ratio and specific leaf area were significantly positively correlated with mean diurnal range. Stem biomass had a strong negative correlation with annual precipitation. Moreover, height of S. canadensis individuals was strongly positively correlated with altitude and precipitation of the wettest quarter. Precipitation seasonality that was associated with longitudinal shift in China likely selected for larger stem biomass in S. canadensis. Overall, these results suggest that longitudinal and altitudinal clines in climate exerted strong selection pressures that shaped phenotypic traits of S. canadensis.ExcelFunding provided by: National Natural Science Foundation of ChinaCrossref Funder Registry ID: http://dx.doi.org/10.13039/501100001809Award Number: 31850410484Funding provided by: National Natural Science Foundation of ChinaCrossref Funder Registry ID: http://dx.doi.org/10.13039/501100001809Award Number: 31270461Funding provided by: Ten Thousand Talent Program of Zhejiang Province*Crossref Funder Registry ID: Award Number: 2019R52043Funding provided by: National Key Research and Development Program of ChinaCrossref Funder Registry ID: http://dx.doi.org/10.13039/501100012166Award Number: 2016YFC1201100Funding provided by: Taizhou city 500 Talent Program*Crossref Funder Registry ID: Award Number:Common garden experiment In April 2013, we established a common-garden experiment on Linhai campus of Taizhou University (121˚17´E, 28˚87´N) to assess quantitative trait variation among S. canadensis individuals from the 14 populations. We obtained three similar-sized (ca. 15 cm) plantlets for each of the 12 families in the 14 populations from the pre-experiment cultivation of plants described above. The plantlets were grown in a common garden in the field in three blocks. Thus, the total number of experimental plants was 504: 14 populations × 12 individuals (representing 12 maternal families) per population × 3 blocks. Within a block, the 12 individuals were planted in 12 separate plots that each measured 1.5 m × 1.5 m. The individual plants were grown 30 cm apart from each other. Throughout, the experimental plants received only rain-fed water and were not fertilized. Phenotypic trait measurement in the common garden One month after transplant (21st - 27th May 2013), we took in situ measurements of photosynthesis on the third fully expanded leaf from the top using a portable photosynthesis meter (LI-6400 XT, Li-COR, Inc., Lincoln, NE, USA). The measurements were taken between 9:00 and 11:00 AM under the following conditions: a photosynthetically active radiation of 1,400 μmol m-2 s-1, leaf temperature of 25°C, CO2 concentration of 400 ppm, and relative humidity of 70%. For each plant, we recorded net photosynthetic rate (Pn), intercellular CO2 concentration (Ci), stomatal conductance (Gs), and transpiration rate (Tr). These measurements were taken six times on different dates, and the average of the six were used in the statistical analyses described below. Seven months after transplant (November 2013), we measured various morphological, physiological, and reproductive traits of the S. canadensis individuals. Individual plant height, leaf length (L), and leaf width (L) for the third leaf per plant were measured to a precision of 0.1 cm. We used the L and W dimensions to compute the L/W ratio. Chlorophyll content of the third leaf was measured using a portable chlorophyll meter. We also measured basal stem diameter to an accuracy of 0.01 cm. After 97% of the plants had flowered, we counted the total number of inflorescences per individual plant. Then after the plants had matured, we obtained a total seed count and 1000-seed weight per plant. We also measured the total leaf area per plant using a WinFOLIA computer image analysis system (Regent Instruments Inc, Quebec, Canada) for use in computation of specific leaf area (SLA). The above-ground plant parts were then harvested and separated into shoots, leaves, and seeds. Biomass of the stem, leaf, and seed were measured after the samples had been dried at 80 °C for 72 h. We measured biomass using an electronic balance to an accuracy of 0.01 g. To obtain total vegetative biomass per S. canadensis individual, we summed up oven-dried stem and leaf biomass of the individuals. We computed SLA as the ratio of leaf area to leaf dry biomass. A year later (November 2014), we repeated measurements of the same traits above from the same S. canadensis individuals. After the measurements had been taken, the whole root system for each S. canadensis individual was dug up and washed free of soil particles under running water. The roots and shoots were then dried at 80 °C for 72 h and then weighed individually to the nearest 0.01 g. Analysis of neutral genetic diversity Leaf sampling In October 2012, we collected leaf tissues from 30 S. canadensis individuals in the same 14 populations that we sampled for rhizomes as described above. We kept the minimum distance between the sampled S. canadensis individuals at 10 m to reduce the chance of sampling the same maternal family more than once. The sampled leaves were immediately immersed in self-sealing plastic bags that contained silica gel and then stored in the laboratory at room temperature until use in DNA extraction. Microsatellite analysis In November 2012, genomic DNA was extracted from 0.1 g of each leaf sample using a modified sodium dodecyl sulfate protocol on a FastPrep-24 Automated Lysis and Homogenization System (MP Biomedicals, Santa Ana, CA, USA). Total DNA concentration was determined with a NanoDrop 2000 Lite Spectrophotometers (ThermoFisher Scientific, Inc. Rockford, IL, USA). The DNA samples were then diluted to 10 ng/μL and stored at -20°C until use in simple sequence repeats (SSR) analyses. Five prime pairs (synthesized by Boshang Biotechonology Co., Ltd. in China) were used in SSR amplifications (Wieczorek and Geber 2002). We ran PCR in a 20 μL volume that was made up of 1 × PCR reaction buffer, 1.5 mM Mg2+, 4 ng template DNA, 0.2 μM each of forward and reverse primers, 0.2 mM 4 × dNTP mixture, and 1 U Taq polymerase (Promega Cooperation, Madison, WI, USA). We performed PCR amplifications using a PTC 220 Peltier Thermal Cycler (Bio-Rad Laboratories, Hercules, CA, USA) with the following settings: denaturation at 95 ºC for 5 min, followed by 34 cycles of 30 s at 95 ºC, 30 s at 58 ºC (﹣1 ºC per cycle), 45 s at 72 ºC, with a final elongation of 7 min at 72 ºC. We analyzed the PCR products using a Fragment Analyzer™ Automated CE System (Advanced Analytical Technologies, Inc, Ankeny, IA, USA) with an 80 cm-long capillary column. We used a DNF-900 35-500 bp ds DNA Reagent Kit in the analysis. We genotyped DNA fragments using PROSize® 2.0 Data Analysis Software based on the elution time compared with a size standard