Systems biology meets stress ecology: linking molecular and organismal stress responses in Daphnia magna

Background: Ibuprofen and other nonsteroidal anti-inflammatory drugs have been designed to interrupt eicosanoid metabolism in mammals, but little is known of how they affect nontarget organisms. Here we report a systems biology study that simultaneously describes the transcriptomic and phenotypic stress responses of the model crustacean Daphnia magna after exposure to ibuprofen. Results: Our findings reveal intriguing similarities in the mode of action of ibuprofen between vertebrates and invertebrates, and they suggest that ibuprofen has a targeted impact on reproduction at the molecular, organismal, and population level in daphnids. Microarray expression and temporal real-time quantitative PCR profiles of key genes suggest early ibuprofen interruption of crustacean eicosanoid metabolism, which appears to disrupt signal transduction affecting juvenile hormone metabolism and oogenesis. Conclusion: Combining molecular and organismal stress responses provides a guide to possible chronic consequences of environmental stress for population health. This could improve current environmental risk assessment by providing an early indication of the need for higher tier testing. Our study demonstrates the advantages of a systems approach to stress ecology, in which Daphnia will probably play a major role.

Amphioxus molecular biology: insights into vertebrate evolution and developmental mechanisms

The cephalochordate amphioxus is the best available proxy for the last common invertebrate ancestor of the vertebrates. During the last decade, the developmental genetics of amphioxus have been extensively examined for insights into the evolutionary origin and early evolution of the vertebrates. Comparisons between expression domains of homologous genes in amphioxus and vertebrates have strengthened proposed homologies between specific body parts. Molecular genetic studies have also highlighted parallels in the developmental mechanisms of amphioxus and vertebrates. In both groups, a similar nested pattern of Hox gene expression is involved in rostrocaudal patterning of the neural tube, and homologous genes also appear to be involved in dorsoventral neural patterning. Studies of amphioxus molecular biology have also hinted that the protochordate ancestor of the vertebrates included cell populations that modified their developmental genetic pathways during early vertebrate evolution to yield definitive neural crest and neurogenic placodes. We also discuss how the application of expressed sequence tag and gene-mapping approaches to amphioxus have combined with developmental studies to advance our understanding of chordate genome evolution. We conclude by considering the potential offered by the sequencing of the amphioxus genome, which was completed in late 2004.

Comparative genome analysis provides insights into the evolution and adaptation of Pseudomonas syringae pv. aesculi on Aesculus hippocastanum.

A recently emerging bleeding canker disease, caused by Pseudomonas syringae pathovar aesculi (Pae), is threatening European horse chestnut in northwest Europe. Very little is known about the origin and biology of this new disease. We used the nucleotide sequences of seven commonly used marker genes to investigate the phylogeny of three strains isolated recently from bleeding stem cankers on European horse chestnut in Britain (E-Pae). On the basis of these sequences alone, the E-Pae strains were identical to the Pae type-strain (I-Pae), isolated from leaf spots on Indian horse chestnut in India in 1969. The phylogenetic analyses also showed that Pae belongs to a distinct clade of P. syringae pathovars adapted to woody hosts. We generated genome-wide Illumina sequence data from the three E-Pae strains and one strain of I-Pae. Comparative genomic analyses revealed pathovar-specific genomic regions in Pae potentially implicated in virulence on a tree host, including genes for the catabolism of plant-derived aromatic compounds and enterobactin synthesis. Several gene clusters displayed intra-pathovar variation, including those encoding type IV secretion, a novel fatty acid biosynthesis pathway and a sucrose uptake pathway. Rates of single nucleotide polymorphisms in the four Pae genomes indicate that the three E-Pae strains diverged from each other much more recently than they diverged from I-Pae. The very low genetic diversity among the three geographically distinct E-Pae strains suggests that they originate from a single, recent introduction into Britain, thus highlighting the serious environmental risks posed by the spread of an exotic plant pathogenic bacterium to a new geographic location. The genomic regions in Pae that are absent from other P. syringae pathovars that infect herbaceous hosts may represent candidate genetic adaptations to infection of the woody parts of the tree.

Comparisons of host mitochondrial, nuclear and endosymbiont bacterial genes reveal cryptic fig wasp species and the effects of Wolbachia on host mtDNA evolution and diversity

Background Figs and fig-pollinating wasp species usually display a highly specific one-to-one association. However, more and more studies have revealed that the "one-to-one" rule has been broken. Co-pollinators have been reported, but we do not yet know how they evolve. They may evolve from insect speciation induced or facilitated by Wolbachia which can manipulate host reproduction and induce reproductive isolation. In addition, Wolbachia can affect host mitochondrial DNA evolution, because of the linkage between Wolbachia and associated mitochondrial haplotypes, and thus confound host phylogeny based on mtDNA. Previous research has shown that fig wasps have the highest incidence of Wolbachia infection in all insect taxa, and Wolbachia may have great influence on fig wasp biology. Therefore, we look forward to understanding the influence of Wolbachia on mitochondrial DNA evolution and speciation in fig wasps. Results We surveyed 76 pollinator wasp specimens from nine Ficus microcarpa trees each growing at a different location in Hainan and Fujian Provinces, China. We found that all wasps were morphologically identified as Eupristina verticillata, but diverged into three clades with 4.22-5.28% mtDNA divergence and 2.29-20.72% nuclear gene divergence. We also found very strong concordance between E. verticillata clades and Wolbachia infection status, and the predicted effects of Wolbachia on both mtDNA diversity and evolution by decreasing mitochondrial haplotypes. Conclusions Our study reveals that the pollinating wasp E. verticillata on F. microcarpa has diverged into three cryptic species, and Wolbachia may have a role in this divergence. The results also indicate that Wolbachia strains infecting E. verticillata have likely resulted in selective sweeps on host mitochondrial DNA.

Intra-tumor heterogeneity: lessons from microbial evolution and clinical implications

Multiple subclonal populations of tumor cells can coexist within the same tumor. This intra-tumor heterogeneity will have clinical implications and it is therefore important to identify factors that drive or suppress such heterogeneous tumor progression. Evolutionary biology can provide important insights into this process. In particular, experimental evolution studies of microbial populations, which exist as clonal populations that can diversify into multiple subclones, have revealed important evolutionary processes driving heterogeneity within a population. There are transferrable lessons that can be learnt from these studies that will help us to understand the process of intra-tumor heterogeneity in the clinical setting. In this review, we summarize drivers of microbial diversity that have been identified, such as mutation rate and environmental influences, and discuss how knowledge gained from microbial experimental evolution studies may guide us to identify and understand important selective factors that promote intra-tumor heterogeneity. Furthermore, we discuss how these factors could be used to direct and optimize research efforts to improve patient care, focusing on therapeutic resistance. Finally, we emphasize the need for longitudinal studies to address the impact of these potential tumor heterogeneity-promoting factors on drug resistance, metastatic potential and clinical outcome.

The mean evolution and variability of the Asian summer monsoon: Comparison of ECMWF and NCEP-NCAR reanalyses

The behavior of the Asian summer monsoon is documented and compared using the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis (ERA) and the National Centers for Environmental Prediction-National Center for Atmospheric Research (NCEP-NCAR) Reanalysis. In terms of seasonal mean climatologies the results suggest that, in several respects, the ERA is superior to the NCEP-NCAR Reanalysis. The overall better simulation of the precipitation and hence the diabatic heating field over the monsoon domain in ERA means that the analyzed circulation is probably nearer reality. In terms of interannual variability, inconsistencies in the definition of weak and strong monsoon years based on typical monsoon indices such as All-India Rainfall (AIR) anomalies and the large-scale wind shear based dynamical monsoon index (DMI) still exist. Two dominant modes of interannual variability have been identified that together explain nearly 50% of the variance. Individually, they have many features in common with the composite flow patterns associated with weak and strong monsoons, when defined in terms of regional AIR anomalies and the large-scale DMI. The reanalyses also show a common dominant mode of intraseasonal variability that describes the latitudinal displacement of the tropical convergence zone from its oceanic-to-continental regime and essentially captures the low-frequency active/break cycles of the monsoon. The relationship between interannual and intraseasonal variability has been investigated by considering the probability density function (PDF) of the principal component of the dominant intraseasonal mode. Based on the DMI, there is an indication that in years with a weaker monsoon circulation, the PDF is skewed toward negative values (i,e., break conditions). Similarly, the PDFs for El Nino and La Nina years suggest that El Nino predisposes the system to more break spells, although the sample size may limit the statistical significance of the results.

A new look at stratospheric sudden warmings. Part III. Polar vortex evolution and vertical structure

The evolution of the Arctic polar vortex during observed major mid-winter stratospheric sudden warmings (SSWs) is investigated for the period 1957-2002, using European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-40 Ertel’s potential vorticity (PV) and temperature fields. Time-lag composites of vertically weighted PV, calculated relative to the SSW onset time, are derived for both vortex displacement SSWs and vortex splitting SSWs, by averaging over the 15 recorded displacement and 13 splitting events. The evolving vertical structure of the polar vortex during a typical SSW of each type is clearly illustrated by plotting an isosurface of the composite PV field, and is shown to be very close to that observed during representative individual events. Results are verified by comparison with an elliptical diagnostic vortex moment technique. For both types of SSW, little variation is found between individual events in the orientation of the developing vortex relative to the underlying topography, i.e. the location of the vortex during SSWs of each type is largely fixed in relation to the Earth’s surface. During each type of SSW, the vortex is found to have a distinctive vertical structure. Vortex splitting events are typically barotropic, with the vortex split occurring near-simultaneously over a large altitude range (20-40 km). In the majority of cases, of the two daughter vortices formed, it is the ‘Siberian’ vortex that dominates over its ‘Canadian’ counterpart. In contrast, displacement events are characterized by a very clear baroclinic structure; the vortex tilts significantly westward with height, so that the top and bottom of the vortex are separated by nearly 180◦ longitude before the upper vortex is sheared away and destroyed.

Temporal evolution and electric potential structure of the auroral acceleration region from multispacecraft measurements

Bright aurorae can be excited by the acceleration of electrons into the atmosphere in violation of ideal magnetohydrodynamics. Modelling studies predict that the accelerating electric potential consists of electric double layers at the boundaries of an acceleration region but observations suggest that particle acceleration occurs throughout this region. Using multi-spacecraft observations from Cluster we have examined two upward current regions on 14 December 2009. Our observations show that the potential difference below C4 and C3 changed by up to 1.7 kV between their respective crossings, which were separated by 150 s. The field-aligned current density observed by C3 was also larger than that observed by C4. The potential drop above C3 and C4 was approximately the same in both crossings. Using a novel technique of quantitatively comparing the electron spectra measured by Cluster 1 and 3, which were separated in altitude, we determine when these spacecraft made effectively magnetically conjugate observations and use these conjugate observations to determine the instantaneous distribution of the potential drop in the AAR. Our observations show that an average of 15% of the potential drop in the AAR was located between C1 at 6235 km and C3 at 4685 km altitude, with a maximum potential drop between the spacecraft of 500~V and that the majority of the potential drop was below C3. By assuming a spatial invariance along the length of the upward current region, we discuss these observations in terms of temporal changes and the vertical structure of the electrostatic potential drop and in the context of existing models and previous observations single- and multi-spacecraft observations.