Integration of biological networks and gene expression data using Cytoscape

Cytoscape is a free software package for visualizing, modeling and analyzing molecular and genetic interaction networks. This protocol explains how to use Cytoscape to analyze the results of mRNA expression profiling, and other functional genomics and proteomics experiments, in the context of an interaction network obtained for genes of interest. Five major steps are described: (i) obtaining a gene or protein network, (ii) displaying the network using layout algorithms, (iii) integrating with gene expression and other functional attributes, (iv) identifying putative complexes and functional modules and (v) identifying enriched Gene Ontology annotations in the network. These steps provide a broad sample of the types of analyses performed by Cytoscape.     

Population persistence under high mutation rate: From evolutionary rescue to lethal mutagenesis

Populations may genetically adapt to severe stress that would otherwise cause their extirpation. Recent theoretical work, combining stochastic demography with Fisher's geometric model of adaptation, has shown how evolutionary rescue becomes unlikely beyond some critical intensity of stress. Increasing mutation rates may however allow adaptation to more intense stress, raising concerns about the effectiveness of treatments against pathogens. This previous work assumes that populations are rescued by the rise of a single resistance mutation. However, even in asexual organisms, rescue can also stem from the accumulation of multiple mutations in a single genome. Here, we extend previous work to study the rescue process in an asexual population where the mutation rate is sufficiently high so that such events may be common. We predict both the ultimate extinction probability of the population and the distribution of extinction times. We compare the accuracy of different approximations covering a large range of mutation rates. Moderate increase in mutation rates favors evolutionary rescue. However, larger increase leads to extinction by the accumulation of a large mutation load, a process called lethal mutagenesis. We discuss how these results could help design “evolution‐proof” antipathogen treatments that even highly mutable strains could not overcome.     

The locus encoding an oligomorphic family of <Emphasis Type="Italic">MHC-A</Emphasis> alleles (<Emphasis Type="Italic">Mane-A</Emphasis>*06/<Emphasis Type="Italic">Mamu-A</Emphasis>*05) is present at high frequency in several macaque species

Several macaques species are used for HIV pathogenesis and vaccine studies, and the characterization of their major histocompatibility complex (MHC) class I genes is required to rigorously evaluate the cellular immune responses induced after immunization and/or infection. In this study, we demonstrate that the gene expressing the Mane-A*06 allele of pig-tailed macaques is an orthologue of the locus encoding the Mamu-A*05 allele family in rhesus macaques. Analysis of the distribution of this locus in a cohort of 63 pig-tailed macaques revealed that it encodes an oligomorphic family of alleles, highly prevalent (90%) in the pig-tailed macaque population. Similarly, this locus was very frequently found (62%) in a cohort of 80 Indian rhesus macaques. An orthologous gene was also detected in cynomolgus monkeys originating from four different geographical locations, but was absent in two African monkey species. Expression analysis in pig-tailed macaques revealed that the Mane-A*06 alleles encoded by this locus are transcribed at 10- to 20-fold lower levels than other MHC-A alleles (Mane-A*03 or Mane-A*10). Despite their conservation and high prevalence among Asian macaque species, the alleles of the Mane-A*06 family and, by extension their orthologues in rhesus and cynomolgus monkeys, may only modestly contribute to cellular immune responses in macaques because of their low level of expression.     

Synteny on mouse chromosome 5 of homologs for human DNA loci linked to the Huntington disease gene

Comparative mapping in man and mouse has revealed frequent conservation of chromosomal segments, offering a potential approach to human disease genes via their murine homologs. Using DNA markers near the Huntington disease gene on the short arm of chromosome 4, we defined a conserved linkage group on mouse chromosome 5. Linkage analyses using recombinant inbred strains, a standard outcross, and an interspecific backcross were used to assign homologs for five human loci, D4S43, D4S62, QDPR, D4S76, and D4S80, to chromosome 5 and to determine their relationships with previously mapped markers for this autosome. The relative order of the conserved loci was preserved in a linkage group that spanned 13% recombination in the interspecific backcross analysis. The most proximal of the conserved markers on the mouse map, D4S43h, showed no recombination with Emv-1, an endogenous ecotropic virus, in 84 outcross progeny and 19 recombinant inbred strains. Hx, a dominant mutation that causes deformities in limb development, maps approximately 2 cM proximal to Emv-1. Since the human D4S43 locus is less than 1 cM proximal to HD near the telomere of chromosome 4, the murine counterpart of the HD gene might lie between Hx and Emv-1 or D4S43h. Cloning of the region between these markers could generate new probes for conserved human sequences in the vicinity of the HD gene or possibly candidates for the murine counterpart of this human disease locus.     

Frequency analysis of atrial fibrillation from the surface electrocardiogram

Atrial fibrillation (AF) is the most common arrhythmia encountered in clinical practice. Neither the natural history of AF nor its response to therapy are sufficiently predictable by clinical and echocardiographic parameters. Atrial fibrillatory frequency (or rate) can reliably be assessed from the surface electrocardiogram (ECG) using digital signal processing (filtering, subtraction of averaged QRST complexes, and power spectral analysis) and shows large inter-individual variability. This measurement correlates well with intraatrial cycle length, a parameter which appears to have primary importance in AF domestication and response to therapy. AF with a low fibrillatory rate is more likely to terminate spontaneously, and responds better to antiarrhythmic drugs or cardioversion while high rate AF is more often persistent and refractory to therapy. In conclusion, frequency analysis of AF seems to be useful for non-invasive assessment of electrical remodeling in AF and may subsequently be helpful for guiding AF therapy.     

Soluble Fas ligand induces epithelial cell apoptosis in humans with acute lung injury (ARDS).

The goals of this study were to determine whether the Fas-dependent apoptosis pathway is active in the lungs of patients with the acute respiratory distress syndrome (ARDS), and whether this pathway can contribute to lung epithelial injury. We found that soluble Fas ligand (sFasL) is present in bronchoalveolar lavage (BAL) fluid of patients before and after the onset of ARDS. The BAL concentration of sFasL at the onset of ARDS was significantly higher in patients who died. BAL from patients with ARDS induced apoptosis of distal lung epithelial cells, which express Fas, and this effect was inhibited by blocking the Fas/FasL system using three different strategies: anti-FasL mAb, anti-Fas mAb, and a Fas-Ig fusion protein. In contrast, BAL from patients at risk for ARDS had no effect on distal lung epithelial cell apoptosis. These data indicate that sFasL is released in the airspaces of patients with acute lung injury and suggest that activation of the Fas/FasL system contributes to the severe epithelial damage that occurs in ARDS. These data provide the first evidence that FasL can be released as a biologically active, death-inducing mediator capable of inducing apoptosis of cells of the distal pulmonary epithelium during acute lung injury.     

Migration of herbaceous plant species across ancient–recent forest ecotones in central Belgium

1 We studied the migration of forest plant species using their percentage cover and frequency in 197 plots distributed over 26 transects across ecotones between ancient and recent deciduous forests in the Meerdaal forest complex in central Belgium. The recent forest stands varied in age between 36 and 132 years, and all occurred on silty, well-drained soils.
2 The total cover, number and diversity of field layer species did not differ significantly between ancient and recent forest stands.
3 The number and cover of the ancient forest plant species and of ant-dispersed species correlated positively with the age of the recent forest and negatively with both the duration of its former agricultural land use and the distance to the ancient forest. This implies a slow colonization of the recent forest stands by these species; all species were, however, able to migrate across the ecotones.
4 The cover of four species ( Anemone nemorosa , Lamium galeobdolon , Convallaria majalis and Polygonatum multiflorum ) declined along the transect, suggesting that they are limited by seed dispersal. Their colonization rates, calculated from the occurrence of the farthest individual, ranged from < 0.05 to 1.15 m year⁻¹ and for other measures from < 0.05 to 0.65 m year⁻¹. Anemone and Lamium appeared to colonize the recent forest by establishment of isolated individuals, while Polygonatum and Convallaria expanded populations from existing patches on the border between ancient and recent forest.
5 Several forest species were able to colonize the recent forest rapidly, where some of them even reached a higher abundance, due to the increased availability of colonization sites with a higher nutrient content and a thinner organic layer.     

Death-associated protein kinase phosphorylates mammalian ribosomal protein S6 and reduces protein synthesis

Death-associated protein kinase (DAPK) is a pro-apoptotic, calcium/calmodulin-regulated protein kinase that is a drug discovery target for neurodegenerative disorders. Despite the potential profound physiological role of DAPK in neuronal function and pathophysiology, the endogenous substrate(s) of this kinase and the mechanisms via which DAPK elicits its biological action remain largely unknown. We report here that the mammalian 40S ribosomal protein S6 is a DAPK substrate. Results from immunoprecipitation experiments are consistent with endogenous DAPK being associated with endogenous S6 in rat brain. When S6 is a component of the 40S ribosomal subunit complex, DAPK selectively phosphorylates it at serine 235, one of the five sites in S6 that are phosphorylated by the S6 kinase family of proteins. The amino acid sequence flanking serine 235 matches the established pattern for DAPK peptide and protein substrates. Kinetic analyses using purified 40S subunits revealed a K(m) value of 9 microM, consistent with S6 being a potential physiological substrate of DAPK. This enzyme-substrate relationship has functional significance. DAPK suppresses translation in rabbit reticulocyte lysate, and treatment of neuroblastoma cells with a stimulator of DAPK reduces protein synthesis. In both cases, suppression of translation correlates with increased phosphorylation of S6 at serine 235. These results demonstrate that DAPK is a S6 kinase and provide evidence for a novel role of DAPK in the regulation of translation.     

Spore and micro-particle capture on an immunosensor surface in an ultrasound standing wave system

The capture of Bacillus subtilis var. niger spores on an antibody-coated surface can be enhanced when that coated surface acts as an acoustic reflector in a quarter wavelength ultrasonic (3 MHz) standing wave resonator. Immunocapture in such a resonator has been characterised here for both spores and 1 microm diameter biotinylated fluorescent microparticles. A mean spatial acoustic pressure amplitude of 460 kPa and a frequency of 2.82 MHz gave high capture efficiencies. It was shown that capture was critically dependent on reflector thickness. The time dependence of particle deposition on a reflector in a batch system was broadly consistent with a calculated time of 35 s to bring 95% of particles to the coated surface. A suspension flow rate of 0.1 ml/min and a reflector thickness of 1.01 mm gave optimal capture in a 2 min assay. The enhancement of particle detection compared with the control (no ultrasound) situation was x 70. The system detects a total of five particles in 15 fields of view in a 2 min assay when the suspending phase concentration was 10(4) particles/ml. A general expression for the dependence of minimum concentration detectable on; number of fields examined, sample volume flowing through the chamber and assay time shows that, for a practical combination of these variables, the threshold detection concentration can be two orders of magnitude lower.     

Adaptive origins of primates revisited

Interpretation of the adaptive profile of ancestral primates is controversial and has been constrained for decades by general acceptance of the premise that the first primates were very small. Here we show that neither the fossil record nor modern species provide evidence that the last common ancestor of living primates was small. Instead, comparative weight distributions of arboreal mammals and a phylogenetic reconstruction of ancestral primate body mass indicate that the reduction of functional claws to nails -- a primate characteristic that had up until now eluded satisfactory explanation - resulted from an increase in body mass to around 1000 g or more in the primate stem lineage. The associated shift to a largely vegetarian diet coincided with increased angiosperm diversity and the evolution of larger fruit size during the Late Cretaceous.