Thermal optima in the hypoxia tolerance of marine ectotherms: Physiological causes and biogeographic consequences.

The minimum O2 needed to fuel the demand of aquatic animals is commonly observed to increase with temperature, driven by accelerating metabolism. However, recent measurements of critical O2 thresholds ("Pcrit") reveal more complex patterns, including those with a minimum at an intermediate thermal "optimum". To discern the prevalence, physiological drivers, and biogeographic manifestations of such curves, we analyze new experimental and biogeographic data using a general dynamic model of aquatic water breathers. The model simulates the transfer of oxygen from ambient water through a boundary layer and into animal tissues driven by temperature-dependent rates of metabolism, diffusive gas exchange, and ventilatory and circulatory systems with O2-protein binding. We find that a thermal optimum in Pcrit can arise even when all physiological rates increase steadily with temperature. This occurs when O2 supply at low temperatures is limited by a process that is more temperature sensitive than metabolism, but becomes limited by a less sensitive process at warmer temperatures. Analysis of published species respiratory traits suggests that this scenario is not uncommon in marine biota, with ventilation and circulation limiting supply under cold conditions and diffusion limiting supply at high temperatures. Using occurrence data, we show that species with these physiological traits inhabit lowest O2 waters near the optimal temperature for hypoxia tolerance and are restricted to higher O2 at temperatures above and below this optimum. Our results imply that hypoxia tolerance can decline under both cold and warm conditions and thus may influence both poleward and equatorward species range limits.

Constraining the oxygen requirements for modern microbial eukaryote diversity.

Eukaryotes originated prior to the establishment of modern marine oxygen (O2) levels. According to the body fossil and lipid biomarker records, modern (crown) microbial eukaryote lineages began diversifying in the ocean no later than ~800 Ma. While it has long been predicted that increasing atmospheric O2 levels facilitated the early diversification of microbial eukaryotes, the O2 levels needed to permit this diversification remain unconstrained. Using time-resolved geochemical parameter and gene sequence information from a model marine oxygen minimum zone spanning a range of dissolved O2 levels and redox states, we show that microbial eukaryote taxonomic richness and phylogenetic diversity remain the same until O2 declines to around 2 to 3% of present atmospheric levels, below which these diversity metrics become significantly reduced. Our observations suggest that increasing O2 would have only directly promoted early crown-eukaryote diversity if atmospheric O2 was below 2 to 3% of modern levels when crown-eukaryotes originated and then later met or surpassed this range as crown-eukaryotes diversified. If atmospheric O2 was already consistently at or above 2 to 3% of modern levels by the time that crown-eukaryotes originated, then the subsequent diversification of modern microbial eukaryotes was not directly driven by atmospheric oxygenation.

Decreasing Phanerozoic extinction intensity as a consequence of Earth surface oxygenation and metazoan ecophysiology.

The decline in background extinction rates of marine animals through geologic time is an established but unexplained feature of the Phanerozoic fossil record. There is also growing consensus that the ocean and atmosphere did not become oxygenated to near-modern levels until the mid-Paleozoic, coinciding with the onset of generally lower extinction rates. Physiological theory provides us with a possible causal link between these two observations-predicting that the synergistic impacts of oxygen and temperature on aerobic respiration would have made marine animals more vulnerable to ocean warming events during periods of limited surface oxygenation. Here, we evaluate the hypothesis that changes in surface oxygenation exerted a first-order control on extinction rates through the Phanerozoic using a combined Earth system and ecophysiological modeling approach. We find that although continental configuration, the efficiency of the biological carbon pump in the ocean, and initial climate state all impact the magnitude of modeled biodiversity loss across simulated warming events, atmospheric oxygen is the dominant predictor of extinction vulnerability, with metabolic habitat viability and global ecophysiotype extinction exhibiting inflection points around 40% of present atmospheric oxygen. Given this is the broad upper limit for estimates of early Paleozoic oxygen levels, our results are consistent with the relative frequency of high-magnitude extinction events (particularly those not included in the canonical big five mass extinctions) early in the Phanerozoic being a direct consequence of limited early Paleozoic oxygenation and temperature-dependent hypoxia responses.

Ediacaran reorganization of the marine phosphorus cycle.

The Ediacaran Period (635 to 541 Ma) marks the global transition to a more productive biosphere, evidenced by increased availability of food and oxidants, the appearance of macroscopic animals, significant populations of eukaryotic phytoplankton, and the onset of massive phosphorite deposition. We propose this entire suite of changes results from an increase in the size of the deep-water marine phosphorus reservoir, associated with rising sulfate concentrations and increased remineralization of organic P by sulfate-reducing bacteria. Simple mass balance calculations, constrained by modern anoxic basins, suggest that deep-water phosphate concentrations may have increased by an order of magnitude without any increase in the rate of P input from the continents. Strikingly, despite a major shift in phosphorite deposition, a new compilation of the phosphorus content of Neoproterozoic and early Paleozoic shows little secular change in median values, supporting the view that changes in remineralization and not erosional P fluxes were the principal drivers of observed shifts in phosphorite accumulation. The trigger for these changes may have been transient Neoproterozoic weathering events whose biogeochemical consequences were sustained by a set of positive feedbacks, mediated by the oxygen and sulfur cycles, that led to permanent state change in biogeochemical cycling, primary production, and biological diversity by the end of the Ediacaran Period.

Calibrating the coevolution of Ediacaran life and environment.

The rise of animals occurred during an interval of Earth history that witnessed dynamic marine redox conditions, potentially rapid plate motions, and uniquely large perturbations to global biogeochemical cycles. The largest of these perturbations, the Shuram carbon isotope excursion, has been invoked as a driving mechanism for Ediacaran environmental change, possibly linked with evolutionary innovation or extinction. However, there are a number of controversies surrounding the Shuram, including its timing, duration, and role in the concomitant biological and biogeochemical upheavals. Here we present radioisotopic dates bracketing the Shuram on two separate paleocontinents; our results are consistent with a global and synchronous event between 574.0 ± 4.7 and 567.3 ± 3.0 Ma. These dates support the interpretation that the Shuram is a primary and synchronous event postdating the Gaskiers glaciation. In addition, our Re-Os ages suggest that the appearance of Ediacaran macrofossils in northwestern Canada is identical, within uncertainty, to similar macrofossils from the Conception Group of Newfoundland, highlighting the coeval appearance of macroscopic metazoans across two paleocontinents. Our temporal framework for the terminal Proterozoic is a critical step for testing hypotheses related to extreme carbon isotope excursions and their role in the evolution of complex life.

Marine sponges in the once and future ocean.

Animals originated under hypoxia-low-oxygen conditions that are currently expanding throughout the global ocean. How marine sponges respond to hypoxia is both relevant to reconstructing early animal evolution and for forecasting the fate of modern marine ecosystems. In a new effort, multiple sponge species from two different oceans were found to tolerate hypoxia in the lab, revealing a more general capacity for hypoxia tolerance across sponges with implications for both the deep past and near future of animal life.

Statistical analysis of iron geochemical data suggests limited late Proterozoic oxygenation.

Sedimentary rocks deposited across the Proterozoic-Phanerozoic transition record extreme climate fluctuations, a potential rise in atmospheric oxygen or re-organization of the seafloor redox landscape, and the initial diversification of animals. It is widely assumed that the inferred redox change facilitated the observed trends in biodiversity. Establishing this palaeoenvironmental context, however, requires that changes in marine redox structure be tracked by means of geochemical proxies and translated into estimates of atmospheric oxygen. Iron-based proxies are among the most effective tools for tracking the redox chemistry of ancient oceans. These proxies are inherently local, but have global implications when analysed collectively and statistically. Here we analyse about 4,700 iron-speciation measurements from shales 2,300 to 360 million years old. Our statistical analyses suggest that subsurface water masses in mid-Proterozoic oceans were predominantly anoxic and ferruginous (depleted in dissolved oxygen and iron-bearing), but with a tendency towards euxinia (sulfide-bearing) that is not observed in the Neoproterozoic era. Analyses further indicate that early animals did not experience appreciable benthic sulfide stress. Finally, unlike proxies based on redox-sensitive trace-metal abundances, iron geochemical data do not show a statistically significant change in oxygen content through the Ediacaran and Cambrian periods, sharply constraining the magnitude of the end-Proterozoic oxygen increase. Indeed, this re-analysis of trace-metal data is consistent with oxygenation continuing well into the Palaeozoic era. Therefore, if changing redox conditions facilitated animal diversification, it did so through a limited rise in oxygen past critical functional and ecological thresholds, as is seen in modern oxygen minimum zone benthic animal communities.

Oxygenated Mesoproterozoic lake revealed through magnetic mineralogy.

Terrestrial environments have been suggested as an oxic haven for eukaryotic life and diversification during portions of the Proterozoic Eon when the ocean was dominantly anoxic. However, iron speciation and Fe/Al data from the ca. 1.1-billion-year-old Nonesuch Formation, deposited in a large lake and bearing a diverse assemblage of early eukaryotes, are interpreted to indicate persistently anoxic conditions. To shed light on these distinct hypotheses, we analyzed two drill cores spanning the transgression into the lake and its subsequent shallowing. While the proportion of highly reactive to total iron (FeHR/FeT) is consistent through the sediments and typically in the range taken to be equivocal between anoxic and oxic conditions, magnetic experiments and petrographic data reveal that iron exists in three distinct mineral assemblages resulting from an oxycline. In the deepest waters, reductive dissolution of iron oxides records an anoxic environment. However, the remainder of the sedimentary succession has iron oxide assemblages indicative of an oxygenated environment. At intermediate water depths, a mixed-phase facies with hematite and magnetite indicates low oxygen conditions. In the shallowest waters of the lake, nearly every iron oxide has been oxidized to its most oxidized form, hematite. Combining magnetics and textural analyses results in a more nuanced understanding of ambiguous geochemical signals and indicates that for much of its temporal duration, and throughout much of its water column, there was oxygen in the waters of Paleolake Nonesuch.

Oxygen, temperature and the deep-marine stenothermal cradle of Ediacaran evolution.

Ediacaran fossils document the early evolution of complex megascopic life, contemporaneous with geochemical evidence for widespread marine anoxia. These data suggest early animals experienced frequent hypoxia. Research has thus focused on the concentration of molecular oxygen (O2) required by early animals, while also considering the impacts of climate. One model, the Cold Cradle hypothesis, proposed the Ediacaran biota originated in cold, shallow-water environments owing to increased O2 solubility. First, we demonstrate using principles of gas exchange that temperature does have a critical role in governing the bioavailability of O2-but in cooler water the supply of O2 is actually lower. Second, the fossil record suggests the Ediacara biota initially occur approximately 571 Ma in deep-water facies, before appearing in shelf environments approximately 555 Ma. We propose an ecophysiological underpinning for this pattern. By combining oceanographic data with new respirometry experiments we show that in the shallow mixed layer where seasonal temperatures fluctuate widely, thermal and partial pressure ( pO2) effects are highly synergistic. The result is that temperature change away from species-specific optima impairs tolerance to low pO2. We hypothesize that deep and particularly stenothermal (narrow temperature range) environments in the Ediacaran ocean were a physiological refuge from the synergistic effects of temperature and low pO2.

Biodiversity response to natural gradients of multiple stressors on continental margins.

Sharp increases in atmospheric CO2 are resulting in ocean warming, acidification and deoxygenation that threaten marine organisms on continental margins and their ecological functions and resulting ecosystem services. The relative influence of these stressors on biodiversity remains unclear, as well as the threshold levels for change and when secondary stressors become important. One strategy to interpret adaptation potential and predict future faunal change is to examine ecological shifts along natural gradients in the modern ocean. Here, we assess the explanatory power of temperature, oxygen and the carbonate system for macrofaunal diversity and evenness along continental upwelling margins using variance partitioning techniques. Oxygen levels have the strongest explanatory capacity for variation in species diversity. Sharp drops in diversity are seen as O2 levels decline through the 0.5-0.15 ml l(-1) (approx. 22-6 µM; approx. 21-5 matm) range, and as temperature increases through the 7-10°C range. pCO2 is the best explanatory variable in the Arabian Sea, but explains little of the variance in diversity in the eastern Pacific Ocean. By contrast, very little variation in evenness is explained by these three global change variables. The identification of sharp thresholds in ecological response are used here to predict areas of the seafloor where diversity is most at risk to future marine global change, noting that the existence of clear regional differences cautions against applying global thresholds.

Oxygen, ecology, and the Cambrian radiation of animals.

The Proterozoic-Cambrian transition records the appearance of essentially all animal body plans (phyla), yet to date no single hypothesis adequately explains both the timing of the event and the evident increase in diversity and disparity. Ecological triggers focused on escalatory predator-prey "arms races" can explain the evolutionary pattern but not its timing, whereas environmental triggers, particularly ocean/atmosphere oxygenation, do the reverse. Using modern oxygen minimum zones as an analog for Proterozoic oceans, we explore the effect of low oxygen levels on the feeding ecology of polychaetes, the dominant macrofaunal animals in deep-sea sediments. Here we show that low oxygen is clearly linked to low proportions of carnivores in a community and low diversity of carnivorous taxa, whereas higher oxygen levels support more complex food webs. The recognition of a physiological control on carnivory therefore links environmental triggers and ecological drivers, providing an integrated explanation for both the pattern and timing of Cambrian animal radiation.

Biotic replacement and mass extinction of the Ediacara biota.

The latest Neoproterozoic extinction of the Ediacara biota has been variously attributed to catastrophic removal by perturbations to global geochemical cycles, 'biotic replacement' by Cambrian-type ecosystem engineers, and a taphonomic artefact. We perform the first critical test of the 'biotic replacement' hypothesis using combined palaeoecological and geochemical data collected from the youngest Ediacaran strata in southern Namibia. We find that, even after accounting for a variety of potential sampling and taphonomic biases, the Ediacaran assemblage preserved at Farm Swartpunt has significantly lower genus richness than older assemblages. Geochemical and sedimentological analyses confirm an oxygenated and non-restricted palaeoenvironment for fossil-bearing sediments, thus suggesting that oxygen stress and/or hypersalinity are unlikely to be responsible for the low diversity of communities preserved at Swartpunt. These combined analyses suggest depauperate communities characterized the latest Ediacaran and provide the first quantitative support for the biotic replacement model for the end of the Ediacara biota. Although more sites (especially those recording different palaeoenvironments) are undoubtedly needed, this study provides the first quantitative palaeoecological evidence to suggest that evolutionary innovation, ecosystem engineering and biological interactions may have ultimately caused the first mass extinction of complex life.

miRNAs: small genes with big potential in metazoan phylogenetics.

microRNAs (miRNAs) are a key component of gene regulatory networks and have been implicated in the regulation of virtually every biological process found in multicellular eukaryotes. What makes them interesting from a phylogenetic perspective is the high conservation of primary sequence between taxa, their accrual in metazoan genomes through evolutionary time, and the rarity of secondary loss in most metazoan taxa. Despite these properties, the use of miRNAs as phylogenetic markers has not yet been discussed within a clear conceptual framework. Here we highlight five properties of miRNAs that underlie their utility in phylogenetics: 1) The processes of miRNA biogenesis enable the identification of novel miRNAs without prior knowledge of sequence; 2) The continuous addition of miRNA families to metazoan genomes through evolutionary time; 3) The low level of secondary gene loss in most metazoan taxa; 4) The low substitution rate in the mature miRNA sequence; and 5) The small probability of convergent evolution of two miRNAs. Phylogenetic analyses using both Bayesian and parsimony methods on a eumetazoan miRNA data set highlight the potential of miRNAs to become an invaluable new tool, especially when used as an additional line of evidence, to resolve previously intractable nodes within the tree of life.

Deep-water first occurrences of Ediacara biota prior to the Shuram carbon isotope excursion in the Wernecke Mountains, Yukon, Canada.

Ediacara-type macrofossils appear as early as ~575 Ma in deep-water facies of the Drook Formation of the Avalon Peninsula, Newfoundland, and the Nadaleen Formation of Yukon and Northwest Territories, Canada. Our ability to assess whether a deep-water origination of the Ediacara biota is a genuine reflection of evolutionary succession, an artifact of an incomplete stratigraphic record, or a bathymetrically controlled biotope is limited by a lack of geochronological constraints and detailed shelf-to-slope transects of Ediacaran continental margins. The Ediacaran Rackla Group of the Wernecke Mountains, NW Canada, represents an ideal shelf-to-slope depositional system to understand the spatiotemporal and environmental context of Ediacara-type organisms' stratigraphic occurrence. New sedimentological and paleontological data presented herein from the Wernecke Mountains establish a stratigraphic framework relating shelfal strata in the Goz/Corn Creek area to lower slope deposits in the Nadaleen River area. We report new discoveries of numerous Aspidella hold-fast discs, indicative of frondose Ediacara organisms, from deep-water slope deposits of the Nadaleen Formation stratigraphically below the Shuram carbon isotope excursion (CIE) in the Nadaleen River area. Such fossils are notably absent in coeval shallow-water strata in the Goz/Corn Creek region despite appropriate facies for potential preservation. The presence of pre-Shuram CIE Ediacara-type fossils occurring only in deep-water facies within a basin that has equivalent well-preserved shallow-water facies provides the first stratigraphic paleobiological support for a deep-water origination of the Ediacara biota. In contrast, new occurrences of Ediacara-type fossils (including juvenile fronds, Beltanelliformis, Aspidella, annulated tubes, and multiple ichnotaxa) are found above the Shuram CIE in both deep- and shallow-water deposits of the Blueflower Formation. Given existing age constraints on the Shuram CIE, it appears that Ediacaran organisms may have originated in the deeper ocean and lived there for up to ~15 million years before migrating into shelfal environments in the terminal Ediacaran. This indicates unique ecophysiological constraints likely shaped the initial habitat preference and later environmental expansion of the Ediacara biota.

Lithium isotopic constraints on the evolution of continental clay mineral factory and marine oxygenation in the earliest Paleozoic Era.

The evolution of oxygen cycles on Earth's surface has been regulated by the balance between molecular oxygen production and consumption. The Neoproterozoic-Paleozoic transition likely marks the second rise in atmospheric and oceanic oxygen levels, widely attributed to enhanced burial of organic carbon. However, it remains disputed how marine organic carbon production and burial respond to global environmental changes and whether these feedbacks trigger global oxygenation during this interval. Here, we report a large lithium isotopic and elemental dataset from marine mudstones spanning the upper Neoproterozoic to middle Cambrian [~660 million years ago (Ma) to 500 Ma]. These data indicate a dramatic increase in continental clay formation after ~525 Ma, likely linked to secular changes in global climate and compositions of the continental crust. Using a global biogeochemical model, we suggest that intensified continental weathering and clay delivery to the oceans could have notably increased the burial efficiency of organic carbon and facilitated greater oxygen accumulation in the earliest Paleozoic oceans.

The identification of microRNAs in calcisponges: independent evolution of microRNAs in basal metazoans.

We present the discovery of microRNAs (miRNAs) in the calcisponges Sycon and Leucosolenia (phylum Calcarea), and potential miRNAs in the homoscleromorph Oscarella carmela (Phylum Homoscleromorpha), expanding the complement of poriferan miRNAs previously known only from the siliceous sponges (demosponges and hexactinellids). Comparison of these miRNAs with those previously described from silicisponges and eumetazoans reveals that these newly described miRNAs are novel, with each metazoan lineage (Silicea, Calcarea, Homoscleromorpha, and Eumetazoa) characterized by a unique and non-overlapping repertoire of miRNAs (or potential miRNAs as in the case of the homoscleromorphs). Because each group is characterized by a unique repertoire of miRNAs, miRNAs cannot be used to help resolve the contentious issue of sponge mono- versus paraphyly. Further, because all sponges are characterized by a similar repertoire of tissue types and body plan organisation, we hypothesize that the lack of conserved miRNAs amongst the three primary sponge lineages is evidence that cellular differentiation and cell type specificity in sponges are not dependent upon conserved miRNAs, contrary to many known cases in eumetazoans. Finally, we suggest that miRNAs evolved multiple times independently not only among eukaryotes, but even within animals, independently evolved miRNAs representing molecular exaptations of RNAi machinery into pre-existing gene regulatory networks. The role(s) miRNAs play though in sponge biology and evolution remains an open question.

Oxygen availability and body mass modulate ectotherm responses to ocean warming.

In an ocean that is rapidly warming and losing oxygen, accurate forecasting of species' responses must consider how this environmental change affects fundamental aspects of their physiology. Here, we develop an absolute metabolic index (ΦA) that quantifies how ocean temperature, dissolved oxygen and organismal mass interact to constrain the total oxygen budget an organism can use to fuel sustainable levels of aerobic metabolism. We calibrate species-specific parameters of ΦA with physiological measurements for red abalone (Haliotis rufescens) and purple urchin (Strongylocentrotus purpuratus). ΦA models highlight that the temperature where oxygen supply is greatest shifts cooler when water loses oxygen or organisms grow larger, providing a mechanistic explanation for observed thermal preference patterns. Viable habitat forecasts are disproportionally deleterious for red abalone, revealing how species-specific physiologies modulate the intensity of a common climate signal, captured in the newly developed ΦA framework.

Mesoproterozoic surface oxygenation accompanied major sedimentary manganese deposition at 1.4 and 1.1 Ga.

Manganese (Mn) oxidation in marine environments requires oxygen (O2 ) or other reactive oxygen species in the water column, and widespread Mn oxide deposition in ancient sedimentary rocks has long been used as a proxy for oxidation. The oxygenation of Earth's atmosphere and oceans across the Archean-Proterozoic boundary are associated with massive Mn deposits, whereas the interval from 1.8-1.0 Ga is generally believed to be a time of low atmospheric oxygen with an apparent hiatus in sedimentary Mn deposition. Here, we report geochemical and mineralogical analyses from 1.1 Ga manganiferous marine-shelf siltstones from the Bangemall Supergroup, Western Australia, which underlie recently discovered economically significant manganese deposits. Layers bearing Mn carbonate microspheres, comparable with major global Mn deposits, reveal that intense periods of sedimentary Mn deposition occurred in the late Mesoproterozoic. Iron geochemical data suggest anoxic-ferruginous seafloor conditions at the onset of Mn deposition, followed by oxic conditions in the water column as Mn deposition persisted and eventually ceased. These data imply there was spatially widespread surface oxygenation ~1.1 Ga with sufficiently oxic conditions in shelf environments to oxidize marine Mn(II). Comparable large stratiform Mn carbonate deposits also occur in ~1.4 Ga marine siltstones hosted in underlying sedimentary units. These deposits are greater or at least commensurate in scale (tonnage) to those that followed the major oxygenation transitions from the Neoproterozoic. Such a period of sedimentary manganogenesis is inconsistent with a model of persistently low O2 throughout the entirety of the Mesoproterozoic and provides robust evidence for dynamic redox changes in the mid to late Mesoproterozoic.

A molecular palaeobiological hypothesis for the origin of aplacophoran molluscs and their derivation from chiton-like ancestors.

Aplacophorans have long been argued to be basal molluscs. We present a molecular phylogeny, including the aplacophorans Neomeniomorpha (Solenogastres) and Chaetodermomorpha (Caudofoveata), which recovered instead the clade Aculifera (Aplacophora + Polyplacophora). Our relaxed Bayesian molecular clock estimates an Early Ordovician appearance of the aculiferan crown group consistent with the presence of chiton-like molluscs with seven or eight dorsal shell plates by the Late Cambrian (approx. 501-490 Ma). Molecular, embryological and palaeontological data indicate that aplacophorans, as well as chitons, evolved from a paraphyletic assemblage of chiton-like ancestors. The recovery of cephalopods as a sister group to aculiferans suggests that the plesiomorphic condition in molluscs might be a morphology similar to that found in monoplacophorans.

MicroRNAs support a turtle + lizard clade.

Despite much interest in amniote systematics, the origin of turtles remains elusive. Traditional morphological phylogenetic analyses place turtles outside Diapsida-amniotes whose ancestor had two fenestrae in the temporal region of the skull (among the living forms the tuatara, lizards, birds and crocodilians)-and allied with some unfenestrate-skulled (anapsid) taxa. Nonetheless, some morphological analyses place turtles within Diapsida, allied with Lepidosauria (tuatara and lizards). Most molecular studies agree that turtles are diapsids, but rather than allying them with lepidosaurs, instead place turtles near or within Archosauria (crocodilians and birds). Thus, three basic phylogenetic positions for turtles with respect to extant Diapsida are currently debated: (i) sister to Diapsida, (ii) sister to Lepidosauria, or (iii) sister to, or within, Archosauria. Interestingly, although these three alternatives are consistent with a single unrooted four-taxon tree for extant reptiles, they differ with respect to the position of the root. Here, we apply a novel molecular dataset, the presence versus absence of specific microRNAs, to the problem of the phylogenetic position of turtles and the root of the reptilian tree, and find that this dataset unambiguously supports a turtle + lepidosaur group. We find that turtles and lizards share four unique miRNA gene families that are not found in any other organisms' genome or small RNA library, and no miRNAs are found in all diapsids but not turtles, or in turtles and archosaurs but not in lizards. The concordance between our result and some morphological analyses suggests that there have been numerous morphological convergences and reversals in reptile phylogeny, including the loss of temporal fenestrae.