Whole-genome sequencing for TB source investigations: principles of ethical precision public health

BACKGROUND: Whole-genome sequencing (WGS) of Mycobacterium tuberculosis allows rapid, accurate inferences about the sources, location and timing of transmission. However, in an era of heightened concern for personal privacy and science distrust, such inferences could result in unintended harm and undermine the public´s trust.METHODS: We held interdisciplinary stakeholder discussions and performed ethical analyses of real-world illustrative cases to identify principles that optimise benefit and mitigate harm of M. tuberculosis WGS-driven TB source investigations.RESULTS: The speed and precision with which real-time WGS can be used to associate M. tuberculosis strains with sensitive information has raised important concerns. While detailed understanding of transmission events could mitigate harm to vulnerable patients and communities when otherwise unfairly blamed for TB outbreaks, the precision of WGS can also identify transmission events resulting in social blame, fear, discrimination, individual or location stigma, and the use of defaming language by the public, politicians and scientists. Public health programmes should balance the need to safeguard privacy with public health goals, transparency and individual rights, including the right to know who infects whom or where.CONCLUSIONS: Ethical challenges raised by real-time WGS-driven TB source investigation requires public health authorities to move beyond their current legal mandate and embrace transparency, privacy and community engagement.

Genomic determinants of speciation and spread of the Mycobacterium tuberculosis complex

Models on how bacterial lineages differentiate increase our understanding of early bacterial speciation events and the genetic loci involved. Here, we analyze the population genomics events leading to the emergence of the tuberculosis pathogen. The emergence is characterized by a combination of recombination events involving core pathogenesis functions and purifying selection on early diverging loci. We identify the phoR gene, the sensor kinase of a two-component system involved in virulence, as a key functional player subject to pervasive positive selection after the divergence of the Mycobacterium tuberculosis complex from its ancestor. Previous evidence showed that phoR mutations played a central role in the adaptation of the pathogen to different host species. Now, we show that phoR mutations have been under selection during the early spread of human tuberculosis, during later expansions, and in ongoing transmission events. Our results show that linking pathogen evolution across evolutionary and epidemiological time scales points to past and present virulence determinants.

Transmission dynamics of HIV-1 subtype B in the Basque Country, Spain

This work was aimed to study the HIV-1 subtype B epidemics in the Basque Country, Spain. 1727 HIV-1 subtype B sequences comprising protease and reverse transcriptase (PR/RT) coding regions, sampled between 2001 and 2008, were analyzed. 156 transmission clusters were detected by means of phylogenetic analyses. Most of them comprised less than 4 individuals and, in total, they included 441 patients. Six clusters comprised 10 or more patients and were further analyzed in order to study their origin and diversification. Four clusters included men who had unprotected homosexual sex (MSM), one group was formed by intravenous drug users (IDUs), and another included both IDUs and people infected through unprotected heterosexual sex (HTs). Most of these clusters originated from the mid-1980s to the mid-1990s. Only one cluster, formed by MSM, originated after 2000. The time between infections was significantly lower in MSM groups than in those containing IDUs (P-value <0.0001). Nucleoside RT and non-nucleoside RT inhibitor (NRTI and NNRTI)-resistance mutations to antiretroviral treatment were found in these six clusters except the most recent MSM group, but only the IDU clusters presented protease inhibitor (PI)-resistance mutations. The most prevalent mutations for each inhibitor class were PI L90M, NRTI T215D/Y/F, and NNRTI K103N, which were also among the most prevalent resistant variants in the whole dataset. In conclusion, while most infections occur as isolated introductions into the population, the number of infections found to be epidemiologically related within the Basque Country is significant. Public health control measures should be reinforced to prevent the further expansion of transmission clusters and resistant mutations occurring within them.

Refined analysis of genetic variability parameters in hepatitis C virus and the ability to predict antiviral treatment response

Hepatitis C virus (HCV) infects approximately 3% of the world population. The chronicity of hepatitis C seems to depend on the level of genetic variability. We have recently (Torres-Puente et al., J Viral Hepat, 2008; 15: 188) reported genetic variability estimates from a large-scale sequence analysis of 67 patients infected with HCV subtypes 1a (23 patients) and 1b (44 patients) and related them to response, or lack of, to alpha-interferon plus ribavirin treatment.. Two HCV genome regions were analysed in samples prior to antiviral therapy, one compressing the three hypervariable regions of the E2 glycoprotein and another one including the interferon sensitive determining region and the V3 domain of the NS5A protein. Haplotype and nucleotide diversity measures showed a clear tendency to higher genetic variability levels in nonresponder than in responder patients. Here, we have refined the analysis of genetic variability (haplotype and nucleotide diversity, number of haplotypes and mutations) by considering their distribution in each of the biologically meaningful subregions mentioned above, as well as in their surrounding and intervening regions. Variability levels are very heterogeneous among the different subregions, being higher for nonresponder patients. Interestingly, significant differences were detected in the biologically relevant regions, but also in the surrounding regions, suggesting that the level of variability of the whole HCV genome, rather than exclusively that from the hypervariable regions, is the main indicator of the treatment response. Finally, the number of haplotypes and mutations seem to be better discriminators than haplotype and nucleotide diversity, especially in the NS5A region.

Genetic variability in hepatitis C virus and its role in antiviral treatment response

Hepatitis C virus (HCV) is a major health problem worldwide, infecting an estimated 170 million people. The high genetic variability of HCV contributes to the chronicity of hepatitis C. Here, we report results from a large-scale sequence analysis of 67 patients infected with HCV genotype 1, 23 with subtype 1a and 44 with subtype 1b. Two regions of the HCV genome were analysed in samples prior to combined therapy with alpha interferon plus ribavirin, one compressing the hypervariable regions (HVR1, HVR2 and HVR3) of the E2 glycoprotein and another one including the interferon-sensitive determining region (ISDR) and the V3 domain of the NS5A protein. Genetic diversity measures showed a clear tendency to higher genetic variability levels in nonresponder patients to antiviral treatment than in responder patients, although highly disperse values were present within each response group for both subtypes. A more detailed analysis of amino acid composition revealed the presence of several subtype-specific variants in a few positions, but no discriminating positions between responder and nonresponder patients were detected. Our results also revealed that most amino acid positions were highly conserved, especially for subtype 1a. We conclude that the outcome of the antiviral treatment might depend not only on the nature of one or a few independent positions, but more likely on the combination of several positions along the HCV genome. Moreover, the own host's ability to generate an appropriate systemic response, in combination with the action of antivirals, is also likely to be essential for treatment outcome.

Historical and biological determinants of genetic diversity in the highly endemic triploid sea lavender Limonium dufourii (Plumbaginaceae)

Microsatellite markers were used to evaluate the genetic diversity and population genetic structure in the critically endangered Limonium dufourii (Plumbaginaceae), a highly endemic triploid species from the coasts of eastern Spain. Sixty-five alleles from 13 microsatellite regions were amplified in a sample of 122 individuals collected from the six extant populations. Microsatellite patterns were consistent with the triploid nature of L. dufourii. Alleles were unambiguously assigned to two different parental subgenomes in this hybrid species and the greater contribution of the diploid parental subgenome was confirmed. Eleven, 25 and 26 multilocus genotypes were recorded from the haploid, diploid and from the combined information of both subgenomes, respectively. Genetic diversity was mostly distributed among populations (72.06% of the total genetic variation). Genotypes from Marjal del Moro populations grouped into two highly structured clusters (88.41% of the total variance). The observed patterns of distribution of genetic diversity are interpreted to result from multiple hybridization events and isolation between populations. Threats to this species are mainly anthropogenic (urbanization and tourism pressure), although stochastic risks cannot be ignored. Therefore, in order to preserve extant genetic variation of L. dufourii, in situ strategies such as the preservation of its habitat are a high priority. Several recommendations in order to assist ex situ measures to guarantee the success of conservation strategies and maintain the relationships between individuals and populations are proposed.

Evolutionary analysis of genetic variation observed in citrus tristeza virus (CTV) after host passage

We have studied the genetic variability in two genes (p18 and p20) from two groups of Citrus tristeza virus (CTV) isolates. One group (isolates T385, T317, T318, and T305) was derived from a Spanish source by successive host passages while the other (isolates T388 and T390) was obtained after aphid transmission from a Japanese source. A total of 274 sequences were obtained for gene p18 and 451 for p20. In the corresponding phylogenetic trees, sequences derived from the severe isolates (T318, T305, and T388) clustered together and separately from those derived from mild or moderate isolates (T385, T317, and T390), regardless of their geographic origin. Hierarchical analyses of molecular variance showed that up to 53% of the total genetic variability in p18 and up to 87% of the variation in p20 could be explained by differences in the pathogenicity features of the isolates. Neutrality tests revealed that different selection forces had been acting between isolates and between genes, with purifying selection being suggested for p18 from isolates T385 and T390 and for p20 from isolates T385, T317, and T388, and balancing selection for p18 from isolates T318, T305, and T388 and for p20 from isolates T318 and T390. Furthermore, several models of codon selection were observed, with purifying selection being the most notable one, compatible with low effective population size of the virus populations resulting from transmission bottlenecks. We found no evidence of recombination playing a significant role during p18 and p20 evolution in these isolates. These results suggest that hosts can be an important evolutionary factor for CTV isolates.

Genetic similarity of hepatitis C virus and fibrosis progression in chronic and recurrent infection after liver transplantation

The effect of hepatitis C virus (HCV) genetic heterogeneity on clinical features of post-transplantation hepatitis C is controversial. Different regions of the HCV genome have been associated with apoptosis, fibrosis, and other pathways leading to liver damage in chronic HCV infection. Besides, differences in immunodominant regions, such as NS3, may influence HCV-specific immune responses and disease outcome. In the liver transplant setting, a recent study has reported a positive association between HCV-1b Core region genetic relatedness 5-year post-transplantation and histological severity of recurrent hepatitis C. We have compared nucleotide sequences of HCV Core, NS3 and NS5b regions in HCV-1b-infected patients 3 years post-transplantation (n = 22). A cohort of nontransplanted patients (n = 22) was used as control of natural chronic HCV-1b infection. Histological evaluation was used to define the rate of fibrosis progression. Molecular variance analysis did not show significant differences in HCV sequences between transplanted and nontransplanted patients, or between those with fast or slow fibrosis progression. The same results were obtained when analysing phylogenetic trees for Core, NS3 and NS5b regions. A more appropriate clustering method (using minimum spanning networks) revealed a significant positive relationship between HCV genetic similarity in Core (r = 0.550, P < 0.01) and NS5b regions (r = 0.847, P < 0.01) and the yearly rate of fibrosis progression in nontransplanted patients which, in contrast, was not observed in transplanted patients. Our results indicate that some strains of HCV-1b might be more pathogenic in the natural course of chronic infection by this virus subtype. In the liver transplant setting, when the immune response is severely compromised, other mechanisms are probably more important in determining hepatitis C progression.

On the verge of extinction: genetics of the critically endangered Iberian plant species, Borderea chouardii (Dioscoreaceae) and implications for conservation management

Borderea chouardii is a relictual and dioecious, strictly sexually reproducing, long-living geophyte of the Dioscoreaceae family. Previous biological and demographic studies have indicated the existence of a uniformly distributed panmictic population of this taxon at the southernmost Spanish pre-Pyrenean mountain ranges where it occurs in rather inaccessible crevices of a single limestone cliff. However, individuals of B. chouardii are spatially subdivided into two subpopulations located, respectively, on the upper and lower parts of the cliff, and vertically separated 150 m. Because of its extreme rarity, B. chouardii was the first Iberian taxon to have a specific conservation plan and has been included in several red lists under the category of critically endangered (CR). However, no previous attempts have been conducted to analyse the fine scale evolutionary mechanisms involved in its present microspatial distribution. Genetic diversity and population structure have been investigated through the analysis of neutral hypervariable markers such as simple sequence repeats (SSRs) and randomly amplified polymorphic DNAs (RAPDs) to unravel the impact of life history traits in the differentiation of the two subpopulations. Both types of molecular markers were unequivocal in distinguishing two genetically distinct groups of individuals corresponding to their spatial separation. However, SSRs detected a higher level of subpopulation differentiation (F(ST) = 0.35, R(ST) = 0.32) than RAPDs (F(ST) = 0.21). SSR data indicated significant deviation from random dispersal of genes and genotypes between the two groups, suggesting that mating occurs mainly among individuals within subpopulations, thus, favouring the divergence between the two groups. This microevolutionary differentiation scenario might have been caused by a coupled effect of past genetic drift and reproductive isolation, as a result of strong glacial age bottlenecks and inefficient dispersal system of pollen and seeds, respectively. The identification of such genetic structure in this narrow endemic prompts a modification of the management strategies of its single extant population.

Characterization of seven (CTT)(n) microsatellite loci in the Pyrenean endemic Borderea pyrenaica (Dioscoreaceae): remarks on ploidy level and hybrid origin assessed through allozymes and microsatellite analyses

We present the identification and characterization of microsatellite loci in the Pyrenean endemic Borderea pyrenaica Miégeville (Dioscoreaceae). Seven microsatellite loci were isolated from a (CTT)(n)-enriched partial genomic library. Electropherograms patterns suggest that B. pyrenaica is a tetraploid species, as is its congener B. chouardii. One microsatellite locus was monomorphic, whereas the remaining ones presented from 2 to 10 alleles when analyzed in a sample of 60 individuals. Microsatellites have revealed higher levels of genetic variability than those in previous studies based on allozymes. Levels of genetic diversity are discussed in terms of tetrasomic (autotetraploidy) or duplicated disomic (allotetraploidy) modes of allele segregation. According to the first hypothesis, mean levels of genetic variability (H(min)-H(max)) range between 0.36 and 0.41, whereas, according to the second hypothesis, the 7 primer pairs amplified 11 chromosomal loci, and mean levels of observed and expected heterozygosities were 0.217 and 0.229, respectively, and did not differ significantly from HW expectations. These results suggest a hybrid allopolyploid origin for the Borderea taxa.

Detecting changes in the functional constraints of paralogous genes

We describe a new procedure to determine whether regional alterations in the evolutionary constraints imposed on paralogous proteins have occurred. We used as models the A and B (alternatively called alpha and beta) subunits of V/F/A-ATPases, originated by a gene duplication more than 3 billion years ago. Changes associated to three major splits (eubacteria versus Archaea-eukaryotes; Archaea versus eukaryotes; and among free-living bacteria and symbiotic mitochondria) were studied. Only in the first case, when we compared eubacterial or mitochondrial F-ATPases versus eukaryotic vacuolar V-ATPases or archaeal A-ATPases, constraint changes were observed. Modifications in the degree of regional constraining were not detected for the other two types of comparisons (V-ATPases versus A-ATPases and within F-ATPases, respectively). When the rates of evolution of the two subunits were compared, it was found that F-ATPases regulatory subunits evolved faster than catalytic subunits, but the opposite was true for A- and V-ATPases. Our results suggest that, even for universal and essential proteins, selective constraints may be occasionally altered. On the other hand, in some cases no changes were detected after periods of more than 2.2 billion years.

Ribotyping of vibrio populations associated with cultured oysters (Ostrea edulis)

The intraspecific variability of Vibrio splendidus, V. harveyi and V. tubiashii recovered from oysters (Ostrea edulis) collected at the Mediterranean coast near Valencia, Spain, was analyzed by ribotyping. The two former species represented the most abundant ones, and the third one was the only species described as pathogenic for oysters. A total of 115 environmental strains were studied, 84 of V. splendidus, 23 of V. harveyi and 8 of V. tubiashii. Chromosomal DNA was digested with KpnI and hybridized with an oligonucleotide probe complementary to a highly conserved sequence in the 23S rRNA gene. Ribotyping among natural populations of the three species rendered 5 to 9 bands, and showed a high genetic diversity, with a ratio no. of strains/no. of ribotypes between 1.1 and 1.5. Cluster analysis of V. splendidus ribotypes suggests a seasonal pattern of incidence, with those ribotypes corresponding to winter and spring samples being maintained in the oysters over the year.

Prephenate dehydratase from the aphid endosymbiont (Buchnera) displays changes in the regulatory domain that suggest its desensitization to inhibition by phenylalanine

Buchnera aphidicola, the prokaryotic endosymbiont of aphids, complements dietary deficiencies with the synthesis and provision of several essential amino acids. We have cloned and sequenced a region of the genome of B. aphidicola isolated from Acyrthosiphon pisum which includes the two-domain aroQ/pheA gene. This gene encodes the bifunctional chorismate mutase-prephenate dehydratase protein, which plays a central role in L-phenylalanine biosynthesis. Two changes involved in the overproduction of this amino acid have been detected. First, the absence of an attenuator region suggests a constitutive expression of this gene. Second, the regulatory domain of the Buchnera prephenate dehydratase shows changes in the ESRP sequence, which is involved in the allosteric binding of phenylalanine and is strongly conserved in prephenate dehydratase proteins from practically all known organisms. These changes suggest the desensitization of the enzyme to inhibition by phenylalanine and would permit the bacterial endosymbiont to overproduce phenylalanine.

Genetic differentiation and structure of Hippocrepis valentina (Leguminosae) populations

We present an analysis of isozyme variability in natural populations of the plant species Hippocrepis valentina (Leguminosae), which is endemic to the eastern Mediterranean coast of Spain and currently has endangered species status. Our results, obtained by starch-gel electrophoresis of 15 loci, show normal levels of variability for species with similar biology. The comparison with the patterns of genetic variability of two closely related species, H. balearica and H. grosii, confirms the taxonomic status of H. valentina as a proper species, independent of H. balearica, as previously suggested. The analysis of population subdivision shows that substantial variation among populations is present, and a hierarchical analysis demonstrates that when zones are defined according to their geographic location, a higher differentiation among populations within zones than among zones is found. Indirect estimates of gene flow indicate that levels of migration per generation are relatively low, except for a group composed of four populations, three of which are in close proximity. For the remaining populations, there is evidence of substantial differentiation. These results have implications for the design of a conservation strategy for this species.

Study of the evolutionary relationships among Limonium species (Plumbaginaceae) using nuclear and cytoplasmic molecular markers

The genus Limonium, due to the patchiness of the natural habitats of its species as well as the high frequency of hybridization and polyploidy and the possibility of reproduction by apomixis, provides an example of all the principal mechanisms of rapid speciation of plants. As an initial study of evolution in this genus, we have analyzed intra- and interspecific variability in 17 species from section Limonium, the largest in the genus, based on RFLPs of cpDNA and nuclear rDNA ITS sequences. In the cpDNA analysis, 21 restriction enzymes were used, resulting in 779 fragments, 490 of which were variable and 339 parsimony informative. L. furfuraceum exhibited two relatively divergent cpDNA haplotypes. The relationships found among the species based on cpDNA restriction fragments were coincident using different methods of phylogenetic analysis. Due to the presumed reticulate evolution in the genus Limonium, the comparison of these results with data from the nuclear DNA was necessary; ITS sequences were analyzed. The final alignment contained 488 characters, of which 198 were variable and 156 parsimony informative. Two relatively divergent ITS types were present at the intraindividual level in L. delicatulum, a triploid species. Each type was related to ITS from different groups of diploid Limonium species, one with a base haploid chromosome number n = 8 (represented by L. cossonianum) and the other with n = 9 (represented by L. minutum). The different phylogenetic inference methods used for the analysis of ITS sequences rendered very similar topologies. In general, the relationships among the species studied were coincident with those obtained with the chloroplast genome. Both nuclear and cytoplasmic markers support the polyphyly of section Limonium, with at least two species, L. narbonense and L. vulgare, clearly divergent from the rest. Moreover, the remaining subsections into which section Limonium is currently divided seem to be artificial.

An 18S rDNA-based molecular phylogeny of aphidiinae (Hymenoptera: braconidae)

We have obtained a molecular phylogeny of the subfamily Aphidiinae (Hymenoptera: Braconidae) by sequencing the 18S rDNA in 37 aphidiine taxa. Approximately 1857 nucleotides were sequenced in each species. Evolutionary relationships were established by comparing the results of maximum-parsimony, maximum-likelihood, and distance analyses. The most variable region of this gene, V4 (approx 403 nucleotides), was employed to establish the basality of the tribe Ephedrini within this subfamily. All phylogenetic reconstructions yielded trees with very similar topologies that confirmed the existence of two of the four traditionally accepted tribes, Ephedrini and Praini, but questioned the existence of Trioxini and Aphidiini. To better ascertain the status of some groups, the same analyses were repeated with a reduced taxonomic sample in which some species that produced systematic errors in the former phylogenetic reconstructions had been removed. The results from this second analysis favor either the three-tribes hypothesis (Ephedrini, Praini, and Aphidiini) or a new classification with at least five tribes (Ephedrini, Praini, Monoctonini, Trioxini, and Aphidiini). The 18S rDNA gene is a useful marker to recover relationships not only at the tribe but also at the subtribe and genus levels in this group. The natural status of some traditionally accepted clusters is also corroborated with the present data whereas the placement of other clusters is questioned or remains unresolved.

Analysis of population genetic structure and variability using RAPD markers in the endemic and endangered Limonium dufourii (Plumbaginaceae)

Limonium dufourii (Plumbaginaceae) is a triploid species, with apomictic reproduction, endemic to the east mediterranean coast of Spain, where it is present in only six populations with a few individuals in most of them. L. dufourii is included in the Red List of Endangered Species by the IUCN. Genetic variation and population structure in this species has been studied using RAPDs. Twelve different primers provided 124 reliable bands of which 33 were polymorphic among the 165 individuals analysed. Those polymorphic bands were able to define 44 different patterns, of which all but six were present in only one population. Several methods for statistical evaluation have been used for intra- and interpopulation analysis of genetic variability. Relationships among patterns have led to the identification of four main clusters. Two of them show a perfect correspondence to the population of origin of those individuals that present them (Cullera and Torreblanca), and the other two (Groups A and B) include patterns found in individuals coexisting in the same populations (Marjal del Moro populations) and in El Saler. Most of the variation found in this species is due to differences among populations as shown by the analysis of molecular variance. This agrees with the expectation for an apomictic species such as L. dufourii. The analysis of homogeneity of variance shows that substantial differences in the amount of genetic variability present in the six populations exist. These results have been used to understand the evolutionary and demographic history of L. dufourii, which is a requisite in order to establish efficient conservation measures for this species.

Evolution of fitness in experimental populations of vesicular stomatitis virus

The evolution of fitness in experimental clonal populations of vesicular stomatitis virus (VSV) has been compared under different genetic (fitness of initial clone) and demographic (population dynamics) regimes. In spite of the high genetic heterogeneity among replicates within experiments, there is a clear effect of population dynamics on the evolution of fitness. Those populations that went through strong periodic bottlenecks showed a decreased fitness in competition experiments with wild type. Conversely, mutant populations that were transferred under the dynamics of continuous population expansions increased their fitness when compared with the same wild type. The magnitude of the observed effect depended on the fitness of the original viral clone. Thus, high fitness clones showed a larger reduction in fitness than low fitness clones under dynamics with included periodic bottleneck. In contrast, the gain in fitness was larger the lower the initial fitness of the viral clone. The quantitative genetic analysis of the trait "fitness" in the resulting populations shows that genetic variation for the trait is positively correlated with the magnitude of the change in the same trait. The results are interpreted in terms of the operation of Muller's ratchet and genetic drift as opposed to the appearance of beneficial mutations.

A method for determining the position and size of optimal sequence regions for phylogenetic analysis

The availability of fast and accurate sequencing procedures along with the use of PCR has led to a proliferation of studies of variability at the molecular level in populations. Nevertheless, it is often impractical to examine long genomic stretches and a large number of individuals at the same time. In order to optimize this kind of study, we suggest a heuristic procedure for detection of the shortest region whose informational content can be considered sufficient for significant phylogenetic reconstruction. The method is based on the comparison of the pairwise genetic distances obtained from a set of sequences of reference to those obtained for different windows of variable size and position by means of a simple index. We also present an approach for testing whether the informative content in the stretches selected in this way is significantly different from the corresponding content shown by the larger genomic regions used as reference. Application of this test to the analysis of the VP1 protein gene of foot-and-mouth-disease type C virus allowed us to define optimal stretches whose informative content is not significantly different from that displayed by the complete VP1 sequence. We showed that the predictions made for type C sequences are valid for type O sequences, indicating that the results of the procedure are consistent.

Molecular characterization and cytonuclear disequilibria of two Drosophila subobscura mitochondrial haplotypes

According to restriction site analyses of mitochondrial DNA, Drosophila subobscura shows a polymorphism that consists of two frequent haplotypes that are evenly distributed all over the Old World and several rare haplotypes never present in more than one locality. To ascertain the causes responsible for such distribution, three different mtDNA fragments from haplotypes I and II sampled in a population from Zürich have been partially sequenced. Only three silent nucleotide changes have been detected in the ND5 gene. One of them implies the loss of the HaeIII restriction site, which differentiates haplotype I from haplotype II. On the basis of these results as well as on others involving the geographic distribution of haplotypes I and II, they can be considered phenotypically equivalent. The sequencing study has been complemented with the analysis of cytonuclear disequilibria between mitochondrial haplotypes and different nuclear loci in four D. subobscura populations. As expected, no significant cytonuclear disequilibria have been found between haplotypes I and II. Moreover, when haplotypes I and II were pooled and compared with the rare, endemic haplotypes similar results were obtained.

The effect of periodic bottlenecks on the competitive ability of Drosophila pseudoobscura lines

Competition experiments between several Drosophila pseudoobscura strains that were previously subjected to periodic bottlenecks of different sizes and D. willistoni have been carried out. Contrary to previous results with these two species, where stable coexistence was detected, populations of D. pseudoobscura displaced D. willistoni in a few generations. By using a relative fitness measure, the control lines (that had not gone through bottlenecks) outcompeted D. willistoni faster than bottlenecked lines. One of these, corresponding to the minimum possible bottleneck size of one pair, also shows significantly lower relative fitness than the other strains. These results are discussed in relation to the effect of periodic bottlenecks in reducing genetic variability and its incidence on the competitive ability of these lines.

Does the VP1 gene of foot-and-mouth disease virus behave as a molecular clock?

We have carried out a phylogenetic study of the evolution of the VP1 gene sequence from different serological types and subtypes of foot-and-mouth disease virus (FMDV). The maximum-likelihood method developed by Hasegawa and co-workers (Hasegawa et al. 1985) for the estimation of evolutionary parameters and branching dates has been used to decide between alternative models of evolution: constant versus variable rates. The results obtained indicate that a constant rate model, i.e., a molecular clock, seems to be the most plausible one. However, additional information suggests the possibility that the appearance of serotype CS has been accompanied by an episode of rapid evolution (Villaverde et al. 1991). We discuss the possibility that this evolution of RNA viruses was due to episodic positive Darwinian selection, which would have helped the new variant to escape the immunogenic pressure from the hosts.

Fitness alteration of foot-and-mouth disease virus mutants: measurement of adaptability of viral quasispecies

We document the rapid alteration of fitness of two foot-and-mouth disease virus (FMDV) mutants resistant to a neutralizing monoclonal antibody. Both mutants showed a selective disadvantage in BHK-21 cells when passaged in competition with their parental FMDV. Upon repeated replication of the mutants alone, they acquired a selective advantage over the parental FMDV and fixed additional genomic substitutions without reversion of the monoclonal antibody-resistant phenotype. Thus, variants that were previously kept at low frequency in the mutant spectrum of a viral quasispecies rapidly became the master sequence of a new genomic distribution and dominated the viral population.

On the analysis of viability data: an example with Drosophila

Larval competition experiments involving two wild type and eight mutant strains of Drosophila melanogaster have been carried out following the substitution procedure proposed by Mather and Caligari (1981). Our main goal has been to compare the competitive abilities of two phenotypically indistinguishable strains (wild and Oregon-R) by means of their responses with eight different mutants. Prior to the analyses of viability data, we have studied the normalizing effect of several transformations in order to determine which was best suited for the analyses. The differences found among the five transformations tested and the untransformed data were not very great. The folded power transformation (Mosteller and Tukey, 1977) was finally chosen. No constant pattern in the responses of the two wild type strains to the mutant competitors was detected. This leads us to conclude that the nature of the competition between the two wild type strains cannot be predicted from a knowledge of their competition with other strains.

Larval competition in Drosophila melanogaster: effects on development time

The effects on development time of larval competition that is interrupted at different times are analyzed in an experiment with constant larval density and varying proportions of two competing strains. Interruption of competition is carried out by means of the overfeeding technique, which allows for the study of individuals that, after several days under competition stress, either complete their development in crowded conditions or migrate to a vial with fresh food where the rest of their development takes place. Results of mean development time and within-vial variance of development time show that individuals under strong competition stress can arrest their development retaining the capability to resume it, once they are exposed to fresh food. The effects of competition are observed sooner on development time than on viabilities, which is in agreement with Ohba's observations (Biol. J. Okayama Univ. 7:87-125, 1961), but the observed effects are not always coincident for both fitness components.