Plastome phylogenomics and new fossil evidence from Dominican amber shed light on the evolutionary history of the Neotropical fern genus Pecluma.

PREMISE: Molecular studies based on chloroplast markers have questioned the monophyly of the fern genus Pecluma (Polypodioideae, Polypodiaceae), which has several species of Polypodium nested within it. We explored the delimitation of Pecluma and its biogeographic pattern by evaluating the phylogenetic position of four Polypodium species not sequenced thus far and integrating the first fossil evidence of Pecluma. METHODS: Using herbarium material, we applied a genome-skimming approach to obtain a phylogenetic hypothesis of Polypodioideae; assessed the combination of character states observed in the fossil from Miocene Dominican amber using a previously published phylogeny of Polypodioideae based on four plastid markers as framework; calculated divergence times; and conducted an ancestral area estimation. RESULTS: Within Polypodioideae, Pecluma was recovered as sister to Phlebodium. Three of the newly sequenced species-Polypodium otites, P. pinnatissimum, and P. ursipes-were recovered with maximum support within the Pecluma clade, whereas P. christensenii remained within Polypodium. The closest combination of character states of the fossil was found within Pecluma. Our biogeographic analyses suggest an Eocene origin of the genus in South America, with several subsequent Oligocene and Miocene colonization events to Mexico-Central America and to the West Indies. CONCLUSIONS: Although the circumscription of Pecluma is still challenging, our results elucidate the origin and age of the genus. The newly described fossil, Pecluma hispaniolae sp. nov., supports the hypothesis that the epiphytic communities of the Greater Antilles exhibit a constant generic composition since the Miocene. We propose new combinations (Pecluma otites, Pecluma pinnatissima, and Pecluma ursipes) to accommodate three species previously classified in Polypodium.

Expression and Functional Studies of Leaf, Floral, and Fruit Developmental Genes in Non-model Species.

Researchers working on evolutionary developmental plant biology are inclined to choose non-model taxa to address how specific features have been acquired during ontogeny and fixed during phylogeny. In this chapter we describe methods to extract RNA, to assemble de-novo transcriptomes, to isolate orthologous genes within gene families, and to evaluate expression and function of target genes. We have successfully optimized these protocols for non-model plant species including ferns, gymnosperms, and a large assortment of angiosperms. In the latter, we have ranged a large number of families including Aristolochiaceae, Apodanthaceae, Chloranthaceae, Orchidaceae, Papaveraceae, Rubiaceae, Solanaceae, and Tropaeolaceae.

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.

DNA barcoding ferns in an unexplored tropical montane cloud forest area of southeast Oaxaca, Mexico.

DNA barcoding can be useful for species identification and phylogenetic analysis, but its effectivity has not been verified in most neotropical cloud forest plants. We tested three plastid barcodes, rbcLa, matK, and trnH-psbA, in selected pteridophytes, a well-represented group in these forests, from a little-explored area in Oaxaca, Mexico, applying the CBOL criteria for barcoding. We used BLASTn, genetic distance, and monophyly tree-based analyses employing neighbor-joining (NJ), maximum likelihood (ML), and Bayesian inference methods. Universal primers for rbcLa and trnH-psbA were successfully amplified and bi-directionally sequenced, but matK could not be amplified for most species. rbcLa showed the highest species discrimination in BLASTn (66.67%). trnH-psbA exhibited higher significant interspecific divergence values than rbcL and rbcLa + trnH-psbA (two-sample sign test, P value < 2.2e-16). Using NJ and ML phylogenetic trees, monophyletic species were successfully resolved (100%), differing only in support values and displaying full agreement with the most recent fern classification. ML trees showed the highest mean support value (80.95%). trnH-psbA was the only barcode that could detect the Elaphoglossoideae subfamily. Species discrimination did not increase using rbcLa + trnH-psbA. rbcLa is useful for fern barcoding, trnH-psbA is most helpful for phylogenetic analyses, and matK may not work as a universal barcoding marker.

State of the art in cytogenetics, insights into chromosome number evolution, and new C-value reports for the fern family Gleicheniaceae.

Studies concerning the cytogenetics of Gleicheniaceae have been scarce, especially those employing evolutionary approaches. Two chromosome number evolutionary models have been hypothesized for Gleicheniaceae. One proposes that ancestral haploid numbers were small and that the chromosome numbers of extant species evolved through polyploidy. The other model proposes that, at the genus level, fern chromosome evolution occurred from ancestors with essentially the same high chromosome numbers seen in living lineages. Neither of those hypotheses has been tested based on phylogenetic frameworks. We sought to (i) present the state of the art of Gleicheniaceae chromosome numbers; (ii) test the two evolutionary models of chromosome numbers within a phylogenetic framework; (iii) test correlations between DNA contents and chromosome numbers in the family. We report here DNA C-values for five species, which increases the number of investigated taxa nearly twofold and report two new genera records. Ancestral state chromosome reconstruction corroborates the hypothesis that ancestral chromosome numbers in Gleicheniaceae were as high as those of extant lineages. Our results demonstrate the possible role of dysploidy in the evolutionary chromosome history of Gleicheniaceae at the genus level and suggest that the relationship between chromosome number and DNA content does not appear to be linear.

A taxonomic and molecular survey of the pteridophytes of the Nectandra Cloud Forest Reserve, Costa Rica.

Floristic surveys are crucial to the conservation of biodiversity, but the vast majority of such surveys are limited to listing species names, and few take into account the evolutionary history of species. Here, we combine classical taxonomic and molecular phylogenetic (DNA barcoding) approaches to catalog the biodiversity of pteridophytes (ferns and lycophytes) of the Nectandra Cloud Forest Reserve, Costa Rica. Surveys were carried out over three field seasons (2008, 2011, and 2013), resulting in 176 species representing 69 genera and 22 families of pteridophytes. Our literature survey of protected areas in Costa Rica shows that Nectandra has an exceptionally diverse pteridophyte flora for its size. Plastid rbcL was selected as a DNA barcode marker and obtained for >95% of pteridophyte taxa at this site. Combined molecular and morphological analyses revealed two previously undescribed taxa that appear to be of hybrid origin. The utility of rbcL for species identification was assessed by calculating minimum interspecific distances and found to have a failure rate of 18%. Finally we compared the distribution of minimum interspecific rbcL distances with two other areas that have been the focus of pteridophyte molecular surveys: Japan and Tahiti. The comparison shows that Nectandra is more similar to Japan than Tahiti, which may reflect the biogeographic history of these floras.

A comparative gene co-expression analysis using self-organizing maps on two congener filmy ferns identifies specific desiccation tolerance mechanisms associated to their microhabitat preference.

BACKGROUND: Filmy-ferns (Hymenophyllaceae) are poikilohydric, homoiochlorophyllous desiccation-tolerant (DT) epiphytes. They can colonize lower and upper canopy environments of humid forest. Filmy-ferns desiccate rapidly (hours), contrasting with DT angiosperms (days/weeks). It has been proposed that desiccation tolerance in filmy-ferns would be associated mainly with constitutive features rather than induced responses during dehydration. However, we hypothesize that the inter-specific differences in vertical distribution would be associated with different dynamics of gene expression within the dehydration or rehydration phases. A comparative transcriptomic analysis with an artificial neural network was done on Hymenophyllum caudiculatum (restricted to lower canopy) and Hymenophyllum dentatum (reach upper canopy) during a desiccation/rehydration cycle. RESULTS: Raw reads were assembled into 69,599 transcripts for H. dentatum and 34,726 transcripts for H. caudiculatum. Few transcripts showed significant changes in differential expression (DE). H. caudiculatum had ca. twice DE genes than H. dentatum and higher proportion of increased-and-decreased abundance of genes occurs during dehydration. In contrast, the abundance of genes in H. dentatum decreased significantly when transitioning from dehydration to rehydration. According to the artificial neural network results, H. caudiculatum enhanced osmotic responses and phenylpropanoid related pathways, whilst H. dentatum enhanced its defense system responses and protection against high light stress. CONCLUSIONS: Our findings provide a deeper understanding of the mechanisms underlying the desiccation tolerance responses of two filmy ferns and the relationship between the species-specific response and the microhabitats these ferns occupy in nature.

Genes coding for transporters showed a rapid and sharp increase in their expression in response to lead, in the aquatic fern (Salvinia minima Baker).

Salvinia minima was assessed for its ability to accumulate lead (Pb) by exposing it to concentrations of 40µM Pb(NO3)2 during 24h. At the same time, the expression levels were quantified, of four genes coding for transporters: SmABCC (ABCC-MRP), SmATPase (ATPase-P3A), SmNhaD (Type-Na+/H+) and SmABCG (ABCG-WBC). In the absence of lead, S. minima had very low expression of those genes, when plants were exposed to the metal however, those genes showed a rapid (in just three hours or less) and sharp increase (up to 60 times) in their expression, particularly the SmNhaD (Type-Na+/H+) gene. This sharp increase in expression levels of the genes studied, occurred at the same time that the plant accumulated the highest content of lead in its tissues. The first two genes, are apparently implicated in detoxification and lead accumulation mechanisms, while the other two genes are apparently involved in maintaining cell balance (homeostatic control) and membrane integrity. Our results confirmed that S. minima is efficient for phytoremediation of water bodies contaminated by lead, as it is efficient in accumulating this metal in its tissues (bioconcentration factor; BCF) values greater than 1000, in short times of exposure. More importantly, our data on the expression profiles of four genes coding for transporters, represent a first sight scenario of the molecular basis for understanding the different mechanism of detoxification, apparently present in this aquatic fern.

Sequencing of small RNAs of the fern Pleopeltis minima (Polypodiaceae) offers insight into the evolution of the microrna repertoire in land plants.

MicroRNAs (miRNAs) are short, single stranded RNA molecules that regulate the stability and translation of messenger RNAs in diverse eukaryotic groups. Several miRNA genes are of ancient origin and have been maintained in the genomes of animal and plant taxa for hundreds of millions of years, playing key roles in development and physiology. In the last decade, genome and small RNA (sRNA) sequencing of several plant species have helped unveil the evolutionary history of land plants. Among these, the fern group (monilophytes) occupies a key phylogenetic position, as it represents the closest extant cousin taxon of seed plants, i.e. gymno- and angiosperms. However, in spite of their evolutionary, economic and ecological importance, no fern genome has been sequenced yet and few genomic resources are available for this group. Here, we sequenced the small RNA fraction of an epiphytic South American fern, Pleopeltis minima (Polypodiaceae), and compared it to plant miRNA databases, allowing for the identification of miRNA families that are shared by all land plants, shared by all vascular plants (tracheophytes) or shared by euphyllophytes (ferns and seed plants) only. Using the recently described transcriptome of another fern, Lygodium japonicum, we also estimated the degree of conservation of fern miRNA targets in relation to other plant groups. Our results pinpoint the origin of several miRNA families in the land plant evolutionary tree with more precision and are a resource for future genomic and functional studies of fern miRNAs.

Rates of ecological divergence and body size evolution are correlated with species diversification in scaly tree ferns.

Variation in species richness across regions and between different groups of organisms is a major feature of evolution. Several factors have been proposed to explain these differences, including heterogeneity in the rates of species diversification and the age of clades. It has been frequently assumed that rapid rates of diversification are coupled to high rates of ecological and morphological evolution, leading to a prediction that remains poorly explored for most species: the positive association between ecological niche divergence, morphological evolution and species diversification. We combined a time-calibrated phylogeny with distribution, ecological and body size data for scaly tree ferns (Cyatheaceae) to test whether rates of species diversification are predicted by the rates at which clades have evolved distinct ecological niches and body sizes. We found that rates of species diversification are positively correlated with rates of ecological and morphological evolution, with rapidly diversifying clades also showing rapidly evolving ecological niches and body sizes. Our results show that rapid diversification of scaly tree ferns is associated with the evolution of species with comparable morphologies that diversified into similar, yet distinct, environments. This suggests parallel evolutionary pathways opening in different tropical regions whenever ecological and geographical opportunities arise. Accordingly, rates of ecological niche and body size evolution are relevant to explain the current patterns of species richness in this 'ancient' fern lineage across the tropics.

Development of microsatellite markers in the tetraploid fern Ceratopteris thalictroides (Parkeriaceae) using RAD tag sequencing.

To understand the genetic variability of the tetraploid fern Ceratopteris thalictroides (Parkeriaceae), we described 30 polymorphic microsatellite markers obtained using the restriction site-associated DNA (RAD) tag sequencing technique. A total of 26 individuals were genotyped for each marker. The number of alleles per locus ranged from 4 to 10, and the expected heterozygosity and the Shannon-Wiener index ranged from 0.264 to 0.852 and 0.676 to 2.032, respectively. Because these 30 microsatellite markers exhibit high degrees of genetic variation, they will be useful tools for studying the adaptive genetic variation and sustainable conservation of C. thalictroides.

Molecular characterization of Platonia insignis Mart. ("bacurizeiro") using inter simple sequence repeat (ISSR) markers.

Platonia insignis Mart. (Clusiaceae) is widespread throughout the Amazon and adjacent areas. The fruits (known locally as "bacuri") have significant commercial potential, but the species is under threat from agro-industrial expansion. The genetic variability within 72 genotypes of P. insignis belonging to ten populations collected in the Brazilian states of Maranhão and Piauí, and maintained in the germplasm collection of Embrapa Meio-Norte, has been determined, and the organization of genetic diversity within populations, between populations and among geographic groups verified. Eighteen selected inter simple sequence repeat primers allowed amplification of 236 loci of which 221 (93.64%) were polymorphic, indicating a high level of genetic diversity. At the population level, the Shannon and Nei diversity indices ranged from 0.082 to 0.323 and from 0.120 to 0.480, respectively. The global coefficient of genetic differentiation (G(ST)) was 0.4730 indicating that differentiation between populations was significant, a finding that was confirmed by analysis of molecular variance (Φ(ST) = 0.28). UPGMA cluster analysis revealed that the genotypes could be stratified into groups that were well defined and consistent with those identified in the dendrogram constructed using pair wise Φ(ST) values. The high genetic diversity established in this study may facilitate the management and conservation of the germplasm of P. insignis.

Effects of 2,4-D on the germination of megaspores and initial development of Regnellidium diphyllum Lindman (Monilophyta, Marsileaceae) / Efeitos do 2,4-D sobre a germinação de megásporos e o desenvolvimento inicial de Regnellidium diphyllum Lindman (Monilophyta, Marsileaceae)

Regnellidium diphyllum is considered as endangered, occurring in the State of Rio Grande do Sul, Brazil, and a few adjoining localities in Uruguay, Argentina and the State of Santa Catarina. It grows in wetlands frequently altered for agricultural activities. Herbicides based on 2,4-dichlorophenoxyacetic acid (2,4-D) are widely used in these fields. The effects of 2,4-D on the germination of megaspores and initial sporophytic development of R. diphyllum were investigated. Six concentrations of 2,4-D (0.32; 0.64; 1.92; 4.80; 9.60 and 19.20 mg.L-1), and the control (0.00 mg.L-1), were tested in vitro, using Meyer's medium. Cultures were maintained in a growth chamber at 24 ± 1 ºC, under artificial light with nominal irradiance of 110 µmol.m-2/s and 16 hours photoperiod. Megaspore germination was lower at 9.60 and 19.20 mg.L-1 of 2,4-D (56 and 48 percent, respectively), compared with the control (68 percent). Herbicide concentrations of up to 1.92 mg.L-1 did not significantly decrease the number of sporophytes formed. At 19.20 mg.L-1, no sporophytes were formed. The lengths of the primary root, primary and secondary leaves were greater at concentrations of 0.32 and 0.64 mg.L-1 of 2,4-D. Low concentrations of 2,4-D do not affect germination rates and initial development of R. diphyllum in a significant way. However, higher concentrations (9.60 and 19.20 mg.L-1) affect substantially the germination of the megaspores and interfere with the establishment of the species.(AU)

Regnellidium diphyllum é considerada ameaçada, ocorrendo no Estado do Rio Grande do Sul, Brasil, e em algumas localidades vizinhas no Uruguai, na Argentina e no Estado de Santa Catarina. Cresce em áreas alagáveis, frequentemente alteradas para atividades agrícolas. Herbicidas baseados em ácido 2,4-diclorofenoxiacético (2,4-D) são largamente utilizados nestas plantações. Os efeitos do 2,4-D sobre a germinação de megásporos e o desenvolvimento esporofítico inicial de R. diphyllum foram investigados. Seis concentrações de 2,4-D (0,32; 0,64; 1,92; 4,80; 9,60 e 19,20 mg.L-1), além do controle (0,00 mg.L-1) foram testadas in vitro, utilizando meio de Meyer. As culturas foram mantidas em câmara de germinação a 24 ± 1 ºC, sob luz artificial, com irradiância nominal de 110 mmol.m-2/s e fotoperíodo de 16 horas. A germinação de megásporos foi menor em 9,60 e 19,20 mg.L-1 de 2,4-D (56 e 48 por cento, respectivamente), comparada ao controle (68 por cento). Concentrações até 1,92 mg.L-1 não diminuíram significativamente o número de esporófitos formados. Em 19,20 mg.L-1, não houve formação de esporófitos. Os comprimentos da raiz primária e das folhas primária e secundária foram maiores em concentrações de 0,32 e 0,64 mg.L-1 de 2,4-D. Baixas concentrações de 2,4-D não afetam significantemente as taxas de germinação e o desenvolvimento inicial de R. diphyllum. Entretanto, maiores concentrações (9,60 e 19,20 mg.L-1) afetam substancialmente a germinação de megásporos e interferem no estabelecimento da espécie.(AU)

The Pb-hyperaccumulator aquatic fern Salvinia minima Baker, responds to Pb(2+) by increasing phytochelatins via changes in SmPCS expression and in phytochelatin synthase activity.

The relationship between accumulation of Pb(2+) and the activation of chelation and metal sequestration mechanisms mediated by phytochelatins (PC) was analyzed in the Pb(2+) hyperaccumulator aquatic fern Salvinia minima, after exposure to 40microM Pb(NO(3))(2). The tissue accumulation pattern of lead and the phytochelatin biosynthesis responses were analyzed in both, S. minima submerged root-like modified fronds (here named "roots"), and in its aerial leaf-like fronds ("leaves"). S. minima roots accumulated a significantly higher concentrations of Pb(+2) than leaves did. Accumulation of Pb(2+) in roots was bi-phasic with a first uptake phase reached after 3h exposure and a second higher uptake phase reached after 24h exposure. In leaves, a single delayed, smaller uptake phase was attained only after 9h of exposure. In roots lead accumulation correlated with an increased phytochelatin synthase (PCS) activity and an enhanced PC production. A higher proportion of polymerized PC(4) was observed in both tissues of exposed S. minima plants relative to unexposed ones, although a higher concentration of PC(4) was found in roots than in leaves. PCS activity and Pb(2+) accumulation was also higher in roots than in leaves. The expression levels of the S. minima PCS gene (SmPCS), in response to Pb(2+) treatment, were also evaluated. In S. minima leaves, the accumulation of Pb(2+) correlated with a marked increase in expression of SmPCS, suggesting a transcriptional regulation in the PCS activation and PC accumulation in this S. minima tissue. However, in roots, the basal expression of SmPCS was down-regulated after Pb(2+) treatment. This fact did not correlate with the later but strong increase in both, PCS activity and PC production; suggesting that the PC biosynthesis activation in S. minima roots occurs only by post-translational activation of PCS. Taken together, our data suggest that the accumulation of PC in S. minima is a direct response to Pb(2+) accumulation, and phytochelatins do participate as one of the mechanism to cope with Pb(2+) of this Pb-hyperaccumulator aquatic fern.