Skip to main content


Birth order modifies the effect of IL13 gene polymorphisms on serum IgE at age 10 and skin prick test at ages 4, 10 and 18: a prospective birth cohort study

Article metrics



Susceptibility to atopy originates from effects of the environment on genes. Birth order has been identified as a risk factor for atopy and evidence for some candidate genes has been accumulated; however no study has yet assessed a birth order-gene interaction.


To investigate the interaction of IL13 polymorphisms with birth order on allergic sensitization at ages 4, 10 and 18 years.


Mother-infant dyads were recruited antenatally and followed prospectively to age 18 years. Questionnaire data (at birth, age 4, 10, 18); skin prick test (SPT) at ages 4, 10, 18; total serum IgE and specific inhalant screen at age 10; and genotyping for IL13 were collected. Three SNPs were selected from IL13: rs20541 (exon 4, nonsynonymous SNP), rs1800925 (promoter region) and rs2066960 (intron 1). Analysis included multivariable log-linear regression analyses using repeated measurements to estimate prevalence ratios (PRs).


Of the 1456 participants, birth order information was available for 83.2% (1212/1456); SPT was performed on 67.4% at age 4, 71.2% at age 10 and 58.0% at age 18. The prevalence of atopy (sensitization to one or more food or aeroallergens) increased from 19.7% at age 4, to 26.7% at 10 and 41.1% at age 18. Repeated measurement analysis indicated interaction between rs20541 and birth order on SPT. The stratified analyses demonstrated that the effect of IL13 on SPT was restricted only to first-born children (p = 0.007; adjusted PR = 1.35; 95%CI = 1.09, 1.69). Similar findings were noted for firstborns regarding elevated total serum IgE at age 10 (p = 0.007; PR = 1.73; 1.16, 2.57) and specific inhalant screen (p = 0.034; PR = 1.48; 1.03, 2.13).


This is the first study to show an interaction between birth order and IL13 polymorphisms on allergic sensitization. Future functional genetic research need to determine whether or not birth order is related to altered expression and methylation of the IL13 gene.


Atopy has been defined as the propensity of an individual to produce IgE in response to allergen and a predisposition to the development of allergic diseases, such as asthma, atopic dermatitis, food allergy or hay fever. It is defined operationally by elevations in serum levels of immunoglobulin E (IgE) reactive with allergens or skin test reactivity to allergens [1]. Thus atopy is strictly linked to IgE production and describes the personal or familial propensity to become sensitized and produce IgE antibodies in response to environmental triggers [2].

The documented increase in childhood asthma and other atopic disorders over the past three decades has necessitated a search for possible underlying mechanisms and mediators [3, 4]. Genetic variation has been documented to play a role; however, when found, gene effects for allergic diseases are typically small and not easily replicable [1]. In addition, recent studies reveal that a change in the genetic pool is insufficient to account for the temporal and spatial correlates of this increase in prevalence, and suggest that modifications of gene expression and function may be more important [5, 6]. Such modifications of gene expression are thought to result from gene-environment interactions occurring before and after birth. Indeed, previous evidence suggests that developmental processes may modify the impact of genetics on atopy development later in life [7, 8]. Thus, current etiological research is focusing on modifiers of gene effects as possible precipitating factors for the increasing prevalence of atopy in childhood.

Previous studies have assessed two related ideas - the hygiene hypothesis and the birth order effect [9]. The hygiene hypothesis suggests that exposure to infections after birth (due to transmission from older siblings or other children), may influence the development of the immune system along a non-allergic (T helper 1; Th1) pathway, leading to a reduced risk of asthma and other allergic diseases. Although still under investigation, the hygiene hypothesis may explain observed associations between family size, birth order, day-care attendance, and the risk of asthma and allergy. According to the hygiene concept, the effect of birth order is through sibling hierarchy, where the younger child is prone to infection from the older sibling and hence is at lower risk of atopy.

Other studies suggest that birth order may act through a different mechanism [9]. In utero sensitization as a possible mechanism of modification of gene effects has been proposed, where birth order, as an indicator of variations in prenatal exposure, acts independent of the number of siblings [10, 11]. A previous report by our group [12] suggested that the birth order effect may result from in utero exposure as indicated by changes in maternal IgE [12]. If a protective effect of birth order is already present at the time of birth of the child, then the effect of hygiene later in childhood may not be as relevant as is currently thought. In support of this explanation, other non-allergic diseases such as preeclampsia are more common in first born children or first pregnancies, when compared to subsequent births or pregnancies, [13] suggesting altered intra uterine immune adaptation as a possible underlying mechanism.

In a recent review, Vercelli [1] suggested that one of the major challenges facing geneticists is to understand how environmental and developmental factors interact with genetic determinants to increase disease susceptibility. Previous studies have assessed the individual effects of birth order [10, 14, 15] and IL13 polymorphisms [1621] on allergic phenotypes. It is noteworthy that none of these reports examined the possible interaction between birth order and IL13 polymorphisms on childhood atopy.

In previous reports by our group and other researchers, birth order has been repeatedly found to be associated with atopic markers. For instance, we found that cord serum IgE was reduced in offspring with higher order, indicating that the sibling effect may have its origin in utero[12], before any hygiene effect in early childhood could affect the offspring. In addition, we found that the number of live offspring significantly reduces maternal IgE [22], which supports the idea that maternal immune tolerance against allergens may increase with increasing order of live offspring; this may transmit a lower risk of developing atopy to children of a higher birth order.

Our a priori hypothesis was that the birth order effect may interact with the genetic predisposition of the offspring in influencing future atopic manifestations. In order to test this hypothesis, we chose one of the most-studied "atopy genes," IL13. This choice seems appropriate because, if a gene-by-birth order interaction exists, it would be best demonstrated with a known gene. Finding such an interaction with a novel gene that does not have the kind of consistent replication so far demonstrated for IL13 may not be as scientifically robust.

This report will be the first study to investigate the interaction between these important risk factors during early childhood, as well as sustained effects in adolescence, testing both the main effects and the interaction effects of the IL13 gene with prenatal factors.

Materials and methods

Study population

The study participants were mother-infant pairs enrolled in the 1989 Isle of Wight birth cohort study. This study represents an unselected whole population birth cohort based on the Isle of Wight, U.K. The Isle of Wight is an island (13 × 23 miles) just off the South Coast of England with a resident population of 130,000. The population is stable to the extent that the majority of participants in the cohort has not moved away and was thus available for follow up. This makes the island particularly attractive for long-term prospective epidemiological studies.

Between January 1989 and February 1990, children born on the Isle of Wight, U.K. were recruited to participate in a longitudinal study (n = 1,456). The local Research Ethics Committee approved the study and informed written parental consent was obtained for all the participants at recruitment and subsequently at each follow up. This whole-population birth cohort was largely Caucasian (99%) and was living in a semi-rural, non-farming environment with no heavy industry.

The Isle of Wight birth cohort has been described in detail elsewhere [23, 24]. Briefly, pregnant women were recruited prenatally and data from birth records and extensive questionnaires were collected, including information on family history of asthma and allergy, as well as maternal smoking habits. At ages 1, 2, 4, 10 and 18 years, the original questionnaire-based information was updated, a study investigator performed physical examinations on the children, and symptoms of asthma and allergic diseases were recorded. Skin prick tests (SPT) to common food and aero-allergens were administered at the 4, 10 and 18 year follow-up visits [25], and at age 10, anti-coagulated blood samples were collected and stored frozen for subsequent DNA analysis. Also at age 10, total serum IgE and inhalant IgE screen were performed in the blood samples of the children. Further information was collected at the 10 and 18 year visits using standardized International Study of Asthma and Allergy in Childhood (ISAAC) questionnaires.

Study variables

This is a population-based association study that tests the interaction effects between IL13 SNPs and birth order on atopic outcomes. Thus, the main exposure variables involved in the interaction terms were IL13 polymorphisms and birth order. Outcome variables include skin prick tests at ages 4, 10, and 18, total IgE at birth (cord serum) and at age 10, and inhalant IgE screen at age 10. The potential confounders we adjusted for were gender, environmental tobacco smoke exposure, family SES cluster, and gestational age. Gender was classified into male and female while gestational age was assessed in weeks. Other variables are described in greater detail below.

IL13 genotyping

Genomic DNA was isolated from blood samples using QIAamp DNA Blood Kits (Qiagen, Valencia, CA) or the ABI PRISM™ 6100 Nucleic Acid PrepStation (Applied Biosystems, Foster City, CA). Polymorphisms in the IL13 gene were examined using the SNPper and Applied Biosystems databases. Genotyping was conducted by fluorogenic 5' nuclease chemistry PCR using Assays on Demand kits cycled on a 7900HT Sequence Detection System (Applied Biosystems, Foster City, CA), or biotin-streptavidin-based pyrosequencing performed on PSQ-96 instrumentation (Biotage AB, Uppsala, Sweden).

IL13 is a small gene (2.9 kb) characterized by extensive linkage disequilibrium, thus genotyping a few polymorphisms provided adequate assessment of genetic associations. Five IL13 SNPs were genotyped (Table 1): rs1800925 at the promoter region, rs2066960 at intron 1 and rs1295686 at intron 3 positions, rs20541 at the coding exon 4 position ([G] Arg to [A] Gln at amino acid position 144) and rs1295685 at the 3' UTR/exon 4 position of the gene. Each of these SNPs had a minor allele frequency ≥ 19%.

Using the Haploview program [26], three of the five single nucleotide polymorphisms (SNPs) of the interleukin 13 gene were in a linkage disequilibrium (LD) block. SNP selection was carried out using the tag-SNP approach. Since rs1295686, rs20541 and rs1295685 were in the same LD block (Figure 1 in [27]), we selected rs20541 for further analysis. In addition to being the non-synonymous SNP among these three, it is also the IL13 SNP that has been most studied [17, 21, 2831].

Table 1 Genotype proportions for IL13 single nucleotide polymorphisms (SNPs)

Birth order

Birth order was obtained from the questionnaire data collected at birth. We defined three groups for this analysis: "first", "second", and "third and higher" order birth positions. Birth order equals the number of older siblings plus one [12].

Skin Prick Testing (SPT)

Skin prick testing to 14 common food and aero-allergens was performed at 4 (n = 981), 10 (n = 1036) and 18 years (n = 845) using a standard battery of food and aeroallergens (ALK, Horsholm, Denmark), which have been previously described [32]. After 15 minutes, a mean wheal diameter ("sum of longest diameter and diameter diagonal to it" divided by 2) of at least 3 mm greater than the negative control was considered evidence of sensitization. Sensitization to at least one food or aeroallergen was recorded as a positive reaction to SPT.

Serum Total IgE Determination

Total IgE was measured in samples of cord serum (n = 1340) and serum collected at age 10 (n = 923). Total IgE in cord serum was measured using Pharmacia IgE EIA® (Pharmacia Diagnostics AB, Uppsala, Sweden) [33], which is designed to measure IgE between 0.2 to 50 kU/L on 0.1 ml of serum or plasma [34, 35]. Maternal IgE and IgE at age 10 were determined using PRIST® (Pharmacia Diagnostics AB, Uppsala, Sweden) designed to measure IgE between 2.0 to 1000 kU/L. For our analysis, maternal IgE was dichotomized into < 100 and ≥ 100 kU/L.

Maternal Atopic Status

Using information from questionnaires on maternal atopic history (Yes vs. No; n = 1213) and data from the measured maternal IgE level at birth (n = 1037), we created a composite variable with four levels: "definite" maternal atopic status (elevated IgE, positive history of atopy, n = 163); "latent" (elevated IgE, negative history, n = 201); "probable" (normal IgE, positive history, n = 121); and "none" (normal IgE, negative history, n = 552).

Serum Specific Inhalant IgE Screening

The inhalant screen was a non-quantitative test for specific allergens. The test was positive if it detected antibodies against one or more of the following allergens; house dust mite (D. pteronyssinus and D. farinae), cat dander, dog dander, horse dander, timothy grass, cladosporium, silver birch, olive, mugwort and nettle. Blood samples were allowed to stand and coagulate in Gel and Clot Activator tubes (Vacutainer Systems, Europe) for at least 10 minutes. They were then centrifuged at 3000 revolutions per minute for a further 15 minutes. Serum was then stored at minus 40°C until analysis for serum IgE. Results were recorded as either positive or negative to inhalant IgE screen. Individuals were classified as positive if they had IgE to one or more of the above tested aero-allergens.

Family Social Status Cluster

"Family social status cluster" is a composite variable that accounts for "socio-economic status" broadly defined [36]. The Isle of Wight population has been characterized as semi-rural, with most families (63%) residing in "owner-occupied" homes that have been owned by their families for decades. In order to correctly classify "social status," we chose to cluster family social status using the following three variables: a) the British socioeconomic classes (1 - 6) derived from parental occupation reported at birth; b) the number of children in the index child's bedroom (collected at age 4); and c) family income at age 10. This composite variable captures the family social class across the entire study period and has been described in more detail elsewhere [36].

Environmental Tobacco Smoke Exposure

Information on tobacco smoking by mothers (during pregnancy and later), by fathers or any other individual inside the home was recorded at recruitment and updated at each follow-up. Exposures to environmental tobacco smoke (ETS) in the household and maternal smoking during pregnancy were combined and classified into three groups. When mothers did not smoke during pregnancy and there was no exposure to household ETS in children up to the age of 10 years, children were categorized as "ETS-0". When mothers did not smoke during pregnancy but household members smoked within the home at some point up to the child's age of 10 years, the exposure status was categorized as "ETS-1". When mothers smoked during pregnancy and the children were also exposed to household ETS at some point up to the age of 10 years, the exposure was categorized as "ETS-2". None of the children had mothers who smoked during pregnancy with no exposure to household tobacco smoke after birth [28].

Data analysis and statistics

Serum total IgE at age 10 (dichotomized into ≤ 200 kU/L and >200 kU/L) and inhalant IgE at age 10 (dichotomized into positive or negative) were analyzed as binary outcomes. Because SPT was measured at ages 4, 10 and 18 years, the three measurements for each individual child were correlated over time. To account for these correlations, we estimated the effect of birth order, IL13 polymorphisms, and their interaction on SPT, using generalized estimating equations (GEE). In addition, since SPTs are prevalent outcomes, we did not estimate odds ratios but prevalence ratios (PROC GENMOD with the REPEATED statement and LOG link function in SAS). The association of birth order with IgE and the inhalant screen measured at age 10 was also determined using prevalence ratios (PRs).

In order to exclude the possible confounding effect of correlated variables, multivariable regression was used to explore the predictive effect of each variable, adjusting for all potential confounders. Assessment of interaction was carried out using backward elimination from the full model (which contained three birth order by IL13 SNP interaction terms). The significance level for the interaction effects was set at an alpha level of 0.1, for interaction on a multiplicative scale (log-linear models), and results were presented by stratification when the "birth order by gene" effect was significant. Stratification assesses interaction between IL13 and birth order on an additive scale, meaning that the combined effect is more than the sum of the single effects, whereas multiplicative interaction implies that the combined effect is more than the product of the individual effects.

Following the pattern of previous publications [27, 37, 38], we estimated the effect of IL13 polymorphisms using the dominant model, i.e. the heterozygous genotype was combined with the homozygous minor allele genotype in one category, while the homozygous common allele served as the referent group. This classification was also necessary because relatively few individuals were homozygous for the risk/minor allele, as has been reported and proposed in previous reports of the IL13 gene [38].


Participant demographic characteristics

Of the available 1536 children born during the recruiting period, 94.8% (1456/1536) consented to participate in this study. Birth order information was available for 83.2% (1212/1456) of the original birth cohort; skin prick test was performed on 67.4% (981/1456) of the original cohort at age 4, 71.2% (1036/1456) at age 10 and 58.0% (845/1456) at age 18 years. Of all the children who were skin-prick tested, the prevalence of atopy (positive reaction to one or more allergens) was 19.7%, 26.7% and 41.1% at ages 4, 10 and 18 years respectively.

Serum total IgE measurements at age 10 were available for 65.5% (953/1456) of the original cohort. Comparison of the demographic data between children with and without IgE measurements at age 10 revealed no evidence of selection bias in the major variables used for this analysis; however, children with no IgE measurements had greater environmental tobacco smoke exposure and were more likely to come from the low family social status cluster (Table 2).

Table 2 Comparison between children with and without total serum IgE measurements at age 10

Association of birth order and IL13 polymorphisms with total serum IgE and inhalant IgE screen at age 10, and with skin prick test positivity at ages 4, 10 and 18 years

Univariable analysis

Univariable analysis showed no significant association between birth order and IgE (serum total IgE and serum specific inhalant IgE screen positivity) at age 10 years (Table 3). Unadjusted bivariable analysis also showed no evidence of an association between birth order and SPT at ages 4, 10 and 18 years (Table 4). In addition, IL13 polymorphisms were not significantly correlated with atopic markers at ages 4, 10 and 18 years, except for rs2066960, which was a significant predictor of elevated serum IgE (>200 kU/L) at age 10 years (p = 0.011) (Table 3).

Other significant predictors of elevated serum IgE at age 10 years in the univariable analysis included maternal IgE (p = 0.009), maternal atopic status (p = 0.006), elevated cord serum IgE (p = 0.004) and prematurity (p = 0.008) (Table 3). Similarly, maternal history of atopy, maternal IgE, maternal atopic status and gender were significantly associated with SPT at ages 10 and 18, but not at age 4 years (Table 4). Only elevated cord serum IgE was significantly associated with skin test positivity at 4 years.

Table 3 Univariate analysis of the association of serum IgE and inhalant IgE screen positivity with IL13 and prenatal factors
Table 4 Univariate analysis of the association of positive skin prick test (SPT) with IL13 and prenatal factors

Multivariable analysis

Elevated serum IgE at age 10 years

Using multivariable log-linear regression analysis and after mutually adjusting for potential confounders, there was a significant interaction between rs1800925 and birth order on elevated serum total IgE at age 10 years (p = 0.023). To unfold this interaction, we stratified the analysis by birth order (Table 5). After stratification, the explanatory model showed the effect of rs1800925 on elevated serum total IgE at age 10 to be restricted only to first-born children (p = 0.007; adjusted Prevalence Ratio (PR) = 1.73; 95% CI = 1.16, 2.57). No other significant associations were found in this model (Table 5).

Serum specific inhalant IgE screen positivity at age 10 years

In the multivariable log-linear regression analysis and after mutually adjusting for potential confounding variables, we identified a statistically significant interaction between rs20541 and birth order on serum specific inhalant IgE positivity at age 10 (p = 0.029). To demonstrate this interaction, we stratified the analysis by birth order (Table 6). The stratified analysis showed the effect of IL13 on positive inhalant IgE screen at age 10 to be restricted only to first-born children (p = 0.034; adjusted PR = 1.48; 95% CI = 1.03, 2.13). No other significant associations were found.

Table 5 Multivariate analysis of the association of elevated total serum IgE at age 10 (>200 kU/L) with IL13 polymorphisms and prenatal factors, stratified by birth order (n = 588) (PR = Prevalence Ratio)
Table 6 Multivariate analysis of the association of inhalant serum IgE at age 10 with IL13 polymorphisms and prenatal factors, overall population, stratified by birth order (n = 588) (PR = Prevalence Ratio)
Skin prick tests at 4, 10 and 18 years: Cross-sectional analyses

We assessed the interaction of IL13 polymorphisms with SPT at the different time-points. Multivariable analysis showed significant interaction on a multiplicative scale between IL13 and birth order on SPT positivity at age 10 (rs20541, p = 0.07) and age 18 (rs20541, p = 0.030 and rs2066960, p = 0.027) but not at age 4. To further demonstrate and compare these effects on an additive scale, we stratified by birth order at each time-point. Among first-born children, minor allele carrier-ship of rs20541 showed a trend towards a higher risk of skin sensitization at age 4 (PR = 1.59; 95% CI = 1.01, 2.50), age 10 (PR = 1.48; 95% CI = 0.99, 2.20) and 18 years (PR = 1.21; 95% CI = 0.86, 1.69). On the other hand, there was no such increased risk of atopy among second-born children at age 4 (PR = 0.81; 95% CI = 0.43, 1.50), age 10 (PR = 0.70; 95% CI = 0.42, 1.15) and 18 years (PR = 0.76; 95% CI = 0.51, 1.14). Compared to the referent group, in children with a birth order of 3 or greater, there was also a higher but non-significant prevalence ratio of SPT in children with minor allele carriership of rs20541.

Skin prick tests at 4, 10 and 18 years: Repeatedmeasurement analysis

After confounder adjustment, the multivariable repeated measurement regression analysis showed a significant interaction between rs20541 and birth order on SPT at ages 4, 10 and 18 years at an alpha level of 0.1 (p = 0.076). To further describe this interaction, we stratified the statistical model by birth order (Table 7). Table 7 focuses on the effect of IL13 polymorphism on SPT within each birth order category. The stratified analysis showed that the effect of IL13 on SPT at ages 4, 10 and 18 was evident only among first-born children (p = 0.007; adjusted PR = 1.35; 95% CI = 1.09, 1.69), with borderline significance for second-born children. Children with a birth order of 3 and higher showed an increased but non-significant prevalence ratio of SPT positivity.

Comparing children from low socio-economic backgrounds with those from the high end of the social status scale, it is also evident from Tables 5 to 7 that there is a decreasing (though not statistically significant) trend of risk for atopy as birth order increases from one to three and above. First born children from low social status backgrounds had increased risk of atopy (PR = 1.51; 95% CI = 0.93, 2.45, not significant; Table 7), while third and higher order birth children were protected from the risk of atopy (PR = 0.34; 95% CI = 0.16, 0.73, p = 0.006; Table 7). In addition, gestational age was a significant predictor of skin test sensitivity for second and higher order birth children (p = 0.03 and 0.01 respectively; Table 7).

Table 7 Multivariate repeated measurement analysis of the association of positive skin prick test (SPT) at ages 4, 10 and 18 with IL13 and prenatal factors, stratified by birth order (n = 825; obs = 1,305) (PR = Prevalence Ratio)


Our analyses showed a statistically significant interaction between IL13 polymorphisms and birth order for elevated serum IgE at age 10, serum inhalant specific IgE positivity at age 10, and for SPT at ages 4, 10 and 18 years. An interaction on an additive scale was found both in the cross-sectional analysis and in the repeated measurement analysis. The predictive value of IL13 genotypes on the atopic markers was restricted only to first born children.

Our findings do not result from a selection bias for several reasons. First, the study cohort maintained high follow-up proportions throughout the entire study period: 83.7% (1,218/1456) at age 4; 94.3% (1,373/1456) at age 10; and 89.6% (1,304/1456) at age 18 years. Secondly, children who had IgE measurement results differed only with regard to environmental tobacco smoke exposure and family social status cluster from those who had no IgE measurements (Table 2). With respect to the genotypes, the presence of a selection bias could result in a violation of the Hardy-Weinberg law. The genotypes of the five IL13 SNPs were in Hardy-Weinberg equilibrium and their allele frequencies were comparable with those of other Caucasian populations [37, 39, 40] (Table 1). Hence, concerning the genetic polymorphisms, a selection bias is unlikely. In addition, due to the commencement of recruitment and assessments pre-natally, and the use of both questionnaires and physical examinations for obtaining information from the participants and their parents, information bias, if present, was minimal and non-differential. Non-differential information bias implies that any bias present would be similar in both affected and unaffected children, leading to a bias of the effect estimates towards the null.

While associations with different SNPs on the same gene are not considered an indication of disagreement [1, 41, 42], we found that two SNPs were related to atopic outcomes assessed in this study. Thus, our finding that different SNPs were associated with elevated total serum IgE (rs1800925), and with inhalant IgE and SPT (rs20541) is not surprising. Indeed, it is possible that total serum IgE is an indicator of general susceptibility (rs1800925 is at the promoter region of IL13), while the more specific reactions to inhalant allergens and skin sensitization are related to the non-synonymous (functional) SNP, rs20541 (exonic SNP).

Birth order [10, 14, 15] and IL13 polymorphisms [1621] have each been previously found to be separately associated with allergy, asthma and atopic markers in childhood. However, this is the first study that shows an interaction between birth order and IL13 polymorphisms. If confirmed by other researchers, this effect modification may in part explain the mechanism of the birth order effect: genetic polymorphisms in the IL13 gene may undergo epigenetic changes in utero due to conditions specific to a first pregnancy compared to subsequent pregnancies. It is now established that the DNA provides the blueprint for the manufacture of all the proteins necessary to create a living organism. Nonetheless, epigenetic modifications provide additional instructions on how, where, and when the genetic information will be used (gene expression). These epigenetic changes generally involve DNA modification such as methylation and acetylation, histone protein modifications, and regulation of gene expression by microRNAs [43]. Such regulatory mechanisms specify which regions of the genome are active in any given cell at any one point in time [44, 45]. Inherited changes in the "epigenome" have been postulated as a possible pathway explaining the differences in gene expression seen in individuals with identical "genomes" [46]. Such in utero epigenetic modifications may result in a higher susceptibility to increased total IgE levels and allergic sensitization in first-born children compared to children of higher-order births.

We have shown that the effect of IL13 polymorphism was restricted to firstborn children, with individuals carrying the homozygous minor and heterozygous genotypes having a higher relative risk (prevalence ratio) among firstborns. In addition, there was a lowered relative risk of the rs20451 SNP in second born, and a higher relative risk in children with a birth order of three or greater, although these were not statistically significant. The recurrence of an increased relative risk in children with 3rd or higher order births was not expected. However, since this group is comprised of children with different birth orders (three and higher), it is possible that a long inter-pregnancy interval may modify the intrauterine environment. In support of this hypothesis, Wegienka et al. [47] found that children born after inter-pregnancy intervals of less than 2 years were less likely to have a positive SPT result compared with children of mothers with no prior pregnancies. Since a successful pregnancy is a fine balance between a mother's tolerance for her fetus' foreign genetic material and the fetus' ability to survive the maternal immunologic defenses, Wegienka et al. suggested that prior pregnancies and their spacing might influence the intrauterine environment. It is possible that closely spaced pregnancies may provide a higher maternal tolerance.

Our findings suggest adverse prenatal programming in firstborn offspring. This impact may be related to immunological differences between a first-born pregnancy and later pregnancies from the same partner, as has been previously documented for the etiology of preeclampsia [13]. However, in order to distinguish the concept of hygiene hypothesis from that of prenatal programming, the greater prevalence of atopy in firstborns may warrant further investigation. To clarify this issue, one needs to investigate whether maternal IgE or cord serum IgE modifies the interaction between birth order and allergic sensitization. If these prenatal markers are shown to be effect modifiers, then the observed effect is more likely to be associated with prenatal exposures than with post-natal infections.

The advantages of a birth cohort with prolonged follow-up, such as the Isle of Wight cohort, has been discussed in the literature as a necessary means to increase our understanding of genetic diseases, especially atopic disorders [48]. Unlike other studies that assessed cord blood IgE at birth as the outcome [49] or that measured child IgE after a few months or years of follow up [50], children in the Isle of Wight birth cohort were recruited prenatally and followed to age 18 years, including repeated assessments at ages 4, 10 and 18 years. Since allergy and atopy are dynamic manifestations that can vary with age due to recoveries and relapses, the long-term predictive value of the majority of previous studies was hampered by the short-term follow up period. The use of data up to age 18 years in this paper addresses this gap.

Finally, future functional genetic research, including gene expression and methylation studies, need to determine whether birth order is related to altered expression or methylation of the IL13 gene in relation to atopic outcomes. If we know that the methylation status is affected, it may be possible to develop interventions during pregnancy that could change global or IL13-specific methylation. Thus, these findings may inform research into allergy prevention by identifying which specific conditions may modify the altered expression of the IL13 gene in firstborns when compared to second borns.


  1. 1.

    Vercelli D: Discovering susceptibility genes for asthma and allergy. Nat Rev Immunol. 2008, 8: 169-182. 10.1038/nri2257.

  2. 2.

    Allam JP, Zivanovic O, Berg C, Gembruch U, Bieber T, Novak N: In search for predictive factors for atopy in human cord blood. Allergy. 2005, 60: 743-750. 10.1111/j.1398-9995.2005.00815.x.

  3. 3.

    Braman SS: The global burden of asthma. Chest. 2006, 130: 4S-12S. 10.1378/chest.130.1_suppl.4S.

  4. 4.

    Global Initative for Asthma (GINA): The Global Burden of Asthma: December 2008. 2008, []

  5. 5.

    Escoubet-Lozach L, Glass CK, Wasserman SI: The role of transcription factors in allergic inflammation. J Allergy Clin Immunol. 2002, 110: 553-564. 10.1067/mai.2002.128076.

  6. 6.

    Gluckman PD, Hanson MA, Cooper C, Thornburg KL: Effect of in utero and early-life conditions on adult health and disease. N Engl J Med. 2008, 359: 61-73. 10.1056/NEJMra0708473.

  7. 7.

    Ege MJ, Bieli C, Frei R, van Strien RT, Riedler J, Ublagger E, Schram-Bijkerk D, Brunekreef B, van Hage M, Scheynius A, Pershagen G, Benz MR, Lauener R, von Mutius E, Braun-Fahrlander C: Prenatal farm exposure is related to the expression of receptors of the innate immunity and to atopic sensitization in school-age children. J Allergy Clin Immunol. 2006, 117: 817-823. 10.1016/j.jaci.2005.12.1307.

  8. 8.

    Martinez FD: Toward asthma prevention--does all that really matters happen before we learn to read?. N Engl J Med. 2003, 349: 1473-1475. 10.1056/NEJMe030041.

  9. 9.

    Karmaus W, Botezan C: Does a higher number of siblings protect against the development of allergy and asthma? A review. J Epidemiol Community Health. 2002, 56: 209-217. 10.1136/jech.56.3.209.

  10. 10.

    Bernsen RM, de Jongste JC, Wouden van der JC: Birth order and sibship size as independent risk factors for asthma, allergy, and eczema. Pediatr Allergy Immunol. 2003, 14: 464-469. 10.1046/j.0905-6157.2003.00108.x.

  11. 11.

    Cullinan P, MacNeill SJ, Harris JM, Moffat S, White C, Mills P, Newman Taylor AJ: Early allergen exposure, skin prick responses, and atopic wheeze at age 5 in English children: a cohort study. Thorax. 2004, 59: 855-861. 10.1136/thx.2003.019877.

  12. 12.

    Karmaus W, Arshad H, Mattes J: Does the sibling effect have its origin in utero? Investigating birth order, cord blood immunoglobulin E concentration, and allergic sensitization at age 4 years. Am J Epidemiol. 2001, 154: 909-915. 10.1093/aje/154.10.909.

  13. 13.

    Dekker GA, Robillard PY, Hulsey TC: Immune maladaptation in the etiology of preeclampsia: a review of corroborative epidemiologic studies. Obstet Gynecol Surv. 1998, 53: 377-382. 10.1097/00006254-199806000-00023.

  14. 14.

    Goldberg S, Israeli E, Schwartz S, Shochat T, Izbicki G, Toker-Maimon O, Klement E, Picard E: Asthma prevalence, family size, and birth order. Chest. 2007, 131: 1747-1752. 10.1378/chest.06-2818.

  15. 15.

    Zekveld C, Bibakis I, Bibaki-Liakou V, Pedioti A, Dimitroulis I, Harris J, Newman Taylor AJ, Cullinan P: The effects of farming and birth order on asthma and allergies. Eur Respir J. 2006, 28: 82-88. 10.1183/09031936.06.00021305.

  16. 16.

    Heinzmann A, Jerkic SP, Ganter K, Kurz T, Blattmann S, Schuchmann L, Gerhold K, Berner R, Deichmann KA: Association study of the IL13 variant Arg110Gln in atopic diseases and juvenile idiopathic arthritis. J Allergy Clin Immunol. 2003, 112: 735-739. 10.1016/S0091-6749(03)01887-6.

  17. 17.

    Hunninghake GM, Soto-Quiros ME, Avila L, Su J, Murphy A, Demeo DL, Ly NP, Liang C, Sylvia JS, Klanderman BJ, Lange C, Raby BA, Silverman EK, Celedon JC: Polymorphisms in IL13, total IgE, eosinophilia, and asthma exacerbations in childhood. J Allergy Clin Immunol. 2007, 120: 84-90. 10.1016/j.jaci.2007.04.032.

  18. 18.

    Liu X, Beaty TH, Deindl P, Huang SK, Lau S, Sommerfeld C, Fallin MD, Kao WH, Wahn U, Nickel R: Associations between total serum IgE levels and the 6 potentially functional variants within the genes IL4, IL13, and IL4RA in German children: the German Multicenter Atopy Study. J Allergy Clin Immunol. 2003, 112: 382-388. 10.1067/mai.2003.1635.

  19. 19.

    Liu X, Beaty TH, Deindl P, Huang SK, Lau S, Sommerfeld C, Fallin MD, Kao WH, Wahn U, Nickel R: Associations between specific serum IgE response and 6 variants within the genes IL4, IL13, and IL4RA in German children: the German Multicenter Atopy Study. J Allergy Clin Immunol. 2004, 113: 489-495. 10.1016/j.jaci.2003.12.037.

  20. 20.

    Liu X, Nickel R, Beyer K, Wahn U, Ehrlich E, Freidhoff LR, Bjorksten B, Beaty TH, Huang SK: An IL13 coding region variant is associated with a high total serum IgE level and atopic dermatitis in the German multicenter atopy study (MAS-90). J Allergy Clin Immunol. 2000, 106: 167-170. 10.1067/mai.2000.107935.

  21. 21.

    Sadeghnejad A, Karmaus W, Hasan Arshad S, Ewart S: IL13 gene polymorphism association with cord serum immunoglobulin E. Pediatr Allergy Immunol. 2007

  22. 22.

    Karmaus W, Arshad SH, Sadeghnejad A, Twiselton R: Does maternal immunoglobulin E decrease with increasing order of live offspring? Investigation into maternal immune tolerance. Clin Exp Allergy. 2004, 34: 853-859. 10.1111/j.1365-2222.2004.01959.x.

  23. 23.

    Kurukulaaratchy RJ, Fenn MH, Waterhouse LM, Matthews SM, Holgate ST, Arshad SH: Characterization of wheezing phenotypes in the first 10 years of life. Clin Exp Allergy. 2003, 33: 573-578. 10.1046/j.1365-2222.2003.01657.x.

  24. 24.

    Arshad SH, Hide DW: Effect of environmental factors on the development pf allergic disorders in infancy. Journal of Allergy and Clinical Immunology. 1992, 90: 235-241. 10.1016/0091-6749(92)90077-F.

  25. 25.

    Kurukulaaratchy RJ, Fenn M, Twiselton R, Matthews S, Arshad SH: The prevalence of asthma and wheezing illnesses amongst 10-year-old schoolchildren. Respir Med. 2002, 96: 163-169. 10.1053/rmed.2001.1236.

  26. 26.

    Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21: 263-265. 10.1093/bioinformatics/bth457.

  27. 27.

    Arshad SH, Karmaus W, Kurukulaaratchy R, Sadeghnejad A, Huebner M, Ewart S: Polymorphisms in the interleukin 13 and GATA binding protein 3 genes and the development of eczema during childhood. Br J Dermatol. 2008, 158: 1315-1322. 10.1111/j.1365-2133.2008.08565.x.

  28. 28.

    Sadeghnejad A, Karmaus W, Arshad SH, Kurukulaaratchy R, Huebner M, Ewart S: IL13 gene polymorphisms modify the effect of exposure to tobacco smoke on persistent wheeze and asthma in childhood, a longitudinal study. Respir Res. 2008, 9: 2-10.1186/1465-9921-9-2.

  29. 29.

    Sadeghnejad A, Meyers DA, Bottai M, Sterling DA, Bleecker ER, Ohar JA: IL13 promoter polymorphism 1112C/T modulates the adverse effect of tobacco smoking on lung function. Am J Respir Crit Care Med. 2007, 176: 748-752. 10.1164/rccm.200704-543OC.

  30. 30.

    Park HW, Lee JE, Kim SH, Kim YK, Min KU, Kim YY, Cho SH: Genetic variation of IL13 as a risk factor of reduced lung function in children and adolescents: A cross-sectional population-based study in Korea. Respir Med. 2008

  31. 31.

    Bottema RW, Reijmerink NE, Kerkhof M, Koppelman GH, Stelma FF, Gerritsen J, Thijs C, Brunekreef B, van Schayck CP, Postma DS: Interleukin 13, CD14, pet and tobacco smoke influence atopy in three Dutch cohorts: the allergenic study. Eur Respir J. 2008, 32: 593-602. 10.1183/09031936.00162407.

  32. 32.

    Sadeghnejad A, Karmaus W, Davis S, Kurukulaaratchy RJ, Matthews S, Arshad SH: Raised cord serum immunoglobulin E increases the risk of allergic sensitisation at ages 4 and 10 and asthma at age 10. Thorax. 2004, 59: 936-942. 10.1136/thx.2004.024224.

  33. 33.

    Arshad SH, Twiselton R, Smith J, Hide DW: Influence of genetic and environmental factors on the level of IgE at birth. Pediatric Allergy and Immunology. 1992, 3: 79-83. 10.1111/j.1399-3038.1992.tb00030.x.

  34. 34.

    Arshad SH, Stevens M, Hide DW: The effect of genetic and environmental factors on the prevalence of allergic disorders at the age of two years. Clinical and Experimental Allergy. 1993, 23: 504-511. 10.1111/j.1365-2222.1993.tb03238.x.

  35. 35.

    Hide DW, Arshad SH, Twiselton R, Stevens M: Cord serum IgE: an insensitive method for prediction of atopy. Clin Exp Allergy. 1991, 21: 739-743. 10.1111/j.1365-2222.1991.tb03204.x.

  36. 36.

    Ogbuanu IU, Karmaus W, Arshad SH, Kurukulaaratchy RJ, Ewart S: Effect of breastfeeding duration on lung function at age 10 years: a prospective birth cohort study. Thorax. 2009, 64: 62-66. 10.1136/thx.2008.101543.

  37. 37.

    Howard TD, Whittaker PA, Zaiman AL, Koppelman GH, Xu J, Hanley MT, Meyers DA, Postma DS, Bleecker ER: Identification and association of polymorphisms in the interleukin-13 gene with asthma and atopy in a Dutch population. Am J Respir Cell Mol Biol. 2001, 25: 377-384.

  38. 38.

    Kang MJ, Lee SY, Kim HB, Yu J, Kim BJ, Choi WA, Jang SO, Hong SJ: Association of IL-13 polymorphisms with leukotriene receptor antagonist drug responsiveness in Korean children with exercise-induced bronchoconstriction. Pharmacogenet Genomics. 2008, 18: 551-558. 10.1097/FPC.0b013e3282fe94c5.

  39. 39.

    SNPper, The Children's Hospital Informatics Program (CHIP) Bioinformatics Tools. Harvard school of medicine (Boston, MA). 2005

  40. 40.

    Graves PE, Kabesch M, Halonen M, Holberg CJ, Baldini M, Fritzsch C, Weiland SK, Erickson RP, von Mutius E, Martinez FD: A cluster of seven tightly linked polymorphisms in the IL-13 gene is associated with total serum IgE levels in three populations of white children. J Allergy Clin Immunol. 2000, 105: 506-513. 10.1067/mai.2000.104940.

  41. 41.

    Neale BM, Sham PC: The future of association studies: gene-based analysis and replication. Am J Hum Genet. 2004, 75: 353-362. 10.1086/423901.

  42. 42.

    Ober C, Hoffjan S: Asthma genetics 2006: the long and winding road to gene discovery. Genes Immun. 2006, 7: 95-100. 10.1038/sj.gene.6364284.

  43. 43.

    Chuang JC, Jones PA: Epigenetics and microRNAs. Pediatr Res. 2007, 61: 24R-29R. 10.1203/pdr.0b013e3180457684.

  44. 44.

    Wade PA, Archer TK: Epigenetics: environmental instructions for the genome. Environ Health Perspect. 2006, 114: A140-141. 10.1289/ehp.114-a140.

  45. 45.

    Bjornsson HT, Fallin MD, Feinberg AP: An integrated epigenetic and genetic approach to common human disease. Trends Genet. 2004, 20: 350-358. 10.1016/j.tig.2004.06.009.

  46. 46.

    Gosden RG, Feinberg AP: Genetics and epigenetics--nature's pen-and-pencil set. N Engl J Med. 2007, 356: 731-733. 10.1056/NEJMe068284.

  47. 47.

    Wegienka G, London SJ, Johnson CC, Ownby DR: Interpregnancy interval might affect the risk of childhood atopy. J Allergy Clin Immunol. 2004, 113: 169-171. 10.1016/j.jaci.2003.10.012.

  48. 48.

    Pembrey M: Genetic epidemiology: some special contributions of birth cohorts. Paediatr Perinat Epidemiol. 2004, 18: 3-7. 10.1111/j.1365-3016.2004.00530.x.

  49. 49.

    Lin YC, Wen HJ, Lee YL, Guo YL: Are maternal psychosocial factors associated with cord immunoglobulin E in addition to family atopic history and mother immunoglobulin E?. Clin Exp Allergy. 2004, 34: 548-554. 10.1111/j.1365-2222.2004.1928.x.

  50. 50.

    Liu CA, Wang CL, Chuang H, Ou CY, Hsu TY, Yang KD: Prenatal prediction of infant atopy by maternal but not paternal total IgE levels. J Allergy Clin Immunol. 2003, 112: 899-904. 10.1016/j.jaci.2003.08.030.

Download references


The authors gratefully acknowledge the cooperation of the children and parents who participated in this study, and appreciate the hard work of Mrs Sharon Matthews and the David Hide Asthma and Allergy Research Centre team in collecting the phenotype data. The authors thank Hans Cheng for use of Pyrosequencing equipment and Dennis Shubitowski for technical assistance.

Author information

Correspondence to Ikechukwu U Ogbuanu.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

IUO conceived of the manuscript idea, drafted the manuscript and performed the statistical analysis. WJK provided supervision and guidance in the design, analysis and writing of the manuscript. HZ provided statistical oversight while TSA assisted with the interpretation of genetic results. SE, GR and SHA participated in the design and coordination of the original cohort study, carried out the molecular genetic studies, provided oversight in interpretation of analytic results and participated in all stages of the manuscript writing process. All authors read and approved the final manuscript.

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article


  • Birth Order
  • Skin Prick Test
  • Prevalence Ratio
  • Environmental Tobacco Smoke Exposure
  • IL13 Polymorphism