Comparison of multiplex ligation-dependent probe amplification and real-time PCR accuracy for gene copy number quantification using the beta-defensin locus

The reliable quantification of gene copy number variations is a precondition for future investigations regarding their functional relevance. To date, there is no generally accepted gold standard method for copy number quantification, and methods in current use have given inconsistent results in selected cohorts. In this study, we compare two methods for copy number quantification. beta-defensin gene copy numbers were determined in parallel in 80 genomic DNA samples by real-time PCR and multiplex ligation-dependent probe amplification (MLPA). The pyrosequencing-based paralog ratio test (PPRT) was used as a standard of comparison in 79 out of 80 samples. Realtime PCR and MPLA results confirmed concordant DEFB4, DEFB103A, and DEFB104A copy numbers within samples. These two methods showed identical results in 32 out of 80 samples; 29 of these 32 samples comprised four or fewer copies. The coefficient of variation of MLPA is lower compared with PCR. In addition, the consistency between MLPA and PPRT is higher than either PCR/MLPA or PCR/PPRT consistency. In summary, these results suggest that MLPA is superior to real-time PCR in beta-defensin copy number quantification.

Changes in species composition in alpine snowbeds with climate change inferred from small-scale spatial patterns

Alpine snowbeds are characterised by a very short growing season. However, the length of the snow-free period is increasingly prolonged due to climate change, so that snowbeds become susceptible to invasions from neighbouring alpine meadow communities. We hypothesised that spatial distribution of species generated by plant interactions may indicate whether snowbed species will coexist with or will be out-competed by invading alpine species – spatial aggregation or segregation will point to coexistence or competitive exclusion, respectively. We tested this hypothesis in snowbeds of the Swiss Alps using the variance ratio statistics. We focused on the relationships between dominant snowbed species, subordinate snowbed species, and potentially invading alpine grassland species. Subordinate snowbed species were generally spatially aggregated with each other, but were segregated from alpine grassland species. Competition between alpine grassland and subordinate snowbed species may have caused this segregation. Segregation between these species groups increased with earlier snowmelt, suggesting an increasing importance of competition with climate change. Further, a dominant snowbed species (Alchemilla pentaphyllea) was spatially aggregated with subordinate snowbed species, while two other dominants (Gnaphalium supinum and Salix herbacea) showed aggregated patterns with alpine grassland species. These dominant species are known to show distinct microhabitat preferences suggesting the existence of hidden microhabitats with different susceptibility to invaders. These results allow us to suggest that alpine snowbed areas are likely to be reduced as a consequence of climate change and that invading species from nearby alpine grasslands could outcompete subordinate snowbed species. On the other hand, microhabitats dominated by Gnaphalium or Salix seem to be particularly prone to invasions by non-snowbed species.

Comparative phylogenies of the housekeeping genes atpD, infB and rpoB and the 16S rRNA gene within the Pasteurellaceae

Phylogenies of housekeeping gene and 16S rRNA gene sequences were compared to improve the classification of the bacterial family Pasteurellaceae and knowledge of the evolutionary relationships of its members. Deduced partial protein sequences of the housekeeping genes atpD, infB and rpoB were compared in 28, 36 and 28 representative taxa of the Pasteurellaceae, respectively. The monophyly of representatives of the genus Gallibacterium was recognized by analysis of all housekeeping genes, while members of Mannheimia, Actinobacillus sensu stricto and the core group of Pasteurella sensu stricto formed monophyletic groups with two out of three housekeeping genes. Representatives of Mannheimia, Actinobacillus sensu stricto, [Haemophilus] ducreyi and [Pasteurella] trehalosi formed a monophyletic unit by analysis of all three housekeeping genes, which was in contrast to the 16S rRNA gene-derived phylogeny, where these taxa occurred at separate positions in the phylogenetic tree. Representatives of the Rodent, Avian and Aphrophilus-Haemophilus 16S rRNA gene groups were weakly supported by phylogenetic analysis of housekeeping genes. Phylogenies derived by comparison of the housekeeping genes diverged significantly from the 16S rRNA gene-derived phylogeny as evaluated by the likelihood ratio test. A low degree of congruence was also observed between the individual housekeeping gene-derived phylogenies. Estimates on speciation derived from 16S rRNA and housekeeping gene sequence comparisons resulted in quite different evolutionary scenarios for members of the Pasteurellaceae. The phylogeny based on the housekeeping genes supported observed host associations between Mannheimia, Actinobacillus sensu stricto and [Pasteurella] trehalosi and animals with paired hooves.

Copy number variation of the Beta-defensin genes in europeans: no supporting evidence for association with lung function, chronic obstructive pulmonary disease or asthma

Lung function measures are heritable, predict mortality and are relevant in diagnosis of chronic obstructive pulmonary disease (COPD). COPD and asthma are diseases of the airways with major public health impacts and each have a heritable component. Genome-wide association studies of SNPs have revealed novel genetic associations with both diseases but only account for a small proportion of the heritability. Complex copy number variation may account for some of the missing heritability. A well-characterised genomic region of complex copy number variation contains beta-defensin genes (DEFB103, DEFB104 and DEFB4), which have a role in the innate immune response. Previous studies have implicated these and related genes as being associated with asthma or COPD. We hypothesised that copy number variation of these genes may play a role in lung function in the general population and in COPD and asthma risk. We undertook copy number typing of this locus in 1149 adult and 689 children using a paralogue ratio test and investigated association with COPD, asthma and lung function. Replication of findings was assessed in a larger independent sample of COPD cases and smoking controls. We found evidence for an association of beta-defensin copy number with COPD in the adult cohort (OR = 1.4, 95%CI:1.02-1.92, P = 0.039) but this finding, and findings from a previous study, were not replicated in a larger follow-up sample(OR = 0.89, 95%CI:0.72-1.07, P = 0.217). No robust evidence of association with asthma in children was observed. We found no evidence for association between beta-defensin copy number and lung function in the general populations. Our findings suggest that previous reports of association of beta-defensin copy number with COPD should be viewed with caution. Suboptimal measurement of copy number can lead to spurious associations. Further beta-defensin copy number measurement in larger sample sizes of COPD cases and children with asthma are needed.

The influence of body mass index on injury pattern in polytrauma: thorax as the main source of complications

Background: Obesity is a growing problem in industrial nations. The aim of this study was to determine the relationship between the body mass index (BMI) and the pattern of injury after polytrauma. Methods: This retrospective study included 651 patients with an injury severity score (ISS) ≥16 and aged ≥16 years who were subdivided into three groups: BMI < 25 kg/m2, BMI 25–30 kg/m2, and BMI > 30 kg/m2. The Abbreviated Injury Scale (AIS) was used to quantify the injuries in the different anatomical regions. The Murray score was assessed at admission and at its maximum during hospitalization to evaluate pulmonary problems. Data are presented as means ± standard errors of the means. One way analysis of variance, χ2 test and Kruskal-Wallis test were used for the analyses and the significance level was set at p < 0.05. Results: The AIS of the thorax was 3.2 ± 0.1 in the BMI < 25 kg/m2 group, 3.3 ± 0.1 in the BMI 25–30 kg/m2 group, and 2.8 ± 0.2 in the BMI > 30 kg/m2 group; p < 0.05. The Murray score at admission increased significantly with increasing BMI (0.8 ± 0.8 for BMI < 25 kg/m2, 0.9 ± 0.9 for BMI 25–30 kg/m2, and 1.0 ± 0.8 for BMI > 30 kg/m2; p < 0.05) as was the maximum Murray score during hospitalization (1.2 ± 0.9 for BMI < 25 kg/m2, 1.6 ± 1.0 for BMI 25–30 kg/m2, and 1.5 ± 0.9 for BMI > 30 kg/m2; p < 0.001). The number of ventilator days was also elevated significantly with increasing BMI (5.9 ± 0.4 for BMI < 25 kg/m2, 7.7 ± 0.8 for BMI 25–30 kg/m2, and 7.9 ± 1.6 for BMI > 30 kg/m2; p < 0.05). Conclusion: Overweight and obesity lead to a higher incidence of thoracic trauma in a polytrauma situation and may additionally handicap ventilation in an obstructive manner.

Adaptive Molecular Evolution in the Opsin Genes of Rapidly Speciating Cichlid Species

Cichlid fish inhabit a diverse range of environments that vary in the spectral content of light available for vision. These differences should result in adaptive selective pressure on the genes involved in visual sensitivity, the opsin genes. This study examines the evidence for differential adaptive molecular evolution in East African cichlid opsin genes due to gross differences in environmental light conditions. First, we characterize the selective regime experienced by cichlid opsin genes using a likelihood ratio test format, comparing likelihood models with different constraints on the relative rates of amino acid substitution, across sites. Second, we compare turbid and clear lineages to determine if there is evidence of differences in relative rates of substitution. Third, we present evidence of functional diversification and its relationship to the photic environment among cichlid opsin genes. We report statistical evidence of positive selection in all cichlid opsin genes, except short wavelength–sensitive 1 and short wavelength–sensitive 2b. In all genes predicted to be under positive selection, except short wavelength–sensitive 2a, we find differences in selective pressure between turbid and clear lineages. Potential spectral tuning sites are variable among all cichlid opsin genes; however, patterns of substitution consistent with photic environment–driven evolution of opsin genes are observed only for short wavelength–sensitive 1 opsin genes. This study identifies a number of promising candidate-tuning sites for future study by site-directed mutagenesis. This work also begins to demonstrate the molecular evolutionary dynamics of cichlid visual sensitivity and its relationship to the photic environment.

Cost and power analysis of linkage disequilibrium methods for locating genes influencing quantitative traits in humans

Linkage disequilibrium methods can be used to find genes influencing quantitative trait variation in humans. Linkage disequilibrium methods can require smaller sample sizes than linkage equilibrium methods, such as the variance component approach to find loci with a specific effect size. The increase in power is at the expense of requiring more markers to be typed to scan the entire genome. This thesis compares different linkage disequilibrium methods to determine which factors influence the power to detect disequilibrium. The costs of disequilibrium and equilibrium tests were compared to determine whether the savings in phenotyping costs when using disequilibrium methods outweigh the additional genotyping costs.^ Nine linkage disequilibrium tests were examined by simulation. Five tests involve selecting isolated unrelated individuals while four involved the selection of parent child trios (TDT). All nine tests were found to be able to identify disequilibrium with the correct significance level in Hardy-Weinberg populations. Increasing linked genetic variance and trait allele frequency were found to increase the power to detect disequilibrium, while increasing the number of generations and distance between marker and trait loci decreased the power to detect disequilibrium. Discordant sampling was used for several of the tests. It was found that the more stringent the sampling, the greater the power to detect disequilibrium in a sample of given size. The power to detect disequilibrium was not affected by the presence of polygenic effects.^ When the trait locus had more than two trait alleles, the power of the tests maximized to less than one. For the simulation methods used here, when there were more than two-trait alleles there was a probability equal to 1-heterozygosity of the marker locus that both trait alleles were in disequilibrium with the same marker allele, resulting in the marker being uninformative for disequilibrium.^ The five tests using isolated unrelated individuals were found to have excess error rates when there was disequilibrium due to population admixture. Increased error rates also resulted from increased unlinked major gene effects, discordant trait allele frequency, and increased disequilibrium. Polygenic effects did not affect the error rates. The TDT, Transmission Disequilibrium Test, based tests were not liable to any increase in error rates.^ For all sample ascertainment costs, for recent mutations ($<$100 generations) linkage disequilibrium tests were less expensive than the variance component test to carry out. Candidate gene scans saved even more money. The use of recently admixed populations also decreased the cost of performing a linkage disequilibrium test. ^

The impact of body mass index and gender on the development of infectious complications in polytrauma patients

Purpose The aim was to test the impact of body mass index (BMI) and gender on infectious complications after polytrauma. Methods A total of 651 patients were included in this retrospective study, with an Injury Severity Score (ISS) C16 and age C16 years. The sample was subdivided into three groups: BMI\25 kg/m2, BMI 25–30 kg/m2, and BMI[30 kg/m2, and a female and a male group. Infectious complications were observed for 31 days after admission. Data are given as mean ± standard errors of the means. Analysis of variance, Kruskal–Wallis test, v2 tests, and Pearson’s correlation were used for the analyses and the significance level was set at P\0.05. Results The overall infection rates were 31.0 % in the BMI\25 kg/m2 group, 29.0 % in the BMI 25–30 kg/m2 group, and 24.5 % in the BMI[30 kg/m2 group (P = 0.519). The female patients developed significantly fewer infectious complications than the male patients (26.8 vs. 73.2 %; P\0.001). The incidence of death was significantly decreased according to the BMI group (8.8 vs. 7.2 vs. 1.5 %; P\0.0001) and the female population had a significantly lower mortality rate (4.1 vs. 13.4 %; P\0.0001). Pearson’s correlations between the Abbreviated Injury Scale (AIS) score and the corresponding infectious foci were not significant. Conclusion Higher BMI seems to be protective against polytrauma-associated death but not polytrauma-associated infections, and female gender protects against both polytrauma- associated infections and death. Understanding gender-specific immunomodulation could improve the outcome of polytrauma patients.

Genome-wide association study of survival from sepsis due to pneumonia: an observational cohort study.

BACKGROUND Sepsis continues to be a major cause of death, disability, and health-care expenditure worldwide. Despite evidence suggesting that host genetics can influence sepsis outcomes, no specific loci have yet been convincingly replicated. The aim of this study was to identify genetic variants that influence sepsis survival. METHODS We did a genome-wide association study in three independent cohorts of white adult patients admitted to intensive care units with sepsis, severe sepsis, or septic shock (as defined by the International Consensus Criteria) due to pneumonia or intra-abdominal infection (cohorts 1-3, n=2534 patients). The primary outcome was 28 day survival. Results for the cohort of patients with sepsis due to pneumonia were combined in a meta-analysis of 1553 patients from all three cohorts, of whom 359 died within 28 days of admission to the intensive-care unit. The most significantly associated single nucleotide polymorphisms (SNPs) were genotyped in a further 538 white patients with sepsis due to pneumonia (cohort 4), of whom 106 died. FINDINGS In the genome-wide meta-analysis of three independent pneumonia cohorts (cohorts 1-3), common variants in the FER gene were strongly associated with survival (p=9·7 × 10(-8)). Further genotyping of the top associated SNP (rs4957796) in the additional cohort (cohort 4) resulted in a combined p value of 5·6 × 10(-8) (odds ratio 0·56, 95% CI 0·45-0·69). In a time-to-event analysis, each allele reduced the mortality over 28 days by 44% (hazard ratio for death 0·56, 95% CI 0·45-0·69; likelihood ratio test p=3·4 × 10(-9), after adjustment for age and stratification by cohort). Mortality was 9·5% in patients carrying the CC genotype, 15·2% in those carrying the TC genotype, and 25·3% in those carrying the TT genotype. No significant genetic associations were identified when patients with sepsis due to pneumonia and intra-abdominal infection were combined. INTERPRETATION We have identified common variants in the FER gene that associate with a reduced risk of death from sepsis due to pneumonia. The FER gene and associated molecular pathways are potential novel targets for therapy or prevention and candidates for the development of biomarkers for risk stratification. FUNDING European Commission and the Wellcome Trust.

Effectiveness of meticillin-resistant Staphylococcus aureus decolonization in long-term haemodialysis patients: a systematic review and meta-analysis.

BACKGROUND Chronic haemodialysis patients are a high-risk population for meticillin-resistant Staphylococcus aureus (MRSA) colonization, which is a precursor of infection. AIM To summarize the effect of nasal (± whole-body wash) MRSA decolonization in haemodialysis patients by means of a systematic review and meta-analysis. METHODS We identified eligible studies using Medline, Embase, the Cochrane database, clinicaltrials.org, and conference abstracts investigating the success of MRSA decolonization in haemodialysis patients. For the statistical analysis, we used Stata 13 to express study-specific proportions with 95% confidence intervals. A likelihood ratio test was used to assess inter-study heterogeneity. FINDINGS Six published prospective cohort studies and one study described in a conference abstract met our inclusion criteria. From 1150 haemodialysis patients enrolled in these studies, MRSA was isolated from nasal swabs of 147 (12.8%) patients. Six of the trials used mupirocin nasal ointment and combined it with chlorhexidine body washes for decolonization. The most widely used protocol was a five-day course of mupirocin nasal ointment application three times a day, and chlorhexidine body wash once daily. The pooled success rate of decolonization was 0.88 (95% confidence interval: 0.75-0.95). A likelihood ratio test of the fixed versus the random-effects model showed significant inter-study heterogeneity (P = 0.047). Four of seven studies determined subsequent MRSA infections in 94 carriers overall, two (2%) of which experienced infection. CONCLUSION The use of mupirocin together with whole-body decolonization is highly effective in eradicating MRSA carriage in haemodialysis patients. The current literature, however, is characterized by a lack of comparative effectiveness studies for this intervention.

DNA nucleotide excision repair gene single nucleotide polymorphisms and hereditary nonpolyposis colorectal cancer

Hereditary nonpolyposis colorectal cancer (HNPCC) is an autosomal dominant disease caused by germline mutations in DNA mismatch repair(MMR) genes. The nucleotide excision repair(NER) pathway plays a very important role in cancer development. We systematically studied interactions between NER and MMR genes to identify NER gene single nucleotide polymorphism (SNP) risk factors that modify the effect of MMR mutations on risk for cancer in HNPCC. We analyzed data from polymorphisms in 10 NER genes that had been genotyped in HNPCC patients that carry MSH2 and MLH1 gene mutations. The influence of the NER gene SNPs on time to onset of colorectal cancer (CRC) was assessed using survival analysis and a semiparametric proportional hazard model. We found the median age of onset for CRC among MMR mutation carriers with the ERCC1 mutation was 3.9 years earlier than patients with wildtype ERCC1(median 47.7 vs 51.6, log-rank test p=0.035). The influence of Rad23B A249V SNP on age of onset of HNPCC is age dependent (likelihood ratio test p=0.0056). Interestingly, using the likelihood ratio test, we also found evidence of genetic interactions between the MMR gene mutations and SNPs in ERCC1 gene(C8092A) and XPG/ERCC5 gene(D1104H) with p-values of 0.004 and 0.042, respectively. An assessment using tree structured survival analysis (TSSA) showed distinct gene interactions in MLH1 mutation carriers and MSH2 mutation carriers. ERCC1 SNP genotypes greatly modified the age onset of HNPCC in MSH2 mutation carriers, while no effect was detected in MLH1 mutation carriers. Given the NER genes in this study play different roles in NER pathway, they may have distinct influences on the development of HNPCC. The findings of this study are very important for elucidation of the molecular mechanism of colon cancer development and for understanding why some mutation carriers of the MSH2 and MLH1 gene develop CRC early and others never develop CRC. Overall, the findings also have important implications for the development of early detection strategies and prevention as well as understanding the mechanism of colorectal carcinogenesis in HNPCC. ^

Exploration of genetic susceptibility to pancreatic cancer: A GWAS approach

Pancreatic cancer is the 4th most common cause for cancer death in the United States, accompanied by less than 5% five-year survival rate based on current treatments, particularly because it is usually detected at a late stage. Identifying a high-risk population to launch an effective preventive strategy and intervention to control this highly lethal disease is desperately needed. The genetic etiology of pancreatic cancer has not been well profiled. We hypothesized that unidentified genetic variants by previous genome-wide association study (GWAS) for pancreatic cancer, due to stringent statistical threshold or missing interaction analysis, may be unveiled using alternative approaches. To achieve this aim, we explored genetic susceptibility to pancreatic cancer in terms of marginal associations of pathway and genes, as well as their interactions with risk factors. We conducted pathway- and gene-based analysis using GWAS data from 3141 pancreatic cancer patients and 3367 controls with European ancestry. Using the gene set ridge regression in association studies (GRASS) method, we analyzed 197 pathways from the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Using the logistic kernel machine (LKM) test, we analyzed 17906 genes defined by University of California Santa Cruz (UCSC) database. Using the likelihood ratio test (LRT) in a logistic regression model, we analyzed 177 pathways and 17906 genes for interactions with risk factors in 2028 pancreatic cancer patients and 2109 controls with European ancestry. After adjusting for multiple comparisons, six pathways were marginally associated with risk of pancreatic cancer ( P < 0.00025): Fc epsilon RI signaling, maturity onset diabetes of the young, neuroactive ligand-receptor interaction, long-term depression (Ps < 0.0002), and the olfactory transduction and vascular smooth muscle contraction pathways (P = 0.0002; Nine genes were marginally associated with pancreatic cancer risk (P < 2.62 × 10−5), including five reported genes (ABO, HNF1A, CLPTM1L, SHH and MYC), as well as four novel genes (OR13C4, OR 13C3, KCNA6 and HNF4 G); three pathways significantly interacted with risk factors on modifying the risk of pancreatic cancer (P < 2.82 × 10−4): chemokine signaling pathway with obesity ( P < 1.43 × 10−4), calcium signaling pathway (P < 2.27 × 10−4) and MAPK signaling pathway with diabetes (P < 2.77 × 10−4). However, none of the 17906 genes tested for interactions survived the multiple comparisons corrections. In summary, our current GWAS study unveiled unidentified genetic susceptibility to pancreatic cancer using alternative methods. These novel findings provide new perspectives on genetic susceptibility to and molecular mechanisms of pancreatic cancer, once confirmed, will shed promising light on the prevention and treatment of this disease. ^

Fuste comercial de clones del género Populus : relación entre volumen y densidad de la madera con nutrientes inorgánicos

El objetivo del trabajo fue determinar sobre 10 clones de álamos, provenientes de cruzamientos intraespecíficos de Populus deltoides e interespecíficos de Populus x canadensis, implantados en Teodelina, Santa Fe, Argentina (34°09' S; 61°15' W), la concentración de nutrientes inorgánicos en su corteza y leño, y sus relaciones con los volúmenes comerciales alcanzados a los 9 años de edad y la densidad básica de la madera. Se apearon ejemplares seleccionados, se calculó su volumen comercial individual y se determinaron la densidad básica y las concentraciones de Ca, Mg, P y K en fuste y corteza. Se realizó análisis de la varianza, test de comparación de medias de Tukey y correlación simple de los diferentes parámetros evaluados. Los resultados arrojaron diferentes contenidos de nutrientes por clon y concentraciones más bajas en la madera que en la corteza, en todos los casos. No hubo correlaciones significativas entre el volumen comercial obtenido, los valores de densidad básica de la madera y los nutrientes inorgánicos, indicando eficiencias individuales en el uso de los recursos inorgánicos.

Homogenization of Portuguese long-term temperature data series: Lisbon, Coimbra and Porto

Three long-term temperature data series measured in Portugal were studied to detect and correct non-climatic homogeneity breaks and are now available for future studies of climate variability. Series of monthly minimum (Tmin) and maximum (Tmax) temperatures measured in the three Portuguese meteorological stations of Lisbon (from 1856 to 2008), Coimbra (from 1865 to 2005) and Porto (from 1888 to 2001) were studied to detect and correct non-climatic homogeneity breaks. These series together with monthly series of average temperature (Taver) and temperature range (DTR) derived from them were tested in order to detect homogeneity breaks, using, firstly, metadata, secondly, a visual analysis and, thirdly, four widely used homogeneity tests: von Neumann ratio test, Buishand test, standard normal homogeneity test and Pettitt test. The homogeneity tests were used in absolute (using temperature series themselves) and relative (using sea-surface temperature anomalies series obtained from HadISST2 close to the Portuguese coast or already corrected temperature series as reference series) modes. We considered the Tmin, Tmax and DTR series as most informative for the detection of homogeneity breaks due to the fact that Tmin and Tmax could respond differently to changes in position of a thermometer or other changes in the instrument's environment; Taver series have been used, mainly, as control. The homogeneity tests show strong inhomogeneity of the original data series, which could have both internal climatic and non-climatic origins. Homogeneity breaks which have been identified by the last three mentioned homogeneity tests were compared with available metadata containing data, such as instrument changes, changes in station location and environment, observing procedures, etc. Significant homogeneity breaks (significance 95% or more) that coincide with known dates of instrumental changes have been corrected using standard procedures. It was also noted that some significant homogeneity breaks, which could not be connected to the known dates of any changes in the park of instruments or stations location and environment, could be caused by large volcanic eruptions. The corrected series were again tested for homogeneity: the corrected series were considered free of non-climatic breaks when the tests of most of monthly series showed no significant (significance 95% or more) homogeneity breaks that coincide with dates of known instrument changes. Corrected series are now available in the frame of ERA-CLIM FP7 project for future studies of climate variability.

Maintenance of pre-mRNA secondary structure by epistatic selection.

Linkage disequilibrium between polymorphisms in a natural population may result from various evolutionary forces, including random genetic drift due to sampling of gametes during reproduction, restricted migration between subpopulations in a subdivided population, or epistatic selection. In this report, we present evidence that the majority of significant linkage disequilibria observed in introns of the alcohol dehydrogenase locus (Adh) of Drosophila pseudoobscura are due to epistatic selection maintaining secondary structure of precursor mRNA (pre-mRNA). Based on phylogenetic-comparative analysis and a likelihood approach, we propose secondary structure models of Adh pre-mRNA for the regions of the adult intron and intron 2 where clustering of linkage disequilibria has been observed. Furthermore, we applied the likelihood ratio test to the phylogenetically predicted secondary structure in intron 1. In contrast to the other two structures, polymorphisms associated with the more conserved stem-loop structure of intron 1 are in low frequency, and linkage disequilibria have not been observed. These findings are qualitatively consistent with a model of compensatory fitness interactions. This model assumes that mutations disrupting pairing in a secondary structural element are individually deleterious if they destabilize a functionally important structure; a second "compensatory" mutation, however, may restabilize the structure and restore fitness.