Syntrichia ruralis: emerging model moss genome reveals a conserved and previously unknown regulator of desiccation in flowering plants

Water scarcity, resulting from climate change, poses a significant threat to ecosystems. Syntrichia ruralis, a dryland desiccation-tolerant moss, provides valuable insights into survival of water-limited conditions. We sequenced the genome of S. ruralis, conducted transcriptomic analyses, and performed comparative genomic and transcriptomic analyses with existing genomes and transcriptomes, including with the close relative S. caninervis. We took a genetic approach to characterize the role of an S. ruralis transcription factor, identified in transcriptomic analyses, in Arabidopsis thaliana. The genome was assembled into 12 chromosomes encompassing 21 169 protein-coding genes. Comparative analysis revealed copy number and transcript abundance differences in known desiccation-associated gene families, and highlighted genome-level variation among species that may reflect adaptation to different habitats. A significant number of abscisic acid (ABA)-responsive genes were found to be negatively regulated by a MYB transcription factor (MYB55) that was upstream of the S. ruralis ortholog of ABA-insensitive 3 (ABI3). We determined that this conserved MYB transcription factor, uncharacterized in Arabidopsis, acts as a negative regulator of an ABA-dependent stress response in Arabidopsis. The new genomic resources from this emerging model moss offer novel insights into how plants regulate their responses to water deprivation.

The dynamics of external water conduction in the dryland moss Syntrichia

Syntrichia relies on external water conduction for photosynthesis, survival, and reproduction, a condition referred to as ectohydry. Capillarity spaces are abundant in Syntrichia, but the link between function and morphology is complex. The aim of this study was to provide a better understanding of species-specific morphological traits underlying the functions of water conduction and storage. We used an environmental scanning electron microscope and confocal microscopy for observing anatomical characters in the leaves of Syntrichia species. We also measured hydration/dehydration curves to understand the rate of conduction and dehydration by experimental approaches. Syntrichia is an ectohydric moss that can externally transport and store water from the base of the stem using capillary action. We propose a new framework to study ectohydric capabilities, which incorporates three morphological scales and the timing of going from completely dehydrated to fully hydrated. Characters of interest in this model include cell anatomy (papillae development, hyaline basal cells and laminar cells), architecture of the stem (concavity and orientation) and whole clump characteristics (density of stems). We report significant variations in the speed of conduction, water holding capacity and hydration associated with each species studied (11 in total). All Syntrichia species are capable of external water conduction and storage, but the relevant traits differ among species. These results help to understand potential evolutionary and ecological trade-offs among speed of water conduction, water holding capacity, ontogeny, and differing habitat requirements. An integrative view of ectohydry in Syntrichia contributes to understanding the water relationships of mosses.

The phylogeny of Syntrichia: An ecologically diverse clade of mosses with an origin in South America

PREMISE

To address the biodiversity crisis, we need to understand the evolution of all organisms and how they fill geographic and ecological space. Syntrichia is one of the most diverse and dominant genera of mosses, ranging from alpine habitats to desert biocrusts, yet its evolutionary history remains unclear.

METHODS

We present a comprehensive phylogenetic analysis of Syntrichia, based on both molecular and morphological data, with most of the named species and closest outgroups represented. In addition, we provide ancestral-state reconstructions of water-related traits and a global biogeographic analysis.

RESULTS

We found 10 major well-resolved subclades of Syntrichia that possess geographical or morphological coherence, in some cases representing previously accepted genera. We infer that the extant species diversity of Syntrichia likely originated in South America in the early Eocene (56.5-43.8 million years ago [Mya]), subsequently expanded its distribution to the neotropics, and finally dispersed to the northern hemisphere. There, the clade experienced a recent diversification (15-12 Mya) into a broad set of ecological niches (e.g., the S. caninervis and S. ruralis complexes). The transition from terricolous to either saxicolous or epiphytic habitats occurred more than once and was associated with changes in water-related traits.

CONCLUSIONS

Our study provides a framework for understanding the evolutionary history of Syntrichia through the combination of morphological and molecular characters, revealing that migration events that shaped the current distribution of the clade have implications for morphological character evolution in relation to niche diversity.

Responses of Biocrust and Associated Soil Bacteria to Novel Climates Are Not Tightly Coupled

Climate change is expanding drylands even as land use practices degrade them. Representing ∼40% of Earth’s terrestrial surface, drylands rely on biological soil crusts (biocrusts) for key ecosystem functions including soil stability, biogeochemical cycling, and water capture. Understanding how biocrusts adapt to climate change is critical to understanding how dryland ecosystems will function with altered climate. We investigated the sensitivity of biocrusts to experimentally imposed novel climates to track changes in productivity and stability under both warming and cooling scenarios. We established three common gardens along an elevational-climate gradient on the Colorado Plateau. Mature biocrusts were collected from each site and reciprocally transplanted intact. Over 20 months we monitored visible species composition and cover, chlorophyll a, and the composition of soil bacterial communities using high throughput sequencing. We hypothesized that biocrusts replanted at their home site would show local preference, and biocrusts transplanted to novel environments would maintain higher cover and stability at elevations higher than their origin, compared to at elevations lower than their origin. We expected responses of the visible biocrust cover and soil bacterial components of the biocrust community to be coupled, with later successional taxa showing higher sensitivity to novel environments. Only high elevation sourced biocrusts maintained higher biocrust cover and community stability at their site of origin. Biocrusts from all sources had higher cover and stability in the high elevation garden. Later successional taxa decreased cover in low elevation gardens, suggesting successional reversal with warming. Visible community composition was influenced by both source and transplant environment. In contrast, soil bacterial community composition was not influenced by transplant environments but retained fidelity to the source. Thus, responses of the visible and soil bacterial components of the biocrust community were not coupled. Synthesis: Our results suggest biocrust communities are sensitive to climate change, and loss of species and function can be expected, while associated soil bacteria may be buffered against rapid change.

Spatial phylogenetics of Japanese ferns: Patterns, processes, and implications for conservation

PREMISE

Biodiversity is often only measured with species richness; however, this metric ignores evolutionary history and is not sufficient for making conservation decisions. Here, we characterize multiple facets and drivers of biodiversity to understand how these relate to bioregions and conservation status in the ferns of Japan.

METHODS

We compiled a community data set of 1239 grid cells (20 × 20 km each) including 672 taxa based on >300,000 specimen records. We combined the community data with a phylogeny and functional traits to analyze taxonomic, phylogenetic, and functional diversity and modeled biodiversity metrics in response to environmental factors and reproductive mode. Hierarchical clustering was used to delimit bioregions. Conservation status and threats were assessed by comparing the overlap of significantly diverse grid cells with conservation zones and range maps of native Japanese deer.

RESULTS

Taxonomic richness was highest at mid-latitudes. Phylogenetic and functional diversity and phylogenetic endemism were highest in small southern islands. Relative phylogenetic and functional diversity were high at high and low latitudes, and low at mid-latitudes. Grid cells were grouped into three (phylogenetic) or four (taxonomic) major bioregions. Temperature and apomixis were identified as drivers of biodiversity patterns. Conservation status was generally high for grid cells with significantly high biodiversity, but the threat due to herbivory by deer was greater for taxonomic richness than other metrics.

CONCLUSIONS

Our integrative approach reveals previously undetected patterns and drivers of biodiversity in the ferns of Japan. Future conservation efforts should recognize that threats can vary by biodiversity metric and consider multiple metrics when establishing conservation priorities.

Genotypic confirmation of a biased phenotypic sex ratio in a dryland moss using restriction fragment length polymorphisms

PREMISE

In dioicous mosses, sex is determined by a single U (female, ♀) or V (male, ♂) chromosome. Although a 1 : 1 sex ratio is expected following meiosis, phenotypic sex ratios based on the production of gametangia are often female-biased. The dryland moss Syntrichia caninervis (Pottiaceae) is notable for its low frequency of sex expression and strong phenotypic female bias. Here we present a technique to determine genotypic sex in a single shoot of S. caninervis, and report results of a case study examining genotypic and phenotypic sex ratios.

METHODS

We reanalyzed 271 non-expressing gametophyte shoots from a previous study on S. caninervis sex expression across microhabitats using a restriction fragment length polymorphism (RFLP) method.

RESULTS

We recovered a genotypic sex ratio in non-expressing shoots of 18.4♀ : 1♂, which exceeds the female bias of the phenotypic ratio (5.3♀ : 1♂; P = 0.013). We also found that the distribution of male and female genotypes across microsites with different levels of sun exposure was not predicted by patterns of sex expression in these microsites.

DISCUSSION

These findings contribute to our understanding of how the environment may modulate sex ratios in S. caninervis, either through its direct influence on sex expression or through selection on genotypes with particular sex expression phenotypes.

Spatial phylogenetics of the Chinese angiosperm flora provides insights into endemism and conservation

The flora of China is well known for its high diversity and endemism. Identifying centers of endemism and designating conservation priorities are essential goals for biodiversity studies. However, there is no comprehensive study from a rigorous phylogenetic perspective to understand patterns of diversity and endemism and to guide biodiversity conservation in China. We conducted a spatial phylogenetic analysis of the Chinese angiosperm flora at the generic level to identify centers of neo- and paleo-endemism. Our results indicate that: (i) the majority of grid cells in China with significantly high phylogenetic endemism (PE) were located in the mountainous regions; (ii) four of the nine centers of endemism recognized, located in northern and western China, were recognized for the first time; (iii) arid and semiarid regions in Northwest China were commonly linked to significant PE, consistent with other spatial phylogenetic studies worldwide; and (iv) six high-priority conservation gaps were detected by overlaying the boundaries of China's nature reserves on all significant PE cells. Overall, we conclude that the mountains of southern and northern China contain both paleo-endemics (ancient relictual lineages) and neo-endemics (recently diverged lineages). The areas we highlight as conservation priorities are important for broad-scale planning, especially in the context of evolutionary history preservation.

Transcriptomic Effects of Acute Ultraviolet Radiation Exposure on Two Syntrichia Mosses

Ultraviolet radiation (UVR) is a major environmental stressor for terrestrial plants. Here we investigated genetic responses to acute broadband UVR exposure in the highly desiccation-tolerant mosses Syntrichia caninervis and Syntrichia ruralis, using a comparative transcriptomics approach. We explored whether UVR protection is physiologically plastic and induced by UVR exposure, addressing the following questions: (1) What is the timeline of changes in the transcriptome with acute UVR exposure in these two species? (2) What genes are involved in the UVR response? and (3) How do the two species differ in their transcriptomic response to UVR? There were remarkable differences between the two species after 10 and 30 min of UVR exposure, including no overlap in significantly differentially abundant transcripts (DATs) after 10 min of UVR exposure and more than twice as many DATs for S. caninervis as there were for S. ruralis. Photosynthesis-related transcripts were involved in the response of S. ruralis to UVR, while membrane-related transcripts were indicated in the response of S. caninervis. In both species, transcripts involved in oxidative stress and those important for desiccation tolerance (such as late embryogenesis abundant genes and early light-inducible protein genes) were involved in response to UVR, suggesting possible roles in UVR tolerance and cross-talk with desiccation tolerance in these species. The results of this study suggest potential UVR-induced responses that may have roles outside of UVR tolerance, and that the response to URV is different in these two species, perhaps a reflection of adaptation to different environmental conditions.

Natural ultraviolet radiation exposure alters photosynthetic biology and improves recovery from desiccation in a desert moss

Plants in dryland ecosystems experience extreme daily and seasonal fluctuations in light, temperature, and water availability. We used an in situ field experiment to uncover the effects of natural and reduced levels of ultraviolet radiation (UV) on maximum PSII quantum efficiency (Fv/Fm), relative abundance of photosynthetic pigments and antioxidants, and the transcriptome in the desiccation-tolerant desert moss Syntrichia caninervis. We tested the hypotheses that: (i) S. caninervis plants undergo sustained thermal quenching of light [non-photochemical quenching (NPQ)] while desiccated and after rehydration; (ii) a reduction of UV will result in improved recovery of Fv/Fm; but (iii) 1 year of UV removal will de-harden plants and increase vulnerability to UV damage, indicated by a reduction in Fv/Fm. All field-collected plants had extremely low Fv/Fm after initial rehydration but recovered over 8 d in lab-simulated winter conditions. UV-filtered plants had lower Fv/Fm during recovery, higher concentrations of photoprotective pigments and antioxidants such as zeaxanthin and tocopherols, and lower concentrations of neoxanthin and Chl b than plants exposed to near natural UV levels. Field-grown S. caninervis underwent sustained NPQ that took days to relax and for efficient photosynthesis to resume. Reduction of solar UV radiation adversely affected recovery of Fv/Fm following rehydration.

Spatial phylogenetics of butterflies in relation to environmental drivers and angiosperm diversity across North America

Broad-scale, quantitative assessments of insect biodiversity and the factors shaping it remain particularly poorly explored. Here we undertook a spatial phylogenetic analysis of North American butterflies to test whether climate stability and temperature gradients have shaped their diversity and endemism. We also performed the first quantitative comparisons of spatial phylogenetic patterns between butterflies and flowering plants. We expected concordance between the two groups based on shared historical environmental drivers and presumed strong butterfly-host plant specializations. We instead found that biodiversity patterns in butterflies are strikingly different from flowering plants, especially warm deserts. In particular, butterflies show different patterns of phylogenetic clustering compared with flowering plants, suggesting differences in habitat conservation between the two groups. These results suggest that shared biogeographic histories and trophic associations do not necessarily assure similar diversity outcomes. The work has applied value in conservation planning, documenting warm deserts as a North American butterfly biodiversity hotspot.

To dry perchance to live: Insights from the genome of the desiccation-tolerant biocrust moss Syntrichia caninervis

With global climate change, water scarcity threatens whole agro/ecosystems. The desert moss Syntrichia caninervis, an extremophile, offers novel insights into surviving desiccation and heat. The sequenced S. caninervis genome consists of 13 chromosomes containing 16 545 protein-coding genes and 2666 unplaced scaffolds. Syntenic relationships within the S. caninervis and Physcomitrella patens genomes indicate the S. caninervis genome has undergone a single whole genome duplication event (compared to two for P. patens) and evidence suggests chromosomal or segmental losses in the evolutionary history of S. caninervis. The genome contains a large sex chromosome composed primarily of repetitive sequences with a large number of Copia and Gypsy elements. Orthogroup analyses revealed an expansion of ELIP genes encoding proteins important in photoprotection. The transcriptomic response to desiccation identified four structural clusters of novel genes. The genomic resources established for this extremophile offer new perspectives for understanding the evolution of desiccation tolerance in plants.

Strategies of desiccation tolerance vary across life phases in the moss Syntrichia caninervis

PREMISE

Desiccation tolerance (DT) is a widespread phenomenon among land plants, and variable ecological strategies for DT are likely to exist. Using Syntrichia caninervis, a dryland moss and model system used in DT studies, we hypothesized that DT is lowest in juvenile (protonemal) tissues, highest in asexual reproductive propagules (gemmae), and intermediate in adults (shoots). We tested the long-standing hypothesis of an inherent constitutive strategy of DT in this species.

METHODS

Plants were rapidly dried to levels of equilibrating relative humidity (RHeq) ranging from 0 to 93%. Postrehydration recovery was assessed using chlorophyll fluorescence, regeneration rates, and visual tissue damage. For each life phase, we estimated the minimum rate of drying (RoD_min ) at RHeq = 42% that did not elicit damage 24 h postrehydration.

RESULTS

DT strategy varied with life phase, with adult shoots having the lowest RoD_min (10-25 min), followed by gemmae (3-10 h) and protonema (14-20 h). Adult shoots exhibited no detectable damage 24 h postrehydration following a rapid-dry only at the highest RHeq used (93%), but when dried to lower RHs the response declined to <50% of control fluorescence values. Notably, immediately following rehydration (0 h postrehydration), shoots were damaged below control levels of fluorescence regardless of the RHeq, thus implicating damage.

CONCLUSIONS

Life phases of the moss S. caninervis had a range of strategies from near constitutive (adult shoots) to demonstrably inducible (protonema). A new response variable for assessing degree of DT is introduced as the minimum rate of drying from which full recovery occurs.

Vascular plant extinction in the continental United States and Canada

Extinction rates are expected to increase during the Anthropocene. Current extinction rates of plants and many animals remain unknown. We quantified extinctions among the vascular flora of the continental United States and Canada since European settlement. We compiled data on apparently extinct species by querying plant conservation databases, searching the literature, and vetting the resulting list with botanical experts. Because taxonomic opinion varies widely, we developed an index of taxonomic uncertainty (ITU). The ITU ranges from A to F, with A indicating unanimous taxonomic recognition and F indicating taxonomic recognition by only a single author. The ITU allowed us to rigorously evaluate extinction rates. Our data suggest that 51 species and 14 infraspecific taxa, representing 33 families and 49 genera of vascular plants, have become extinct in our study area since European settlement. Seven of these taxa exist in cultivation but are extinct in the wild. Most extinctions occurred in the west, but this outcome may reflect the timing of botanical exploration relative to settlement. Sixty-four percent of extinct plants were single-site endemics, and many occurred outside recognized biodiversity hotspots. Given the paucity of plant surveys in many areas, particularly prior to European settlement, the actual extinction rate of vascular plants is undoubtedly much higher than indicated here.

Morphological and physiological traits in relation to carbon balance in a diverse clade of dryland mosses

Plant functional trait analyses have focused almost exclusively on vascular plants, but bryophytes comprise ancient and diverse plant lineages that have widespread global distributions and important ecological functions in terrestrial ecosystems. We examined a diverse clade of dryland mosses, Syntrichia, and studied carbon balance during a precipitation event (C-balance), a functional trait related to physiological functioning, desiccation tolerance, survival, and ecosystem carbon and nitrogen cycling. We examined variability in C-balance among 14 genotypes of Syntrichia and measured an additional 10 physiological and 13 morphological traits at the cell, leaf, shoot, and clump level. C-balance varied 20-fold among genotypes, and highest C-balances were associated with long, narrow leaves with awns, and small cells with thick cell walls, traits that may influence water uptake and retention during a precipitation event. Ordination analyses revealed that the axis most strongly correlated with C-balance included the maximum chlorophyll fluorescence, F_m , indicating the importance of photosystem II health for C exchange. C-balance represents a key functional trait in bryophytes, but its measurement is time intensive and not feasible to measure on large scales. We propose two models (using physiological and morphological traits) to predict C-balance, whereby identifying simpler to measure traits for trait databases.

Facets of phylodiversity: evolutionary diversification, divergence and survival as conservation targets

Biodiversity is often described as having multiple facets, including species richness, functional diversity and phylogenetic diversity. In this paper, we argue that phylogenetic diversity itself has three distinct facets-lineage diversification, character divergence and survival time-that can be quantified using distinct branch length metrics on an evolutionary tree. Each dimension is related to different processes of macroevolution, has different spatial patterns and is tied to distinct goals for conserving biodiversity and protecting its future resilience and evolutionary potential. We compared the landscapes identified as top conservation priorities by each of these three metrics in a conservation gap analysis for California, a world biodiversity hotspot, using herbarium data on the biogeography and evolutionary relationships of more than 5000 native plant species. Our analysis incorporated a novel continuous metric of current land protection status, fine-scale data on landscape intactness and an optimization algorithm used to identify complementary priority sites containing concentrations of taxa that are evolutionarily unique, vulnerable due to small range size and/or poorly protected across their ranges. Top conservation priorities included pockets of coastal and northern California that ranked highly for all three phylodiversity dimensions and for species richness, as well as sites uniquely identified by each metric whose value may depend on whether properties such as genetic divergence, high net diversification or independent survival experience are most desirable in an Anthropocene flora.This article is part of the theme issue 'Biological collections for understanding biodiversity in the Anthropocene'.

Digitization protocol for scoring reproductive phenology from herbarium specimens of seed plants

PREMISE OF THE STUDY

Herbarium specimens provide a robust record of historical plant phenology (the timing of seasonal events such as flowering or fruiting). However, the difficulty of aggregating phenological data from specimens arises from a lack of standardized scoring methods and definitions for phenological states across the collections community.

METHODS AND RESULTS

To address this problem, we report on a consensus reached by an iDigBio working group of curators, researchers, and data standards experts regarding an efficient scoring protocol and a data-sharing protocol for reproductive traits available from herbarium specimens of seed plants. The phenological data sets generated can be shared via Darwin Core Archives using the Extended MeasurementOrFact extension.

CONCLUSIONS

Our hope is that curators and others interested in collecting phenological trait data from specimens will use the recommendations presented here in current and future scoring efforts. New tools for scoring specimens are reviewed.

Spatial phylogenetics of the native California flora

BACKGROUND

California is a world floristic biodiversity hotspot where the terms neo- and paleo-endemism were first applied. Using spatial phylogenetics, it is now possible to evaluate biodiversity from an evolutionary standpoint, including discovering significant areas of neo- and paleo-endemism, by combining spatial information from museum collections and DNA-based phylogenies. Here we used a distributional dataset of 1.39 million herbarium specimens, a phylogeny of 1083 operational taxonomic units (OTUs) and 9 genes, and a spatial randomization test to identify regions of significant phylogenetic diversity, relative phylogenetic diversity, and phylogenetic endemism (PE), as well as to conduct a categorical analysis of neo- and paleo-endemism (CANAPE).

RESULTS

We found (1) extensive phylogenetic clustering in the South Coast Ranges, southern Great Valley, and deserts of California; (2) significant concentrations of short branches in the Mojave and Great Basin Deserts and the South Coast Ranges and long branches in the northern Great Valley, Sierra Nevada foothills, and the northwestern and southwestern parts of the state; (3) significant concentrations of paleo-endemism in Northwestern California, the northern Great Valley, and western Sonoran Desert, and neo-endemism in the White-Inyo Range, northern Mojave Desert, and southern Channel Islands. Multiple analyses were run to observe the effects on significance patterns of using different phylogenetic tree topologies (uncalibrated trees versus time-calibrated ultrametric trees) and using different representations of OTU ranges (herbarium specimen locations versus species distribution models).

CONCLUSIONS

These analyses showed that examining the geographic distributions of branch lengths in a statistical framework adds a new dimension to California floristics that, in comparison with climatic data, helps to illuminate causes of endemism. In particular, the concentration of significant PE in more arid regions of California extends previous ideas about aridity as an evolutionary stimulus. The patterns seen are largely robust to phylogenetic uncertainty and time calibration but are sensitive to the use of occurrence data versus modeled ranges, indicating that special attention toward improving geographic distributional data should be top priority in the future for advancing understanding of spatial patterns of biodiversity.

Species richness and endemism in the native flora of California

PREMISE OF THE STUDY

California's vascular flora is the most diverse and threatened in temperate North America. Previous studies of spatial patterns of Californian plant diversity have been limited by traditional metrics, non-uniform geographic units, and distributional data derived from floristic descriptions for only a subset of species.

METHODS

We revisited patterns of sampling intensity, species richness, and relative endemism in California based on equal-area spatial units, the full vascular flora, and specimen-based distributional data. We estimated richness, weighted endemism (inverse range-weighting of species), and corrected weighted endemism (weighted endemism corrected for richness), and performed a randomization test for significantly high endemism.

KEY RESULTS

Possible biases in herbarium data do not obscure patterns of high richness and endemism at the spatial resolution studied. High species richness was sometimes associated with significantly high endemism (e.g., Klamath Ranges) but often not. In Stebbins and Major's (1965) main endemism hotspot, Southwestern California, species richness is high across much of the Peninsular and Transverse ranges but significantly high endemism is mostly localized to the Santa Rosa and San Bernardino mountains. In contrast, species richness is low in the Channel Islands, where endemism is significantly high, as also found for much of the Death Valley region.

CONCLUSIONS

Measures of taxonomic richness, even with greater weighting of range-restricted taxa, are insufficient for identifying areas of significantly high endemism that warrant conservation attention. Differences between our findings and those in previous studies appear to mostly reflect the source and scale of distributional data, and recent analytical refinements.

Assessing biodiversity and endemism using phylogenetic methods across multiple taxonomic groups

Identifying geographical areas with the greatest representation of the tree of life is an important goal for the management and conservation of biodiversity. While there are methods available for using a single phylogenetic tree to assess spatial patterns of biodiversity, there has been limited exploration of how separate phylogenies from multiple taxonomic groups can be used jointly to map diversity and endemism. Here, we demonstrate how to apply different phylogenetic approaches to assess biodiversity across multiple taxonomic groups. We map spatial patterns of phylogenetic diversity/endemism to identify concordant areas with the greatest representation of biodiversity across multiple taxa and demonstrate the approach by applying it to the Murray–Darling basin region of southeastern Australia. The areas with significant centers of phylogenetic diversity and endemism were distributed differently for the five taxonomic groups studied (plant genera, fish, tree frogs, acacias, and eucalypts); no strong shared patterns across all five groups emerged. However, congruence was apparent between some groups in some parts of the basin. The northern region of the basin emerges from the analysis as a priority area for future conservation initiatives focused on eucalypts and tree frogs. The southern region is particularly important for conservation of the evolutionary heritage of plants and fishes.

Continental scale patterns and predictors of fern richness and phylogenetic diversity

Because ferns have a wide range of habitat preferences and are widely distributed, they are an ideal group for understanding how diversity is distributed. Here we examine fern diversity on a broad-scale using standard and corrected richness measures as well as phylogenetic indices; in addition we determine the environmental predictors of each diversity metric. Using the combined records of Australian herbaria, a dataset of over 60,000 records was obtained for 89 genera to infer richness. A molecular phylogeny of all the genera was constructed and combined with the herbarium records to obtain phylogenetic diversity patterns. A hotspot of both taxic and phylogenetic diversity occurs in the Wet Tropics of northeastern Australia. Although considerable diversity is distributed along the eastern coast, some important regions of diversity are identified only after sample-standardization of richness and through the phylogenetic metric. Of all of the metrics, annual precipitation was identified as the most explanatory variable, in part, in agreement with global and regional fern studies. However, precipitation was combined with a different variable for each different metric. For corrected richness, precipitation was combined with temperature seasonality, while correlation of phylogenetic diversity to precipitation plus radiation indicated support for the species-energy hypothesis. Significantly high and significantly low phylogenetic diversity were found in geographically separate areas. These separate areas correlated with different climatic conditions such as seasonality in precipitation. The phylogenetic metrics identified additional areas of significant diversity, some of which have not been revealed using traditional taxonomic analyses, suggesting that different ecological and evolutionary processes have operated over the continent. Our study demonstrates that it is possible and vital to incorporate evolutionary metrics when inferring biodiversity hotspots from large compilations of data.

Evolution of oleosin in land plants

Oleosins form a steric barrier surface on lipid droplets in cytoplasm, preventing them from contacting and coalescing with adjacent droplets. Oleosin genes have been detected in numerous plant species. However, the presence of oleosin genes in the most basally diverging lineage of land plants, liverworts, has not been reported previously. Thus we explored whether liverworts have an oleosin gene. In Marchantia polymorpha L., a thalloid liverwort, one predicted sequence was found that could encode oleosin, possessing the hallmark of oleosin, a proline knot (-PX5SPX3P-) motif. The phylogeny of the oleosin gene family in land plants was reconstructed based on both nucleotide and amino acid sequences of oleosins, from 31 representative species covering almost all the main lineages of land plants. Based on our phylogenetic trees, oleosin genes were classified into three groups: M-oleosins (defined here as a novel group distinct from the two previously known groups), low molecular weight isoform (L-oleosin), and high molecular weight isoform (H-oleosin), according to their amino-acid organization, phylogenetic relationships, expression tissues, and immunological characteristics. In liverworts, mosses, lycophytes, and gymnosperms, only M-oleosins have been described. In angiosperms, however, while this isoform remains and is highly expressed in the gametophyte pollen tube, two other isoforms also occur, L-oleosins and H-oleosins. Phylogenetic analyses suggest that the M-oleosin isoform is the precursor to the ancestor of L-oleosins and H-oleosins. The later two isoforms evolved by successive gene duplications in ancestral angiosperms. At the genomic level, most oleosins possess no introns. If introns are present, in both the L-isoform and the M-isoform a single intron inserts behind the central region, while in the H-isoform, a single intron is located at the 5'-terminus. This study fills a major gap in understanding functional gene evolution of oleosin in land plants, shedding new light on evolutionary transitions of lipid storage strategies.

Organellar genomes of the four-toothed moss, Tetraphis pellucida

BACKGROUND

Mosses are the largest of the three extant clades of gametophyte-dominant land plants and remain poorly studied using comparative genomic methods. Major monophyletic moss lineages are characterised by different types of a spore dehiscence apparatus called the peristome, and the most important unsolved problem in higher-level moss systematics is the branching order of these peristomate clades. Organellar genome sequencing offers the potential to resolve this issue through the provision of both genomic structural characters and a greatly increased quantity of nucleotide substitution characters, as well as to elucidate organellar evolution in mosses. We publish and describe the chloroplast and mitochondrial genomes of Tetraphis pellucida, representative of the most phylogenetically intractable and morphologically isolated peristomate lineage.

RESULTS

Assembly of reads from Illumina SBS and Pacific Biosciences RS sequencing reveals that the Tetraphis chloroplast genome comprises 127,489 bp and the mitochondrial genome 107,730 bp. Although genomic structures are similar to those of the small number of other known moss organellar genomes, the chloroplast lacks the petN gene (in common with Tortula ruralis) and the mitochondrion has only a non-functional pseudogenised remnant of nad7 (uniquely amongst known moss chondromes).

CONCLUSIONS

Structural genomic features exist with the potential to be informative for phylogenetic relationships amongst the peristomate moss lineages, and thus organellar genome sequences are urgently required for exemplars from other clades. The unique genomic and morphological features of Tetraphis confirm its importance for resolving one of the major questions in land plant phylogeny and for understanding the evolution of the peristome, a likely key innovation underlying the diversity of mosses. The functional loss of nad7 from the chondrome is now shown to have occurred independently in all three bryophyte clades as well as in the early-diverging tracheophyte Huperzia squarrosa.

Quantifying phytogeographical regions of Australia using geospatial turnover in species composition

The largest digitized dataset of land plant distributions in Australia assembled to date (750,741 georeferenced herbarium records; 6,043 species) was used to partition the Australian continent into phytogeographical regions. We used a set of six widely distributed vascular plant groups and three non-vascular plant groups which together occur in a variety of landscapes/habitats across Australia. Phytogeographical regions were identified using quantitative analyses of species turnover, the rate of change in species composition between sites, calculated as Simpson's beta. We propose six major phytogeographical regions for Australia: Northern, Northern Desert, Eremaean, Eastern Queensland, Euronotian and South-Western. Our new phytogeographical regions show a spatial agreement of 65% with respect to previously defined phytogeographical regions of Australia. We also confirm that these new regions are in general agreement with the biomes of Australia and other contemporary biogeographical classifications. To assess the meaningfulness of the proposed phytogeographical regions, we evaluated how they relate to broad scale environmental gradients. Physiographic factors such as geology do not have a strong correspondence with our proposed regions. Instead, we identified climate as the main environmental driver. The use of an unprecedentedly large dataset of multiple plant groups, coupled with an explicit quantitative analysis, makes this study novel and allows an improved historical bioregionalization scheme for Australian plants. Our analyses show that: (1) there is considerable overlap between our results and older biogeographic classifications; (2) phytogeographical regions based on species turnover can be a powerful tool to further partition the landscape into meaningful units; (3) further studies using phylogenetic turnover metrics are needed to test the taxonomic areas.

Desiccation tolerance in bryophytes: a reflection of the primitive strategy for plant survival in dehydrating habitats?

Bryophytes are a non-monophyletic group of three major lineages (liverworts, hornworts, and mosses) that descend from the earliest branching events in the phylogeny of land plants. We postulate that desiccation tolerance is a primitive trait, thus mechanisms by which the first land plants achieved tolerance may be reflected in how extant desiccation-tolerant bryophytes survive drying. Evidence is consistent with extant bryophytes employing a tolerance strategy of constitutive cellular protection coupled with induction of a recovery/repair mechanism upon rehydration. Cellular structures appear intact in the desiccated state but are disrupted by rapid uptake of water upon rehydration, but cellular integrity is rapidly regained. The photosynthetic machinery appears to be protected such that photosynthetic activity recovers quickly. Gene expression responds following rehydration and not during drying. Gene expression is translationally controlled and results in the synthesis of a number of proteins, collectively called rehydrins. Some prominent rehydrins are similar to Late Embryogenesis Abundant (LEA) proteins, classically ascribed a protection function during desiccation. The role of LEA proteins in a rehydrating system is unknown but data indicates a function in stabilization and reconstitution of membranes. Phylogenetic studies using a Tortula ruralis LEA-like rehydrin led to a re-examination of the evolution of desiccation tolerance. A new phylogenetic analysis suggests that: (i) the basic mechanisms of tolerance seen in modern day bryophytes have changed little from the earliest manifestations of desiccation tolerance in land plants, and (ii) vegetative desiccation tolerance in the early land plants may have evolved from a mechanism present first in spores.

Complete plastome sequences of Equisetum arvense and Isoetes flaccida: implications for phylogeny and plastid genome evolution of early land plant lineages

Background

Despite considerable progress in our understanding of land plant phylogeny, several nodes in the green tree of life remain poorly resolved. Furthermore, the bulk of currently available data come from only a subset of major land plant clades. Here we examine early land plant evolution using complete plastome sequences including two previously unexamined and phylogenetically critical lineages. To better understand the evolution of land plants and their plastomes, we examined aligned nucleotide sequences, indels, gene and nucleotide composition, inversions, and gene order at the boundaries of the inverted repeats.

Results

We present the plastome sequences of Equisetum arvense, a horsetail, and of Isoetes flaccida, a heterosporous lycophyte. Phylogenetic analysis of aligned nucleotides from 49 plastome genes from 43 taxa supported monophyly for the following clades: embryophytes (land plants), lycophytes, monilophytes (leptosporangiate ferns + Angiopteris evecta + Psilotum nudum + Equisetum arvense), and seed plants. Resolution among the four monilophyte lineages remained moderate, although nucleotide analyses suggested that P. nudum and E. arvense form a clade sister to A. evecta + leptosporangiate ferns. Results from phylogenetic analyses of nucleotides were consistent with the distribution of plastome gene rearrangements and with analysis of sequence gaps resulting from insertions and deletions (indels). We found one new indel and an inversion of a block of genes that unites the monilophytes.

Conclusions

Monophyly of monilophytes has been disputed on the basis of morphological and fossil evidence. In the context of a broad sampling of land plant data we find several new pieces of evidence for monilophyte monophyly. Results from this study demonstrate resolution among the four monilophytes lineages, albeit with moderate support; we posit a clade consisting of Equisetaceae and Psilotaceae that is sister to the "true ferns," including Marattiaceae.

Chloroplast genome sequence of the moss Tortula ruralis: gene content, polymorphism, and structural arrangement relative to other green plant chloroplast genomes

Background

Tortula ruralis, a widely distributed species in the moss family Pottiaceae, is increasingly used as a model organism for the study of desiccation tolerance and mechanisms of cellular repair. In this paper, we present the chloroplast genome sequence of T. ruralis, only the second published chloroplast genome for a moss, and the first for a vegetatively desiccation-tolerant plant.

Results

The Tortula chloroplast genome is ~123,500 bp, and differs in a number of ways from that of Physcomitrella patens, the first published moss chloroplast genome. For example, Tortula lacks the ~71 kb inversion found in the large single copy region of the Physcomitrella genome and other members of the Funariales. Also, the Tortula chloroplast genome lacks petN, a gene found in all known land plant plastid genomes. In addition, an unusual case of nucleotide polymorphism was discovered.

Conclusions

Although the chloroplast genome of Tortula ruralis differs from that of the only other sequenced moss, Physcomitrella patens, we have yet to determine the biological significance of the differences. The polymorphisms we have uncovered in the sequencing of the genome offer a rare possibility (for mosses) of the generation of DNA markers for fine-level phylogenetic studies, or to investigate individual variation within populations.

Generational differences in response to desiccation stress in the desert moss Tortula inermis

BACKGROUND AND AIMS

Active growth in post-embryonic sporophytes of desert mosses is restricted to the cooler, wetter months. However, most desert mosses have perennial gametophytes. It is hypothesized that these life history patterns are due, in part, to a reduced desiccation tolerance for sporophytes relative to gametophytes.

METHODS

Gametophytes with attached post-embryonic sporophytes of Tortula inermis (early seta elongation phenophase) were exposed to two levels of desiccation stress, one rapid-dry cycle and two rapid-dry cycles, then moistened and allowed to recover, resume development, and/or regenerate for 35 d in a growth chamber.

KEY RESULTS

Gametophytes tolerated the desiccation treatments well, with 93 % survival through regenerated shoot buds and/or protonemata. At the high stress treatment, a significantly higher frequency of burned leaves and browned shoots occurred. Sporophytes were far more sensitive to desiccation stress, with only 23 % surviving after the low desiccation stress treatment, and 3 % surviving after the high desiccation stress treatment. While the timing of protonemal production and sporophytic phenophases was relatively unaffected by desiccation stress, shoots exposed to one rapid-dry cycle produced shoots more rapidly than shoots exposed to two rapid-dry cycles.

CONCLUSIONS

It is concluded that sporophytes of Tortula inermis are more sensitive to rapid drying than are maternal gametophytes, and that sporophyte abortion in response to desiccation results from either reduced desiccation tolerance of sporophytes relative to gametophytes, or from a termination of the sporophyte on the part of the gametophyte in response to stress.

Taking the first steps towards a standard for reporting on phylogenies: Minimum Information About a Phylogenetic Analysis (MIAPA)

In the eight years since phylogenomics was introduced as the intersection of genomics and phylogenetics, the field has provided fundamental insights into gene function, genome history and organismal relationships. The utility of phylogenomics is growing with the increase in the number and diversity of taxa for which whole genome and large transcriptome sequence sets are being generated. We assert that the synergy between genomic and phylogenetic perspectives in comparative biology would be enhanced by the development and refinement of minimal reporting standards for phylogenetic analyses. Encouraged by the development of the Minimum Information About a Microarray Experiment (MIAME) standard, we propose a similar roadmap for the development of a Minimal Information About a Phylogenetic Analysis (MIAPA) standard. Key in the successful development and implementation of such a standard will be broad participation by developers of phylogenetic analysis software, phylogenetic database developers, practitioners of phylogenomics, and journal editors.

The first complete chloroplast genome sequence of a lycophyte, Huperzia lucidula (Lycopodiaceae)

We used a unique combination of techniques to sequence the first complete chloroplast genome of a lycophyte, Huperzia lucidula. This plant belongs to a significant clade hypothesized to represent the sister group to all other vascular plants. We used fluorescence-activated cell sorting (FACS) to isolate the organelles, rolling circle amplification (RCA) to amplify the genome, and shotgun sequencing to 8x depth coverage to obtain the complete chloroplast genome sequence. The genome is 154,373 bp, containing inverted repeats of 15,314 bp each, a large single-copy region of 104,088 bp, and a small single-copy region of 19,657 bp. Gene order is more similar to those of mosses, liverworts, and hornworts than to gene order for other vascular plants. For example, the Huperzia chloroplast genome possesses the bryophyte gene order for a previously characterized 30 kb inversion, thus supporting the hypothesis that lycophytes are sister to all other extant vascular plants. The lycophyte chloroplast genome data also enable a better reconstruction of the basal tracheophyte genome, which is useful for inferring relationships among bryophyte lineages. Several unique characters are observed in Huperzia, such as movement of the gene ndhF from the small single copy region into the inverted repeat. We present several analyses of evolutionary relationships among land plants by using nucleotide data, inferred amino acid sequences, and by comparing gene arrangements from chloroplast genomes. The results, while still tentative pending the large number of chloroplast genomes from other key lineages that are soon to be sequenced, are intriguing in themselves, and contribute to a growing comparative database of genomic and morphological data across the green plants.

The cost of realized sexual reproduction: assessing patterns of reproductive allocation and sporophyte abortion in a desert moss

The desert moss Syntrichia caninervis exhibits one of the most skewed sex ratios in the plant kingdom, with female individuals far outnumbering male individuals (exceeding 14♀:1♂). The "cost of sex hypothesis" derives from allocational theory and predicts that the sex which is most expensive should be the rarer sex. This hypothesis, which, as considered here represents the realized cost of sexual reproduction, is contingent upon two assumptions that are explored: (1) that male sex expression is more expensive than female sex expression, and (2) that sexual reproduction is resource limited. Using inflorescence biomass and discounting sperm, male sex expression was found to be in the neighborhood of one order of magnitude more expensive than female sex expression, and this difference is reflected in higher numbers of gametangia per male inflorescence, presence of paraphyses in male inflorescences, and a much longer developmental time for male inflorescences. The realized cost of female reproduction from two communities dominated by S. caninervis was found to be lower than the realized cost of male sexual reproduction. Resource-limited reproduction was assessed by determining the frequency of sporophyte abortion, the age distribution of sporophyte abortions, and patterns of sporophyte abortion that may be density dependent. Among ten sexually reproducing populations, abortive sporophytes occurred at a frequency of 0.64. Abortive sporophytes averaged 8% the mass of mature sporophytes, and cohort sporophytes from the same individual female were found to abort in a density-dependent pattern. We conclude that the two assumptions, upon which the cost of sex hypothesis depends, are supported.

Sex expression, skewed sex ratios, and microhabitat distribution in the dioecious desert moss Syntrichia caninervis (Pottiaceae)

The moss Syntrichia caninervis is the dominant soil bryophyte in a blackbrush (Coleogyne ramosissima) community in the southern Nevada Mojave Desert, with a mean cover of 6.3%. A survey of the 10-ha study site revealed an expressed ramet sex ratio of 14♀ : 1♂ (N = 890), with 85% of ramets not expressing sex over their life span, and an expressed population sex ratio of 40♀ : 2♂ : 1♀♂ (female : male : mixed-sex, N = 89), with 52% of populations not expressing sex. A greater incidence of sex expression was associated with shaded microsites, higher soil moisture content, and taller ramets. Shaded microsites had higher surface soil moisture levels than exposed microsites. In the exposed microhabitat, surface soil moisture was positively correlated with ramet height but not with sex expression. Male ramets and populations were restricted to shaded microhabitats, whereas female ramets and populations were found in both shaded and exposed microhabitats, suggesting gender specialization. The rarity of mature sporophytes, found in 0% of the ramets sampled and in only 3% of the populations, is probably due to the rarity of mixed-sex populations. We hypothesize that mixed-sex populations are rare because of factors relating to male rarity and that the differential cost of sex expression reduces the clonal growth capacity of male individuals.

Phylogenetic relationships of the liverworts (Hepaticae), a basal embryophyte lineage, inferred from nucleotide sequence data of the chloroplast gene rbcL

Sequence data from the chloroplast-encoded gene rbcL were obtained for 24 liverworts, a basal group of embryophytes. Maximum likelihood and parsimony analyses of these data, along with data from other major green plant lineages, confirm hypotheses based on morphological data, such as the paraphyly of bryophytes, and the basal position of liverworts. Molecular data corroborate the deep separation between the complex thalloid and leafy/simple thalloid liverworts implied by morphological data, but the monophyly of liverworts could not be rejected. The effects of accounting for site-to-site rate heterogeneity in these data were examined using maximum likelihood methods. Comparison of trees obtained with and without rate heterogeneity showed that simply allowing for heterogeneity had a greater improvement on likelihood score than optimization of transition/transversion bias. Incorporation of site-to-site rate heterogeneity in the larger analysis, however, did not necessarily change which topology was favored. Properties of rbcL sequences from the two liverwort groups were compared. Significantly different substitution rates were found between leafy/simple thalloid and complex thalloid liverwort taxa, with rates of rbcL sequence evolution in leafy/simple thalloid taxa being higher and more indicative of those of vascular plants, and with those of complex thalloid taxa (such as Marchantia) being slower. Codon usage in rbcL in complex thalloid liverworts was biased toward NNU and NNA, compared to the leafy/simple thalloid liverworts. Although base composition and relative substitution rates differed between the two groups, no significant differences were detected within each of the two groups of liverworts. The signal present in first and second codon sites versus third codon sites was compared. While the third codon positions in rbcL across this taxon sampling are highly variable (with only 15 constant sites of 439), the trees obtained were in general agreement with trees from the entire data set and with trees obtained from independent sources of data. The presence of signal in third codon positions across greater than 400 MY of plant evolution means that definitions of saturation based on pair-wise comparisons of sequences inadequately assess phylogenetic signal.

Cladistic analysis of molecular and morphological data

Considerable progress has been made recently in phylogenetic reconstruction in a number of groups of organisms. This progress coincides with two major advances in systematics: new sources have been found for potentially informative characters (i.e., molecular data) and (more importantly) new approaches have been developed for extracting historical information from old or new characters (i.e., Hennigian phylogenetic systematics or cladistics). The basic assumptions of cladistics (the existence and splitting of lineages marked by discrete, heritable, and independent characters, transformation of which occurs at a rate slower than divergence of lineages) are discussed and defended. Molecular characters are potentially greater in quantity than (and usually independent of) more traditional morphological characters, yet their great simplicity (i.e., fewer potential character states; problems with determining homology), and difficulty of sufficient sampling (particularly from fossils) can lead to special difficulties. Expectations of the phylogenetic behavior of different types of data are investigated from a theoretical standpoint, based primarily on variation in the central parameter lambda (branch length in terms of expected number of character changes per segment of a tree), which also leads to possibilities for character and character state weighting. Also considered are prospects for representing diverse yet clearly monophyletic clades in larger-scale cladistic analyses, e.g., the exemplar method vs. "compartmentalization" (a new approach involving substituting an inferred "archetype" for a large clade accepted as monophyletic based on previous analyses). It is concluded that parsimony is to be preferred for synthetic, "total evidence" analyses because it appears to be a robust method, is applicable to all types of data, and has an explicit and interpretable evolutionary basis.