Substitution rate comparisons between grasses and palms: synonymous rate differences at the nuclear gene Adh parallel rate differences at the plastid gene rbcL.

A number of studies have noted that nucleotide substitution rates at the chloroplast-encoded rbcL locus violate the molecular clock principle. Substitution rate variation at this plastid gene is particularly pronounced between palms and grasses; for example, a previous study estimated that substitution rates in rbcL sequences are approximately 5-fold faster in grasses than in palms. To determine whether a proportionate change in substitution rates also occurs in plant nuclear genes, we characterized nucleotide substitution rates in palm and grass sequences for the nuclear gene Adh. In this article, we report that palm sequences evolve at a rate of 2.61 x 10(-9) substitution per synonymous site per year, a rate which is slower than most plant nuclear genes. Grass Adh sequences evolve approximately 2.5-fold faster than palms at synonymous sites. Thus, synonymous rates in nuclear Adh genes show a marked decrease in palms relative to grasses, paralleling the pattern found at the plastid rbcL locus. This shared pattern indicates that synonymous rates are correlated between a nuclear and a plastid gene. Remarkably, nonsynonymous rates do not show this correlation. Nonsynonymous rates vary between two duplicated grass Adh loci, and nonsynonymous rates at the palm Adh locus are not markedly reduced relative to grasses.

Middle East respiratory syndrome coronavirus neutralising serum antibodies in dromedary camels: a comparative serological study

A new betacoronavirus-Middle East respiratory syndrome coronavirus (MERS-CoV)-has been identified in patients with severe acute respiratory infection. Although related viruses infect bats, molecular clock analyses have been unable to identify direct ancestors of MERS-CoV. Anecdotal exposure histories suggest that patients had been in contact with dromedary camels or goats. We investigated possible animal reservoirs of MERS-CoV by assessing specific serum antibodies in livestock. METHODS: We took sera from animals in the Middle East (Oman) and from elsewhere (Spain, Netherlands, Chile). Cattle (n=80), sheep (n=40), goats (n=40), dromedary camels (n=155), and various other camelid species (n=34) were tested for specific serum IgG by protein microarray using the receptor-binding S1 subunits of spike proteins of MERS-CoV, severe acute respiratory syndrome coronavirus, and human coronavirus OC43. Results were confirmed by virus neutralisation tests for MERS-CoV and bovine coronavirus. FINDINGS: 50 of 50 (100%) sera from Omani camels and 15 of 105 (14%) from Spanish camels had protein-specific antibodies against MERS-CoV spike. Sera from European sheep, goats, cattle, and other camelids had no such antibodies. MERS-CoV neutralising antibody titres varied between 1/320 and 1/2560 for the Omani camel sera and between 1/20 and 1/320 for the Spanish camel sera. There was no evidence for cross-neutralisation by bovine coronavirus antibodies. INTERPRETATION: MERS-CoV or a related virus has infected camel populations. Both titres and seroprevalences in sera from different locations in Oman suggest widespread infection. FUNDING: European Union, European Centre For Disease Prevention and Control, Deutsche Forschungsgemeinschaft.

Sampling information, length measurements, DNA alignments, models and phylogenetics trees for Oithona similis s.l.

Traditionally, many small-sized copepod species are considered to be widespread, bipolar or cosmopolitan. However, these large-scale distribution patterns need to be re-examined in view of increasing evidence of cryptic and pseudo-cryptic speciation in pelagic copepods. Here, we present a phylogeographic study of Oithona similis s.l. populations from the Arctic Ocean, the Southern Ocean and its northern boundaries, the North Atlantic and the Mediterrranean Sea. O. similis s.l. is considered as one of the most abundant species in temperate to polar oceans and acts as an important link in the trophic network between the microbial loop and higher trophic levels such as fish larvae. Two gene fragments were analysed: the mitochondrial cytochrome oxidase c subunit I (COI), and the nuclear ribosomal 28S genetic marker. Seven distinct, geographically delimitated, mitochondrial lineages could be identified, with divergences among the lineages ranging from 8 to 24 %, thus representing most likely cryptic or pseudocryptic species within O. similis s.l. Four lineages were identified within or close to the borders of the Southern Ocean, one lineage in the Arctic Ocean and two lineages in the temperate Northern hemisphere. Surprisingly the Arctic lineage was more closely related to lineages from the Southern hemisphere than to the other lineages from the Northern hemisphere, suggesting that geographic proximity is a rather poor predictor of how closely related the clades are on a genetic level. Molecular clock application revealed that the evolutionary history of O. similis s.l. is possibly closely associated with the reorganization of the ocean circulation in the mid Miocene and may be an example of allopatric speciation in the pelagic zone.

Diversity in a honey bee pathogen: first report of a third master variant of the Deformed Wing Virus quasispecies

Treatment of emerging RNA viruses is hampered by the high mutation and replication rates that enable these viruses to operate as a quasispecies. Declining honey bee populations have been attributed to the ectoparasitic mite Varroa destructor and its affiliation with Deformed Wing Virus (DWV). In the current study we use next-generation sequencing to investigate the DWV quasispecies in an apiary known to suffer from overwintering colony losses. We show that the DWV species complex is made up of three master variants. Our results indicate that a new DWV Type C variant is distinct from the previously described types A and B, but together they form a distinct clade compared with other members of the Iflaviridae. The molecular clock estimation predicts that Type C diverged from the other variants ~319 years ago. The discovery of a new master variant of DWV has important implications for the positive identification of the true pathogen within global honey bee populations.

Diversity in a honey bee pathogen: first report of a third master variant of the Deformed Wing Virus quasispecies

Treatment of emerging RNA viruses is hampered by the high mutation and replication rates that enable these viruses to operate as a quasispecies. Declining honey bee populations have been attributed to the ectoparasitic mite Varroa destructor and its affiliation with Deformed Wing Virus (DWV). In the current study we use next-generation sequencing to investigate the DWV quasispecies in an apiary known to suffer from overwintering colony losses. We show that the DWV species complex is made up of three master variants. Our results indicate that a new DWV Type C variant is distinct from the previously described types A and B, but together they form a distinct clade compared with other members of the Iflaviridae. The molecular clock estimation predicts that Type C diverged from the other variants ~319 years ago. The discovery of a new master variant of DWV has important implications for the positive identification of the true pathogen within global honey bee populations.

Host jumps shaped the diversity of extant rust fungi (Pucciniales)

* The aim of this study was to determine the evolutionary time line for rust fungi and date key speciation events using a molecular clock. Evidence is provided that supports a contemporary view for a recent origin of rust fungi, with a common ancestor on a flowering plant. * Divergence times for > 20 genera of rust fungi were studied with Bayesian evolutionary analyses. A relaxed molecular clock was applied to ribosomal and mitochondrial genes, calibrated against estimated divergence times for the hosts of rust fungi, such as Acacia (Fabaceae), angiosperms and the cupressophytes. * Results showed that rust fungi shared a most recent common ancestor with a mean age between 113 and 115 million yr. This dates rust fungi to the Cretaceous period, which is much younger than previous estimations. Host jumps, whether taxonomically large or between host genera in the same family, most probably shaped the diversity of rust genera. Likewise, species diversified by host shifts (through coevolution) or via subsequent host jumps. This is in contrast to strict coevolution with their hosts. * Puccinia psidii was recovered in Sphaerophragmiaceae, a family distinct from Raveneliaceae, which were regarded as confamilial in previous studies.

Concordance in phylogeography and ecological niche modelling identify dispersal corridors for reptiles in arid Australia

Aim: Using the rock-specialist agamid Ctenophorus caudicinctus as a model, we test hypothesized biogeographical dispersal corridors for lizards in the Australian arid zone (across the western sand deserts), and assess how these dispersal routes have shaped phylogeographical structuring. Location: Arid and semi-arid Australia. Methods: We sequenced a c. 1400 bp fragment of mtDNA (ND2) for 134 individuals of C. caudicinctus as well as a subset of each of the mtDNA clades for five nuclear loci (BDNF, BACH1, GAPD, NTF3, and PRLR). We used phylogenetic methods to assess biogeographical patterns within C. caudicinctus, including relaxed molecular clock analyses to estimate divergence times. Ecological niche modelling (Maxent) was employed to estimate the current distribution of suitable climatic envelopes for each lineage. Results: Phylogenetic analyses identified two deeply divergent mtDNA clades within C. caudicinctus - an eastern and western clade - separated by the Western Australian sand deserts. However, divergences pre-date the Pleistocene sand deserts. Phylogenetic analyses of the nuclear DNA data sets generally support major mtDNA clades, suggesting past connections between the western C. c. caudicinctus populations in far eastern Pilbara (EP) and the lineages to the east of the sand deserts. Ecological niche modelling supports the continued suitability of climatic conditions between the Central Ranges and the far EP for C. c. graafi. Main conclusions: Estimates of lineage ages provide evidence of divergence between eastern and western clades during the Miocene with subsequent secondary contact during the Pliocene. Our results suggest that this secondary contact occurred via dispersal between the Central Ranges and the far EP, rather than the more southerly Giles Corridor. These events precede the origins of the western sand deserts and divergence patterns instead appear associated with Miocene and Pliocene climate change.

The impact and control of malignant catarrhal fever in Tanzania

Malignant Catarrhal Fever (MCF), an often-lethal infectious disease, presents as a variable complex of lesions in susceptible ungulate species. The disease is caused by a -herpesvirus following transmission from an inapparent carrier host. Two major epidemiological forms exist: wildebeest-associated MCF (WA-MCF), in which the virus is transmitted to susceptible species by wildebeest calves less than approximately four months of age, and sheepassociated MCF (SA-MCF) in which the virus is spread by sheep (primarily adolescents). Due to the lack of an in-vitro propagation system for the causative agent of the more economically significant SA-MCF, and with the expectation that cross-protective immunity may be provided, vaccine development has focused on the more easily propagated alcelaphine herpesvirus-1 (AlHV-1) that causes WA-MCF. In 2008 a direct viral challenge trial showed that a novel vaccine, employing an attenuated AlHV-1 (atAlHV-1) `C5000 virus strain, protected British Friesian-Holstein (FH) cattle against an intranasal challenge with virulent AlHV-1 `C5000 virus. For cattle keeping people living near wildebeest calving areas in sub-Saharan Africa an effective vaccine would have value as it would release them from the costly annual disease avoidance strategy of having to move their herds away from the oncoming wildebeest. On the other hand, an effective vaccine will release herd owners from the need to avoid MCF, allowing them to graze their cattle alongside wildebeest on the highly nutritious pastures of the calving areas. As such conservationists have raised concerns that the development of a vaccine might lead to detrimental grazing competition. The principle objective of this study was to test the novel vaccine on Tanzanian shorthorn zebu cross cattle (SZC).We did this firstly using a natural challenge field trial (Chapter Two) which demonstrated that immunisation with the atAlHV-1 vaccine was well tolerated and induced an oro-nasopharyngeal AlHV-1-specific and -neutralising antibody response. This resulted in an immunity in SZC cattle that was partially protective and reduced naturally transmitted infection by 56%. We also demonstrated that non-fatal infections occurred with a much higher frequency than previously thought. Because the calculated efficacy of the vaccine was less than that seen in British FH cattle we wanted to determine whether host factors, particular to SZC cattle, had impacted the outcomes of the field trial. To do this we repeated the 2008 direct viral challenge trial using SZC cattle (Chapter Four). During this trial we also investigated whether the recombinant bacterial flagellin monomer (FliC), when used as an adjuvant, might improve the vaccine’s efficacy. The findings from this trial indicated that direct challenge with pathogenic AlHV-1 is effective at inducing MCF in SZC cattle and that FliC is not an appropriate adjuvant for this vaccine. Furthermore, with less control group cattle dying of MCF than expected we speculate that SZC cattle may have a degree of resistance to MCF that affords them protection from infection and developing fatal disease. In Chapter Three we investigated aspects of the epidemiology of MCF, specifically whether wildebeest placenta, long implicated by Maasai cattle owners as a source of MCF, might play a role in viral transmission. Additionally, through comparative sequence analysis, at two specific genes (A9.5 and ORF50) of wild-type and atAlHV-1, we investigated whether the `C5000 strain, the source of which was taken from Africa more than 40 years ago, was appropriate for vaccine development. The detection of AlHV-1 virus in approximately 50% of placentae indicated that infection can occur in-utero and that this tissue might play a role in disease transmission. And, despite describing three new alleles of the A9.5 gene (supporting previous evidence that this gene is polymorphic and encodes a secretory protein with interleukin-4 as the major homologue), the observation that the most frequently detected haplotypes, in both wild-type and attenuated AlHV-1, were identical suggests that AlHV-1 has a slow molecular clock and that the attenuated strain was appropriate for vaccine development. In Chapter Five we present the first quantitative assessment of the annual MCF avoidance costs that Maasai pastoralists incur. In particular we estimated that as a result of MCF avoidance 64% of the total daily milk yield during the MCF season was not available to be used by the 81% of the family unit remaining at the permanent boma. This represents an upper-bound loss of approximately 8% of a household0s annual income. Despite these considerable losses we concluded that, given an incidence of fatal MCF in cattle living in wildebeest calving areas of 5% to 10%, if herd owners were to stop trying to avoid MCF by allowing their cattle to graze alongside wildebeest, any gains made through increased availability of milk, improved body condition and reduced energy demands would be offset by an increase in MCF-incidence. With the development of an effective vaccine, however, this alternative strategy might become optimal. The overall conclusion we draw therefore is that, despite the substantial costs incurred each year avoiding MCF, the partial protection afforded by the novel vaccine strategy is not sufficient to warrant a wholesale change in disease avoidance strategy. Nonetheless, even the partial protection provided by this vaccine could be of value to protect animals that cannot be moved, for example where some of the herd remain at the boma to provide milk or where land-use changes make traditional disease avoidance difficult. Furthermore, the vaccine may offer a feasible solution to some of the current land-use challenges and conflicts, providing a degree of protection to valuable livestock where avoidance strategies are not possible, but with less risk of precipitating the potentially damaging environmental consequences, such as overgrazing of highly nutritious seasonal pastures, that might result if herd owners decide they no longer need to avoid wildebeest.

Evolutionary History of the PER3 Variable Number of Tandem Repeats (VNTR): Idiosyncratic Aspect of Primate Molecular Circadian Clock

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Molecular coevolution within a Drosophila clock gene

The period (per) gene in Drosophila melanogaster provides an integral component of biological rhythmicity and encodes a protein that includes a repetitive threonine-glycine (Thr-Gly) tract. Similar repeats are found in the frq and wc2 clock genes of Neurospora crassa and in the mammalian per homologues, but their circadian functions are unknown. In Drosophilids, the length of the Thr-Gly repeat varies widely between species, and sequence comparisons have suggested that the repeat length coevolves with the immediately flanking amino acids. A functional test of the coevolution hypothesis was performed by generating several hybrid per transgenes between Drosophila pseudoobscura and D. melanogaster, whose repetitive regions differ in length by about 150 amino acids. The positions of the chimeric junctions were slightly altered in each transgene. Transformants carrying per constructs in which the repeat of one species was juxtaposed next to the flanking region of the other were almost arrhythmic or showed a striking temperature sensitivity of the circadian period. In contrast, transgenes in which the repeat and flanking regions were conspecific gave wild-type levels of circadian rescue. These results support the coevolutionary interpretation of the interspecific sequence changes in this region of the PER molecule and reveal a functional dimension to this process related to the clock’s temperature compensation.

The molecular basis of temperature compensation in the Arabidopsis circadian clock

Circadian clocks maintain robust and accurate timing over a broad range of physiological temperatures, a characteristic termed temperature compensation. In Arabidopsis thaliana, ambient temperature affects the rhythmic accumulation of transcripts encoding the clock components TIMING OF CAB EXPRESSION1 (TOC1), GIGANTEA (GI), and the partially redundant genes CIRCADIAN CLOCK ASSOCIATED1 (CCA1) and LATE ELONGATED HYPOCOTYL (LHY). The amplitude and peak levels increase for TOC1 and GI RNA rhythms as the temperature increases (from 17 to 27 degrees C), whereas they decrease for LHY. However, as temperatures decrease ( from 17 to 12 degrees C), CCA1 and LHY RNA rhythms increase in amplitude and peak expression level. At 27 degrees C, a dynamic balance between GI and LHY allows temperature compensation in wild-type plants, but circadian function is impaired in Ihy and gi mutant plants. However, at 12 degrees C, CCA1 has more effect on the buffering mechanism than LHY, as the cca1 and gi mutations impair circadian rhythms more than Ihy at the lower temperature. At 17 degrees C, GI is apparently dispensable for free-running circadian rhythms, although partial GI function can affect circadian period. Numerical simulations using the interlocking-loop model show that balancing LHY/CCA1 function against GI and other evening-expressed genes can largely account for temperature compensation in wild-type plants and the temperature-specific phenotypes of gi mutants.

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.

Four mammal fossil calibrations: Balancing competing palaeontological and molecular considerations

With the introduction of relaxed-clock molecular dating methods, the role of fossil calibration has expanded from providing a timescale, to also informing the models for molecular rate variation across the phylogeny. Here I suggest fossil calibration bounds for four mammal clades, Monotremata (platypus and echidnas), Macropodoidea (kangaroos and potoroos), Caviomorpha-Phiomorpha (South American and African hystricognath rodents), and Chiroptera (bats). In each case I consider sources of uncertainty in the fossil record and provide a molecular dating analysis to examine how the suggested calibration priors are further informed by other mammal fossil calibrations and molecular data.

The circadian oscillator gene GIGANTEA mediates a long-term response of the Arabidopsis thaliana circadian clock to sucrose.

Circadian clocks are 24-h timing devices that phase cellular responses; coordinate growth, physiology, and metabolism; and anticipate the day-night cycle. Here we report sensitivity of the Arabidopsis thaliana circadian oscillator to sucrose, providing evidence that plant metabolism can regulate circadian function. We found that the Arabidopsis circadian system is particularly sensitive to sucrose in the dark. These data suggest that there is a feedback between the molecular components that comprise the circadian oscillator and plant metabolism, with the circadian clock both regulating and being regulated by metabolism. We used also simulations within a three-loop mathematical model of the Arabidopsis circadian oscillator to identify components of the circadian clock sensitive to sucrose. The mathematical studies identified GIGANTEA (GI) as being associated with sucrose sensing. Experimental validation of this prediction demonstrated that GI is required for the full response of the circadian clock to sucrose. We demonstrate that GI acts as part of the sucrose-signaling network and propose this role permits metabolic input into circadian timing in Arabidopsis.