Is the hemochromatosis gene a modifier locus for cystic fibrosis?

Interaction models matter: an efficient, flexible computational framework for model-specific investigation of epistasis

PURPOSE

Epistasis, the interaction between two or more genes, is integral to the study of genetics and is present throughout nature. Yet, it is seldom fully explored as most approaches primarily focus on single-locus effects, partly because analyzing all pairwise and higher-order interactions requires significant computational resources. Furthermore, existing methods for epistasis detection only consider a Cartesian (multiplicative) model for interaction terms. This is likely limiting as epistatic interactions can evolve to produce varied relationships between genetic loci, some complex and not linearly separable.

METHODS

We present new algorithms for the interaction coefficients for standard regression models for epistasis that permit many varied models for the interaction terms for loci and efficient memory usage. The algorithms are given for two-way and three-way epistasis and may be generalized to higher order epistasis. Statistical tests for the interaction coefficients are also provided. We also present an efficient matrix based algorithm for permutation testing for two-way epistasis. We offer a proof and experimental evidence that methods that look for epistasis only at loci that have main effects may not be justified. Given the computational efficiency of the algorithm, we applied the method to a rat data set and mouse data set, with at least 10,000 loci and 1,000 samples each, using the standard Cartesian model and the XOR model to explore body mass index.

RESULTS

This study reveals that although many of the loci found to exhibit significant statistical epistasis overlap between models in rats, the pairs are mostly distinct. Further, the XOR model found greater evidence for statistical epistasis in many more pairs of loci in both data sets with almost all significant epistasis in mice identified using XOR. In the rat data set, loci involved in epistasis under the XOR model are enriched for biologically relevant pathways.

CONCLUSION

Our results in both species show that many biologically relevant epistatic relationships would have been undetected if only one interaction model was applied, providing evidence that varied interaction models should be implemented to explore epistatic interactions that occur in living systems.

Meconium ileus caused by mutations in GUCY2C, encoding the CFTR-activating guanylate cyclase 2C

Meconium ileus, intestinal obstruction in the newborn, is caused in most cases by CFTR mutations modulated by yet-unidentified modifier genes. We now show that in two unrelated consanguineous Bedouin kindreds, an autosomal-recessive phenotype of meconium ileus that is not associated with cystic fibrosis (CF) is caused by different homozygous mutations in GUCY2C, leading to a dramatic reduction or fully abrogating the enzymatic activity of the encoded guanlyl cyclase 2C. GUCY2C is a transmembrane receptor whose extracellular domain is activated by either the endogenous ligands, guanylin and related peptide uroguanylin, or by an external ligand, Escherichia coli (E. coli) heat-stable enterotoxin STa. GUCY2C is expressed in the human intestine, and the encoded protein activates the CFTR protein through local generation of cGMP. Thus, GUCY2C is a likely candidate modifier of the meconium ileus phenotype in CF. Because GUCY2C heterozygous and homozygous mutant mice are resistant to E. coli STa enterotoxin-induced diarrhea, it is plausible that GUCY2C mutations in the desert-dwelling Bedouin kindred are of selective advantage.

Physiology and pathology of proteostasis in the early secretory compartment

The endoplasmic reticulum (ER), the port of entry for proteins into the secretory pathway, is a multifunctional organelle emerging as a central integrator of numerous signalling pathways. The mechanisms that control proteostasis are integral part of this signalling network, providing cues for morphological and functional cell remodelling, proliferation, inflammation and cell death. The complexity of ER responses is exploited during physiological and pathological tissue development, cell differentiation and lifespan control. This essay outlines some of the mechanisms that link proteostasis within the early secretory compartment to signalling in development and disease.

Update on gene modifiers in cystic fibrosis

PURPOSE OF REVIEW

Cystic fibrosis (CF) is a common, life-limiting monogenic disease, which typically manifests as progressive bronchiectasis, exocrine pancreatic dysfunction, and recurrent sinopulmonary infections. Although the gene responsible for CF (CFTR) was described in 1989, it has become increasingly evident that modifier genes and environmental factors play substantial roles in determining the severity of disease, particularly lung disease. Identifying these factors is crucial in devising therapies and other interventions to decrease the morbidity and mortality associated with this disorder.

RECENT FINDINGS

Although many genes have been proposed as potential modifiers of CF, only a handful have withstood the test of replication. Several of the replicated findings reveal that genes affecting inflammation and infection response play a key role in modifying CF lung disease severity. Interactions between CFTR genotype, modifier genes, and environmental factors have been documented to influence lung function measures and infection status in CF patients.

SUMMARY

Several genes have been demonstrated to affect disease severity in CF. Furthermore, it is likely that gene-gene and gene-environment interactions can explain a substantial portion of the variation of lung disease. Ongoing genome-wide studies are likely to identify novel genetic modifiers. Continued exploration of the role of genetic and nongenetic modifiers of CF is likely to yield new options for combating this debilitating disease.

A family with atypical cystic fibrosis: brother and sister with heterozygosity for both G542X and R117H

Clinical manifestations of cystic fibrosis (CF) are variable. Genetic and environmental factors that determine whether an individual will develop associated complications are still under investigation. The present study reports the genetic analysis of a family with different clinical forms of CF and addresses the difficulty of CF diagnosis in an individual with mutant alleles G542X and R117H because of the variable phenotype associated with R117H mutation. Both children in this family were heterozygous for G542X/R117H with the same thymine sequence (7T/9T) in intron 8 of CF transmembrane conductance regulator. The girl was diagnosed with CF, whereas the boy was diagnosed with azoospermia as the sole clinical manifestation. The possible implication of the hemochromatosis gene as a CF modifier locus was analyzed because the 2 children had the same genotype. No genetic differences were detected between brother and sister that explained the different clinical manifestations of CF.

Differential expression of genes related to HFE and iron status in mouse duodenal epithelium

Iron absorption, distribution, use, and storage are thought to be tightly regulated since altered iron stores may lead to cellular damage and disease. HFE, the hereditary hemochromatosis gene product, is expressed in the crypts of the duodenum, but the molecular mechanism by which it contributes to the inhibition of iron absorption is still unknown. In this study we aimed to identify transcriptional profiles in the duodenal epithelium of Hfe(-/-) mice. We used dedicated microarrays to compare gene expression among the duodenum of Hfe(-/-) mice, induced iron overload mice, and control mice. We found 151 differentially expressed genes and unknown sequences between Hfe(-/-) mice and normal littermates. Gene profiling revealed a gene subset more specific for Hfe inactivation. The functional annotation of upregulated genes highlighted that mucus production and cell maintenance may account for the influence of Hfe on epithelium integrity and luminal iron uptake.

Epistasis and the genetics of human diseases

Epistasis or modifier genes, that is, gene-gene interactions of non-allelic partners, play a major role in susceptibility to common human diseases. This old genetic concept has experienced a major renaissance recently. Interestingly, epistatic genes can make the disease less severe, or make it more severe. Hence, most diseases are of different intensities in different individuals and in different ethnicities. This phenomenon affects sickle-cell anemia carriers and other hemoglobinopathies, systemic lupus erythematosus, cystic fibrosis, complex autoimmune diseases, venous thromboembolism, and many others. It is likely, and fortunate, than 20 years form now, patients entering a medical facility will be subjected to a genomic scanning, including pathogenic genes as well as epistatic genes.

Disease modifying genes in cystic fibrosis

The variation in cystic fibrosis (CF) lung disease and development of CF related complications correlates poorly with the genotype of the CF transmembrane regulator (CFTR) and with environmental factors. Increasing evidence suggests that phenotypic variation in CF can be attributed to genetic variation in genes other than the CFTR gene, so-called modifier genes. In recent years, multiple candidate modifier genes have been investigated in CF, especially genes that are involved in the control of infection, immunity and inflammation. Some of these genes have been rather conclusively identified as modifiers of the CF phenotype, whereas associations found in other genes have not been confirmed or are conflicting. Identification of genetic variation in modifier genes, obtained by genotype-phenotype studies in well-defined patient populations, may be used as an aid to prognosis and may provide the possibility of new therapeutic interventions.

Hepatobiliary complications of cystic fibrosis

Cystic fibrosis (CF) is the most common potentially lethal genetic disease in the Caucasian population. The disease results from mutations in the gene for the cystic fibrosis transmembrane conductance regulator (CFTR), a cAMP-activated Cl(-) channel in the apical membrane of most secretory epithelia. In the liver, CFTR is located in biliary epithelial cells or cholangiocytes and gallbladder epithelia, where it appears to play a role in normal bile formation. However, how a defective CFTR protein leads to associated liver and biliary disease in a subset of patients with CF is unknown. Improvements in life expectancy have led to an increasing recognition of hepatobiliary complications from CF. Whereas the biliary tract disease is usually clinically evident, the liver involvement may progress silently, only manifesting as end-stage liver disease and portal hypertension. Unlike the pancreatic involvement in CF, a genotype-phenotype correlation is not apparent in the expression of liver disease, suggesting the presence of as yet unidentifiable "genetic modifiers" influencing disease expression. This review focuses on the pathogenesis, clinical manifestations, screening, diagnosis, and treatment of CF hepatobiliary disease.

HFE alleles in an Irish cystic fibrosis population

The variable clinical manifestations of cystic fibrosis (CF) suggest the influence of modifier genes. Genetic and environmental factors that determine whether an individual will develop associated complications are still being determined. It has been proposed that the gene for hemochromatosis, HFE, may be a modifier locus for CF disease phenotype. Recent research has suggested a relationship between mutations to the HFE gene and the development of meconium ileus (MI) and liver disease in CF. This study aims to expand our knowledge of the HFE mutations C282Y and H63D carrier rate in an Irish population of CF allele carriers. PCR restriction enzyme analysis was performed on blood samples from CF patients to identify the C282Y and H63D mutations. HFE status of CF allele carriers and CF patients (Delta F508) homozygotes with and without meconium ileus was determined. The carrier frequency for C282Y was 30.8% for the Delta F508 homozygote MI positive group, as compared to 12.5% for the non-Delta F508 MI positive group but did not reach statistical significance (p = 0.27). Interestingly, no Delta F508 homozygote patients were homozygous for the C282Y mutation.

Analysis of the two common alpha-1-antitrypsin deficiency alleles PiMS and PiMZ as modifiers of Pseudomonas aeruginosa susceptibility in cystic fibrosis

Lung disease is the direct cause of death in more than 90% of cystic fibrosis (CF) patients. Proteinase-antiproteinase imbalances are common in CF and alpha-1-antitrypsin (AAT) deficiency. We investigated the hypothesis that the AAT deficiency alleles PiS and PiZ contribute to pulmonary prognosis in CF. Two hundred and sixty-nine CF patients from Southern Germany were included in this study. The serum concentrations of AAT and C-reactive protein (CRP) were determined by nephelometry, and patients were screened by polymerase chain reaction (PCR) and restriction enzyme digest for the common AAT deficiency alleles PiS and PiZ. The onset of chronic bacterial colonization by Pseudomonas aeruginosa (Pae) was correlated with the AAT phenotypes PiMM, PiMS and PiMZ. Only three out of nine CF patients (33%) diagnosed with either PiMS or PiMZ had developed chronic Pae lung infection earlier in their lives. The remaining six patients showing a PiMS or PiMZ phenotype showed a later onset of chronic Pae lung infection. Our results indicate that PiMS and PiMZ are not associated with worse pulmonary prognosis in CF. These data need to be confirmed in studies with a much larger number of cases.

Genotype-phenotype correlation in cystic fibrosis: the role of modifier genes

More than 1,000 mutations have been identified in the cystic fibrosis (CF) transmembrane regulator (CFTR) disease gene. The impact of these mutations on the protein and the wide spectrum of CF phenotypes prompted a series of Genotype-Phenotype correlation studies. The CFTR genotype is invariably correlated with pancreatic status-in about 85% of cases with pancreatic insufficiency and in about 15% of cases with pancreatic sufficiency. The correlations between the CFTR genotype and pulmonary, liver, and gastrointestinal expression are debatable. The heterogeneous phenotype in CF patients bearing the same genotype or homozygotes for nonsense mutations implicated environmental and/or genetic factors in the disease. However, the discordant phenotype observed in CF siblings argued against a major role of environmental factors and suggested that genes other than CFTR modulate the CF phenotype. A locus that modulates gastrointestinal expression was identified in mice and subsequently in humans. By analyzing nine CF patients discordant for meconium ileus we were able to show that this locus had a dominant effect. Moreover, in a collaborative study we found a higher rate of polymorphisms in beta-defensin genes 1 and 2 in CF patients and in controls. In another multicenter study mutations in alpha-1 antitrypsin (A1AT) and mannose binding lectin genes were found to be independent risk factors for liver disease in CF patients. The body of evidence available suggests that the variegated CF phenotype results from complex interactions between numerous gene products.

A novel CFTR frame-shift mutation, 935delA, in two Hispanic cystic fibrosis patients

The currently available mutation analysis panel detects about 50-60% of CFTR mutations in Hispanic patients. In order to search for Hispanic CF mutations, we developed a temporal temperature gradient gel electrophoresis (TTGE) method to screen for unknown mutations. Using TTGE to study the CFTR gene has lead to the discovery of many novel mutations in Hispanic patients. A novel frame-shift mutation, 935delA, was found in two unrelated patients. One was heterozygous for two novel frame-shift mutations, 663delT and 935delA, and the other was heterozygous for DeltaF508 and 935delA. Both patients showed severe phenotype with meconium ileus, pancreatic insufficiency, and early pulmonary microbial colonization with Pseudomonas aeruginosa. Patient 1 died at 4 years of age. Patient 2 had an upper lobectomy. The 935delA mutation produces a truncated polypeptide with only 21% of the full-length protein. The severe course of clinical manifestation is consistent with two oppressively truncated mutant polypeptides encoded by both mutant alleles in patient 1 and the compound heterozygosity truncation and DeltaF508 mutations in patient 2.