Placental telomere length decreases with gestational age and is influenced by parity: a study of third trimester live-born twins.

BACKGROUND: In contrast to the postnatal period, little is known about telomere length (TL) during prenatal life. The decrease in placental TL remains unknown, although intra uterine growth retardation and preeclampsia are associated with shorter placental TL. The aim of this study is to assess the decrease of placental TL during the third trimester of gestation and to explore the role of potential "growth influencing factors". METHODS: The study sample consisted of 329 live-born twins from the East Flanders Prospective Twin Survey. TL was determined using a multiplex quantitative PCR method. Gestational age, sex, birth order, placental characteristics, parity, maternal and paternal age, diabetes, hypertension, smoking, alcohol use, and socio economic status (SES) were considered "growth influencing factors". Bivariable multilevel regression analysis with "growth influencing factors" was performed. RESULTS: Placental TL ranged from 4.3 kbp to 84.4 kbp with a median of 10.8 kbp. Ln(TL) decreased in a linear fashion with an estimated TL decreasing from 13.98 kbp at 28 weeks to 10.56 kbp at 42 weeks. The regression coefficient of gestational age became smaller if considered together with SES (b = -0.017; p = 0.08) or diabetes (b = -0.018; p = 0.07) and bigger if considered together with parity (b = -0.022; p = 0.02), indicating that part of the association between gestational age and telomere length is explained by these three confounding factors. CONCLUSION: Placental TL decreases during the third trimester of gestation of live-born twins with approximately 25% indicating that telomere shortening may play a role in aging of the placenta.

Heritability and genetic etiology of habitual physical activity: a twin study with objective measures.

Twin studies with objective measurements suggest habitual physical activity (HPA) are modestly to highly heritable, depending on age. We aimed to confirm or refute this finding and identify relevant genetic variants using a candidate gene approach. HPA was measured for 14 days with a validated triaxial accelerometer (Tracmor) in two populations: (1) 28 monozygotic and 24 dizygotic same-sex twin pairs (aged 22 ± 5 years, BMI 21.8 ± 3.4 kg/m(2), 21 male, 31 female pairs); (2) 52 and 65 unrelated men and women (aged 21 ± 2 years, BMI 22.0 ± 2.5 kg/m(2)). Single nucleotide polymorphisms (SNPs) in PPARD, PPARGC1A, NRF1 and MTOR were considered candidates. Association analyses were performed for both groups separately followed by meta-analysis. Structural equation modeling shows significant familiality for HPA, consistent with a role for additive genetic factors (heritability 57 %, 95 % CI 32-74 %, AE model) or common environmental factors (47 %, 95 % CI 23-65 %, CE model). A moderate heritability was observed for the time spent on low- and high-intensity physical activity (P ≤ 0.05), but could not be confirmed for the time spent on moderate-intensity physical activity. For PPARD, each additional effect allele was inversely associated with HPA (P ≤ 0.01; rs2076168 allele C) or tended to be associated with HPA (P ≤ 0.05; rs2267668 allele G). Linkage disequilibrium existed between those two SNPs (alleles A/G and A/C, respectively) and meta-analysis showed that carriers of the AA GC haplotype were less physically active than carriers of the AA AA and AA AC haplotypes combined (P = 0.017). For PPARGC1A, carriers of AA in rs8192678 spent more time on high-intensity physical activity than GG carriers (P = 0.001). No associations were observed with SNPs in NRF1 and MTOR. In conclusion, HPA may be modestly heritable, which is confirmed by an association with variants in PPARD.

Is PPARα intron 7 G/C polymorphism associated with muscle strength characteristics in nonathletic young men?

Peroxisome proliferator-activated receptor alpha (PPARα), a ligand-dependent transcription factor, regulates fatty acid metabolism in heart and skeletal muscle. The intron 7 G/C polymorphism (rs4253778) has been associated with athletic performance. The rare C-allele was predominant in power athletes, whereas the G-allele was more frequent in endurance athletes. In the present study, we investigated the association between this polymorphism and strength characteristics in nonathletic, healthy young adults (n = 500; age 24.2 ± 4.4 years). Knee torque was measured during concentric knee flexion and extension movements at 60°/s, 120°/s, and 240°/s during 3, 25, and 5 repetitions, respectively. Also, resistance to muscle fatigue (i.e. work last 20% repetitions/work first 20% repetitions *100) was calculated. Differences in knee strength phenotypes between GG homozygous individuals and C-allele carriers were analyzed. The polymorphism did not influence the ability to produce isometric or dynamic knee flexor or extensor peak torque during static or dynamic conditions in this population (0.23 < P < 0.95). Similar results were found for the endurance ratio, a measure for resistance to muscle fatigue. In conclusion, the PPARα intron 7 G/C polymorphism does not seem to influence strength characteristics in a nonathletic population.

Intrauterine environment and cognitive development in young twins.

Intrauterine factors important for cognitive development, such as birth weight, chorionicity and umbilical cord characteristics were investigated. A total of 663 twin pairs completed the Wechsler Intelligence Scale for Children-Revised and scores were available for Performance, Verbal and Total Intelligence Quotient (IQ). The intrauterine factors examined were birth weight, placental weight and morphology, cord knots, cord length and cord insertion. IQ scores for the varying levels of the intrauterine markers adjusting for gender and gestational age were calculated. The heritability of IQ and the association between IQ and intrauterine environment were examined. Twins with lower birth weight and cord knots had lower IQ scores. The aetiology of IQ is largely distinct from that of birth weight and cord knots, and non-shared environment may influence the observed relationships.

Secular trends in gestational age and birthweight in twins.

BACKGROUND: In recent decades, the overall rate of preterm births has increased. The aim of the present study was to examine whether this trend is also seen for multiple gestations. More specifically, we examined if there has been a decrease in gestational age for live born monozygotic (MZ) and dizygotic (DZ) twins and if there has been a simultaneous change in birthweight. The contributions of fertility treatments and Caesarean sections were taken into consideration. All analyses were carried out in two large European twin cohorts. METHODS: Cross-sectional study of 6310 live born twin pairs, born between 1964-2007, from the Belgian East Flanders Prospective Twin Survey and 14,712 twin pairs, born between 1990-2006, from the Netherlands Twin Register. Multiple regression analyses were performed with gestational age as outcome variable, and multilevel analysis with birthweight as outcome variable. All analyses were performed with and without adjustment for zygosity, parity, maternal age, mode of conception and delivery and, for the analyses of birthweight, gestational age. RESULTS: Gestational age decreased in a linear fashion from 1964 to 2007 with a decrease of 0.25 days per year in a similar way for MZ and DZ twins. Changes in birthweight depended on gestational age: up to 32 weeks, birthweight decreased and after 32 weeks birthweight increased. The frequency of infertility treatment and Caesarean sections, primiparity and advanced maternal age increased over the years, but none of these factors influenced the secular trends in gestational age and birthweight. CONCLUSIONS: The decrease in gestational age and change in birthweight in twins are sources of concern, especially for very preterm twins, for whom birthweight decreased. For twins born after 32 weeks, an increase in birthweight was observed and this is very likely the explanation for the decrease in gestational age.

No major role for X-inactivation in variations of intelligence and behavioral problems at middle childhood.

Although members of monozygotic (MZ) twin pairs are identical in genomic sequence, epigenetic mechanisms may occasion difference in gene expression and, consequently, twin discordance in complex traits. Recent work suggests that the epigenetic process of X-inactivation in female individuals may impact on intelligence and child behavioral problems. The timing of X-inactivation has been linked to chorionic splitting in MZ twins. Dichorionic monozygotic (DC-MZ) twinning, unlike monochorionic monozygotic (MC-MZ) twinning, occurs prior to the time of X-inactivation in female organisms. Therefore, the hypothesis of a causal role of X-inactivation in intelligence and behavioral problems can be analyzed by modeling the statistical interaction between sex and chorion type for within-pair differences in these traits in MZ twins. In this study, the effect of X-inactivation on childhood behavioral problems, measured with the CBCL, was studied in a sample of 324 MZ twin pairs from the EFPTS and the effect of X-inactivation on IQ was studied in a sample of 272 twin pairs from the same twin survey. Information on chorion type, gestational age, and birth weight was additionally collated. No significant statistical interaction was found between sex and chorion type, indicating that X-inactivation is not likely involved in variations in intelligence or behavioral problems in middle childhood. Further studies are required to replicate these findings and may explore the role of X-inactivation at different ages or at the extreme scores in the spectrum of intelligence and behavioral problems or may focus on other epigenetic mechanisms.

Parent-of-origin specific linkage and association of the IGF2 gene region with birth weight and adult metabolic risk factors.

OBJECTIVE: The maternally imprinted insulin-like growth factor 2 (IGF2) gene is an important fetal growth factor and is also suggested to have postnatal metabolic effects. In this study, we examined whether common polymorphisms in IGF2 (6815_6819delAGGGC, 1156T>C and 820G>A (ApaI)) and a microsatellite marker in the close vicinity of IGF2 were linked to or associated with birth weight and adult metabolic risk factors. DESIGN AND PARTICIPANTS: Polymorphisms were genotyped in 199 monozygotic complete twin pairs, 109 dizygotic complete twin pairs, 15 single twins, 231 mothers and 228 fathers recruited from the East Flanders Prospective Twin Survey. Conventional and parent-of-origin specific linkage and association analyses were carried out with birth weight, adult body height and parameters quantifying obesity, insulin sensitivity and dyslipidaemia measured at adult age (mean age 25 years). RESULTS: In the parent-of-origin specific association analysis, in which only the paternally inherited allele was incorporated, the 1156T>C SNP (single nucleotide polymorphism) showed significant association with IGF-binding protein 1 (IGFBP1) levels (T and C (mean (95% CI)): 13.2 (12.1-14.3) and 16.2 (14.6-18.0) ng ml(-1), P=0.002). No linkage was observed in either the conventional or in the parent-of-origin specific linkage analysis. CONCLUSION: This study suggests that paternally inherited alleles of a common polymorphism in the IGF2 gene affect IGFBP1 levels.

Common SNPs in LEP and LEPR associated with birth weight and type 2 diabetes-related metabolic risk factors in twins.

OBJECTIVE: Children born small for gestational age are at increased risk of developing type 2 diabetes in adulthood. The satiety signal leptin that regulates food intake and energy expenditure might be a possible molecular link, as umbilical cord leptin levels are positively correlated with birth weight. In the present study, we examined whether common single nucleotide polymorphisms (SNPs) in the leptin (LEP; 19G>A) gene and its receptor (LEPR; Q223R and K109R) are associated with birth weight and adult metabolic risk factors for type 2 diabetes in twins. DESIGN: SNPs were genotyped in 396 monozygotic and 232 dizygotic twins (286 men and 342 women, mean age 25 years) recruited from the East Flanders Prospective Twin Survey. Data were analysed using linear mixed models. RESULTS: The LEPR K109R SNP was associated with birth weight (KK, KR and RR (95% confidence interval, CI): 2511 (2465-2557), 2575 (2516-2635) and 2726 (2606-2845) gram; P(additive)=0.001). Also the LEPR Q223R SNP showed a significant association with weight at birth (QQ, QR and RR (95% CI): 2492 (2431-2554), 2545 (2495-2595) and 2655 (2571-2740) gram; P(additive)=0.003). Furthermore, an interaction between the LEPR K109R and the Q223R SNP on birth weight was observed (P=0.014). G allele carriers of the LEP 19G>A SNP had higher high-density lipoprotein (HDL) cholesterol levels compared to 19A homozygotes (GX vs AA (95% CI): 1.62 (1.58-1.66) vs 1.49 (1.40-1.58) mmol l(-1); P(recessive)=0.013). CONCLUSIONS: This study indicates that leptin may act as a growth-promoting signal during fetal development, and suggests a possible role for the LEPR in explaining the inverse relationship between birth weight and the development of metabolic diseases in adulthood. Additionally, these results suggest that the LEP 19G>A SNP affect HDL cholesterol levels.

Genome-wide linkage scan for maximum and length-dependent knee muscle strength in young men: significant evidence for linkage at chromosome 14q24.3.

BACKGROUND: Maintenance of high muscular fitness is positively related to bone health, functionality in daily life and increasing insulin sensitivity, and negatively related to falls and fractures, morbidity and mortality. Heritability of muscle strength phenotypes ranges between 31% and 95%, but little is known about the identity of the genes underlying this complex trait. As a first attempt, this genome-wide linkage study aimed to identify chromosomal regions linked to muscle and bone cross-sectional area, isometric knee flexion and extension torque, and torque-length relationship for knee flexors and extensors. METHODS: In total, 283 informative male siblings (17-36 years old), belonging to 105 families, were used to conduct a genome-wide SNP-based multipoint linkage analysis. RESULTS: The strongest evidence for linkage was found for the torque-length relationship of the knee flexors at 14q24.3 (LOD = 4.09; p<10(-5)). Suggestive evidence for linkage was found at 14q32.2 (LOD = 3.00; P = 0.005) for muscle and bone cross-sectional area, at 2p24.2 (LOD = 2.57; p = 0.01) for isometric knee torque at 30 degrees flexion, at 1q21.3, 2p23.3 and 18q11.2 (LOD = 2.33, 2.69 and 2.21; p<10(-4) for all) for the torque-length relationship of the knee extensors and at 18p11.31 (LOD = 2.39; p = 0.0004) for muscle-mass adjusted isometric knee extension torque. CONCLUSIONS: We conclude that many small contributing genes rather than a few important genes are involved in causing variation in different underlying phenotypes of muscle strength. Furthermore, some overlap in promising genomic regions were identified among different strength phenotypes.

A novel stepwise analysis procedure of genome-wide expression profiles identifies transcript signatures of thiamine genes as classifiers of mitochondrial mutants.

To extract functional information on genes and processes from large expression datasets, analysis methods are required that can computationally deal with these amounts of data, are tunable to specific research questions, and construct classifiers that are not overspecific to the dataset at hand. To satisfy these requirements, a stepwise procedure that combines elements from principal component analysis and discriminant analysis, was developed to specifically retrieve genes involved in processes of interest and classify samples based upon those genes. In a global expression dataset of 300 gene knock-outs in Saccharomyces cerevisiae, the procedure successfully classified samples with similar 'cellular component' Gene Ontology annotations of the knock-out gene by expression signatures of limited numbers of genes. The genes discriminating 'mitochondrion' from the other subgroups were evaluated in more detail. The thiamine pathway turned out to be one of the processes involved and was successfully evaluated in a logistic model to predict whether yeast knock-outs were mitochondrial or not. Further, this pathway is biologically related to the mitochondrial system. Hence, this strongly indicates that our approach is effective and efficient in extracting meaningful information from large microarray experiments and assigning functions to yet uncharacterized genes.

Anthropometry, carbohydrate and lipid metabolism in the East Flanders Prospective Twin Survey: heritabilities.

AIMS/HYPOTHESIS: We determined the genetic contribution of 18 anthropometric and metabolic risk factors of type 2 diabetes using a young healthy twin population. METHODS: Traits were measured in 240 monozygotic (MZ) and 138 dizygotic (DZ) twin pairs aged 18 to 34 years. Twins were recruited from the Belgian population-based East Flanders Prospective Twin Survey, which is characterised by its accurate zygosity determination and extensive collection of perinatal and placental data, including information on chorionicity. Heritability was estimated using structural equation modelling implemented in the Mx software package. RESULTS: Intra-pair correlations of the anthropometric and metabolic characteristics did not differ between MZ monochorionic and MZ dichorionic pairs; consequently heritabilities were estimated using the classical twin approach. For body mass, BMI and fat mass, quantitative sex differences were observed; genetic variance explained 84, 85 and 81% of the total variation in men and 74, 75 and 70% in women, respectively. Heritability estimates of the waist-to-hip ratio, sum of four skinfold thicknesses and lean body mass were 70, 74 and 81%, respectively. The heritability estimates of fasting glucose, fasting insulin, homeostasis model assessment of insulin resistance and beta cell function, as well as insulin-like growth factor binding protein-1 levels were 67, 49, 48, 62 and 47%, in that order. Finally, for total cholesterol, LDL-cholesterol, HDL-cholesterol, total cholesterol:HDL-cholesterol ratio, triacylglycerol, NEFA and leptin levels, genetic factors explained 75, 78, 76, 79, 58, 37 and 53% of the total variation, respectively. CONCLUSIONS/INTERPRETATION: Genetic factors explain the greater part of the variation in traits related to obesity, glucose intolerance/insulin resistance and dyslipidaemia.

Heritability of somatotype components: a multivariate analysis.

OBJECTIVE: To study the genetic and environmental determination of variation in Heath-Carter somatotype (ST) components (endomorphy, mesomorphy and ectomorphy). DESIGN: Multivariate path analysis on twin data. SUBJECTS: Eight hundred and three members of 424 adult Flemish twin pairs (18-34 years of age). RESULTS: The results indicate the significance of sex differences and the significance of the covariation between the three ST components. After age-regression, variation of the population in ST components and their covariation is explained by additive genetic sources of variance (A), shared (familial) environment (C) and unique environment (E). In men, additive genetic sources of variance explain 28.0% (CI 8.7-50.8%), 86.3% (71.6-90.2%) and 66.5% (37.4-85.1%) for endomorphy, mesomorphy and ectomorphy, respectively. For women, corresponding values are 32.3% (8.9-55.6%), 82.0% (67.7-87.7%) and 70.1% (48.9-81.8%). For all components in men and women, more than 70% of the total variation was explained by sources of variance shared between the three components, emphasising the importance of analysing the ST in a multivariate way. CONCLUSIONS: The findings suggest that the high heritabilities for mesomorphy and ectomorphy reported in earlier twin studies in adolescence are maintained in adulthood. For endomorphy, which represents a relative measure of subcutaneous adipose tissue, however, the results suggest heritability may be considerably lower than most values reported in earlier studies on adolescent twins. The heritability is also lower than values reported for, for example, body mass index (BMI), which next to the weight of organs and adipose tissue also includes muscle and bone tissue. Considering the differences in heritability between musculoskeletal robustness (mesomorphy) and subcutaneous adipose tissue (endomorphy) it may be questioned whether studying the genetics of BMI will eventually lead to a better understanding of the genetics of fatness, obesity and overweight.

A battery of DNA effect biomarkers to evaluate environmental exposure of Flemish adolescents.

The present paper deals with the evaluation of a battery of genotoxicity biomarkers in healthy Flemish adolescents and their relation with common pollutants occurring in their life environment. DNA damage as reflected by the comet assay appeared to be most sensitive to ozone (partial r(2) = 0.102, p < 0.00001), and to a lesser extent to ortho-cresol (partial r(2) = 0.055; p = 0.001) and 1-hydroxy-pyrene (1-OH-pyrene, partial r(2) = 0.031; p = 0.013). 8-hydroxy-deoxyguanosine (8-OHdG) was only related to ortho-cresol (r(2) = 0.069; p < 0.007). Interestingly, the comet assay results and urinary 8-OHdG concentrations were positively correlated with a Pearson r = 0.21 (p = 0.003, N = 200). Logistic regression models revealed significant relations between chromatid breaks and 1-OH-pyrene (relative risk (RR): 1.58; p = 0.008), and t,t-muconic acid (RR: 1.71; p = 0.014). There was no correlation between micronucleus formation or occurrence of chromosomal or chromatid breaks on the one hand and comet or 8-OHdG results on the other hand. Thus, in this study the comet assay on whole blood samples and urine 8-OHdG measurements especially appeared sensitive biomarkers for assessing the genetic effects of environmental pollutants to which adolescents may be exposed.

Deconstructing the familiality of the emotive component of psychotic experiences in the general population.

OBJECTIVE: Genetic and environmental influences on variation in distress associated with subclinical psychotic experiences were examined. METHOD: A total of 289 twin pairs filled in the Community Assessment of Psychic Experiences, a self-report instrument assessing subclinical positive and negative psychotic experiences and associated distress (distresspos and distressneg). Using structural equation modelling, univariate and bivariate models were fitted. RESULTS: Univariate model fitting showed genetic and non-shared environmental influences on both distresspos and distressneg. Bivariate model fitting showed that 52% of the correlation between the two phenotypes (r=0.46) was because of shared genes and that non-shared environmental factors accounted for 48% of the correlation. CONCLUSION: Liability to psychosis not only refers to the development of psychosis per se, but also to the liability to develop dysfunctional emotional appraisals. The emotive component of psychosis liability involves genetic transmission of a general, non-symptom-specific distress factor that may be a target for molecular genetic research.

Quantitative trait loci for human muscle strength: linkage analysis of myostatin pathway genes.

This study reports the results of a multipoint linkage study that aims to unravel the genetic basis of muscle strength and muscle mass in humans. Myostatin (GDF8) is known to be a strong inhibitor of muscle growth in animals. However, studies examining human myostatin polymorphisms are rare and are limited to the GDF8 gene itself. Here, the contribution to isometric and concentric knee strength of nine key proteins involved in the myostatin pathway is studied in a nonparametric multipoint linkage analysis by means of a variance components and regression method. A sample of 367 healthy young male siblings was phenotyped on an isokinetic dynamometer and genotyped for markers of the myostatin pathway genes. Three of the loci were found significantly linked with a quantitative trait locus (QTL) for knee muscle strength. First, D13S1303 showed replication of an explorative single-point linkage study with a maximum LOD score of 2.7 (P = 0.0002). Second, maximum LOD scores of 3.4 (P = 0.00004) and 3.3 (P = 0.00005) were observed for markers D12S1042 and D12S85, respectively, at 12q12-14. Finally, marker D12S78 showed an LOD score of 2.7 at 12q22-23. We conclude that several genes involved in the myostatin pathway, but not the myostatin gene itself, are important QTLs for human muscle strength. An additional set of valuable candidate genes that were not part of the myostatin pathway was found in the chromosome 12 and 13 genomic regions.

Association of tumour necrosis factor alpha variants with the CF pulmonary phenotype.

BACKGROUND: The pulmonary phenotype in patients with cystic fibrosis (CF), even in those with the same CF transmembrane conductance regulator (CFTR) genotype, is variable and must therefore be influenced by secondary genetic factors as well as environmental factors. Possible candidate genes that modulate the CF lung phenotype may include proinflammatory cytokines. One such protein is tumour necrosis factor alpha (TNFalpha), a member of the immune system. METHODS: Three polymorphic loci in the promoter (-851c/t, -308g/a, -238g/a) and one polymorphic locus in intron 1 (+691g ins/del) of the TNFalpha gene were typed by a single nucleotide primer extension assay in CF patients and healthy controls. Spirometric data and first age of infection with Pseudomonas aeruginosa were collected retrospectively from patients' medical records. RESULTS: An association was found between the TNFalpha +691g ins/del polymorphic locus and severity of CF lung disease. Patients heterozygous for +691g ins and +691g del were more likely to have better pulmonary function (mean (SD) forced expiratory volume in 1 second (FEV1) 79.7 (12.8)% predicted) than patients homozygous for +691g ins (mean (SD) FEV1 67.5 (23.0)% predicted; p = 0.008, mean difference 12.2%, 95% CI 3.5 to 21.0). Also, patients heterozygous for +691g ins and +691g del were more likely to have an older first age of infection with P aeruginosa (mean (SD) 11.4 (6.0) years) than patients homozygous for +691g ins (mean (SD) 8.3 (4.6) years; p = 0.018, mean difference 3.1 years, 95% CI 0.5 to 5.6). An association was also found with the -851c/t polymorphic locus. In the group of patients with more severe FEV1% predicted, a higher proportion of patients were homozygous for the -851c allele than in the other group of patients (p = 0.04, likelihood ratio chi2, odds ratio = 2.4). CONCLUSION: TNFalpha polymorphisms are associated with the severity of CF lung disease in Czech and Belgian patients with CF.

Epithelioid granulomas, pattern recognition receptors, and phenotypes of Crohn's disease.

INTRODUCTION: Crohn's disease is a chronic inflammatory disorder of the gut. It is assumed that a defective interaction between the bacterial flora of the gut and the innate immune system plays a key role in the pathogenesis of the disease. This may lead to specific histological lesions. The epithelioid granuloma is particularly interesting in this regard as it is also observed in several bacterial infections of the gut. AIMS AND METHODS: We hypothesised that genetic or environmental factors with a known influence on inflammation or immunity would lead to an increased prevalence of granulomas. Therefore, surgical specimens from 161 patients were evaluated for the presence of granulomas. Patients were genotyped for the three single nucleotide polymorphisms in caspase recruitment domain 15 (CARD15)/NOD2 associated with CD and for Asp299Gly in Toll-like receptor 4 (TLR4). RESULTS: The overall prevalence of granulomas was 68.9%. We did not find a significant correlation between granulomas and TLR4 or CARD15 variants. The frequency of granulomas increased with more distal disease (63% small bowel, 72% right colon, 88% left colon, 90% rectum; p=0.01). Granulomas were more frequent in younger patients (odds ratio 0.95 (95% confidence interval 0.92-0.98) p=0.007). CONCLUSION: In this study of 161 well documented CD patients, we found no significant association between CARD15 and TLR4 variants and granulomas. This finding seems to refute our initial hypothesis. However, it may be that additional factors are needed for granuloma development. Granulomas may develop only when specific bacterial components are present. Therefore, future research on granuloma pathogenesis should be orientated towards detection and identification of bacterial components in these lesions.