Statistical methods to detect mother-father genetic interaction effects on risk of infertility: A genome-wide approach.

Infertility is a heterogeneous phenotype, and for many couples, the causes of fertility problems remain unknown. One understudied hypothesis is that allelic interactions between the genotypes of the two parents may influence the risk of infertility. Our aim was, therefore, to investigate how allelic interactions can be modeled using parental genotype data linked to 15,789 pregnancies selected from the Norwegian Mother, Father, and Child Cohort Study. The newborns in 1304 of these pregnancies were conceived using assisted reproductive technologies (ART), and the remainder were conceived naturally. Treating the use of ART as a proxy for infertility, different parameterizations were implemented in a genome-wide screen for interaction effects between maternal and paternal alleles at the same locus. Some of the models were more similar in the way they were parameterized, and some produced similar results when implemented on a genome-wide scale. The results showed near-significant interaction effects in genes relevant to the phenotype under study, such as Dynein axonemal heavy chain 17 (DNAH17) with a recognized role in male infertility. More generally, the interaction models presented here are readily adaptable to the study of other phenotypes in which maternal and paternal allelic interactions are likely to be involved.

The X-factor in ART: does the use of assisted reproductive technologies influence DNA methylation on the X chromosome?

BACKGROUND: Assisted reproductive technologies (ART) may perturb DNA methylation (DNAm) in early embryonic development. Although a handful of epigenome-wide association studies of ART have been published, none have investigated CpGs on the X chromosome. To bridge this knowledge gap, we leveraged one of the largest collections of mother-father-newborn trios of ART and non-ART (natural) conceptions to date to investigate sex-specific DNAm differences on the X chromosome. The discovery cohort consisted of 982 ART and 963 non-ART trios from the Norwegian Mother, Father, and Child Cohort Study (MoBa). To verify our results from the MoBa cohort, we used an external cohort of 149 ART and 58 non-ART neonates from the Australian 'Clinical review of the Health of adults conceived following Assisted Reproductive Technologies' (CHART) study. The Illumina EPIC array was used to measure DNAm in both datasets. In the MoBa cohort, we performed a set of X-chromosome-wide association studies ('XWASs' hereafter) to search for sex-specific DNAm differences between ART and non-ART newborns. We tested several models to investigate the influence of various confounders, including parental DNAm. We also searched for differentially methylated regions (DMRs) and regions of co-methylation flanking the most significant CpGs. Additionally, we ran an analogous model to our main model on the external CHART dataset. RESULTS: In the MoBa cohort, we found more differentially methylated CpGs and DMRs in girls than boys. Most of the associations persisted after controlling for parental DNAm and other confounders. Many of the significant CpGs and DMRs were in gene-promoter regions, and several of the genes linked to these CpGs are expressed in tissues relevant for both ART and sex (testis, placenta, and fallopian tube). We found no support for parental DNAm-dependent features as an explanation for the observed associations in the newborns. The most significant CpG in the boys-only analysis was in UBE2DNL, which is expressed in testes but with unknown function. The most significant CpGs in the girls-only analysis were in EIF2S3 and AMOT. These three loci also displayed differential DNAm in the CHART cohort. CONCLUSIONS: Genes that co-localized with the significant CpGs and DMRs associated with ART are implicated in several key biological processes (e.g., neurodevelopment) and disorders (e.g., intellectual disability and autism). These connections are particularly compelling in light of previous findings indicating that neurodevelopmental outcomes differ in ART-conceived children compared to those naturally conceived.

Isolated oral clefts and school grades: population-based cohort study from Norway.

OBJECTIVE: To compare school grades of adolescents in Norway born with isolated cleft with those of their unaffected peers. DESIGN: Population-based cohort study. SETTING: Norway. PATIENTS: A total of 347 419 individuals born in Norway between 1986 and 1992, including 523 isolated cleft cases which were identified using data from Norway's two treatment centres. Individuals were followed from birth through compulsory school. MAIN OUTCOME MEASURES: Grade point average (GPA) from middle school graduation (around the age of 16). Specific subject grades were also investigated. RESULTS: Using a grade scale from 1-6, the observed mean GPA for the reference group was 3.99. Both cleft lip only (CLO) and cleft lip with cleft palate (CLP) had a mean GPA similar to the reference group (adjusted GPA differences from the reference with 95% CIs of 0.06 (-0.04 to 0.16) and -0.08 (-0.19 to 0.03), respectively). Cleft palate only (CPO) had a marginally lower GPA (adjusted GPA difference: -0.18 (-0.28 to -0.08)). These comparisons were consistent across specific subjects. Overall, the evidence suggests a larger difference in GPA between cases and controls in males compared with females. Females with CLO even had a higher estimated GPA than females in the reference group (adjusted GPA difference: 0.19 (0.013 to 0.36)). Grades were similar regardless of laterality of cleft lip (CLO or CLP). CONCLUSION: In Norway, individuals born with isolated CLO or CLP did not have lower average school grades when graduating from middle school. Individuals born with isolated CPO had marginally lower grades.

Acute exacerbations of COPD and risk of lung cancer in COPD patients with and without a history of asthma.

RATIONALE: There is limited knowledge on the effect of acute exacerbations in chronic obstructive pulmonary disease (AECOPD) on lung cancer risk in COPD patients with and without a history of asthma. This study aims to examine whether AECOPD is associated with risk of lung cancer, and whether the effect depends on a history of asthma. METHODS: In the GenKOLS study of 2003-2005, 852 subjects with COPD performed spirometry, and filled out questionnaires on smoking habits, symptoms and disease history. These data were linked to lung cancer data from the Cancer Registry of Norway through 2013. AECOPD, measured at baseline was the main predictor. To quantify differences in lung cancer risk, we performed Cox-proportional hazards regression. We adjusted for sex, age, smoking variables, body mass index, and lung function. MEASUREMENTS AND RESULTS: During follow-up, 8.8% of the subjects with, and 5.9% of the subjects without exacerbations were diagnosed with lung cancer. Cox regression showed a significant increased risk of lung cancer with one or more exacerbations in COPD patients without a history of asthma, HRR = 2.77 (95% CI 1.39-5.52). We found a significant interaction between a history of asthma and AECOPD on lung cancer. CONCLUSIONS: AECOPD is associated with an increased risk of lung cancer in COPD patients without a history of asthma.

Gene-methylation interactions: discovering region-wise DNA methylation levels that modify SNP-associated disease risk.

BACKGROUND: Current technology allows rapid assessment of DNA sequences and methylation levels at a single-site resolution for hundreds of thousands of sites in the human genome, in thousands of individuals simultaneously. This has led to an increase in epigenome-wide association studies (EWAS) of complex traits, particularly those that are poorly explained by previous genome-wide association studies (GWAS). However, the genome and epigenome are intertwined, e.g., DNA methylation is known to affect gene expression through, for example, genomic imprinting. There is thus a need to go beyond single-omics data analyses and develop interaction models that allow a meaningful combination of information from EWAS and GWAS. RESULTS: We present two new methods for genetic association analyses that treat offspring DNA methylation levels as environmental exposure. Our approach searches for statistical interactions between SNP alleles and DNA methylation (G ×Me) and between parent-of-origin effects and DNA methylation (PoO ×Me), using case-parent triads or dyads. We use summarized methylation levels over nearby genomic region to ease biological interpretation. The methods were tested on a dataset of parent-offspring dyads, with EWAS data on the offspring. Our results showed that methylation levels around a SNP can significantly alter the estimated relative risk. Moreover, we show how a control dataset can identify false positives. CONCLUSIONS: The new methods, G ×Me and PoO ×Me, integrate DNA methylation in the assessment of genetic relative risks and thus enable a more comprehensive biological interpretation of genome-wide scans. Moreover, our strategy of condensing DNA methylation levels within regions helps overcome specific disadvantages of using sparse chip-based measurements. The methods are implemented in the freely available R package Haplin ( https://cran.r-project.org/package=Haplin ), enabling fast scans of multi-omics datasets.

Design efficiency in genetic association studies.

Selecting the best design for genetic association studies requires careful deliberation; different study designs can be used to scan for different genetic effects, and each design has its own set of strengths and limitations. A variety of family and unrelated control configurations are amenable to genetic association analyses, including the case-control design, case-parent triads, and case-parent triads in combination with unrelated controls or control-parent triads. Ultimately, the goal is to choose the design that achieves the highest statistical power using the lowest cost. For given parameter values and genotyped individuals, designs can be compared directly by computing the power. However, a more informative and general design comparison can be achieved by studying the relative efficiency, defined as the ratio of variances of two different parameter estimators, corresponding to two separate designs. Using log-linear modeling, we derive the relative efficiency from the asymptotic variance of the parameter estimators and relate it to the concept of Pitman efficiency. The relative efficiency takes into account the fact that different designs impose different costs relative to the number of genotyped individuals. We show that while optimal efficiency for analyses of regular autosomal effects is achieved using the standard case-control design, the case-parent triad design without unrelated controls is efficient when searching for parent-of-origin effects. Due to the potential loss of efficiency, maternal genes should generally not be adjusted for in an initial genome-wide association study scan of offspring genes but instead checked post hoc. The relative efficiency calculations are implemented in our R package Haplin.

A genome-wide scan of cleft lip triads identifies parent-of-origin interaction effects between ANK3 and maternal smoking, and between ARHGEF10 and alcohol consumption.

Background: Although both genetic and environmental factors have been reported to influence the risk of isolated cleft lip with or without cleft palate (CL/P), the exact mechanisms behind CL/P are still largely unaccounted for. We recently developed new methods to identify parent-of-origin (PoO) interactions with environmental exposures (PoOxE) and applied them to families with children born with isolated cleft palate only. Here, we used the same genome-wide association study (GWAS) dataset and methodology to screen for PoOxE effects in the larger sample of CL/P triads. Methods: Genotypes from 1594 complete triads and 314 dyads (1908 nuclear families in total) with CL/P were available for the current analyses. Of these families, 1024 were Asian, 825 were European and 59 had other ancestries. After quality control, 341,191 SNPs remained from the original 569,244. The exposures were maternal cigarette smoking, use of alcohol, and use of vitamin supplements in the periconceptional period. The methodology applied in the analyses is implemented in the R-package Haplin. Results: Among Europeans, there was evidence of a PoOxSmoke effect for ANK3 with three SNPs (rs3793861, q=0.20, p=2.6e-6; rs7087489, q=0.20, p=3.1e-6; rs4310561, q=0.67, p=4.0e-5) and a PoOxAlcohol effect for ARHGEF10 with two SNPs (rs2294035, q=0.32, p=2.9e-6; rs4876274, q=0.76, p=1.3e-5). Conclusion: Our results indicate that the detected PoOxE effects have a plausible biological basis, and thus warrant replication in other independent cleft samples. Our demonstration of the feasibility of identifying complex interactions between relevant environmental exposures and PoO effects offers new avenues for future research aimed at unravelling  the complex etiology of cleft lip defects.

Haplin power analysis: a software module for power and sample size calculations in genetic association analyses of family triads and unrelated controls.

BACKGROUND: Log-linear and multinomial modeling offer a flexible framework for genetic association analyses of offspring (child), parent-of-origin and maternal effects, based on genotype data from a variety of child-parent configurations. Although the calculation of statistical power or sample size is an important first step in the planning of any scientific study, there is currently a lack of software for genetic power calculations in family-based study designs. Here, we address this shortcoming through new implementations of power calculations in the R package Haplin, which is a flexible and robust software for genetic epidemiological analyses. Power calculations in Haplin can be performed analytically using the asymptotic variance-covariance structure of the parameter estimator, or else by a straightforward simulation approach. Haplin performs power calculations for child, parent-of-origin and maternal effects, as well as for gene-environment interactions. The power can be calculated for both single SNPs and haplotypes, either autosomal or X-linked. Moreover, Haplin enables power calculations for different child-parent configurations, including (but not limited to) case-parent triads, case-mother dyads, and case-parent triads in combination with unrelated control-parent triads. RESULTS: We compared the asymptotic power approximations to the power of analysis attained with Haplin. For external validation, the results were further compared to the power of analysis attained by the EMIM software using data simulations from Haplin. Consistency observed between Haplin and EMIM across various genetic scenarios confirms the computational accuracy of the inference methods used in both programs. The results also demonstrate that power calculations in Haplin are applicable to genetic association studies using either log-linear or multinomial modeling approaches. CONCLUSIONS: Haplin provides a robust and reliable framework for power calculations in genetic association analyses for a wide range of genetic effects and etiologic scenarios, based on genotype data from a variety of child-parent configurations.

Comparison of two lung cancer screening scores among patients with chronic obstructive pulmonary disease: A community study.

INTRODUCTION: Based on the National Lung Cancer Screening Trial (NLST), guidelines on screening programs for lung cancer have recommended low-dose computed tomography (LDCT). De Torres et al made a score for COPD patients (COPD-LUCSS) to improve their selection criteria. OBJECTIVE: To examine and compare the discriminating value of both scores in a community-based cohort of COPD patients. METHODS: Four hundred and twenty-two ever-smokers with COPD from the GenKOLS study in Bergen were merged with the Cancer Registry of Norway. We divided the patients into groups of high and low risk according to the COPD-LUCSS and the NLST criteria. Cox regression and logistic regression were used to analyse the associations between the scores and lung cancer. We used Harrell's C and area under the curve (AUC) to estimate discriminating values and to compare the models. RESULTS: Hazard ratio for the high risk vs the low risk in the COPD-LUCSS was 3.0 (1.4-6.5 95% CI), P < 0.01. Hazard ratio for the NLST criteria was 2.2 (95% CI 1.1-4.5), P < 0.05. Harrell's C was 0.63 for the COPD-LUCSS and 0.59 for the NLST selection criteria. AUC was 0.61 for COPD-LUCSS and 0.59 for NLST criteria. Comparing tests showed no differences (P = 0.76). CONCLUSION: Although the COPD-LUCSS and the NLST criteria were associated with increased risk of lung cancer, the AUC and Harrell's C values showed that these models have poor discriminating abilities in our cohort of COPD patients. The COPD-LUCSS was not significantly better than the NLST criteria.

A Genome-Wide Search for Gene-Environment Effects in Isolated Cleft Lip with or without Cleft Palate Triads Points to an Interaction between Maternal Periconceptional Vitamin Use and Variants in ESRRG.

Background: It is widely accepted that cleft lip with or without cleft palate (CL/P) results from the complex interplay between multiple genetic and environmental factors. However, a robust investigation of these gene-environment (GxE) interactions at a genome-wide level is still lacking for isolated CL/P. Materials and Methods: We used our R-package Haplin to perform a genome-wide search for GxE effects in isolated CL/P. From a previously published GWAS, genotypes and information on maternal periconceptional cigarette smoking, alcohol intake, and vitamin use were available on 1908 isolated CL/P triads of predominantly European or Asian ancestry. A GxE effect is present if the relative risk estimates for gene-effects in the offspring are different across exposure strata. We tested this using the relative risk ratio (RRR). Besides analyzing all ethnicities combined ("pooled analysis"), separate analyses were conducted on Europeans and Asians to investigate ethnicity-specific effects. To control for multiple testing, q-values were calculated from the p-values. Results: We identified significant GxVitamin interactions with three SNPs in "Estrogen-related receptor gamma" (ESRRG) in the pooled analysis. The RRRs (95% confidence intervals) were 0.56 (0.45-0.69) with rs1339221 (q = 0.011), 0.57 (0.46-0.70) with rs11117745 (q = 0.011), and 0.62 (0.50-0.76) with rs2099557 (q = 0.037). The associations were stronger when these SNPs were analyzed as haplotypes composed of two-SNP and three-SNP combinations. The strongest effect was with the "t-t-t" haplotype of the rs1339221-rs11117745-rs2099557 combination [RRR = 0.50 (0.40-0.64)], suggesting that the effects observed with the other SNP combinations, including those in the single-SNP analyses, were mainly driven by this haplotype. Although there were potential GxVitamin effects with rs17734557 and rs1316471 and GxAlcohol effects with rs9653456 and rs921876 in the European sample, respectively, none of the SNPs was located in or near genes with strong links to orofacial clefts. GxAlcohol and GxSmoke effects were not assessed in the Asian sample because of a lack of observations for these exposures. Discussion/Conclusion: We identified significant interactions between vitamin use and variants in ESRRG in the pooled analysis. These GxE effects are novel and warrant further investigations to elucidate their roles in orofacial clefting. If validated, they could provide prospects for exploring the impact of estrogens and vitamins on clefting, with potential translational applications.

Analysis of Parent-of-Origin Effects on the X Chromosome in Asian and European Orofacial Cleft Triads Identifies Associations with DMD, FGF13, EGFL6, and Additional Loci at Xp22.2.

Background: Although both the mother's and father's alleles are present in the offspring, they may not operate at the same level. These parent-of-origin (PoO) effects have not yet been explored on the X chromosome, which motivated us to develop new methods for detecting such effects. Orofacial clefts (OFCs) exhibit sex-specific differences in prevalence and are examples of traits where a search for various types of effects on the X chromosome might be relevant. Materials and Methods: We upgraded our R-package Haplin to enable genome-wide analyses of PoO effects, as well as power simulations for different statistical models. 14,486 X-chromosome SNPs in 1,291 Asian and 1,118 European case-parent triads of isolated OFCs were available from a previous GWAS. For each ethnicity, cleft lip with or without cleft palate (CL/P) and cleft palate only (CPO) were analyzed separately using two X-inactivation models and a sliding-window approach to haplotype analysis. In addition, we performed analyses restricted to female offspring. Results: Associations were identified in "Dystrophin" (DMD, Xp21.2-p21.1), "Fibroblast growth factor 13" (FGF13, Xq26.3-q27.1) and "EGF-like domain multiple 6" (EGFL6, Xp22.2), with biologically plausible links to OFCs. Unlike EGFL6, the other associations on chromosomal region Xp22.2 had no apparent connections to OFCs. However, the Xp22.2 region itself is of potential interest because it contains genes for clefting syndromes [for example, "Oral-facial-digital syndrome 1" (OFD1) and "Midline 1" (MID1)]. Overall, the identified associations were highly specific for ethnicity, cleft subtype and X-inactivation model, except for DMD in which associations were identified in both CPO and CL/P, in the model with X-inactivation and in Europeans only. Discussion/Conclusion: The specificity of the associations for ethnicity, cleft subtype and X-inactivation model underscores the utility of conducting subanalyses, despite the ensuing need to adjust for additional multiple testing. Further investigations are needed to confirm the associations with DMD, EGF16, and FGF13. Furthermore, chromosomal region Xp22.2 appears to be a hotspot for genes implicated in clefting syndromes and thus constitutes an exciting direction to pursue in future OFCs research. More generally, the new methods presented here are readily adaptable to the study of X-linked PoO effects in other outcomes that use a family-based design.

Parent-of-origin-environment interactions in case-parent triads with or without independent controls.

With case-parent triad data, one can frequently deduce parent of origin of the child's alleles. This allows a parent-of-origin (PoO) effect to be estimated as the ratio of relative risks associated with the alleles inherited from the mother and the father, respectively. A possible cause of PoO effects is DNA methylation, leading to genomic imprinting. Because environmental exposures may influence methylation patterns, gene-environment interaction studies should be extended to allow for interactions between PoO effects and environmental exposures (i.e., PoOxE). One should thus search for loci where the environmental exposure modifies the PoO effect. We have developed an extensive framework to analyze PoOxE effects in genome-wide association studies (GWAS), based on complete or incomplete case-parent triads with or without independent control triads. The interaction approach is based on analyzing triads in each exposure stratum using maximum likelihood estimation in a log-linear model. Interactions are then tested applying a Wald-based posttest of parameters across strata. Our framework includes a complete setup for power calculations. We have implemented the models in the R software package Haplin. To illustrate our PoOxE test, we applied the new methodology to top hits from our previous GWAS, assessing whether smoking during the periconceptional period modifies PoO effects on cleft palate only.

A new approach to chromosome-wide analysis of X-linked markers identifies new associations in Asian and European case-parent triads of orofacial clefts.

BACKGROUND: GWAS discoveries on the X-chromosome are underrepresented in the literature primarily because the analytical tools that have been applied were originally designed for autosomal markers. Our objective here is to employ a new robust and flexible tool for chromosome-wide analysis of X-linked markers in complex traits. Orofacial clefts are good candidates for such analysis because of the consistently observed excess of females with cleft palate only (CPO) and excess of males with cleft lip with or without cleft palate (CL/P). METHODS: Genotypes for 14,486 X-chromosome SNPs in 1,291 Asian and 1,118 European isolated cleft triads were available from a previously published GWAS. The R-package HAPLIN enables genome-wide-level analyses as well as statistical power simulations for a range of biologic scenarios. We analyzed isolated CL/P and isolated CPO for each ethnicity in HAPLIN, using a sliding-window approach to haplotype analysis and two different statistical models, with and without X-inactivation in females. RESULTS: There was a larger number of associations in the Asian versus the European sample, and similar to previous reports that have analyzed the same GWAS dataset using different methods, we identified associations with EFNB1/PJA1 and DMD. In addition, new associations were detected with several other genes, among which KLHL4, TBX22, CPXCR1 and BCOR were noteworthy because of their roles in clefting syndromes. A few of the associations were only detected by one particular X-inactivation model, whereas a few others were only detected in one sex. DISCUSSION/CONCLUSION: We found new support for the involvement of X-linked variants in isolated clefts. The associations were specific for ethnicity, sex and model parameterization, highlighting the need for flexible tools that are capable of detecting and estimating such effects. Further efforts are needed to verify and elucidate the potential roles of EFNB1/PJA1, KLHL4, TBX22, CPXCR1 and BCOR in isolated clefts.

Genome-wide analysis of parent-of-origin interaction effects with environmental exposure (PoOxE): An application to European and Asian cleft palate trios.

Cleft palate only is a common birth defect with high heritability. Only a small fraction of this heritability is explained by the genetic variants identified so far, underscoring the need to investigate other disease mechanisms, such as gene-environment (GxE) interactions and parent-of-origin (PoO) effects. Furthermore, PoO effects may vary across exposure levels (PoOxE effects). Such variation is the focus of this study. We upgraded the R-package Haplin to enable direct tests of PoOxE effects at the genome-wide level. From a previous GWAS, we had genotypes for 550 case-parent trios, of mainly European and Asian ancestry, and data on three maternal exposures (smoking, alcohol, and vitamins). Data were analyzed for Europeans and Asians separately, and also for all ethnicities combined. To account for multiple testing, a false discovery rate method was used, where q-values were generated from the p-values. In the Europeans-only analyses, interactions with maternal smoking yielded the lowest q-values. Two SNPs in the 'Interactor of little elongation complex ELL subunit 1' (ICE1) gene had a q-value of 0.14, and five of the 20 most significant SNPs were in the 'N-acetylated alpha-linked acidic dipeptidase-like 2' (NAALADL2) gene. No evidence of PoOxE effects was found in the other analyses. The connections to ICE1 and NAALADL2 are novel and warrant further investigation. More generally, the new methodology presented here is easily applicable to other traits and exposures in which a family-based study design has been implemented.

Incidence of non-pulmonary cancer and lung cancer by amount of emphysema and airway wall thickness: a community-based cohort.

There is limited knowledge about the prognostic value of quantitative computed tomography (CT) measures of emphysema and airway wall thickness in cancer.The aim of this study was to investigate if using CT to quantitatively assess the amount of emphysema and airway wall thickness independently predicts the subsequent incidence of non-pulmonary cancer and lung cancer.In the GenKOLS study of 2003-2005, 947 ever-smokers performed spirometry and underwent CT examination. The main predictors were the amount of emphysema measured by the percentage of low attenuation areas (%LAA) on CT and standardised measures of airway wall thickness (AWT-PI10). Cancer data from 2003-2013 were obtained from the Norwegian Cancer Register. The hazard ratio associated with emphysema and airway wall thickness was assessed using Cox proportional hazards regression for cancer diagnoses.During 10 years of follow-up, non-pulmonary cancer was diagnosed in 11% of the subjects with LAA <3%, in 19% of subjects with LAA 3-10%, and in 17% of subjects with LAA ≥10%. Corresponding numbers for lung cancer were 2%, 3% and 11%, respectively. After adjustment, the baseline amount of emphysema remained a significant predictor of the incidence of non-pulmonary cancer and lung cancer. Airway wall thickness did not predict cancer independently.This study offers a strong argument that emphysema is an independent risk factor for both non-pulmonary cancer and lung cancer.

Pneumothorax size measurements on digital chest radiographs: Intra- and inter- rater reliability.

PURPOSE: Detailed and reliable methods may be important for discussions on the importance of pneumothorax size in clinical decision-making. Rhea's method is widely used to estimate pneumothorax size in percent based on chest X-rays (CXRs) from three measure points. Choi's addendum is used for anterioposterior projections. The aim of this study was to examine the intrarater and interrater reliability of the Rhea and Choi method using digital CXR in the ward based PACS monitors. MATERIALS AND METHODS: Three physicians examined a retrospective series of 80 digital CXRs showing pneumothorax, using Rhea and Choi's method, then repeated in a random order two weeks later. We used the analysis of variance technique by Eliasziw et al. to assess the intrarater and interrater reliability in altogether 480 estimations of pneumothorax size. RESULTS: Estimated pneumothorax sizes ranged between 5% and 100%. The intrarater reliability coefficient was 0.98 (95% one-sided lower-limit confidence interval C 0.96), and the interrater reliability coefficient was 0.95 (95% one-sided lower-limit confidence interval 0.93). CONCLUSION: This study has shown that the Rhea and Choi method for calculating pneumothorax size has high intrarater and interrater reliability. These results are valid across gender, side of pneumothorax and whether the patient is diagnosed with primary or secondary pneumothorax.

The Relationship of Educational Attainment with Pulmonary Emphysema and Airway Wall Thickness.

RATIONALE: Low educational attainment is a risk factor of chronic obstructive pulmonary disease (COPD). There is limited knowledge on the relationship between educational level and computed tomography measures of emphysema and airway wall thickness (AWT). OBJECTIVES: We hypothesized that low educational attainment is associated with increased emphysema and AWT in ever-smokers with and without COPD. METHODS: We included 462 and 485 ever-smokers with and without COPD in a cross-sectional study, aged 40-86 years. The sample was divided into groups reflecting educational attainment: primary, secondary, and university. We performed linear regression to examine associations between educational attainment and both emphysema and AWT separately for those with and without COPD. We adjusted for sex, age, smoking status, age of onset of smoking, pack-years, height, and body mass index. MEASUREMENTS AND MAIN RESULTS: Compared with university education, in subjects with COPD, primary education was associated with a 68.1% (95% confidence interval = 14.2-147.6%; P = 0.01) relative increase in emphysema and secondary education was associated with a 50.6% (95% confidence interval = 5.7-114.6%; P = 0.02) relative increase. There was a nonsignificant trend toward an association between lower educational attainment and increased emphysema among those without COPD (P = 0.18), yet greater age appeared to modify this association (P = 0.01). We did not detect significant linear relationships between educational attainment and AWT in subjects with or without COPD. CONCLUSIONS: Lower educational attainment was associated with increased emphysema among adults with COPD. Among those without COPD, this association was more pronounced with increasing age. No significant linear relationship between educational attainment and AWT was found. Clinicians treating adults with emphysema should keep in mind that factors related to low education beyond that of smoking and occupational dust exposure might be of importance to the disease.

Improving the error rates of the Begg and Mazumdar test for publication bias in fixed effects meta-analysis.

BACKGROUND: The rank correlation test introduced by Begg and Mazumdar is extensively used in meta-analysis to test for publication bias in clinical and epidemiological studies. It is based on correlating the standardized treatment effect with the variance of the treatment effect using Kendall's tau as the measure of association. To our knowledge, the operational characteristics regarding the significance level of the test have not, however, been fully assessed. METHODS: We propose an alternative rank correlation test to improve the error rates of the original Begg and Mazumdar test. This test is based on the simulated distribution of the estimated measure of association, conditional on sampling variances. Furthermore, Spearman's rho is suggested as an alternative rank correlation coefficient. The attained level and power of the tests are studied by simulations of meta-analyses assuming the fixed effects model. RESULTS: The significance levels of the original Begg and Mazumdar test often deviate considerably from the nominal level, the null hypothesis being rejected too infrequently. It is proven mathematically that the assumptions for using the rank correlation test are not strictly satisfied. The pairs of variables fail to be independent, and there is a correlation between the standardized effect sizes and sampling variances under the null hypothesis of no publication bias. In the meta-analysis setting, the adverse consequences of a false negative test are more profound than the disadvantages of a false positive test. Our alternative test improves the error rates in fixed effects meta-analysis. Its significance level equals the nominal value, and the Type II error rate is reduced. In small data sets Spearman's rho should be preferred to Kendall's tau as the measure of association. CONCLUSIONS: As the attained significance levels of the test introduced by Begg and Mazumdar often deviate greatly from the nominal level, modified rank correlation tests, improving the error rates, should be preferred when testing for publication bias assuming fixed effects meta-analysis.