A Modified Entropy-Based Approach for Identifying Gene-Gene Interactions in Case-Control Study

Gene-gene interactions may play an important role in the genetics of a complex disease. Detection and characterization of gene-gene interactions is a challenging issue that has stimulated the development of various statistical methods to address it. In this study, we introduce a method to measure gene interactions using entropy-based statistics from a contingency table of trait and genotype combinations. We also developed an exploration procedure by using graphs. We propose a standardized relative information gain (RIG) measure to evaluate the interactions between single nucleotide polymorphism (SNP) combinations. To identify the k ^th order interactions, contingency tables of trait and genotype combinations of k SNPs are constructed, with which RIGs are calculated. The RIGs are standardized using the mean and standard deviation from the permuted datasets. SNP combinations yielding high standardized RIG are chosen for gene-gene interactions. Detection of high-order interactions and comparison of interaction strengths between different orders are made possible by using standardized RIG. We have applied the proposed standardized entropy-based method to two types of data sets from a simulation study and a real genetic association study. We have compared our method and the multifactor dimensionality reduction (MDR) method through power analysis of eight different genetic models with varying penetrance rates, number of SNPs, and sample sizes. Our method shows successful identification of genetic associations and gene-gene interactions both in simulation and real genetic data. Simulation results suggest that the proposed entropy-based method is better able to detect high-order interactions and is superior to the MDR method in most cases. The proposed method is well suited for detecting interactions without main effects as well as for models including main effects.     

Synthetic biology for biofuels: Building designer microbes from the scratch

The ultimate goal in the production of biofuels is to produce fuels identical or similar to petroleum-derived transportation fuels more efficiently and in commercial quantities. Synthetic biologists have been engineering microbes to synthesize biofuels, such as butanol and fatty acid- or isoprenoid-based fuels, which are nearly identical to gasoline and diesel. One of the most urgent demands along this direction is to attain a solid framework for characterizing and standardizing the biological parts and devices. It seems quite promising because biotechnologies specially based on miniaturizations have been making a big contribution to this work. Therefore, in this review, recent advances and difficulties in the biofuel field are discussed, along with the advances of synthetic biology, which will make it possible to create designer microorganisms that produce economically viable next generation biofuels, aside from bioethanol, from corn or sugar cane, and biodiesel from plant or animal oils.     

Quantitative proteomic analysis of ribosomal protein L35b mutant of Saccharomyces cerevisiae

Recent studies have revealed that in higher eukaryotes, several ribosomal proteins are involved in some pathological events or developmental defects, indicating that ribosomal proteins perform unconventional functions other than protein biosynthesis. To obtain an insight into the novel roles of ribosomal proteins, we aimed to analyze the changes in proteome expression in ribosomal protein mutants by using Saccharomyces cerevisiae as a model system. We introduced the rpl35bΔ mutation into the 4159 green fluorescent protein (GFP)-tagged yeast strains by using the synthetic genetic array (SGA) method, and performed quantitative proteomic analysis by using a multilabel microplate reader and flow cytometer. We identified 22 upregulated and 20 downregulated proteins in the rpl35bΔ mutant. These proteins were primarily classified into the Gene Ontology (GO) categories of cellular biosynthetic process, translation, protein or nucleotide metabolic process, cell wall organization and biogenesis, and hyperosmotic response. We also investigated the correlation between the mRNA and protein levels of the identified proteins. Our results show that a ribosomal protein mutation can lead to perturbation in the expression of several proteins, including some other ribosomal proteins. Furthermore, our approach of combining a library of GFP-tagged yeast strains and the SGA method provides an effective and highly sensitive method for dynamic analysis of the effects of various mutations on proteome expression.     

Error-pooling-based statistical methods for identifying novel temporal replication profiles of human chromosomes observed by DNA tiling arrays

Statistical analysis on tiling array data is extremely challenging due to the astronomically large number of sequence probes, high noise levels of individual probes and limited number of replicates in these data. To overcome these difficulties, we first developed statistical error estimation and weighted ANOVA modeling approaches to high-density tiling array data, especially the former based on an advanced error-pooling method to accurately obtain heterogeneous technical error of small-sample tiling array data. Based on these approaches, we analyzed the high-density tiling array data of the temporal replication patterns during cell-cycle S phase of synchronized HeLa cells on human chromosomes 21 and 22. We found many novel temporal replication patterns, identifying about 26% of over 1 million tiling array sequence probes with significant differential replication during the four 2-h time periods of S phase. Among these differentially replicated probes, 126941 sequence probes were matched to 417 known genes. The majority of these genes were found to be replicated within one or two consecutive time periods, while the others were replicated at two non-consecutive time periods. Also, coding regions found to be more differentially replicated in particular time periods than noncoding regions in the gene-poor chromosome 21 (25% differentially replicated among genic probes versus 18.6% among intergenic probes), while such a phenomenon was less prominent in gene-rich chromosome 22. A rigorous statistical testing for local proximity of differentially replicated genic and intergenic probes was performed to identify significant stretches of differentially replicated sequence regions. From this analysis, we found that adjacent genes were frequently replicated at different time periods, potentially implying the existence of quite dense replication origins. Evaluating the conditional probability significance of identified gene ontology terms on chromosomes 21 and 22, we detected some over-represented molecular functions and biological processes among these differentially replicated genes, such as the ones relevant to hydrolase, transferase and receptor-binding activities. Some of these results were confirmed showing >70% consistency with cDNA microarray data that were independently generated in parallel with the tiling arrays. Thus, our improved analysis approaches specifically designed for high-density tiling array data enabled us to reliably and sensitively identify many novel temporal replication patterns on human chromosomes.     

Molecular and ecological analyses of microbial community structures in biofilms of a full-scale Aerated Up-Flow Biobead process

Molecular and cultivation techniques were used to characterize the bacterial communities of biobead reactor biofilms in a sewage treatment plant to which an Aerated Up-Flow Biobead process was applied. With this biobead process, the monthly average values of various chemical parameters in the effluent were generally kept under the regulation limits of the effluent quality of the sewage treatment plant during the operation period. Most probable number (MPN) analysis revealed that the population of denitrifying bacteria was abundant in the biobead #1 reactor, denitrifying and nitrifying bacteria coexisted in the biobead #2 reactor, and nitrifying bacteria prevailed over denitrifying bacteria in the biobead #3 reactor. The results of the MPN test suggested that the biobead #2 reactor was a transition zone leading to acclimated nitrifying biofilms in the biobead #3 reactor. Phylogenetic analysis of 16S rDNA sequences cloned from biofilms showed that the biobead #1 reactor, which received a high organic loading rate, had much diverse microorganisms, whereas the biobead #2 and #3 reactors were dominated by the members of Proteobacteria. DGGE analysis with the ammonia monooxygenase (amoA) gene supported the observation from the MPN test that the biofilms of September were fully developed and specialized for nitrification in the biobead reactor #3. All of the DNA sequences of the amoA DGGE bands were very similar to the sequence of the amoA gene of Nitrosomonas species, the presence of which is typical in the biological aerated filters. The results of this study showed that organic and inorganic nutrients were efficiently removed by both denitrifying microbial populations in the anaerobic tank and heterotrophic and nitrifying bacterial biofilms well-formed in the three functional biobead reactors in the Aerated Up-Flow Biobead process.     

Log-linear model-based multifactor dimensionality reduction method to detect gene gene interactions

MOTIVATION: The identification and characterization of susceptibility genes that influence the risk of common and complex diseases remains a statistical and computational challenge in genetic association studies. This is partly because the effect of any single genetic variant for a common and complex disease may be dependent on other genetic variants (gene-gene interaction) and environmental factors (gene-environment interaction). To address this problem, the multifactor dimensionality reduction (MDR) method has been proposed by Ritchie et al. to detect gene-gene interactions or gene-environment interactions. The MDR method identifies polymorphism combinations associated with the common and complex multifactorial diseases by collapsing high-dimensional genetic factors into a single dimension. That is, the MDR method classifies the combination of multilocus genotypes into high-risk and low-risk groups based on a comparison of the ratios of the numbers of cases and controls. When a high-order interaction model is considered with multi-dimensional factors, however, there may be many sparse or empty cells in the contingency tables. The MDR method cannot classify an empty cell as high risk or low risk and leaves it as undetermined. RESULTS: In this article, we propose the log-linear model-based multifactor dimensionality reduction (LM MDR) method to improve the MDR in classifying sparse or empty cells. The LM MDR method estimates frequencies for empty cells from a parsimonious log-linear model so that they can be assigned to high-and low-risk groups. In addition, LM MDR includes MDR as a special case when the saturated log-linear model is fitted. Simulation studies show that the LM MDR method has greater power and smaller error rates than the MDR method. The LM MDR method is also compared with the MDR method using as an example sporadic Alzheimer's disease.     

Odds ratio based multifactor-dimensionality reduction method for detecting gene-gene interactions

MOTIVATION: The identification and characterization of genes that increase the susceptibility to common complex multifactorial diseases is a challenging task in genetic association studies. The multifactor dimensionality reduction (MDR) method has been proposed and implemented by Ritchie et al. (2001) to identify the combinations of multilocus genotypes and discrete environmental factors that are associated with a particular disease. However, the original MDR method classifies the combination of multilocus genotypes into high-risk and low-risk groups in an ad hoc manner based on a simple comparison of the ratios of the number of cases and controls. Hence, the MDR approach is prone to false positive and negative errors when the ratio of the number of cases and controls in a combination of genotypes is similar to that in the entire data, or when both the number of cases and controls is small. Hence, we propose the odds ratio based multifactor dimensionality reduction (OR MDR) method that uses the odds ratio as a new quantitative measure of disease risk. RESULTS: While the original MDR method provides a simple binary measure of risk, the OR MDR method provides not only the odds ratio as a quantitative measure of risk but also the ordering of the multilocus combinations from the highest risk to lowest risk groups. Furthermore, the OR MDR method provides a confidence interval for the odds ratio for each multilocus combination, which is extremely informative in judging its importance as a risk factor. The proposed OR MDR method is illustrated using the dataset obtained from the CDC Chronic Fatigue Syndrome Research Group. AVAILABILITY: The program written in R is available.     

Combining multiple microarrays in the presence of controlling variables

MOTIVATION: Microarray technology enables the monitoring of expression levels for thousands of genes simultaneously. When the magnitude of the experiment increases, it becomes common to use the same type of microarrays from different laboratories or hospitals. Thus, it is important to analyze microarray data together to derive a combined conclusion after accounting for the differences. One of the main objectives of the microarray experiment is to identify differentially expressed genes among the different experimental groups. The analysis of variance (ANOVA) model has been commonly used to detect differentially expressed genes after accounting for the sources of variation commonly observed in the microarray experiment. RESULTS: We extended the usual ANOVA model to account for an additional variability resulting from many confounding variables such as the effect of different hospitals. The proposed model is a two-stage ANOVA model. The first stage is the adjustment for the effects of no interests. The second stage is the detection of differentially expressed genes among the experimental groups using the residuals obtained from the first stage. Based on these residuals, we propose a permutation test to detect the differentially expressed genes. The proposed model is illustrated using the data from 133 microarrays collected at three different hospitals. The proposed approach is more flexible to use, and it is easier to accommodate the individual covariates in this model than using the meta-analysis approach. AVAILABILITY: A set of programs written in R will be electronically sent upon request.