PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Cancer Epidemiol Biomarkers Prev. Author manuscript; available in PMC 2011 March 3.
Published in final edited form as:
PMCID: PMC2838211
NIHMSID: NIHMS167989

A Segregation Analysis of Barrett’s Esophagus and Associated Adenocarcinomas

Abstract

Familial aggregation of esophageal adenocarcinomas, esophagogastric junction adenocarcinomas, and their precursor Barrett’s esophagus has been termed Familial Barrett’s Esophagus (FBE). Numerous studies documenting increased familial risk for these diseases raise the hypothesis that there may be an inherited susceptibility to the development of BE and its associated cancers. In this study, using segregation analysis for a binary trait as implemented in S.A.G.E. 6.0.1, we analyzed data on 881singly ascertained pedigrees in order to determine whether FBE is caused by a common environmental or genetic agent and, if genetic, to identify the mode of inheritance of FBE. The inheritance models were compared by likelihood ratio tests and Akaike’s A Information Criterion. Results indicated that random environmental and/or multifactorial components were insufficient to fully explain the familial nature of FBE, but rather there is segregation of a major type transmitted from one generation to the next (p-value < 10−10). An incompletely dominant inheritance model together with a polygenic component fits the data best. For this dominant model, the estimated penetrance of the dominant allele is 0.1005 (95% confidence interval, CI: 0.0587 to 0.1667) and the sporadic rate is 0.0012 (95% CI: 0.0004 to 0.0042), corresponding to a relative risk of 82.53 (95% CI: 28.70 to 237.35), or odds ratio of 91.63 (95% CI: 32.01 to 262.29). This segregation analysis provides epidemiological evidence in support of one or more rare autosomally inherited dominant susceptibility allele(s) in FBE families, and hence motivates linkage analyses.

Keywords: familial esophageal adenocarcinomas, complex segregation analysis, dominant major gene inheritance, polygenic component, likelihood, AIC, unified model

Introduction

The striking increase in the incidence of Barrett’s esophagus (BE), esophageal adenocarcinoma (EAC) and esophagogastric junctional adenocarcinoma (EGJAC) over the last few decades implicates in the pathogenesis of these diseasesa potential major environmental factor, such as gastroesophageal reflux (GERD), obesity, and/or smoking. Yet, in the past two decades there have been numerous case reports, case series, and cross sectional studies documenting increased familial risk for these diseases (18) raising the hypothesis that there may be an inherited susceptibility to the development of BE and its associated cancers. Indeed, there is strong evidence of a genetic component in the pathogenesis of GERD (911) and obesity (1214). Just as a change in environmental factors within the past several decades has interacted with undiscovered genetic factors to contribute to the increasing epidemic of obesity in the United States, it is quite plausible that changes in environmental factors are acting together with genetic susceptibility factors to contribute to the rising incidence of BE, EAC and EGJAC. Of course, the reported increased familial risk of BE and its associated cancers could be related to either shared environmental factors within families or to transmitted genetic susceptibility factors. The purpose of this study was to perform a complex segregation analysis on pedigrees of probands with BE, EAC, and EGJAC to determine whether the familial aggregation of these diseases is consistent with the transmission of autosomal genes. An incompletely dominant inheritance model together with a polygenic component was identified. Such a finding on a large sample would provide strong motivation for a linkage analysis of multiplex families, as well as provide parameter estimates that would be expected to maximize the power of a model-based linkage analysis.

Materials and Methods

Recruitment of Probands

The familial Barrett’s esophagus (FBE) pedigrees were collected as part of an active multi-center study registered at clinicaltrials.gov, identifier NCT00288119, where the ultimate goal of this accrual is to identify putative susceptibility genes. The methodology used in this study to recruit probands and develop the family structures to identify FBE pedigrees has been reported before (8, 15). Briefly, families were ascertained through probands with established or newly diagnosed BE, EAC, or GEJAC at 5 tertiary care academic hospitals in the United States as previously described (15). Recruitment periods varied among the hospitals, depending on available personnel and institutional review board approval, and ranged from 6 months to 7 years in duration. All eligible patients seen in the endoscopy suite during the active recruitment period were approached for study entry.

Determination of Family History

A FBE questionnaire that elicits information about reflux symptoms (patterned on the validated Mayo GERQ)(16), relevant covariates and a detailed family history of BE, esophageal cancer, and other cancers, was used in this ongoing study. Probands were defined as the affected member who brought the family to the study and each family was assumed to be singly ascertained via the sole proband in each family. A FBE questionnaire was administered to these probands by a study nurse or sent them via mail. Mailed questionnaires were followed up with a phone call by the study nurse. Probands were requested to obtain permission for study personnel to contact their first degree family members, and similarly for their affected relatives. Questionnaires were also administered to relatives who contacted the study centers and consented to participate. For deceased relatives, the next of kin was asked to complete relevant portions of the questionnaire and provide consent for review of the medical records. Screening endoscopy was offered to family members who had not had previous upper endoscopy. The endoscopic and histological diagnosis for family members reported to have BE or esophageal cancer was confirmed by reviewing endoscopy and pathology records. For all reported diagnoses that could not be confirmed by review of medical records, the diagnosis of BE or EAC was classified as possible but not confirmed. Each institutional review board for human investigation at participating hospitals approved the study protocol.

Definition of Trait and Other Variables

EAC was defined as adenocarcinoma from a mass that predominantly involves the tubular esophagus, and EGJAC was defined as adenocarcinoma on histology from a mass that predominantly involves the gastro-esophageal junction but does have some esophageal involvement. The definition of BE required clear and rigorous documentation of measurements of the affected segment of tubular esophagus in the endoscopy report and documented presence of intestinal metaplasia on biopsy. Biopsies from the gastroesophageal junction or an irregular Z-line containing intestinal metaplasia were classified as intestinal metaplasia of the cardia and not considered to be BE. Institutional biopsies and outside biopsies, when obtained, were reviewed by a gastrointestinal pathologist. BE, EAC, and EGJAC have been considered part of a single complex trait termed FBE because these conditions are epidemiologically similar and there is strong evidence that BE is the precursor of nearly all EACs and a substantial proportion of EGJACs (1720). The disease trait for the analysis performed here was defined as BE, EAC, or EGJAC. For the purpose of this segregation analysis, an individual with a confirmed diagnosis of BE, EAC, or EGJAC was considered affected. Individuals with no reported history, or individuals who did not have BE or adenocarcinoma upon review of medical records, were considered unaffected. The trait status for individuals with a diagnosis that could not be confirmed was treated as unknown.

Statistical Analysis

Complex segregation analysis of FBE was performed as implemented in the program SEGREG within the S.A.G.E. 6.0.1 program package(21). The models implemented in SEGREG are briefly described here, with special emphasis on the details that are specific to modeling complex segregation of a binary trait.

The segregation model for a binary trait assumes that susceptibility to disease, γ, defined as the probability that an individual is affected with disease, depends on an unobserved latent factor termed type, designated as u, which can take on one of the three values AA, AB, or BB. If the segregation is Mendelian (in which case γ is interpreted as a penetrance function), the type u represents a putative genotype that underlies the distribution of the observed phenotype. In the transmission probability model, types are characterized by two parameter sets, type frequencies and transmission parameters (22). If one assumes Hardy-Weinberg equilibrium (HWE), the type frequencies in the population can be defined by a single parameter qA, the frequency of A. The set of transmission parameters, τAA, τAB, and τBB, represent the probabilities respectively that individuals of types AA, AB, and BB will transmit the component A (allele, if the type is a genotype) to offspring. Assuming that mating is random, the transmissions from each parent are independent. For a Mendelian locus, HWE is assumed, and τAA, τAB, and τBB are equal to 1, 0.5, and 0, respectively. In the specific case that the disease segregation is caused purely by a random environmental factor, there is no transmission from generation to generation, and the three unobserved types transmit equally (τAA = τAB = τBB). In addition, a “general” transmission model allows the transmission probabilities τAA, τAB and τBB to take on any arbitrary values between 0 and 1. A more restricted general transmission model assumes homogeneity of the phenotypic distribution across generations and must satisfy two conditions: the type frequencies must follow HWE proportions and τAB must be equal to a specific function of the frequency of A, τAA and τBB; thus only two of the transmission probabilities can be freely estimated (23).

For a binary trait, the susceptibility γ is defined as the cumulative logistic function,

equation M1

where, conditional on type u, the logit of the susceptibility for the i-th individual, θu(i), can depend on both major type u and mean centered covariate values xi1, xi2, …, xip:

equation M2

where βu is the intercept corresponding to type u and ξ1, ξ2, … &xi;p are the covariate regression coefficients. This is analogous to estimating covariate coefficients in a linear model simultaneously with all the other parameters in the segregation analysis of a quantitative trait (24).

Two ways to allow for multifactorial components in the analysis of a binary trait are implemented in SEGREG. The first is a multivariate logistic model(MLM) that allows the incorporation of first order nuclear family residual association parameters into the logit of susceptibility(25). In this model, the residual association parameters include δFM for father-mother, δFO for father-offspring, δMO for mother-offspring, and δSS for sib-sib residual associations. These association parameters can be converted to residual correlations that depend on the expected value of the logit. In multigenerational pedigrees, the class D assumption (26) is used, i.e. that all sib-sib residual associations are equal but not necessarily due to common parentage alone; nor does it necessarily assume that they are equal to the parent-offspring residual association. The second model is the finite polygenic mixed model (FPMM) (27, 28), which assumes that the logit of susceptibility is influenced by a small number, v, of additive diallelic loci in addition to possible segregation at a single major locus with large effect. In this model the effect of the finite number of polygenic loci is represented by additive polygenic variance equation M3 of the logit, i.e. within nuclear families the sibling correlation is assumed to be the same as the parent-offspring correlation and the correlation between k-th degree relatives is (½)k−1 times that of first degree relatives. The incorporation of a polygenic component into the transmission probability model results in the unified model (29), but in SEGREG it is extended to multigenerational pedigrees by assuming a finite number of polygenic loci.

The segregation models studied can be classified according to the number of underlying susceptibility types in the model: one major susceptibility type, and therefore no segregation; major type models that assume a mixture of two susceptibility types (dominant or recessive, if there is Mendelian segregation); and major type models that assume a mixture of three susceptibility types. In the analysis presented here, because of difficulties in maximizing a general likelihood in the general codominant case (three unrestricted susceptibilities), we restricted our attention to three susceptibilities that are additive on the logistic scale. In each case, additional familial/multifactorial components were investigated. For the one-susceptibility case, both a purely random environmental model (no transmission, τAA = τAB = τBB = qA) without any multifactorial component, and the same model but with familial/multifactorial components (sometimes called a “sporadic model”), were fitted. For the major type models that have two or three susceptibility types, a homogeneous model with no transmission (τAA = τAB = τBB = qA), a Mendelian transmission model (τAA = 1, τAB = 0.5, τBB = 0), and two “general” transmission models were fitted. These general transmission models subsume both the environmental no transmission and the Mendelian transmission models as special cases, and either assume homogeneity of the phenotype distribution across generations (τAB is a function of qAAA and τBB) or they do not. In the latter case, all three transmission probabilities are freely estimated in the interval [0,1] without any further restriction.

SEGREG also allows computing the likelihood of a no transmission model (τAA = τAB = τBB) that is not homogeneous across generations, but results for this model were not reported here because they did not change any of the conclusions. It should be noted that, in contrast to what occurs for a trait that can take on more than two values, for a binary trait the no transmission model, either with or without the assumption of homogeneity across generations (corresponding to τAA = τAB = τBB = qA and τAA = τAB = τBB, respectively), is the same whatever the number of susceptibility types if there is no multifactorial component incorporated into the model. This is so because the phenotypic distribution is then a mixture of Bernoulli random variables, which is simply a Bernoulli random variable.

All parameters were estimated by numerical maximization of the likelihood, with standard errors (SEs) calculated by numerical double differentiation of the log likelihood evaluated at the maximum likelihood estimates of all the model parameters. The 95% confidence interval (CI) of an estimate was calculated as estimate±1.96×SE, and the CIs for susceptibility, γ, relative risk and odds ratio were calculated by the delta method. Various hypotheses about mode of transmission are incorporated into models that impose different restrictions on the unrestricted general model, and cause the likelihood to be smaller than that for the general model. Likelihood ratio tests (LRTs) were used to test the significance of the departure from specified null hypothesis models using the asymptotic properties of the LRT: when the null hypothesis is not on a boundary of the unrestricted model, twice the difference in ln(likelihood) between the two models is asymptotically distributed as χ2, with the number of degrees of freedom equal to the difference in the number of independent parameters estimated. In some cases, where the null hypothesis is on the boundary of the unrestricted model, the asymptotic distribution is a mixture of χ2 distributions (30). In other cases the asymptotic distribution is unknown, and then Akaike’s A Information Criterion (AIC) (31) can be used to select the better model. For a fixed number of susceptibility types, the SEGREG output includes p-values based on the appropriate asymptotic distribution of the likelihood ratio criterion, either a chi-square distribution or a mixture of chi-square distributions, as well as AIC values, which can be compared across results for different numbers of susceptibility types. All tests of statistical significance except for variances were two-sided.

To allow for ascertainment, the likelihood of each pedigree was conditioned on the affection status of all individuals in the proband sampling frame, which for single ascertainment comprise the probands (32).

Results

The Population Sample

A total of 881 pedigrees were included in the study and the average pedigree size was 10.03 with standard deviation 8.47. There were 81 (9.2%) pedigrees with 1 generation, 152 (17.2%) pedigrees with 2 generations, 474 (53.8%) pedigrees with 3 generations and 174 (19.8%) pedigrees with 4 or more generations. Of the total 8835 individuals in the data, 992 (11.2%) were affected with BE, EAC or EGJAC. As shown in Table 1, the percentage of males affected is about four times that of females (17.8% vs. 4.2%; assuming independent observations, equation M4, p-value < 2.2×10−16), and the percentage of affected among non-founders is about five times that among founders (15.4% vs. 3.2%; assuming independent observations, equation M5, p-value < 2.2×10−6). These imbalances of affected individuals between males vs. females and founders vs. non-founders were accounted for in the segregation analysis using (0,1) coded covariates.

Table 1
Affection status by gender and founder status for a total of 8835 individuals from FBE families for segregation analysis

One Susceptibility Type Models

First, the effect of covariates on the probability of being affected was evaluated in a one susceptibility type model (which assumes no segregation), denoted “no transmission”. In addition to sex, founder status was also investigated as a covariate because of the noted heterogeneity of the proportions affected between founders and non-founders.

Results in Table 2 indicate that the estimated regression coefficients ξsex and ξfounder are highly significant (more than four times their standard errors) for any model with sex and/or founder as covariate(s);and also by comparing Akaike’s AIC values, the full model, including both sex and founder as covariates, fits the data best. Similar results were obtained by likelihood ratio tests. So, in the following segregation analyses, sex and founder status were always included as covariates in the logit of FBE susceptibility.

Table 2
Parameter estimate in the one susceptibility type no transmission models, with and without (centered) covariates

Next, the effect of incorporating different familial/multifactorial components in these non-segregating models was investigated. The results of familial associations in the MLM model are shown in Table 3. Because there were insufficient data to estimate the father–mother association (it was frequently not possible to maximize the likelihood), we set it equal to 0 for all these models. Although the AIC suggests that including separate parent-offspring and sib-sib associations (δFO = δMOSS) in the model is best, likelihood ratio tests indicated no significant differences among the three associations: the model with three equal familial associations (δFO = δMO = δSS) was not rejected when compared to the model with three different association coefficients (δFO, δMO, δSS) (p-value = 0.269), or with separate parent-offspring and sib-sib associations (δFO = δMOSS) (p-value = 0.106), which suggests that the model with three equal familial multifactorial components fits the data adequately.

Table 3
Parameter estimate in the one susceptibility type (sporadic)model, incorporating various residual familial associations by MLM (δFM = 0)

The results of the one susceptibility type models incorporating various numbers of additive polygenic loci in the FPMM are shown in Supplementary Table S1. They indicate that the likelihood increased with the number of polygenic loci, but the increase was small, especially after including more than three polygenic loci. Furthermore, the AIC was always much smaller than seen in the MLM models in Table 3, which indicates that including a multifactorial component via the FPMM model fits the data better.

In summary, results of analyses with one susceptibility models indicate that (1) sex and founder status were significant covariates of susceptibility to FBE, (2) Including three equal familial/multifactorial components significantly improved the fitting, (3) the multifactorial familial components can be best summarized by an additive polygenic variance, and three additive loci are sufficient to capture this polygenic variance.

Two Susceptibility Type Models

Our next step is to investigate whether the residual variability in the logit of susceptibility, after allowing for the covariates, is significantly better explained by the incorporation of a further susceptibility type (with unequal transmission probabilities) to decide if there may be a major latent factor – whether environmentally or genetically caused – segregating. The FPMM model incorporating three polygenic loci was used to allow for a multifactorial/polygenic component. The transmission mode was described by the transmission parameters τAA, τAB, τBB, and qA, the frequency of susceptibility allele A. Four modes of transmissions were compared: homogeneous no transmission (τAA = τAB = τBB = qA), Mendelian transmission (τAA = 1, τAB = 0.5, τBB = 0), a homogeneous general transmission model (τAA, τBB), and a similar non-homogeneous general model in which the three transmission probabilities (τAA, τAB, τBB) are freely estimated in the interval [0,1]. If there is a polygenic component incorporated into the model, for two and three susceptibility types we denote the homogeneous no transmission model “environmental plus polygenic transmission” because there is a latent major environmental factor “segregating” in addition to polygenic transmission.

First, homogeneity across generations was tested by comparing the homogeneous general transmission (τAA, τBB)model with the general (τAA, τAB, τBB)model. Then, the homogeneous no transmission model and the Mendelian model were each compared with both the homogeneous general (τAA, τBB) model and the general model (τAA, τAB, τBB), to test first if there is a major type segregating, and then, if so, whether this type could be Mendelian.

For two susceptibility types, a purely random environmental model was rejected on comparing it with either the general (τAA, τAB, τBB) model or the homogeneous general (τAA, τBB) model (Supplementary Table S2A, p-value = 1.12×10−13 and 2.70×10−12 respectively). Similarly the environmental model incorporating a polygenic component was also rejected (Supplementary Table S2B, p-value = 8.76×10−13 and 3.11×10−13, respectively). Moreover, purely polygenic transmission together with random environment, represented by the one susceptibility type no transmission model incorporating three polygenic loci, has a much larger AIC value (1090.49, Supplementary Table S1) than any of the genetic models incorporating three polygenic loci that also incorporate two susceptibility types (Supplementary Table S2B, AIC values <1042). These results indicate that neither random environment, nor polygenic transmission, nor a major random environmental effect together with polygenic transmission, are sufficient to fit the data; rather there is a major type being transmitted in some manner. Furthermore, at the 0.05 level, without a polygenic component (Supplementary Table S2A), homogeneity across generations was rejected on comparing the homogeneous general (τAA, τBB) model with the general (τAA, τAB, τBB) model (p-value = 0.002), and both purely dominant and recessive models were rejected on comparing with the general (τAA, τAB, τBB)model (p-values = 0.003). Only under the erroneous homogeneity assumption could dominant and recessive inheritance not be rejected in a model that did not include a polygenic component (p-value = 0.110 and 0.126 respectively for dominant and recessive inheritance). After including three polygenic loci (Supplementary Table S2B), however, homogeneity across generations was not rejected (p-value = 0.206), and Mendelian dominant inheritance was also not rejected – whether compared to the homogenerous general (τAA, τBB)model or to the general (τAA, τAB, τBB) model (p-values = 0.750 and 0.441, respectively) – whereas recessive inheritance was rejected (p-values = 0.001 and 0.002).

In summary, these results indicate that two susceptibility types with a polygenic component support homogeneity across generations, and Mendelian dominant inheritance with a polygenic component fits the data best.

Three Susceptibility Type Models

Because the data were not sufficiently informative to consider a general codominant model, only three susceptibilities that are additive on the logit scale were studied. Results for three additive-susceptibility types in Supplementary Table S3 indicate that there is a major type transmitted (p-values < 10−12), but heterogeneity across generations exists no matter whether a polygenic component is incorporated into the model or not (p-values = 0.035 and 0.030, respectively). Mendelian additive transmission was not rejected when including a polygenic component in the model (p-values = 0.750 and 0.118, comparing with the homogeneous general and the general model, respectively). These results indicate that homogeneity across generations was not supported by three additive susceptibility types, although Mendelian additive transmission was supported after including a polygenic component. Moreover, the additive inheritance model is not better than the dominant inheritance model on comparing their AIC values (Table 4, 1032.10 vs. 1031.98). Table 4, collated from Supplementary Table S2 and S3, gives detailed results of the major genetic models to be compared. Of particular note is the much better fit of dominant over recessive inheritance once a polygenic component is included in the model.

Table 4
Parameter estimates and model AICs under Mendelian transmission of FBE, without and with a polygenic component in the model

In summary, these results indicate that three additive susceptibility types do not fit better than two susceptibility types, and dominant transmission is sufficient to fit the data.

Disease Prevalence

Assuming the estimates for Mendelian dominant inheritance together with a polygenic component, the mode of inheritance found to best fit the data (Table 4), the logit of the susceptibility for types AA and AB, −2.19, corresponds to a penetrance of γAA = γAB = 0.1005 (95% CI: 0.0587, 0.1667), while the logit of the susceptibility for BB, −6.71, corresponds to a sporadic rate of γBB = 0.0012 (95% CI: 0.0004, 0.0042). Because all covariates were centered, this penetrance and sporadic rate are those at average values of the covariates, i.e. the proportions of females, 48.5%, and founders, 34.2%, in these pedigrees. These results correspond, again for average covariate values, to a relative risk of 82.53 (95% CI: 28.70 to 237.35), or odds ratio of 91.63 (95% CI: 32.01 to 262.29). Incorporating the covariates into the logits, and weighting the penetrance and sporadic rate by the estimated genotype frequencies, the disease prevalence was calculated among males and females, and founders and non-founders; this was done both with and without a polygenic component in the model (Table 5). With a polygenic component, the estimated prevalence was 0.48% in males (95% CI: 0.41%, 0.54%) and 0.14% in females(95% CI: 0.11%, 0.16%); and 0.40% in non-founders (95% CI: 0.34%, 0.45%) and 0.12% in founders(95% CI: 0.10%, 0.14%). Combining the sex and founder status, the estimated prevalence was 0.71% in male non-founders (95% CI: 0.63%, 0.80%), 0.21% in male founders (95% CI: 0.18%, 0.25%), 0.21% in female non-founders (95% CI: 0.18%, 0.24%), and 0.06% in female founders (95% CI: 0.05%, 0.07%). Without a polygenic component, the estimated prevalences were higher, 0.60% in males (95% CI: 0.54%, 0.66%) and 0.19% in females (95% CI: 0.17%, 0.21%), and 0.55% in non-founders (95% CI: 0.49%, 0.61%) and 0.14% in founders (95% CI: 0.12%, 0.15%).

Table 5
The estimatedprevalence (and standard error) of FBE in males and females, founders and non-founders under Mendelian dominant transmission

Discussion

The present segregation study analyzed the factors influencing the phenotype transmission step by step. We first confirmed that sex and founder status influence the susceptibility to FBE significantly when included as covariates in the logit of susceptibility. Familial correlations (representing arbitrary multifactorial components) or a polygenic component were also found to contribute to the susceptibility of FBE, and we found that incorporating the latter was better, it being adequate to allow for three additive polygenic loci in the model. Then various transmission models were fitted for two susceptibility types and three susceptibility types. After including the significant polygenic component, dominant inheritance is seen to be the best model.

This segregation analysis of families of patients with BE, EAC, or EGJAC provides strong epidemiological evidence for the segregation of autosomal dominant inheritance at a Mendelian locus that achieves homogeneity across generations. Including a polygenic component also influences the susceptibility to the development of BE or its associated cancers, but it is not sufficient by itself to explain the familial transmission of the disease. In the dominant inheritance model with a polygenic component, the estimated “genetic” susceptibility to develop FBE is limited, in that the overall penetrance of the dominant genotypes, γAA and γAB, is only 0.1005 (95% CI: 0.0587, 0.1667). However, this penetrance is 82.53(95% CI: 28.70, 237.35) times the sporadic rate (γBB = 0.0012, 95% CI: 0.0004, 0.0042) and, using the estimated allele frequency (0.0018) it is found that this corresponds to an odds ratio of 91.63 (95% CI: 32.01, 262.29). It is this difference in magnitude that is key for a successful linkage study. Although the estimated sporadic rate is very low, it is nevertheless significantly different from 0. By comparing the likelihood of this best model to the dominant model with a sporadic rate γBB fixed at 0 (effectively accomplished by fixing βBB = −500), the p-value for departure of γBB from 0 is 1.11×10−13.

Under our best model, the estimated prevalence of disease is higher in non-founders than in founders and higher in males than in females (0.71% in male non-founders, 0.21% in male founders, 0.21% in female non-founders, 0.06% in female founders, Table 5). There has been an increase in obesity over time, which probably accounts for the change in susceptibility to FBE. Moreover, the incidence of EAC has increased dramatically within the U.S. in the past three decades (33, 34). This increase is partly attributable to an epidemic of obesity (35, 36) and might reflect other unrecognized environmental factors. It is likely that the incidence of BE is also increasing along with the increased incidence of its related cancers. The non-homogeneity in disease prevalence between founders and non-founders shown in this study likely reflects the effect of obesity and other environmental factors that have increased the prevalence of disease each generation. Ideally we would have regressed on year of birth and age to allow for such non-homogeneity, but founder status was used as a surrogate because year of birth and/or age were missing on 42% of the pedigree members in the data. Similarly, obesity and smoking status were not considered as covariates because they are only available on 16% of the subjects. The difference in FBE susceptibility between males and females is consistent with the known male predominance of BE and its related cancers (3739).

Our analysis estimated an overall prevalence of BE and its associated cancers to be 0.48% for men and 0.14% in women. Since the prevalence of EAC or EGJAC is only a small fraction of that of BE, this overall prevalence estimate in our study population is essentially an estimate of the prevalence of BE. The estimate is lower than the 1.6% prevalence of BE reported in a Swedish population study (39) and the 6.8% prevalence of BE reported in an American study of patients undergoing colonoscopy (40). The differences in prevalence estimates are likely related to several factors. The prevalence of BE is known to increase with age (3941). The study by Ronkainen et al. (39) measured prevalence in an adult population, the study by Rex et al. (40) examined an even more selected older population undergoing colonoscopy, whereas our study included individuals of all ages. Furthermore, the prevalence of BE has likely undergone a rapid increase in the past several decades along with the increasing incidence of EAC (34). Thus, our study population, which included founders from several decades ago and included a sizable proportion of individuals who never had endoscopy, is expected to have a lower prevalence of BE than the more comprehensive and more recent studies by Ronkainen et al. (39) and Rex et al. (40). Among those individuals we classified as affected, 66.7% had BE, the remainder had EAC or EGJAC. On classifying as unknown those individuals with confirmed EAC or EGJAC, as well as those with a diagnosis that could not be confirmed, the parameter estimates of the dominant model were similar, with (as expected) lower penetrance of the dominant genotypes. In addition, the estimated polygenic variance was 11% larger.

In summary, this segregation analysis provides the first epidemiological evidence in support of a genetic etiology for the familial aggregation of BE, EAC, and EGJAC. Furthermore, the pattern of disease in FBE families is consistent with a rare autosomally inherited dominant susceptibility allele. Although the segregation model assumed only one such allele, this kind of analysis of a binary trait cannot detect allelic heterogeneity, nor can it detect locus heterogeneity if the penetrances are similar across loci. Indeed, the significant polygenic component could be nothing more than an indication of locus heterogeneity, rather than of polygenic inheritance, and that its estimate increased on classifying as unknown those with EAC or EGAC suggests that FBE as defined here would be a better phenotype for linkage analysis. Whereas the power of a linkage analysis is not affected by allelic heterogeneity, it can be seriously diminished by locus heterogeneity. However, with an odds ratio of more than 91, it is difficult to believe that no locus exists for which the odds ratio is high enough to make linkage detectable. Studies need to be performed to discover these loci with putative susceptibility alleles.

Supplementary Material

Acknowledgments

The authors appreciate the accrual efforts of Denise Buonocore-Sassano, Anna Haas, and Kasey Orlowski. The authors thank Mike Warfe for securing the data in a web based data repository.

Grant support: This work was supported by the U.S. Public Health Service research grants R01 DK070863 and K24 DK002800 from the National Institute of Diabetes and Digestive and Kidney Diseases, and R37 GM028356 from National Institute of General Medical Sciences; resource grant P41 RR003655 from the National Center for Research Resources; and Cancer Center Support Grant P30 CA043703 from the National Cancer Institute.

References

1. Gelfand MD. Barrett’s esophagus in sexagenerian identical twins. J Clin Gastroenterol. 1983;5:251–3. [PubMed]
2. Crabb DW, Berk MA, Hall TR, Conneally PM, Biegel AA, Lehman GA. Familial gastroesophageal reflux and development of Barrett’s esophagus. Ann Intern Med. 1985;103:52–4. [PubMed]
3. Prior A, Whorwell PJ. Familial Barrett’s esophagus? Hepatogastroenterology. 1986;33:86–7. [PubMed]
4. Jochem VJ, Fuerst PA, Fromkes JJ. Familial Barrett’s esophagus associated with adenocarcinoma. Gastroenterology. 1992;102:1400–2. [PubMed]
5. Eng C, Spechler SJ, Ruben R, Li FP. Familial Barrett’s esophagus and adenocarcinoma of the gastroesophageal junction. Cancer Epidemiol Biomarkers Prev. 1993;2:397–9. [PubMed]
6. Fahmy N, King JF. Barrett’s esophagus: An acquired condition with a genetic predisposition. Am J Gastroenterol. 1993;88:1262–5. [PubMed]
7. Poynton AR, Walsh TN, O’Sullivan G, Hennessy TP. Carcinoma arising in familial Barrett’s esophagus. Am J Gastroenterol. 1996;91:1855–6. [PubMed]
8. Chak A, Lee T, Kinnard MF, et al. Familial aggregation of Barrett’s esophagus, esophageal adenocarcinoma, and esophagogastric junctional adenocarcinoma in Caucasian adults. Gut. 2002;51:323–8. [PMC free article] [PubMed]
9. Cameron AJ, Lagergren J, Henriksson C, Nyren O, Locke GR, 3rd, Pedersen NL. Gastroesophageal reflux disease in monozygotic and dizygotic twins. Gastroenterology. 2002;122:55–9. [PubMed]
10. Mohammed I, Cherkas LF, Riley SA, Spector TD, Trudgill NJ. Genetic influences in gastro-oesophageal reflux disease: a twin study. Gut. 2003;52:1085–9. [PMC free article] [PubMed]
11. Zheng Z, Nordenstedt H, Pedersen NL, Lagergren J, Ye W. Lifestyle factors and risk for symptomatic gastroesophageal reflux in monozygotic twins. Gastroenterology. 2007;132:87–95. [PMC free article] [PubMed]
12. Farooqi S, O’Rahilly S. Genetics of obesity in humans. Endocr Rev. 2006;27:710–8. [PubMed]
13. Farooqi IS. Monogenic human obesity. Front Horm Res. 2008;36:1–11. [PubMed]
14. Ranadive SA, Vaisse C. Lessons from extreme human obesity: monogenic disorders. Endocrinol Metab Clin North Am. 2008;37:733–51. [PubMed]
15. Chak A, Ochs-Balcom H, Falk G, et al. Familiality in Barrett’s esophagus, adenocarcinoma of the esophagus, and adenocarcinoma of the gastroesophageal junction. Cancer Epidemiol Biomarkers Prev. 2006;15:1668–73. [PubMed]
16. Locke GR, Talley NJ, Weaver AL, Zinsmeister AR. A new questionnaire for gastroesophageal reflux disease. Mayo Clin Proc. 1994;69:539–47. [PubMed]
17. Cameron AJ, Lomboy CT, Pera M, Carpenter HA. Adenocarcinoma of the esophagogastric junctionand Barrett’s esophagus. Gastroenterology. 1995;109:1541–6. [PubMed]
18. Haggitt RC, Tryzelaar J, Ellis FH, Colcher H. Adenocarcinoma complicating columnar epithelium-lined (Barrett’s) esophagus. Am J Clin Pathol. 1978;70:1–5. [PubMed]
19. Hirota WK, Loughney TM, Lazas DJ, Maydonovitch CL, Rholl V, Wong RK. Specialized intestinal metaplasia, dysplasia, and cancer of the esophagus and esophagogastric junction: prevalence and clinical data. Gastroenterology. 1999;116:277–85. [PubMed]
20. Ruol A, Parenti A, Zaninotto G, et al. Intestinal metaplasia is the probable common precursor of adenocarcinoma in barrett esophagus and adenocarcinoma of the gastric cardia. Cancer. 2000;88:2520–8. [PubMed]
21. S.A.G.E 6.0.1. Statistical Analysis for Genetic Epidemiology. 2009. http://darwin.cwru.edu/sage/
22. Elston RC, Stewart J. A general model for the genetic analysis of pedigree data. Hum Hered. 1971;21:523–42. [PubMed]
23. Demenais FM, Elston RC. A general transmission probability model for pedigree data. Hum Hered. 1981;31:93–9. [PubMed]
24. Elston RC, Namboodiri KK, Glueck CJ, Fallat R, Tsang R, Leuba V. Study of the genetic transmission of hypocholesterolemia and hypertriglyceridemia in a 195 member kindred. Ann Hum Genet. 1975;39:67–87. [PubMed]
25. Karunaratne PM, Elston RC. A multivariate logistic model (MLM) for analyzing binary family data. Am J Med Genet. 1998;76:428–37. [PubMed]
26. Bonney GE. On the statistical determination of major gene mechanisms in continuous human traits: regressive models. Am J Med Genet. 1984;18:731–49. [PubMed]
27. Fernando RL, Stricker C, Elston RC. The finite polygenic mixed model: An alternative formulation for the mixed model of inheritance. Theor Appl Genet. 1994;88:573–80. [PubMed]
28. Lange K. An Approximate Model of Polygenic Inheritance. Genetics. 1997;147:1423–30. [PubMed]
29. Lalouel JM, Rao DC, Morton NE, Elston RC. A unified model for complex segregation analysis. Am J Hum Genet. 1983;35:816–26. [PubMed]
30. Self S, Liang K. Asymptotic properties of maximum likelihood estimators and likelihood ratio tests under nonstandard conditions. J Am Stat Assoc. 1987;82:605–10.
31. Akaike H. A new lookat the statistical model identification. IEEE Trans Automat Contr. 1974;AC-19:716–3.
32. Elston RC, Sobel E. Sampling considerations in the gathering and analysis of pedigree data. Am J Hum Genet. 1979;31:62–9. [PubMed]
33. Pera M, Trastek VF, Pairolero PC, Cardesa A, Allen MS, Deschamps C. Barrett’s disease: pathophysiology of metaplasia and adenocarcinoma. Ann Thorac Surg. 1993;56:1191–7. [PubMed]
34. Pohl H, Welch HG. The role of overdiagnosis and reclassification in the marked increase of esophageal adenocarcinoma incidence. J Natl Cancer Inst. 2005;97:142–6. [PubMed]
35. Chow WH, Blot WJ, Vaughan TL, et al. Body mass index and risk of adenocarcinomas of the esophagus and gastric cardia. J Natl Cancer Inst. 1998;90:150–5. [PubMed]
36. Engel LS, Chow WH, Vaughan TL, et al. Population attributable risks of esophageal and gastric cancers. J Natl Cancer Inst. 2003;95:1404–13. [PubMed]
37. Drewitz DJ, Sampliner RE, Garewal HS. The incidence of adenocarcinoma in Barrett’s esophagus: a prospective study of 170 patients followed 4.8 years. Am J Gastroenterol. 1997;92:212–5. [PubMed]
38. Avidan B, Sonnenberg A, Schnell TG, Chejfec G, Metz A, Sontag SJ. Hiatal hernia size, Barrett’s length, and severity of acid reflux are all risk factors for esophageal adenocarcinoma. Am J Gastroenterol. 2002;97:1930–6. [PubMed]
39. Ronkainen J, Aro P, Storskrubb T, et al. Prevalence of Barrett’s esophagus in the general population: an endoscopic study. Gastroenterology. 2005;129:1825–31. [PubMed]
40. Rex DK, Cummings OW, Shaw M, et al. Screening for Barrett’s esophagus in colonoscopy patients with and without heartburn. Gastroenterology. 2003;125(6):1670–7. [PubMed]
41. Cameron AJ, Lomboy CT. Barrett’s esophagus: age, prevalence, and extent of columnar epithelium. Gastroenterology. 1992;103:1241–5. [PubMed]