150 Years of Coevolution Research: Evolution and Ecology of Yucca Moths (Prodoxidae) and Their Hosts

Yucca moths (Tegeticula and Parategeticula) are specialized pollinators of yucca plants, possessing unique, tentacle-like mouthparts used to actively collect pollen and deposit it onto the flowers of their hosts. The moths' larvae feed on the developing seeds and fruit tissue. First described in 1873, the yucca-yucca moth pollination system is now considered the archetypical example of a coevolved intimate mutualism. Research conducted over the past three decades has transformed our understanding of yucca moth diversity and host plant interactions. We summarize the current understanding of the diversity, ecology, and evolution of this group, review evidence for coevolution of the insects and their hosts, and describe how the nature of the interaction varies across evolutionary time and ecological contexts. Finally, we identify unresolved questions and areas for future research.

Genome-scale data resolves the timing of divergence in Joshua trees

PREMISE

Joshua trees (Yucca brevifolia and Y. jaegeriana) and their yucca moth pollinators (Tegeticula synthetica and T. antithetica) are a model system for studies of plant-pollinator coevolution and, they are thought to be one of the only cases in which there is compelling evidence for cospeciation driven by coevolution. Previous work attempted to evaluate whether divergence between the plant and their pollinators was contemporaneous. That work concluded that the trees diverged more than 5 million years ago-well before the pollinators. However, clear inferences were hampered by a lack of data from the nuclear genome and low genetic variation in chloroplast genes. As a result, divergence times in the trees could not be confidently estimated.

METHODS

We present an analysis of whole chloroplast genome sequence data and RADseq data from >5000 loci in the nuclear genome. We developed a molecular clock for the Asparagales and the Agavoideae using multiple fossil calibration points. Using Bayesian inference, we produced new estimates for the age of the genus Yucca and for Joshua trees. We used calculated summary statistics describing genetic variation and used coalescent-based methods to estimate population genetic parameters.

RESULTS

We find that the Joshua trees are moderately genetically differentiated, but that they diverged quite recently (~100-200 kya), and much more recently than their pollinators.

CONCLUSIONS

The results argue against the notion that coevolution directly contributed to speciation in this system, suggesting instead that coevolution with pollinators may have reinforced reproductive isolation following initial divergence in allopatry.

Strong Selection Against Early Generation Hybrids in Joshua Tree Hybrid Zone Not Explained by Pollinators Alone

Coevolution frequently plays an important role in diversification, but the role of obligate pollination mutualisms in the maintenance of hybrid zones has rarely been investigated. Like most members of the genus Yucca, the two species of Joshua tree (Yucca brevifolia and Yucca jaegeriana) are involved in a tightly coevolved mutualism with yucca moths. There is strong evidence of a history of coevolution between Joshua trees and their moth pollinators. We use a geographic clines approach in the Joshua tree hybrid zone to ask if selection by the moths may currently contribute to maintaining separation between these species. We compare genomic, phenotypic, and pollinator frequency clines to test whether pollinators maintain the hybrid zone or follow it as passive participants. The results reveal dramatic overlapping genomic and pollinator clines, consistent with a narrow hybrid zone maintained by strong selection. Wider phenotypic clines and a chloroplast genomic cline displaced opposite the expected direction suggest that pollinators are not the main source of selection maintaining the hybrid zone. Rather, it seems that high levels of reproductive isolation, likely acting through multiple barriers and involving many parts of the genome, keep the hybrid zone narrow.

Mitochondrial Divergence between Western and Eastern Great Bustards: Implications for Conservation and Species Status

The great bustard is the heaviest bird capable of flight and an iconic species of the Eurasian steppe. Populations of both currently recognized subspecies are highly fragmented and critically small in Asia. We used DNA sequence data from the mitochondrial cytochrome b gene and the mitochondrial control region to estimate the degree of mitochondrial differentiation and rates of female gene flow between the subspecies. We obtained genetic samples from 51 individuals of Otis tarda dybowskii representing multiple populations, including the first samples from Kazakhstan and Mongolia and samples from near the Altai Mountains, the proposed geographic divide between the subspecies, allowing for better characterization of the boundary between the 2 subspecies. We compared these with existing sequence data (n = 66) from Otis tarda tarda. Our results suggest, though do not conclusively prove, that O. t. dybowskii and O. t. tarda may be distinct species. The geographic distribution of haplotypes, phylogenetic analysis, analyses of molecular variance, and coalescent estimation of divergence time and female migration rates indicate that O. t. tarda and O. t. dybowskii are highly differentiated in the mitochondrial genome, have been isolated for approximately 1.4 million years, and exchange much less than 1 female migrant per generation. Our findings indicate that the 2 forms should at least be recognized and managed as separate evolutionary units. Populations in Xinjiang, China and Khövsgöl and Bulgan, Mongolia exhibited the highest levels of genetic diversity and should be prioritized in conservation planning.

First Recorded Observations of Pollination and Oviposition Behavior in Tegeticula antithetica (Lepidoptera: Prodoxidae) Suggest a Functional Basis for Coevolution With Joshua Tree (Yucca) Hosts

Yucca moths (Tegeticula spp.) are the exclusive pollinators of Joshua trees (Yucca brevifolia s. l.). The moths actively pollinate the Joshua tree flowers and lay their eggs in the style. Recent studies have revealed that the plants commonly known as Joshua trees include two distinct, sister-species of plant: Yucca brevifolia Engelm. and Yucca jaegeriana McKelvey, each pollinated by two sister-species of yucca moth Tegeticula synthetica Riley and Tegeticula antithetica Pellmyr, respectively. A number of studies have argued that the moths have coevolved with their hosts, producing a pattern of phenotype matching between moth ovipositor length and floral style length. However, the only known descriptions of yucca moth pollination and oviposition behavior on Joshua trees are observations of T. synthetica made in 1893. The behavior of T. antithetica has never been observed before. We produced the first video recordings of the behavior of T. antithetica, and measured the points of oviposition and egg placement within the floral style. We found a number of differences between the behaviors of T. antithetica and T. synthetica, which appear to be a consequence of differences in floral morphology between Y. jaegeriana and Y. brevifolia. We also found that variation in floral style length strongly influences the placement of eggs within the flower, which may explain patterns of phenotype matching described previously. However, unlike in other yucca moths, we find that the mode of oviposition is unlikely to wound the floral ovules, and thus that oviposition by T. antithetica is unlikely to prompt floral abscission.

Population genomics of divergence within an obligate pollination mutualism: Selection maintains differences between Joshua tree species

PREMISE OF THE STUDY

Speciation is a complex process that can be shaped by many factors, from geographic isolation to interspecific interactions. In Joshua trees, selection from pollinators on style length has been hypothesized to contribute to the maintenance of differentiation between two hybridizing sister species. We used population genomics approaches to measure the extent of genetic differentiation between these species, test whether selection maintains differences between them, and determine whether genetic variants associated with style length show signatures of selection.

METHODS

Using restriction-site-associated DNA (RAD)-sequencing, we identified 9516 single nucleotide polymorphisms (SNPs) across the Joshua tree genome. We characterized the genomic composition of trees in a narrow hybrid zone and used genomic scans to search for signatures of selection acting on these SNPs. We used a genome-wide association study to identify SNPs associated with variation in phenotypic traits, including style length, and asked whether those SNPs were overrepresented among the group under selection.

KEY RESULTS

The two species were highly genetically differentiated (F_ST = 0.25), and hybrids were relatively rare in the hybrid zone. Approximately 20% of SNPs showed evidence of selection maintaining divergence. While SNPs associated with style length were overrepresented among those under selection (P << 0.0001), the same was true for SNPs associated with highly differentiated vegetative traits.

CONCLUSIONS

The two species of Joshua tree are clearly genetically distinct, and selection is maintaining differences between them. We found that loci associated with differentiated traits were likely to be under selection. However, many traits other than style length appeared to be under selection. Together with the dearth of intermediate hybrids, these findings reveal that these taxa are more strongly diverged than previously suspected and that selection, likely on many targets, is maintaining separation where the two species meet and hybridize.

Anion replacement at Au(110)/electrolyte interfaces

A characteristic reflection anisotropy spectrum (RAS) is observed from a Au(110) surface in a wide range of electrolytes and combinations of pH and applied potentials. It is suggested that this common RAS profile arises from an interaction between the potential applied to the Au(110) electrode and the dipole moments of oxidized species that locates the Fermi level at a common position with respect to the electronic band structure of Au. Rapid changes in this RAS profile are observed for Au(110)/H2SO4 as the potential is switched between 0.3 V and 0.6 V, a potential range in which the surface is not reconstructed and below the potential range of surface oxidation. The spectral changes are completed in less than 10 ms, are reversible and are attributed to the replacement of adsorbed anions by an oxygenated species.

The stability of the Au(1 1 0)-(1 × 3) surface reconstruction in electrochemical environments

Changes in the reflection anisotropy (RAS) profile of the Au(1 1 0)-(1  ×  3)/Na2SO4 interface over 25 h are attributed to the slow accumulation of impurities on the Au(1 1 0) surface which reduce the intensity of a transition involving a surface state that makes a positive contribution to the RAS profile at 1.8 eV. The growth in the intensity of a feature that makes a negative contribution to the RAS profile at 2.6 eV and the reduction in the intensity of contributions to higher energy is attributed to shifts in the energy of the surface band structure relative to the Fermi level caused by the accumulation of impurities. There is no clear explanation of the subsequent decay of the 2.6 eV feature or the long term reduction in intensity to high energy of the RAS profile.

The reflection anisotropy spectroscopy of the Au(1 1 0) surface structures in liquid environments

The reflection anisotropy (RAS) profiles of the Au(1 1 0)-(1  ×  1), (1  ×  2) and (1  ×  3) surface structures in electrochemical environments are shown to arise mainly from surface dipole transitions directed along the principal axes of the Au(1 1 0) surface. There are weak contributions to the RAS profiles of the Au(1 1 0)-(1  ×  1) and (1  ×  3) surfaces in the region of 4.0 eV which probably arise from (1 1 1) facets that are either intrinsic to the surface structures or are associated with steps. A transition involving a surface state just above the Fermi level, E F, contributes to the RAS profiles of the (1  ×  2) and (1  ×  3) surfaces but not to the RAS profile of the (1  ×  1) surface. A strong feature at 2.5 eV in the RAS profiles of the Au(1 1 0)-(1  ×  1) and (1  ×  2) surfaces is attributed to a transition in the vicinity of the L point of the Brillouin zone between the 5d band and the [Formula: see text] band at E F. It is argued that the applied potential of  -0.6 V, which creates the Au(1 1 0)-(1  ×  3) surface, lifts E F above the [Formula: see text] band causing it to become occupied and quenching this contribution to the RAS profile.

Conformational change in cytochrome P450 reductase adsorbed at a Au(110)-phosphate buffer interface induced by interaction with nicotinamide adenine dinucleotide phosphate

Changes observed in the reflection anisotropy spectroscopy (RAS) profiles of monolayers of cytochrome P450 reductase adsorbed at Au(110)-electrolyte interfaces at 0.056 V following the addition of nicotinamide adenine dinucleotide phosphate (NADP(+)) are explained in terms of a simple model as arising from changes in the orientation of an isoalloxazine ring located in the flavin mononucleotide binding domain of the protein. The model also accounts for the changes observed in the RAS as the potential applied to the Au(110) surface is varied and suggests that differences in the dependence of the RAS profile of the adsorbed protein on the potential applied to the electrode in the absence and presence of NADP(+) are explicable as arising from a competition between the applied potential acting to reduce the protein and the NADP(+) to oxidize it.

Conformational change induced by electron transfer in a monolayer of cytochrome P450 reductase adsorbed at the Au(110)-phosphate buffer interface

The reflection anisotropy spectroscopy profiles of a variant of cytochrome P450 reductase adsorbed at the Au(110)-phosphate buffer interface depend on the sequence of potentials applied to the Au(110) electrode. It is suggested that this dependence arises from changes in the orientation of the isoalloxazine ring structures in the protein with respect to the Au(110) surface. This offers a method of monitoring conformational change in this protein by measuring variations in the reflection anisotropy spectrum arising from changes in the redox potential.

Long-term clearance of hepatitis C virus following interferon α-2b or peginterferon α-2b, alone or in combination with ribavirin

Sustained virologic response (SVR) is the standard measure for evaluating response to therapy in patients with chronic hepatitis C (CHC). The aim of this study was to prospectively assess the durability of SVR in the pivotal studies of peginterferon (PEG-IFN) α-2b or IFN α-2b. We conducted two phase 3b long-term follow-up studies of patients previously treated for CHC in eight prospective randomized studies of IFN α-2b and/or PEG-IFN α-2b. Patients who achieved SVR [undetectable hepatitis C virus (HCV) RNA 24 weeks after completion of treatment] were eligible for inclusion in these follow-up studies. In total, 636 patients with SVR following treatment with IFN α-2b and 366 with SVR following treatment with PEG-IFN α-2b were enrolled. Definite relapse (quantifiable serum HCV RNA with no subsequent undetectable HCV RNA) was reported in six patients treated with IFN α-2b and three patients treated with PEG-IFN α-2b. Based on these relapses, the point estimate for the likelihood of maintaining response after 5 years was 99.2% [95% confidence interval (CI), 98.1-99.7%] for IFN α-2b and 99.4% (95% CI, 97.7-99.9%) for PEG-IFN α-2b. Successful treatment of hepatitis C with PEG-IFN α-2b or IFN α-2b leads to clinical cure of hepatitis C in the vast majority of cases.

The nature and stability of the Au(110)/electrochemical interface produced by flame annealing

It is shown using reflection anisotropy spectroscopy (RAS) that following flame annealing and immersion in pure water, the Au(110) surface adopts a (1×1) structure and that this structure is preserved in a 0.1 M H(2)SO(4) environment. The surface transforms to the (1×2) reconstruction following the application of a potential of 0.0 V versus SCE (a saturated calomel electrode). This surface is unstable and the RAS profile changes over periods of 15 min and 1 h in a manner which suggests that changes are occurring in the structure and distribution of [11(-)0] steps. Over longer periods the RAS transforms towards a profile attributed to a surface associated with the specific adsorption of anions.

Controlling the formation of a monolayer of cytochrome P450 reductase onto Au surfaces

The conditions necessary for the formation of a monolayer and a bilayer of a mutated form (P499C) of human cytochrome P450 reductase on a Au(110)/electrolyte interface have been determined using a quartz crystal microbalance with dissipation, atomic force microscopy, and reflection anisotropy spectroscopy (RAS). The molecules adsorb through a Au-S linkage and, for the monolayer, adopt an ordered structure on the Au(110) substrate in which the optical axes of the dipoles contributing to the RAS signal are aligned roughly along the optical axes of the Au(110) substrate. Differences between the absorption spectrum of the molecules in a solution and the RAS profile of the adsorbed monolayer are attributed to surface order in the orientation of dipoles that contribute in the low energy region of the spectrum, a roughly vertical orientation on the surface of the long axes of the isoalloxazine rings and the lack of any preferred orientation in the molecular structure of the dipoles in the aromatic amino acids. Our studies establish an important proof of principle for immobilizing large biological macromolecules to gold surfaces. This opens up detailed studies of the dynamics of biological macromolecules by RAS, which have general applications in studies of biological redox chemistry that are coupled to protein dynamics.

Comparative phylogeography of a coevolved community: concerted population expansions in Joshua trees and four yucca moths

Comparative phylogeographic studies have had mixed success in identifying common phylogeographic patterns among co-distributed organisms. Whereas some have found broadly similar patterns across a diverse array of taxa, others have found that the histories of different species are more idiosyncratic than congruent. The variation in the results of comparative phylogeographic studies could indicate that the extent to which sympatrically-distributed organisms share common biogeographic histories varies depending on the strength and specificity of ecological interactions between them. To test this hypothesis, we examined demographic and phylogeographic patterns in a highly specialized, coevolved community--Joshua trees (Yucca brevifolia) and their associated yucca moths. This tightly-integrated, mutually interdependent community is known to have experienced significant range changes at the end of the last glacial period, so there is a strong a priori expectation that these organisms will show common signatures of demographic and distributional changes over time. Using a database of >5000 GPS records for Joshua trees, and multi-locus DNA sequence data from the Joshua tree and four species of yucca moth, we combined paleaodistribution modeling with coalescent-based analyses of demographic and phylgeographic history. We extensively evaluated the power of our methods to infer past population size and distributional changes by evaluating the effect of different inference procedures on our results, comparing our palaeodistribution models to Pleistocene-aged packrat midden records, and simulating DNA sequence data under a variety of alternative demographic histories. Together the results indicate that these organisms have shared a common history of population expansion, and that these expansions were broadly coincident in time. However, contrary to our expectations, none of our analyses indicated significant range or population size reductions at the end of the last glacial period, and the inferred demographic changes substantially predate Holocene climate changes.

Spectral signatures of the surface reconstructions of Au(110)/electrolyte interfaces

It is demonstrated that the (1 × 1) structure and the (1 × 2) and (1 × 3) surface reconstructions that occur at Au(110)/electrolyte interfaces have unique optical fingerprints. The optical fingerprints are potential, pH and anion dependent and have potential for use in monitoring dynamic changes at this interface. We also observe a specific reflection anisotropy spectroscopy signature that may arise from anions adsorbed on the (1 × 1) structure of Au(110).

How to become a yucca moth: Minimal trait evolution needed to establish the obligate pollination mutualism

The origins of obligate pollination mutualisms, such as the classic yucca-yucca moth association, appear to require extensive trait evolution and specialization. To understand the extent to which traits truly evolved as part of establishing the mutualistic relationship, rather than being preadaptations, we used an expanded phylogenetic estimate with improved sampling of deeply-diverged groups to perform the first formal reconstruction of trait evolution in pollinating yucca moths and their non-pollinating relatives. Our analysis demonstrates that key life history traits of yucca moths, including larval feeding in the floral ovary and the associated specialized cutting ovipositor, as well as colonization of woody monocots in xeric habitats, may have been established before the obligate mutualism with yuccas. Given these preexisting traits, novel traits in the mutualist moths are limited to the active pollination behaviors and the tentacular appendages that facilitate pollen collection and deposition. These results suggest that a highly specialized obligate mutualism was built on the foundation of preexisting interactions between early Prodoxidae and their host plants, and arose with minimal trait evolution.

Divergence in an obligate mutualism is not explained by divergent climatic factors

Adaptation to divergent environments creates and maintains biological diversity, but we know little about the importance of different agents of ecological divergence. Coevolution in obligate mutualisms has been hypothesized to drive divergence, but this contention has rarely been tested against alternative ecological explanations. Here, we use a well-established example of coevolution in an obligate pollination mutualism, Yucca brevifolia and its two pollinating yucca moths, to test the hypothesis that divergence in this system is the result of mutualists adapting to different abiotic environments as opposed to coevolution between mutualists. We used a combination of principal component analyses and ecological niche modeling to determine whether varieties of Y. brevifolia associated with different pollinators specialize on different environments. Yucca brevifolia occupies a diverse range of climates. When the two varieties can disperse to similar environments, they occupy similar habitats. This suggests that the two varieties have not specialized on distinct habitats. In turn, this suggests that nonclimatic factors, such as the biotic interaction between Y. brevifolia and its pollinators, are responsible for evolutionary divergence in this system.

Pattern and timing of diversification in Yucca (Agavaceae): specialized pollination does not escalate rates of diversification

The yucca-yucca moth interaction is one of the most well-known and remarkable obligate pollination mutualisms, and is an important study system for understanding coevolution. Previous research suggests that specialist pollinators can promote rapid diversification in plants, and theoretical work has predicted that obligate pollination mutualism promotes cospeciation between plants and their pollinators, resulting in contemporaneous, parallel diversification. However, a lack of information about the age of Yucca has impeded efforts to test these hypotheses. We used analyses of 4322 AFLP markers and cpDNA sequence data representing six non-protein-coding regions (trnT-trnL, trnL, trnL intron, trnL-trnF, rps16 and clpP intron 2) from all 34 species to recover a consensus organismal phylogeny, and used penalized likelihood to estimate divergence times and speciation rates in Yucca. The results indicate that the pollination mutualism did not accelerate diversification, as Yucca diversity (34 species) is not significantly greater than that of its non-moth-pollinated sister group, Agave sensu latissimus (240 species). The new phylogenetic estimates also corroborate the suggestion that the plant-moth pollination mutualism has at least two origins within the Agavaceae. Finally, age estimates show significant discord between the age of Yucca (ca 6-10Myr) and the current best estimates for the age of their pollinators (32-40Myr).

Historical Biogeography: The New Synthesis

Since Darwin first noted the impact of past ice ages on the distribution of organisms, biogeographers have debated whether Pleistocene glaciations shaped evolutionary patterns. A new synthesis of population genetics and palaeoclimatology promises unprecedented insights into Pleistocene evolutionary history.

Dos and don'ts of testing the geographic mosaic theory of coevolution

The geographic mosaic theory of coevolution is stimulating much new research on interspecific interactions. We provide a guide to the fundamental components of the theory, its processes and main predictions. Our primary objectives are to clarify misconceptions regarding the geographic mosaic theory of coevolution and to describe how empiricists can test the theory rigorously. In particular, we explain why confirming the three main predicted empirical patterns (spatial variation in traits mediating interactions among species, trait mismatching among interacting species and few species-level coevolved traits) does not provide unequivocal support for the theory. We suggest that strong empirical tests of the geographic mosaic theory of coevolution should focus on its underlying processes: coevolutionary hot and cold spots, selection mosaics and trait remixing. We describe these processes and discuss potential ways each can be tested.

Orientation of ordered structures of cytosine and cytidine -monophosphate adsorbed at Au(110)/liquid interfaces

It is demonstrated using reflection anisotropy spectroscopy that the adsorption of cytosine and cytidine -monophosphate at the Au(110) 1 x 2/electrolyte interface gives rise to ordered structures in which the base is oriented vertical to the surface and parallel to the [110] axis of the Au(110) plane.

Evolutionary Consequences of Dispersal Ability in Cactus-feeding Insects

Although gene flow is an important determinant of evolutionary change, the role of ecological factors such as specialization in determining migration and gene flow has rarely been explored empirically. To examine the consequences of dispersal ability and habitat patchiness on gene flow, migration rates were compared in three cactophagous longhorn beetles using coalescent analyses of mtDNA sequences. Analyses of covariance were used to identify the roles of dispersal ability and habitat distribution in determining migration patterns. Dispersal ability was a highly significant predictor of gene flow (p< 0.001), and was more important than any other factor. These findings predict that dispersal ability may be an import factor shaping both microevolutionary and macroevolutionary patterns; this prediction is borne out by comparisons of species diversity in cactus-feeding groups.

Phylogeography of the longhorn cactus beetle Moneilema appressum LeConte (Coleoptera: Cerambycidae): was the differentiation of the Madrean sky islands driven by Pleistocene climate changes?

Although it has been suggested that Pleistocene climate changes drove population differentiation and speciation in many groups of organisms, population genetic evidence in support of this scenario has been ambiguous, and it has often been difficult to distinguish putative vicariance from simple isolation by distance. The sky island communities of the American Southwest present an ideal system in which to compare late Pleistocene range fragmentations documented by palaeoenvironmental studies with population genetic data from organisms within these communities. In order to elucidate the impact of Pleistocene climate fluctuations on these environments, biogeographic patterns in the flightless longhorn cactus beetle, Moneilema appressum were examined using mitochondrial DNA sequence data. Gene tree relationships between haplotypes were inferred using parsimony, maximum-likelihood, and Bayesian analysis. Nested clade analysis, Mantel tests, and coalescent modelling were employed to examine alternative biogeographic scenarios, and to test the hypothesis that Pleistocene climate changes drove population differentiation in this species. The program mdiv was used to estimate migration and divergence times between populations, and to measure the statistical support for isolation over ongoing migration. These analyses showed significant geographic structure in genetic relationships, and implicated topography as a key determinant of isolation. However, although the coalescent analyses suggested that a history of past habitat fragmentation underlies the observed geographic patterns, the nested clade analysis indicated that the pattern was consistent with isolation by distance. Estimated divergence times indicated that range fragmentation in M. appressum is considerably older than the end of the most recent glacial, but coincided with earlier interglacial warming events and with documented range expansions in other, desert-dwelling species of Moneilema.

Range expansions in the flightless longhorn cactus beetles, Moneilema gigas and Moneilema armatum, in response to Pleistocene climate changes

Pollen cores and plant and animal fossils suggest that global climate changes at the end of the last glacial period caused range expansions in organisms indigenous to the North American desert regions, but this suggestion has rarely been investigated from a population genetic perspective. In order to investigate the impact of Pleistocene climate changes and glacial/interglacial cycling on the distribution and population structure of animals in North American desert communities, biogeographical patterns in the flightless, warm-desert cactus beetles, Moneilema gigas and Moneilema armatum, were examined using mitochondrial DNA (mtDNA) sequence data from the cytochrome oxidase I (COI) gene. Gene tree relationships between haplotypes were inferred using parsimony, maximum-likelihood, and Bayesian analysis. Nested clade analysis and coalescent modelling using the programs mdiv and fluctuate were used to identify demographically independent populations, and to test the hypothesis that Pleistocene climate changes caused recent range expansions in these species. A sign test was used to evaluate the probability of observing concerted population growth across multiple, independent populations. The phylogeographical and nested clade analyses reveal a history of northward expansion in both of these species, as well as a history of past range fragmentation, followed by expansion from refugia. The coalescent analyses provide highly significant evidence for independent range expansions from multiple refugia, but also identify biogeographical patterns that predate the most recent glacial period. The results indicate that widespread desert environments are more ancient than has been suggested in the past.

Phosphorylation of Bruton's tyrosine kinase by c-Abl

Bruton's tyrosine kinase (Btk) is necessary for B-lymphocyte development. Mutation in the gene coding for Btk causes X-linked agammaglobulinemia (XLA) in humans. Similar to Btk, c-Abl is a tyrosine kinase shuttling between the cytoplasm and the nucleus where it is involved in different functions depending on the localization. In this report we describe for the first time that c-Abl and Btk physically interact and that c-Abl can phosphorylate tyrosine 223 in the SH3 domain of Btk. Interestingly, the Btk sequence matched a v-Abl substrate [correction] identified from a randomized peptide library and was also highly related to a number of previously found c-Abl substrates.

Both proline-rich sequences in the TH region of Bruton's tyrosine kinase stabilize intermolecular interactions with the SH3 domain

The Tec homology (TH) region located N-terminal to the Src homology 3 (SH3) domain of Bruton's tyrosine kinase (Btk) contains two proline-rich SH3-binding sequences (PRRs). We have previously demonstrated that the TH region acts to stabilize intermolecular interactions in N-terminally extended SH3 (PRR-SH3) fragments. Here, we analyze six PRR-SH3 fragments with different proline-to-alanine substitutions in the two PRRs. Gel permeation chromatography and nuclear magnetic resonance spectroscopy show that both PRRs can stabilize self-association. This observation provides an explanation to why the TH region of Btk makes intermolecular interactions, whereas the corresponding interaction in the related Itk kinase with only one PRR, is intramolecular.

Variations in the human phospholipase Cgamma2 gene in patients with B-cell defects of unknown etiology

Our recent studies using targeted gene disruption have shown that defects in phospholipase Cgamma2 (PLCgamma2) result in a B-cell abnormality that is very similar to that seen in Btk-deficient mice. Null mutations in either PLCG2 or BTK are associated with decreased numbers of mature B cells, failure to make antibodies to some T cell-independent antigens and the absence of CD5+ peritoneal B cells. Mutations in BTK in humans cause a more severe defect in B-cell development characterized by almost complete absence of B cells in the peripheral circulation, profound hypogammaglobulinemia and an inability to produce antibodies to any antigens. However, not all patients with severe defects in B-cell development have mutations in BTK or the components of the B-cell signal transduction complex. To explore the possibility that some patients with defects in B-cell development of unknown etiology might have mutations in PLCG2, we determined the genomic structure of this gene and established conditions to analyze the 32 exons of the gene and the flanking sequences by single-strand conformation polymorphism. Although 24 polymorphic variants of this gene were found in 35 patients, we did not identify any alterations that were likely to be the cause of disease.

A mutation in Bruton's tyrosine kinase as a cause of selective anti-polysaccharide antibody deficiency

Children and adults can have recurrent infection with invasive encapsulated bacterial pathogens such as Streptococcus pneumoniae as a result of a selective inability to respond to polysaccharide antigens. We have identified a mutation in the gene encoding Bruton's tyrosine kinase in a male patient with selective anti-polysaccharide antibody deficiency.

Microphthalmia, facial anomalies, microcephaly, thumb and hallux hypoplasia, and agammaglobulinemia

We report a boy, born to consanguineous patients, with agammaglobulinemia associated with multiple physical anomalies: mild intrauterine growth retardation, extreme microphthalmia (clinical anophthalmia), severe microcephaly, blepharophimosis, long face with temporal narrowing, scaphocephalic skull shape, posterior cleft palate, hypoplastic, adducted thumbs with small nails, and short, inward turned halluces with absent distal phalanges and nails. Psychomotor development was moderately delayed. No mutations were found in exons of BKT and PAX-5 genes. This unreported constellation could represent a novel, autosomal recessive syndrome.

Mapping of herpes simplex virus-1 VP22 functional domains for inter- and subcellular protein targeting

The herpes simplex virus 1 (HSV-1) tegument protein VP22 has been utilised as a vehicle for trafficking proteins. It has a remarkable property of exiting the cell that is producing it and entering the neighbouring cells, which has been used to deliver therapeutic proteins, p53 and herpes simplex virus thymidine kinase (tk). It has a complex pattern of expression and subcellular localisation. Functions of VP22 include intercellular transport, binding to and bundling of microfilaments, inducing cytoskeleton collapse, nuclear translocation during mitosis, and binding to chromatin and nuclear membrane. The regions of VP22 which contain each of these functions have not been characterised. Finding the region carrying the property of intercellular spread would facilitate enhancement of transport function. By constructing a series of deletion constructs of VP22 tagged by the green fluorescent protein (GFP) we have mapped the functions of VP22 to specific regions in the polypeptide as follows: intercellular transport - aa 81-195; binding and reorganisation of cytoskeleton - aa 159-267; nuclear targeting, inhibition of cytoskeleton collapse - aa 81-121; and nuclear targeting and facilitation of intercellular transport - aa 267-301. Separation of VP22 functions enables focus on the mechanism of VP22-mediated transport and improve the transportation efficiency of VP22.

Expression and release of stable and active forms of murine granulocyte-macrophage colony-stimulating factor (mGM-CSF) targeted to different subcellular compartments

Cytokines have been used for several years as immunomodulators. However, one of the main drawbacks of systemically applied cytokines is their high toxicity. In addition, cytokines work in a paracrine form and frequently after cell-to-cell interaction. Therefore, a very restricted release of cytokines-in time and space-could be desired for most of their therapeutic applications. The murine granulocyte-macrophage colony-stimulating factor (mGM-CSF) is one of the most promising cytokine candidates for cancer immunotherapy and as an adjuvant of DNA vaccines. With the aim of improving the administration and release of cytokines in a very restricted area, we have designed vectors expressing the mGM-CSF cDNA with different localization signals. Using this strategy we have shown that cytokines can be expressed and targeted to different subcellular compartments (i.e. the cytoplasm, the endoplasmic reticulum and the nucleus), stored inside the cells and released after cell lysis as stable active proteins. Moreover, a plasma membrane targeted form of mGM-CSF displayed substantial amount of biological activity. These vectors could have potential applications in immunotherapy for tumours and DNA vaccination protocols.

The Tec family of cytoplasmic tyrosine kinases: mammalian Btk, Bmx, Itk, Tec, Txk and homologs in other species

Cytoplasmic protein-tyrosine kinases (PTKs) are enzymes involved in transducing a vast number of signals in metazoans. The importance of the Tec family of kinases was immediately recognized when, in 1993, mutations in the gene encoding Bruton's tyrosine kinase (Btk) were reported to cause the human disease X-linked agammaglobulinemia (XLA). Since then, additional kinases belonging to this family have been isolated, and the availability of full genome sequences allows identification of all members in selected species enabling phylogenetic considerations. Tec kinases are endowed with Pleckstrin homology (PH) and Tec homology (TH) domains and are involved in diverse biological processes related to the control of survival and differentiation fate. Membrane translocation resulting in the activation of Tec kinases with subsequent Ca2+ release seems to be a general feature. However, nuclear translocation may also be of importance. The purpose of this essay is to characterize members of the Tec family and discuss their involvement in signaling. The three-dimensional structure, expression pattern and evolutionary aspects will also be considered.

Proximal promoter of the murine syndecan-1 gene is not sufficient for the developmental pattern of syndecan expression in B lineage cells

Syndecan-1 (CD138) is a cell membrane proteoglycan that binds extracellular matrix components and various growth factors. The role of syndecan-1 in the control of cell growth and morphology has been illustrated by its altered expression in hematological malignancies such as multiple myeloma as well as some solid tumors. It has been reported that the expression of syndecan-1 in cells of the B lineage is developmentally regulated such that pre-B cells and plasma cells express syndecan-1 while mature B cells do not. Thus, we investigated whether the proximal promoter region of the murine syndecan-1 promoter was able to confer the observed on-off-on expression of syndecan-1 in cells of the B lineage as they develop from pre-B cells to plasma cells. Experiments carried out using deletion mutants of the proximal promoter cloned upstream of the CAT reporter gene transfected into murine cell lines, representing the above stages of B-cell development, such as BA/F3 (pro-B cell), 70Z/3 (pre-B cell), 2PK3 (late mature B cell), and MPC-11 (plasma cell), showed detectable levels of CAT expression. The WEHI-231 (mature B cell) cell lines did not show detectable levels of CAT reporter activity. The strong levels of expression were observed with a fragment of the proximal promoter spanning the region from -365 to -95 (from the translation start point). However, Northern analysis of RNA obtained from the five murine B-cell lines, representing various stages of B-cell development, showed that the 70Z/3, MPC-11 but not BA/F3, and 2PK3 cells expressed detectable levels of syndecan-1 mRNA. By FACS analysis, using a rat anti mouse syndecan-1 antibody, syndecan-1 expression on the cell surface was found to correlate with the observed mRNA expression patterns in these cell lines. Our results indicate that the proximal promoter of the murine syndecan-1 promoter is not sufficient for the observed developmental pattern of syndecan expression in B cells.

High levels of interferon gamma in the plasma of children with complete interferon gamma receptor deficiency

We have found that children with complete interferon gamma (IFNgamma) receptor deficiency, unlike patients with other genetic defects predisposing them to mycobacterial diseases, have very high levels of IFNgamma in their plasma. This unexpected observation provides a simple and accurate diagnostic method for complete IFNgamma receptor deficiency in children with clinical disease caused by bacille Calmette-Guérin vaccines or environmental nontuberculous mycobacteria.

Primary immunodeficiency mutation databases

Primary immunodeficiencies are intrinsic defects of immune systems. Mutations in a large number of cellular functions can lead to impaired immune responses. More than 80 primary immunodeficiencies are known to date. During the last years genes for several of these disorders have been identified. Here, mutation information for 23 genes affected in 14 immunodefects is presented. The proteins produced are employed in widely diverse functions, such as signal transduction, cell surface receptors, nucleotide metabolism, gene diversification, transcription factors, and phagocytosis. Altogether, the genetic defect of 2,140 families has been determined. Diseases with X-chromosomal origin constitute about 70% of all the cases, presumably due to full penetrance and because the single affected allele causes the phenotype. All types of mutations have been identified; missense mutations are the most common mutation type, and truncation is the most common effect on the protein level. Mutational hotspots in many disorders appear in CPG dinucleotides. The mutation data for the majority of diseases are distributed on the Internet with a special database management system, MUTbase. Despite large numbers of mutations, it has not been possible to make genotype-phenotype correlations for many of the diseases.

Intermolecular interactions between the SH3 domain and the proline-rich TH region of Bruton's tyrosine kinase

The SH3 domain of Bruton's tyrosine kinase (Btk) is preceded by the Tec homology (TH) region containing proline-rich sequences. We have studied a protein fragment containing both the Btk SH3 domain and the proline-rich sequences of the TH region (PRR-SH3). Intermolecular NMR cross-relaxation measurements, gel permeation chromatography profiles, titrations with proline-rich peptides, and (15)N NMR relaxation measurements are all consistent with a monomer-dimer equilibrium with a dissociation constant on the order of 60 microM. The intermolecular interactions do, at least in part, involve proline-rich sequences in the TH region. This behavior of Btk PRR-SH3 may have implications for the functional action of Btk.

In vivo nuclear delivery of oligonucleotides via hybridizing bifunctional peptides

Linking proteins directly to nucleic acids has been a complex task. By hybridizing a bifunctional peptide nucleic acid (PNA) consisting of a nucleic acid binding moiety and a nuclear localization signal (NLS) we have previously demonstrated that it is possible to link protein functions directly to nucleic acids containing a PNA target site. By hybridizing fluorescently labeled oligonucleotides to PNA-NLS molecules and subsequently transfecting different organs in vivo we demonstrate an active nuclear translocation of the PNA-NLS/oligonucleotide complex in different mouse organs.

Nucleocytoplasmic shuttling of Bruton's tyrosine kinase

Bruton's tyrosine kinase (Btk), a nonreceptor cytoplasmic tyrosine kinase belonging to the Tec family of kinases, has been shown to be critical for B cell proliferation, differentiation, and signaling. Loss-of-function mutations in the Btk gene lead to X-linked agammaglobulinemia (XLA), a primary immunodeficiency in humans, and the less severe condition xid in mice. Although Btk is mainly localized in the cytoplasm under steady state conditions, it translocates to the plasma membrane upon growth factor stimulation and cross-linking of the B cell receptor. Nevertheless, in ectopically as well as endogenously Btk-expressing cells, it can also translocate to the nucleus. Deletion of the pleckstrin homology (PH) domain (DeltaPH1) leads, however, to an even redistribution of Btk within the nucleus and cytoplasm in the majority of transfected cells. In contrast, an SH3-deleted (DeltaSH3) mutant of Btk has been found to be predominantly nuclear. We also demonstrate that the nuclear accumulation of DeltaPH1 is dependent on Src expression. This nucleocytoplasmic shuttling is sensitive to the exportin 1/CRM1-inactivating drug, leptomycin B, indicating that Btk utilizes functional nuclear export signals. In addition, while the DeltaPH1 mutant of Btk was found to be active and tyrosine-phosphorylated in vivo, DeltaSH3 displayed decreased autokinase activity and was not phosphorylated. Our findings indicate that the nucleocytoplasmic shuttling of Btk has implications regarding potential targets inside the nucleus, which may be critical in gene regulation during B cell development and differentiation.

The cellular phenotype conditions Btk for cell survival or apoptosis signaling

Bruton's agammaglobulinemia tyrosine kinase (Btk) was identified as the gene mutated in X-linked agammaglobulinemia. Btk is involved in B-cell receptor signaling and B cell ontogeny as the disease is characterized by a block in B-cell development. Cell proliferation and apoptosis are integral parts of B-cell development. We have demonstrated that overexpression of Btk in HeLa cells led to apoptosis, whereas in B cells, en dogenous levels of Btk protected the cells from apoptosis. We suggest a dual role for Btk in apoptosis and cell survival. We further propose that the phenotype of the cell may direct Btk for either cell survival or apoptosis. Our model is in line with the general feature of mammalian cells that have the inherent property to die unless the survival signals are triggered. The interplay between survival and apoptotic signaling is regulated by cell sur face receptors, cytoplasmic and nuclear regulatory molecules. These regu latory molecules may be simple unidirectional regulators or bidirectional regulators such as Btk. The different molecules involved in these pathways bring about an orchestrated signal depending on the phenotype of the cell, and a cell type-specific biological response is achieved that decides the fate of the cell.