The effect of genotype-by-environment interaction on longitudinal triglyceride profiles in the Bogalusa Heart Study

Triglyceride levels are a component of plasma lipids that are thought to be an important risk factor for coronary heart disease and are influenced by genetic and environmental factors, such as single nucleotide polymorphisms (SNPs), alcohol intake, and smoking. This study used longitudinal data from the Bogalusa Heart Study, a biracial community-based survey of cardiovascular disease risk factors. A sample of 1191 individuals, 4 to 38 years of age, was measured multiple times from 1973 to 2000. The study sample consisted of 730 white and 461 African American participants. Individual growth models were developed in order to assess gene-environment interactions affecting plasma triglycerides over time. After testing for inclusion of significant covariates and interactions, final models, each accounting for the effects of a different SNP, were assessed for fit and normality. After adjustment for all other covariates and interactions, LIPC -514C/T was found to interact with age3, age2, and age and a non-significant interaction of CETP -971G/A genotype with smoking status was found (p = 0.0812). Ever-smokers had higher triglyceride levels than never smokers, but persons heterozygous at this locus, about half of both races, had higher triglyceride levels after smoking cessation compared to current smokers. Since tobacco products increase free fatty acids circulating in the bloodstream, smoking cessation programs have the potential to ultimately reduce triglyceride levels for many persons. However, due to the effect of smoking cessation on the triglyceride levels of CETP -971G/A heterozygotes, the need for smoking prevention programs is also demonstrated. Both smoking cessation and prevention programs would have a great public health impact on minimizing triglyceride levels and ultimately reducing heart disease. ^

Cell survival in early embryogenesis requires proper dosage of GCN5 and p300

Histone acetylation is a central event in transcriptional activation. The importance of this modification in mammalian development is highlighted by knockout studies that revealed loss of the histone acetyltransferases GCN5, p300, or CBP results in embryonic lethality. Furthermore, early embryogenesis is sensitive to the dosage of p300 and CBP since double p300 +/−CBP+/− heterozygotes die in utero, although either single heterozygote survives. PCAF and GCN5 physically interact with p300 and CBP in vitro. To determine whether these two groups of HATs interact functionally in vivo, we created mice lacking one or more allele of p300, GCN5 or PCAF. As expected, we found that mice heterozygous for any one of these null alleles are viable. The majority of GCN5 p300 double heterozygotes also survive to adulthood with no apparent abnormalities. However, a portion of these mice die prior to birth. These embryos are developmentally stunted and exhibit increased apoptosis compared to wild type or single GCN5 or p300 heterozygous littermates at E8.5. Tissue specification is unaffected in these embryos but organ formation is compromised. In contrast, no abnormalities were observed in mice harboring mutations in both PCAF and p300 , emphasizing the specificity of HAT functions in mammalian development. ^ Since GCN5 null embryos die early in embryogenesis because of a marked increase in apoptosis, studies of its function and mechanism in late development and in tissue specific differentiation are precluded. Here, we also report the establishment of a GCN5 null embryonic stem cell line and a conditional floxGCN5 mouse line, which will serve as powerful genetic tools to examine in depth the function of GCN5 in mammalian development and in adult tissues. ^

Disease sequence from mutant rhodopsin allele to rod and cone photoreceptor degeneration in man

Mutations in the gene encoding rhodopsin, the visual pigment in rod photoreceptors, lead to retinal degeneration in species from Drosophila to man. The pathogenic sequence from rod cell-specific mutation to degeneration of rods and cones remains unclear. To understand the disease process in man, we studied heterozygotes with 18 different rhodopsin gene mutations by using noninvasive tests of rod and cone function and retinal histopathology. Two classes of disease expression were found, and there was allele-specificity. Class A mutants lead to severely abnormal rod function across the retina early in life; topography of residual cone function parallels cone cell density. Class B mutants are compatible with normal rods in adult life in some retinal regions or throughout the retina, and there is a slow stereotypical disease sequence. Disease manifests as a loss of rod photoreceptor outer segments, not singly but in microscopic patches that coalesce into larger irregular areas of degeneration. Cone outer segment function remains normal until >75% of rod outer segments are lost. The topography of cone loss coincides with that of rod loss. Most class B mutants show an inferior-nasal to superior-temporal retinal gradient of disease vulnerability associated with visual cycle abnormalities. Class A mutant alleles behave as if cytotoxic; class B mutants can be relatively innocuous and epigenetic factors may play a major role in the retinal degeneration.

Xeroderma pigmentosum and trichothiodystrophy are associated with different mutations in the XPD (ERCC2) repair/transcription gene

The xeroderma pigmentosum group D (XPD) protein has a dual function, both in nucleotide excision repair of DNA damage and in basal transcription. Mutations in the XPD gene can result in three distinct clinical phenotypes, XP, trichothiodystrophy (TTD), and XP with Cockayne syndrome. To determine if the clinical phenotypes of XP and TTD can be attributed to the sites of the mutations, we have identified the mutations in a large group of TTD and XP-D patients. Most sites of mutations differed between XP and TTD, but there are three sites at which the same mutation is found in XP and TTD patients. Since the corresponding patients were all compound heterozygotes with different mutations in the two alleles, the alleles were tested separately in a yeast complementation assay. The mutations which are found in both XP and TTD patients behaved as null alleles, suggesting that the disease phenotype was determined by the other allele. If we eliminate the null mutations, the remaining mutagenic pattern is consistent with the site of the mutation determining the phenotype.

Gilbert syndrome and glucose-6-phosphate dehydrogenase deficiency: A dose-dependent genetic interaction crucial to neonatal hyperbilirubinemia

Severe jaundice leading to kernicterus or death in the newborn is the most devastating consequence of glucose-6-phosphate dehydrogenase (EC 1.1.1.49; G-6-PD) deficiency. We asked whether the TA repeat promoter polymorphism in the gene for uridinediphosphoglucuronate glucuronosyltransferase 1 (EC 2.4.1.17; UDPGT1), associated with benign jaundice in adults (Gilbert syndrome), increases the incidence of neonatal hyperbilirubinemia in G-6-PD deficiency. DNA from term neonates was analyzed for UDPGT1 polymorphism (normal homozygotes, heterozygotes, variant homozygotes), and for G-6-PD Mediterranean deficiency. The variant UDPGT1 promoter allele frequency was similar in G-6-PD-deficient and normal neonates. Thirty (22.9%) G-6-PD deficient neonates developed serum total bilirubin ≥ 257 μmol/liter, vs. 22 (9.2%) normals (P = 0.0005). Of those with the normal homozygous UDPGT1 genotype, the incidence of hyperbilirubinemia was similar in G-6-PD-deficients and controls (9.7% and 9.9%). In contrast, in the G-6-PD-deficient neonates, those with the heterozygous or homozygous variant UDPGT1 genotype had a higher incidence of hyperbilirubinemia than corresponding controls (heterozygotes: 31.6% vs. 6.7%, P < 0.0001; variant homozygotes: 50% vs. 14.7%, P = 0.02). Among G-6-PD-deficient infants the incidence of hyperbilirubinemia was greater in those with the heterozygous (31.6%, P = 0.006) or variant homozygous (50%, P = 0.003) UDPGT1 genotype than in normal homozygotes. In contrast, among those normal for G-6-PD, the UDPGT1 polymorphism had no significant effect (heterozygotes: 6.7%; variant homozygotes: 14.7%). Thus, neither G-6-PD deficiency nor the variant UDPGT1 promoter, alone, increased the incidence of hyperbilirubinemia, but both in combination did. This gene interaction may serve as a paradigm of the interaction of benign genetic polymorphisms in the causation of disease.

Candida albicans strains heterozygous and homozygous for mutations in mitogen-activated protein kinase signaling components have defects in hyphal development

The Candida albicans genes, CST20 and HST7, were cloned by their ability to suppress the mating defects of Saccharomyces cerevisiae mutants in the ste20 and ste7 genes, which code for elements of the mating mitogen-activated protein (MAP) kinase pathway. These Candida genes are both structural and functional homologs of the cognate Saccharomyces genes. The pattern of suppression in Saccharomyces is related to their presumptive position in the MAP kinase cascade. Null alleles of these genes were constructed in Candida. The Candida homozygous null mutants are defective in hyphal formation on some media, but are still induced to form hyphae by serum, showing that serum induction of hyphae is independent of the MAP kinase cascade. The Candida heterozygotes CST20/cst20 and HST7/hst7 are also defective in hyphal formation. This lack of dominance of the wild-type allele suggests that gene dosage is important in Candida.

Disruption of the murine gene encoding phosphatidylethanolamine N-methyltransferase

All nucleated cells make phosphatidylcholine via the CDP-choline pathway. Liver has an alternative pathway in which phosphatidylcholine is made by methylation of phosphatidylethanolamine catalyzed by phosphatidylethanolamine N-methyltransferase (PEMT). We investigated the function of PEMT and its role in animal physiology by targeted disruption of its gene, Pempt2. A targeting vector that interrupts exon 2 was constructed and introduced into mice yielding three genotypes: normal (+/+), heterozygotes (+/−), and homozygotes (−/−) for the disrupted PEMT gene. Only a trace of PE methylation activity remained in Pempt2(−/−) mice. Antibody to one form of the enzyme, PEMT2, indicated complete loss of this protein from Pempt2(−/−) mice and a decrease in Pempt2(+/−) mice, compared with Pempt2(+/+) mice. The levels of hepatic phosphatidylethanolamine and phosphatidylcholine were minimally affected. The active form of CTP:phosphocholine cytidylyltransferase, the regulated enzyme in the CDP-choline pathway, was increased 60% in the PEMT-deficient mice. Injection of [l-methyl-3H]methionine demonstrated that the in vivo PEMT activity was eliminated in the Pempt2(−/−) mice and markedly decreased in the Pempt2(+/−) mice. This experiment also demonstrated that the choline moiety derived from PEMT in the liver can be distributed via the plasma throughout the mouse where it is found as phosphatidylcholine, lysophosphatidylcholine, and sphingomyelin. Mice homozygous for the disrupted Pempt2 gene displayed no abnormal phenotype, normal hepatocyte morphology, normal plasma lipid levels and no differences in bile composition. This is the first application of the “knockout mouse” technique to a gene for phospholipid biosynthesis.

Rump white inversion in the mouse disrupts dipeptidyl aminopeptidase-like protein 6 and causes dysregulation of Kit expression

The mouse rump white (Rw) mutation causes a pigmentation defect in heterozygotes and embryonic lethality in homozygotes. At embryonic day (E) 7.5, Rw/Rw embryos are retarded in growth, fail to complete neurulation and die around E 9.5. The Rw mutation is associated with a chromosomal inversion spanning 30 cM of the proximal portion of mouse chromosome 5. The Rw embryonic lethality is complemented by the W19H deletion, which spans the distal boundary of the Rw inversion, suggesting that the Rw lethality is not caused by the disruption of a gene at the distal end of the inversion. Here, we report the molecular characterization of sequences disrupted by both inversion breakpoints. These studies indicate that the distal breakpoint of the inversion is associated with ectopic Kit expression and therefore may be responsible for the dominant pigmentation defect in Rw/+ mice; whereas the recessive lethality of Rw is probably due to the disruption of the gene encoding dipeptidyl aminopeptidase-like protein 6, Dpp6 [Wada, K., Yokotani, N., Hunter, C., Doi, K., Wenthold, R. J. & Shimasaki, S. (1992) Proc. Natl. Acad. Sci. USA 89, 197–201] located at the proximal inversion breakpoint.

Strong balancing selection at HLA loci: Evidence from segregation in South Amerindian families

The genotypic proportions for major histocompatibility complex loci, HLA-A and HLA-B, of progeny in families in 23 South Amerindian tribes in which segregation for homozygotes and heterozygotes could occur are examined. Overall, there is a large deficiency of homozygotes compared with Mendelian expectations (for HLA-A, 114 observed and 155.50 expected and for HLA-B 110 observed and 144.75 expected), consistent with strong balancing selection favoring heterozygotes. There is no evidence that these deficiencies were associated with particular alleles or with the age of the individuals sampled. When these families were divided into four mating types, there was strong selection against homozygotes, averaging 0.462 for three of the mating types over the two loci. For the other mating type in which the female parent is homozygous and shares one allele with the heterozygous male parent, there was no evidence of selection against homozygotes. A theoretical model incorporating these findings surprisingly does not result in a stable polymorphism for two alleles but does result in an excess of heterozygotes and a minimum fitness at intermediate allele frequencies. However, for more than two alleles, balancing selection does occur and the model approaches the qualities of the symmetrical heterozygote advantage model as the number of alleles increases.

α+-Thalassemia protects children against disease caused by other infections as well as malaria

In the South West Pacific region, the striking geographical correlation between the frequency of α+-thalassemia and the endemicity of Plasmodium falciparum suggests that this hemoglobinopathy provides a selective advantage against malaria. In Vanuatu, paradoxically, α+-thalassemia increases the incidence of contracting mild malaria in the first 2 years of life, but severe disease was too uncommon to assess adequately. Therefore, we undertook a prospective case-control study of children with severe malaria on the north coast of Papua New Guinea, where malaria transmission is intense and α+-thalassemia affects more than 90% of the population. Compared with normal children, the risk of having severe malaria was 0.40 (95% confidence interval 0.22–0.74) in α+-thalassemia homozygotes and 0.66 (0.37–1.20) in heterozygotes. Unexpectedly, the risk of hospital admission with infections other than malaria also was reduced to a similar degree in homozygous (0.36; 95% confidence interval 0.22–0.60) and heterozygous (0.63; 0.38–1.07) children. This clinical study demonstrates that a malaria resistance gene protects against disease caused by infections other than malaria. The mechanism of the remarkable protective effect of α+-thalassemia against severe childhood disease remains unclear but must encompass the clear interaction between this hemoglobinopathy and both malarial and nonmalarial infections.

An α-Tubulin Mutant Destabilizes the Heterodimer: Phenotypic Consequences and Interactions with Tubulin-binding Proteins

Many effectors of microtubule assembly in vitro enhance the polymerization of subunits. However, several Saccharomyces cerevisiae genes that affect cellular microtubule-dependent processes appear to act at other steps in assembly and to affect polymerization only indirectly. Here we use a mutant α-tubulin to probe cellular regulation of microtubule assembly. tub1-724 mutant cells arrest at low temperature with no assembled microtubules. The results of several assays reported here demonstrate that the heterodimer formed between Tub1-724p and β-tubulin is less stable than wild-type heterodimer. The unstable heterodimer explains several conditional phenotypes conferred by the mutation. These include the lethality of tub1-724 haploid cells when the β-tubulin–binding protein Rbl2p is either overexpressed or absent. It also explains why the TUB1/tub1-724 heterozygotes are cold sensitive for growth and why overexpression of Rbl2p rescues that conditional lethality. Both haploid and heterozygous tub1-724 cells are inviable when another microtubule effector, PAC2, is overexpressed. These effects are explained by the ability of Pac2p to bind α-tubulin, a complex we demonstrate directly. The results suggest that tubulin-binding proteins can participate in equilibria between the heterodimer and its components.

Characterization of serum amyloid A protein mRNA expression and secondary amyloidosis in the domestic duck

Secondary amyloidosis is a common disease of water fowl and is characterized by the deposition of extracellular fibrils of amyloid A (AA) protein in the liver and certain other organs. Neither the normal role of serum amyloid A (SAA), a major acute phase response protein, nor the causes of secondary amyloidosis are well understood. To investigate a possible genetic contribution to disease susceptibility, we cloned and sequenced SAA cDNA derived from livers of domestic ducks. This revealed that the three C-terminal amino acids of SAA are removed during conversion to insoluble AA fibrils. Analysis of SAA cDNA sequences from several animals identified a distinct genetic dimorphism that may be relevant to susceptibility to secondary amyloid disease. The duck genome contained a single copy of the SAA gene that was expressed in liver and lung tissue of ducklings, even in the absence of induction of acute phase response. Genetic analysis of heterozygotes indicated that only one SAA allele is expressed in livers of adult birds. Immunofluorescence staining of livers from adult ducks displaying early symptoms of amyloidosis revealed what appear to be amyloid deposits within hepatocytes that are expressing unusually high amounts of SAA protein. This observation suggests that intracellular deposition of AA may represent an early event during development of secondary amyloidosis in older birds.

Mutations in the nebulin gene associated with autosomal recessive nemaline myopathy

The congenital nemaline myopathies are rare hereditary muscle disorders characterized by the presence in the muscle fibers of nemaline bodies consisting of proteins derived from the Z disc and thin filament. In a single large Australian family with an autosomal dominant form of nemaline myopathy, the disease is caused by a mutation in the α-tropomyosin gene TPM3. The typical form of nemaline myopathy is inherited as an autosomal recessive trait, the locus of which we previously assigned to chromosome 2q21.2-q22. We show here that mutations in the nebulin gene located within this region are associated with the disease. The nebulin protein is a giant protein found in the thin filaments of striated muscle. A variety of nebulin isoforms are thought to contribute to the molecular diversity of Z discs. We have studied the 3′ end of the 20.8-kb cDNA encoding the Z disc part of the 800-kDa protein and describe six disease-associated mutations in patients from five families of different ethnic origins. In two families with consanguineous parents, the patients were homozygous for point mutations. In one family with nonconsanguineous parents, the affected siblings were compound heterozygotes for two different mutations, and in two further families with one detected mutation each, haplotypes are compatible with compound heterozygosity. Immunofluorescence studies with antibodies specific to the C-terminal region of nebulin indicate that the mutations may cause protein truncation possibly associated with loss of fiber-type diversity, which may be relevant to disease pathogenesis.

Functional and physiological consequences of genetic variation at phosphoglucose isomerase: Heat shock protein expression is related to enzyme genotype in a montane beetle

Allele frequency variation at the phosphoglucose isomerase (PGI) locus in Californian populations of the beetle Chrysomela aeneicollis suggests that PGI may be undergoing natural selection. We quantified (i) apparent Michaelis-Menten constant (Km) of fructose 6-phosphate at different temperatures and (ii) thermal stability for three common PGI genotypes (1–1, 1–4, and 4–4). We also measured air temperature (Ta) and beetle body temperature (Tb) in three montane drainages in the Sierra Nevada, California. Finally, we measured 70-kDa heat shock protein (Hsp70) expression in field-collected and laboratory-acclimated beetles. We found that PGI allele 1 predominated in the northernmost drainage, Rock Creek (RC), which was also significantly cooler than the southernmost drainage, Big Pine Creek (BPC), where PGI allele 4 predominated. Allele frequencies and air temperatures were intermediate in the middle drainage, Bishop Creek (BC). Differences among genotypes in Km (1–1 > 1–4 > 4–4) and thermal stability (4–4 > 1–4 > 1–1) followed a pattern consistent with temperature adaptation. In nature, Tb was closely related to Ta. Hsp70 expression in adult beetles decreased with elevation and differed among drainages (BPC > BC > RC). After laboratory acclimation (8 days, 20°C day, 4°C night) and heat shock (4 h, 28–36°C), Hsp70 expression was greater for RC than BPC beetles. In RC, field-collected beetles homozygous for PGI 1–1 had higher Hsp70 levels than heterozygotes or a 4–4 homozygote. These results reveal functional and physiological differences among PGI genotypes, which suggest that montane populations of this beetle are locally adapted to temperature.

Nociception in cyclooxygenase isozyme-deficient mice

Prostaglandins formed by cyclooxygenase-1 (COX-1) or COX-2 produce hyperalgesia in sensory nerve endings. To assess the relative roles of the two enzymes in pain processing, we compared responses of COX-1- or COX-2-deficient homozygous and heterozygous mice with wild-type controls in the hot plate and stretching tests for analgesia. Preliminary observational studies determined that there were no differences in gross parameters of behavior between the different groups. Surprisingly, on the hot plate (55°C), the COX-1-deficient heterozygous groups showed less nociception, because mean reaction time was longer than that for controls. All other groups showed similar reaction times. In the stretching test, there was less nociception in COX-1-null and COX-1-deficient heterozygotes and also, unexpectedly, in female COX-2-deficient heterozygotes, as shown by a decreased number of writhes. Measurements of mRNA levels by reverse transcription–PCR demonstrated a compensatory increase of COX-1 mRNA in spinal cords of COX-2-null mice but no increase in COX-2 mRNA in spinal cords of COX-1-null animals. Thus, compensation for the absence of COX-1 may not involve increased expression of COX-2, whereas up-regulation of COX-1 in the spinal cord may compensate for the absence of COX-2. The longer reaction times on the hot plate of COX-1-deficient heterozygotes are difficult to explain, because nonsteroid anti-inflammatory drugs have no analgesic action in this test. Reduction in the number of writhes of the COX-1-null and COX-1-deficient heterozygotes may be due to low levels of COX-1 at the site of stimulation with acetic acid. Thus, prostaglandins made by COX-1 mainly are involved in pain transmission in the stretching test in both male and female mice, whereas those made by COX-2 also may play a role in the stretching response in female mice.