Research article | Open | Open Peer Review | Published:
Association between fetal growth restriction and polymorphisms at sites -1 and +3 of pituitary growth hormone: a case-control study
BMC Pregnancy and Childbirthvolume 5, Article number: 2 (2005)
Fetal growth restriction is associated with significantly increased risks of neonatal death and morbidity and with susceptibility to hypertension, cardiovascular disease and NIDDM later in life. Human birth weight has a substantial genetic component, with at least a quarter of the variation attributable to additive genetic effects.
One hundred twenty-five subjects (83 control and 42 case) were selected using stringent inclusion/exclusion criteria. DNA sequencing was used to identify 26 single nucleotide polymorphisms in the pituitary growth hormone gene (GH1) at which all subjects were genotyped. Association with fetal growth restriction was tested by logistic regression for all sites with minor allele frequencies greater than 5%.
Logistic regression identified significant association with fetal growth restriction of C alleles at sites -1 and +3 (relative to the start of transcription) that are in complete linkage disequilibrium. These alleles are present at higher frequency (6% vs. 0.4%) in fetal growth restricted subjects and are associated with an average reduction in birth weight of 152 g in normal birth weight and 97 g in low birth weight subjects.
There is suggestive association between fetal growth restriction and the presence of C alleles at sites -1 and +3 of the pituitary growth hormone gene.
Fetal growth restriction (FGR) is a major risk factor for illness in the perinatal period and throughout life, with the smallest 7.5 percent of infants accounting for two-thirds of infant deaths . Term, low birth weight infants are at least five times more likely to die in the first year [2, 3] and are second only to premature infants in their rates of morbidity and mortality . FGR infants have an increased frequency of hypoglycemia, hypothermia, polycythemia, neurodevelopmental deficits, and cerebral palsy . Later in life, individuals born FGR are at elevated risk of hypertension, cardiovascular disease, and non-insulin dependent diabetes . For example, FGR increases the risk for adult onset of non-insulin dependent diabetes two- to three-fold [7, 8]. The ability to both diagnose and treat FGR early in gestation has enormous potential to reduce childhood and adult illness.
It is difficult to distinguish the genetic and environmental components of human birth weight variation, but recent studies support a major genetic component to birth weight variation. Clausson et al.'s  study of the offspring of female dizygotic and monozygotic twins (2,009 twin pairs) estimated heritability for human birth weight of 42%, although the confounding influence of shared environmental effects must be considered. In a study of 3,562 captive macaques that minimized environmental heterogeneity, Ha et al.  estimated a total heritability for birth weight of 51%, with an additive genetic component of 23%. These findings demonstrate that comparatively simple and readily identifiable genetic factors influence birth weight. In concert with this recent research, FGR tends both to cluster in families and to recur in successive generations [11–14].
The five genes of the human growth hormone locus reside within about 45 kb on chromosome 17 . Pituitary growth hormone (GH1) is by far the most thoroughly studied of the genes and lies at the 5' end of the cluster. The remaining four genes, placental growth hormone (GH2) and three chorionic somatomammotropins (CS1, CS2, and pseudogene CS5 or CSHP1), are expressed only from the placenta. The promoter region of GH1 is unusually polymorphic, with 16 SNPs having been identified in a span of 535 bp [16–18]. Most of these SNPs occur at the comparatively small number of sites that exhibit sequence differences among the five genes of the GH locus, and this has been interpreted as evidence of gene conversion [16, 19].
Here we use DNA sequencing to identify and to determine the frequencies of both 12 newly-identified single nucleotide polymorphisms (SNPs) in the promoter and coding region of the GH1 gene and 15 previously reported SNPs [16, 18]. Using a case-control design, we identify two SNPs in complete linkage disequilibrium near the start of transcription of the GH1 gene that may predispose to reduced birth weight.
DNA was extracted from placental tissue from 125 live births (83 normal birth weight, 42 low birth weight) at Baystate Medical Center (Springfield, MA) in a case-control study of the genetic predisposition to fetal growth restriction. All subjects were Caucasian, both Hispanic and non-Hispanic. Classification of newborns as fetal growth restricted (or IUGR, intrauterine growth retarded) followed the definition of the American College of Obstetricians and Gynecologists as those newborns below the 10th percentile of size for gestational age . Our cut-off of 2,500 g at term (>37 weeks gestation) corresponds to the lowest 7.5 percentile of all US Caucasian deliveries, or the lowest 10th percentile of female and lowest 6.5 percentile of male deliveries . Stringent inclusion/exclusion criteria (Table 1) were employed by a placental pathologist (TKB) to reduce nongenetic contributors to birth weight variation. Ethical approval to conduct this study was obtained from the Human Subjects Institutional Review Board of the University of Massachusetts.
Polymerase chain reaction (PCR) and sequencing
The region from -624 (relative to the start of transcription; GenBank accession J03071) to +1,726 (197 nucleotides after the termination codon) of GH1 was amplified in two overlapping fragments: -624 to +541 (GHN-1F 5' AGGGACCTGGGGGAGCCCCAGCAA 3', GHN-1R 5' TCACCCCTTCCTGCCACCCCTGAT 3') and +450 to +1,726 (GHN-2F 5' CCATCGTCTGCACCAGCTGGCCTT 3'; GHN-2R 5' GCCCTACAGGTTGTCTTCCCAACT 3'). Approximately 50 ng of DNA was used as a template in a polymerase chain reaction with 30 cycles of 95°C (1 minute), 62°C (1 minute), and 72°C (2 minutes 30 seconds). PCR products were purified from agarose using a QIAquick PCR Purification kit (Qiagen) and sequenced directly with BigDye v2.0 chemistry (Applied Biosystems) and run on either an ABI Prism 377 or 3100 automated DNA sequencer. PCR products were sequenced with the PCR primers and additional internal primers: 5' AAGCACAGCCAATAGATTG 3', -459 to -441; 5' GCACAAGCCCGTCAGTGGCC 3', -108 to -89; 5' GGATTTTAGGGGCGCTTACC 3', +71 to +90; 5' CATCTCCCTGCTGCTCATC 3', +931 to+949; 5' GCGCTTGGGYACTGTTCCCT 3', +1280 to +1299.
Single nucleotide polymorphism (SNP) genotyping
Sequence traces were aligned and assembled into contigs by the program Polyphred . Contigs were viewed in either the program Consed  or Sequencher (Gene Codes Corp.) and polymorphisms confirmed visually. Twenty-six polymorphic sites were identified (Table 2), including all of the sites identified by Horan et al.  in 154 British military recruits with the exception of site -339 which had a minor allele (deletion of G) frequency of 3.6% in their study.
The five genes of the human growth hormone locus exhibit high sequence similarity, and the paralogous regions corresponding to the portion of GH1 sequenced in this study (-624 to +1,726) were multiply aligned. Nucleotide positions (Table 2) were designated as invariant if all five genes had the same nucleotide or the four paralogs of GH1 were identical and matched the major allele at that site in GH1. This categorization explicitly assumes that only the minor alleles in GH1 are the product of gene conversion and that minor alleles not observed in paralogs of GH1 are the result of unique mutations. The proportion of GH1 polymorphisms at invariant versus variant sites was compared by Fisher's Exact Test to determine if there was an over-representation of polymorphic sites among variant sites.
Logistic regression, using FGR status as the outcome, was performed on gestational age and genotypes at sites with a minor allele frequency above 5% (Table 3). Sites -301 and -308 (relative to the start of transcription) are in complete linkage disequilibrium and the minor C alleles at sites -1 and +3 are in complete linkage disequilibrium. Therefore, sites -301 and +3 were excluded from logistic regression to avoid multicollinearity. Based on the results from logistic regression, separate ANOVA (Table 4) was performed on gestational age and the AA versus AC genotypes at site -1 within low birth weight and within normal birth weight subjects. Empirical p values for the F statistic for the genotypic effect, corrected for multiple comparisons, were determined by 2,000 random permutations of the genotypic data . All statistical analyses were performed using the Stata program (Stata Corp., College Station, TX).
Polymorphism in the GH1 gene
Among the 125 subjects sequenced from -624 to +1,726 of the GH1 gene, 26 polymorphic sites were identified (Table 2). These included all but one of the 15 sites characterized by Horan et al. (2003). In the region of overlap between the two studies, we failed to detect variation at site -339, where there is a minor allele deletion of a single nucleotide with frequency 3.6% in British army recruits, and identified two additional variants at sites -360 and -352 with minor allele frequencies 2.8% and 1.2%, respectively. Therefore, the discrepancies between the two studies in the identification of SNPs can most likely be ascribed to sampling error. Outside the region surveyed by Horan et al. (2003), we detected 10 additional polymorphisms. All of these had minor allele frequencies ≤ 3.2%, except an intron four polymorphism at site +1169 with a minor allele frequency of 33.1%. In general, polymorphisms in the promoter of GH1 tend to be more densely clustered and exhibit higher minor allele frequencies than in the transcribed region.
Evidence of gene conversion
An alignment of the region of GH1 sequenced in this study with the paralogous sequences of placental growth hormone and chorionic somatomammotropins indicated that there are 1,979 invariant sites and 293 variant sites (excluding 78 sites that could not be unambiguously aligned), as defined in the methods. Of the invariant sites, 5 are polymorphic in GH1, while 21 of the variant sites are polymorphic. A comparison of the proportion of polymorphic sites at invariant and variant sites by Fisher's Exact Test is highly significant (P <<0.001). This result indicates that there is a strong bias for polymorphisms in GH1 to occur at the minority of sites that exhibit sequence divergence among the paralogous genes. The high correspondence between the sequence of minor alleles in GH1 and nucleotides present in paralogs of GH1 supports previous assertions that the unusually high polymorphism of the GH1 gene is driven by gene conversion [16, 18].
A proportion of this bias may be explained by selective constraint, in that sites that are polymorphic within GH1 may be under less selective constraint and thus more free to exhibit sequence divergence among paralogs. However, to entirely explain the bias towards polymorphism at sites of divergence, one must assume that about 1,804 (91%) of the invariant sites are selectively constrained and not free to vary. Given that a substantial proportion (814 bp, ~36%) of the sequence surveyed in this study is composed of introns, this assumption seems unreasonable.
Association with fetal growth restriction
Logistic regression (Table 3) was performed on gestational age and SNPs with minor allele frequencies greater than 5%, excluding sites -301 and +3 because they are in complete linkage disequilibrium with other sites, to identify associations with fetal growth restriction. Gestational age was significant because the FGR subjects exhibit a slightly younger average estimate of gestational age (38.1 vs. 39.1 weeks, t = 4.9, p < 0.0001; gestational ages rounded to nearest week). However, even accounting for the effect of gestational age, the C allele at site -1 was significantly associated with FGR in the combined regression. However, this allele did not retain significance in a regression on only gestational age (p < 0.001) and the A/C polymorphism at site -1 (p = 0.242).
Although the A/C polymorphism at site -1 was not significant in a reduced model, we decided to investigate this site for three reasons. First, the C allele at site -1 is in complete linkage disequilibrium with C at site +3. Second, both of these sites are located at the start of transcription, making them good candidates for affecting the level of transcription of GH1. Third, the only paralog that shares C at these sites is CS-5, a pseudogene, consistent with the possibility that C at sites -1 and +3 is disadvantageous. The C allele at both sites exhibits a much higher frequency in low birth weight (6%) versus normal birth weight (0.4%) subjects. Restricting examination only within the normal or low birth weight subjects, the AC genotype is associated with an average reduction of birth weight of 152 g and 97 g in normal and low birth weight subjects, respectively. However, this difference does not achieve significant in an ANOVA on the A/C polymorphism and gestational age (Table 4).
Birth weight in humans and other primates exhibits substantial heritability [9–13]. Although a suite of environmental and genetic/karyological insults are known to cause fetal growth restriction, perhaps as many as 40% of FGR cases have no known etiology . Therefore, identification of the underlying genetic variants that predispose to FGR could have significant medical significance if it allows us to identify early in gestation those pregnancies that are at increased risk of growth retardation. A logical place to begin such a search is among those genes that are known to be major regulators of fetal growth and to exhibit significant differences in circulating protein concentrations between normal and FGR pregnancies, such as the members of the growth hormone-prolactin and insulin-IGF families of hormones, receptors, and binding proteins. Here, we use a stringently selected set of subjects to report suggestive association of SNPs in GH1 with fetal growth restriction.
Adjusted for gestational age, the C alleles at sites -1 and +3 of the GH1 gene appear to be associated with reduction in birth weight. The marginal significance of these results may be the result of several factors. First, we examined 125 total subjects, roughly one-third of them FGR, and this number may give inadequate statistical power to identify weak genetic effects. Second, the C allele is low frequency, providing a small number of heterozygotes for the allele and no homozygotes. It is worth noting that among normal birth weight subjects we observed a very similar frequency for the C allele as Horan et al. (2003) did among British army recruits (0.4% vs. 0.3%) but that among the low birth weight subjects the frequency of the C allele is substantially higher (6%). Comparing all the other allele frequencies (Table 2) between the two studies of Caucasian populations, no other allele shows such a large magnitude of difference, although the proportional difference in frequency of C at -1 and +3 may be somewhat distorted by sampling error at low allele frequencies. Third, FGR undoubtedly has a heterogeneous genetic etiology, and it may be unlikely to find any genetic variant that accounts for more than a small proportion of cases.
If the C allele at -1 and +3 exhibits true association with FGR, it unfortunately may be very difficult to determine if the effect is due to one or both sites because they are in complete linkage disequilibrium. Nevertheless, the presence of both alleles at the start of GH1 transcription provides both variants with a biologically plausible possibility to affect the level of transcription of GH1. Future work should be devoted to examining the effect of these alleles on transcription levels both alone and in combination. Site -1 is somewhat unique among single nucleotide polymorphisms in that three alternate alleles exist at that site (the major allele A and the minor alleles C and T). There are two possible explanations for this observation. First, site -1 could be a hotspot for nucleotide mutation, with transversions from A to C/T occurring often. Second, GH1 may be the recipient of gene conversion events from more than one paralog within the GH locus. Among the four other genes of the GH locus, all three alleles observed at site -1 occur (Table 2). Given the predominant effect that gene conversion [16, 18] appears to have on the patterns of nucleotide polymorphism in the GH1 gene, the latter explanation may be more plausible.
It must be pointed out that the common wisdom is that GH1 plays no role in regulating fetal growth because the GH receptor is expressed fairly late in gestation and because anencephalic infants or those born without a pituitary achieve nearly normal length . Importantly, we restricted our study to full term deliveries. It is possible that late in gestation, when fetal GH1 is expressed and GHR receptors are present in a wide variety of fetal tissues, pituitary growth hormone begins to have a growth stimulatory role sufficient to account for the 90–150 g difference in birth weight between genotypes for the A/C polymorphism at -1 and +3 that we observed.
In a stringently selected set of subjects, C alleles at sites -1 and +3 relative to the start of transcription of GH1 have a higher frequency (0.4% vs. 6%) in fetal growth restricted newborns. These two alleles are in complete linkage disequilibrium and their presence is associated with a reduction in birth weight of 152 g in term, normal birth weight subjects and 97 g in term, low birth weight (<2,500 g) subjects. In combination with environmental, behavioural and other genetic factors, these alleles may contribute to fetal growth restriction.
MacDorman MF, Atkinson JO: Infant mortality statistics from the 1997 period linked birth/infant death data set. Natl Vital Stat Rep. 1999, 47: 1-23.
McCormick MC: The contribution of low birth weight to infant mortality and childhood morbidity. N Engl J Med. 1985, 312: 82-90.
McIntire DD, Bloom SL, Casey BM, Leveno KJ: Birth weight in relation to morbidity and mortality among newborn infants. N Engl J Med. 1999, 340: 1234-1238. 10.1056/NEJM199904223401603.
Bernstein I, Gabbe SG: Intrauterine growth retardation. Obstetrics: normal and problem pregnancies. Edited by: Gabbe SG, Niebyl JR, Simpson JL and Annas GJ. 1989, New York, Churchill Livingstone, 863-886. 3rd
Hay WWJ, Catz CS, Grave GD, Yaffe SJ: Workshop summary: fetal growth: its regulation and disorders. Pediatrics. 1997, 99: 585-591. 10.1542/peds.99.4.585.
Barker DJP: Fetal and infant origins of adult disease. 1992, London, British Medical Association, 343-
Rich-Edwards JW, Colditz GA, Stampfer MJ, Willett WC, Gillman MW, Hennekens CH, Speizer FE, Manson JE: Birthweight and the risk for type 2 diabetes mellitus in adult women. Ann Intern Med. 1999, 130: 278-284.
Lithell HO, McKeigue PM, Berglund L, Mohsen R, Lithell UB, Leon DA: Relation of size at birth to non-insulin dependent diabetes and insulin concentrations in men aged 50-60 years. Bmj. 1996, 312: 406-410.
Clausson B, Lichtenstein P, Cnattingius S: Genetic influence on birthweight and gestational length determined by studies in offspring of twins. Bjog. 2000, 107: 375-381.
Ha JC, Ha RR, Almasy L, Dyke B: Genetics and caging type affect birth weight in captive pigtailed macaques (Macaca nemestrina). Am J Primatol. 2002, 56: 207-213. 10.1002/ajp.1075.
Magnus P, Bakketeig LS, Hoffman H: Birth weight of relatives by maternal tendency to repeat small-for-gestational-age (SGA) births in successive pregnancies. Acta Obstet Gynecol Scand Suppl. 1997, 165: 35-38.
Klebanoff MA, Meirik O, Berendes HW: Second-generation consequences of small-for-dates birth. Pediatrics. 1989, 84: 343-347.
Wang X, Zuckerman B, Coffman GA, Corwin MJ: Familial aggregation of low birth weight among whites and blacks in the United States. N Engl J Med. 1995, 333: 1744-1749. 10.1056/NEJM199512283332606.
Strauss RS, Dietz WH: Growth and development of term children born with low birth weight: effects of genetic and environmental factors. J Pediatr. 1998, 133: 67-72.
Chen EY, Liao YC, Smith DH, Barrera-Saldana HA, Gelinas RE, Seeburg PH: The human growth hormone locus: nucleotide sequence, biology, and evolution. Genomics. 1989, 4: 479-497. 10.1016/0888-7543(89)90271-1.
Giordano M, Marchetti C, Chiorboli E, Bona G, Momigliano Richiardi P: Evidence for gene conversion in the generation of extensive polymorphism in the promoter of the growth hormone gene. Hum Genet. 1997, 100: 249-255. 10.1007/s004390050500.
Wagner JK, Eble A, Cogan JD, Prince MA, Phillips JA, Mullis PE: Allelic variations in the human growth hormone-1 gene promoter of growth hormone-deficient patients and normal controls. Eur J Endocrinol. 1997, 137: 474-481. 10.1530/eje.0.1370474.
Horan M, Millar DS, Hedderich J, Lewis G, Newsway V, Mo N, Fryklund L, Procter AM, Krawczak M, Cooper DN: Human growth hormone 1 (GH1) gene expression: complex haplotype-dependent influence of polymorphic variation in the proximal promoter and locus control region. Hum Mutat. 2003, 21: 408-423. 10.1002/humu.10167.
Krawczak M, Chuzhanova NA, Cooper DN: Evolution of the proximal promoter region of the mammalian growth hormone gene. Gene. 1999, 237: 143-151. 10.1016/S0378-1119(99)00313-3.
ACOG: Compendium of Selected Publications. 2003, , American College of Obstetricians and Gynecologists
Oken E, Kleinman KP, Rich-Edwards J, Gillman MW: A nearly continuous measure of birth weight for gestational age using a United States national reference. BMC Pediatr. 2003, 3: 6-10.1186/1471-2431-3-6.
Nickerson DA, Tobe VO, Taylor SL: PolyPhred: automating the detection and genotyping of single nucleotide substitutions using fluorescence-based resequencing. Nucleic Acids Res. 1997, 25: 2745-2751. 10.1093/nar/25.14.2745.
Gordon D, Abajian C, Green P: Consed: a graphical tool for sequence finishing. Genome Res. 1998, 8: 195-202.
Churchill GA, Doerge RW: Empirical threshold values for quantitative trait mapping. Genetics. 1994, 138: 963-971.
Wollmann HA: Intrauterine growth restriction: definition and etiology. Horm Res. 1998, 49 Suppl 2: 1-6. 10.1159/000053079.
Gluckman PD: Growth hormone deficiency diagnosed and treated in the first 2 years of life: evidence of the role of growth hormone in human perinatal growth. Progress in growth hormone therapy. Edited by: Ranke MB and Gunnarsson R. 1994, Mannheim, J & J Verlag
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2393/5/2/prepub
This work was supported by funding from the Children's Foundation Research Center of Memphis at Le Bonheur Children's Medical Center, the Center of Genomics and Bioinformatics at the University of Tennessee Health Science Center (both to RMA), and a University of Massachusetts-Baystate Medical Center Collaborative Biomedical Research Grant (to RMA and TKB). We thank Drs. Julia Krushkal and David Hosmer for analytical advice and Jeannette Peeples and Eric Gelke for their technical work.
The author(s) declare that they have no competing interests.
RMA conceived the project and directed its design and execution. CC and RV performed the molecular genetic work and participated in preliminary analyses. TKB performed the subject selection.