A General Model for Multilocus Epistatic Interactions in Case-Control Studies

Background

Epistasis, i.e., the interaction of alleles at different loci, is thought to play a central role in the formation and progression of complex diseases. The complexity of disease expression should arise from a complex network of epistatic interactions involving multiple genes.

Methodology

We develop a general model for testing high-order epistatic interactions for a complex disease in a case-control study. We incorporate the quantitative genetic theory of high-order epistasis into the setting of cases and controls sampled from a natural population. The new model allows the identification and testing of epistasis and its various genetic components.

Conclusions

Simulation studies were used to examine the power and false positive rates of the model under different sampling strategies. The model was used to detect epistasis in a case-control study of inflammatory bowel disease, in which five SNPs at a candidate gene were typed, leading to the identification of a significant three-locus epistasis.     

OpenDMAP: An open source, ontology-driven concept analysis engine, with applications to capturing knowledge regarding protein transport, protein interactions and cell-type-specific gene expression

Background

Microarray technology provides an efficient means for globally exploring physiological processes governed by the coordinated expression of multiple genes. However, identification of genes differentially expressed in microarray experiments is challenging because of their potentially high type I error rate. Methods for large-scale statistical analyses have been developed but most of them are applicable to two-sample or two-condition data.

Results

We developed a large-scale multiple-group F-test based method, named ranking analysis of F-statistics (RAF), which is an extension of ranking analysis of microarray data (RAM) for two-sample t-test. In this method, we proposed a novel random splitting approach to generate the null distribution instead of using permutation, which may not be appropriate for microarray data. We also implemented a two-simulation strategy to estimate the false discovery rate. Simulation results suggested that it has higher efficiency in finding differentially expressed genes among multiple classes at a lower false discovery rate than some commonly used methods. By applying our method to the experimental data, we found 107 genes having significantly differential expressions among 4 treatments at <0.7% FDR, of which 31 belong to the expressed sequence tags (ESTs), 76 are unique genes who have known functions in the brain or central nervous system and belong to six major functional groups.

Conclusion

Our method is suitable to identify differentially expressed genes among multiple groups, in particular, when sample size is small.     

A complete classification of epistatic two-locus models

Background

The study of epistasis is of great importance in statistical genetics in fields such as linkage and association analysis and QTL mapping. In an effort to classify the types of epistasis in the case of two biallelic loci Li and Reich listed and described all models in the simplest case of 0/1 penetrance values. However, they left open the problem of finding a classification of two-locus models with continuous penetrance values.

Results

We provide a complete classification of biallelic two-locus models. In addition to solving the classification problem for dichotomous trait disease models, our results apply to any instance where real numbers are assigned to genotypes, and provide a complete framework for studying epistasis in QTL data. Our approach is geometric and we show that there are 387 distinct types of two-locus models, which can be reduced to 69 when symmetry between loci and alleles is accounted for. The model types are defined by 86 circuits, which are linear combinations of genotype values, each of which measures a fundamental unit of interaction.

Conclusion

The circuits provide information on epistasis beyond that contained in the additive × additive, additive × dominance, and dominance × dominance interaction terms. We discuss the connection between our classification and standard epistatic models and demonstrate its utility by analyzing a previously published dataset.     

Investigation of non-uniform airflow signal oscillation during high frequency chest compression

Background  

High frequency chest compression (HFCC) is a useful and popular therapy for clearing bronchial airways of excessive or thicker mucus. Our observation of respiratory airflow of a subject during use of HFCC showed the airflow oscillation by HFCC was strongly influenced by the nonlinearity of the respiratory system. We used a computational model-based approach to analyse the respiratory airflow during use of HFCC.     

Natural selection governs local, but not global, evolutionary gene coexpression networks in Caenorhabditis elegans

Background

Signal duration (e.g. the time over which an active signaling intermediate persists) is a key regulator of biological decisions in myriad contexts such as cell growth, proliferation, and developmental lineage commitments. Accompanying differences in signal duration are numerous downstream biological processes that require multiple steps of biochemical regulation.

Results

Here we present an analysis that investigates how simple biochemical motifs that involve multiple stages of regulation can be constructed to differentially process signals that persist at different time scales. We compute the dynamic, frequency dependent gain within these networks and resulting power spectra to better understand how biochemical networks can integrate signals at different time scales. We identify topological features of these networks that allow for different frequency dependent signal processing properties.

Conclusion

We show that multi-staged cascades are effective in integrating signals of long duration whereas multi-staged cascades that operate in the presence of negative feedback are effective in integrating signals of short duration. Our studies suggest principles for why signal duration in connection with multiple steps of downstream regulation is a ubiquitous motif in biochemical systems.     

Expression, purification and structural analysis of the Pyrococcus abyssi RNA binding protein PAB1135

Background

Epistatic interactions of multiple single nucleotide polymorphisms (SNPs) are now believed to affect individual susceptibility to common diseases. The detection of such interactions, however, is a challenging task in large scale association studies. Ant colony optimization (ACO) algorithms have been shown to be useful in detecting epistatic interactions.

Findings

AntEpiSeeker, a new two-stage ant colony optimization algorithm, has been developed for detecting epistasis in a case-control design. Based on some practical epistatic models, AntEpiSeeker has performed very well.

Conclusions

AntEpiSeeker is a powerful and efficient tool for large-scale association studies and can be downloaded from http://nce.ads.uga.edu/~romdhane/AntEpiSeeker/index.html.     

QTLRel: an R Package for Genome-wide Association Studies in which Relatedness is a Concern

Background

Existing software for quantitative trait mapping is either not able to model polygenic variation or does not allow incorporation of more than one genetic variance component. Improperly modeling the genetic relatedness among subjects can result in excessive false positives. We have developed an R package, QTLRel, to enable more flexible modeling of genetic relatedness as well as covariates and non-genetic variance components.

Results

We have successfully used the package to analyze many datasets, including F₃₄ body weight data that contains 688 individuals genotyped at 3105 SNP markers and identified 11 QTL. It took 295 seconds to estimate variance components and 70 seconds to perform the genome scan on an Linux machine equipped with a 2.40GHz Intel(R) Core(TM)2 Quad CPU.

Conclusions

QTLRel provides a toolkit for genome-wide association studies that is capable of calculating genetic incidence matrices from pedigrees, estimating variance components, performing genome scans, incorporating interactive covariates and genetic and non-genetic variance components, as well as other functionalities such as multiple-QTL mapping and genome-wide epistasis.     

BglBrick vectors and datasheets: A synthetic biology platform for gene expression

Background

As engineered biological systems become more complex, it is increasingly common to express multiple operons from different plasmids and inducible expression systems within a single host cell. Optimizing such systems often requires screening combinations of origins of replication, expression systems, and antibiotic markers. This procedure is hampered by a lack of quantitative data on how these components behave when more than one origin of replication or expression system are used simultaneously. Additionally, this process can be time consuming as it often requires the creation of new vectors or cloning into existing but disparate vectors.

Results

Here, we report the development and characterization of a library of expression vectors compatible with the BglBrick standard (BBF RFC 21). We have designed and constructed 96 BglBrick-compatible plasmids with a combination of replication origins, antibiotic resistance genes, and inducible promoters. These plasmids were characterized over a range of inducer concentrations, in the presence of non-cognate inducer molecules, and with several growth media, and their characteristics were documented in a standard format datasheet. A three plasmid system was used to investigate the impact of multiple origins of replication on plasmid copy number.

Conclusions

The standardized collection of vectors presented here allows the user to rapidly construct and test the expression of genes with various combinations of promoter strength, inducible expression system, copy number, and antibiotic resistance. The quantitative datasheets created for these vectors will increase the predictability of gene expression, especially when multiple plasmids and inducers are utilized.     

Assessment of algorithms for high throughput detection of genomic copy number variation in oligonucleotide microarray data

Background

Genomic deletions and duplications are important in the pathogenesis of diseases, such as cancer and mental retardation, and have recently been shown to occur frequently in unaffected individuals as polymorphisms. Affymetrix GeneChip whole genome sampling analysis (WGSA) combined with 100 K single nucleotide polymorphism (SNP) genotyping arrays is one of several microarray-based approaches that are now being used to detect such structural genomic changes. The popularity of this technology and its associated open source data format have resulted in the development of an increasing number of software packages for the analysis of copy number changes using these SNP arrays.

Results

We evaluated four publicly available software packages for high throughput copy number analysis using synthetic and empirical 100 K SNP array data sets, the latter obtained from 107 mental retardation (MR) patients and their unaffected parents and siblings. We evaluated the software with regards to overall suitability for high-throughput 100 K SNP array data analysis, as well as effectiveness of normalization, scaling with various reference sets and feature extraction, as well as true and false positive rates of genomic copy number variant (CNV) detection.

Conclusion

We observed considerable variation among the numbers and types of candidate CNVs detected by different analysis approaches, and found that multiple programs were needed to find all real aberrations in our test set. The frequency of false positive deletions was substantial, but could be greatly reduced by using the SNP genotype information to confirm loss of heterozygosity.     

Signals of demographic expansion in Drosophila virilis

Background

The antagonistic co-evolution of hosts and their parasites is considered to be a potential driving force in maintaining host genetic variation including sexual reproduction and recombination. The examination of this hypothesis calls for information about the genetic basis of host-parasite interactions – such as how many genes are involved, how big an effect these genes have and whether there is epistasis between loci. We here examine the genetic architecture of quantitative resistance in animal and plant hosts by concatenating published studies that have identified quantitative trait loci (QTL) for host resistance in animals and plants.

Results

Collectively, these studies show that host resistance is affected by few loci. We particularly show that additional epistatic interactions, especially between loci on different chromosomes, explain a majority of the effects. Furthermore, we find that when experiments are repeated using different host or parasite genotypes under otherwise identical conditions, the underlying genetic architecture of host resistance can vary dramatically – that is, involves different QTLs and epistatic interactions. QTLs and epistatic loci vary much less when host and parasite types remain the same but experiments are repeated in different environments.

Conclusion

This pattern of variability of the genetic architecture is predicted by strong interactions between genotypes and corroborates the prevalence of varying host-parasite combinations over varying environmental conditions. Moreover, epistasis is a major determinant of phenotypic variance for host resistance. Because epistasis seems to occur predominantly between, rather than within, chromosomes, segregation and chromosome number rather than recombination via cross-over should be the major elements affecting adaptive change in host resistance.     

Multiple sclerosis and pregnancy: what does the patient think? a questionnaire study

Background

The first objective of a DNA microarray experiment is typically to generate a list of genes or probes that are found to be differentially expressed or represented (in the case of comparative genomic hybridizations and/or copy number variation) between two conditions or strains. Rank Products analysis comprises a robust algorithm for deriving such lists from microarray experiments that comprise small numbers of replicates, for example, less than the number required for the commonly used t-test. Currently, users wishing to apply Rank Products analysis to their own microarray data sets have been restricted to the use of command line-based software which can limit its usage within the biological community.

Findings

Here we have developed a web interface to existing Rank Products analysis tools allowing users to quickly process their data in an intuitive and step-wise manner to obtain the respective Rank Product or Rank Sum, probability of false prediction and p-values in a downloadable file.

Conclusions

The online interactive Rank Products analysis tool RankProdIt, for analysis of any data set containing measurements for multiple replicated conditions, is available at: http://strep-microarray.sbs.surrey.ac.uk/RankProducts     

Characterization of antibodies elicited by XMRV infection and development of immunoassays useful for epidemiologic studies

Background

An essential event during the replication cycle of HIV-1 is the integration of the reverse transcribed viral DNA into the host cellular genome. Our former report revealed that HIV-1 integrase (IN), the enzyme that catalyzes the integration reaction, is positively regulated by acetylation mediated by the histone acetyltransferase (HAT) p300.

Results

In this study we demonstrate that another cellular HAT, GCN5, acetylates IN leading to enhanced 3'-end processing and strand transfer activities. GCN5 participates in the integration step of HIV-1 replication cycle as demonstrated by the reduced infectivity, due to inefficient provirus formation, in GCN5 knockdown cells. Within the C-terminal domain of IN, four lysines (K258, K264, K266, and K273) are targeted by GCN5 acetylation, three of which (K264, K266, and K273) are also modified by p300. Replication analysis of HIV-1 clones carrying substitutions at the IN lysines acetylated by both GCN5 and p300, or exclusively by GCN5, demonstrated that these residues are required for efficient viral integration. In addition, a comparative analysis of the replication efficiencies of the IN triple- and quadruple-mutant viruses revealed that even though the lysines targeted by both GCN5 and p300 are required for efficient virus integration, the residue exclusively modified by GCN5 (K258) does not affect this process.

Conclusions

The results presented here further demonstrate the relevance of IN post-translational modification by acetylation, which results from the catalytic activities of multiple HATs during the viral replication cycle. Finally, this study contributes to clarifying the recent debate raised on the role of IN acetylated lysines during HIV-1 infection.     

Predicting the sensitivity and specificity of published real-time PCR assays

Background

In recent years real-time PCR has become a leading technique for nucleic acid detection and quantification. These assays have the potential to greatly enhance efficiency in the clinical laboratory. Choice of primer and probe sequences is critical for accurate diagnosis in the clinic, yet current primer/probe signature design strategies are limited, and signature evaluation methods are lacking.

Methods

We assessed the quality of a signature by predicting the number of true positive, false positive and false negative hits against all available public sequence data. We found real-time PCR signatures described in recent literature and used a BLAST search based approach to collect all hits to the primer-probe combinations that should be amplified by real-time PCR chemistry. We then compared our hits with the sequences in the NCBI taxonomy tree that the signature was designed to detect.

Results

We found that many published signatures have high specificity (almost no false positives) but low sensitivity (high false negative rate). Where high sensitivity is needed, we offer a revised methodology for signature design which may designate that multiple signatures are required to detect all sequenced strains. We use this methodology to produce new signatures that are predicted to have higher sensitivity and specificity.

Conclusion

We show that current methods for real-time PCR assay design have unacceptably low sensitivities for most clinical applications. Additionally, as new sequence data becomes available, old assays must be reassessed and redesigned. A standard protocol for both generating and assessing the quality of these assays is therefore of great value. Real-time PCR has the capacity to greatly improve clinical diagnostics. The improved assay design and evaluation methods presented herein will expedite adoption of this technique in the clinical lab.     

Gelatinase (MMP-2 and -9) expression profiles during gestation in the bovine endometrium

Background

The angiogenic and invasive properties of the cytotrophoblast are crucial to provide an adequate area for feto-maternal exchange. The present study aimed at identifying the localization of interrelated angiogenic, hyperpermeability and vasodilator factors in the feto-maternal interface in pregnant guinea-pigs.

Methods

Utero-placental units were collected from early to term pregnancy. VEGF, Flt-1, KDR, B2R and eNOS were analyzed by immunohistochemistry, and the intensity of the signals in placenta and syncytial streamers was digitally analysed. Flt1 and eNOS content of placental homogenates was determined by western blotting. Statistical analysis used one-way analysis of variance and Tukey's Multiple Comparison post-hoc test.

Results

In the subplacenta, placental interlobium and labyrinth VEGF, Flt-1, KDR, B2R and eNOS were expressed in all stages of pregnancy. Syncytial streamers in all stages of gestation, and cytotrophoblasts surrounding myometrial arteries in early and mid pregnancy – and replacing the smooth muscle at term – displayed immunoreactivity for VEGF, Flt-1, KDR, eNOS and B2R. In partly disrupted mesometrial arteries in late pregnancy cytotrophoblasts and endothelial cells expressed VEGF, Flt-1, KDR, B2R and eNOS. Sections incubated in absence of the first antibody, or in presence of rabbit IgG fraction and mouse IgG serum, yielded no staining. According to the digital analysis, Flt-1 increased in the placental interlobium in days 40 and 60 as compared to day 20 (P = 0.016), and in the labyrinth in day 60 as compared to days 20 and 40 (P = 0.026), while the signals for VEGF, KDR, B2R, and eNOS showed no variations along pregnancy. In syncytial streamers the intensity of VEGF immunoreactivity was increased in day 40 in comparison to day 20 (P = 0.027), while that of B2R decreased in days 40 and 60 as compared to day 20 (P = 0.011); VEGF, Flt-1, KDR, B2R and eNOS expression showed no variations. Western blots for eNOS and Flt-1 in placental homogenates showed no significant temporal differences along pregnancy.

Conclusion

The demonstration of different angiogenic, hyperpermeability and vasodilator factors in the same cellular protagonists of angiogenesis and invasion in the pregnant guinea-pig, supports the presence of a functional network, and strengthens the argument that this species provides an adequate model to understand human pregnancy.     

Midkine is a NF-κB-inducible gene that supports prostate cancer cell survival

Background

Atherosclerotic peripheral arterial disease (PAD) affects 8–10 million people in the United States and is associated with a marked impairment in quality of life and an increased risk of cardiovascular events. Noninvasive assessment of PAD is performed by measuring the ankle-brachial index (ABI). Complex traits, such as ABI, are influenced by a large array of genetic and environmental factors and their interactions. We attempted to characterize the genetic architecture of ABI by examining the main and interactive effects of individual single nucleotide polymorphisms (SNPs) and conventional risk factors.

Methods

We applied linear regression analysis to investigate the association of 435 SNPs in 112 positional and biological candidate genes with ABI and related physiological and biochemical traits in 1046 non-Hispanic white, hypertensive participants from the Genetic Epidemiology Network of Arteriopathy (GENOA) study. The main effects of each SNP, as well as SNP-covariate and SNP-SNP interactions, were assessed to investigate how they contribute to the inter-individual variation in ABI. Multivariable linear regression models were then used to assess the joint contributions of the top SNP associations and interactions to ABI after adjustment for covariates. We reduced the chance of false positives by 1) correcting for multiple testing using the false discovery rate, 2) internal replication, and 3) four-fold cross-validation.

Results

When the results from these three procedures were combined, only two SNP main effects in NOS3, three SNP-covariate interactions (ADRB2 Gly 16 – lipoprotein(a) and SLC4A5 – diabetes interactions), and 25 SNP-SNP interactions (involving SNPs from 29 different genes) were significant, replicated, and cross-validated. Combining the top SNPs, risk factors, and their interactions into a model explained nearly 18% of variation in ABI in the sample. SNPs in six genes (ADD2, ATP6V1B1, PRKAR2B, SLC17A2, SLC22A3, and TGFB3) were also influencing triglycerides, C-reactive protein, homocysteine, and lipoprotein(a) levels.

Conclusion

We found that candidate gene SNP main effects, SNP-covariate and SNP-SNP interactions contribute to the inter-individual variation in ABI, a marker of PAD. Our findings underscore the importance of conducting systematic investigations that consider context-dependent frameworks for developing a deeper understanding of the multidimensional genetic and environmental factors that contribute to complex diseases.     

Interaction among apoptosis-associated sequence variants and joint effects on aggressive prostate cancer

Background

Molecular and epidemiological evidence demonstrate that altered gene expression and single nucleotide polymorphisms in the apoptotic pathway are linked to many cancers. Yet, few studies emphasize the interaction of variant apoptotic genes and their joint modifying effects on prostate cancer (PCA) outcomes. An exhaustive assessment of all the possible two-, three- and four-way gene-gene interactions is computationally burdensome. This statistical conundrum stems from the prohibitive amount of data needed to account for multiple hypothesis testing.

Methods

To address this issue, we systematically prioritized and evaluated individual effects and complex interactions among 172 apoptotic SNPs in relation to PCA risk and aggressive disease (i.e., Gleason score ≥ 7 and tumor stages III/IV). Single and joint modifying effects on PCA outcomes among European-American men were analyzed using statistical epistasis networks coupled with multi-factor dimensionality reduction (SEN-guided MDR). The case-control study design included 1,175 incident PCA cases and 1,111 controls from the prostate, lung, colo-rectal, and ovarian (PLCO) cancer screening trial. Moreover, a subset analysis of PCA cases consisted of 688 aggressive and 488 non-aggressive PCA cases. SNP profiles were obtained using the NCI Cancer Genetic Markers of Susceptibility (CGEMS) data portal. Main effects were assessed using logistic regression (LR) models. Prior to modeling interactions, SEN was used to pre-process our genetic data. SEN used network science to reduce our analysis from > 36 million to < 13,000 SNP interactions. Interactions were visualized, evaluated, and validated using entropy-based MDR. All parametric and non-parametric models were adjusted for age, family history of PCA, and multiple hypothesis testing.

Results

Following LR modeling, eleven and thirteen sequence variants were associated with PCA risk and aggressive disease, respectively. However, none of these markers remained significant after we adjusted for multiple comparisons. Nevertheless, we detected a modest synergistic interaction between AKT3 rs2125230-PRKCQ rs571715 and disease aggressiveness using SEN-guided MDR (p = 0.011).

Conclusions

In summary, entropy-based SEN-guided MDR facilitated the logical prioritization and evaluation of apoptotic SNPs in relation to aggressive PCA. The suggestive interaction between AKT3-PRKCQ and aggressive PCA requires further validation using independent observational studies.