Homozygous missense variants in YKT6 result in loss of function and are associated with developmental delay, with or without severe infantile liver disease and risk for hepatocellular carcinoma.

PURPOSE: YKT6 plays important roles in multiple intracellular vesicle trafficking events but has not been associated with Mendelian diseases. METHODS: We report 3 unrelated individuals with rare homozygous missense variants in YKT6 who exhibited neurological disease with or without a progressive infantile liver disease. We modeled the variants in Drosophila. We generated wild-type and variant genomic rescue constructs of the fly ortholog dYkt6 and compared their ability in rescuing the loss-of-function phenotypes in mutant flies. We also generated a dYkt6KozakGAL4 allele to assess the expression pattern of dYkt6. RESULTS: Two individuals are homozygous for YKT6 [NM_006555.3:c.554A>G p.(Tyr185Cys)] and exhibited normal prenatal course followed by failure to thrive, developmental delay, and progressive liver disease. Haplotype analysis identified a shared homozygous region flanking the variant, suggesting a common ancestry. The third individual is homozygous for YKT6 [NM_006555.3:c.191A>G p.(Tyr64Cys)] and exhibited neurodevelopmental disorders and optic atrophy. Fly dYkt6 is essential and is expressed in the fat body (analogous to liver) and central nervous system. Wild-type genomic rescue constructs can rescue the lethality and autophagic flux defects, whereas the variants are less efficient in rescuing the phenotypes. CONCLUSION: The YKT6 variants are partial loss-of-function alleles, and the p.(Tyr185Cys) is more severe than p.(Tyr64Cys).

A disorder-related variant (E420K) of a PP2A-regulatory subunit (PPP2R5D) causes constitutively active AKT-mTOR signaling and uncoordinated cell growth.

Functional genomic approaches have facilitated the discovery of rare genetic disorders and improved efforts to decipher their underlying etiology. PPP2R5D-related disorder is an early childhood onset condition characterized by intellectual disability, hypotonia, autism-spectrum disorder, macrocephaly, and dysmorphic features. The disorder is caused by de novo single nucleotide changes in PPP2R5D, which generate heterozygous dominant missense variants. PPP2R5D is known to encode a B'-type (B'56δ) regulatory subunit of a PP2A-serine/threonine phosphatase. To help elucidate the molecular mechanisms altered in PPP2R5D-related disorder, we used a CRISPR-single-base editor to generate HEK-293 cells in which a single transition (c.1258G>A) was introduced into one allele, precisely recapitulating a clinically relevant E420K variant. Unbiased quantitative proteomic and phosphoproteomic analyses of endogenously expressed proteins revealed heterozygous-dominant changes in kinase/phosphatase signaling. These data combined with orthogonal validation studies revealed a previously unrecognized interaction of PPP2R5D with AKT in human cells, leading to constitutively active AKT-mTOR signaling, increased cell size, and uncoordinated cellular growth in E420K-variant cells. Rapamycin reduced cell size and dose-dependently reduced RPS6 phosphorylation in E420K-variant cells, suggesting that inhibition of mTOR1 can suppress both the observed RPS6 hyperphosphorylation and increased cell size. Together, our findings provide a deeper understanding of PPP2R5D and insight into how the E420K-variant alters signaling networks influenced by PPP2R5D. Our comprehensive approach, which combines precise genome editing, isobaric tandem mass tag labeling of peptides generated from endogenously expressed proteins, and concurrent liquid chromatography-mass spectrometry (LC-MS3), also provides a roadmap that can be used to rapidly explore the etiologies of additional genetic disorders.

Bi-allelic PAGR1 variants are associated with microcephaly and a severe neurodevelopmental disorder: Genetic evidence from two families.

Exome and genome sequencing were used to identify the genetic etiology of a severe neurodevelopmental disorder in two unrelated Ashkenazi Jewish families with three affected individuals. The clinical findings included a prenatal presentation of microcephaly, polyhydramnios and clenched hands while postnatal findings included microcephaly, severe developmental delay, dysmorphism, neurologic deficits, and death in infancy. A shared rare homozygous, missense variant (c.274A > G; p.Ser92Gly, NM_024516.4) was identified in PAGR1, a gene currently not associated with a Mendelian disease. PAGR1 encodes a component of the histone methyltransferase MLL2/MLL3 complex and may function in the DNA damage response pathway. Complete knockout of the murine Pagr1a is embryonic-lethal. Given the available evidence, PAGR1 is a strong candidate gene for a novel autosomal recessive severe syndromic neurodevelopmental disorder.

A homozygous splice variant in ATP5PO, disrupts mitochondrial complex V function and causes Leigh syndrome in two unrelated families.

Mitochondrial complex V plays an important role in oxidative phosphorylation by catalyzing the generation of ATP. Most complex V subunits are nuclear encoded and not yet associated with recognized Mendelian disorders. Using exome sequencing, we identified a rare homozygous splice variant (c.87+3A>G) in ATP5PO, the complex V subunit which encodes the oligomycin sensitivity conferring protein, in three individuals from two unrelated families, with clinical suspicion of a mitochondrial disorder. These individuals had a similar, severe infantile and often lethal multi-systemic disorder that included hypotonia, developmental delay, hypertrophic cardiomyopathy, progressive epileptic encephalopathy, progressive cerebral atrophy, and white matter abnormalities on brain MRI consistent with Leigh syndrome. cDNA studies showed a predominant shortened transcript with skipping of exon 2 and low levels of the normal full-length transcript. Fibroblasts from the affected individuals demonstrated decreased ATP5PO protein, defective assembly of complex V with markedly reduced amounts of peripheral stalk proteins, and complex V hydrolytic activity. Further, expression of human ATP5PO cDNA without exon 2 (hATP5PO-∆ex2) in yeast cells deleted for yATP5 (ATP5PO homolog) was unable to rescue growth on media which requires oxidative phosphorylation when compared to the wild type construct (hATP5PO-WT), indicating that exon 2 deletion leads to a non-functional protein. Collectively, our findings support the pathogenicity of the ATP5PO c.87+3A>G variant, which significantly reduces but does not eliminate complex V activity. These data along with the recent report of an affected individual with ATP5PO variants, add to the evidence that rare biallelic variants in ATP5PO result in defective complex V assembly, function and are associated with Leigh syndrome.

A novel homozygous variant in TRAPPC2L results in a neurodevelopmental disorder and disrupts TRAPP complex function.

BACKGROUND: Next-generation sequencing has facilitated the diagnosis of neurodevelopmental disorders with variable and non-specific clinical findings. Recently, a homozygous missense p.(Asp37Tyr) variant in TRAPPC2L, a core subunit of TRAPP complexes which function as tethering factors during membrane trafficking, was reported in two unrelated individuals with neurodevelopmental delay, post-infectious encephalopathy-associated developmental arrest, tetraplegia and accompanying rhabdomyolysis. METHODS: We performed whole genome sequencing on members of an Ashkenazi Jewish pedigree to identify the underlying genetic aetiology of global developmental delay/intellectual disability in three affected siblings. To assess the effect of the identified TRAPPC2L variant, we performed biochemical and cell biological functional studies on the TRAPPC2L protein. RESULTS: A rare homozygous predicted deleterious missense variant, p.(Ala2Gly), in TRAPPC2L was identified in the affected siblings and it segregated with the neurodevelopmental phenotype within the family. Using a yeast two-hybrid assay and in vitro binding, we demonstrate that the p.(Ala2Gly) variant, but not the p.(Asp37Tyr) variant, disrupted the interaction between TRAPPC2L and another core TRAPP protein, TRAPPC6a. Size exclusion chromatography suggested that this variant affects the assembly of TRAPP complexes. Employing two different membrane trafficking assays using fibroblasts from one of the affected siblings, we found a delay in traffic into and out of the Golgi. Similar to the p.(Asp37Tyr) variant, the p.(Ala2Gly) variant resulted in an increase in the levels of active RAB11. CONCLUSION: Our data fill in a gap in the knowledge of TRAPP architecture with TRAPPC2L interacting with TRAPPC6a, positioning it as a putative adaptor for other TRAPP subunits. Collectively, our findings support the pathogenicity of the TRAPPC2L p.(Ala2Gly) variant.

Peroxisomes can oxidize medium- and long-chain fatty acids through a pathway involving ABCD3 and HSD17B4.

Peroxisomes are essential organelles for the specialized oxidation of a wide variety of fatty acids, but they are also able to degrade fatty acids that are typically handled by mitochondria. Using a combination of pharmacological inhibition and clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR associated protein 9 genome editing technology to simultaneously manipulate peroxisomal and mitochondrial fatty acid ß-oxidation (FAO) in HEK-293 cells, we identified essential players in the metabolic crosstalk between these organelles. Depletion of carnitine palmitoyltransferase (CPT)2 activity through pharmacological inhibition or knockout (KO) uncovered a significant residual peroxisomal oxidation of lauric and palmitic acid, leading to the production of peroxisomal acylcarnitine intermediates. Generation and analysis of additional single- and double-KO cell lines revealed that the D-bifunctional protein (HSD17B4) and the peroxisomal ABC transporter ABCD3 are essential in peroxisomal oxidation of lauric and palmitic acid. Our results indicate that peroxisomes not only accept acyl-CoAs but can also oxidize acylcarnitines in a similar biochemical pathway. By using an Hsd17b4 KO mouse model, we demonstrated that peroxisomes contribute to the plasma acylcarnitine profile after acute inhibition of CPT2, proving in vivo relevance of this pathway. We summarize that peroxisomal FAO is important when mitochondrial FAO is defective or overloaded.-Violante, S., Achetib, N., van Roermund, C. W. T., Hagen, J., Dodatko, T., Vaz, F. M., Waterham, H. R., Chen, H., Baes, M., Yu, C., Argmann, C. A., Houten, S. M. Peroxisomes can oxidize medium- and long-chain fatty acids through a pathway involving ABCD3 and HSD17B4.

Mild inborn errors of metabolism in commonly used inbred mouse strains.

Inbred mouse strains are a cornerstone of translational research but paradoxically many strains carry mild inborn errors of metabolism. For example, α-aminoadipic acidemia and branched-chain ketoacid dehydrogenase deficiency are known in C57BL/6J mice. Using RNA sequencing, we now reveal the causal variants in Dhtkd1 and Bckdhb, and the molecular mechanism underlying these metabolic defects. C57BL/6J mice have decreased Dhtkd1 mRNA expression due to a solitary long terminal repeat (LTR) in intron 4 of Dhtkd1. This LTR harbors an alternate splice donor site leading to a partial splicing defect and as a consequence decreased total and functional Dhtkd1 mRNA, decreased DHTKD1 protein and α-aminoadipic acidemia. Similarly, C57BL/6J mice have decreased Bckdhb mRNA expression due to an LTR retrotransposon in intron 1 of Bckdhb. This transposable element encodes an alternative exon 1 causing aberrant splicing, decreased total and functional Bckdhb mRNA and decreased BCKDHB protein. Using a targeted metabolomics screen, we also reveal elevated plasma C5-carnitine in 129 substrains. This biochemical phenotype resembles isovaleric acidemia and is caused by an exonic splice mutation in Ivd leading to partial skipping of exon 10 and IVD protein deficiency. In summary, this study identifies three causal variants underlying mild inborn errors of metabolism in commonly used inbred mouse strains.

A next generation sequencing based approach to identify extracellular vesicle mediated mRNA transfers between cells.

BACKGROUND: Exosomes and other extracellular vesicles (EVs) have emerged as an important mechanism of cell-to-cell communication. However, previous studies either did not fully resolve what genetic materials were shuttled by exosomes or only focused on a specific set of miRNAs and mRNAs. A more systematic method is required to identify the genetic materials that are potentially transferred during cell-to-cell communication through EVs in an unbiased manner. RESULTS: In this work, we present a novel next generation of sequencing (NGS) based approach to identify EV mediated mRNA exchanges between co-cultured adipocyte and macrophage cells. We performed molecular and genomic profiling and jointly considered data from RNA sequencing (RNA-seq) and genotyping to track the "sequence varying mRNAs" transferred between cells. We identified 8 mRNAs being transferred from macrophages to adipocytes and 21 mRNAs being transferred in the opposite direction. These mRNAs represented biological functions including extracellular matrix, cell adhesion, glycoprotein, and signal peptides. CONCLUSIONS: Our study sheds new light on EV mediated RNA communications between adipocyte and macrophage cells, which may play a significant role in developing insulin resistance in diabetic patients. This work establishes a new method that is applicable to examining genetic material exchanges in many cellular systems and has the potential to be extended to in vivo studies as well.

Germline deletion of Krüppel-like factor 14 does not increase risk of diet induced metabolic syndrome in male C57BL/6 mice.

OBJECTIVE: The transcription factor Krüppel-like factor 14 (KLF14) has been associated with type 2 diabetes and high-density lipoprotein-cholesterol (HDL-C) through genome-wide association studies. The mechanistic underpinnings of KLF14's control of metabolic processes remain largely unknown. We studied the physiological roles of KLF14 in a knockout (KO) mouse model. METHODS: Male whole body Klf14 KO mice were fed a chow or high fat diet (HFD) and diet induced phenotypes were analyzed. Additionally, tissue-specific expression of Klf14 was determined using RT-PCR, RNA sequencing, immunoblotting and whole mount lacZ staining. Finally, the consequences of KLF14 loss-of-function were studied using RNA sequencing in tissues with relatively high Klf14 expression levels. RESULTS: KLF14 loss-of-function did not affect HFD-induced weight gain or insulin resistance. Fasting plasma concentrations of glucose, insulin, cholesterol, HDL-C and ApoA-I were also comparable between Klf14+/+ and Klf14-/- mice on chow and HFD. We found that in mice expression of Klf14 was the highest in the anterior pituitary (adenohypophysis), lower but detectable in white adipose tissue and undetectable in liver. Loss of KLF14 function impacted on the pituitary transcriptome with extracellular matrix organization as the primary affected pathway and a predicted link to glucocorticoid receptor signaling. CONCLUSIONS: Whole body loss of KLF14 function in male mice does not result in metabolic abnormalities as assessed under chow and HFD conditions. Mostly likely there is redundancy for the role of KLF14 in the mouse and a diverging function in humans.

Novel, compound heterozygous, single-nucleotide variants in MARS2 associated with developmental delay, poor growth, and sensorineural hearing loss.

Novel, single-nucleotide mutations were identified in the mitochondrial methionyl amino-acyl tRNA synthetase gene (MARS2) via whole exome sequencing in two affected siblings with developmental delay, poor growth, and sensorineural hearing loss.We show that compound heterozygous mutations c.550C>T:p.Gln 184* and c.424C>T:p.Arg142Trp in MARS2 lead to decreased MARS2 protein levels in patient lymphoblasts. Analysis of respiratory complex enzyme activities in patient fibroblasts revealed decreased complex I and IV activities. Immunoblotting of patient fibroblast and lymphoblast samples revealed reduced protein levels of NDUFB8 and COXII, representing complex I and IV, respectively. Additionally, overexpression of wild-type MARS2 in patient fibroblasts increased NDUFB8 and COXII protein levels. These findings suggest that recessive single-nucleotide mutations in MARS2 are causative for a new mitochondrial translation deficiency disorder with a primary phenotype including developmental delay and hypotonia. Identification of additional patients with single-nucleotide mutations in MARS2 is necessary to determine if pectus carinatum is also a consistent feature of this syndrome.

Genetic basis of alpha-aminoadipic and alpha-ketoadipic aciduria.

Alpha-aminoadipic and alpha-ketoadipic aciduria is an autosomal recessive inborn error of lysine, hydroxylysine, and tryptophan degradation. To date, DHTKD1 mutations have been reported in two alpha-aminoadipic and alpha-ketoadipic aciduria patients. We have now sequenced DHTKD1 in nine patients diagnosed with alpha-aminoadipic and alpha-ketoadipic aciduria as well as one patient with isolated alpha-aminoadipic aciduria, and identified causal mutations in eight. We report nine novel mutations, including three missense mutations, two nonsense mutations, two splice donor mutations, one duplication, and one deletion and insertion. Two missense mutations, one of which was reported before, were observed in the majority of cases. The clinical presentation of this group of patients was inhomogeneous. Our results confirm that alpha-aminoadipic and alpha-ketoadipic aciduria is caused by mutations in DHTKD1, and further establish that DHTKD1 encodes the E1 subunit of the alpha-ketoadipic acid dehydrogenase complex.

Tissue specific expression of myosin IC isoforms.

BACKGROUND: Myosin IC is a single headed member of the myosin superfamily that localizes to the cytoplasm and the nucleus and is implicated in a variety of processes in both compartments. We recently identified a novel isoform of myosin IC and showed that the MYOIC gene in mammalian cells encodes three isoforms (isoforms A, B, and C) that differ only in the addition of short isoform-specific N-terminal peptides. The expression pattern of the isoforms and the mechanisms of expression regulation remain unknown. RESULTS: To determine the expression patterns of myosin IC isoforms, we performed a comprehensive expression analysis of the two myosin IC isoforms (isoform A and B) that contain isoform-specific sequences. By immunoblotting with isoform-specific antibodies and by qRT-PCR with isoform-specific primer we demonstrate that myosin IC isoforms A and B have distinct expression patterns in mouse tissues. Specifically, we show that myosin IC isoform A is expressed in a tissue specific pattern, while myosin IC isoform B is ubiquitously expressed at comparable levels in mouse tissues. CONCLUSIONS: The differences in the expression profile of the myosin IC isoforms indicate a tissue-specific MYOIC gene regulation and further suggest that the myosin IC isoforms, despite their high sequence homology, might have tissue-specific and isoform-specific functions.

Functional dissection of human cardiac enhancers and noncoding de novo variants in congenital heart disease.

Rare coding mutations cause ∼45% of congenital heart disease (CHD). Noncoding mutations that perturb cis-regulatory elements (CREs) likely contribute to the remaining cases, but their identification has been problematic. Using a lentiviral massively parallel reporter assay (lentiMPRA) in human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs), we functionally evaluated 6,590 noncoding de novo variants (ncDNVs) prioritized from the whole-genome sequencing of 750 CHD trios. A total of 403 ncDNVs substantially affected cardiac CRE activity. A majority increased enhancer activity, often at regions with undetectable reference sequence activity. Of ten DNVs tested by introduction into their native genomic context, four altered the expression of neighboring genes and iPSC-CM transcriptional state. To prioritize future DNVs for functional testing, we used the MPRA data to develop a regression model, EpiCard. Analysis of an independent CHD cohort by EpiCard found enrichment of DNVs. Together, we developed a scalable system to measure the effect of ncDNVs on CRE activity and deployed it to systematically assess the contribution of ncDNVs to CHD.

A temporal classifier predicts histopathology state and parses acute-chronic phasing in inflammatory bowel disease patients.

Previous studies have conducted time course characterization of murine colitis models through transcriptional profiling of differential expression. We characterize the transcriptional landscape of acute and chronic models of dextran sodium sulfate (DSS) and adoptive transfer (AT) colitis to derive temporal gene expression and splicing signatures in blood and colonic tissue in order to capture dynamics of colitis remission and relapse. We identify sub networks of patient-derived causal networks that are enriched in these temporal signatures to distinguish acute and chronic disease components within the broader molecular landscape of IBD. The interaction between the DSS phenotype and chronological time-point naturally defines parsimonious temporal gene expression and splicing signatures associated with acute and chronic phases disease (as opposed to ordinary time-specific differential expression/splicing). We show these expression and splicing signatures are largely orthogonal, i.e. affect different genetic bodies, and that using machine learning, signatures are predictive of histopathological measures from both blood and intestinal data in murine colitis models as well as an independent cohort of IBD patients. Through access to longitudinal multi-scale profiling from disease tissue in IBD patient cohorts, we can apply this machine learning pipeline to generation of direct patient temporal multimodal regulatory signatures for prediction of histopathological outcomes.

Identification and validation of candidate risk genes in endocytic vesicular trafficking associated with esophageal atresia and tracheoesophageal fistulas.

Esophageal atresias/tracheoesophageal fistulas (EA/TEF) are rare congenital anomalies caused by aberrant development of the foregut. Previous studies indicate that rare or de novo genetic variants significantly contribute to EA/TEF risk, and most individuals with EA/TEF do not have pathogenic genetic variants in established risk genes. To identify the genetic contributions to EA/TEF, we performed whole genome sequencing of 185 trios (probands and parents) with EA/TEF, including 59 isolated and 126 complex cases with additional congenital anomalies and/or neurodevelopmental disorders. There was a significant burden of protein-altering de novo coding variants in complex cases (p = 3.3 × 10-4), especially in genes that are intolerant of loss-of-function variants in the population. We performed simulation analysis of pathway enrichment based on background mutation rate and identified a number of pathways related to endocytosis and intracellular trafficking that as a group have a significant burden of protein-altering de novo variants. We assessed 18 variants for disease causality using CRISPR-Cas9 mutagenesis in Xenopus and confirmed 13 with tracheoesophageal phenotypes. Our results implicate disruption of endosome-mediated epithelial remodeling as a potential mechanism of foregut developmental defects. Our results suggest significant genetic heterogeneity of EA/TEF and may have implications for the mechanisms of other rare congenital anomalies.

Erratum: Identification and validation of candidate risk genes in endocytic vesicular trafficking associated with esophageal atresia and tracheoesophageal fistulas.

[This corrects the article DOI: 10.1016/j.xhgg.2022.100107.].

Genome-Wide De Novo Variants in Congenital Heart Disease Are Not Associated With Maternal Diabetes or Obesity.

BACKGROUND: Congenital heart disease (CHD) is the most common anomaly at birth, with a prevalence of ≈1%. While infants born to mothers with diabetes or obesity have a 2- to 3-fold increased incidence of CHD, the cause of the increase is unknown. Damaging de novo variants (DNV) in coding regions are more common among patients with CHD, but genome-wide rates of coding and noncoding DNVs associated with these prenatal exposures have not been studied in patients with CHD. METHODS: DNV frequencies were determined for 1812 patients with CHD who had whole-genome sequencing and prenatal history data available from the Pediatric Cardiac Genomics Consortium's CHD GENES study (Genetic Network). The frequency of DNVs was compared between subgroups using t test or linear model. RESULTS: Among 1812 patients with CHD, the number of DNVs per patient was higher with maternal diabetes (76.5 versus 72.1, t test P=3.03×10-11), but the difference was no longer significant after including parental ages in a linear model (paternal and maternal correction P=0.42). No interaction was observed between diabetes risk and parental age (paternal and maternal interaction P=0.80 and 0.68, respectively). No difference was seen in DNV count per patient based on maternal obesity (72.0 versus 72.2 for maternal body mass index <25 versus maternal body mass index >30, t test P=0.86). CONCLUSIONS: After accounting for parental age, the offspring of diabetic or obese mothers have no increase in DNVs compared with other children with CHD. These results emphasize the role for other mechanisms in the cause of CHD associated with these prenatal exposures. REGISTRATION: URL: https://clinicaltrials.gov; NCT01196182.

Rare variant analysis of 4241 pulmonary arterial hypertension cases from an international consortium implicates FBLN2, PDGFD, and rare de novo variants in PAH.

BACKGROUND: Pulmonary arterial hypertension (PAH) is a lethal vasculopathy characterized by pathogenic remodeling of pulmonary arterioles leading to increased pulmonary pressures, right ventricular hypertrophy, and heart failure. PAH can be associated with other diseases (APAH: connective tissue diseases, congenital heart disease, and others) but often the etiology is idiopathic (IPAH). Mutations in bone morphogenetic protein receptor 2 (BMPR2) are the cause of most heritable cases but the vast majority of other cases are genetically undefined. METHODS: To identify new risk genes, we utilized an international consortium of 4241 PAH cases with exome or genome sequencing data from the National Biological Sample and Data Repository for PAH, Columbia University Irving Medical Center, and the UK NIHR BioResource - Rare Diseases Study. The strength of this combined cohort is a doubling of the number of IPAH cases compared to either national cohort alone. We identified protein-coding variants and performed rare variant association analyses in unrelated participants of European ancestry, including 1647 IPAH cases and 18,819 controls. We also analyzed de novo variants in 124 pediatric trios enriched for IPAH and APAH-CHD. RESULTS: Seven genes with rare deleterious variants were associated with IPAH with false discovery rate smaller than 0.1: three known genes (BMPR2, GDF2, and TBX4), two recently identified candidate genes (SOX17, KDR), and two new candidate genes (fibulin 2, FBLN2; platelet-derived growth factor D, PDGFD). The new genes were identified based solely on rare deleterious missense variants, a variant type that could not be adequately assessed in either cohort alone. The candidate genes exhibit expression patterns in lung and heart similar to that of known PAH risk genes, and most variants occur in conserved protein domains. For pediatric PAH, predicted deleterious de novo variants exhibited a significant burden compared to the background mutation rate (2.45×, p = 2.5e-5). At least eight novel pediatric candidate genes carrying de novo variants have plausible roles in lung/heart development. CONCLUSIONS: Rare variant analysis of a large international consortium identified two new candidate genes-FBLN2 and PDGFD. The new genes have known functions in vasculogenesis and remodeling. Trio analysis predicted that ~ 15% of pediatric IPAH may be explained by de novo variants.

Development and analytical validation of a next-generation sequencing based microsatellite instability (MSI) assay.

BACKGROUND: We have developed and analytically validated a next-generation sequencing (NGS) assay to classify microsatellite instability (MSI) in formalin-fixed paraffin-embedded (FFPE) tumor specimens. METHODOLOGY: The assay relies on DNA-seq evaluation of insertion/deletion (indel) variability at 29 highly informative genomic loci to estimate MSI status without the requirement for matched-normal tissue. The assay has a clinically relevant five-day turnaround time and can be conducted on as little as 20 ng genomic DNA with a batch size of up to forty samples in a single run. RESULTS: The MSI detection method was developed on a training set (n = 94) consisting of 22 MSI-H, 24 MSS, and 47 matched normal samples and tested on an independent test set of 24 MSI-H and 24 MSS specimens. Assay performance with respect to accuracy, reproducibility, precision as well as control sample performance was estimated across a wide range of FFPE samples of multiple histologies to address pre-analytical variability (percent tumor nuclei), and analytical variability (batch size, run, day, operator). Analytical precision studies demonstrated that the assay is highly reproducible and accurate as compared to established gold standard PCR methodology and has been validated through NYS CLEP. SIGNIFICANCE: This assay provides clinicians with robust and reproducible NGS-based MSI testing without the need of matched normal tissue to inform clinical decision making for patients with solid tumors.

Proliferative potential and resistance to immune checkpoint blockade in lung cancer patients.

BACKGROUND: Resistance to immune checkpoint inhibitors (ICIs) has been linked to local immunosuppression independent of major ICI targets (e.g., PD-1). Clinical experience with response prediction based on PD-L1 expression suggests that other factors influence sensitivity to ICIs in non-small cell lung cancer (NSCLC) patients. METHODS: Tumor specimens from 120 NSCLC patients from 10 institutions were evaluated for PD-L1 expression by immunohistochemistry, and global proliferative profile by targeted RNA-seq. RESULTS: Cell proliferation, derived from the mean expression of 10 proliferation-associated genes (namely BUB1, CCNB2, CDK1, CDKN3, FOXM1, KIAA0101, MAD2L1, MELK, MKI67, and TOP2A), was identified as a marker of response to ICIs in NSCLC. Poorly, moderately, and highly proliferative tumors were somewhat equally represented in NSCLC, with tumors with the highest PD-L1 expression being more frequently moderately proliferative as compared to lesser levels of PD-L1 expression. Proliferation status had an impact on survival in patients with both PD-L1 positive and negative tumors. There was a significant survival advantage for moderately proliferative tumors compared to their combined highly/poorly counterparts (p = 0.021). Moderately proliferative PD-L1 positive tumors had a median survival of 14.6 months that was almost twice that of PD-L1 negative highly/poorly proliferative at 7.6 months (p = 0.028). Median survival in moderately proliferative PD-L1 negative tumors at 12.6 months was comparable to that of highly/poorly proliferative PD-L1 positive tumors at 11.5 months, but in both instances less than that of moderately proliferative PD-L1 positive tumors. Similar to survival, proliferation status has impact on disease control (DC) in patients with both PD-L1 positive and negative tumors. Patients with moderately versus those with poorly or highly proliferative tumors have a superior DC rate when combined with any classification schema used to score PD-L1 as a positive result (i.e., TPS ≥ 50% or ≥ 1%), and best displayed by a DC rate for moderately proliferative tumors of no less than 40% for any classification of PD-L1 as a negative result. While there is an over representation of moderately proliferative tumors as PD-L1 expression increases this does not account for the improved survival or higher disease control rates seen in PD-L1 negative tumors. CONCLUSIONS: Cell proliferation is potentially a new biomarker of response to ICIs in NSCLC and is applicable to PD-L1 negative tumors.