|Home | About | Journals | Submit | Contact Us | Français|
To determine haplotype background of common mutations in the genes encoding surfactant proteins B and C (SFTPB and SFTPC) and to assess recombination in SFTPC.
Using comprehensive resequencing of SFTPC and SFTPB, we assessed linkage disequilibrium (LD) (D’), and computationally inferred haplotypes. We computed average recombination rates and Bayes factors (BFs) within SFTPC in a population cohort and near SFTPC (±50kb) in HapMap cohorts. We then biochemically confirmed haplotypes in families with sporadic SFTPC mutations (n=11) and in individuals with the common SFTPB mutation (121ins2, n=34).
We detected strong evidence (weak LD and BFs > 1400) for an intragenic recombination hot spot in both genes. The 121ins2 SFTPB mutation occurred predominantly (89%) on 2 common haplotypes. In contrast, no consistent haplotypes were associated with mutated SFTPC alleles. Sporadic SFTPC mutations arose on the paternal allele in 4 of 5 families; the remaining child had evidence for somatic recombination on the mutated allele.
In contrast to SFTPB, disease alleles at SFTPC do not share a common haplotype background. Most sporadic mutations in SFTPC occurred on the paternal allele, but somatic recombination may be an important mechanism of mutation in SFTPC.
Pulmonary surfactant is a phospholipid-protein complex synthesized by pulmonary alveolar type II cells that lines alveoli at the air-liquid interface, maintains alveolar patency at end expiration, and is required for successful fetal-neonatal pulmonary transition. Mutations in genes encoding two hydrophobic surfactant associated proteins, surfactant proteins-B and -C (SP-B and SP-C), result in surfactant dysfunction and lung disease in children 1. The 9.7kb SP-B gene (SFTPB) encodes a 381 amino acid proprotein that undergoes sequential proteolytic cleavages and glycosylation to yield a 79 amino acid mature peptide. The 3.5 kb SP-C gene (SFTPC) encodes an alternatively spliced 191 or 197 amino acid proprotein (proSP-C) which undergoes sequential proteolytic cleavages and palmitoylation to yield a 35 amino acid mature peptide 2. Rare, recessive loss of function mutations in SFTPB are completely penetrant and result in lethal respiratory distress syndrome in the newborn period 3,4 while dominantly expressed mutations in SFTPC are variably penetrant and may be silent or result in acute or chronic lung disease in individuals ranging from the newborn to adult 5-7.
The most common mutation in SFTPB (>60% of mutated alleles) results from a frameshift due to net insertion of 2 basepairs at codon 121 (121ins2) 3,8. The frequency of the 121ins2 mutation in the general population is 0.1%, and its frequency is enriched >100 fold in cohorts of term infants with severe neonatal respiratory distress 4,9. In a European population, the 121ins2 mutation occurred on 2 common SFTPB haplotypes, suggesting a founder effect and an older mutation despite evidence for high intragenic recombination and low linkage disequilibrium (LD) 10,11. No sporadic mutations in SFTPB in affected families studied have been identified to date.
In contrast, sporadic SFTPC mutations with variable phenotypes occur in approximately 55% of families studied to date. The most common SFTPC mutation results from a nonsynonymous T>C transition at g.1295, resulting in a threonine substitution for isoleucine at codon 73 (I73T) and accounts for approximately 25% of mutations in SFTPC 12. A previous study found no common haplotype background for I73T, suggesting that SFTPC may also have a high degree of recombination 13. To determine whether haplotype background or recombination contributes to sporadic mutation in SFTPC, we assessed linkage disequilibrium (LD) (D’) in families with sporadic SFTPC mutations and in individuals with the common SFTPB mutation (121ins2). We also computationally inferred and biochemically confirmed haplotypes and computed average recombination rates and Bayes factors (BFs) within SFTPC in a population-based cohort and near SFTPC (±50kb) in HapMap cohorts. We report that in contrast to SFTPB, we observed a male germline mutation predominance in SFTPC as well as an example of somatic recombination associated with high sporadic mutation frequency.
We bidirectionally resequenced SFTPC from lymphocyte-derived DNA obtained from two cohorts of individuals: 1) a “population cohort” of 278 infants ≥ 34 weeks gestation [117 African-American individuals (AA) and 161 European-American individuals (EA)], and 2) a “family collection” comprised of 11 children of European descent with lung disease due to mutations in SFTPC and their families (34 individuals from 11 families, 3 of whom have been included in previous publications13-15). We bidirectionally resequenced SFTPB in individuals carrying the 121ins2 mutation including newborns with SP-B deficiency (n=11 homozygous for 121ins2, and 5 compound heterozygous for 121ins2 and another loss of function mutation) and in 14 heterozygous individuals of predominantly of European descent identified through population based screening of DNA samples obtained from the newborn screening programs of Missouri, Norway and South Africa, as previously described 10,16. In all, 41 alleles from 30 individuals were analyzed. All SFTPB amplification primers are available at http://genome.wustl.edu/activity/med_seq/primers. Amplification primers for SFTPC can be located in the on-line supplementary material.
We used sequence data from the population cohort to estimate pairwise LD (as defined by D’) between all SNPs in SFTPC. Haplotypes were inferred computationally in the population cohort using PHASE and all loci with minor allele frequency (MAF) >5% (n= 22 and 30 loci in AA and EA, respectively) 17 18. Haplotypes in the family collection were inferred for all single nucleotide polymorphisms (SNPs) (including the mutations) using both PHASE and pedigree information. We also inferred haplotypes using common polymorphisms within 50 kb of SFTPC from the HapMap Yoruba (YRI) and European descent (CEU) populations (HapMap Project data release #21 http://www.hapmap.org/ , a catalog of common genetic variants that provides data about haplotypes in populations). To estimate background recombination rate and determine hot spot location, we used PHASE to compute average recombination rates and Bayes factors (BFs), where a BF>10 is evidence of a recombination hot spot 10,19.
To infer haplotypes associated with 121ins2 mutation in SFTPB, we ran PHASE with loci with MAF >5% (n=11) and the 121ins2 mutation on SFTPB sequence data from the 34 individuals with at least one 121ins2 allele, combined with sequence data from 875 European American individuals obtained through the State of Missouri Newborn Screening Program 10.
To screen for evidence of gene conversion or recombination events associated with de novo mutations, we directly ascertained the haplotypes in both the affected child and both parents by cloning long range PCR (LR-PCR) products spanning the mutated SNP from members of the family collection with sporadic SFTPC mutations. After amplification and gel purification (Montage gel extraction kit, Millipore, Bedford, MA) of a single 6.9 kb SFTPC fragment, we ligated the fragment into pCR-XL-TOPO vector (Invitrogen, Carlsbad, CA) according to the manufacturer’s instructions. Multiple clones were picked from each individual and DNA was extracted using the High Pure Plasmid Isolation Kit (Roche Diagnostics, Indianapolis, IN). We sequenced multiple clones from each individual across all heterozygous sites to directly resolve haplotypes across SFTPC in that individual. Computationally inferred 121ins2 bearing haplotypes at SFTPB were also confirmed by amplification and cloning of a single 8.5 kb SFTPB LR-PCR.
To evaluate the role of demethylation transitions as a mechanism for mutations, we compared the incidence of demethylation transitions from the ancestral allele at CG sites (...CR... where R=A/G base change or ...YG... where Y=C/T base change) in SFTPB and SFTPC as compared to SFTPD in the Seattle SNPs database (http://pga.gs.washington.edu/ ). We determined information about the ancestral allele at each SNP from dbSNP or from a BLAST search of sequence around the SNP against the chimp genome when not available in dbSNP (http://18.104.22.168/BLAST/ ).
The Washington University Human Research Protection Office approved this study. Anonymized blood samples were obtained under a protocol to evaluate the epidemiology of surfactant-associated mutations in diverse ethnic groups. We assembled, reviewed, and edited all sequence data, and we formatted all confirmed polymorphic sites for further statistical and haplotype analysis as previously described 10. We analyzed all data using Statistical Analysis System (SAS, v. 9.3.1) (SAS, Inc., Cary NC).
In a population cohort of 278 infants ≥ 34 weeks gestation [117 African-American individuals (AA) and 161 European-American individuals (EA)], we identified 90 SNPs in SFTPC, including 5 SNPs in coding regions (2 synonymous and 3 non-synonymous), 34 SNPs in intronic regions, and 51 SNPs in the promoter region. We detected no known deleterious variants in this cohort. In a family cohort of 11 children of European descent with lung disease due to mutations in SFTPC and their families (34 individuals from 11 families), we found 8 potentially deleterious variants in children with interstitial lung disease (Table I): these variants were non-synonymous and were predicted to be damaging by 2 homology-based software tools (SIFT or Polyphen; http://blocks.fhcrc.org/sift/SIFT.html; http://genetics.bwh.harvard.edu/pph/), or involved splice sites (donor/acceptor) that could alter exon splicing. Our failure to find these variants in 556 alleles in the population cohort and their locations in evolutionarily conserved regions suggest they may functionally account for the interstitial lung disease in these children.
Despite its small genomic size, we found weak LD within and near SFTPC in both African-American and European-American populations, which suggests that SFTPC includes a recombination hot spot. Computationally inferred haplotype diversity using the common SNPs in the population cohort revealed 104 and 136 haplotypes in the AA and EA populations, respectively, among 117 and 195 individuals from each population, which suggests significant amounts of intragenic recombination. PHASE inferred a very high background recombination rate in the resequencing population, so there was little evidence of a hotspot considering this 6 kb region (low BF of 5.3 and 0.24 respectively). However, when we expanded the genomic context by using HapMap data to examine recombination across 100 kb surrounding SFTPC we obtained BF of 1446 and 1416 in the CEU and YRI populations, respectively (Figure I). Thus, a strong recombination hotspot spans most, if not all, of the SFTPC transcript. We found no common haplotype associated with either the mutated or nonmutated alleles in either the children with known SFTPC mutations or their family members (Table II), suggesting absence of a founder effect.
Among children with interstitial lung disease due to mutations in SFTPC, we found 5 of 11 to have alleles not inherited from either parent that presumably represent sporadic mutations. Biochemical haplotype determination revealed that the sporadic mutation occurred on the paternal allele in 4 of 5 families. Because the transmitted paternal haplotype could be completely phased, we conclude that none of these mutations appears to have occurred coincident with a meiotic recombination event in the paternal germ line. Neither haplotype found in the fifth child occurred in either parent (Table III), and both haplotypes were confirmed in multiple independent clones (mutated haplotype in two clones and nonmutated in three). However, the observed haplotypes are consistent with a somatic recombination event between the inherited maternal and paternal chromosomes. There were no clinical characteristics that distinguished this child’s presentation from those of the other children (Table I).
We identified five SFTPB haplotypes associated with the 121ins2 mutation, two of which represented 90.2% of all 121ins2 haplotypes in this cohort (37/41). A third haplotype accounted for 4.9% of the haplotypes (2/41) and was seen in 2 unrelated families. The other two haplotypes were each seen once in the same individual. We confirmed all 5 haplotypes biochemically. Using 6 SNPs from a study in a Western European population by Tredano et al. (g.-18, 1013, 1580, 4546/4550, 4559/4563, 4564/4568), we found the same two common 121ins2 bearing haplotypes that they reported plus two unique haplotypes in this study 11 (Table IV).
Neither SFTPB nor SFTPC showed a significant excess of demethylation transitions from the ancestral allele, when compared to SFTPD.
SFTPB and SFTPC have important biologic and genomic similarities. Both genes are small (<10 kb) and encode hydrophobic proteins that are critical for function of the pulmonary surfactant. Comprehensive resequencing of SFTPB and SFTPC revealed excess low frequency variation and evidence of intragenic recombination hot spots in both genes 10,20. These similarities raise the possibility that intragenic recombination hotspots are mechanistically related to increased rare variant frequency.
However, mutation inheritance and expression are different in SFTPB and SFTPC. Mutations in SFTPB are inherited in an autosomal recessive manner and are completely penetrant. Heterozygous carriers are asymptomatic during infancy and adulthood 21. Because SFTPB disease alleles can persist silently in the population, affected individuals are more likely to inherit two disease alleles (one from each parent) than to inherit one disease allele and spontaneously mutate the other allele. Consistent with this finding is the limited number of haplotypes upon which the most common SFTPB mutation, 121ins2, occurs, an observation that suggests a founder effect for this mutation, with all observed 121ins2 alleles sharing a single common ancestor in recent history.
In contrast, SFTPC mutations are dominant with variable penetrance and approximately 50% of affected patients have sporadic mutations. Because a single mutated SFTPC allele may result in a lung disease phenotype, we observed a number of sporadic mutations in the clinically affected population but no mutations in the population cohort. In contrast to the founder effect observed in SFTPB, the most common SFTPC mutation (I73T) likely represents a recurrent mutation, because it is observed on a variety of haplotype backgrounds 13.
We sought to determine whether the mechanism of sporadic mutations in SFTPC is related to methylation or intragenic recombination. We found no evidence that either SFTPC or SFTPB is associated with an elevated rate of demethylation associated transitions at CpG dinucleotides in either SFTPB or SFTPC, a common mechanism of mutations across the genome 22. In examining families with sporadic mutations in SFTPC, we began with the hypothesis that the mutations arose as a consequence of meiotic recombination or gene conversion events in the parental germline. In meiosis, cells must orient chromosomes accurately on the meiotic spindle via physical connection by deliberately causing DNA damage and relying on a subsequent repair process 23. It is through this DNA repair process, during the resolution of Halladay junctions, that the potential for mutation events arises. Because we could resolve the parental alleles in 4 of the 5 children with a sporadic mutation, we conclude that meiotic recombination events within SFTPC were not associated with the mutation. Thus, neither CpG demethylation nor meiotic recombination appears to account for the high de novo mutation rate in SFTPC.
The fact that the SFTPC sporadic mutations were observed on the paternal allele in four of the five families is consistent with male mutation bias, first proposed by Haldane in 1935. A male bias in the origin of new mutations should occur if the number of cell divisions between fertilized egg and gamete is higher in males than in females, as observed in humans 24. This bias has been previously demonstrated in Multiple Endocrine Neoplasia 2B patients in whom Carlson et al. demonstrated that 25 of 25 sporadic cases occurred due to mutations on the paternally derived chromosome 25. Although we also observe a trend toward male bias, our sample of sporadic mutations is too small to reject the null hypothesis of equal mutation rates in both sexes.
Our data for the four paternally derived sporadic alleles are not consistent with a recombination or large gene conversion event in the area of the mutation, because paternal haplotypes of markers flanking the mutation bearing haplotype are intact. Thus, paternal meiotic recombination/gene conversion is not obviously driving the recurrent mutations. Gene conversion frequently occurs on a scale of 500bp-2kb 26, so the possibility remains that these mutations arose during undetected gene conversion events, but the trend toward paternal bias suggests that the mutations arose during mitotic divisions.
The del91-93 mutation in the fifth child appears to be a somatic mutation because the two haplotypes recovered from the child, although complementary, are not consistent with either parent. Rather, the haplotypes carried by the child appear to be the products of recombination between a maternal and a paternal haplotype, which could only occur after fertilization. Somatic (or mitotic) recombination also involves repair of double strand breaks and is a frequent occurrence in the mammalian genome as homologous recombination occurs 10-6 to 10-5 times per cell cycle between repeated DNA sequences. If recombination occurs in roughly 1-5kb segments, approximately 10 spontaneous homologous recombination events occur per cell in each cell cycle 27. It is not clear whether the meiotic recombination hotspot in SFTPC is also a somatic recombination hotspot, but if it is, then somatic mutations in SFTPC may contribute to pediatric interstitial lung disease more than previously recognized and may explain some of the variability in disease expression.
Somatic recombination as a mechanism for mutation in non-cancerous diseases has been reported in neurofibromatosis 1 and 2, McCune-Albright syndrome, Paroxysmal nocturnal hemoglobinuria, and males with incontinentia pigmenti. Disorders with proven somatic and germ-line mosaicism include osteogenesis imperfecta II, Duchenne muscular dystrophy, Hunter syndrome, retinoblastoma, neurofibromatosis 2, congenital contractural arachnodactyly, hemophilia B, and hemophilia A 28. In 15 patients with idiopathic atrial fibrillation, Gollob et al. found four mutations in the connexin 40 gene, three of which were present in cardiac tissue only, consistent with somatic mutation. The fourth mutation was found in both heart tissue and lymphocyte DNA, consistent with a germ-line mutation 29.
We believe that the somatic mutation in our patient occurred at a very early stage of development because the majority of both lung and blood cells must have arisen from the mutated cell. Because the child came to clinical attention with interstitial lung disease, both endodermal-derived lung tissue and mesodermal-derived peripheral blood lymphocytes must carry the mutation. The evidence of somatic mutation in this family raises the possibility that mutations at a later stage of development than the one in this family may be detectable to a different extent in lung tissue and peripheral blood. For example, Beck et al. identified a patient with spontaneous early-onset familial Alzheimer’s Disease due to a somatic mutation in the presenilin-1 gene that was detectable in 8% of peripheral lymphocytes and 14% of the cerebral cortex, but enough to cause phenotypic disease 30. As peripheral blood lymphocytes are the most common DNA source used for resequencing, individuals without detectable mutations by resequencing of lymphocyte-derived DNA may carry somatic mutations that could be identified through the resequencing of tissue-specific DNA. Furthermore, this evidence suggests that mutations previously thought to be sporadic may actually be inherited from a parent who carried the mutated allele in the germline but not in lymphocytes. Thus, in the absence of an identifiable mutation in lymphocyte-derived DNA, sequencing SFPTC directly from lung tissue of patients with both sporadic and familial ILD may provide specific insight into the etiology of disease. Studies are underway to determine the concordance of SFTPC genotype between lung tissue and lymphocyte-derived DNA.
Further studies may help to elucidate whether mitosis associated mutation could contribute to other diseases with high sporadic mutation rates, including Tuberous Sclerosis Complex (autosomal dominant with mutation rate approximately 2/3) 31, Multiple Endocrine Neoplasia Type 2B (autosomal dominant, 50% sporadic mutations) 32, Cornelia de Lange syndrome (autosomal dominant, 99% sporadic mutations) 33, and epidermolytic hyperkeratosis (autosomal dominant, 50% sporadic) 34. Better understanding of the mechanisms of mutation may aid in diagnosis of patients with previously negative screening results and will allow for better prospective guidance and reproductive counseling for affected families.
The authors would like to thank Yumi Kasai, PhD of the Genome Sequencing Center at Washington University for overseeing the resequencing of SFTPC in the population cohort, Kimberly Fitzroy for sequencing the 121ins2 carriers, and Peter van Asperen, MD of the University of Sydney, Australia, James Chmiel, MD of Case Western Reserve University, Cleveland, OH, and Jan Kouwenberg, MD, of HAGA Juliana Children’s Hospital, The Hague, Netherlands, for providing DNA samples and clinical data.
This work was supported by grants from the National Heart, Lung, and Blood Institute (HL 82747 and HL 65174 to F.S.C., HL 65385 to A.H.), the Saigh Foundation (FSC and AH), NIH Training Grant T32-HD043010 (ADM), Discovery Labs Pulmonary Research Grant (ADM), NIH Training Grant T32-HD041925 (JW), and Washington University Genome Sequencing Center Pilot Scale Sequencing Program (AH).
This work was presented, in part, at the Pediatric Academic Societies’ Annual Meeting, Toronto, Canada, May 2007 (E-PAS2007:616292.10).
Conflict of interest statement: Dr. McBee’s support from Discovery Labs, Inc. did not introduce any bias into this work. None of the authors has a conflict of interest.