Integrating functional genomics data using maximum likelihood based simultaneous component analysis

Background  

In contemporary biology, complex biological processes are increasingly studied by collecting and analyzing measurements of the same entities that are collected with different analytical platforms. Such data comprise a number of data blocks that are coupled via a common mode. The goal of collecting this type of data is to discover biological mechanisms that underlie the behavior of the variables in the different data blocks. The simultaneous component analysis (SCA) family of data analysis methods is suited for this task. However, a SCA may be hampered by the data blocks being subjected to different amounts of measurement error, or noise. To unveil the true mechanisms underlying the data, it could be fruitful to take noise heterogeneity into consideration in the data analysis. Maximum likelihood based SCA (MxLSCA-P) was developed for this purpose. In a previous simulation study it outperformed normal SCA-P. This previous study, however, did not mimic in many respects typical functional genomics data sets, such as, data blocks coupled via the experimental mode, more variables than experimental units, and medium to high correlations between variables. Here, we present a new simulation study in which the usefulness of MxLSCA-P compared to ordinary SCA-P is evaluated within a typical functional genomics setting. Subsequently, the performance of the two methods is evaluated by analysis of a real life Escherichia coli metabolomics data set.     

Detection of MET mRNA in gastric cancer in situ. Comparison with immunohistochemistry and sandwich immunoassays

Determination of predictive biomarkers by immunohistochemistry (IHC) relies on antibodies with high selectivity. RNA in situ hybridization (RNA ISH) may be used to confirm IHC and may potentially replace it if suitable antibodies are not available or are insufficiently selective to discriminate closely related protein isoforms. We validated RNA ISH as specificity control for IHC and as a potential alternative method for selecting patients for treatment with MET inhibitors. MET, the HGF receptor, is encoded by the MET proto-oncogene that may be activated by mutation or amplification. MET expression and activity were tested in a panel of control cell lines. MET could be detected in formalin fixed paraffin, embedded (FFPE) samples by IHC and RNA ISH, and this was confirmed by sandwich immunoassays of fresh frozen samples. Gastric cancer cell lines with high MET expression and phosphorylation of tyrosine-1349 respond to the MET inhibitor, BAY-853474. High expression and phosphorylation of MET is a predictive biomarker for response to MET inhibitors. We then analyzed MET expression and activity in a matched set of FFPE vs. fresh frozen tumor samples consisting of 20 cases of gastric cancer. Two of 20 clinical samples investigated exhibited high MET expression with RNA ISH and IHC. Both cases were shown by sandwich immunoassays to exhibits strong functional activity. Expression levels and functional activity in these two cases were in a range that predicted response to treatment. Our findings indicate that owing to its high selectivity, RNA ISH can be used to confirm findings obtained by IHC and potentially may replace IHC for certain targets if no suitable antibodies are available. RNA ISH is a valid platform for testing predictive biomarkers for patient selection.     

Nonlinear relationships among evolutionary rates identify regions of functional divergence in heat-shock protein 70 genes

The neutral theory predicts that, in comparisons among related genes, the number of amino acid replacements per site in a given gene region should be a linear function of that in another region of the same gene, unless the genes have diverged functionally in one region. Therefore, nonlinearity of this relationship can be used to identify regions of possible functional divergence among members of a multigene family. This method of analysis was applied to members of the heat-shock protein 70 (HSP70) gene family, which encode highly conserved ATP- dependent chaperone proteins found in all organisms. A nonlinear relationship was found between the rate of amino acid replacement in the conserved IA domain of the ATPase portion of the molecule and that in other ATPase domains and the peptide-binding domain. These results suggest that genes in the HSP70 subfamily C (dnaK of bacteria and SSC1 of yeast) may have diverged functionally from other subfamilies in the ATPase domains, especially IIB, whereas SSB1 of yeast has diverged markedly in the peptide-binding domain. Functional divergence within these regions is consistent with what is known about functional differences between the HSP70 subfamilies in yeast.      

Blebs produced by actin–myosin contraction during apoptosis release damage-associated molecular pattern proteins before secondary necrosis occurs

Apoptosis is a fundamental homeostatic mechanism essential for the normal growth, development and maintenance of every tissue and organ. Dying cells have been defined as apoptotic by distinguishing features, including cell contraction, nuclear fragmentation, blebbing, apoptotic body formation and maintenance of intact cellular membranes to prevent massive protein release and consequent inflammation. We now show that during early apoptosis limited membrane permeabilization occurs in blebs and apoptotic bodies, which allows release of proteins that may affect the proximal microenvironment before the catastrophic loss of membrane integrity during secondary necrosis. Blebbing, apoptotic body formation and protein release during early apoptosis are dependent on ROCK and myosin ATPase activity to drive actomyosin contraction. We identified 231 proteins released from actomyosin contraction-dependent blebs and apoptotic bodies by adapted SILAC (stable isotope labeling with amino acids in cell culture) combined with mass spectrometry analysis. The most enriched proteins released were the nucleosomal histones, which have previously been identified as damage-associated molecular pattern proteins (DAMPs) that can initiate sterile inflammatory responses. These results indicate that limited membrane permeabilization occurs in blebs and apoptotic bodies before secondary necrosis, leading to acute and localized release of immunomodulatory proteins during the early phase of active apoptotic membrane blebbing. Therefore, the shift from apoptosis to secondary necrosis is more graded than a simple binary switch, with the membrane permeabilization of apoptotic bodies and consequent limited release of DAMPs contributing to the transition between these states.     

Drivers of beta diversity of macroinvertebrate communities in tropical forest streams

1. There has recently been increasing interest in patterns of beta diversity but we still lack a comprehensive understanding of these patterns in various regions (e.g. the tropics), ecosystems (e.g. streams) and organism groups (e.g. invertebrates). 2. Our aim was to investigate the patterns of beta diversity of stream macroinvertebrates in relation to key environmental (i.e. stream size, pH and habitat degradation) and geographical variables (i.e. latitude, longitude, altitude) in a tropical region. We surveyed a total of 8–10 riffle sites in each of 34 streams (altogether 337 riffle sites were sampled) in Peninsular Malaysia to examine variation in macroinvertebrate community composition at within‐stream and among‐stream scales. 3. Based on test of homogeneity of dispersion, we found that the streams studied differed significantly in within‐stream variation in community composition (i.e. among‐site variation of within stream beta diversity). The patterns were similar based on Bray–Curtis coefficient on abundance data, Sorensen coefficient on presence–absence data and Simpson coefficient on presence–absence data. We also found that within‐stream beta diversity was significantly related to stream size, pH and latitude, with each of these variables individually accounting for around 20% of the variation in beta diversity in simple regressions, while the total variation explained by the three significant variables amounted to around 50% in multiple regressions. By contrast, habitat degradation, longitude and altitude were not significantly related to beta diversity. We also found that the factor drainage basin accounted for much of the variation in beta diversity in general linear models, suppressing the effects of environmental variables. 4. We concluded that within‐stream beta diversity is mainly related to a combination of the identity of a drainage basin and stream environmental factors. Our findings provide important background for stream environmental assessment and conservation planning by emphasising that (i) macroinvertebrate communities within streams are not homogeneous, but show considerable beta diversity, (ii) streams differ in their degree of within‐stream beta diversity, (iii) stream size and water pH should be considered in applied contexts related to within‐stream beta diversity and (iv) historical effects may be different in different drainage basins and may affect present‐day patterns of within‐stream beta diversity.     

Studies on sodium bypass flow in lateral rootless mutants lrt1 and lrt2, and crown rootless mutant crl1 of rice (Oryza sativa L.)

An apoplastic pathway, the so‐called bypass flow, is important for Na⁺ uptake in rice (Oryza sativa L.) under saline conditions; however, the precise site of entry is not yet known. We report the results of our test of the hypothesis that bypass flow of Na⁺ in rice occurs at the site where lateral roots emerge from the main roots. We investigated Na⁺ uptake and bypass flow in lateral rootless mutants (lrt1, lrt2), a crown rootless mutant (crl1), their wild types (Oochikara, Nipponbare and Taichung 65, respectively) and in seedlings of rice cv. IR36. The results showed that shoot Na⁺ concentration in lrt1, lrt2 and crl1 was lower (by 20–23%) than that of their wild types. In contrast, the bypass flow quantified using trisodium‐8‐hydroxy‐1,3,6‐pyrenetrisulphonic acid (PTS) was significantly increased in the mutants, from an average of 1.1% in the wild types to 3.2% in the mutants. Similarly, bypass flow in shoots of IR36 where the number of lateral and crown roots had been reduced through physical and hormonal manipulations was dramatically increased (from 5.6 to 12.5%) as compared to the controls. The results suggest that the path of bypass flow in rice is not at the sites of lateral root emergence.     

促胰液素和胆囊收缩素族激素对豚鼠肝胆汁分泌的影响

用具备胃瘘和胆瘘的豚鼠于人工维持胆汁酸池恒定的条件下,观察促胰液素(SEC)和胆囊收缩素(CCK)族激素[包括雨蛙肽(CAE)、五肽胃泌素(G5)和内源性CCK]对肝胆汁分泌的影响及其相互作用。结果表明:静脉灌注SEC、CAE或肠内灌注左旋苯丙氨酸(L-PHE,促内源性CCK释放剂)后,胆汁流量、胆汁HCO~-_3和Cl~-排出量均显著增多,并呈剂量-效应关系,但静脉注射G5则无利胆效应。在恒速灌注SEC的背景下,CAE或CCK对胆汁HCO~-_3排出的效应分别大于它们单独给予时的效应(P<0.05或P<0.01)。这些激素对胆汁酸的排出量均无影响。上述结果表明,SEC,CAE和内源性CCK均有利胆作用,所刺激的肝胆汁属于不依赖胆汁酸部分。G5则无利胆效应。对胆汁中HCO~-_3的排出,SEC与CAE或内源性CCK间有相互加强作用。     

Phylogeny of Six Genera of the Subclass Haptoria (Ciliophora,Litostomatea) Inferred from Sequences of the Gene Coding for Small Subunit Ribosomal RNA

ABSTRACT. The small subunit ribosomal RNA genes of nine species belonging to six genera of litostome ciliates, namely Amphileptus aeschtae, Chaenea teres, Chaenea vorax, Lacrymaria marina, Litonotus paracygnus, Loxophyllum sp.‐GD‐070419, Loxophyllum jini, Loxophyllum rostratum, and Phialina salinarum, were sequenced for the first time. Phylogenetic trees were constructed using different methods to assess the inter‐ and intra‐generic relationships of haptorians, of which Chaenea, Lacrymaria, Litonotus, and Phialina were analyzed for the first time based on molecular data. Monophyly of the order Pleurostomatida was strongly confirmed, and the two existing families of pleurostomatids, created on the basis of morphology, were confirmed by molecular evidence. Within the Pleurostomatida, Siroloxophyllum utriculariae occupied a well‐supported position basal to the Loxophyllum clade, supporting the separation of these genera from one another. Both the subclass Haptoria and the order Haptorida were partially unresolved, possibly paraphyletic assemblages of taxa in all analyses, creating doubts about the traditional placement of some haptorid taxa. The existing sequence of L. rostratum in GenBank (DQ411864) was conspicuously different from that of the isolate from Qingdao, China sequenced in the present work, indicating that they are different species. The isolate from Qingdao was verified as L. rostratum by morphological analysis, and the published morphology of existing GenBank record of L. rostratum is different from it. Based on both morphological and molecular evidence, the latter may be congeneric with an undescribed species of Loxophyllum from Guangdong Province, China.