A linear programming approach for estimating the structure of a sparse linear genetic network from transcript profiling data

Background: A genetic network can be represented as a directed graph in which a node corresponds to a gene and a directed edge specifies the direction of influence of one gene on another. The reconstruction of such networks from transcript profiling data remains an important yet challenging endeavor. A transcript profile specifies the abundances of many genes in a biological sample of interest. Prevailing strategies for learning the structure of a genetic network from high-dimensional transcript profiling data assume sparsity and linearity. Many methods consider relatively small directed graphs, inferring graphs with up to a few hundred nodes. This work examines large undirected graphs representations of genetic networks, graphs with many thousands of nodes where an undirected edge between two nodes does not indicate the direction of influence, and the problem of estimating the structure of such a sparse linear genetic network (SLGN) from transcript profiling data. Results: The structure learning task is cast as a sparse linear regression problem which is then posed as a LASSO (l1-constrained fitting) problem and solved finally by formulating a Linear Program (LP). A bound on the Generalization Error of this approach is given in terms of the Leave-One-Out Error. The accuracy and utility of LP-SLGNs is assessed quantitatively and qualitatively using simulated and real data. The Dialogue for Reverse Engineering Assessments and Methods (DREAM) initiative provides gold standard data sets and evaluation metrics that enable and facilitate the comparison of algorithms for deducing the structure of networks. The structures of LP-SLGNs estimated from the INSILICO1, INSILICO2 and INSILICO3 simulated DREAM2 data sets are comparable to those proposed by the first and/or second ranked teams in the DREAM2 competition. The structures of LP-SLGNs estimated from two published Saccharomyces cerevisae cell cycle transcript profiling data sets capture known regulatory associations. In each S. cerevisiae LP-SLGN, the number of nodes with a particular degree follows an approximate power law suggesting that its degree distributions is similar to that observed in real-world networks. Inspection of these LP-SLGNs suggests biological hypotheses amenable to experimental verification. Conclusion: A statistically robust and computationally efficient LP-based method for estimating the topology of a large sparse undirected graph from high-dimensional data yields representations of genetic networks that are biologically plausible and useful abstractions of the structures of real genetic networks. Analysis of the statistical and topological properties of learned LP-SLGNs may have practical value; for example, genes with high random walk betweenness, a measure of the centrality of a node in a graph, are good candidates for intervention studies and hence integrated computational – experimental investigations designed to infer more realistic and sophisticated probabilistic directed graphical model representations of genetic networks. The LP-based solutions of the sparse linear regression problem described here may provide a method for learning the structure of transcription factor networks from transcript profiling and transcription factor binding motif data.

Structural studies of analgesics and their interactions. VII. Stereoisomerism and disorder in the structure of oxyphenbutazone monohydrate

Oxyphenbutazone, C19H20N203, a metabolite and perhaps the active form of phenylbutazone, is a widely used non-narcotic analgesic and anti-inflammatory pyrazolidinedione derivative. The monohydrate of the compound crystallizes in the triclinic space group Pi with two molecules in a unit cell of dimensions a -- 9.491 (4), b = 10.261 (5), c = 11.036 (3)A and ¢~ = 72.2 (1), fl = 64.3 (1), 7 = 73.0 (1) °. The structure was solved by direct methods and refined to an R value of 0.107 for 1498 observed reflections. The butyl group in the molecule is disordered. The hydroxyl group occupies two sites with unequal occupancies. On account of the asymmetry at the two N atoms and one of the C atoms in the central five-membered ring, the molecule can exist in eight isomeric states, of which four are sterically unfavourable. The disorder in the position of the hydroxyl group can be readily explained on the basis of the existence, with unequal abundances, of all four sterically favourable isomers.The bond lengths and angles in the molecule are similar to those in phenylbutazone. The crystal structure is stabilized by van der Waals interactions, and O-H... O hydrogen bonds involving the carbonyl and the hydroxyl groups as well as a water molecule.

Proteomics of Trypanosoma evansi Infection in Rodents

Background: Trypanosoma evansi infections, commonly called 'surra', cause significant economic losses to livestock industry. While this infection is mainly restricted to large animals such as camels, donkeys and equines, recent reports indicate their ability to infect humans. There are no World Animal Health Organization (WAHO) prescribed diagnostic tests or vaccines available against this disease and the available drugs show significant toxicity. There is an urgent need to develop improved methods of diagnosis and control measures for this disease. Unlike its related human parasites T. brucei and T. cruzi whose genomes have been fully sequenced T. evansi genome sequence remains unavailable and very little efforts are being made to develop improved methods of prevention, diagnosis and treatment. With a view to identify potential diagnostic markers and drug targets we have studied the clinical proteome of T. evansi infection using mass spectrometry (MS).Methodology/Principal Findings: Using shot-gun proteomic approach involving nano-lc Quadrupole Time Of Flight (QTOF) mass spectrometry we have identified over 160 proteins expressed by T. evansi in mice infected with camel isolate. Homology driven searches for protein identification from MS/MS data led to most of the matches arising from related Trypanosoma species. Proteins identified belonged to various functional categories including metabolic enzymes; DNA metabolism; transcription; translation as well as cell-cell communication and signal transduction. TCA cycle enzymes were strikingly missing, possibly suggesting their low abundances. The clinical proteome revealed the presence of known and potential drug targets such as oligopeptidases, kinases, cysteine proteases and more.Conclusions/Significance: Previous proteomic studies on Trypanosomal infections, including human parasites T. brucei and T. cruzi, have been carried out from lab grown cultures. For T. evansi infection this is indeed the first ever proteomic study reported thus far. In addition to providing a glimpse into the biology of this neglected disease, our study is the first step towards identification of diagnostic biomarkers, novel drug targets as well as potential vaccine candidates to fight against T. evansi infections.

Long term monitering of vegetation in a tropical deciduous forest in Mudumalai,southern India

As part of an international network of large plots to study tropical vegetation dynamics on a long-term basis, a 50-hectare permanent plot was set up during 1988-89 in the deciduous forests of Mudumalai, southern India. Within this plot 25,929 living woody plants (71 species) above 1 cm DBH (diameter at breast height) were identified, measured, tagged and mapped. Species abundances corresponded to the characteristic log-normal distribution. The four most abundant species (Kydia calycina, Lagerstroemia microcarpa, Terminalia crenulata and Helicteres isora) constituted nearly 56% of total stems, while seven species were represented by only one individual each in the plot. Variance/mean ratios of density showed most species to have clumped distributions. The population declined overall by 14% during the first two years, largely due to elephant and fire-mediated damage to Kydia calycina and Helicteres isora. In this article we discuss the need for large plots to study vegetation dynamics.

A Multi-layer Perceptron based Non-linear Mixture Model to estimate class abundance from mixed pixels

Sub-pixel classification is essential for the successful description of many land cover (LC) features with spatial resolution less than the size of the image pixels. A commonly used approach for sub-pixel classification is linear mixture models (LMM). Even though, LMM have shown acceptable results, pragmatically, linear mixtures do not exist. A non-linear mixture model, therefore, may better describe the resultant mixture spectra for endmember (pure pixel) distribution. In this paper, we propose a new methodology for inferring LC fractions by a process called automatic linear-nonlinear mixture model (AL-NLMM). AL-NLMM is a three step process where the endmembers are first derived from an automated algorithm. These endmembers are used by the LMM in the second step that provides abundance estimation in a linear fashion. Finally, the abundance values along with the training samples representing the actual proportions are fed to multi-layer perceptron (MLP) architecture as input to train the neurons which further refines the abundance estimates to account for the non-linear nature of the mixing classes of interest. AL-NLMM is validated on computer simulated hyperspectral data of 200 bands. Validation of the output showed overall RMSE of 0.0089±0.0022 with LMM and 0.0030±0.0001 with the MLP based AL-NLMM, when compared to actual class proportions indicating that individual class abundances obtained from AL-NLMM are very close to the real observations.

Hybrid Bayesian Classifier for Improved Classification Accuracy

The widely used Bayesian classifier is based on the assumption of equal prior probabilities for all the classes. However, inclusion of equal prior probabilities may not guarantee high classification accuracy for the individual classes. Here, we propose a novel technique-Hybrid Bayesian Classifier (HBC)-where the class prior probabilities are determined by unmixing a supplemental low spatial-high spectral resolution multispectral (MS) data that are assigned to every pixel in a high spatial-low spectral resolution MS data in Bayesian classification. This is demonstrated with two separate experiments-first, class abundances are estimated per pixel by unmixing Moderate Resolution Imaging Spectroradiometer data to be used as prior probabilities, while posterior probabilities are determined from the training data obtained from ground. These have been used for classifying the Indian Remote Sensing Satellite LISS-III MS data through Bayesian classifier. In the second experiment, abundances obtained by unmixing Landsat Enhanced Thematic Mapper Plus are used as priors, and posterior probabilities are determined from the ground data to classify IKONOS MS images through Bayesian classifier. The results indicated that HBC systematically exploited the information from two image sources, improving the overall accuracy of LISS-III MS classification by 6% and IKONOS MS classification by 9%. Inclusion of prior probabilities increased the average producer's and user's accuracies by 5.5% and 6.5% in case of LISS-III MS with six classes and 12.5% and 5.4% in IKONOS MS for five classes considered.

Features of a unique intronless cluster of class I small heat shock protein genes in tandem with box C/D snoRNA genes on chromosome 6 in tomato (Solanum lycopersicum)

Physical clustering of genes has been shown in plants; however, little is known about gene clusters that have different functions, particularly those expressed in the tomato fruit. A class I 17.6 small heat shock protein (Sl17.6 shsp) gene was cloned and used as a probe to screen a tomato (Solanum lycopersicum) genomic library. An 8.3-kb genomic fragment was isolated and its DNA sequence determined. Analysis of the genomic fragment identified intronless open reading frames of three class I shsp genes (Sl17.6, Sl20.0, and Sl20.1), the Sl17.6 gene flanked by Sl20.1 and Sl20.0, with complete 5' and 3' UTRs. Upstream of the Sl20.0 shsp, and within the shsp gene cluster, resides a box C/D snoRNA cluster made of SlsnoR12.1 and SlU24a. Characteristic C and D, and C' and D', boxes are conserved in SlsnoR12.1 and SlU24a while the upstream flanking region of SlsnoR12.1 carries TATA box 1, homol-E and homol-D box-like cis sequences, TM6 promoter, and an uncharacterized tomato EST. Molecular phylogenetic analysis revealed that this particular arrangement of shsps is conserved in tomato genome but is distinct from other species. The intronless genomic sequence is decorated with cis elements previously shown to be responsive to cues from plant hormones, dehydration, cold, heat, and MYC/MYB and WRKY71 transcription factors. Chromosomal mapping localized the tomato genomic sequence on the short arm of chromosome 6 in the introgression line (IL) 6-3. Quantitative polymerase chain reaction analysis of gene cluster members revealed differential expression during ripening of tomato fruit, and relatively different abundances in other plant parts.

The Contribution of Rare Species to Community Phylogenetic Diversity across a Global Network of Forest Plots

Niche differentiation has been proposed as an explanation for rarity in species assemblages. To test this hypothesis requires quantifying the ecological similarity of species. This similarity can potentially be estimated by using phylogenetic relatedness. In this study, we predicted that if niche differentiation does explain the co-occurrence of rare and common species, then rare species should contribute greatly to the overall community phylogenetic diversity (PD), abundance will have phylogenetic signal, and common and rare species will be phylogenetically dissimilar. We tested these predictions by developing a novel method that integrates species rank abundance distributions with phylogenetic trees and trend analyses, to examine the relative contribution of individual species to the overall community PD. We then supplement this approach with analyses of phylogenetic signal in abundances and measures of phylogenetic similarity within and between rare and common species groups. We applied this analytical approach to 15 long-term temperate and tropical forest dynamics plots from around the world. We show that the niche differentiation hypothesis is supported in six of the nine gap-dominated forests but is rejected in the six disturbance-dominated and three gap-dominated forests. We also show that the three metrics utilized in this study each provide unique but corroborating information regarding the phylogenetic distribution of rarity in communities.

Nucleosynthesis in the accretion disks of type II collapsars

We investigate nucleosynthesis inside the gamma-ray burst (GRB) accretion disks formed by the Type II collapsars. In these collapsars, the core collapse of massive stars first leads to the formation of a proto-neutron star. After that, an outward moving shock triggers a successful supernova. However, the supernova ejecta lacks momentum and within a few seconds the newly formed neutron star gets transformed to a stellar mass black hole via massive fallback. The hydrodynamics of such an accretion disk formed from the fallback material of the supernova ejecta has been studied extensively in the past. We use these well-established hydrodynamic models for our accretion disk in order to understand nucleosynthesis, which is mainly advection dominated in the outer regions. Neutrino cooling becomes important in the inner disk where the temperature and density are higher. The higher the accretion rate (M) over dot is, the higher the density and temperature are in the disks. We deal with accretion disks with relatively low accretion rates: 0.001 M-circle dot s(-1) less than or similar to (M) over dot less than or similar to 0.01 M-circle dot s(-1) and hence these disks are predominantly advection dominated. We use He-rich and Si-rich abundances as the initial condition of nucleosynthesis at the outer disk, and being equipped with the disk hydrodynamics and the nuclear network code, we study the abundance evolution as matter inflows and falls into the central object. We investigate the variation in the nucleosynthesis products in the disk with the change in the initial abundance at the outer disk and also with the change in the mass accretion rate. We report the synthesis of several unusual nuclei like P-31, K-39, Sc-43, Cl-35 and various isotopes of titanium, vanadium, chromium, manganese and copper. We also confirm that isotopes of iron, cobalt, nickel, argon, calcium, sulphur and silicon get synthesized in the disk, as shown by previous authors. Much of these heavy elements thus synthesized are ejected from the disk via outflows and hence they should leave their signature in observed data.

Chromium isotope variations (delta Cr-53/52) in mantle-derived sources and their weathering products: Implications for environmental studies and the evolution of delta Cr-53/52 in the Earth's mantle over geologic time

Here we report chromium isotope compositions, expressed as delta Cr-53/ 52 in per mil (&) relative to NIST 979, measured in selected Cr-rich minerals and rocks formed by the primary magmatic as well as the secondary metamorphic and weathering processes. The main objectives of this study were: (i) to further constrain the isotope composition of the Earth's mantle Cr inventory and its possible variation during geological history, based on the analysis of globally distributed and stratigraphically constrained mantle-derived chromites; and (ii) to investigate the magnitude and systematics of Cr isotope fractionation during oxidative weathering and secondary alteration (i. e., hydration, serpentinization) of the magmatic Cr sources. Specifically, we analyzed delta Cr-53/ 52 in a set of globally distributed mantle-derived chromites (FeMgCr2O4, n = 30) collected from various locations in Europe, Asia, Africa and South America, and our results confirm that a chromite-hosted Earth's mantle Cr inventory is uniform at - 0.079 +/- 0.129& (2SD), which we named here as a ` canonical' mantle d 53/ 52 Cr signature. Furthermore our dataset of stratigraphically constrained chromites, whose crystallization ages cover most of the Earth's geological history, indicate that the bulk Cr isotope composition of the chromite-hosted mantle inventory has remained uniform, within about +/- 0.100&, since at least the Early Archean times (similar to 3500 million years ago, Ma). To investigate the systematics of Cr isotope fractionation associated with alteration processes we analyzed a number of secondary Cr-rich minerals and variably altered ultramafic rocks (i. e., serpentinized harzburgites, lherzolites) that revealed large positive delta Cr-53/ 52 anomalies that are systematically shifted to higher values with an increasing degree of alteration and serpentinization. The degree of aqueous alteration and serpentinization was quantified by the abundances of fluid-mobile (Rb, K) elements, and by the Loss On Ignition (LOI) parameter, which determines the amount of structurally bound water (OH/ H2O) present in secondary hydrated minerals like serpentine. Overall, we observed that altered ultramafic rocks that yielded the highest LOI values, and the lowest amounts of fluid mobile elements, also yielded the heaviest delta Cr-53/ 52 signatures. Therefore, we conclude that secondary alteration (i.e., hydration, serpentinization) of ultramafic rocks in near-surface oxidative environments tend to shift the bulk Cr isotope composition of the weathered products to isotopically heavier values, pointing to a dynamic redox cycling of Cr in the Earth's crustal and near-surface environments. Hence, if validated by future

Loss of frugivore seed dispersal services under climate change

The capacity of species to track shifting climates into the future will strongly influence outcomes for biodiversity under a rapidly changing climate. However, we know remarkably little about the dispersal abilities of most species and how these may be influenced by climate change. Here we show that climate change is projected to substantially reduce the seed dispersal services provided by frugivorous vertebrates in rainforests across the Australian Wet Tropics. Our model projections show reductions in both median and long-distance seed dispersal, which may markedly reduce the capacity of many rainforest plant species to track shifts in suitable habitat under climate change. However, our analyses suggest that active management to maintain the abundances of a small set of important frugivores under climate change could markedly reduce the projected loss of seed dispersal services and facilitate shifting distributions of rainforest plant species.

Contrasting geochemical signatures on land from the Middle and Late Permian extinction events

The end of the Palaeozoic is marked by two mass-extinction events during the Middle Permian (Capitanian) and the Late Permian (Changhsingian). Given similarities between the two events in geochemical signatures, such as large magnitude negative C-13 anomalies, sedimentological signatures such as claystone breccias, and the approximate contemporaneous emplacement of large igneous provinces, many authors have sought a common causal mechanism. Here, a new high-resolution continental record of the Capitanian event from Portal Mountain, Antarctica, is compared with previously published Changhsingian records of geochemical signatures of weathering intensity and palaeoclimatic change. Geochemical means of discriminating sedimentary provenance (Ti/Al, U/Th and La/Ce ratios) all indicate a common provenance for the Portal Mountain sediments and associated palaeosols, so changes spanning the Capitanian extinction represent changes in weathering intensity rather than sediment source. Proxies for weathering intensity chemical index of alteration, W and rare earth element accumulation all decline across the Capitanian extinction event at Portal Mountain, which is in contrast to the increased weathering recorded globally at the Late Permian extinction. Furthermore, palaeoclimatic proxies are consistent with unchanging or cooler climatic conditions throughout the Capitanian event, which contrasts with Changhsingian records that all indicate a significant syn-extinction and post-extinction series of greenhouse warming events. Although both the Capitanian and Changhsingian event records indicate significant redox shifts, palaeosol geochemistry of the Changhsingian event indicates more reducing conditions, whereas the new Capitanian record of reduced trace metal abundances (Cr, Cu, Ni and Ce) indicates more oxidizing conditions. Taken together, the differences in weathering intensity, redox and the lack of evidence for significant climatic change in the new record suggest that the Capitanian mass extinction was not triggered by dyke injection of coal-beds, as in the Changhsingian extinction, and may instead have been triggered directly by the Emeishan large igneous province or by the interaction of Emeishan basalts with platform carbonates.

Off-site impacts of agricultural composting: role of terrestrially derived organic matter in structuring aquatic microbial communities and their metabolic potential

While considered as sustainable and low-cost agricultural amendments, the impacts of organic fertilizers on downstream aquatic microbial communities remain poorly documented. We investigated the quantity and quality of the dissolved organic matter leaching from agricultural soil amended with compost, vermicompost or biochar and assessed their effects on lake microbial communities, in terms of viral and bacterial abundances, community structure and metabolic potential. The addition of compost and vermicompost significantly increased the amount of dissolved organic carbon in the leachate compared with soil alone. Leachates from these additions, either with or without biochar, were highly bioavailable to aquatic microbial communities, although reducing the metabolic potential of the community and harbouring more specific communities. Although not affecting bacterial richness or taxonomic distributions, the specific addition of biochar affected the original lake bacterial communities, resulting in a strongly different community. This could be partly explained by viral burst and converging bacterial abundances throughout the samples. These results underline the necessity to include off-site impacts of agricultural amendments when considering their cascading effect on downstream aquatic ecosystems.