Subgenome evolutionary dynamics in allotetraploid ferns: insights from the gene expression patterns in the allotetraploid species Phegopteris decursivepinnata (Thelypteridacea, Polypodiales)

Allopolyploidization often leads to disruptive conflicts among more than two sets of subgenomes, leading to genomic modifications and changes in gene expression. Although the evolutionary trajectories of subgenomes in allopolyploids have been studied intensely in angiosperms, the dynamics of subgenome evolution remain poorly understood in ferns, despite the prevalence of allopolyploidization. In this study, we have focused on an allotetraploid fern-Phegopteris decursivepinnata-and its diploid parental species, P. koreana (K) and P. taiwaniana (T). Using RNA-seq analyses, we have compared the gene expression profiles for 9,540 genes among parental species, synthetic F1 hybrids, and natural allotetraploids. The changes in gene expression patterns were traced from the F1 hybrids to the natural allopolyploids. This study has revealed that the expression patterns observed in most genes in the F1 hybrids are largely conserved in the allopolyploids; however, there were substantial differences in certain genes between these groups. In the allopolyploids compared with the F1 hybrids, the number of genes showing a transgressive pattern in total expression levels was increased. There was a slight reduction in T-dominance and a slight increase in K-dominance, in terms of expression level dominance. Interestingly, there is no obvious bias toward the T- or K-subgenomes in the number and expression levels overall, showing the absence of subgenome dominance. These findings demonstrated the impacts of the substantial transcriptome change after hybridization and the moderate modification during allopolyploid establishment on gene expression in ferns and provided important insights into subgenome evolution in polyploid ferns.

Global patterns and climatic determinants of phylogenetic structure of regional fern floras

Knowledge of relationships between phylogenetic structure of a biological assemblage and ecological factors that drive the variation of phylogenetic structure among regions is crucial for understanding the causes of variation in taxonomic composition and richness among regions, but this knowledge is lacking for the global flora of ferns. Here, we fill this critical knowledge gap. We divided the globe into 392 geographic units on land, collated species lists of ferns for each geographic unit, and used different phylogenetic metrics (tip- vs basal-weighted) reflecting different evolutionary depths to quantify phylogenetic structure. We then related taxonomic and phylogenetic structure metrics to six climatic variables for ferns as a whole and for two groups of ferns (old clades vs polypods) reflecting different evolutionary histories across the globe and within each continental region. We found that when old clades and polypods were considered separately, temperature-related variables explained more variation in these metrics than did precipitation-related variables in both groups. When analyses were conducted for continental regions separately, this pattern holds in most cases. Climate extremes have a stronger relationship with phylogenetic structure of ferns than does climate seasonality. Climatic variables explained more variation in phylogenetic structure at deeper evolutionary depths.

An open and continuously updated fern tree of life

Ferns, with about 12,000 species, are the second most diverse lineage of vascular plants after angiosperms. They have been the subject of numerous molecular phylogenetic studies, resulting in the publication of trees for every major clade and DNA sequences from nearly half of all species. Global fern phylogenies have been published periodically, but as molecular systematics research continues at a rapid pace, these become quickly outdated. Here, we develop a mostly automated, reproducible, open pipeline to generate a continuously updated fern tree of life (FTOL) from DNA sequence data available in GenBank. Our tailored sampling strategy combines whole plastomes (few taxa, many loci) with commonly sequenced plastid regions (many taxa, few loci) to obtain a global, species-level fern phylogeny with high resolution along the backbone and maximal sampling across the tips. We use a curated reference taxonomy to resolve synonyms in general compliance with the community-driven Pteridophyte Phylogeny Group I classification. The current FTOL includes 5,582 species, an increase of ca. 40% relative to the most recently published global fern phylogeny. Using an updated and expanded list of 51 fern fossil constraints, we find estimated ages for most families and deeper clades to be considerably older than earlier studies. FTOL and its accompanying datasets, including the fossil list and taxonomic database, will be updated on a regular basis and are available via a web portal (https://fernphy.github.io) and R packages, enabling immediate access to the most up-to-date, comprehensively sampled fern phylogeny. FTOL will be useful for anyone studying this important group of plants over a wide range of taxonomic scales, from smaller clades to the entire tree. We anticipate FTOL will be particularly relevant for macroecological studies at regional to global scales and will inform future taxonomic systems with the most recent hypothesis of fern phylogeny.

Antheridiogen controls spatial dynamics of sex expression in naturally occurring gametophytes of the tree fern Cyathea multiflora

PREMISE

Antheridiogen systems are a set of pheromonal mechanisms that control sex expression in fern gametophytes. However, antheridiogen has rarely been studied outside of the laboratory, and little is known about its function in natural settings. Combining predictions based on field and laboratory study, we tested whether the sexual structure of gametophytic colonies of a tree fern were attributable to antheridiogen.

METHODS

Gametophytic colonies of the antheridiogen-producing tree fern Cyathea multiflora were collected at La Selva Biological Station in Costa Rica in January 2019. The sex of each gametophyte was determined, mapped, and spatial statistic approaches were used to examine the distribution of sex in each colony.

RESULTS

In all gametophytic colonies, males were most common, representing 62-68% of individuals. No hermaphroditic gametophytes were identified in any colony. A quadrat-based method showed female gametophytes were not clustered in each colony, while male gametophytes were clustered. In two of the colonies, the K(r) test statistic for males was greater than expected compared to random simulations of sex expression, indicating male sex expression was spatially associated with females.

CONCLUSIONS

This study provides the first documentation of spatial sex expression in natural settings of gametophytes of an antheridiogen-producing tree fern species. The profound impact of antheridiogen on gametophytic sex expression in field settings suggests this system is intimately tied to mating system, fitness, and genetic diversity in Cyathea multiflora.

An experimental assessment of competitive interactions between sexual and apomictic fern gametophytes using Easy Leaf Area

PREMISE

Few studies have explored competition in fern gametophyte populations. One limiting factor is the tedious measurement of gametophyte size as a proxy for biomass in these small plants. Here, an alternative approach of estimating the number of green pixels from photos was employed to measure the competitive interactions among apomictic and sexual Dryopteris gametophytes.

METHODS

We cultivated the gametophytes of two apomictic (diploid and triploid) and one sexual (tetraploid) Dryopteris species in monocultures and in two-species mixtures in the ratios 1 : 1 and 1 : 3. The total gametophyte cover of each population originating from 20 spores was assessed using Easy Leaf Area. Assessments were performed weekly between weeks 4 and 10 of cultivation. Additionally, during week 5, the cover of each species in each mixture was estimated separately.

RESULTS

We identified a positive correlation between gametophyte size and ploidy level as well as sexual reproduction. The performance of the tested species in mixtures was dependent on the competitor species identity, indicating the importance of competition between gametophytes.

DISCUSSION

The methods outlined can be used for a rapid assessment of fern gametophyte cover in large gametophyte populations. Ploidy level and reproduction type seem to play a major role in the competitive abilities of fern gametophytes, but more research is needed on this topic.

Antinociceptive and Anti-Inflammatory Effects of Saline Extract and Lectin-Rich Fraction from Microgramma vacciniifolia Rhizome in Mice

Previous studies have characterized a saline extract from Microgramma vacciniifolia rhizome and its lectin (MvRL)-rich fraction with low acute toxicity. In the present study, we evaluated these preparations for acute toxicity (1,000 mg/kg) and antinociceptive and anti-inflammatory activities (100-400 mg/kg for the extract and 25-50 mg/kg for the fraction). No signs of toxicity were observed. Both the extract and fraction increased the latency period for nociception in the hot plate assay, decreased writhing induced by acetic acid, and promoted analgesic effects in phases 1 and 2 of the formalin test. The antinociceptive mechanism was attributed to interactions with opioid receptors and K⁺ ATPase channels. The extract and fraction decreased carrageenan-induced paw edema in 46.15 % and 77.22 %, respectively, at the highest doses evaluated. Furthermore, the fraction was shown to act on the bradykinin pathway. The ability to decrease leukocyte migration after treatment was also verified in the peritonitis and air pouch models. In exudates collected from air pouches, decreased tumor necrosis factor (TNF)-α and increased interleukin (IL)-10 levels were noted. Both the extract and fraction also effectively inhibited the development of granulomatous tissue. In conclusion, the substances investigated in this study can be used for the development of novel therapeutic options for pain and inflammatory processes.

Effects of graywater on the growth and survival of ornamental plants in nature-based systems

The current paper investigates the development of two ornamental plants, canna lily (Canna x generalis) and giant horsetail (Equisetum giganteum), at both bench and pilot scale. Combinations of gravel-filled mesocosm, planted and unplanted (control), irrigated with light greywater (GW_L) or tap water (W_T), were used. Both species were able to grow under the tested conditions with no indication of toxicity that could affect the development. Irrigation with GW_L, resulted in higher evapotranspiration (2.2 mm-2.8 mm) in canna lily than giant horsetail (1.7 mm-2.3 mm) in mesocosm system. When the plants were mature and the season was more humid and warmer, canna lily and giant horsetail irrigated with GW_L evapotranspirated 69.23% and 30.77%, respectively as compared to the unplanted GW_L-irrigated-mesocosm. Principal components and cluster analysis identified similarity between evapotranspiration (ET) and the characteristics of the plants. Both species can thus be used in constructed wetlands taking into consideration elements such as the space available, level of water and solar incidence so as to allow the full development of the plants. The roots of giant horsetail require high water availability. Low solar incidence is indicated for giant horsetail, and the opposite for canna lily, if flowering is desired.

Functional diversity and convergence in the evolution of plant reproductive structures

BACKGROUND AND AIMS

Structures that simultaneously perform many functional roles are likely to show a variety of morphological solutions to these demands, and thus probably exhibit high morphological disparity. In contrast, specialization for a few simple functions should result in a more limited suite of morphologies. We explore this idea using lycopsid reproductive structures, which, throughout their history, have performed a limited set of functional roles compared with the reproductive structures of other plant groups such as seed plants.

METHODS

We scored living and fossil lycopsid taxa for 18 discrete character measurements and several continuous traits, including sporangium size, supporting axis diameter, and strobilus length and width. We used the discrete characters to construct a multivariate morphospace for lycopsid reproductive morphology through time, and the continuous characters to test whether fossil and extant lycopsids show similar patterns of tissue allocation within reproductive structures.

RESULTS

Lycopsids occupy similar areas of reproductive morphospace and show similar patterns of tissue allocation over most of their history, alternating between diffuse fertile zones with leaf-like sporophylls and compact strobili with specialized sporophylls that allow sporangia to be closely packed while also protected during their development. Growth habit also plays an important role in lycopsid reproductive evolution, broadly influencing the size and shape of reproductive structures.

CONCLUSIONS

Lycopsid reproductive structures are primarily specialized for densely packaging sporangia, and are consistent with the idea that performing limited functional roles is associated with reduced morphological disparity. Morphologies similar to lycopsid strobili are also found in other groups with simple, wind-dispersed propagules, suggesting that the same processes occur across plant lineages.

Seasonal changes in tissue-water relations for eight species of ferns during historic drought in California

PREMISE OF THE STUDY

California experienced severe drought between 2012 and 2016. During this period, we compared seasonal changes in tissue-water relations among eight fern species in the Santa Monica Mountains of southern California to elucidate differential mechanisms of drought survival and physiological performance during extreme water deficits.

METHODS

We monitored seasonal changes in water potential (Ψmd) and dark-adapted chlorophyll fluorescence (Fv/Fm), assessed tissue-water relations including osmotic potential at saturation and the turgor loss point (Ψπ, sat and Ψπ, tlp), and measured, for two evergreen species, xylem-specific and leaf-specific hydraulic conductivity (Ks and Kl) and vulnerability of stem xylem to water stress-induced embolism (water potential at 50% loss hydraulic conductivity, Ψ50).

KEY RESULTS

Species grew in either riparian or chaparral understory. The five chaparral species had a wider range of seasonal water potentials, root depths, and frond phenological traits, including one evergreen, two summer-deciduous, and two desiccation-tolerant (resurrection) species. Evergreen species were especially diverse, with an evergreen riparian species maintaining seasonal water potentials above -1.3 MPa, while an evergreen chaparral species had seasonal water potentials below -8 MPa. In those two species the Ψ50 values were -2.5 MPa and -4.3 MPa, respectively.

CONCLUSIONS

Observed differences in physiological performance among eight fern species reflected niche partitioning in water utilization and habitat preference associated with distinct phenological traits. We predict differential survival among fern species as future drought events in California intensify, with desiccation-tolerant resurrection ferns being the most resistant.

Air humidity as key determinant of morphogenesis and productivity of the rare temperate woodland fern Polystichum braunii

(1) Most ferns are restricted to moist and shady habitats, but it is not known whether soil moisture or atmospheric water status are decisive limiting factors, or if both are equally important. (2) Using the rare temperate woodland fern Polystichum braunii, we conducted a three-factorial climate chamber experiment (soil moisture (SM) × air humidity (RH) × air temperature (T)) to test the hypotheses that: (i) atmospheric water status (RH) exerts a similarly large influence on the fern's biology as soil moisture, and (ii) both a reduction in RH and an increase in air temperature reduce vigour and growth. (3) Nine of 11 morphological, physiological and growth-related traits were significantly influenced by an increase in RH from 65% to 95%, leading to higher leaf conductance, increased above- and belowground productivity, higher fertility, more epidermal trichomes and fewer leaf deformities under high air humidity. In contrast, soil moisture variation (from 66% to 70% in the moist to ca. 42% in the dry treatment) influenced only one trait (specific leaf area), and temperature variation (15 °C versus 19 °C during daytime) only three traits (leaf conductance, root/shoot ratio, specific leaf area); RH was the only factor affecting productivity. (4) This study is the first experimental proof for a soil moisture-independent air humidity effect on the growth of terrestrial woodland ferns. P. braunii appears to be an air humidity hygrophyte that, whithin the range of realistic environmental conditions set in this study, suffers more from a reduction in RH than in soil moisture. A climate warming-related increase in summer temperatures, however, seems not to directly threaten this endangered species.

Development of microsatellite markers for the apomictic triploid fern Myriopteris lindheimeri (Pteridaceae)

PREMISE OF THE STUDY

Microsatellite markers were developed for investigating the population dynamics of Myriopteris lindheimeri (Pteridaceae), an apomictic triploid fern endemic to deserts of the southwestern United States and Mexico.

METHODS AND RESULTS

Using 454 sequencing, 21 microsatellite markers were developed. Of these, 14 were polymorphic with up to five alleles per locus and eight markers amplified in one or more congeneric close relatives (M. covillei, M. fendleri, M. aurea, and M. rufa). To demonstrate marker utility, M. lindheimeri samples from three Arizona populations were genotyped at nine loci. For each population, diversity measures including percent polymorphic loci, frequency of heterozygotes across all loci, and genotypic diversity were calculated. Across the three populations, on average, 63% of loci were polymorphic, the average frequency of heterozygotes (across all loci) was 0.32, and average genotypic diversity was 0.34.

CONCLUSIONS

These markers provide a foundation for future studies exploring polyploidy and apomixis in myriopterid ferns.

Mapping nutrient resorption efficiencies of subarctic cryptogams and seed plants onto the Tree of Life

Nutrient resorption from senescing photosynthetic organs is a powerful mechanism for conserving nitrogen (N) and phosphorus (P) in infertile environments. Evolution has resulted in enhanced differentiation of conducting tissues to facilitate transport of photosynthate to other plant parts, ultimately leading to phloem. Such tissues may also serve to translocate N and P to other plant parts upon their senescence. Therefore, we hypothesize that nutrient resorption efficiency (RE, % of nutrient pool exported) should correspond with the degree of specialization of these conducting tissues across the autotrophic branches of the Tree of Life. To test this hypothesis, we had to compare members of different plant clades and lichens within a climatic region, to minimize confounding effects of climatic drivers on nutrient resorption. Thus, we compared RE among wide-ranging basal clades from the principally N-limited subarctic region, employing a novel method to correct for mass loss during senescence. Even with the limited numbers of species available for certain clades in this region, we found some consistent patterns. Mosses, lichens, and lycophytes generally showed low REN (<20%), liverworts and conifers intermediate (40%) and monilophytes, eudicots, and monocots high (>70%). REP appeared higher in eudicots and liverworts than in mosses. Within mosses, taxa with more efficient conductance also showed higher REN. The differences in REN among clades broadly matched the degree of specialization of conducting tissues. This novel mapping of a physiological process onto the Tree of Life broadly supports the idea that the evolution of conducting tissues toward specialized phloem has aided land plants to optimize their internal nitrogen recycling. The generality of evolutionary lines in conducting tissues and nutrient resorption efficiency needs to be tested across different floras in different climatic regions with different levels of N versus P availability.

Genetic consequence of restricted habitat and population decline in endangered Isoetes sinensis (Isoetaceae)

BACKGROUND AND AIMS

Isoetes sinensis (Isoeteaceae) is a critically endangered aquatic quillwort in eastern China. Rapid decline of extant population size and local population extinction have occurred in recent years and have raised great concerns among conservationists.

METHODS

Amplified fragment length polymorphisms (AFLPs) were used to investigate the genetic variation and population structure of seven extant populations of the species.

KEY RESULTS

Eight primer combinations produced a total of 343 unambiguous bands of which 210 (61.2 %) were polymorphic. Isoetes sinensis exhibited a high level of intra-population genetic diversity (H(E) = 0.118; hs = 0.147; I = 0.192; P = 35.2 %). The genetic variation within each of the populations was not positively correlated with their size, suggesting recent population decline, which is well in accordance with field data of demographic surveys. Moreover, a high degree of genetic differentiation (F(ST) = 0.535; G(ST) = 0.608; theta(B) = 0.607) was detected among populations and no correlation was found between geographical and genetic distance, suggesting that populations were in disequilibrium of migration-drift. Genetic drift played a more important role than gene flow in the current population genetic structure of I. sinensis because migration of I. sinensis is predominantly water-mediated and habitat range was highly influenced by environment changes.

CONCLUSIONS

Genetic information obtained in the present study provides useful baseline data for formulating conservation strategies. Conservation management, including both reinforcement for in situ populations and ex situ conservation programmes should be carefully designed to avoid the potential risk of outbreeding depression by admixture of individuals from different regions. However, translocation within the same regional population should be considered as a measure of genetic enhancement to rehabilitate local populations. An ex situ conservation strategy for conserving all extant populations to maximize genomic representation of the species is also recommended.