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.

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.

Molecular phylogeny of wood mice (Apodemus, Muridae) in East Asia

Sequences of the mitochondrial cytochrome b (1140 bp) and nuclear IRBP (1152 bp) genes were used to assess the evolutionary history of Apodemus, using the complete set of Asian species. Our results indicate that speciation in Asia involved three radiations, which supports an earlier study. The initial radiation yielded A. argenteus (Japanese endemic), A. gurkha (Nepalese endemic), and the ancestral lineage of the remaining Asian species. This lineage subsequently diverged into four groups: agrarius-chevrieri (agrarius group), draco-latronum-semotus (draco group), A. peninsulae, and A. speciosus (Japanese endemic). The final step consisted of divergence within two species groups as a consequence of the geography of the Yunnan-Guizhou plateau and Taiwan. The ecological ability of two Apodemus-species to inhabit one locality via niche partitioning likely drove the second radiation and shaped the basic geographical pattern seen today: A. argenteus and A. speciosus in Japan, A. agrarius and A. peninsulae in northern China, and the A. agrarius and A. draco groups in southern China. The three radiations are estimated to have occurred 7.5, 6.6, and 1.8-0.8 Mya respectively, using the IRBP clock, based on rat-mouse divergence 12 Mya. (C) 2003 The Linnean Society of London.

Variance-Constrained Capacity of the Molecular Timing Channel with Synchronization Error

Molecular communication is set to play an important role in the design of complex biological and chemical systems. An important class of molecular communication systems is based on the timing channel, where information is encoded in the delay of the transmitted molecule - a synchronous approach. At present, a widely used modeling assumption is the perfect synchronization between the transmitter and the receiver. Unfortunately, this assumption is unlikely to hold in most practical molecular systems. To remedy this, we introduce a clock into the model - leading to the molecular timing channel with synchronization error. To quantify the behavior of this new system, we derive upper and lower bounds on the variance-constrained capacity, which we view as the step between the mean-delay and the peak-delay constrained capacity. By numerically evaluating our bounds, we obtain a key practical insight: the drift velocity of the clock links does not need to be significantly larger than the drift velocity of the information link, in order to achieve the variance-constrained capacity with perfect synchronization.

CLOCK-controlled polyphonic regulation of circadian rhythms through canonical and noncanonical E-boxes.

In mammalian circadian clockwork, the CLOCK-BMAL1 complex binds to DNA enhancers of target genes and drives circadian oscillation of transcription. Here we identified 7,978 CLOCK-binding sites in mouse liver by chromatin immunoprecipitation-sequencing (ChIP-Seq), and a newly developed bioinformatics method, motif centrality analysis of ChIP-Seq (MOCCS), revealed a genome-wide distribution of previously unappreciated noncanonical E-boxes targeted by CLOCK. In vitro promoter assays showed that CACGNG, CACGTT, and CATG(T/C)G are functional CLOCK-binding motifs. Furthermore, we extensively revealed rhythmically expressed genes by poly(A)-tailed RNA-Seq and identified 1,629 CLOCK target genes within 11,926 genes expressed in the liver. Our analysis also revealed rhythmically expressed genes that have no apparent CLOCK-binding site, indicating the importance of indirect transcriptional and posttranscriptional regulations. Indirect transcriptional regulation is represented by rhythmic expression of CLOCK-regulated transcription factors, such as Krüppel-like factors (KLFs). Indirect posttranscriptional regulation involves rhythmic microRNAs that were identified by small-RNA-Seq. Collectively, CLOCK-dependent direct transactivation through multiple E-boxes and indirect regulations polyphonically orchestrate dynamic circadian outputs.