Using Linkage Analysis to Detect Gene-Gene Interactions. 2. Improved Reliability and Extension to More-Complex Models

Detecting gene-gene interaction in complex diseases has become an important priority for common disease genetics, but most current approaches to detecting interaction start with disease-marker associations. These approaches are based on population allele frequency correlations, not genetic inheritance, and therefore cannot exploit the rich information about inheritance contained within families. They are also hampered by issues of rigorous phenotype definition, multiple test correction, and allelic and locus heterogeneity. We recently developed, tested, and published a powerful gene-gene interaction detection strategy based on conditioning family data on a known disease-causing allele or a disease-associated marker allele⁴. We successfully applied the method to disease data and used computer simulation to exhaustively test the method for some epistatic models. We knew that the statistic we developed to indicate interaction was less reliable when applied to more-complex interaction models. Here, we improve the statistic and expand the testing procedure. We computer-simulated multipoint linkage data for a disease caused by two interacting loci. We examined epistatic as well as additive models and compared them with heterogeneity models. In all our models, the at-risk genotypes are “major” in the sense that among affected individuals, a substantial proportion has a disease-related genotype. One of the loci (A) has a known disease-related allele (as would have been determined from a previous analysis). We removed (pruned) family members who did not carry this allele; the resultant dataset is referred to as “stratified.” This elimination step has the effect of raising the “penetrance” and detectability at the second locus (B). We used the lod scores for the stratified and unstratified data sets to calculate a statistic that either indicated the presence of interaction or indicated that no interaction was detectable. We show that the new method is robust and reliable for a wide range of parameters. Our statistic performs well both with the epistatic models (false negative rates, i.e., failing to detect interaction, ranging from 0 to 2.5%) and with the heterogeneity models (false positive rates, i.e., falsely detecting interaction, ≤1%). It works well with the additive model except when allele frequencies at the two loci differ widely. We explore those features of the additive model that make detecting interaction more difficult. All testing of this method suggests that it provides a reliable approach to detecting gene-gene interaction.     

Ectomycorrhizal fungal exoenzyme activity differs on spruce seedlings planted in forests versus clearcuts

Key message

Ectomycorrhizal (ECM) fungal community structure and potential exoenzymatic activity change after clearcut harvesting, but functional complementarity and redundancy among those ECM fungal species remaining support growth of regenerating seedlings.

Abstract

Ectomycorrhizal (ECM) fungal community composition is altered by forest harvesting, but it is not clear if this shift in structure influences ECM fungal physiological function at the community level. In this study, we characterized activities of extracellular enzymes in the ectomycorrhizospheres of Picea engelmannii seedlings grown in forest and clearcut plots. These exoenzymes are critical for the breakdown of large organic molecules, from which nutrients are subsequently absorbed and translocated by ECM fungi to host plants. We found that ectomycorrhizae on seedlings planted in forests had different exoenzyme activity profiles than those on seedlings planted in clearcuts. Specifically, the activities of glucuronidase, laccase, and acid phosphatase were higher on forest seedlings (P ≤ 0.006). These differences may have been partly driven by soil properties. Total carbon, total nitrogen (N), extractable phosphorus, extractable ammonium-N, and mineralizable N were higher, while pH was lower in forest plots (P ≤ 0.01). However, we also found that enzyme activity only shifted where community composition also changed. Functional complementarity can be inferred within ECM fungal communities in both forests and clearcuts because ectomycorrhizae formed by different species in the same environment had distinct enzyme profiles (P < 0.0001). However, ectomycorrhizae of Thelephora terrestris exhibited high levels of N- and P-mobilizing exoenzyme activities. Seedling biomass did not differ between forest and clearcut environments, so the high abundance of T. terrestris ectomycorrhizae in the clearcuts may have sustained nutrient acquisition by clearcut seedlings even in soils with lower N and P and with reduced ECM fungal species richness.

     

The response of cerambycid beetles (Coleoptera: Cerambycidae) to long‐term fire frequency regimes in subtropical eucalypt forest

Fire has a varied influence on plant and animal species through direct (e.g. fire‐induced mortality) and indirect (e.g. modification of habitat) effects. Our understanding of the influence of fire regime on invertebrates and their response to fire‐induced modifications to habitat is poor. We aimed to determine the response of a beetle family (Coleoptera: Cerambycidae) to varying fire treatments and hypothesised that the abundance of cerambycid beetles is influenced by fire frequency due to modifications in habitat associated with the fire treatments. Arthropods were sampled across 3 months in annually and triennially burnt areas (treatments starting in 1952 and 1973 respectively), an area unburnt since 1946, and a former unburnt treatment, burnt by wildfire in 2006. Eleven different cerambycid taxa were collected using flight intercept panel traps, dominated by three species (Ipomoria tillides, Adrium sp. and Bethelium signiferum) which made up 99% of individuals collected. Over the sampling period the long unburnt treatment had significantly lower species richness than the triennial and wildfire treatments. Cerambycid abundance was significantly higher in the triennially burnt treatment than in all other fire treatments. Ipomoria tillides was more abundant in both frequently burnt treatments, Adrium sp. was more common in triennially burnt areas, whereas B. signiferum, was more common in the wildfire affected treatment. Some, but not all, cerambycid beetles were more common in areas with a more open understorey (i.e. resulting from frequent burning), and lower tree basal area, as this likely influences their ability to fly easily between food sources. Cerambycid abundance was positively related to the volume of coarse woody debris and healthy tree crowns. Cerambycid beetles were clearly influenced by historic fire regime, suggesting that changes in fire regime can potentially have a profound influence on arthropod assemblages, and subsequent influences on ecosystem processes, which are currently poorly understood.     

Redefining the role of Ca2+-permeable channels in photoreceptor degeneration using diltiazem

Hereditary degeneration of photoreceptors has been linked to over-activation of Ca²⁺-permeable channels, excessive Ca²⁺-influx, and downstream activation of Ca²⁺-dependent calpain-type proteases. Unfortunately, after more than 20 years of pertinent research, unequivocal evidence proving significant and reproducible photoreceptor protection with Ca²⁺-channel blockers is still lacking. Here, we show that both D- and L-cis enantiomers of the anti-hypertensive drug diltiazem were very effective at blocking photoreceptor Ca²⁺-influx, most probably by blocking the pore of Ca²⁺-permeable channels. Yet, unexpectedly, this block neither reduced the activity of calpain-type proteases, nor did it result in photoreceptor protection. Remarkably, application of the L-cis enantiomer of diltiazem even led to a strong increase in photoreceptor cell death. These findings shed doubt on the previously proposed links between Ca²⁺ and retinal degeneration and are highly relevant for future therapy development as they may serve to refocus research efforts towards alternative, Ca²⁺-independent degenerative mechanisms.Subject terms: Cell death in the nervous system, Ion channels in the nervous system, Molecular neuroscience     

Influence of pH on the human prion protein: insights into the early steps of misfolding

Transmissible spongiform encephalopathies, or prion diseases, are caused by misfolding and aggregation of the prion protein PrP. Conversion from the normal cellular form (PrP^C) or recombinant PrP (recPrP) to a misfolded form is pH-sensitive, in that misfolding and aggregation occur more readily at lower pH. To gain more insight into the influence of pH on the dynamics of PrP and its potential to misfold, we performed extensive molecular-dynamics simulations of the recombinant PrP protein (residues 90-230) in water at three different pH regimes: neutral (or cytoplasmic) pH (∼7.4), middle (or endosomal) pH (∼5), and low pH (<4). We present five different simulations of 50 ns each for each pH regime, amounting to a total of 750 ns of simulation time. A detailed analysis and comparison with experiment validate the simulations and lead to new insights into the mechanism of pH-induced misfolding. The mobility of the globular domain increases with decreasing pH, through displacement of the first helix and instability of the hydrophobic core. At middle pH, conversion to a misfolded (PrP^Sc-like) conformation is observed. The observed changes in conformation and stability are consistent with experimental data and thus provide a molecular basis for the initial steps in the misfolding process.     

Calcium-independent phospholipase A2 and apoptosis

Apoptosis or programmed cell death is associated with changes in glycerophospholipid metabolism. Cells undergoing apoptosis generally release free fatty acids including arachidonic acid, which parallels the reduction in cell viability. The involvement of cytosolic group IVA phospholipase A(2)alpha (cPLA(2)alpha) in apoptosis has been the subject of numerous studies but a clear picture of the role(s) played by this enzyme is yet to emerge. More recently, the importance of lipid products generated by the action of a second phospholipase A(2), the group VIA calcium-independent phospholipase A(2) (iPLA(2)-VIA) in apoptosis has begun to be unveiled. Current evidence suggests that iPLA(2)-VIA-derived lysophosphatidylcholine may play a prominent role in mediating the chemoattractant and recognition/engulfment signals that accompany the process of apoptotic cell death, and gives possibility to the efficient clearance of dying cells by circulating phagocytes. Other lines of evidence suggest that perturbations in the control of free arachidonic acid levels within the cells, a process that may implicate iPLA(2)-VIA as well, may provide important cellular signals for the onset of apoptosis.     

Tandem affinity purification tagging of fatty acid biosynthetic enzymes in Synechocystis sp. PCC6803 and Arabidopsis thaliana

De novo fatty acid synthesis in plants occurs primarily in the plastids and is catalysed by a type-II fatty acid synthase (FAS) in which separate enzymes catalyse sequential reactions. Genes encoding all of the plant FAS components have been identified, following enzyme purification or by homology to Escherichia coli genes, and the structure of a number of the individual proteins determined. There are several lines of biochemical evidence indicating that FAS enzymes form a multi-protein complex and both in vitro and in vivo strategies can be used to investigate the association and interactions between them. To investigate protein interactions in vivo, tandem affinity purification-tagged FAS components are being used to purify complexes from both Arabidopsis thaliana and Synechocystis PCC6803. Here, the development of the tandem affinity purification method, its modification, and its use in plants is described and the experimental results achieved so far are reported.