Cystic fibrosis-related diabetes onset can be predicted using biomarkers measured at birth

Purpose

Cystic fibrosis (CF), caused by pathogenic variants in the CF transmembrane conductance regulator (CFTR), affects multiple organs including the exocrine pancreas, which is a causal contributor to cystic fibrosis–related diabetes (CFRD). Untreated CFRD causes increased CF-related mortality whereas early detection can improve outcomes.

Methods

Using genetic and easily accessible clinical measures available at birth, we constructed a CFRD prediction model using the Canadian CF Gene Modifier Study (CGS; n = 1,958) and validated it in the French CF Gene Modifier Study (FGMS; n = 1,003). We investigated genetic variants shown to associate with CF disease severity across multiple organs in genome-wide association studies.

Results

The strongest predictors included sex, CFTR severity score, and several genetic variants including one annotated to PRSS1, which encodes cationic trypsinogen. The final model defined in the CGS shows excellent agreement when validated on the FGMS, and the risk classifier shows slightly better performance at predicting CFRD risk later in life in both studies.

Conclusion

We demonstrated clinical utility by comparing CFRD prevalence rates between the top 10% of individuals with the highest risk and the bottom 10% with the lowest risk. A web-based application was developed to provide practitioners with patient-specific CFRD risk to guide CFRD monitoring and treatment.

Positive epistasis between disease-causing missense mutations and silent polymorphism with effect on mRNA translation velocity

Epistasis refers to the dependence of a mutation on other mutation(s) and the genetic context in general. In the context of human disorders, epistasis complicates the spectrum of disease symptoms and has been proposed as a major contributor to variations in disease outcome. The nonadditive relationship between mutations and the lack of complete understanding of the underlying physiological effects limit our ability to predict phenotypic outcome. Here, we report positive epistasis between intragenic mutations in the cystic fibrosis transmembrane conductance regulator (CFTR)-the gene responsible for cystic fibrosis (CF) pathology. We identified a synonymous single-nucleotide polymorphism (sSNP) that is invariant for the CFTR amino acid sequence but inverts translation speed at the affected codon. This sSNP in cis exhibits positive epistatic effects on some CF disease-causing missense mutations. Individually, both mutations alter CFTR structure and function, yet when combined, they lead to enhanced protein expression and activity. The most robust effect was observed when the sSNP was present in combination with missense mutations that, along with the primary amino acid change, also alter the speed of translation at the affected codon. Functional studies revealed that synergistic alteration in ribosomal velocity is the underlying mechanism; alteration of translation speed likely increases the time window for establishing crucial domain-domain interactions that are otherwise perturbed by each individual mutation.

LocusFocus: Web-based colocalization for the annotation and functional follow-up of GWAS

Genome-wide association studies (GWAS) have primarily identified trait-associated loci in the non-coding genome. Colocalization analyses of SNP associations from GWAS with expression quantitative trait loci (eQTL) evidence enable the generation of hypotheses about responsible mechanism, genes and tissues of origin to guide functional characterization. Here, we present a web-based colocalization browsing and testing tool named LocusFocus (https://locusfocus.research.sickkids.ca). LocusFocus formally tests colocalization using our established Simple Sum method to identify the most relevant genes and tissues for a particular GWAS locus in the presence of high linkage disequilibrium and/or allelic heterogeneity. We demonstrate the utility of LocusFocus, following up on a genome-wide significant locus from a GWAS of meconium ileus (an intestinal obstruction in cystic fibrosis). Using LocusFocus for colocalization analysis with eQTL data suggests variation in ATP12A gene expression in the pancreas rather than intestine is responsible for the GWAS locus. LocusFocus has no operating system dependencies and may be installed in a local web server. LocusFocus is available under the MIT license, with full documentation and source code accessible on GitHub at https://github.com/naim-panjwani/LocusFocus.

Genetic Modifiers of Cystic Fibrosis-Related Diabetes Have Extensive Overlap With Type 2 Diabetes and Related Traits

CONTEXT

Individuals with cystic fibrosis (CF) develop a distinct form of diabetes characterized by β-cell dysfunction and islet amyloid accumulation similar to type 2 diabetes (T2D), but generally have normal insulin sensitivity. CF-related diabetes (CFRD) risk is determined by both CFTR, the gene responsible for CF, and other genetic variants.

OBJECTIVE

To identify genetic modifiers of CFRD and determine the genetic overlap with other types of diabetes.

DESIGN AND PATIENTS

A genome-wide association study was conducted for CFRD onset on 5740 individuals with CF. Weighted polygenic risk scores (PRSs) for type 1 diabetes (T1D), T2D, and diabetes endophenotypes were tested for association with CFRD.

RESULTS

Genome-wide significance was obtained for variants at a novel locus (PTMA) and 2 known CFRD genetic modifiers (TCF7L2 and SLC26A9). PTMA and SLC26A9 variants were CF-specific; TCF7L2 variants also associated with T2D. CFRD was strongly associated with PRSs for T2D, insulin secretion, postchallenge glucose concentration, and fasting plasma glucose, and less strongly with T1D PRSs. CFRD was inconsistently associated with PRSs for insulin sensitivity and was not associated with a PRS for islet autoimmunity. A CFRD PRS comprising variants selected from these PRSs (with a false discovery rate < 0.1) and the genome-wide significant variants was associated with CFRD in a replication population.

CONCLUSIONS

CFRD and T2D have more etiologic and mechanistic overlap than previously known, aligning along pathways involving β-cell function rather than insulin sensitivity. Two CFRD risk loci are unrelated to T2D and may affect multiple aspects of CF. An 18-variant PRS stratifies risk of CFRD in an independent population.

Correlating Cystic Fibrosis Transmembrane Conductance Regulator Function with Clinical Features to Inform Precision Treatment of Cystic Fibrosis

Rationale: The advent of precision treatment for cystic fibrosis using small-molecule therapeutics has created a need to estimate potential clinical improvements attributable to increases in cystic fibrosis transmembrane conductance regulator (CFTR) function. Objectives: To derive CFTR function of a variety of CFTR genotypes and correlate with key clinical features (sweat chloride concentration, pancreatic exocrine status, and lung function) to develop benchmarks for assessing response to CFTR modulators. Methods: CFTR function assigned to 226 unique CFTR genotypes was correlated with the clinical data of 54,671 individuals enrolled in the Clinical and Functional Translation of CFTR (CFTR2) project. Cross-sectional FEV₁% predicted measurements were plotted by age at which measurement was obtained. Shifts in sweat chloride concentration and lung function reported in CFTR modulator trials were compared with function-phenotype correlations to assess potential efficacy of therapies. Measurements and Main Results: CFTR genotype function exhibited a logarithmic relationship with each clinical feature. Modest increases in CFTR function related to differing genotypes were associated with clinically relevant improvements in cross-sectional FEV₁% predicted over a range of ages (6-82 yr). Therapeutic responses to modulators corresponded closely to predictions from the CFTR2-derived relationship between CFTR genotype function and phenotype. Conclusions: Increasing CFTR function in individuals with severe disease will have a proportionally greater effect on outcomes than similar increases in CFTR function in individuals with mild disease and should reverse a substantial fraction of the disease process. This study provides reference standards for clinical outcomes that may be achieved by increasing CFTR function.

Mammographic density defined by higher than conventional brightness thresholds better predicts breast cancer risk

Background: Mammographic density defined by the conventional pixel brightness threshold, and adjusted for age and body mass index (BMI), is a well-established risk factor for breast cancer. We asked if higher thresholds better separate women with and without breast cancer.

Methods: We studied Australian women, 354 with breast cancer over-sampled for early-onset and family history, and 944 unaffected controls frequency-matched for age at mammogram. We measured mammographic dense area and percent density using the CUMULUS software at the conventional threshold, which we call Cumulus, and at two increasingly higher thresholds, which we call Altocumulus and Cirrocumulus, respectively. All measures were Box–Cox transformed and adjusted for age and BMI. We estimated the odds per adjusted standard deviation (OPERA) using logistic regression and the area under the receiver operating characteristic curve (AUC).

Results:Altocumulus and Cirrocumulus were correlated with Cumulus (r ∼ 0.8 and 0.6, respectively). For dense area, the OPERA was 1.62, 1.74 and 1.73 for Cumulus, Altocumulus and Cirrocumulus, respectively (all P < 0.001). After adjusting for Altocumulus and Cirrocumulus, Cumulus was not significant (P > 0.6). The OPERAs for percent density were less but gave similar findings. The mean of the standardized adjusted Altocumulus and Cirrocumulus dense area measures was the best predictor; OPERA = 1.87 [95% confidence interval (CI): 1.64–2.14] and AUC = 0.68 (0.65–0.71).

Conclusions: The areas of higher mammographically dense regions are associated with almost 30% stronger breast cancer risk gradient, explain the risk association of the conventional measure and might be more aetiologically important. This has substantial implications for clinical translation and molecular, genetic and epidemiological research.

Sources of Variation in Sweat Chloride Measurements in Cystic Fibrosis

RATIONALE

Expanding the use of cystic fibrosis transmembrane conductance regulator (CFTR) potentiators and correctors for the treatment of cystic fibrosis (CF) requires precise and accurate biomarkers. Sweat chloride concentration provides an in vivo assessment of CFTR function, but it is unknown the degree to which CFTR mutations account for sweat chloride variation.

OBJECTIVES

To estimate potential sources of variation for sweat chloride measurements, including demographic factors, testing variability, recording biases, and CFTR genotype itself.

METHODS

A total of 2,639 sweat chloride measurements were obtained in 1,761 twins/siblings from the CF Twin-Sibling Study, French CF Modifier Gene Study, and Canadian Consortium for Genetic Studies. Variance component estimation was performed by nested mixed modeling.

MEASUREMENTS AND MAIN RESULTS

Across the tested CF population as a whole, CFTR gene mutations were found to be the primary determinant of sweat chloride variability (56.1% of variation) with contributions from variation over time (e.g., factors related to testing on different days; 13.8%), environmental factors (e.g., climate, family diet; 13.5%), other residual factors (e.g., test variability; 9.9%), and unique individual factors (e.g., modifier genes, unique exposures; 6.8%) (likelihood ratio test, P < 0.001). Twin analysis suggested that modifier genes did not play a significant role because the heritability estimate was negligible (H2 = 0; 95% confidence interval, 0.0-0.35). For an individual with CF, variation in sweat chloride was primarily caused by variation over time (58.1%) with the remainder attributable to residual/random factors (41.9%).

CONCLUSIONS

Variation in the CFTR gene is the predominant cause of sweat chloride variation; most of the non-CFTR variation is caused by testing variability and unique environmental factors. If test precision and accuracy can be improved, sweat chloride measurement could be a valuable biomarker for assessing response to therapies directed at mutant CFTR.

Cystic fibrosis gene modifier SLC26A9 modulates airway response to CFTR-directed therapeutics

Cystic fibrosis is realizing the promise of personalized medicine. Recent advances in drug development that target the causal CFTR directly result in lung function improvement, but variability in response is demanding better prediction of outcomes to improve management decisions. The genetic modifier SLC26A9 contributes to disease severity in the CF pancreas and intestine at birth and here we assess its relationship with disease severity and therapeutic response in the airways. SLC26A9 association with lung disease was assessed in individuals from the Canadian and French CF Gene Modifier consortia with CFTR-gating mutations and in those homozygous for the common Phe508del mutation. Variability in response to a CFTR-directed therapy attributed to SLC26A9 genotype was assessed in Canadian patients with gating mutations. A primary airway model system determined if SLC26A9 shows modification of Phe508del CFTR function upon treatment with a CFTR corrector. In those with gating mutations that retain cell surface-localized CFTR we show that SLC26A9 modifies lung function while this is not the case in individuals homozygous for Phe508del where cell surface expression is lacking. Treatment response to ivacaftor, which aims to improve CFTR-channel opening probability in patients with gating mutations, shows substantial variability in response, 28% of which can be explained by rs7512462 in SLC26A9 (P = 0.0006). When homozygous Phe508del primary bronchial cells are treated to restore surface CFTR, SLC26A9 likewise modifies treatment response (P = 0.02). Our findings indicate that SLC26A9 airway modification requires CFTR at the cell surface, and that a common variant in SLC26A9 may predict response to CFTR-directed therapeutics.

Deficiency of the ribosome biogenesis gene Sbds in hematopoietic stem and progenitor cells causes neutropenia in mice by attenuating lineage progression in myelocytes

Shwachman-Diamond syndrome is a congenital bone marrow failure disorder characterized by debilitating neutropenia. The disease is associated with loss-of-function mutations in the SBDS gene, implicated in ribosome biogenesis, but the cellular and molecular events driving cell specific phenotypes in ribosomopathies remain poorly defined. Here, we established what is to our knowledge the first mammalian model of neutropenia in Shwachman-Diamond syndrome through targeted downregulation of Sbds in hematopoietic stem and progenitor cells expressing the myeloid transcription factor CCAAT/enhancer binding protein α (Cebpa). Sbds deficiency in the myeloid lineage specifically affected myelocytes and their downstream progeny while, unexpectedly, it was well tolerated by rapidly cycling hematopoietic progenitor cells. Molecular insights provided by massive parallel sequencing supported cellular observations of impaired cell cycle exit and formation of secondary granules associated with the defect of myeloid lineage progression in myelocytes. Mechanistically, Sbds deficiency activated the p53 tumor suppressor pathway and induced apoptosis in these cells. Collectively, the data reveal a previously unanticipated, selective dependency of myelocytes and downstream progeny, but not rapidly cycling progenitors, on this ubiquitous ribosome biogenesis protein, thus providing a cellular basis for the understanding of myeloid lineage biased defects in Shwachman-Diamond syndrome.

A Joint Location-Scale Test Improves Power to Detect Associated SNPs, Gene Sets, and Pathways

Gene-based, pathway, and other multivariate association methods are motivated by the possibility of GxG and GxE interactions; however, accounting for such interactions is limited by the challenges associated with adequate modeling information. Here we propose an easy-to-implement joint location-scale (JLS) association testing framework for single-variant and multivariate analysis that accounts for interactions without explicitly modeling them. We apply the JLS method to a gene-set analysis of cystic fibrosis (CF) lung disease, which is influenced by multiple environmental and genetic factors. We identify and replicate an association between the constituents of the apical plasma membrane and CF lung disease (p = 0.0099 and p = 0.0180, respectively) and highlight a role for the SLC9A3-SLC9A3R1/2-EZR complex in contributing to CF lung disease. Many association studies could benefit from re-analysis with the JLS method that leverages complex genetic architecture for SNP, gene, and pathway identification. Analytical verification, simulation, and additional proof-of-principle applications support our approach.

In Vivo Senescence in the Sbds-Deficient Murine Pancreas: Cell-Type Specific Consequences of Translation Insufficiency

Genetic models of ribosome dysfunction show selective organ failure, highlighting a gap in our understanding of cell-type specific responses to translation insufficiency. Translation defects underlie a growing list of inherited and acquired cancer-predisposition syndromes referred to as ribosomopathies. We sought to identify molecular mechanisms underlying organ failure in a recessive ribosomopathy, with particular emphasis on the pancreas, an organ with a high and reiterative requirement for protein synthesis. Biallelic loss of function mutations in SBDS are associated with the ribosomopathy Shwachman-Diamond syndrome, which is typified by pancreatic dysfunction, bone marrow failure, skeletal abnormalities and neurological phenotypes. Targeted disruption of Sbds in the murine pancreas resulted in p53 stabilization early in the postnatal period, specifically in acinar cells. Decreased Myc expression was observed and atrophy of the adult SDS pancreas could be explained by the senescence of acinar cells, characterized by induction of Tgfβ, p15^Ink4b and components of the senescence-associated secretory program. This is the first report of senescence, a tumour suppression mechanism, in association with SDS or in response to a ribosomopathy. Genetic ablation of p53 largely resolved digestive enzyme synthesis and acinar compartment hypoplasia, but resulted in decreased cell size, a hallmark of decreased translation capacity. Moreover, p53 ablation resulted in expression of acinar dedifferentiation markers and extensive apoptosis. Our findings indicate a protective role for p53 and senescence in response to Sbds ablation in the pancreas. In contrast to the pancreas, the Tgfβ molecular signature was not detected in fetal bone marrow, liver or brain of mouse models with constitutive Sbds ablation. Nevertheless, as observed with the adult pancreas phenotype, disease phenotypes of embryonic tissues, including marked neuronal cell death due to apoptosis, were determined to be p53-dependent. Our findings therefore point to cell/tissue-specific responses to p53-activation that include distinction between apoptosis and senescence pathways, in the context of translation disruption.

Growth of all living things relies on protein synthesis. Failure of components of the complex protein synthesis machinery underlies a growing list of inherited and acquired multi—organ syndromes referred to as ribosomopathies. While ribosomes, the critical working components of the protein synthesis machinery, are required in all cell types to translate the genetic code, only certain organs manifest clinical symptoms in ribosomopathies, indicating specific cell-type features of protein synthesis control. Further, many of these diseases result in cancer despite an inherent deficit in growth. Here we report a range of consequences of protein synthesis insufficiency with loss of a broadly expressed ribosome factor, leading to growth impairment and cell cycle arrest at different stages. Apparent induction of p53-dependent cell death and arrest pathways included apoptosis in the fetal brain and senescence in the mature exocrine pancreas. The senescence, considered a tumour suppression mechanism, was accompanied by the expression of biomarkers associated with early stages of malignant transformation. These findings inform how cancer may initiate when growth is compromised and provide new insights into cell-type specific consequences of protein synthesis insufficiency.

Variants in Solute Carrier SLC26A9 Modify Prenatal Exocrine Pancreatic Damage in Cystic Fibrosis

OBJECTIVES

To test the hypothesis that multiple constituents of the apical plasma membrane residing alongside the causal cystic fibrosis (CF) transmembrane conductance regulator protein, including known CF modifiers SLC26A9, SLC6A14, and SLC9A3, would be associated with prenatal exocrine pancreatic damage as measured by newborn screened (NBS) immunoreactive trypsinogen (IRT) levels.

STUDY DESIGN

NBS IRT measures and genome-wide genotype data were available on 111 subjects from Colorado, 37 subjects from Wisconsin, and 80 subjects from France. Multiple linear regression was used to determine whether any of 8 single nucleotide polymorphisms (SNPs) in SLC26A9, SLC6A14, and SLC9A3 were associated with IRT and whether other constituents of the apical plasma membrane contributed to IRT.

RESULTS

In the Colorado sample, 3 SLC26A9 SNPs were associated with NBS IRT (min P=1.16×10(-3); rs7512462), but no SLC6A14 or SLC9A3 SNPs were associated (P>.05). The rs7512462 association replicated in the Wisconsin sample (P=.03) but not in the French sample (P=.76). Furthermore, rs7512462 was the top-ranked apical membrane constituent in the combined Colorado and Wisconsin sample.

CONCLUSIONS

NBS IRT is a biomarker of prenatal exocrine pancreatic disease in patients with CF, and a SNP in SLC26A9 accounts for significant IRT variability. This work suggests SLC26A9 as a potential therapeutic target to ameliorate exocrine pancreatic disease.

Evidence for a causal relationship between early exocrine pancreatic disease and cystic fibrosis-related diabetes: a Mendelian randomization study

Circulating immunoreactive trypsinogen (IRT), a biomarker of exocrine pancreatic disease in cystic fibrosis (CF), is elevated in most CF newborns. In those with severe CF transmembrane conductance regulator (CFTR) genotypes, IRT declines rapidly in the first years of life, reflecting progressive pancreatic damage. Consistent with this progression, a less elevated newborn IRT measure would reflect more severe pancreatic disease, including compromised islet compartments, and potentially increased risk of CF-related diabetes (CFRD). We show in two independent CF populations that a lower newborn IRT estimate is associated with higher CFRD risk among individuals with severe CFTR genotypes, and we provide evidence to support a causal relationship. Increased loge(IRT) at birth was associated with decreased CFRD risk in Canadian and Colorado samples (hazard ratio 0.30 [95% CI 0.15-0.61] and 0.39 [0.18-0.81], respectively). Using Mendelian randomization with the SLC26A9 rs7512462 genotype as an instrumental variable since it is known to be associated with IRT birth levels in the CF population, we provide evidence to support a causal contribution of exocrine pancreatic status on CFRD risk. Our findings suggest CFRD risk could be predicted in early life and that maintained ductal fluid flow in the exocrine pancreas could delay the onset of CFRD.

Genetic modifiers of cystic fibrosis-related diabetes

Diabetes is a common age-dependent complication of cystic fibrosis (CF) that is strongly influenced by modifier genes. We conducted a genome-wide association study in 3,059 individuals with CF (644 with CF-related diabetes [CFRD]) and identified single nucleotide polymorphisms (SNPs) within and 5' to the SLC26A9 gene that associated with CFRD (hazard ratio [HR] 1.38; P = 3.6 × 10(-8)). Replication was demonstrated in 694 individuals (124 with CFRD) (HR, 1.47; P = 0.007), with combined analysis significant at P = 9.8 × 10(-10). SLC26A9 is an epithelial chloride/bicarbonate channel that can interact with the CF transmembrane regulator (CFTR), the protein mutated in CF. We also hypothesized that common SNPs associated with type 2 diabetes also might affect risk for CFRD. A previous association of CFRD with SNPs in TCF7L2 was replicated in this study (P = 0.004; combined analysis P = 3.8 × 10(-6)), and type 2 diabetes SNPs at or near CDKAL1, CDKN2A/B, and IGF2BP2 were associated with CFRD (P < 0.004). These five loci accounted for 8.3% of the phenotypic variance in CFRD onset and had a combined population-attributable risk of 68%. Diabetes is a highly prevalent complication of CF, for which susceptibility is determined in part by variants at SLC26A9 (which mediates processes proximate to the CF disease-causing gene) and at four susceptibility loci for type 2 diabetes in the general population.

Breast cancer in a case of Shwachman Diamond syndrome

Shwachman Diamond syndrome (SDS) is a rare inherited bone marrow failure syndrome (IBMFS) characterized by neutropenia, exocrine pancreatic dysfunction, and cancer predisposition. Patients are at risk for myelodysplastic syndrome (MDS) and acute myelogenous leukemia (AML) but, unlike other IBMFS, there have been no reported cases of solid tumors. We report a novel case of a solid tumor in a patient with SDS and biallelic mutations in the Shwachman Bodian Diamond Syndrome gene (SBDS). Whether the development of breast cancer in this patient is due to SDS or an isolated case due to unknown factors requires further study.

Deficiency of Sbds in the mouse pancreas leads to features of Shwachman-Diamond syndrome, with loss of zymogen granules

BACKGROUND & AIMS

Shwachman-Diamond syndrome (SDS) is the second leading cause of hereditary exocrine pancreatic dysfunction. More than 90% of patients with SDS have biallelic loss-of-function mutations in the Shwachman-Bodian Diamond syndrome (SBDS) gene, which encodes a factor involved in ribosome function. We investigated whether mutations in Sbds lead to similar pancreatic defects in mice.

METHODS

Pancreas-specific knock-out mice were generated using a floxed Sbds allele and bred with mice carrying a null or disease-associated missense Sbds allele. Cre recombinase, regulated by the pancreatic transcription factor 1a promoter, was used to disrupt Sbds specifically in the pancreas. Models were assessed for pancreatic dysfunction and growth impairment.

RESULTS

Disruption of Sbds in the mouse pancreas was sufficient to recapitulate SDS phenotypes. Pancreata of mice with Sbds mutations had decreased mass, fat infiltration, but general preservation of ductal and endocrine compartments. Pancreatic extracts from mutant mice had defects in formation of the 80S ribosomal complex. The exocrine compartment of mutant mice was hypoplastic and individual acini produced few zymogen granules. The null Sbds allele resulted in an earlier onset of phenotypes as well as endocrine impairment. Mutant mice had reduced serum levels of digestive enzymes and overall growth impairment.

CONCLUSIONS

We developed a mouse model of SDS with pancreatic phenotypes similar to those of the human disease. This model could be used to investigate organ-specific consequences of Sbds-associated ribosomopathy. Sbds genotypes correlated with phenotypes. Defects developed specifically in the pancreata of mice, reducing growth of mice and production of digestive enzymes. SBDS therefore appears to be required for normal pancreatic development and function.

Multiple apical plasma membrane constituents are associated with susceptibility to meconium ileus in individuals with cystic fibrosis

Variants associated with meconium ileus in cystic fibrosis were identified in 3,763 affected individuals by genome-wide association study (GWAS). Five SNPs at two loci near SLC6A14 at Xq23-24 (minimum P = 1.28 × 10(-12) at rs3788766) and SLC26A9 at 1q32.1 (minimum P = 9.88 × 10(-9) at rs4077468) accounted for ~5% of phenotypic variability and were replicated in an independent sample of affected individuals (n = 2,372; P = 0.001 and 0.0001, respectively). By incorporating the knowledge that disease-causing mutations in CFTR alter electrolyte and fluid flux across surface epithelium into a hypothesis-driven GWAS (GWAS-HD), we identified associations with the same SNPs in SLC6A14 and SLC26A9 and established evidence for the involvement of SNPs in a third solute carrier gene, SLC9A3. In addition, GWAS-HD provided evidence of association between meconium ileus and multiple genes encoding constituents of the apical plasma membrane where CFTR resides (P = 0.0002; testing of 155 apical membrane genes jointly and in replication, P = 0.022). These findings suggest that modulating activities of apical membrane constituents could complement current therapeutic paradigms for cystic fibrosis.

A genome-wide linkage study of mammographic density, a risk factor for breast cancer

Introduction

Mammographic breast density is a highly heritable (h² > 0.6) and strong risk factor for breast cancer. We conducted a genome-wide linkage study to identify loci influencing mammographic breast density (MD).

Methods

Epidemiological data were assembled on 1,415 families from the Australia, Northern California and Ontario sites of the Breast Cancer Family Registry, and additional families recruited in Australia and Ontario. Families consisted of sister pairs with age-matched mammograms and data on factors known to influence MD. Single nucleotide polymorphism (SNP) genotyping was performed on 3,952 individuals using the Illumina Infinium 6K linkage panel.

Results

Using a variance components method, genome-wide linkage analysis was performed using quantitative traits obtained by adjusting MD measurements for known covariates. Our primary trait was formed by fitting a linear model to the square root of the percentage of the breast area that was dense (PMD), adjusting for age at mammogram, number of live births, menopausal status, weight, height, weight squared, and menopausal hormone therapy. The maximum logarithm of odds (LOD) score from the genome-wide scan was on chromosome 7p14.1-p13 (LOD = 2.69; 63.5 cM) for covariate-adjusted PMD, with a 1-LOD interval spanning 8.6 cM. A similar signal was seen for the covariate adjusted area of the breast that was dense (DA) phenotype. Simulations showed that the complete sample had adequate power to detect LOD scores of 3 or 3.5 for a locus accounting for 20% of phenotypic variance. A modest peak initially seen on chromosome 7q32.3-q34 increased in strength when only the 513 families with at least two sisters below 50 years of age were included in the analysis (LOD 3.2; 140.7 cM, 1-LOD interval spanning 9.6 cM). In a subgroup analysis, we also found a LOD score of 3.3 for DA phenotype on chromosome 12.11.22-q13.11 (60.8 cM, 1-LOD interval spanning 9.3 cM), overlapping a region identified in a previous study.

Conclusions

The suggestive peaks and the larger linkage signal seen in the subset of pedigrees with younger participants highlight regions of interest for further study to identify genes that determine MD, with the goal of understanding mammographic density and its involvement in susceptibility to breast cancer.

Genome-wide association and linkage identify modifier loci of lung disease severity in cystic fibrosis at 11p13 and 20q13.2

A combined genome-wide association and linkage study was used to identify loci causing variation in cystic fibrosis lung disease severity. We identified a significant association (P = 3.34 × 10(-8)) near EHF and APIP (chr11p13) in p.Phe508del homozygotes (n = 1,978). The association replicated in p.Phe508del homozygotes (P = 0.006) from a separate family based study (n = 557), with P = 1.49 × 10(-9) for the three-study joint meta-analysis. Linkage analysis of 486 sibling pairs from the family based study identified a significant quantitative trait locus on chromosome 20q13.2 (log(10) odds = 5.03). Our findings provide insight into the causes of variation in lung disease severity in cystic fibrosis and suggest new therapeutic targets for this life-limiting disorder.

Bone progenitor dysfunction induces myelodysplasia and secondary leukaemia

Mesenchymal cells contribute to the 'stroma' of most normal and malignant tissues, with specific mesenchymal cells participating in the regulatory niches of stem cells. By examining how mesenchymal osteolineage cells modulate haematopoiesis, here we show that deletion of Dicer1 specifically in mouse osteoprogenitors, but not in mature osteoblasts, disrupts the integrity of haematopoiesis. Myelodysplasia resulted and acute myelogenous leukaemia emerged that had acquired several genetic abnormalities while having intact Dicer1. Examining gene expression altered in osteoprogenitors as a result of Dicer1 deletion showed reduced expression of Sbds, the gene mutated in Schwachman-Bodian-Diamond syndrome-a human bone marrow failure and leukaemia pre-disposition condition. Deletion of Sbds in mouse osteoprogenitors induced bone marrow dysfunction with myelodysplasia. Therefore, perturbation of specific mesenchymal subsets of stromal cells can disorder differentiation, proliferation and apoptosis of heterologous cells, and disrupt tissue homeostasis. Furthermore, primary stromal dysfunction can result in secondary neoplastic disease, supporting the concept of niche-induced oncogenesis.

The behavioral phenotype of school-age children with shwachman diamond syndrome indicates neurocognitive dysfunction with loss of Shwachman-Bodian-Diamond syndrome gene function

OBJECTIVE

To investigate the cognitive, behavioral and adaptive functioning of children with Shwachman-Diamond syndrome (SDS).

STUDY DESIGN

Thirty-two children with SDS (6-17 years) were evaluated by use of standardized neuropsychological tests. Results were compared with normative data, unaffected siblings (n = 13), and age-and sex-matched children with cystic fibrosis (CF; n = 20).

RESULTS

Although intragroup variability in performance was evident, children with SDS displayed weaker overall intellectual reasoning, higher-order language skills, perceptual reasoning, visual-motor processing speed, visual motor- integration, visual executive problem-solving, attention, and aspects of academic achievement, as well as lower functional level of independence relative to the general population. Significant issues with behavior were also identified, including prior formal diagnoses and social problems. Lower abilities were found relative to sibling and CF control subjects and were not associated with secondary complications of SDS, age, or sex.

CONCLUSION

Neurocognitive deficits in subjects with SDS are largely independent of family environment and having a chronic illness and are likely the consequences of Shwachman-Bodian-Diamond syndrome gene dysfunction. There is a need for a broad-based approach to the assessment of cognitive function and appropriate remediation of individuals with SDS.

Shwachman-Bodian Diamond syndrome is a multi-functional protein implicated in cellular stress responses

Shwachman-Diamond syndrome (SDS; OMIM 260400) results from loss-of-function mutations in the Shwachman-Bodian Diamond syndrome (SBDS) gene. It is a multi-system disorder with clinical features of exocrine pancreatic dysfunction, skeletal abnormalities, bone marrow failure and predisposition to leukemic transformation. Although the cellular functions of SBDS are still unclear, its yeast ortholog has been implicated in ribosome biogenesis. Using affinity capture and mass spectrometry, we have developed an SBDS-interactome and report SBDS binding partners with diverse molecular functions, notably components of the large ribosomal subunit and proteins involved in DNA metabolism. Reciprocal co-immunoprecipitation confirmed the interaction of SBDS with the large ribosomal subunit protein RPL4 and with DNA-PK and RPA70, two proteins with critical roles in DNA repair. Function for SBDS in response to cellular stresses was implicated by demonstrating that SBDS-depleted HEK293 cells are hypersensitive to multiple types of DNA damage as well as chemically induced endoplasmic reticulum stress. Furthermore, using multiple routes to impair translation and mimic the effect of SBDS-depletion, we show that SBDS-dependent hypersensitivity of HEK293 cells to UV irradiation can be distinguished from a role of SBDS in translation. These results indicate functions of SBDS beyond ribosome biogenesis and may provide insight into the poorly understood cancer predisposition of SDS patients.

A systematic analysis of intronic sequences downstream of 5' splice sites reveals a widespread role for U-rich motifs and TIA1/TIAL1 proteins in alternative splicing regulation

To identify human intronic sequences associated with 5' splice site recognition, we performed a systematic search for motifs enriched in introns downstream of both constitutive and alternative cassette exons. Significant enrichment was observed for U-rich motifs within 100 nucleotides downstream of 5' splice sites of both classes of exons, with the highest enrichment between positions +6 and +30. Exons adjacent to U-rich intronic motifs contain lower frequencies of exonic splicing enhancers and higher frequencies of exonic splicing silencers, compared with exons not followed by U-rich intronic motifs. These findings motivated us to explore the possibility of a widespread role for U-rich motifs in promoting exon inclusion. Since cytotoxic granule-associated RNA binding protein (TIA1) and TIA1-like 1 (TIAL1; also known as TIAR) were previously shown in vitro to bind to U-rich motifs downstream of 5' splice sites, and to facilitate 5' splice site recognition in vitro and in vivo, we investigated whether these factors function more generally in the regulation of splicing of exons followed by U-rich intronic motifs. Simultaneous knockdown of TIA1 and TIAL1 resulted in increased skipping of 36/41 (88%) of alternatively spliced exons associated with U-rich motifs, but did not affect 32/33 (97%) alternatively spliced exons that are not associated with U-rich motifs. The increase in exon skipping correlated with the proximity of the first U-rich motif and the overall "U-richness" of the adjacent intronic region. The majority of the alternative splicing events regulated by TIA1/TIAL1 are conserved in mouse, and the corresponding genes are associated with diverse cellular functions. Based on our results, we estimate that approximately 15% of alternative cassette exons are regulated by TIA1/TIAL1 via U-rich intronic elements.

Exon skipping through the creation of a putative exonic splicing silencer as a consequence of the cystic fibrosis mutation R553X

Nonsense mutations that occur more than 50 bases upstream of terminal spliced junctions are generally thought to lead to degradation of the corresponding transcripts by the process of nonsense-mediated mRNA decay. It has also been proposed that some nonsense mutations may affect splicing by the process of nonsense-associated altered splicing (NAS), or by the disruption of a splicing regulatory element. In this study, the effect of the R553X mutation on the splicing of exon 11 of the cystic fibrosis transmembrane conductance regulator gene was investigated. Evidence that R553X causes exon 11 to skip through the creation of a putative exonic splicing silencer (ESS) was provided. The putative ESS appears to be active when located immediately upstream of a 5' splice site. These findings argue against the possibility that R553X-associated exon 11 skipping is caused by NAS. The study further suggests that aminoglycoside antibiotic treatment would not be effective for patients with the R553X mutation, owing to the skipping of exon 11, and further emphasises the need for detailed mechanistic characterisation of the consequences of nonsense disease mutations.

Loss of the mouse ortholog of the shwachman-diamond syndrome gene (Sbds) results in early embryonic lethality

Mutations in SBDS are responsible for Shwachman-Diamond syndrome (SDS), a disorder with clinical features of exocrine pancreatic insufficiency, bone marrow failure, and skeletal abnormalities. SBDS is a highly conserved protein whose function remains largely unknown. We identified and investigated the expression pattern of the murine ortholog. Variation in levels was observed, but Sbds was found to be expressed in all embryonic stages and most adult tissues. Higher expression levels were associated with rapid proliferation. A targeted disruption of Sbds was generated in order to understand the consequences of its loss in an in vivo model. Consistent with recessive disease inheritance for SDS, Sbds(+/-) mice have normal phenotypes, indistinguishable from those of their wild-type littermates. However, the development of Sbds(-/-) embryos arrests prior to embryonic day 6.5, with muted epiblast formation leading to early lethality. This finding is consistent with the absence of patients who are homozygous for early truncating mutations. Sbds is an essential gene for early mammalian development, with an expression pattern consistent with a critical role in cell proliferation.

Human homologs of Ubc6p ubiquitin-conjugating enzyme and phosphorylation of HsUbc6e in response to endoplasmic reticulum stress

Ubiquitin-conjugating enzyme Ubc6p is a tail-anchored protein that is localized to the cytoplasmic face of the endoplasmic reticulum (ER) membrane and has been implicated in the degradation of many misfolded membrane proteins in yeast. We have undertaken characterization studies of two human homologs, hsUbc6 and hsUbc6e. Both possess tail-anchored protein motifs, display high conservation in their catalytic domains, and are functional ubiquitin-conjugating enzymes as determined by in vitro thiol-ester assay. Both also display induction by the unfolded protein response, a feature of many ER-associated degradation (ERAD) components. Post-translational modification involving phosphorylation of hsUbc6e was observed to be ER-stress-related and dependent on signaling of the PRK-like ER kinase (PERK). The phosphorylation site was mapped to Ser-184, which resides within the uncharacterized region linking the highly conserved catalytic core and the C-terminal transmembrane domain. Phosphorylation of hsUbc6e also did not alter stability, subcellular localization, or interaction with a partner ubiquitin-protein isopeptide ligase. Assays of hsUbc6e(S184D) and hsUbc6e(S184E), which mimic the phosphorylated state, suggest that phosphorylation may reduce capacity for forming ubiquitin-enzyme thiol-esters. The occurrence of two distinct Ubc6p homologs in vertebrates, including one with phosphorylation modification in response to ER stress, emphasizes diversity in function between these Ub-conjugating enzymes during ERAD processes.

Phylogeny, sequence conservation, and functional complementation of the SBDS protein family

The Shwachman-Bodian-Diamond syndrome (SBDS) protein family occurs widely in nature, although its function has not been determined. Comprehensive database searches revealed SBDS homologues from 159 species, including examples from all sequenced archaeal and eukaryotic genomes and all eukaryotic kingdoms. Sequence alignment with ClustalX and MUSCLE algorithms led to the identification of conserved residues that occurred predominantly in the amino-terminal FYSH domain where they appeared to contribute to protein folding or stability. Only SBDS residue Gly91 was invariant in all species. Four distantly related protists were found to have two divergent SBDS genes in their genomes. In each case, phylogenetic analyses and the identification of shared sequence features suggested that one gene was derived from lateral gene transfer. We also identified a shared C-terminal zinc finger domain fusion in flowering plants and chromalveolates that may shed light on the function of the protein family and the evolutionary histories of these kingdoms. To assess the extent of SBDS functional conservation, we carried out complementation studies of SBDS homologues and interspecies chimeras in Saccharomyces cerevisiae. We determined that the FYSH domain was widely interchangeable among eukaryotes, while domain 2 imparted species specificity to protein function. Domain 3 was largely dispensable for function in our yeast complementation assay. Overall, the phylogeny of SBDS was shared with a group of proteins that were markedly enriched for RNA metabolism and/or ribosome-associated functions. These findings link Shwachman-Diamond syndrome to other bone marrow failure syndromes with defects in nucleolus-associated processes, including Diamond-Blackfan anemia, cartilage-hair hypoplasia, and dyskeratosis congenita.

Identification of the gene responsible for methylmalonic aciduria and homocystinuria, cblC type

Methylmalonic aciduria and homocystinuria, cblC type (OMIM 277400), is the most common inborn error of vitamin B(12) (cobalamin) metabolism, with about 250 known cases. Affected individuals have developmental, hematological, neurological, metabolic, ophthalmologic and dermatologic clinical findings. Although considered a disease of infancy or childhood, some individuals develop symptoms in adulthood. The cblC locus was mapped to chromosome region 1p by linkage analysis. We refined the chromosomal interval using homozygosity mapping and haplotype analyses and identified the MMACHC gene. In 204 individuals, 42 different mutations were identified, many consistent with a loss of function of the protein product. One mutation, 271dupA, accounted for 40% of all disease alleles. Transduction of wild-type MMACHC into immortalized cblC fibroblast cell lines corrected the cellular phenotype. Molecular modeling predicts that the C-terminal region of the gene product folds similarly to TonB, a bacterial protein involved in energy transduction for cobalamin uptake.

Mammographic breast density as an intermediate phenotype for breast cancer

The amount of radiologically dense breast-tissue appearing on a mammogram varies between women because of differences in the composition of breast tissue, and is referred to here as mammographic density. This review presents evidence that mammographic density is a strong risk factor for breast cancer, and that risk of breast cancer is four to five times greater in women with density in more than 75% of the breast than in women with little or no density in the breast. Density in more than 50% of the breast could account for about a third of breast cancers. The epidemiology of mammographic density is consistent with its being a marker of susceptibility to breast cancer. Twin studies have shown that the proportion of the breast occupied by density, at a given age, is highly heritable, and inherited factors explain 63% of the variance. Mammographic breast density has the characteristics of a quantitative trait and might be determined by genes that are easier to identify than those for breast cancer itself. The genes that determine breast density might also be associated with risk of breast cancer, and their identification is also likely to provide insights into the biology of the breast and identify potential targets for preventive strategies.

The Shwachman-Bodian-Diamond Syndrome Protein Family Is Involved in RNA Metabolism

A combination of structural, biochemical, and genetic studies in model organisms was used to infer a cellular role for the human protein (SBDS) responsible for Shwachman-Bodian-Diamond syndrome. The crystal structure of the SBDS homologue in Archaeoglobus fulgidus, AF0491, revealed a three domain protein. The N-terminal domain, which harbors the majority of disease-linked mutations, has a novel three-dimensional fold. The central domain has the common winged helix-turn-helix motif, and the C-terminal domain shares structural homology with known RNA-binding domains. Proteomic analysis of the SBDS sequence homologue in Saccharomyces cerevisiae, YLR022C, revealed an association with over 20 proteins involved in ribosome biosynthesis. NMR structural genomics revealed another yeast protein, YHR087W, to be a structural homologue of the AF0491 N-terminal domain. Sequence analysis confirmed them as distant sequence homologues, therefore related by divergent evolution. Synthetic genetic array analysis of YHR087W revealed genetic interactions with proteins involved in RNA and rRNA processing including Mdm20/Nat3, Nsr1, and Npl3. Our observations, taken together with previous reports, support the conclusion that SBDS and its homologues play a role in RNA metabolism.

Structural and Mutational Analysis of the SBDS Protein Family

Shwachman-Diamond Syndrome (SDS) is an autosomal recessive disorder characterized by bone marrow failure with significant predisposition to the development of poor prognosis myelodysplasia and leukemia, exocrine pancreatic failure and metaphyseal chondrodysplasia. Although the SBDS gene mutated in this disorder is highly conserved in Archaea and all eukaryotes, the function is unknown. To interpret the molecular consequences of SDS-associated mutations, we have solved the crystal structure of the Archaeoglobus fulgidus SBDS protein orthologue at a resolution of 1.9 angstroms, revealing a three domain architecture. The N-terminal (FYSH) domain is the most frequent target for disease mutations and contains a novel mixed alpha/beta-fold identical to the single domain yeast protein Yhr087wp that is implicated in RNA metabolism. The central domain consists of a three-helical bundle, whereas the C-terminal domain has a ferredoxin-like fold. By genetic complementation analysis of the essential Saccharomyces cerevisiae SBDS orthologue YLR022C, we demonstrate an essential role in vivo for the FYSH domain and the central three-helical bundle. We further show that the common SDS-related K62X truncation is non-functional. Most SDS-related missense mutations that alter surface epitopes do not impair YLR022C function, but mutations affecting residues buried in the hydrophobic core of the FYSH domain severely impair or abrogate complementation. These data are consistent with absence of homozygosity for the common K62X truncation mutation in individuals with SDS, indicating that the SDS disease phenotype is a consequence of expression of hypomorphic SBDS alleles and that complete loss of SBDS function is likely to be lethal.

Skeletal phenotype in patients with Shwachman-Diamond syndrome and mutations in SBDS

Pancreatic exocrine and bone marrow dysfunctions are considered to be universal features of Shwachman-Diamond syndrome (SDS) whereas the associated skeletal dysplasia is variable and not consistently observed. The genetic defect in SDS has recently been identified; causative mutations have been shown in the SBDS gene. The aims of this study were to characterize the nature, frequency, and age-related changes of radiographic skeletal abnormalities in patients with SBDS mutations and to assess genotype-phenotype correlation. Fifteen patients (mean age 9.7 years) with a clinical diagnosis of SDS and documented SBDS gene mutations were included. Review of their skeletal radiographs showed abnormalities in all patients. The skeletal changes were variable, even in patients with identical genotypes. The typical features were (1) delayed appearance of secondary ossification centers, (2) variable widening and irregularity of the metaphyses in early childhood, followed by progressive thickening and irregularity of the growth plates, and (3) generalized osteopenia. There was a tendency towards normalization of the epiphyseal maturation defect and progression of the metaphyseal changes with age. The results suggest that the characteristic skeletal changes are present in all patients with SDS and SBDS mutations, but their severity and localization varies with age. No phenotype-genotype correlation was observed.

Human chromosome 7: DNA sequence and biology

DNA sequence and annotation of the entire human chromosome 7, encompassing nearly 158 million nucleotides of DNA and 1917 gene structures, are presented. To generate a higher order description, additional structural features such as imprinted genes, fragile sites, and segmental duplications were integrated at the level of the DNA sequence with medical genetic data, including 440 chromosome rearrangement breakpoints associated with disease. This approach enabled the discovery of candidate genes for developmental diseases including autism.

Mutations in SBDS are associated with Shwachman-Diamond syndrome

Shwachman-Diamond syndrome (SDS; OMIM 260400) is an autosomal recessive disorder with clinical features that include pancreatic exocrine insufficiency, hematological dysfunction and skeletal abnormalities. Here, we report identification of disease-associated mutations in an uncharacterized gene, SBDS, in the interval of 1.9 cM at 7q11 previously shown to be associated with the disease. We report that SBDS has a 1.6-kb transcript and encodes a predicted protein of 250 amino acids. A pseudogene copy (SBDSP) with 97% nucleotide sequence identity resides in a locally duplicated genomic segment of 305 kb. We found recurring mutations resulting from gene conversion in 89% of unrelated individuals with SDS (141 of 158), with 60% (95 of 158) carrying two converted alleles. Converted segments consistently included at least one of two pseudogene-like sequence changes that result in protein truncation. SDBS is a member of a highly conserved protein family of unknown function with putative orthologs in diverse species including archaea and eukaryotes. Archaeal orthologs are located within highly conserved operons that include homologs of RNA-processing genes, suggesting that SDS may be caused by a deficiency in an aspect of RNA metabolism that is essential for development of the exocrine pancreas, hematopoiesis and chrondrogenesis.

Serum pancreatic enzymes define the pancreatic phenotype in patients with Shwachman-Diamond syndrome

OBJECTIVE

To evaluate the role of serum enzymes for defining the pancreatic phenotype in Shwachman-Diamond syndrome (SDS), an inherited multisystem condition.

STUDY DESIGN

Serum pancreatic trypsinogen and isoamylase were measured in 164 patients known or presumed to have SDS. The diagnosis was confirmed in 90 patients. Among 74 unconfirmed cases, 35 ("probable SDS") had hematologic dysfunction but lacked documented pancreatic dysfunction, whereas 39 patients ("improbable SDS") lacked both documented pancreatic and hematologic dysfunction. Classification and regression tree (CART) analysis was performed in 90 patients with SDS and 134 control patients to establish a rule for defining the pancreatic phenotype of SDS; the rule was then applied to the patients with unconfirmed diagnosis.

RESULTS

In the control patients, serum trypsinogen showed little variation with age, whereas serum isoamylase values rose from birth on, attaining adult values by 3 years. For patients with SDS, serum trypsinogen values were low in young patients and tended to increase with age, whereas serum isoamylase values remained low at all ages. The CART rule combined results from both enzymes and classified the pancreatic phenotype in all but one SDS patient, who was <3 years of age. Excluding patients <3 years of age, CART identified the pancreatic phenotype in 82% and 7% of the "probable SDS" and "improbable SDS" cases, respectively.

CONCLUSIONS

Serum pancreatic enzymes are useful for determining the pancreatic phenotype and confirming the diagnosis of SDS.

Non-CFTR chloride channels likely contribute to secretion in the murine small intestine

While most cystic fibrosis (CF) transmembrane conductance regulator (CFTR)-knockout animals die due to intestinal obstruction before or at the time of weaning, a subpopulation of these animals are long living and exhibit a milder phenotype. The decreased severity of intestinal disease in these mildly affected CF mice is related to the expression of non-CFTR genetic modifiers. The identity of these genetic modifiers is not known, but we hypothesize that they may complement CFTR function as a chloride channel in this tissue. To assess the contribution of non-CFTR chloride channels to chloride secretion across the small intestine of CF mice with mild disease, we measured the basal transepithelial potential difference across this tissue as well as the secretory response to agonists of the cAMP and the calcium-mediated signaling pathways. Chloride secretion across the small intestine of mildly affected CF mice was not stimulated by forskolin or by carbachol. The absence of CFTR is thus not compensated by the activity of a distinct, cAMP- or calcium-activated chloride channel at the apical surface of the intestinal epithelium. On the other hand, a basal chloride secretion across the intestinal epithelium was present in these animals, and we hypothesize that this activity may be linked to improved survival of these animals.

Fine mapping of the locus for Shwachman-Diamond syndrome at 7q11, identification of shared disease haplotypes, and exclusion of TPST1 as a candidate gene

Shwachman-Diamond syndrome (SDS) is an autosomal recessive disorder characterised by exocrine pancreatic dysfunction, haematological and skeletal abnormalities. We have previously defined the SDS locus as a 2.7 cM interval spanning the centromere of chromosome 7. To facilitate additional analysis of this complex and poorly characterised region, a framework of ordered genetic markers at 7p11-q11, including six newly identified, has been constructed using somatic cell, radiation hybrid and STS-content mapping. We have identified shared disease haplotypes, that recur in unrelated families of common ethnic origin, and extend across the SDS locus. Detection of ancestral and intrafamilial recombination events in patients refined the SDS locus to a 1.9 cM interval at 7q11, which contains the tyrosylprotein sulfotransferase 1 (TPST1) gene. Patients with SDS were screened for mutations in TPST1 by sequencing of exons and intron-exon junctions. Two single nucleotide polymorphisms, but no disease-causing mutations, were identified. In addition, Southern blot analysis yielded no evidence of large-scale mutations, and RT-PCR analysis failed to detect alterations in expression. These results exclude TPST1 as the causative gene for SDS. The established map of the refined SDS locus will assist in the identification and characterisation of other candidate genes for SDS.

Shwachman-Diamond syndrome with exocrine pancreatic dysfunction and bone marrow failure maps to the centromeric region of chromosome 7

Shwachman-Diamond syndrome (SDS) is an autosomal recessive disorder characterized by exocrine pancreatic insufficiency and hematologic and skeletal abnormalities. A genomewide scan of families with SDS was terminated at approximately 50% completion, with the identification of chromosome 7 markers that showed linkage with the disease. Finer mapping revealed significant linkage across a broad interval that included the centromere. The maximum two-point LOD score was 8.7, with D7S473, at a recombination fraction of 0. The maximum multipoint LOD score was 10, in the interval between D7S499 and D7S482 (5.4 cM on the female map and 0 cM on the male map), a region delimited by recombinant events detected in affected children. Evidence from all 15 of the multiplex families analyzed provided support for the linkage, consistent with a single locus for SDS. However, the presence of several different mutations is suggested by the heterogeneity of disease-associated haplotypes in the candidate region.

Determinants of the nuclear localization of the heterodimeric DNA fragmentation factor (ICAD/CAD)

Programmed cell death or apoptosis leads to the activation of the caspase-activated DNase (CAD), which degrades chromosomal DNA into nucleosomal fragments. Biochemical studies revealed that CAD forms an inactive heterodimer with the inhibitor of caspase-activated DNase (ICAD), or its alternatively spliced variant, ICAD-S, in the cytoplasm. It was initially proposed that proteolytic cleavage of ICAD by activated caspases causes the dissociation of the ICAD/CAD heterodimer and the translocation of active CAD into the nucleus in apoptotic cells. Here, we show that endogenous and heterologously expressed ICAD and CAD reside predominantly in the nucleus in nonapoptotic cells. Deletional mutagenesis and GFP fusion proteins identified a bipartite nuclear localization signal (NLS) in ICAD and verified the function of the NLS in CAD. The two NLSs have an additive effect on the nuclear targeting of the CAD-ICAD complex, whereas ICAD-S, lacking its NLS, appears to have a modulatory role in the nuclear localization of CAD. Staurosporine-induced apoptosis evoked the proteolysis and disappearance of endogenous and exogenous ICAD from the nuclei of HeLa cells, as monitored by immunoblotting and immunofluorescence microscopy. Similar phenomenon was observed in the caspase-3-deficient MCF7 cells upon expressing procaspase-3 transiently. We conclude that a complex mechanism, involving the recognition of the NLSs of both ICAD and CAD, accounts for the constitutive accumulation of CAD/ICAD in the nucleus, where caspase-3-dependent regulation of CAD activity takes place.

Cystic fibrosis mutations lead to carboxyl-terminal fragments that highlight an early biogenesis step of the cystic fibrosis transmembrane conductance regulator

Inefficient delivery of the cystic fibrosis transmembrane conductance regulator (CFTR) to the surface of cells contributes to disease in the majority of cystic fibrosis patients. Analysis of cystic fibrosis-associated missense mutations in the first nucleotide binding domain (NBD1), including A455E, S549R, Y563N, and P574H, revealed reduced levels of mature CFTR with elevated levels of carboxyl-terminal polypeptide fragments of 105 and 90 kDa. These fragments appear early in biogenesis and degrade rapidly in four distinct cell types tested including the bronchial epithelial IB3-1 cell line. They were detected at highest levels with CFTRA455E where the 105-kDa fragment accounted for 40% of newly synthesized polypeptide but for only 20 and 7% of nascent wild type and mutant DeltaF508 proteins, respectively. The bands represent core- and unglycosylated forms of the same CFTR fragment supporting that precursor forms are correctly inserted into the membrane of the endoplasmic reticulum. Proteolytic cleavage would be predicted to occur on the cytosolic face of the endoplasmic reticulum within the NBD1-R domain segment, but pharmacological testing did not support involvement of the 26 S proteasome. The examined missense mutations in NBD1 manifest differently than the major mutant, DeltaF508, and highlight a critical conformational aspect of biogenesis of CFTR.

Segregation analysis in Shwachman-Diamond syndrome: evidence for recessive inheritance

Shwachman-Diamond syndrome is a rare disorder of unknown cause. Reports have indicated the occurrence of affected siblings, but formal segregation analysis has not been performed. In families collected for genetic studies, the mean paternal age and mean difference in parental ages were found to be consistent with the general population. We determined estimates of segregation proportion in a cohort of 84 patients with complete sibship data under the assumption of complete ascertainment, using the Li and Mantel estimator, and of single ascertainment with the Davie modification. A third estimate was also computed with the expectation-maximization (EM) algorithm. All three estimates supported an autosomal recessive mode of inheritance, but complete ascertainment was found to be unlikely. Although there are no overt signs of disease in adult carriers (parents), the use of serum trypsinogen levels to indicate exocrine pancreatic dysfunction was evaluated as a potential measure for heterozygote expression. No consistent differences were found in levels between parents and a normal control population. Although genetic heterogeneity cannot be excluded, our results indicate that simulation and genetic analyses of Shwachman-Diamond syndrome should consider a recessive model of inheritance.

Exclusion of linkage of Shwachman-Diamond syndrome to chromosome regions 6q and 12q implicated by a de novo translocation

Shwachman-Diamond syndrome is a rare genetic disorder of unknown pathogenesis involving exocrine pancreatic insufficiency and hematological and skeletal abnormalities. There is broad clinical variability; the extent of heterogeneity is unknown but comparisons within a large cohort of patients show no striking differences between patients of families with single or multiple affected offspring. Segregation analysis of a cohort of 69 families has suggested an autosomal recessive mode of inheritance. A single constitutional de novo chromosome rearrangement was reported in a Japanese patient involving a balanced translocation, t(6;12)(q16.2;q21.2), thereby suggesting possible loci for a genetic defect. Evenly spaced microsatellite markers spanning 26-32 cM intervals from D6S1056 to D6S304 and D12S375 to D12S346 were analyzed for linkage in members of 13 Shwachman-Diamond syndrome families with two or three affected children. Two-point lod scores were calculated for each marker under assumptions of recessive inheritance and complete penetrance. Negative lod scores indicated exclusion of both chromosome regions. Further, affected sibs were discordant for inheritance of chromosomes in most families based on constructed haplotypes. The cytogenetic abnormality is not associated with most cases of Shwachman-Diamond syndrome.

Shwachman syndrome: phenotypic manifestations of sibling sets and isolated cases in a large patient cohort are similar

OBJECTIVES

With the use of clinical data from a large international cohort, we evaluated and compared affected siblings and isolated cases.

STUDY DESIGN

Data from 116 families were collected, and patients conforming to our predetermined diagnostic criteria were analyzed. Phenotypic manifestations of affected siblings and singletons were compared with the use of t tests, Wilcoxon scores, and chi2 analysis.

RESULTS

Eighty-eight patients (33 female, 55 male; median age 5.20 years) fulfilled our predetermined diagnostic criteria for Shwachman syndrome; 63 patients were isolated cases, and 25 affected siblings were from 12 multiplex families. Steatorrhea was present in 86% (57 of 66), and 91% (78 of 86) displayed a low serum trypsinogen concentration. Patients older than 4 years more often had pancreatic sufficiency. Neutropenia occurred in 98%, anemia in 42%, and thrombocytopenia in 34%. Myelodysplasia or cytogenetic abnormalities were reported in 7 patients. Short stature with normal nutritional status was a prominent feature.

CONCLUSIONS

Clinical features among patients with Shwachman syndrome varied between patients and with age. Similarities in phenotype between isolated cases and affected sibling sets support the hypothesis that Shwachman syndrome is a single disease entity.

A physical and transcriptional map of the preaxial polydactyly locus on chromosome 7q36

Preaxial polydactyly is a congenital hand malformation that includes duplicated thumbs, various forms of triphalangeal thumbs, and duplications of the index finger. A locus for preaxial polydactyly has been mapped to a region of 1.9 cM on chromosome 7q36 between polymorphic markers D7S550 and D7S2423. We constructed a detailed physical map of the preaxial polydactyly candidate region. With a combination of methods we identified and positioned 11 transcripts within this map. By recombination analysis on families with preaxial polydactyly, using newly developed polymorphic markers, we were able to reduce the candidate region to approximately 450 kb. The homeobox gene HLXB9, a putative receptor C7orf2, and two transcripts of unknown function, C7orf3 and C7orf4, map in the refined candidate region and have been subjected to mutation analysis in individuals with preaxial polydactyly.