Proteome-Wide Assessment of Clustering of Missense Variants in Neurodevelopmental Disorders Versus Cancer

Missense de novo variants (DNVs) and missense somatic variants contribute to neurodevelopmental disorders (NDDs) and cancer, respectively. Proteins with statistical enrichment based on analyses of these variants exhibit convergence in the differing NDD and cancer phenotypes. Herein, the question of why some of the same proteins are identified in both phenotypes is examined through investigation of clustering of missense variation at the protein level. Our hypothesis is that missense variation is present in different protein locations in the two phenotypes leading to the distinct phenotypic outcomes. We tested this hypothesis in 1D protein space using our software CLUMP. Furthermore, we newly developed 3D-CLUMP that uses 3D protein structures to spatially test clustering of missense variation for proteome-wide significance. We examined missense DNVs in 39,883 parent-child sequenced trios with NDDs and missense somatic variants from 10,543 sequenced tumors covering five TCGA cancer types and two COSMIC pan-cancer aggregates of tissue types. There were 57 proteins with proteome-wide significant missense variation clustering in NDDs when compared to cancers and 79 proteins with proteome-wide significant missense clustering in cancers compared to NDDs. While our main objective was to identify differences in patterns of missense variation, we also identified a novel NDD protein BLTP2. Overall, our study is innovative, provides new insights into differential missense variation in NDDs and cancer at the protein-level, and contributes necessary information toward building a framework for thinking about prognostic and therapeutic aspects of these proteins.

High Coverage Highly Accurate Long-Read Sequencing of a Mouse Neuronal Cell Line Using the PacBio Revio Sequencer

Recently, Pacific Biosciences released a new highly accurate long-read sequencer called the Revio System that is projected to generate 30× HiFi whole-genome sequencing for the human genome within one sequencing SMRT Cell. Mouse and human genomes are similar in size. In this study, we sought to test this new sequencer by characterizing the genome and epigenome of the mouse neuronal cell line Neuro-2a. We generated long-read HiFi whole-genome sequencing on three Revio SMRT Cells, achieving a total coverage of 98×, with 30×, 32×, and 36× coverage respectively for each of the three Revio SMRT Cells. We performed several tests on these data including single-nucleotide variant and small insertion detection using GPU-accelerated DeepVariant, structural variant detection with pbsv, methylation detection with pb-CpG-tools, and generating de novo assemblies with the HiCanu and hifiasm assemblers. Overall, we find consistency across SMRT Cells in coverage, detection of variation, methylation, and de novo assemblies for each of the three SMRT Cells.

Biologic characterization of ABCA3 variants in lung tissue from infants and children with ABCA3 deficiency

ABCA3 is a phospholipid transporter protein required for surfactant assembly in lamellar bodies of alveolar type II cells. Biallelic pathogenic ABCA3 variants cause severe neonatal respiratory distress syndrome or childhood interstitial lung disease. However, ABCA3 genotype alone does not explain the diversity in disease presentation, severity, and progression. Additionally, monoallelic ABCA3 variants have been reported in infants and children with ABCA3-deficient phenotypes. The effects of most ABCA3 variants identified in patients have not been characterized at the RNA level. ABCA3 allele-specific expression occurs in some cell types due to epigenetic regulation. We obtained lung tissue at transplant or autopsy from 16 infants and children with ABCA3 deficiency due to compound heterozygous ABCA3 variants for biologic characterization of the predicted effects of ABCA3 variants at the RNA level and determination of ABCA3 allele expression. We extracted DNA and RNA from frozen lung tissue and reverse-transcribed cDNA from mRNA. We performed Sanger sequencing to assess allele-specific expression by comparing the heights of variant nucleotide peaks in amplicons from genomic DNA and cDNA. We found similar genomic and cDNA variant nucleotide peak heights and no evidence of allele-specific expression among explant or autopsy samples with biallelic missense ABCA3 variants (n = 6). We observed allele-specific expression of missense alleles in trans with frameshift (n = 4) or nonsense (n = 1) variants, attributable to nonsense-mediated decay. The missense variant c.53 A > G;p.Gln18Arg, located near an exon-intron junction, encoded abnormal splicing with skipping of exon 4. Biologic characterization of ABCA3 variants can inform discovery of variant-specific disease mechanisms.

Functional Genomics of ABCA3 Variants

Rare or private, biallelic variants in the ABCA3 (ATP-binding cassette transporter A3) gene are the most common monogenic cause of lethal neonatal respiratory failure and childhood interstitial lung disease. Functional characterization of fewer than 10% of over 200 disease-associated ABCA3 variants (majority missense) suggests either disruption of ABCA3 protein trafficking (type I) or of ATPase-mediated phospholipid transport (type II). Therapies remain limited and nonspecific. A scalable platform is required for functional characterization of ABCA3 variants and discovery of pharmacologic correctors. To address this need, we first silenced the endogenous ABCA3 locus in A549 cells with CRISPR/Cas9 genome editing. Next, to generate a parent cell line (A549/ABCA3^-/-) with a single recombination target site for genomic integration and stable expression of individual ABCA3 missense variant cDNAs, we used lentiviral-mediated integration of a LoxFAS cassette, FACS, and dilutional cloning. To assess the fidelity of this cell-based model, we compared functional characterization (ABCA3 protein processing, ABCA3 immunofluorescence colocalization with intracellular markers, ultrastructural vesicle phenotype) of two individual ABCA3 mutants (type I mutant, p.L101P; type II mutant, p.E292V) in A549/ABCA3^-/- cells and in both A549 cells and primary, human alveolar type II cells that transiently express each cDNA after adenoviral-mediated transduction. We also confirmed pharmacologic rescue of ABCA3 variant-encoded mistrafficking and vesicle diameter in A549/ABCA3^-/- cells that express p.G1421R (type I mutant). A549/ABCA3^-/- cells provide a scalable, genetically versatile, physiologically relevant functional genomics platform for discovery of variant-specific mechanisms that disrupt ABCA3 function and for screening of potential ABCA3 pharmacologic correctors.

Functional characterization of four ATP-binding cassette transporter A3 gene (ABCA3) variants

ABCA3 transports phospholipids across lamellar body membranes in pulmonary alveolar type II cells and is required for surfactant assembly. Rare, biallelic, pathogenic ABCA3 variants result in lethal neonatal respiratory distress syndrome and childhood interstitial lung disease. Qualitative functional characterization of ABCA3 missense variants suggests two pathogenic classes: disrupted intracellular trafficking (type I mutant) or impaired ATPase-mediated phospholipid transport into the lamellar bodies (type II mutant). We qualitatively compared wild-type (WT-ABCA3) with four uncharacterized ABCA3 variants (c.418A>C;p.Asn140His, c.3609_3611delCTT;p.Phe1203del, c.3784A>G;p.Ser1262Gly, and c.4195G>A;p.Val1399Met) in A549 cells using protein processing, colocalization with intracellular organelles, lamellar body ultrastructure, and ATPase activity. We quantitatively measured lamellar body-like vesicle diameter and intracellular ABCA3 trafficking using fluorescence-based colocalization. Three ABCA3 variants (p.Asn140His, p.Ser1262Gly, and p.Val1399Met) were processed and trafficked normally and demonstrated well-organized lamellar body-like vesicles, but had reduced ATPase activity consistent with type II mutants. P.Phe1203del was processed normally, had reduced ATPase activity, and well-organized lamellar body-like vesicles, but quantitatively colocalized with both endoplasmic reticulum and lysosomal markers, an intermediate phenotype suggesting disruption of both intracellular trafficking and phospholipid transport. All ABCA3 mutants demonstrated mean vesicle diameters smaller than WT-ABCA3. Qualitative and quantitative functional characterization of ABCA3 variants informs mechanisms of pathogenicity.

Functional Characterization of ATP-Binding Cassette Transporter A3 Mutations from Infants with Respiratory Distress Syndrome

Mutations in the ATP-binding cassette transporter A3 gene (ABCA3) result in severe neonatal respiratory distress syndrome and childhood interstitial lung disease. As most ABCA3 mutations are rare or private, determination of mutation pathogenicity is often based on results from in silico prediction tools, identification in unrelated diseased individuals, statistical association studies, or expert opinion. Functional biologic studies of ABCA3 mutations are needed to confirm mutation pathogenicity and inform clinical decision making. Our objective was to functionally characterize two ABCA3 mutations (p.R288K and p.R1474W) identified among term and late-preterm infants with respiratory distress syndrome with unclear pathogenicity in a genetically versatile model system. We performed transient transfection of HEK293T cells with wild-type or mutant ABCA3 alleles to assess protein processing with immunoblotting. We used transduction of A549 cells with adenoviral vectors, which concurrently silenced endogenous ABCA3 and expressed either wild-type or mutant ABCA3 alleles (p.R288K and p.R1474W) to assess immunofluorescent localization, ATPase activity, and organelle ultrastructure. Both ABCA3 mutations (p.R288K and p.R1474W) encoded proteins with reduced ATPase activity but with normal intracellular localization and protein processing. Ultrastructural phenotypes of lamellar body-like vesicles in A549 cells transduced with mutant alleles were similar to wild type. Mutant proteins encoded by ABCA3 mutations p.R288K and p.R1474W had reduced ATPase activity, a biologically plausible explanation for disruption of surfactant metabolism by impaired phospholipid transport into the lamellar body. These results also demonstrate the usefulness of a genetically versatile, human model system for functional characterization of ABCA3 mutations with unclear pathogenicity.

Measurement of human surfactant protein-B turnover in vivo from tracheal aspirates using targeted proteomics

We describe a method to measure protein synthesis and catabolism in humans without prior purification and use the method to measure the turnover of surfactant protein-B (SP-B). SP-B, a lung-specific, hydrophobic protein essential for fetal-neonatal respiratory transition, is present in only picomolar quantities in tracheal aspirate samples and difficult to isolate for dynamic turnover studies using traditional in vivo tracer techniques. Using infusion of [5,5,5-(2)H(3)] leucine and a targeted proteomics method, we measured both the quantity and kinetics of SP-B tryptic peptides in tracheal aspirate samples of symptomatic newborn infants. The fractional synthetic rate (FSR) of SP-B measured using the most abundant proteolytic fragment, a 10 amino acid peptide from the carboxy-terminus of proSP-B (SPTGEWLPR), from the circulating leucine pool was 0.035 +/- 0.005 h(-1), and the fractional catabolic rate was 0.044 +/- 0.003 h(-1). This technique permits high-throughput and sensitive measurement of turnover of low abundance proteins with minimal sample preparation.

Developmental and genetic regulation of human surfactant protein B in vivo

BACKGROUND

Genetic and developmental disruption of surfactant protein B (SP-B) expression causes neonatal respiratory distress syndrome (RDS).

OBJECTIVES

To assess developmental and genetic regulation of SP-B expression in vivo.

METHODS

To evaluate in vivo developmental regulation of SP-B, we used immunoblotting to compare frequency of detection of mature and pro-SP-B peptides in developmentally distinct cohorts: 24 amniotic fluid samples, unfractionated tracheal aspirates from 101 infants >or=34 weeks' gestation with (75) and without (26) neonatal RDS, and 6 nonsmoking adults. To examine genetic regulation, we used univariate and logistic regression analyses to detect associations between common SP-B (SFTPB) genotypes and SP-B peptides in the neonatal RDS cohort.

RESULTS

We found pro-SP-B peptides in 24/24 amniotic fluid samples and in 100/101 tracheal aspirates from newborn infants but none in bronchoalveolar lavage from normal adults (0/6) (p < 0.001). We detected an association (p = 0.0011) between pro-SP-B peptides (M(r) 40 and 42 kDa) and genotype of a nonsynonymous single nucleotide polymorphism at genomic position 1580 that regulates amino-terminus glycosylation.

CONCLUSIONS

Pro-SP-B peptides are more common in developmentally less mature humans. Association of genotype at genomic position 1580 with pro-SP-B peptides (M(r) 40 and 42 kDa) suggests genetic regulation of amino terminus glycosylation in vivo.

Population and disease-based prevalence of the common mutations associated with surfactant deficiency

The prevalence of the common mutations in the surfactant protein-B (121ins2), surfactant protein-C (I73T), and ATP-binding cassette member A3 (E292V) genes in population-based or case-control cohorts of newborn respiratory distress syndrome (RDS) is unknown. We determined the frequencies of these mutations in ethnically diverse population and disease-based cohorts using restriction enzyme analysis (121ins2 and E292V) and a 5' nuclease assay (I73T) in DNA samples from population-based cohorts in Missouri, Norway, South Korea, and South Africa, and from a case-control cohort of newborns with and without RDS (n = 420). We resequenced the ATP-binding cassette member A3 gene (ABCA3) in E292V carriers and computationally inferred ABCA3 haplotypes. The population-based frequencies of 121ins2, E292V, and I73T were rare (<0.4%). E292V was present in 3.8% of newborns with RDS, a 10-fold greater prevalence than in the Missouri cohort (p < 0.001). We did not identify other loss of function mutations in ABCA3 among patients with E292V that would account for their RDS. E292V occurred on a unique haplotype that was derived from a recombination of two common ABCA3 haplotypes. E292V was over-represented in newborns with RDS suggesting that E292V or its unique haplotype impart increased genetic risk for RDS.

Recombination as a mechanism for sporadic mutation in the surfactant protein-C gene

OBJECTIVE

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.

STUDY DESIGN

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 (+/-50 kb) 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 = 30).

RESULTS

We detected strong evidence (weak LD and BFs > 1,400) 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 four of five families; the remaining child had evidence for somatic recombination on the mutated allele.

CONCLUSIONS

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.

A major deletion in the surfactant protein-B gene causing lethal respiratory distress

BACKGROUND

Loss of function mutations in the surfactant protein-B gene (SFTPB) cause lethal neonatal respiratory distress due to reduced or absent expression of mature surfactant protein B (SP-B, encoded in exons 6 and 7). No large deletions in SFTPB have been previously identified.

AIM

Genomic, proteomic and immunohistochemical characterization of a 3 kb deletion in SFTPB.

METHODS

A full-term newborn presented with refractory respiratory failure. We amplified and sequenced SFTPB from the infant and both parents, determined SP-B protein expression in tracheal aspirate samples using Western-blot analysis, and performed immunohistochemical staining and electron microscopy of lung biopsy tissue.

RESULTS

The infant was homozygous for a 2958 bp deletion in SFTPB that included exons 7 and 8. Both asymptomatic parents were heterozygous for the deletion. A truncated mature SP-B peptide was detected on Western blotting of tracheal aspirate. Amino acid sequence specific to that encoded in exon 5 was present, but that encoded by exon 7 was absent. ProSP-B expression was robust within alveolar type II cells and lamellar body structure was disrupted.

CONCLUSIONS

This deletion in SFTPB resulted in SP-B deficiency due to absence of elements in mature SP-B that are critical for appropriate peptide folding, trafficking and processing.

Informed consent for genetic research

BACKGROUND

Rapid technological advances in genetic research and public concern about genetic discrimination have led to anticipatory safeguards in the informed consent process in the absence of legal examples of proven discrimination. Despite federal and state regulations to restrict access to personal health information, including genetic information, institutional review boards have required the addition of language to informed consent documents that warns about the risks of discrimination with participation in genetic research.

OBJECTIVE

To determine the reasons that families refused consent for their infant's participation in a study evaluating a genetic cause of respiratory distress syndrome.

DESIGN

Survey conducted between February 1, 2002, and March 31, 2003.

SETTING

Academic, tertiary free-standing children's hospital.

PARTICIPANTS

A convenience sample of 465 families were approached for consent. The 135 families who refused consent were surveyed.

MAIN OUTCOME MEASURES

Reasons for refusal.

RESULTS

Of the nonconsenting families, 79% spontaneously and specifically identified institutionally required language in our consent form concerning the risk of denial of access to health insurance and employment as the primary reason for refusal; 97% indicated that their fears resulted directly from language in our consent form. Only 20% of families who refused consent cited inadequate time to consider the study.

CONCLUSIONS

The institutionally required description of risk of genetic discrimination due solely to participation in genetic research was the primary reason for refusal to consent in this cohort. Information about federally and institutionally mandated protections for confidentiality of participants in genetic research should be included in the informed consent document to balance the description of hypothetical risks and more accurately inform subjects.