The Accumulation of Reproductive Incompatibilities in African Cichlid Fish

The rate at which different components of reproductive isolation accumulate with divergence time between species has only been studied in a limited, but growing, number of species. We measured premating isolation and hybrid inviability at four different ontogenetic stages from zygotes to adults in interspecific hybrids of 26 pairs of African cichlid species, spanning the entire East African haplochromine radiation. We then used multiple relaxed molecular clock calibrations to translate genetic distances into absolute ages to compare evolutionary rates of different components of reproductive isolation. We find that premating isolation accumulates fast initially but then changes little with increasing genetic distance between species. In contrast, postmating isolation between closely related species is negligible but then accumulates rapidly, resulting in complete hybrid inviability after 4.4/8.5/18.4 million years (my). Thus, the rate at which complete intrinsic incompatibilities arise in this system is orders of magnitude lower than rates of speciation within individual lake radiations. Together these results suggest divergent ecological adaptations may prevent populations from interbreeding and help maintain cichlid species diversity, which may be vulnerable to environmental degradation. By quantifying the capacity to produce viable hybrids between allopatric, distantly related lineages our results also provide an upper divergence time limit for the "hybrid swarm origin" model of adaptive radiation.

Models of DNA sequence evolution

Models of DNA sequence evolution and methods for estimating evolutionary distances are needed for studying the rate and pattern of molecular evolution and for inferring the evolutionary relationships of organisms or genes. In this dissertation, several new models and methods are developed.^ The rate variation among nucleotide sites: To obtain unbiased estimates of evolutionary distances, the rate heterogeneity among nucleotide sites of a gene should be considered. Commonly, it is assumed that the substitution rate varies among sites according to a gamma distribution (gamma model) or, more generally, an invariant+gamma model which includes some invariable sites. A maximum likelihood (ML) approach was developed for estimating the shape parameter of the gamma distribution $(\alpha)$ and/or the proportion of invariable sites $(\theta).$ Computer simulation showed that (1) under the gamma model, $\alpha$ can be well estimated from 3 or 4 sequences if the sequence length is long; and (2) the distance estimate is unbiased and robust against violations of the assumptions of the invariant+gamma model.^ However, this ML method requires a huge amount of computational time and is useful only for less than 6 sequences. Therefore, I developed a fast method for estimating $\alpha,$ which is easy to implement and requires no knowledge of tree. A computer program was developed for estimating $\alpha$ and evolutionary distances, which can handle the number of sequences as large as 30.^ Evolutionary distances under the stationary, time-reversible (SR) model: The SR model is a general model of nucleotide substitution, which assumes (i) stationary nucleotide frequencies and (ii) time-reversibility. It can be extended to SRV model which allows rate variation among sites. I developed a method for estimating the distance under the SR or SRV model, as well as the variance-covariance matrix of distances. Computer simulation showed that the SR method is better than a simpler method when the sequence length $L>1,000$ bp and is robust against deviations from time-reversibility. As expected, when the rate varies among sites, the SRV method is much better than the SR method.^ The evolutionary distances under nonstationary nucleotide frequencies: The statistical properties of the paralinear and LogDet distances under nonstationary nucleotide frequencies were studied. First, I developed formulas for correcting the estimation biases of the paralinear and LogDet distances. The performances of these formulas and the formulas for sampling variances were examined by computer simulation. Second, I developed a method for estimating the variance-covariance matrix of the paralinear distance, so that statistical tests of phylogenies can be conducted when the nucleotide frequencies are nonstationary. Third, a new method for testing the molecular clock hypothesis was developed in the nonstationary case. ^

Different population dynamics of human T cell lymphotropic virus type II in intravenous drug users compared with endemically infected tribes

The phylogeny of human T cell lymphotropic virus type II (HTLV-II) was investigated by using strains isolated from Amerindian and Pygmy tribes, in which the virus is maintained primarily through mother-to-child transmission via breast-feeding, and strains from intravenous drug users (IDUs), in which spread is mainly blood-borne via needle sharing. Molecular clock analysis showed that HTLV-II has two different evolutionary rates with the molecular clock for the virus in IDUs ticking 150–350 times faster than the one in endemically infected tribes: 2.7 × 10−4 compared with 1.71/7.31 × 10−7 nucleotide substitutions per site per year in the long terminal repeat region. This dramatic acceleration of the evolutionary rate seems to be related with the mode of transmission. Mathematical models showed the correlation of these two molecular clocks with an endemic spread of HTLV-II in infected tribes compared with the epidemic spread in IDUs. We also noted a sharp increase in the population size of the virus among IDUs during the last decades probably caused by the worldwide increase in intravenous drug use.

Nonsynonymous substitution in abalone sperm fertilization genes exceeds substitution in introns and mitochondrial DNA

Strong positive Darwinian selection acts on two sperm fertilization proteins, lysin and 18-kDa protein, from abalone (Haliotis). To understand the phylogenetic context for this dramatic molecular evolution, we obtained sequences of mitochondrial cytochrome c oxidase subunit I (mtCOI), and genomic sequences of lysin, 18-kDa, and a G protein subunit. Based on mtDNA differentiation, four north Pacific abalone species diverged within the past 2 million years (Myr), and remaining north Pacific species diverged over a period of 4–20 Myr. Between-species nonsynonymous differences in lysin and 18-kDa exons exceed nucleotide differences in introns by 3.5- to 24-fold. Remarkably, in some comparisons nonsynonymous substitutions in lysin and 18-kDa genes exceed synonymous substitutions in mtCOI. Lysin and 18-kDa intron/exon segments were sequenced from multiple red abalone individuals collected over a 1,200-km range. Only two nucleotide changes and two sites of slippage variation were detected in a total of >29,000 nucleotides surveyed. However, polymorphism in mtCOI and a G protein intron was found in this species. This finding suggests that positive selection swept one lysin allele and one 18-kDa allele to fixation. Similarities between mtCOI and lysin gene trees indicate that rapid adaptive evolution of lysin has occurred consistently through the history of the group. Comparisons with mtCOI molecular clock calibrations suggest that nonsynonymous substitutions accumulate 2–50 times faster in lysin and 18-kDa genes than in rapidly evolving mammalian genes.

Estimation of divergence times from multiprotein sequences for a few mammalian species and several distantly related organisms

When many protein sequences are available for estimating the time of divergence between two species, it is customary to estimate the time for each protein separately and then use the average for all proteins as the final estimate. However, it can be shown that this estimate generally has an upward bias, and that an unbiased estimate is obtained by using distances based on concatenated sequences. We have shown that two concatenation-based distances, i.e., average gamma distance weighted with sequence length (d2) and multiprotein gamma distance (d3), generally give more satisfactory results than other concatenation-based distances. Using these two distance measures for 104 protein sequences, we estimated the time of divergence between mice and rats to be approximately 33 million years ago. Similarly, the time of divergence between humans and rodents was estimated to be approximately 96 million years ago. We also investigated the dependency of time estimates on statistical methods and various assumptions made by using sequence data from eubacteria, protists, plants, fungi, and animals. Our best estimates of the times of divergence between eubacteria and eukaryotes, between protists and other eukaryotes, and between plants, fungi, and animals were 3, 1.7, and 1.3 billion years ago, respectively. However, estimates of ancient divergence times are subject to a substantial amount of error caused by uncertainty of the molecular clock, horizontal gene transfer, errors in sequence alignments, etc.

Estimation of evolutionary distances under stationary and nonstationary models of nucleotide substitution

Estimation of evolutionary distances has always been a major issue in the study of molecular evolution because evolutionary distances are required for estimating the rate of evolution in a gene, the divergence dates between genes or organisms, and the relationships among genes or organisms. Other closely related issues are the estimation of the pattern of nucleotide substitution, the estimation of the degree of rate variation among sites in a DNA sequence, and statistical testing of the molecular clock hypothesis. Mathematical treatments of these problems are considerably simplified by the assumption of a stationary process in which the nucleotide compositions of the sequences under study have remained approximately constant over time, and there now exist fairly extensive studies of stationary models of nucleotide substitution, although some problems remain to be solved. Nonstationary models are much more complex, but significant progress has been recently made by the development of the paralinear and LogDet distances. This paper reviews recent studies on the above issues and reports results on correcting the estimation bias of evolutionary distances, the estimation of the pattern of nucleotide substitution, and the estimation of rate variation among the sites in a sequence.

Accurate reconstruction of a known HIV-1 transmission history by phylogenetic tree analysis.

Phylogenetic analyses are increasingly used in attempts to clarify transmission patterns of human immunodeficiency virus type 1 (HIV-1), but there is a continuing discussion about their validity because convergent evolution and transmission of minor HIV variants may obscure epidemiological patterns. Here we have studied a unique HIV-1 transmission cluster consisting of nine infected individuals, for whom the time and direction of each virus transmission was exactly known. Most of the transmissions occurred between 1981 and 1983, and a total of 13 blood samples were obtained approximately 2-12 years later. The p17 gag and env V3 regions of the HIV-1 genome were directly sequenced from uncultured lymphocytes. A true phylogenetic tree was constructed based on the knowledge about when the transmissions had occurred and when the samples were obtained. This complex, known HIV-1 transmission history was compared with reconstructed molecular trees, which were calculated from the DNA sequences by several commonly used phylogenetic inference methods [Fitch-Margoliash, neighbor-joining, minimum-evolution, maximum-likelihood, maximum-parsimony, unweighted pair group method using arithmetic averages (UPGMA), and a Fitch-Margoliash method assuming a molecular clock (KITSCH)]. A majority of the reconstructed trees were good estimates of the true phylogeny; 12 of 13 taxa were correctly positioned in the most accurate trees. The choice of gene fragment was found to be more important than the choice of phylogenetic method and substitution model. However, methods that are sensitive to unequal rates of change performed more poorly (such as UPGMA and KITSCH, which assume a constant molecular clock). The rapidly evolving V3 fragment gave better reconstructions than p17, but a combined data set of both p17 and V3 performed best. The accuracy of the phylogenetic methods justifies their use in HIV-1 research and argues against convergent evolution and selective transmission of certain virus variants.

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.

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.