Data from: Phylogeny of Megalastrum based on plastid nucleotide DNA sequences
Data files
Jun 04, 2026 version files 1.86 MB
-
Concat_2025_export.nex
1.85 MB
-
IQ_tree_output_(1).contree.tre
11.69 KB
-
README.md
3.23 KB
Abstract
This dataset is an alignment of seven plastid markers. These were generated from 156 specimens of Megalastrum, and represent 72 of the 93 species that are currently recognized. The alignment also includes 46 species from 14 other Dryopteridaceae genera. The alignment was generated using MAFFT. Maximium likelihood phylogenies were were generated using IQ-tree v2.1.3 and MrBayes v.3.2.6. In agreement with previous findings, results indicate that Megalastrum is monophyletic and sister to Rumohra. We recognize six main clades: 1) The Mascarene Clade, with two species endemic to the Mascarenes and having echinate-fenestrate spores; 2) the Pulverulentum Clade, consisting of ten species typically with large, 3–4-pinnate leaves and strongly denticulate scales often darkened apically; 3) the Cristate-spored Clade, harboring 57 species with cristate spores; 4) the Villosum Clade, with a single species endemic to the Greater Antilles (M. villosum, the type of the genus); 5) the Circumaustral Clade, consisting of seven species occurring at high (33°–45°) southern latitudes; and 6) the Southeastern Brazil Clade, containing 16 species mostly or entirely restricted to the region.
Dataset DOI: 10.5061/dryad.8931zcs5t
Description of the data and file structure
Files:
1. Concat_2025_export.nex
Final concatenated nucleotide alignment in NEXUS format.
2. IQ_tree_output_(1).contree.tre
Consensus phylogenetic tree inferred from the concatenated alignment.
Description:
This dataset contains the final concatenated molecular sequence alignment and corresponding phylogenetic tree used in the Megalastrum study. The alignment contains 202 taxa and 9103 aligned nucleotide positions.
Notes:
- Taxon names use underscores instead of spaces.
- Missing data are represented by question marks (?).
- Alignment gaps are represented by hyphens (-).
- The alignment and tree terminal labels were curated to correspond directly.
Contact:
Michael Sundue
Royal Botanic Garden Edinburgh
Date:
24 May 2026
Files and variables
File: Concat_2025_export.nex
Description: FILES AND VARIABLES
File 1:
Concat_2025_export.nex
Description:
Concatenated nucleotide sequence alignment in NEXUS format used for phylogenetic analyses of Megalastrum and related Dryopteridaceae.
Variables / structure:
- Rows correspond to taxa (terminal labels).
- Columns correspond to aligned nucleotide positions.
- Character states are standard DNA nucleotide codes:
A, C, G, T
- Missing data are represented by question marks (?).
- Alignment gaps are represented by hyphens (-).
- Taxon names use underscores (_) instead of spaces.
- Terminal labels generally consist of:
Genus_species_collectionCode
Example:
Megalastrum_marginatum_RG232
Units:
Not applicable.
File 2:
IQ_tree_output (1).contree.tre
Description:
Consensus phylogenetic tree inferred from the concatenated alignment.
Variables / structure:
- Tree format: Newick
- Branch lengths are proportional to inferred sequence divergence.
- Terminal labels correspond directly to taxon names in the alignment.
Abbreviations:
- RG, ES, etc. = collector or DNA sample codes associated with voucher specimens.
Missing values:
- Missing nucleotide data are represented by "?" in the alignment.
- Alignment gaps are represented by "-".
Code/software
CODE / SOFTWARE
The data files can be viewed and analyzed using standard phylogenetic software that supports NEXUS and Newick formats.
Software used:
- IQ-TREE v1 or v2 (phylogenetic inference)
- FigTree v1.4 or later (tree visualization)
- Geneious (sequence assembly, alignment inspection, and matrix export)
- Mesquite (alignment and tree inspection)
- R (optional downstream analyses and visualization)
Common R packages used for phylogenetic analyses include:
- ape
- phytools
- ggtree
File relationships:
- Concat_2025_export.nex was used as the input alignment for phylogenetic analyses.
- IQ_tree_output (1).contree.tre is the resulting consensus phylogenetic tree inferred from the alignment.
File formats:
- .nex = NEXUS formatted nucleotide alignment
- .tre = Newick formatted phylogenetic tree
No custom scripts or proprietary code are included in this submission.
Taxon sampling—We included sequences from 156 Megalastrum specimens . These represent 72 of the 93 species (77%) recognized in the genus. As outgroups, we used 46 species from 14 other Dryopteridaceae genera (Appendix 1), including Lastreopsis, Parapolystichum, and Rumohra, which are most closely related to Megalastrum and classified together in the Lastreopsidoideae (Zuo et al. 2025). For the ingroup and outgroup samples sequenced and included in this study, voucher information and GenBank accession numbers are provided in Appendix 1 of the associated manuscript.
For five specimens, the species names used here differ from those provided in our monographs, where they were incorrectly cited. The first is Megalastrum rupicola (Kessler et al. 11621) from Bolivia, which was incorrectly cited by Moran et al. (2014b) as “M. praetermissum”; the second is M. organense (Labiak et al. 4371) from Brazil, which was incorrectly cited by Moran et al. (2009a) as “M. umbrinum”; the last three are M. oppositum from the Mascarenes, which were incorrectly cited by Rouhan & Moran (2011) as “M. canacae” (Duan & Zhang 2014). The correct names are used in Fig. 2 and Appendix 1.
DNA extraction, amplification, and sequencing—DNA was isolated from either herbarium specimens or leaves collected in the field and dried in silica-gel (Appendix 1). To avoid contamination, we collected young leaves free of epiphyllous bryophytes and obvious fungal or insect parasites. For silica-gel dried samples, total genomic DNA was extracted from approximately 1 cm2 of leaf tissue using a modified CTAB protocol based on Doyle & Doyle (1987; see Struwe et al. 1998). For herbarium material, the DNEasy Plant Mini kit Qiagen (Valencia, CA) was used following the manufacturer's protocol, but with a proteinase K digestion during the lysis step: 30 μL of proteinase K was added per tube, and tubes were incubated on a tipping plate at 42 °C for 24h.
The polymerase chain reaction (PCR) was used to amplify seven plastid regions: two coding genes (atpA and rbcL) and five non-coding intergenic spacers (psbA-trnH, atpB-rbcL, trnG-trnR, trnL-trnF, and rps4-trnS). PCR reactions were typically carried out in a 25 μL volume containing 1x PCR buffer, 2.5 mM MgCl2, 250 μM of each dNTP, 1 M betaine, 0.4 μM of each primer, 0.75 U Taq polymerase (Taq CORE kit; MP Biomedicals, Illkirch, France), and 0.75 μL of template DNA. Primer sequences and thermal cycling conditions are reported in Table 1.
The resulting PCR products were visualized on a 1 % agarose gel with ethidium bromide and purified using QIAquick PCR purification kit (Qiagen). They were directly sequenced with a Perkin Elmer ABI 377XL automated sequencer using the amplification primers and additional internal primers for some markers (Table 1). Forward and reverse sequences obtained were edited and assembled using Sequencher 4.9 (Gene Codes Corporation, Ann Arbor, MI, USA). The newly obtained consensus sequences were submitted to GenBank, and the alignments to Dryad (Moran et al. 2026); the accession numbers for both repositories are listed in Appendix 1.
Sequence alignment, indel coding, and phylogenetic analyses—We assembled a data matrix of two gene regions (atpA, rbcL) and five intergenic spacers (psbA-trnH, trnG-trnR, trnL-trnF, psbA-trnH, and rps4-trnS). Sequences for each region were aligned using the MAFFT (Katoh et al. 2002) plugin in Geneious 2024.07. Maximum likelihood (ML) phylogenetic analyses were performed using IQ-tree v2.1.3 (Nguyen et al. 2015; Minh et al. 2020). We implemented ModelFinder in IQ-tree (Kalyaanamoorthy et al. 2017) to select the best-fitting model of sequence evolution via Bayesian Information Criterion (BIC). Node support values were generated from 1000 ultrafast bootstrap (UFBoot) replicates (Hoang et al. 2018; Minh et al. 2013). Bayesian inference (BI) of phylogeny was performed using MrBayes v.3.2.6 (Ronquist & Huelsenbeck 2003) and implemented in the CIPRES portal (Miller et al. 2010), with four chains run for 6 million generations, sampling every 1000 generations. The resulting log files were inspected for convergence and adequate sampling using Tracer v.1.7.2 (Rambaut et al. 2018). The first 25% of trees was discarded as burn-in, and a majority-rule consensus tree was generated from the remaining trees. Support values shown are BI posterior probabilities (PP) and ML bootstrap (BS) percentages. The aligned data matrix is available in the Dryad Digital Repository (Moran et al. 2026). In the tree presented (Fig. 2), outgroups were omitted to save space and make the species names easier to read.
