General Epistatic Models of the Risk of Complex Diseases

The range of possible gene interactions in a multilocus model of a complex inherited disease is studied by exploring genotype-specific risks subject to the constraint that the allele frequencies and marginal risks are known. We quantify the effect of gene interactions by defining the interaction ratio, , where K_R is the recurrence risk to relatives with relationship R for the true model and is the recurrence risk to relatives for a multiplicative model with the same marginal risks. We use a Markov chain Monte Carlo (MCMC) procedure to sample from the space of possible models. We find that the average of C_R increases with the number of loci for both low frequency (p = 0.03) and higher frequency (p = 0.25) causative alleles. Furthermore, the probability that C_R > 1 is nearly 1. Similar results are obtained when more weight is given to risk models that are closer to the comparable multiplicative model. These results imply that, in general, gene interactions will result in greater heritability of a complex inherited disease than is expected on the basis of a multiplicative model of interactions and hence may provide a partial explanation for the problem of missing heritability of complex diseases.ALTHOUGH many genome-wide association studies (GWAS) have been performed and have found hundreds of SNPs associated with higher risk of complex inherited diseases, those SNPs so far account for only a small fraction of the inherited risk of those diseases (Altshuler et al. 2008). Several not mutually exclusive explanations have been proposed for the “missing heritability,” i.e., the heritability that is not yet accounted for by SNPs found in GWAS (Manolio et al. 2009): (i) common alleles of small effect that have not been found because GWAS done so far have been underpowered, (ii) low-frequency alleles of moderate effect that are difficult to find using HapMap SNPs, (iii) rare copy-number variants that are not in strong linkage disequilibrium (LD) with HapMap SNPs, (iv) inherited epigenetic factors that are not in strong LD with HapMap SNPs, and (v) interactions among causative alleles that conceal their true contribution to heritability. In this article we investigate the last possibility and determine the extent to which interactions may account for missing heritability.Our analysis is in the same spirit as that of Culverhouse et al. (2002). We assume that the risk of being affected by a complex disease is determined by an individual''s genotype at two or more loci and that the frequencies of causative alleles and the average risks for each one-locus genotype (the marginal risks) are known. Culverhouse et al. (2002) assumed the marginal risks were the same for all genotypes and all loci. In that case, causative alleles have odds ratios of 1; they contribute to risk only through their interactions. Culverhouse et al. found the risk function that maximized the heritability and showed that the maximum possible heritability attributable to interactions increased with the number of loci. They concluded that it is quite possible that interactions among loci that have no main effect could contribute substantially to the heritability of a complex disease and indeed could account for “virtually all the variation in affection status for diseases with any prevalence” (Culverhouse et al. 2002, p. 468).We generalize the analysis of Culverhouse et al. in three ways. First, we allow causative alleles to have odds ratios >1. Second, we explore the entire space of models instead of focusing only on the risk model that maximizes heritability. Third, we examine how the importance of gene interactions depends on the “distance” between a risk model and a comparable multiplicative model. We show that gene interactions can substantially increase the heritability of risk as measured by recurrence risk, K_R, and that the effect increases with the number of loci carrying causative alleles. Furthermore, we show that these results are true even if more weight is given to models that are closer to a comparable multiplicative model.Geometrically, the space of feasible genotype-specific risks subject to the aforementioned constraints (i.e., that the allele frequencies and marginal risks are known) corresponds to a high-dimensional convex polytope, and the computational problem of interest involves integrating a quadratic function over the polytope. The dimension of the polytope grows exponentially with the number of loci, and, therefore, analytic computation is intractable for more than two loci. Hence, we devise a Monte Carlo approach to tackle the problem. Note that, because of high dimensionality, rejection algorithms are not appropriate for this kind of problem. We instead employ a Markov chain Monte Carlo (MCMC) algorithm based on a random walk that always stays inside the polytope. We present empirical results for up to five loci and obtain a closed-form formula for the minimum of K_R over the polytope; the latter result applies to an arbitrary number of loci. Interestingly, the minimum of K_R decreases as the number L of loci increases, but the average of K_R over the polytope increases with L.     

eNOS-β-Actin Interaction Contributes to Increased Peroxynitrite Formation during Hyperoxia in Pulmonary Artery Endothelial Cells and Mouse Lungs

Oxygen toxicity is the most severe side effect of oxygen therapy in neonates and adults. Pulmonary damage of oxygen toxicity is related to the overproduction of reactive oxygen species (ROS). In the present study, we investigated the effect of hyperoxia on the production of peroxynitrite in pulmonary artery endothelial cells (PAEC) and mouse lungs. Incubation of PAEC under hyperoxia (95% O₂) for 24 h resulted in an increase in peroxynitrite formation. Uric acid, a peroxynitrite scavenger, prevented hyperoxia-induced increase in peroxynitrite. The increase in peroxynitrite formation is accompanied by increases in nitric oxide (NO) release and endothelial NO synthase (eNOS) activity. We have previously reported that association of eNOS with β-actin increases eNOS activity and NO production in lung endothelial cells. To study whether eNOS-β-actin association contributes to increased peroxynitrite production, eNOS-β-actin interaction were inhibited by reducing β-actin availability or by using a synthetic peptide (P326TAT) containing a sequence corresponding to the actin binding site on eNOS. We found that disruption of eNOS-β-actin interaction prevented hyperoxia-induced increases in eNOS-β-actin association, eNOS activity, NO and peroxynitrite production, and protein tyrosine nitration. Hyperoxia failed to induce the increases in eNOS activity, NO and peroxynitrite formation in COS-7 cells transfected with plasmids containing eNOS mutant cDNA in which amino acids leucine and tryptophan were replaced with alanine in the actin binding site on eNOS. Exposure of mice to hyperoxia resulted in significant increases in eNOS-β-actin association, eNOS activity, and protein tyrosine nitration in the lungs. Our data indicate that increased association of eNOS with β-actin in PAEC contributes to hyperoxia-induced increase in the production of peroxynitrite which may cause nitrosative stress in pulmonary vasculature.     

Managing and mining protein crystallization data

The crystallization of macromolecules remains a major bottleneck in structural biology. The routine screening of more than one thousand crystallization conditions and subsequent optimization by fine screening presents a challenge to conventional laboratory notebook keeping. In addition, the development of high-throughput robotic crystallization and imaging systems presents a pressing need for low-cost laboratory information management system (LIMS). Here we describe CLIMS2, a crystallization LIMS that features a simple, user-friendly graphical interface, allowing the storage, management, retrieval and mining of crystallization data. The CLIMS2 executable and documentation is freely available at http://clims.med.monash.edu.au.     

Molecular genotype identification of the Gallus gallus major histocompatibility complex

The chicken major histocompatibility complex (MHC) is commonly defined by serologic reactions of erythrocytes with antibodies specific to the highly polymorphic MHC class I (BF) and MHC class IV (BG) antigens. The microsatellite marker LEI0258 is known to be physically located within the MHC, between the BG and BF regions. DNA from various serologically defined MHC haplotypes was amplified by polymerase chain reaction with primers surrounding this marker. Twenty-six distinctive allele sizes were identified. Some serologically well-defined MHC haplotypes shared a common LEI0258 allele size but could be distinguished either by the addition of information from another nearby marker (MCW0371) or by small indels or single nucleotide polymorphism (SNP) differences between the alleles. The association between LEI0258 allele and serologically defined MHC haplotype was very consistent for the same haplotype from multiple sources. Sequence information for the region defined by LEI0258 was obtained for 51 different haplotypes. Two internal repeats whose lengths were 13 and 12 bp, respectively, are the primary basis for allelic variability. Allele size variation ranges from 182 to 552 bp. Four indels and five SNPs in the surrounding sequence provide additional means for distinguishing alleles. Typing with LEI0258 and MCW0371 will be useful in identifying MHC haplotypes in outbred populations of chickens particularly for the initial development of serological reagents.     

The distribution and evolutionary history of the PRP8 intein

Background  

We recently described a mini-intein in the PRP8 gene of a strain of the basidiomycete Cryptococcus neoformans, an important fungal pathogen of humans. This was the second described intein in the nuclear genome of any eukaryote; the first nuclear encoded intein was found in the VMA gene of several saccharomycete yeasts. The evolution of eukaryote inteins is not well understood. In this report we describe additional PRP8 inteins (bringing the total of these to over 20). We compare and contrast the phylogenetic distribution and evolutionary history of the PRP8 intein and the saccharomycete VMA intein, in order to derive a broader understanding of eukaryote intein evolution. It has been suggested that eukaryote inteins undergo horizontal transfer and the present analysis explores this proposal.     

Nitric oxide reduces NADPH oxidase 5 (Nox5) activity by reversible S-nitrosylation

The NADPH oxidases (Noxs) are a family of transmembrane oxidoreductases that produce superoxide and other reactive oxygen species (ROS). Nox5 was the last of the conventional Nox isoforms to be identified and is a calcium-dependent enzyme that does not depend on accessory subunits for activation. Recently, Nox5 was shown to be expressed in human blood vessels and therefore the goal of this study was to determine whether nitric oxide (NO) can modulate Nox5 activity. Endogenously produced NO potently inhibited basal and stimulated Nox5 activity and this inhibition was reversible with chronic, but not acute, exposure to L-NAME. Nox5 activity was reduced by NO donors, iNOS, and eNOS and in endothelial cells and LPS-stimulated smooth muscle cells in a manner dependent on NO concentration. ROS production was diminished by NO in an isolated enzyme activity assay replete with surplus calcium and NADPH. There was no evidence for NO-dependent changes in tyrosine nitration, glutathiolation, or phosphorylation of Nox5. In contrast, there was evidence for the increased nitrosylation of Nox5 as determined by the biotin-switch assay and mass spectrometry. Four S-nitrosylation sites were identified and of these, mutation of C694 dramatically lowered Nox5 activity, NO sensitivity, and biotin labeling. Furthermore, coexpression of the denitrosylation enzymes thioredoxin 1 and GSNO reductase prevented NO-dependent inhibition of Nox5. The potency of NO against other Nox enzymes was in the order Nox1 ≥ Nox3 > Nox5 > Nox2, whereas Nox4 was refractory. Collectively, these results suggest that endogenously produced NO can directly S-nitrosylate and inhibit the activity of Nox5.     

Interleukin-10 protects cultured fetal rat type II epithelial cells from injury induced by mechanical stretch

Mechanical ventilation plays a central role in the pathogenesis of bronchopulmonary dysplasia. However, the mechanisms by which excessive stretch of fetal or neonatal type II epithelial cells contributes to lung injury are not well defined. In these investigations, isolated embryonic day 19 fetal rat type II epithelial cells were cultured on substrates coated with fibronectin and exposed to 5% or 20% cyclic stretch to simulate mechanical forces during lung development or lung injury, respectively. Twenty percent stretch of fetal type II epithelial cells increased necrosis, apoptosis, and proliferation compared with control, unstretched samples. By ELISA and real-time PCR (qRT-PCR), 20% stretch increased secretion of IL-8 into the media and IL-8 gene expression and inhibited IL-10 release. Interestingly, administration of recombinant IL-10 before 20% stretch did not affect cell lysis but significantly reduced apoptosis and IL-8 release compared with stretched samples without IL-10. Collectively, our studies suggest that IL-10 may play an important role in protection of fetal type II epithelial cells from injury secondary to stretch.     

Diet-induced obesity is linked to marked but reversible alterations in the mouse distal gut microbiome

We have investigated the interrelationship between diet, gut microbial ecology, and energy balance using a mouse model of obesity produced by consumption of a prototypic Western diet. Diet-induced obesity (DIO) produced a bloom in a single uncultured clade within the Mollicutes class of the Firmicutes, which was diminished by subsequent dietary manipulations that limit weight gain. Microbiota transplantation from mice with DIO to lean germ-free recipients promoted greater fat deposition than transplants from lean donors. Metagenomic and biochemical analysis of the gut microbiome together with sequencing and metabolic reconstructions of a related human gut-associated Mollicute (Eubacterium dolichum) revealed features that may provide a competitive advantage to members of the bloom in the Western diet nutrient milieu, including import and processing of simple sugars. Our study illustrates how combining comparative metagenomics with gnotobiotic mouse models and specific dietary manipulations can disclose the niches of previously uncharacterized members of the gut microbiota.