Perspective: Reproductive isolation caused by natural selection against immigrants from divergent habitats

The classification of reproductive isolating barriers laid out by Dobzhansky and Mayr has motivated and structured decades of research on speciation. We argue, however, that this classification is incomplete and that the unique contributions of a major source of reproductive isolation have often been overlooked. Here, we describe reproductive barriers that derive from the reduced survival of immigrants upon reaching foreign habitats that are ecologically divergent from their native habitat. This selection against immigrants reduces encounters and thus mating opportunities between individuals from divergently adapted populations. It also reduces the likelihood that successfully mated immigrant females will survive long enough to produce their hybrid offspring. Thus, natural selection against immigrants results in distinctive elements of premating and postmating reproductive isolation that we hereby dub "immigrant inviability." We quantify the contributions of immigrant inviability to total reproductive isolation by examining study systems where multiple components of reproductive isolation have been measured and demonstrate that these contributions are frequently greater than those of traditionally recognized reproductive barriers. The relevance of immigrant inviability is further illustrated by a consideration of population-genetic theory, a review of selection against immigrant alleles in hybrid zone studies, and an examination of its participation in feedback loops that influence the evolution of additional reproductive barriers. Because some degree of immigrant inviability will commonly exist between populations that exhibit adaptive ecological divergence, we emphasize that these barriers play critical roles in ecological modes of speciation. We hope that the formal recognition of immigrant inviability and our demonstration of its evolutionary importance will stimulate more explicit empirical studies of its contributions to speciation.     

The 5-lipoxygenase pathway in arterial wall biology and atherosclerosis

Leukotrienes (LTs) are powerful inflammatory lipid mediators derived from the 5-lipoxygenase (5-LO) cascade of arachidonic acid. Recent clinical, population genetic, cell biological, and mouse studies indicate participation of the 5-LO pathway in atherogenesis and arterial wall remodeling. 5-LO is expressed by leukocytes including blood monocytes, tissue macrophages, dendritic cells, neutrophils, and mast cells. LTB4 and the cysteinyl LTs LTC4, LTD4, and LTE4, act through two BLT and two cysLT receptors that are differentially expressed on hematopoietic and arterial wall cells. The precise roles of LTs or the LT receptors in cardiovascular physiology remain largely to be explored. In this review, we will discuss what is currently known about the 5-LO atherosclerosis connection. We will attempt to propose strategies to further explore potential links between the 5-LO pathway and blood vessel physiology and disease progression.     

Macrophage-specific expression of group IIA sPLA2 results in accelerated atherogenesis by increasing oxidative stress

Group IIA secretory phospholipase A2 (sPLA2) is an acute-phase protein mediating decreased plasma HDL cholesterol and increased atherosclerosis. This study investigated the impact of macrophage-specific sPLA2 overexpression on lipoprotein metabolism and atherogenesis. Macrophages from sPLA2 transgenic mice have 2.5 times increased rates of LDL oxidation (thiobarbituric acid-reactive substances formation) in vitro (59 +/- 5 vs. 24 +/- 4 nmol malondialdehyde/mg protein; P < 0.001) dependent on functional 12/15-lipoxygenase (12/15-LO). Low density lipoprotein receptor-deficient (LDLR-/-) mice were transplanted with bone marrow from either sPLA2 transgenic mice (sPLA2--> LDLR-/-; n = 19) or wild-type C57BL/6 littermates (C57 BL/6-->LDLR-/-; n = 19) and maintained for 8 weeks on chow and then for 9 weeks on a Western-type diet. Plasma sPLA2 activity and plasma lipoprotein profiles were not significantly different between sPLA2-->LDLR-/- and C57BL/6-->LDLR-/- mice. Aortic root atherosclerosis was increased by 57% in sPLA2-->LDLR-/- mice compared with C57BL/6-->LDLR-/- controls (P < 0.05). Foam cell formation in vitro and in vivo was increased significantly. Urinary, plasma, and aortic levels of the isoprostane 8,12-iso-iPF2alpha-VI and aortic levels of 12/15-LO reaction products were each significantly higher (P < 0.001) in sPLA2-->LDLR-/- compared with C57BL/6-->LDLR-/- mice, indicating significantly increased in vivo oxidative stress in sPLA2--> LDLR-/-. These data demonstrate that macrophage-specific overexpression of human sPLA2 increases atherogenesis by directly modulating foam cell formation and in vivo oxidative stress without any effect on systemic sPLA2 activity and lipoprotein metabolism.     

Resting state conformation of the MsbA homodimer as studied by site-directed spin labeling

MsbA is the ABC transporter for lipid A and is found in the inner membranes of Gram-negative bacteria such as Escherichia coli. Without MsbA present, bacterial cells accumulate a toxic amount of lipid A within their inner membranes. A crystal structure of MsbA was recently obtained that provides an excellent starting point for functional dynamics studies in membranes [Chang, and Roth (2001) Science 293, 1793-1800]. Although a structure of MsbA is now available, many questions remain concerning its mechanism of transport. Site-directed spin labeling (SDSL) electron paramagnetic resonance (EPR) spectroscopy is a powerful approach for characterizing local areas within a large protein structure in addition to detecting and following changes in local structure due to dynamic interactions within a protein. The quaternary structure of the resting state of the MsbA homodimer reconstituted into lipid membranes has been evaluated by SDSL EPR spectroscopy and chemical cross-linking techniques. SDSL and cross-linking results are consistent with the controversial resting state conformation of the MsbA homodimer found in the crystal structure, with the tips of the transmembrane helices forming a dimer interface. The position of MsbA in the membrane bilayer along with the relative orientation of the transmembrane helical bundles with respect to one another has been determined. Characterization of the resting state of the MsbA homodimer is essential for future studies on the functional dynamics of this membrane transporter.     

Mutation Exposed: A Neutral Explanation for Extreme Base Composition of an Endosymbiont Genome

The influence of neutral mutation pressure versus selection on base composition evolution is a subject of considerable controversy. Yet the present study represents the first explicit population genetic analysis of this issue in prokaryotes, the group in which base composition variation is most dramatic. Here, we explore the impact of mutation and selection on the dynamics of synonymous changes in Buchnera aphidicola, the AT-rich bacterial endosymbiont of aphids. Specifically, we evaluated three forms of evidence. (i) We compared the frequencies of directional base changes (ATGC vs. GCAT) at synonymous sites within and between Buchnera species, to test for selective preference versus effective neutrality of these mutational categories. Reconstructed mutational changes across a robust intraspecific phylogeny showed a nearly 1:1 ATGC:GCAT ratio. Likewise, stationarity of base composition among Buchnera species indicated equal rates of ATGC and GCAT substitutions. The similarity of these patterns within and between species supported the neutral model. (ii) We observed an equivalence of relative per-site AT mutation rate and current AT content at synonymous sites, indicating that base composition is at mutational equilibrium. (iii) We demonstrated statistically greater equality in the frequency of mutational categories in Buchnera than in parallel mammalian studies that documented selection on synonymous sites. Our results indicate that effectively neutral mutational pressure, rather than selection, represents the major force driving base composition evolution in Buchnera. Thus they further corroborate recent evidence for the critical role of reduced N_e in the molecular evolution of bacterial endosymbionts.Reviewing Editor: Dr. J. William Ballard     

RNA Structures Required for Production of Subgenomic Flavivirus RNA

Flaviviruses are a group of single-stranded, positive-sense RNA viruses causing ∼100 million infections per year. We have recently shown that flaviviruses produce a unique, small, noncoding RNA (∼0.5 kb) derived from the 3′ untranslated region (UTR) of the genomic RNA (gRNA), which is required for flavivirus-induced cytopathicity and pathogenicity (G. P. Pijlman et al., Cell Host Microbe, 4: 579-591, 2008). This RNA (subgenomic flavivirus RNA [sfRNA]) is a product of incomplete degradation of gRNA presumably by the cellular 5′-3′ exoribonuclease XRN1, which stalls on the rigid secondary structure stem-loop II (SL-II) located at the beginning of the 3′ UTR. Mutations or deletions of various secondary structures in the 3′ UTR resulted in the loss of full-length sfRNA (sfRNA1) and production of smaller and less abundant sfRNAs (sfRNA2 and sfRNA3). Here, we investigated in detail the importance of West Nile virus Kunjin (WNV_KUN) 3′ UTR secondary structures as well as tertiary interactions for sfRNA formation. We show that secondary structures SL-IV and dumbbell 1 (DB1) downstream of SL-II are able to prevent further degradation of gRNA when the SL-II structure is deleted, leading to production of sfRNA2 and sfRNA3, respectively. We also show that a number of pseudoknot (PK) interactions, in particular PK1 stabilizing SL-II and PK3 stabilizing DB1, are required for protection of gRNA from nuclease degradation and production of sfRNA. Our results show that PK interactions play a vital role in the production of nuclease-resistant sfRNA, which is essential for viral cytopathicity in cells and pathogenicity in mice.Arthropod-borne flaviviruses such as West Nile virus (WNV), dengue virus (DENV), and Japanese encephalitis virus (JEV) cause major outbreaks of potentially fatal disease and affect over 50 million people every year. The highly pathogenic North American strain of WNV (WNV_NY99) has already claimed more than 1,000 lives with over 27,000 cases reported since its emergence in New York in 1999 and has raised global public health concerns (9). In contrast, the closely related Australian strain of WNV, WNV_KUN, is highly attenuated and does not cause overt disease in humans and animals (11). WNV_KUN has been used extensively as a model virus to study flavivirus replication and flavivirus-host interactions (13, 14, 16-19, 26, 38, 39).The ∼11-kb positive-stranded, capped WNV genomic RNA (gRNA) lacks a poly(A) tail and consists of 5′ and 3′ untranslated regions (UTRs) flanking one open reading frame, which encodes the viral proteins required for the viral life cycle (6, 15, 38, 39). Flavivirus UTRs are involved in translation and initiation of RNA replication and likely determine genome packaging (13, 14, 16, 21, 30, 39-41). Both the 5′ UTR (∼100 nucleotides [nt] in size) and the 3′ UTR (from ∼400 to 700 nucleotides) can form secondary and tertiary structures which are highly conserved among mosquito-borne flaviviruses (1, 8, 10, 14, 29, 32, 34). More specifically, the WNV_KUN 3′ UTR consists of several conserved regions and secondary structures (Fig. ​(Fig.1A)1A) which were previously predicted or shown to exist in various flaviviruses by computational and chemical analyses, respectively (4, 10, 25, 26, 29-32). The 5′ end of the 3′ UTR starts with an AU-rich region which can form stem-loop structure I (SL-I) followed by SL-II, which we previously showed to be vitally important for subgenomic flavivirus RNA (sfRNA) production (26; see also below). SL-II is followed by a short, repeated conserved hairpin (RCS3) and SL-III (26). Further downstream of SL-III are the SL-IV and CS3 structures, which are remarkably similar to the preceding SL-II-RCS3 structure (26, 29). Further downstream of the SL-IV-CS3 structure are dumbbells 1 and 2 (DB1 and DB2, respectively) followed by a short SL and the 3′ SL (25, 26).Open in a separate windowFIG. 1.(A) Model of the WNV_KUN 3′ UTR RNA structure. Highlighted in bold are the secondary structures investigated here. Dashed lines indicate putative PKs. The two sites of the putative PK interactions are shown in open boxes. sfRNA1, -2, -3, and -4 start sites are indicated by arrows. (R)CS, (repeated) conserved sequence; DB, dumbbell structure; PK, pseudoknot; SL, stem-loop. (B) Structural model of PK1 in SL-II with disruptive mutations. Nucleotide numbering is from the end of the 3′ UTR. The sfRNA1 start is indicated by an arrow. Nucleotides forming PK1 are on a gray background, and mutated nucleotides are white on a black background. (C) Sequences mutated in the different constructs. Nucleotides in the wt PK sequences used for mutations are bold and underlined. Introduced mutations are shown under the corresponding nucleotides in the wt sequence.The described structures have been investigated in some detail for their requirement in RNA replication and translation. Generally, a progressive negative effect on viral growth was shown with progressive deletions into the 3′-proximal region of the JEV 3′ UTR (41). However, only a relatively short region of the JEV 3′ UTR, consisting of the 3′-terminal 193 nt, was shown to be absolutely essential for gRNA replication (41). The minimal region for DENV replication was reported to be even shorter (23). Extensive analysis has shown that the most 3′-terminal, essential regions of the 3′ UTR include the cyclization sequence and 3′ SL, which are required for efficient RNA replication (2, 14, 16, 23, 35). As we showed, deletion of SL-II or SL-I did not overtly affect WNV_KUN replication (26). However, deletion of CS2, RCS2, CS3, or RCS3 in WNV replicon RNA significantly reduced RNA replication but not translation (20), indicating that these elements facilitate but are not essential for RNA replication. In addition, it was shown that deletion of DB1 or DB2 resulted in a viable mutant virus that was reduced in growth efficiency, while deletion of both DB structures resulted in a nonviable mutant (23).In addition to the above-mentioned secondary stem-loop structures, computational and chemical analysis of the flavivirus 3′ UTR suggested the presence of 5 pseudoknot (PK) interactions (Fig. ​(Fig.1A)1A) (25, 26, 32). A PK is a structure formed upon base pairing of a single-stranded region of RNA in the loop of a hairpin to a stretch of complementary nucleotides elsewhere in the RNA chain (Fig. ​(Fig.1B).1B). These structures are referred to as hairpin type (H-type) PKs (3), and they usually stabilize secondary RNA structures. Typically, the final tertiary structure does not significantly alter the preformed secondary structure (5). In general, PK interactions have been shown to be important in biological processes such as initiation and/or elongation of translation, initiation of gRNA replication, and ribosomal frameshifting for a number of different viruses, including flaviviruses (reviewed in references 3 and 22). The first PK in the WNV 3′ UTR was predicted to form in SL-II, followed by a similar PK in SL-IV (26) (PK1 and PK2 in Fig. ​Fig.1A).1A). For the DENV, yellow fever virus (YFV), and JEV subgroup of flaviviruses, two PKs further downstream were predicted to form between DB1 and DB2 and corresponding single-stranded RNA regions located further downstream (25) (PK3 and PK4 in Fig. ​Fig.1A).1A). The formation of these structures is supported by covariations in the WNV RNAs. In addition, a PK was proposed to form between a short SL and the 3′ SL at the 3′ terminus of the viral genome (32) (PK5 in Fig. ​Fig.1A1A).Importantly, in addition to its role in viral replication and translation, we have shown that the WNV_KUN 3′ UTR is important for the production of a small noncoding RNA fragment designated sfRNA (26). This short RNA fragment of ∼0.5 kb is derived from the 3′ UTR of the gRNA and exclusively produced by the members of the Flavivirus genus of the Flaviviridae family, where it is required for efficient viral replication, cytopathicity, and pathogenicity (26). Our studies suggested that sfRNA is a product of incomplete degradation of the gRNA presumably by the cellular 5′-3′ exoribonuclease XRN1, resulting from XRN1 stalling on the rigid secondary/tertiary structures located at the beginning of the 3′ UTR (26). XRN1 is an exoribonuclease which usually degrades mRNA from the 5′ to the 3′ end as part of cellular mRNA decay and turnover (33), and it was shown previously that XRN1 can be stalled by SL structures (28). Mutations or deletions of WNV 3′ UTR secondary structures resulted in the loss of full-length sfRNA (sfRNA1) and production of smaller and less abundant sfRNAs (sfRNA2 and sfRNA3) (26). In particular, SL-II (Fig. ​(Fig.1A)1A) was shown to be important for sfRNA1 production; deletion of this structure either alone or in conjunction with other structures located downstream of SL-II abolished sfRNA1 production, leading to the production of the smaller RNA fragments sfRNA2 and sfRNA3.Here, we extended our investigation and studied the importance of several predicted 3′ UTR secondary structures and PK interactions for the production of sfRNA. To further understand the generation mechanism of sfRNA and its requirements, we deleted or mutated a number of RNA structures in the WNV_KUN 3′ UTR and investigated the size and amount of sfRNA generated from these mutant RNAs. The results show that not only SLs but also PK interactions play a vital role in stabilizing the 3′ UTR RNA and preventing complete degradation of viral gRNA to produce nuclease-resistant sfRNA, which is required for efficient virus replication and cytopathicity in cells and virulence in mice.     

Genetic model of selective COX2 inhibition reveals novel heterodimer signaling

Selective inhibitors of cyclooxygenase-2 (COX2) have attracted widespread media attention because of evidence of an elevated risk of cardiovascular complications in placebo-controlled trials, resulting in the market withdrawal of some members of this class. These drugs block the cyclooxygenase activity of prostaglandin H synthase-2 (PGHS2), but do not affect the associated peroxidase function. They were developed with the rationale of conserving the anti-inflammatory and analgesic actions of traditional nonsteroidal anti-inflammatory drugs (tNSAIDs) while sparing the ability of PGHS1-derived prostaglandins to afford gastric cytoprotection. PGHS1 and PGHS2 coexist in the vasculature and in macrophages, and are upregulated together in inflammatory tissues such as rheumatoid synovia and atherosclerotic plaque. They are each believed to function as homodimers. Here, we developed a new genetic mouse model of selective COX2 inhibition using a gene-targeted point mutation, resulting in a Y385F substitution. Structural modeling and biochemical assays showed the ability of PGHS1 and PGHS2 to heterodimerize and form prostaglandins. The heterodimerization of PGHS1-PGHS2 may explain how the ductus arteriosus closes normally at birth in mice expressing PGHS2 Y385F, but not in PGHS2-null mice.     

High levels of cryptic species diversity uncovered in Amazonian frogs

One of the greatest challenges for biodiversity conservation is the poor understanding of species diversity. Molecular methods have dramatically improved our ability to uncover cryptic species, but the magnitude of cryptic diversity remains unknown, particularly in diverse tropical regions such as the Amazon Basin. Uncovering cryptic diversity in amphibians is particularly pressing because amphibians are going extinct globally at an alarming rate. Here, we use an integrative analysis of two independent Amazonian frog clades, Engystomops toadlets and Hypsiboas treefrogs, to test whether species richness is underestimated and, if so, by how much. We sampled intensively in six countries with a focus in Ecuador (Engystomops: 252 individuals from 36 localities; Hypsiboas: 208 individuals from 65 localities) and combined mitochondrial DNA, nuclear DNA, morphological, and bioacoustic data to detect cryptic species. We found that in both clades, species richness was severely underestimated, with more undescribed species than described species. In Engystomops, the two currently recognized species are actually five to seven species (a 150-250% increase in species richness); in Hypsiboas, two recognized species represent six to nine species (a 200-350% increase). Our results suggest that Amazonian frog biodiversity is much more severely underestimated than previously thought.     

Modeling perimenopause in Sprague-Dawley rats by chemical manipulation of the transition to ovarian failure

Various age-related diseases increase in incidence during perimenopause. However, our understanding of the effects of aging compared with hormonal changes of perimenopause in mediating these disease risks is incomplete, in part due to the lack of an experimental perimenopause model. We therefore aimed to determine whether manipulation of the transition to ovarian failure in rats via the use of 4-vinylcyclohexene diepoxide (VCD) could be used to model and accelerate hormonal changes characteristic of perimenopause. We examined long-term (11 to 20 mo), dose-dependent effects of VCD on reproductive function in 1- and 3-mo-old female Sprague-Dawley rats. Twenty-five daily doses of VCD (80 or 160 mg/kg daily compared with vehicle alone) depleted ovarian follicles in a dose-dependent fashion in rats of both ages, accelerated the onset of acyclicity, and caused dose-dependent increases in follicle-stimulating hormone that exceeded those naturally occurring with age in control rats but left serum levels of 17β-estradiol unchanged, with continued ovarian production of androstenedione. High-dose VCD caused considerable nonovarian toxicities in 3-mo-old Sprague-Dawley rats, making this an unsuitable model. In contrast, 1-mo-old rats had more robust dose-dependent increases in follicle-stimulating hormone without evidence of systemic toxicity in response to either VCD dose. Because perimenopause is characterized by an increase in follicle-stimulating hormone with continued secretion of ovarian steroids, VCD acceleration of an analogous hormonal milieu in 1-mo-old Sprague-Dawley rats may be useful for probing the hormonal effects of perimenopause on age-related disease risk.