Sean Graham


Research Interests

Plant phylogenetics
Plant systematics
Plant evolution
Land plants

Relevant Degree Programs



Master's students
Doctoral students
Postdoctoral Fellows
Any time / year round

Plant phylogenetics and molecular evolution

I am open to hosting Visiting International Research Students (non-degree, up to 12 months).

Complete these steps before you reach out to a faculty member!

Check requirements
  • Familiarize yourself with program requirements. You want to learn as much as possible from the information available to you before you reach out to a faculty member. Be sure to visit the graduate degree program listing and program-specific websites.
  • Check whether the program requires you to seek commitment from a supervisor prior to submitting an application. For some programs this is an essential step while others match successful applicants with faculty members within the first year of study. This is either indicated in the program profile under "Admission Information & Requirements" - "Prepare Application" - "Supervision" or on the program website.
Focus your search
  • Identify specific faculty members who are conducting research in your specific area of interest.
  • Establish that your research interests align with the faculty member’s research interests.
    • Read up on the faculty members in the program and the research being conducted in the department.
    • Familiarize yourself with their work, read their recent publications and past theses/dissertations that they supervised. Be certain that their research is indeed what you are hoping to study.
Make a good impression
  • Compose an error-free and grammatically correct email addressed to your specifically targeted faculty member, and remember to use their correct titles.
    • Do not send non-specific, mass emails to everyone in the department hoping for a match.
    • Address the faculty members by name. Your contact should be genuine rather than generic.
  • Include a brief outline of your academic background, why you are interested in working with the faculty member, and what experience you could bring to the department. The supervision enquiry form guides you with targeted questions. Ensure to craft compelling answers to these questions.
  • Highlight your achievements and why you are a top student. Faculty members receive dozens of requests from prospective students and you may have less than 30 seconds to pique someone’s interest.
  • Demonstrate that you are familiar with their research:
    • Convey the specific ways you are a good fit for the program.
    • Convey the specific ways the program/lab/faculty member is a good fit for the research you are interested in/already conducting.
  • Be enthusiastic, but don’t overdo it.
Attend an information session

G+PS regularly provides virtual sessions that focus on admission requirements and procedures and tips how to improve your application.


Graduate Student Supervision

Doctoral Student Supervision (Jan 2008 - May 2021)
Phylogenetics and evolution of monocot mycoheterotrophs and a newly demonstrated lineage of carnivorous monocots (2020)

No abstract available.

Phylogenomic studies of the monocot sister orders Pandanales and Dioscoreales (2020)

No abstract available.

Systematics of deep and recent lineages of bryophytes using phylogenomic approaches (2020)

No abstract available.

Phylogenomics and Comparative Plastome Analysis of Mycoheterotrophic Plants (2016)

Plastid genomes (plastomes) of fully mycoheterotrophic plants (which obtain nutrition from fungi and have lost photosynthesis) may exhibit accelerated substitution rates, gene losses and structural rearrangements compared to their more stable photosynthetic relatives. Distantly related lineages provide independent data points to study plastome degradation. I used Sanger sequencing to assess the utility of three nonphotosynthetic plastid genes in phylogenetic inference of seven monocot families that include mycoheterotrophic taxa. I also assembled full plastome genomes for multiple mycoheterotrophic monocots, a heterotrophic conifer (Parasitaxus, Podocarpaceae) and autotrophic relatives for comparative analysis. Phylogenomic inferences are robust to different likelihood approaches and often extensive gene loss, are generally congruent with the few-gene analyses, and are insensitive to long branches, in contrast to parsimony. Patterns of gene loss and retention are largely in agreement with hypothesized trajectories, starting with plastid NAD(P)H dehydrogenase, followed by the loss of other photosynthesis-related genes, and ending in gradual loss of transcriptional apparatus and other non-photosynthesis related genes. I observed retention (delayed loss?) of genes encoding subunits of plastid-encoded RNA polymerase (Parasitaxus and some species in Petrosavia, Petrosaviaceae), plastid ATP synthase (Petrosavia, perhaps Parasitaxus in modified form) and Rubisco (Petrosavia), consistent with secondary non-photosynthetic functions of the latter two complexes. Some group IIA introns appear to be retained despite the loss of the plastid intron maturase, matK. Retained open reading frames are generally under strong purifying selection in Sciaphila (Triuridaceae). Genome contraction is the major mode of genome rearrangement, with severe reduction seen in some lineages (e.g., Apteria in Burmanniaceae is reduced to ~16 kb). Some mycoheterotrophs are nearly or completely colinear with autotrophic lineages (Geosiris in Iridaceae, at ~123 kb). Others have multiple minor or major rearrangements, which may be unrelated to the presence or absence of an inverted repeat (IR). Four independent IR losses were observed (in Burmanniaceae, Corsiaceae, Petrosaviaceae and Triuridaceae), an extra IR copy evolved in Campylosiphon (Burmanniaceae), and an entire IR re-evolved in Parasitaxus. Shifts in IR boundaries were also found in all mycoheterotrophs. Within-taxon comparisons (e.g., in Corsiaceae and Petrosavia) also underline that idiosyncratic evolutionary changes may occur following each loss of photosynthesis.

View record

The phylogeny and evolution of two ancient lineages of aquatic plants (2013)

In my thesis I aim to improve our phylogenetic and evolutionary knowledge of two ancient and distantly related groups of aquatic flowering plants, Hydatellaceae and Alismatales. While the phylogeny of monocots has received fairly intense scrutiny for two decades, some parts of its diversification have been less frequently investigated. One such lineage is the order Alismatales, which defines one of the deepest splits in monocot evolution. Many families of Alismatales are aquatic or semi-aquatic, and they have been implicated in historical discussions of monocot origins. I evaluate inter-familial relationships in the order, considering a suite of 17 plastid genes for 31 Alismatales taxa for all 13 recognized families. This study improves on our understanding of, and confidence in, higher-order Alismatales relationships. I also uncovered convergent gene loss of plastid-encoded subunits for the NADH dehydrogenase complex. I then expand monocot coverage outside Alismatales by including unpublished and newly sequenced data for otherorders. This large-scale sample facilitated a re-evaluation of monocot phylogeny and molecular dating, the latter using 25 fossil constraints. Previously included in the monocot order Poales, Hydatellaceae are a small family of ephemeral aquatics relatively recently found to be the sister group of water lilies (Cabombaceae and Nymphaeaceae). I present the first molecularphylogeny of the family and evaluate aspects of the family's morphological evolution. I show how sexual system shifts are associated with shifts in otherreproductive traits. I also infer a temporal scale for Hydatellaceae diversification using a two-step Bayesian approach. I use the resulting dated tree to address biogeographic patterns and aspects of niche evolution. I show that its "Gondwanan" distribution is the result of long-distance dispersal and not continental rifting, and demonstrate strong phylogenetic niche conservatism in the family. These studies expand our understanding of evolution in Hydatellaceae, and provide a substantial update to our understanding of Alismatales (and more generally monocot) phylogeny and divergence times.

View record

Effects of adaxial-abaxial signalling on leaf polarity (2012)

The unifying theme of this thesis is adaxial-abaxial or dorsiventral patterning in leaves. The adaxial-abaxial axis sets the thickness of a leaf and without the appropriate juxtaposition of the adaxial and abaxial domains, radialized leaves develop. The underlying genetic mechanisms of the development of these polarity defects started to be elucidated only over the past 20 years in the model Arabidopsis, in particular. I investigated this patterning in a variety of non-model species. Firstly, I investigated the variability of dorsiventral polarity in plants with naturally occurring radialized leaves including Allium, Nepenthes, Krishna fig, Pelargonium, several Cactaceae species, and popREVOLUTA mutant of a poplar hybrid. Subsequent chapters aimed to incorporate morphology and anatomy with molecular genetics in order to elucidate the underlying basis of the phenotype of interest in species that have not been used as model systems for leaf development, including canola and poplar. A novel mutant (lamina epiphylla, lip) was identified in canola, which has adaxialized leaves and leaf-derived organs. Some of the HD-ZIPIII candidate genes were sequenced in canola, but I was unable to determine the location of the LIP mutation. The rest of this thesis focuses on the abaxial greening and unifacial petiole phenotypes seen in some species of poplar that have isobilateral leaves (others in the genus have bifacial leaves). YABBY, KANADI, and HD-ZIPIII genes are some of the major contributors to setting proper adaxial-abaxial polarity and I investigated the relationships of these genes by identifying the orthologs in Arabidopsis, poplar, and eucalyptus (a genus that shares the abaxial greening phenotype with poplar). Further, I studied the species relationships within the genus Populus in order to establish the ancestral state of leaf type. I determined that bifacial leaves are likely derived within the genus. Finally, two poplar species (black cottonwood with bifacial leaves and hybrid aspen with isobilateral leaves) were compared on the basis of morphology, anatomy, and molecular genetics in order to determine the underlying basis of the abaxial greening and unifacial petiole phenotypes in hybrid aspen. I identified a subset of genes that may be involved in determining these phenotypes, but further investigation is needed.

View record

Molecular phylogenetics of mosses and relatives (2011)

Substantial ambiguities still remain concerning the broad backbone of moss phylogeny. I surveyed 17 slowly evolving plastid genes from representative taxa to reconstruct phylogenetic relationships among the major lineages of mosses in the overall context of land-plant phylogeny. I first designed 78 bryophyte-specific primers and demonstrated that they permit straightforward amplification and sequencing of 14 core genes across a broad range of bryophytes (three of the 17 genes required more effort). In combination, these genes can generate sturdy and well-resolved phylogenetic inferences of higher-order moss phylogeny, with little evidence of conflict among different data partitions or analyses. Liverworts are strongly supported as the sister group of the remaining land plants, and hornworts as sister to vascular plants. Within mosses, besides confirming some previously published findings based on other markers, my results substantiallyimprove support for major branching patterns that were ambiguous before. The monogeneric classes Takakiopsida and Sphagnopsida likely represent the first and second split within moss phylogeny, respectively. However, this result is shown to be sensitive to the strategy used to estimate DNA substitution model parameter values and to different data partitioning methods. Regarding the placement of remaining nonperistomate lineages, the [[[Andreaeobryopsida, Andreaeopsida], Oedipodiopsida], peristomate mosses] arrangement receives moderate to strong support. Among peristomate mosses, relationships among Polytrichopsida, Tetraphidopsida andBryopsida remain unclear, as do the earliest splits within sublcass Bryidae. A Funariidae, [Timmiidae, [Dicranidae, Bryidae]]] arrangement is strongly supported, as are major relationships within subclasses Funariidae and Dicranidae. I also reconstructed the phylogeny of the nonperistomate moss family Andreaeaceae, with a focus on costate taxa, using two complementary sets of plastid markers and taxa. The major subgenera (Andreaea and Chasmocalyx) and sections of Andreaea (Andreaea and Nerviae) are rejected as monophyletic. Well-supported lineages include clades comprising: (1) Andreaea nivalis and A. rigida (northern hemisphere members of subgenus Chasmocalyx) and A. blyttii (section Nerviae); (2) most of the remainder of Nerviae; (3) a mixture of costate and ecostate species from Chasmocalyx, Nerviae,all sampled members of section Andreaea, and subgenus Acroschisma. Relationships among the major lineages, including the root of the family, are all well supported.

View record

Molecular Phylogenetic Studies of the Vascular Plants (2009)

No abstract available.

Master's Student Supervision (2010 - 2020)
DNA barcoding the vascular plant flora of southern British Columbia (2017)

DNA barcoding is a tool for rapidly identifying species based on short, standardized sequences of DNA, for example in situations where this may be difficult using morphology alone. I assembled a core DNA barcode reference library for southern British Columbia, home to ~54% of the vascular plant species in Canada, using the core plastid barcode loci rbcL and matK, and assessed its utility for identifying species in this region. The library comprises 4,812 sequences obtained from field-collected and herbarium tissue samples, supplemented with sequences downloaded from BOLD and GenBank, with at least one sequence for 75.4% of the vascular plant species occurring below 50°N in British Columbia. Sequence recoverability was significantly higher for rbcL than for matK (93.5% and 80.2%, respectively), and only marginally lower for both markers when using herbarium specimens (90.5% for rbcL and 77.8% for matK), which demonstrates the future feasibility of using museum specimens for completing a southern BC barcode reference library. As a proxy for assessing marker effectiveness, I scored resolution at the level of species and genus using tests of monophyly for Neighbour Joining trees, and performed sequence similarity searches with BLASTn analyses, both for each locus separately and for a dual-locus marker system (rbcL+matK; scored as a cumulative percentage in the BLAST analyses). Ignoring species represented by singleton sequences, the highest overall level of discrimination (66.9% of species and 91.6% of genera) was achieved for BLASTn analysis of rbcL+matK together. This work represents a significant contribution to a nation-wide barcode database, and provides a preliminary platform for ecological and other applications requiring species identification, where traditional methods are not feasible.

View record

Plastid genome evolution in partially and fully mycoheterotrophic eudicots (2016)

Plastid-genome evolution following photosynthesis loss is characterized by substantial change, contrasting with strong conservation in most photosynthetic land plants. Common features of reduced plastid genomes across diverse heterotrophic lineages point to a predictable trajectory of genome degradation, but this has been only partly tested. Here I document the molecular evolution of plastid genomes belonging to several mycoheterotroph lineages in Ericaceae, Gentianaceae and Polygalaceae, which include several independent origins of mycoheterotrophy in eudicot angiosperms that span different time scales since photosynthesis loss. I used next-generation and Sanger sequencing techniques to assemble complete plastomes or gene sets for comparative analyses of gene content and genome structure, and phylogenomic inference. I also sequenced several partially mycoheterotrophic and fully autotrophic relatives. Patterns of gene loss in mycoheterotroph plastomes are generally consistent with a previously hypothesized trajectory of change, starting with the loss of plastid NAD(P)H dehydrogenase before full loss of photosynthesis, and ending (here) with substantial reduction in genes involved in the translation apparatus and other nonphotosynthetic functions. Several retentions (delayed losses) of subunit genes for plastid-encoded polymerase, plastid ATP synthase and Rubisco are also consistent with hypothesized secondary (nonphotosynthetic) functions for these complexes. Two within-genus comparisons (for Epirixanthes in Polygalaceae and Voyria in Gentianaceae) demonstrate substantially different levels of genome degradation, consistent with heterogeneity in rates of genome change after a given origin of full mycoheterotrophy. Mycoheterotrophs in two families (Ericaceae, Polygalaceae) have extensive genome rearrangement compared to most land plants, contrasting with near colinearity in mycoheterotrophic members of Gentianaceae (despite sometimes extensive genome reduction in the latter). However, these contrasting patterns are apparently not associated with transitions to mycoheterotrophy, as photosynthetic relatives in Ericaceae and Polygalaceae are also substantially rearranged—or with inverted repeat loss (evident in Epirixanthes pallida, Polygalaceae), as autotrophic Polygala retains its inverted repeats. Phylogenomic inferences of core eudicot phylogeny made using the retained genes are generally well supported and robust to a variety of phylogenetic approaches, and are also congruent with recent phylogenetic studies in each mycoheterotrophic family.

View record

Phylogenetics and molecular evolution of Alismatales based on whole plastid genomes (2015)

The order Alismatales is a mostly aquatic group of monocots that displays substantial morphological and life history diversity, including the seagrasses, the only land plants that have re-colonized marine environments. Past phylogenetic studies of the order have either considered a single gene with dense taxonomic sampling, or several genes with thinner sampling. Despite substantial progress based on these studies, multiple phylogenetic uncertainties still remain concerning higher-order phylogenetic relationships. To address these issues, I completed a near-genus level sampling of the core alismatid families and the phylogenetically isolated family Tofieldiaceae, adding these new data to published sequences of Araceae and other monocots, eudicots and ANITA-grade angiosperms. I recovered whole plastid genomes (plastid gene sets representing up to 83 genes per taxa) and analyzed them using maximum likelihood and parsimony approaches. I recovered a well supported phylogenetic backbone for most of the order, with all families supported as monophyletic, and with strong support for most inter- and intrafamilial relationships. A major exception is the relative arrangement of Araceae, core alismatids and Tofieldiaceae; although most analyses recovered Tofieldiaceae as the sister-group of the rest of the order, this result was not well supported. Different partitioning schemes used in the likelihood analyses had little effect on patterns of clade support across the order, and the parsimony and likelihood results were generally highly congruent. I also used the inferred phylogeny of Alismatales to study the loss of the mostly plastid-encoded NADH dehydrogenase enzyme complex in the order. This enzyme is hypothesized to be involved in mitigating photooxidative stress by inducing chlororespiration. The inclusion or exclusion of ndh pseudogenes had little impact on the main phylogenetic results. Previous work hypothesized three independent losses/pseudogenization events within the core alismatids, which I confirmed here. I also inferred an additional loss in Tofieldiaceae, the first example in unsubmerged species of Alismatales. The repeated loss of plastid NADH dehydrogenase may spur future research into the physiological bases of the loss.

View record

Using a low-copy nuclear gene (phosphoglycerate kinase; PGK) to explore the phylogeny of the aquatic plant family Hydatellaceae (Nymphaeales) (2014)

Hydatellaceae are a small aquatic family of 12 species related to water lilies, part of the ANITA grade of angiosperms. Our current understanding of phylogenetic relationships in the family comes from several plastid genes and nuclear ITS data. These data sets are generally highly congruent, and lend support to the monophyly of multiple species. However, the published nuclear ITS tree was unrooted (outgroups were too distant to align), and there were several minor phylogenetic conflicts between plastid and ITS gene trees for three closely related species, Trithuria bibracteata, T. occidentalis, and T. submersa; two of these species were also not reciprocally monophyletic in individual gene trees. The position of T. occidentalis was also based on very limited plastid data, and there was no molecular evidence to link staminate and pistillate individuals in this species. To further clarify phylogenetic relationships and species boundaries, I recovered two copies of nuclear-encoded phosphoglycerate kinase (PGK) gene from taxa in Hydatellaceae and several water lilies. I reconstructed the history of the PGK duplication in angiosperms as a whole. I also added plastid data from additional populations of several species, and estimated the dated species tree using a Bayesian multispecies coalescent approach to reconcile different gene trees. The angiosperm-level PGK gene tree indicated that the duplication of PGK gene may have happened around the origin of angiosperms. The root of Hydatellaceae implied by concatenated nuclear PGK matches that inferred from plastid data. Trithuria occidentalis is clearly placed in sect. Trithuria, and staminate and pistillate individuals of this species are linked together using new evidence from the plastid and PGK genes. Phylogenetic relationships inferred using each PGK copy are consistent with the sectional relationships inferred using plastid and ITS data, with less sharply defined species boundaries. I also explore the possibility here that some of the incongruence that I observed between individual genes trees and in inferred species trees is a consequence of additional minor gene duplications or polyploidization/introgression events.

View record



If this is your researcher profile you can log in to the Faculty & Staff portal to update your details and provide recruitment preferences.


Discover the amazing research that is being conducted at UBC!