Simultaneous Non-Negative Matrix Factorization for Multiple Large Scale Gene Expression Datasets in Toxicology

Non-negative matrix factorization is a useful tool for reducing the dimension of large datasets. This work considers simultaneous non-negative matrix factorization of multiple sources of data. In particular, we perform the first study that involves more than two datasets. We discuss the algorithmic issues required to convert the approach into a practical computational tool and apply the technique to new gene expression data quantifying the molecular changes in four tissue types due to different dosages of an experimental panPPAR agonist in mouse. This study is of interest in toxicology because, whilst PPARs form potential therapeutic targets for diabetes, it is known that they can induce serious side-effects. Our results show that the practical simultaneous non-negative matrix factorization developed here can add value to the data analysis. In particular, we find that factorizing the data as a single object allows us to distinguish between the four tissue types, but does not correctly reproduce the known dosage level groups. Applying our new approach, which treats the four tissue types as providing distinct, but related, datasets, we find that the dosage level groups are respected. The new algorithm then provides separate gene list orderings that can be studied for each tissue type, and compared with the ordering arising from the single factorization. We find that many of our conclusions can be corroborated with known biological behaviour, and others offer new insights into the toxicological effects. Overall, the algorithm shows promise for early detection of toxicity in the drug discovery process.     

Sexual signalling in female crested macaques and the evolution of primate fertility signals

Background

Host-parasite coevolution can lead to local adaptation of either parasite or host if there is specificity (GxG interactions) and asymmetric evolutionary potential between host and parasite. This has been demonstrated both experimentally and in field studies, but a substantial proportion of studies fail to detect such clear-cut patterns. One explanation for this is that adaptation can be masked by counter-adaptation by the antagonist. Additionally, genetic architecture underlying the interaction is often highly complex thus preventing specific adaptive responses. Here, we have employed a reciprocal cross-infection experiment to unravel the adaptive responses of two components of fitness affecting both parties with different complexities of the underlying genetic architecture (i.e. mortality and spore load). Furthermore, our experimental coevolution of hosts (Tribolium castaneum) and parasites (Nosema whitei) included paired replicates of naive hosts from identical genetic backgrounds to allow separation between host- and parasite-specific responses.

Results

In hosts, coevolution led to higher resistance and altered resistance profiles compared to paired control lines. Host genotype × parasite genotype interactions (G_H × G_P) were observed for spore load (the trait of lower genetic complexity), but not for mortality. Overall parasite performance correlated with resistance of its matching host coevolution background reflecting a directional and unspecific response to strength of selection during coevolution. Despite high selective pressures exerted by the obligatory killing parasite, and host- and parasite-specific mortality profiles, no general pattern of local adaptation was observed, but one case of parasite maladaptation was consistently observed on both coevolved and control host populations. In addition, the use of replicate control host populations in the assay revealed one case of host maladaptation and one case of parasite adaptation that was masked by host counter-adaptation, suggesting the presence of complex and probably dynamically changing fitness landscapes.

Conclusions

Our results demonstrate that the use of replicate naive populations can be a useful tool to differentiate between host and parasite adaptation in complex and dynamic fitness landscapes. The absence of clear local adaptation patterns during coevolution with a sexual host showing a complex genetic architecture for resistance suggests that directional selection for generality may be more important attributes of host-parasite coevolution than commonly assumed.     

Development of 15 polymorphic microsatellite markers in the Atlantic tarpon (Megalops atlanticus) for capture-recapture studies

Fifteen polymorphic microsatellite DNA loci for Atlantic tarpon, Megalops atlanticus, were isolated by using PIMA, a polymerase chain reaction‐based technique. The number of alleles at each locus ranged from two to 24 (mean = 7.7) in 65 specimens from Tampa Bay, Florida. Observed and expected heterozygosities ranged from 0.27 to 0.92 (mean = 0.60) and from 0.28 to 0.95 (mean = 0.62), respectively. Genotypes at one locus deviated significantly from Hardy–Weinberg equilibrium. In exact tests for genotypic disequilibrium, there was no evidence of associations between any pair of loci. Overall, loci were well resolved and highly polymorphic, confirming their suitability for DNA fingerprinting applications and other genetic studies.     

Why are birds’ eggs speckled?

Birds are unique in laying eggs with pigmented shells, but for most species (e.g. most passerines, which lay white eggs speckled with reddish spots of protoporphyrin) the pigments’ function is unknown. We studied a bird population at a geologically variable site, and considered a hitherto untested hypothesis: that protoporphyrin pigments might compensate for reduced eggshell‐thickness (caused partly by calcium deficiency), which is known to reduce eggshell‐strength and increase eggshell‐permeability. We found that pigment spots specifically demarcated thinner areas of shell, with darker spots marking yet thinner shell than paler spots. Variation in pigmentation was thus associated with variation in shell thickness both within and between clutches, so accounting for the eggshell's characteristic spot patterns. Geological variability at this site has resulted in a great range of calcium availability and, as predicted by the hypothesis, variation in calcium availability was found to affect between‐clutch variation in both eggshell‐mass (+) and pigmentation characteristics (?). We suggest a physiological mechanism and some important implications of these findings.     

Rhesus macaques build new social connections after a natural disaster

1. Download : Download high-res image (239KB)
2. Download : Download full-size image

     

First Direct Evidence of Chalcolithic Footwear from the Near Eastern Highlands

In 2008, a well preserved and complete shoe was recovered at the base of a Chalcolithic pit in the cave of Areni-1, Armenia. Here, we discuss the chronology of this find, its archaeological context and its relevance to the study of the evolution of footwear. Two leather samples and one grass sample from the shoe were dated at the Oxford Radiocarbon Accelerator Unit (ORAU). A third leather sample was dated at the University of California-Irvine Accelerator Mass Spectrometry Facility (UCIAMS). The R_Combine function for the three leather samples provides a date range of 3627–3377 Cal BC (95.4% confidence interval) and the calibrated range for the straw is contemporaneous (3627–3377 Cal BC). The shoe was stuffed with loose, unfastened grass (Poaceae) without clear orientation which was more than likely used to maintain the shape of the shoe and/or prepare it for storage. The shoe is 24.5 cm long (European size 37), 7.6 to 10 cm wide, and was made from a single piece of leather that wrapped around the foot. It was worn and shaped to the wearer''s right foot, particularly around the heel and hallux where the highest pressure is exerted in normal gait. The Chalcolithic shoe provides solid evidence for the use of footwear among Old World populations at least since the Chalcolithic. Other 4^th millennium discoveries of shoes (Italian and Swiss Alps), and sandals (Southern Israel) indicate that more than one type of footwear existed during the 4^th millennium BC, and that we should expect to discover more regional variations in the manufacturing and style of shoes where preservation conditions permit.     

Evidence for diet partitioning among three coexisting native freshwater fishes in South Africa's Cape Fold Ecoregion

The partitioning of limited resources commonly explains how different species can coexist within the same ecological community. In this 2010 study, the diets of three coexisting freshwater fishes (Cape galaxias Galaxias zebratus, n = 27; Cape kurper Sandelia capensis, n = 60; Breede River redfin Pseudobarbus burchelli, n = 77) were characterised and compared in three headwater streams in South Africa's Cape Fold Ecoregion using gut contents and stable isotope analyses. These data were analysed to ascertain whether the three species exploit distinct trophic niches. Both approaches provided evidence that these species occupy different trophic niches, though with some overlap. However, dietary differences among sites were not consistent and were probably influenced by site-specific factors like resource availability. Pseudobarbus burchelli had a broader niche breadth at Tierkloof Stream than the other two species, but not at Waaihoek or Tierstel Streams. Our results also suggest that P. burchelli consumed a more omnivorous diet than do the other two species, whereas S. capensis occupied a higher trophic position than the other two species and consumed vertebrates. Our findings suggest that these species occupy non-equivalent feeding niches in Cape Fold Ecoregion headwater streams, and that diet partitioning might facilitate their coexistence in these systems.     

Formation of amyloid fibrils by peptides derived from the bacterial cold shock protein CspB.

Three peptides covering the sequence regions corresponding to the first two (CspB-1), the first three (CspB-2), and the last two (CspB-3) beta-strands of CspB, the major cold shock protein of Bacillus subtilis, have been synthesized and analyzed for their conformations in solution and for their precipitation behavior. The peptides are nearly insoluble in water, but highly soluble in aqueous solutions containing 50% acetonitrile (pH 4.0). Upon shifts of the solvent condition toward lower or higher acetonitrile concentrations, the peptides all form fibrils resembling those observed in amyloid associated diseases. These fibrils have been identified and characterized by electron microscopy, binding of the dye congo red, and X-ray fiber diffraction. Characterization of the peptides in solution by circular dichroism and NMR spectroscopy shows that the formation of these fibrils does not require specific preformed secondary structure in the solution state species. While the majority of the soluble fraction of each peptide is monomeric and unstructured, different types of structures including alpha-helical, beta-sheet, and random coil conformations are observed under conditions that eventually lead to fibril formation. We conclude that the absence of tertiary contacts under solution conditions where binding interactions between peptide units are still favorable is a crucial requirement for amyloid formation. Thus, fragmentation of a sequence, like partial chemical denaturation or mutation, can enhance the capacity of specific protein sequences to form such fibrils.