Evolution of the diatoms: major steps in their evolution and a review of the supporting molecular and morphological evidence.

Over the years, many reviews of different aspects of diatom biology, ecology and evolution have appeared. Since 1993 many molecular trees have been produced to infer diatom phylogeny. In 2004, Medlin & Kaczmarska revised the systematics of the diatoms based on more than 20 years of consistent recovery of two major clades of diatoms that did not correspond to a traditional concept of centrics and pennates and established three classes of diatoms: Clade 1 = Coscinodiscophyceae (radial centrics) and Clade 2 = Mediophyceae (polar centrics + radial Thalassiosirales) and Bacillariophyceae (pennates). However, under certain analytical conditions, an alternative view of diatom evolution, a grades of clades, has been recovered that suggests a gradual evolution from centric to pennate symmetry. These two schemes of diatom evolution are evaluated in terms of whether or not the criteria advocated by Medlin & Kaczmarska that should be met to recover monophyletic classes have been used. The monophyly of the three diatom classes can only be achieved if (1) a secondary structure of the small subunit (SSU) rRNA gene was used to construct the alignment and not an alignment based on primary structure and (2) multiple outgroups were used. These requirements have not been met in each study of diatom evolution; hence, the grade of clades, which is useful in reconstructing the sequence of evolution, is not useful for accepting the new classification of the diatoms. Evidence for how these two factors affect the recovery of the three monophyletic classes is reviewed here. The three classes have been defined by clear morphological differences primarily based on gametangia and auxospore ontogeny and envelope structure, the presence or absence of a structure (tube process or sternum) associated with the annulus and the location of the cribrum in those genera with loculate areolae. New evidence supporting the three clades is reviewed. Other features of the cell are examined to determine whether they can also be used to support the monophyly of the three classes.

Evolution of the diatoms: major steps in their evolution and a review of the supporting molecular and morphological evidence.

Over the years, many reviews of different aspects of diatom biology, ecology and evolution have appeared. Since 1993 many molecular trees have been produced to infer diatom phylogeny. In 2004, Medlin & Kaczmarska revised the systematics of the diatoms based on more than 20 years of consistent recovery of two major clades of diatoms that did not correspond to a traditional concept of centrics and pennates and established three classes of diatoms: Clade 1 = Coscinodiscophyceae (radial centrics) and Clade 2 = Mediophyceae (polar centrics + radial Thalassiosirales) and Bacillariophyceae (pennates). However, under certain analytical conditions, an alternative view of diatom evolution, a grades of clades, has been recovered that suggests a gradual evolution from centric to pennate symmetry. These two schemes of diatom evolution are evaluated in terms of whether or not the criteria advocated by Medlin & Kaczmarska that should be met to recover monophyletic classes have been used. The monophyly of the three diatom classes can only be achieved if (1) a secondary structure of the small subunit (SSU) rRNA gene was used to construct the alignment and not an alignment based on primary structure and (2) multiple outgroups were used. These requirements have not been met in each study of diatom evolution; hence, the grade of clades, which is useful in reconstructing the sequence of evolution, is not useful for accepting the new classification of the diatoms. Evidence for how these two factors affect the recovery of the three monophyletic classes is reviewed here. The three classes have been defined by clear morphological differences primarily based on gametangia and auxospore ontogeny and envelope structure, the presence or absence of a structure (tube process or sternum) associated with the annulus and the location of the cribrum in those genera with loculate areolae. New evidence supporting the three clades is reviewed. Other features of the cell are examined to determine whether they can also be used to support the monophyly of the three classes.

Molecular evolutionary dynamics of Ross River virus and implications for vaccine efficacy

Ross River virus (RRV) is a mosquito-borne member of the genus Alphavirus that causes epidemic polyarthritis in humans, costing the Australian health system at least US$10 million annually. Recent progress in RRV vaccine development requires accurate assessment of RRV genetic diversity and evolution, particularly as they may affect the utility of future vaccination. In this study, we provide novel RRV genome sequences and investigate the evolutionary dynamics of RRV from time-structured E2 gene datasets. Our analysis indicates that, although RRV evolves at a similar rate to other alphaviruses (mean evolutionary rate of approx. 8x10(-4) nucleotide substitutions per site year(-1)), the relative genetic diversity of RRV has been continuously low through time, possibly as a result of purifying selection imposed by replication in a wide range of natural host and vector species. Together, these findings suggest that vaccination against RRV is unlikely to result in the rapid antigenic evolution that could compromise the future efficacy of current RRV vaccines.

Multiple evolutionary rate classes in animal genome evolution

The proportion of functional sequence in the human genome is currently a subject of debate. The most widely accepted figure is that approximately 5% is under purifying selection. In Drosophila, estimates are an order of magnitude higher, though this corresponds to a similar quantity of sequence. These estimates depend on the difference between the distribution of genomewide evolutionary rates and that observed in a subset of sequences presumed to be neutrally evolving. Motivated by the widening gap between these estimates and experimental evidence of genome function, especially in mammals, we developed a sensitive technique for evaluating such distributions and found that they are much more complex than previously apparent. We found strong evidence for at least nine well-resolved evolutionary rate classes in an alignment of four Drosophila species and at least seven classes in an alignment of four mammals, including human. We also identified at least three rate classes in human ancestral repeats. By positing that the largest of these ancestral repeat classes is neutrally evolving, we estimate that the proportion of nonneutrally evolving sequence is 30% of human ancestral repeats and 45% of the aligned portion of the genome. However, we also question whether any of the classes represent neutrally evolving sequences and argue that a plausible alternative is that they reflect variable structure-function constraints operating throughout the genomes of complex organisms.

Iron biogeochemistry and associated greenhouse gas evolution in a forested subtropical Australian coastal catchment : Poona Creek, Southeast Queensland

Iron (Fe) is the fourth most abundant element in the Earth’s crust. Excess Fe mobilization from terrestrial into aquatic systems is of concern for deterioration of water quality via biofouling and nuisance algal blooms in coastal and marine systems. Substantial Fe dissolution and transport involve alternate Fe(II) oxidation followed by Fe(III) reduction, with a diversity of Bacteria and Archaea acting as the key catalyst. Microbially-mediated Fe cycling is of global significance with regard to cycles of carbon (C), sulfur (S) and manganese (Mn). However, knowledge regarding microbial Fe cycling in circumneutral-pH habitats that prevail on Earth has been lacking until recently. In particular, little is known regarding microbial function in Fe cycling and associated Fe mobilization and greenhouse (CO2 and CH4, GHG) evolution in subtropical Australian coastal systems where microbial response to ambient variations such as seasonal flooding and land use changes is of concern. Using the plantation-forested Poona Creek catchment on the Fraser Coast of Southeast Queensland (SEQ), this research aimed to 1) study Fe cycling-associated bacterial populations in diverse terrestrial and aquatic habitats of a representative subtropical coastal circumneutral-pH (4–7) ecosystem; and 2) assess potential impacts of Pinus plantation forestry practices on microbially-mediated Fe mobilization, organic C mineralization and associated GHG evolution in coastal SEQ. A combination of wet-chemical extraction, undisturbed core microcosm, laboratory bacterial cultivation, microscopy and 16S rRNA-based molecular phylogenetic techniques were employed. The study area consisted primarily of loamy sands, with low organic C and dissolved nutrients. Total reactive Fe was abundant and evenly distributed within soil 0–30 cm profiles. Organic complexation primarily controlled Fe bioavailability and forms in well-drained plantation soils and water-logged, native riparian soils, whereas tidal flushing exerted a strong “seawater effect” in estuarine locations and formed a large proportion of inorganic Fe(III) complexes. There was a lack of Fe(II) sources across the catchment terrestrial system. Mature, first-rotation plantation clear-felling and second-rotation replanting significantly decreased organic matter and poorly crystalline Fe in well-drained soils, although variations in labile soil organic C fractions (dissolved organic C, DOC; and microbial biomass C, MBC) were minor. Both well-drained plantation soils and water-logged, native-vegetation soils were inhabited by a variety of cultivable, chemotrophic bacterial populations capable of C, Fe, S and Mn metabolism via lithotrophic or heterotrophic, (micro)aerobic or anaerobic pathways. Neutrophilic Fe(III)-reducing bacteria (FeRB) were most abundant, followed by aerobic, heterotrophic bacteria (heterotrophic plate count, HPC). Despite an abundance of FeRB, cultivable Fe(II)-oxidizing bacteria (FeOB) were absent in associated soils. A lack of links between cultivable Fe, S or Mn bacterial densities and relevant chemical measurements (except for HPC correlated with DOC) was likely due to complex biogeochemical interactions. Neither did variations in cultivable bacterial densities correlate with plantation forestry practices, despite total cultivable bacterial densities being significantly lower in estuarine soils when compared with well-drained plantation soils and water-logged, riparian native-vegetation soils. Given that bacterial Fe(III) reduction is the primary mechanism of Fe oxide dissolution in soils upon saturation, associated Fe mobilization involved several abiotic and biological processes. Abiotic oxidation of dissolved Fe(II) by Mn appeared to control Fe transport and inhibit Fe dissolution from mature, first-rotation plantation soils post-saturation. Such an effect was not observed in clear-felled and replanted soils associated with low SOM and potentially low Mn reactivity. Associated GHG evolution post-saturation mainly involved variable CO2 emissions, with low, but consistently increasing CH4 effluxes in mature, first-rotation plantation soil only. In comparison, water-logged soils in the riparian native-vegetation buffer zone functioned as an important GHG source, with high potentials for Fe mobilization and GHG, particularly CH4 emissions in riparian loam soils associated with high clay and crystalline Fe fractions. Active Fe–C cycling was unlikely to occur in lower-catchment estuarine soils associated with low cultivable bacterial densities and GHG effluxes. As a key component of bacterial Fe cycling, neutrophilic FeOB widely occurred in diverse aquatic, but not terrestrial, habitats of the catchment study area. Stalked and sheathed FeOB resembling Gallionella and Leptothrix were limited to microbial mat material deposited in surface fresh waters associated with a circumneutral-pH seep, and clay-rich soil within riparian buffer zones. Unicellular, Sideroxydans-related FeOB (96% sequence identity) were ubiquitous in surface and subsurface freshwater environments, with highest abundance in estuary-adjacent shallow coastal groundwater water associated with redox transition. The abundance of dissolved C and Fe in the groundwater-dependent system was associated with high numbers of cultivable anaerobic, heterotrophic FeRB, microaerophilic, putatively lithotrophic FeOB and aerobic, heterotrophic bacteria. This research represents the first study of microbial Fe cycling in diverse circumneutral-pH environments (terrestrial–aquatic, freshwater–estuarine, surface–subsurface) of a subtropical coastal ecosystem. It also represents the first study of its kind in the southern hemisphere. This work highlights the significance of bacterial Fe(III) reduction in terrestrial, and bacterial Fe(II) oxidation in aquatic catchment Fe cycling. Results indicate the risk of promotion of Fe mobilization due to plantation clear-felling and replanting, and GHG emissions associated with seasonal water-logging. Additional significant outcomes were also achieved. The first direct evidence for multiple biomineralization patterns of neutrophilic, microaerophilic, unicellular FeOB was presented. A putatively pure culture, which represents the first cultivable neutrophilic FeOB from the southern hemisphere, was obtained as representative FeOB ubiquitous in diverse catchment aquatic habitats.

Review ["Telling the Evolutionary Time : Molecular Clocks and the Fossil Record" edited by Philip C. J. Donoghue and M. Paul Smith]

Determining the temporal scale of biological evolution has traditionally been the preserve of paleontology, with the timing of species originations and major diversifications all being read from the fossil record. However, the ages of the earliest (correctly identified) records will underestimate actual origins due to the incomplete nature of the fossil record and the necessity for lineages to have evolved sufficiently divergent morphologies in order to be distinguished. The possibility of inferring divergence times more accurately has been promoted by the idea that the accumulation of genetic change between modern lineages can be used as a molecular clock (Zuckerkandl and Pauling, 1965). In practice, though, molecular dates have often been so old as to be incongruent even with liberal readings of the fossil record. Prominent examples include inferred diversifications of metazoan phyla hundreds of millions of years before their Cambrian fossil record appearances (e.g., Nei et al., 2001) and a basal split between modern birds (Neoaves) that is almost double the age of their earliest recognizable fossils (e.g., Cooper and Penny, 1997).

Evolution of the extinct Sabretooths and the American cheetah-like cat

The sabretooths (Smilodon and Homotherium) and the American cheetah-like cat (Miracinonyx) were the top predators in Late Pleistocene America, but became extinct about 13 thousand years ago. As the evolutionary history of these taxa remains poorly understood , we analysed their phylogenetic relationship to extant felids. In contrast to previous molecular studies , our results show that the sabretooths diverge early and are not closely related to any living cats. This supports their morphological placement in a separate subfamily (Machairodontinae). Despite its remarkable morphological similarity to the African cheetah (Acinonyx jubatus), Miracinonyx appears to have evolved from a puma-like ancestor, presumably in response to similar ecological pressures.

Time dependency of molecular rate estimates and systematic overestimation of recent divergence times

Studies of molecular evolutionary rates have yielded a wide range of rate estimates for various genes and taxa. Recent studies based on population-level and pedigree data have produced remarkably high estimates of mutation rate, which strongly contrast with substitution rates inferred in phylogenetic (species-level) studies. Using Bayesian analysis with a relaxed-clock model, we estimated rates for three groups of mitochondrial data: avian protein-coding genes, primate protein-coding genes, and primate d-loop sequences. In all three cases, we found a measurable transition between the high, short-term (<1–2 Myr) mutation rate and the low, long-term substitution rate. The relationship between the age of the calibration and the rate of change can be described by a vertically translated exponential decay curve, which may be used for correcting molecular date estimates. The phylogenetic substitution rates in mitochondria are approximately 0.5% per million years for avian protein-coding sequences and 1.5% per million years for primate protein-coding and d-loop sequences. Further analyses showed that purifying selection offers the most convincing explanation for the observed relationship between the estimated rate and the depth of the calibration. We rule out the possibility that it is a spurious result arising from sequence errors, and find it unlikely that the apparent decline in rates over time is caused by mutational saturation. Using a rate curve estimated from the d-loop data, several dates for last common ancestors were calculated: modern humans and Neandertals (354 ka; 222–705 ka), Neandertals (108 ka; 70–156 ka), and modern humans (76 ka; 47–110 ka). If the rate curve for a particular taxonomic group can be accurately estimated, it can be a useful tool for correcting divergence date estimates by taking the rate decay into account. Our results show that it is invalid to extrapolate molecular rates of change across different evolutionary timescales, which has important consequences for studies of populations, domestication, conservation genetics, and human evolution.

Branch-length estimation bias misleads molecular dating for a vertebrate mitochondrial phylogeny

Despite recent methodological advances in inferring the time-scale of biological evolution from molecular data, the fundamental question of whether our substitution models are sufficiently well specified to accurately estimate branch-lengths has received little attention. I examine this implicit assumption of all molecular dating methods, on a vertebrate mitochondrial protein-coding dataset. Comparison with analyses in which the data are RY-coded (AG → R; CT → Y) suggests that even rates-across-sites maximum likelihood greatly under-compensates for multiple substitutions among the standard (ACGT) NT-coded data, which has been subject to greater phylogenetic signal erosion. Accordingly, the fossil record indicates that branch-lengths inferred from the NT-coded data translate into divergence time overestimates when calibrated from deeper in the tree. Intriguingly, RY-coding led to the opposite result. The underlying NT and RY substitution model misspecifications likely relate respectively to “hidden” rate heterogeneity and changes in substitution processes across the tree, for which I provide simulated examples. Given the magnitude of the inferred molecular dating errors, branch-length estimation biases may partly explain current conflicts with some palaeontological dating estimates.

Bird evolution: testing the Metaves clade with six new mitochondrial genomes

Background Evolutionary biologists are often misled by convergence of morphology and this has been common in the study of bird evolution. However, the use of molecular data sets have their own problems and phylogenies based on short DNA sequences have the potential to mislead us too. The relationships among clades and timing of the evolution of modern birds (Neoaves) has not yet been well resolved. Evidence of convergence of morphology remain controversial. With six new bird mitochondrial genomes (hummingbird, swift, kagu, rail, flamingo and grebe) we test the proposed Metaves/Coronaves division within Neoaves and the parallel radiations in this primary avian clade. Results Our mitochondrial trees did not return the Metaves clade that had been proposed based on one nuclear intron sequence. We suggest that the high number of indels within the seventh intron of the β-fibrinogen gene at this phylogenetic level, which left a dataset with not a single site across the alignment shared by all taxa, resulted in artifacts during analysis. With respect to the overall avian tree, we find the flamingo and grebe are sister taxa and basal to the shorebirds (Charadriiformes). Using a novel site-stripping technique for noise-reduction we found this relationship to be stable. The hummingbird/swift clade is outside the large and very diverse group of raptors, shore and sea birds. Unexpectedly the kagu is not closely related to the rail in our analysis, but because neither the kagu nor the rail have close affinity to any taxa within this dataset of 41 birds, their placement is not yet resolved. Conclusion Our phylogenetic hypothesis based on 41 avian mitochondrial genomes (13,229 bp) rejects monophyly of seven Metaves species and we therefore conclude that the members of Metaves do not share a common evolutionary history within the Neoaves.

A molecular phylogeny for the Tribe Dacini (Diptera: Tephritidae) : systematic and biogeographic implications

With well over 700 species, the Tribe Dacini is one of the most species-rich clades within the dipteran family Tephritidae, the true fruit flies. Nearly all Dacini belong to one of two very large genera, Dacus Fabricius and Bactrocera Macquart. The distribution of the genera overlap in or around the Indian subcontinent, but the greatest diversity of Dacus is in Africa and the greatest diversity of Bactrocera is in south-east Asia and the Pacific. The monophyly of these two genera has not been rigorously established, with previous phylogenies only including a small number of species and always heavily biased to one genus over the other. Moreover, the subgeneric taxonomy within both genera is complex and the monophyly of many subgenera has not been explicitly tested. Previous hypotheses about the biogeography of the Dacini based on morphological reviews and current distributions of taxa have invoked an out-of-India hypothesis; however this has not been tested in a phylogenetic framework. We attempted to resolve these issues with a dated, molecular phylogeny of 125 Dacini species generated using 16S, COI, COII and white eye genes. The phylogeny shows that Bactrocera is not monophyletic, but rather consists of two major clades: Bactrocera s.s. and the ‘Zeugodacus group of subgenera’ (a recognised, but informal taxonomic grouping of 15 Bactrocera subgenera). This ‘Zeugodacus’ clade is the sister group to Dacus, not Bactrocera and, based on current distributions, split from Dacus before that genus moved into Africa. We recommend that taxonomic consideration be given to raising Zeugodacus to genus level. Supportive of predictions following from the out-of-India hypothesis, the first common ancestor of the Dacini arose in the mid-Cretaceous approximately 80 mya. Major divergence events occurred during the Indian rafting period and diversification of Bactrocera apparently did not begin until after India docked with Eurasia (50–35 mya). In contrast, diversification in Dacus, at approximately 65 mya, apparently began much earlier than predicted by the out-of-India hypothesis, suggesting that, if the Dacini arose on the Indian plate, then ancestral Dacus may have left the plate in the mid to late Cretaceous via the well documented India–Madagascar–Africa migration route. We conclude that the phylogeny does not disprove the predictions of an out-of-India hypothesis for the Dacini, although modification of the original hypothesis is required.

A mitochondrial genome phylogeny of termites (Blattodea: Termitoidae) : robust support for interfamilial relationships and molecular synapomorphies define major clades

Despite their ecological significance as decomposers and their evolutionary significance as the most speciose eusocial insect group outside the Hymenoptera, termite (Blattodea: Termitoidae or Isoptera) evolutionary relationships have yet to be well resolved. Previous morphological and molecular analyses strongly conflict at the family level and are marked by poor support for backbone nodes. A mitochondrial (mt) genome phylogeny of termites was produced to test relationships between the recognised termite families, improve nodal support and test the phylogenetic utility of rare genomic changes found in the termite mt genome. Complete mt genomes were sequenced for 7 of the 9 extant termite families with additional representatives of each of the two most speciose families Rhinotermitidae (3 of 7 subfamilies) and Termitidae (3 of 8 subfamilies). The mt genome of the well supported sister group of termites, the subsocial cockroach Cryptocercus, was also sequenced. A highly supported tree of termite relationships was produced by all analytical methods and data treatment approaches, however the relationship of the termites + Cryptocercus clade to other cockroach lineages was highly affected by the strong nucleotide compositional bias found in termites relative to other dictyopterans. The phylogeny supports previously proposed suprafamilial termite lineages, the Euisoptera and Neoisoptera, a later derived Kalotermitidae as sister group of the Neoisoptera and a monophyletic clade of dampwood (Stolotermitidae, Archotermopsidae) and harvester termites (Hodotermitidae). In contrast to previous termite phylogenetic studies, nodal supports were very high for family-level relationships within termites. Two rare genomic changes in the mt genome control region were found to be molecular synapomorphies for major clades. An elongated stem-loop structure defined the clade Polyphagidae + (Cryptocercus + termites), and a further series of compensatory base changes in this stem loop is synapomorphic for the Neoisoptera. The complicated repeat structures first identified in Reticulitermes, composed of short (A-type) and long (B-type repeats) defines the clade Heterotermitinae + Termitidae, while the secondary loss of A-type repeats is synapomorphic for the non-macrotermitine Termitidae.

Continuum mechanics modelling of microfilament networks with different architectures based on molecular investigation of single F-actin

The actin microfilament plays a critical role in many cellular processes including embryonic development, wound healing, immune response, and tissue development. It is commonly organized in the form of networks whose mechanical properties change with changes in their architecture due to cell evolution processes. This paper presents a new nonlinear continuum mechanics model of single filamentous actin (F-actin) that is based on nanoscale molecular simulations. Following this continuum model of the single F-actin, mechanical properties of differently architected lamellipodia are studied. The results provide insight that can contribute to the understanding of the cell edge motions of living cells.

Molecular characterisation and expression analysis of interferon gamma in response to natural chlamydia infection in the koala, phascolarctos cinereus

Interferon gamma (IFNγ) is a key Th1 cytokine, with a principal role in the immune response against intracellular organisms such as Chlamydia. Along with being responsible for significant morbidity in human populations, Chlamydia is also responsible for wide spread infection and disease in many animal hosts, with reports that many Australian koala subpopulations are endemically infected. An understanding of the role played by IFNγ in koala chlamydial diseases is important for the establishment of better prophylactic and therapeutic approaches against chlamydial infection in this host. A limited number of IFNγ sequences have been published from marsupials and no immune reagents to measure expression have been developed. Through preliminary analysis of the koala transcriptome, we have identified the full coding sequence of the koala IFNγ gene. Transcripts were identified in spleen and lymph node tissue samples. Phylogenetic analysis demonstrated that koala IFNγ is closely related to other marsupial IFNγ sequences and more distantly related to eutherian mammals. To begin to characterise the role of this important cytokine in the koala's response to chlamydial infection, we developed a quantitative real time PCR assay and applied it to a small cohort of koalas with and without active chlamydial disease, revealing significant differences in expression patterns between the groups. Description of the IFNγ sequence from the koala will not only assist in understanding this species' response to its most important pathogen but will also provide further insight into the evolution of the marsupial immune system