libSRES: a C library for stochastic ranking evolution strategy for parameter estimation

SUMMARY: Estimation of kinetic parameters in a biochemical pathway or network represents a common problem in systems studies of biological processes. We have implemented a C library, named libSRES, to facilitate a fast implementation of computer software for study of non-linear biochemical pathways. This library implements a (mu, lambda)-ES evolutionary optimization algorithm that uses stochastic ranking as the constraint handling technique. Considering the amount of computing time it might require to solve a parameter-estimation problem, an MPI version of libSRES is provided for parallel implementation, as well as a simple user interface. libSRES is freely available and could be used directly in any C program as a library function. We have extensively tested the performance of libSRES on various pathway parameter-estimation problems and found its performance to be satisfactory. AVAILABILITY: The source code (in C) is free for academic users at http://csbl.bmb.uga.edu/~jix/science/libSRES/     

天敌对油菜蚜虫的综合控制作用分析

本研究对不同肥力水平两个品系油菜上蚜虫自然种群和阻隔天敌条件下的蚜虫种群采用蚜虫增长指数的超越函数定积分方法评价天敌的综合控制作用。这是评价自然条件下天敌控制作用的一种好方法。     

The lifestyle transition of Arthrobotrys oligospora is mediated by microRNA-like RNAs

The lifestyle transition of fungi, defined as switching from taking organic material as nutrients to pathogens, is a fundamental phenomenon in nature. However, the mechanisms of such transition remain largely unknown. Here we show microRNA-like RNAs(milRNAs) play a key role in fungal lifestyle transition for the first time. We identified milRNAs by small RNA sequencing in Arthrobotrys oligospora, a known nematode-trapping fungus. Among them, 7 highly expressed milRNAs were confirmed by northern-blot analysis. Knocking out two milRNAs significantly decreased A. oligospora's ability to switch lifestyles. We further identified that two of these milRNAs were associated with argonaute protein QDE-2 by RNA-immunoprecipitation(RIP) analysis. Three of the predicted target genes of milRNAs were found in immunoprecipitation(IP)products of QDE-2. Disruption of argonaute gene qde-2 also led to serious defects in lifestyle transition. Interestingly, knocking out individual milRNAs or qde-2 lead to diverse responses under different conditions, and qde-2 itself may be targeted by the milRNAs. Collectively, it indicates the lifestyle transition of fungi is mediated by milRNAs through RNA interference(RNAi)machinery, revealing the wide existence of miRNAs in fungi kingdom and providing new insights into understanding the adaptation of fungi from scavengers to predators and the mechanisms underlying fungal infections.     

DA and Xiao-two giant and composite LTR-retrotransposon-like elements identified in the human genome

We discovered two new complex elements while studying large genomic rearrangements and segmental duplications in the human genome. Both resemble bacterial composite DNA transposon Tn9, consisting of a core flanked by mobile elements, except that the flanking element is not a DNA transposon but instead is long terminal repeat retrotransposon-like with human endogenous retrovirus and satellite sequences. Based on the core size, we named them Xiao ( approximately 30 kb) and DA ( approximately 280 kb), meaning small and big, respectively, in Chinese. Xiao originated from a 19p region encoding olfactory receptor 7E members after the human/ape divergence from Old World monkeys, while DA likely evolved from a Xiao by inserting approximately 200 kb of chimeric sequence from 16p and 21q into the Xiao core, resulting in a target site duplication of 3.4 kb. DA/Xiao was identified in 30 loci on 12 chromosomes, and only DAs mediated intrachromosomal rearrangements, based on our reconstructed human-mouse-rat ancestral genome and the rhesus macaque genome.     

HMMGEP: clustering gene expression data using hidden Markov models

SUMMARY: The package HMMGEP performs cluster analysis on gene expression data using hidden Markov models. AVAILABILITY: HMMGEP, including the source code, documentation and sample data files, is available at http://www.bioinfo.tsinghua.edu.cn:8080/~rich/hmmgep_download/index.html.     

TWO INNEXINS OF Spodoptera litura INFLUENCES HEMICHANNEL AND GAP JUNCTION FUNCTIONS IN CELLULAR IMMUNE RESPONSES

Insect cellular immune responses include encapsulation, nodule formation, and phagocytosis. Hemichannels and gap junctions are involved in these cellular actions. Innexins (Inxs: analogous to the vertebrate connexins) form hemichannels and gap junctions, but the molecular mechanisms underlying their biology is still unclear. In this article, we reported a steady‐state level of Inxs (SpliInxs) in hemocytes of Spodoptera litura, which formed nonfunctional hemichannels on the cell surface to maintain normal metabolism. We also reported that two innnexins (SpliInx2 and SpliInx3) were expressed significantly higher in hemocytes compared to other tissues, suggesting that they play important roles in hemocytes. Amino acid analysis found that two cysteine residues in two extracellular loops provided the capability for SpliInx2 and SpliInx3 hemichannels to dock into gap junctions. Western blotting demonstrated that both extracellular and intracellular loops of SpliInx3 and the extracellular loops of SpliInx2 might undergo posttranslational modification during the formation of a steady‐state hemichannel. During hemichannel formation, SpliInx2 presented as one isoform, while SpliInx3 presented as three isoforms. These results provide fundamental knowledge for further study of how steady‐state levels of SpliInxs are dynamically adjusted to perform cellular immune responses under immune challenge.     

Frequent Alteration of the Tumor Suppressor Gene APC in Sporadic Canine Colorectal Tumors

Sporadic canine colorectal cancers (CRCs) should make excellent models for studying the corresponding human cancers. To molecularly characterize canine CRC, we investigated exonic sequence mutations of adenomatous polyposis coli (APC), the best known tumor suppressor gene of human CRC, in 23 sporadic canine colorectal tumors, including 8 adenomas and 15 adenocarcinomas, via exon-resequencing analysis. As a comparison, we also performed the same sequencing analysis on 10 other genes, either located at human 5q22 (the same locus as APC) or 18q21 (also frequently altered in human CRC), or known to play a role in human carcinogenesis. We noted that APC was the most significantly mutated gene in both canine adenomas and adenocarcinomas among the 11 genes examined. Significantly, we detected large deletions of ≥10 bases, many clustered near the mutation cluster region, as well as single or two base deletions in ∼70% canine tumors of both subtypes. These observations indicate that like in the human, APC is also frequently altered in sporadic colorectal tumors in the dog and its alteration is an early event in canine colorectal tumorigenesis. Our study provides further evidence demonstrating the molecular similarity in pathogenesis between sporadic human and canine CRCs. This work, along with our previous copy number abnormality study, supports that sporadic canine CRCs are valid models of human CRCs at the molecular level.     

Mining gene expression data using a novel approach based on hidden Markov models

In this work we have developed a new framework for microarray gene expression data analysis. This framework is based on hidden Markov models. We have benchmarked the performance of this probability model-based clustering algorithm on several gene expression datasets for which external evaluation criteria were available. The results showed that this approach could produce clusters of quality comparable to two prevalent clustering algorithms, but with the major advantage of determining the number of clusters. We have also applied this algorithm to analyze published data of yeast cell cycle gene expression and found it able to successfully dig out biologically meaningful gene groups. In addition, this algorithm can also find correlation between different functional groups and distinguish between function genes and regulation genes, which is helpful to construct a network describing particular biological associations. Currently, this method is limited to time series data. Supplementary materials are available at http://www.bioinfo.tsinghua.edu.cn/~rich/hmmgep_supp/.