VarSifter: visualizing and analyzing exome-scale sequence variation data on a desktop computer

Clinical and molecular heterogeneity of syndromic hypothalamic hamartoma

Two children are presented who have a distinct syndrome of multiple buccolingual frenula, a stiff and short fifth finger with small nails, a hypothalamic hamartoma, mild to moderate neurological impairment, and mild endocrinological symptoms. No variant assessed to be pathogenic or likely pathogenic was detected in the GLI3 gene in either child. This syndrome appears to be distinct from the inherited Pallister-Hall syndrome associated with GLI3 variants, which is characterized by hypothalamic hamartoma, mesoaxial polydactyly, and other anomalies. In the individuals described here, manifestations outside of the central nervous system were milder and the mesoaxial polydactyly, which is common in individuals with Pallister-Hall syndrome, was absent. Instead, these children had multiple buccolingual frenula together with the unusual appearance of the fifth digit. It remains unclear whether these two individuals represent a separate nosologic entity or if they represent a milder manifestation of one of the more severe syndromes associated with a hypothalamic hamartoma.

A pilot study to evaluate tissue- and plasma-based DNA driver mutations in a cohort of patients with pancreatic intraductal papillary mucinous neoplasms

Intraductal papillary mucinous neoplasms (IPMNs) are precursor lesions to pancreatic ductal adenocarcinoma that are challenging to manage due to limited imaging, cytologic, and molecular markers that accurately classify lesions, grade of dysplasia, or focus of invasion preoperatively. The objective of this pilot study was to determine the frequency and type of DNA mutations in a cohort of surgically resected, pathologically confirmed IPMN, and to determine if concordant mutations are detectable in paired pretreatment plasma samples. Formalin-fixed paraffin-embedded (FFPE) tissue from 46 surgically resected IPMNs (31 low-grade, 15 high-grade) and paired plasma from a subset of 15 IPMN cases (10 low-grade, 5 high-grade) were subjected to targeted mutation analysis using a QIAseq Targeted DNA Custom Panel. Common driver mutations were detected in FFPE from 44 of 46 (95.6%) IPMN cases spanning all grades; the most common DNA mutations included: KRAS (80%), RNF43 (24%), and GNAS (43%). Of note, we observed a significant increase in the frequency of RNF43 mutations from low-grade to high-grade IPMNs associated or concomitant with invasive carcinoma (trend test, P = 0.01). Among the subset of cases with paired plasma, driver mutations identified in the IPMNs were not detected in circulation. Overall, our results indicate that mutational burden for IPMNs is a common occurrence, even in low-grade IPMNs. Furthermore, although blood-based biopsies are an attractive, noninvasive method for detecting somatic DNA mutations, the QIAseq panel was not sensitive enough to detect driver mutations that existed in IPMN tissue using paired plasma in the volume we were able to retrieve for this retrospective study.

A neomorphic variant in SP7 alters sequence specificity and causes a high-turnover bone disorder

SP7/Osterix is a transcription factor critical for osteoblast maturation and bone formation. Homozygous loss-of-function mutations in SP7 cause osteogenesis imperfecta type XII, but neomorphic (gain-of-new-function) mutations of SP7 have not been reported in humans. Here we describe a de novo dominant neomorphic missense variant (c.926 C > G:p.S309W) in SP7 in a patient with craniosynostosis, cranial hyperostosis, and long bone fragility. Histomorphometry shows increased osteoblasts but decreased bone mineralization. Mice with the corresponding variant also show a complex skeletal phenotype distinct from that of Sp7-null mice. The mutation alters the binding specificity of SP7 from AT-rich motifs to a GC-consensus sequence (typical of other SP family members) and produces an aberrant gene expression profile, including increased expression of Col1a1 and endogenous Sp7, but decreased expression of genes involved in matrix mineralization. Our study identifies a pathogenic mechanism in which a mutation in a transcription factor shifts DNA binding specificity and provides important in vivo evidence that the affinity of SP7 for AT-rich motifs, unique among SP proteins, is critical for normal osteoblast differentiation.

ZBTB11 dysfunction: spectrum of brain abnormalities, biochemical signature and cellular consequences

Bi-allelic pathogenic variants in ZBTB11 have been associated with intellectual developmental disorder, autosomal recessive 69 (MRT69; OMIM 618383). We report five patients from three families with novel, bi-allelic variants in ZBTB11. We have expanded the clinical phenotype of MRT69, documenting varied severity of atrophy affecting different brain regions and described combined malonic and methylmalonic aciduria as a biochemical manifestation. As ZBTB11 encodes for a transcriptional regulator, we performeded chromatin immunoprecipitation-sequencing targeting ZBTB11 in fibroblasts from patients and controls. Chromatin immunoprecipitation-sequencing revealed binding of wild-type ZBTB11 to promoters in 238 genes, among which genes encoding proteins involved in mitochondrial functions and RNA processing are over-represented. Mutated ZBTB11 showed reduced binding to 61 of the targeted genes, indicating that the variants act as loss of function. Most of these genes are related to mitochondrial functions. Transcriptome analysis of the patient fibroblasts revealed dysregulation of mitochondrial functions. In addition, we uncovered that reduced binding of the mutated ZBTB11 to ACSF3 leads to decreased ACSF3 transcript level, explaining combined malonic and methylmalonic aciduria. Collectively, these results expand the clinical spectrum of ZBTB11-related neurological disease and give insight into the pathophysiology in which the dysfunctional ZBTB11 affect mitochondrial functions and RNA processing contributing to the neurological and biochemical phenotypes.

KLF3 and PAX6 are candidate driver genes in late-stage, MSI-hypermutated endometrioid endometrial carcinomas

Endometrioid endometrial carcinomas (EECs) are the most common histological subtype of uterine cancer. Late-stage disease is an adverse prognosticator for EEC. The purpose of this study was to analyze EEC exome mutation data to identify late-stage-specific statistically significantly mutated genes (SMGs), which represent candidate driver genes potentially associated with disease progression. We exome sequenced 15 late-stage (stage III or IV) non-ultramutated EECs and paired non-tumor DNAs; somatic variants were called using Strelka, Shimmer, SomaticSniper and MuTect. Additionally, somatic mutation calls were extracted from The Cancer Genome Atlas (TCGA) data for 66 late-stage and 270 early-stage (stage I or II) non-ultramutated EECs. MutSigCV (v1.4) was used to annotate SMGs in the two late-stage cohorts and to derive p-values for all mutated genes in the early-stage cohort. To test whether late-stage SMGs are statistically significantly mutated in early-stage tumors, q-values for late-stage SMGs were re-calculated from the MutSigCV (v1.4) early-stage p-values, adjusting for the number of late-stage SMGs tested. We identified 14 SMGs in the combined late-stage EEC cohorts. When the 14 late-stage SMGs were examined in the TCGA early-stage data, only Krüppel-like factor 3 (KLF3) and Paired box 6 (PAX6) failed to reach significance as early-stage SMGs, despite the inclusion of enough early-stage cases to ensure adequate statistical power. Within TCGA, nonsynonymous mutations in KLF3 and PAX6 were, respectively, exclusive or nearly exclusive to the microsatellite instability (MSI)-hypermutated molecular subgroup and were dominated by insertions-deletions at homopolymer tracts. In conclusion, our findings are hypothesis-generating and suggest that KLF3 and PAX6, which encode transcription factors, are MSI target genes and late-stage-specific SMGs in EEC.

BET1 variants establish impaired vesicular transport as a cause for muscular dystrophy with epilepsy

BET1 is required, together with its SNARE complex partners GOSR2, SEC22b, and Syntaxin-5 for fusion of endoplasmic reticulum-derived vesicles with the ER-Golgi intermediate compartment (ERGIC) and the cis-Golgi. Here, we report three individuals, from two families, with severe congenital muscular dystrophy (CMD) and biallelic variants in BET1 (P1 p.(Asp68His)/p.(Ala45Valfs*2); P2 and P3 homozygous p.(Ile51Ser)). Due to aberrant splicing and frameshifting, the variants in P1 result in low BET1 protein levels and impaired ER-to-Golgi transport. Since in silico modeling suggested that p.(Ile51Ser) interferes with binding to interaction partners other than SNARE complex subunits, we set off and identified novel BET1 interaction partners with low affinity for p.(Ile51Ser) BET1 protein compared to wild-type, among them ERGIC-53. The BET1/ERGIC-53 interaction was validated by endogenous co-immunoprecipitation with both proteins colocalizing to the ERGIC compartment. Mislocalization of ERGIC-53 was observed in P1 and P2's derived fibroblasts; while in the p.(Ile51Ser) P2 fibroblasts specifically, mutant BET1 was also mislocalized along with ERGIC-53. Thus, we establish BET1 as a novel CMD/epilepsy gene and confirm the emerging role of ER/Golgi SNAREs in CMD.

Biallelic mutations in ELFN1 gene associated with developmental and epileptic encephalopathy and joint laxity

ELFN1, a transmembrane leucine rich repeat protein, is involved in signal transduction in both neural cells and ROD ON-bipolar synaptogenesis. We present three siblings with developmental and epileptic encephalopathy and co-morbidities due to ELFN1 gene mutation; this is the first report in literature defining the human phenotype of ELFN1 gene mutation. Clinical, electrophysiological, and radiological findings along with comprehensive genetic studies of the patients and their family members are presented. Developmental and epileptic encephalopathy, autistic features, pyramidal signs, joint laxity, and dysmorphic features are the characteristic findings of this new clinical entity, involving mainly nervous system and possibly connective tissue. Whole exome sequence analysis followed by Sanger sequencing in all family members revealed disease-causing 8 bp frameshift mutation depicted as NM_001128636.2: c.42_49delGGCCGCCA; p. (Ala15Profs*241) in ELFN1. The variant, located in the signal peptide domain in the ELFN1 gene, was found to be homozygous in three patients, and heterozygous in the parents and three healthy siblings. Segregation analysis in family members together with pathogenicity assessment tools strongly supported the damaging effect of the frameshift variant on the function of the ELFN1 protein. Mutations in ELFN1 gene may be considered in patients with neonatal and infantile-onset epileptic encephalopathy before the full clinical picture is apparent.

Evidence That the Etiology of Congenital Hypopituitarism Has a Major Genetic Component but Is Infrequently Monogenic

Purpose

Congenital hypopituitarism usually occurs sporadically. In most patients, the etiology remains unknown.

Methods

We studied 13 children with sporadic congenital hypopituitarism. Children with non-endocrine, non-familial idiopathic short stature (NFSS) (n = 19) served as a control group. Exome sequencing was performed in probands and both unaffected parents. A burden testing approach was used to compare the number of candidate variants in the two groups.

Results

First, we assessed the frequency of rare, predicted-pathogenic variants in 42 genes previously reported to be associated with pituitary gland development. The average number of variants per individual was greater in probands with congenital hypopituitarism than those with NFSS (1.1 vs. 0.21, mean variants/proband, P = 0.03). The number of probands with at least 1 variant in a pituitary-associated gene was greater in congenital hypopituitarism than in NFSS (62% vs. 21%, P = 0.03). Second, we assessed the frequency of rare, predicted-pathogenic variants in the exome (to capture undiscovered causes) that were inherited in a fashion that could explain the sporadic occurrence of the proband's condition with a monogenic etiology (de novo mutation, autosomal recessive, or X-linked recessive) with complete penetrance. There were fewer monogenic candidates in the probands with congenital hypopituitarism than those with NFSS (1.3 vs. 2.5 candidate variants/proband, P = 0.024). We did not find any candidate variants (0 of 13 probands) in genes previously reported to explain the phenotype in congenital hypopituitarism, unlike NFSS (8 of 19 probands, P = 0.01).

Conclusion

Our findings provide evidence that the etiology of sporadic congenital hypopituitarism has a major genetic component but may be infrequently monogenic with full penetrance, suggesting a more complex etiology.

Pathogenic variants in TNNC2 cause congenital myopathy due to an impaired force response to calcium

Troponin C (TnC) is a critical regulator of skeletal muscle contraction; it binds Ca2+ to activate muscle contraction. Surprisingly, the gene encoding fast skeletal TnC (TNNC2) has not yet been implicated in muscle disease. Here, we report 2 families with pathogenic variants in TNNC2. Patients present with a distinct, dominantly inherited congenital muscle disease. Molecular dynamics simulations suggested that the pathomechanisms by which the variants cause muscle disease include disruption of the binding sites for Ca2+ and for troponin I. In line with these findings, physiological studies in myofibers isolated from patients' biopsies revealed a markedly reduced force response of the sarcomeres to [Ca2+]. This pathomechanism was further confirmed in experiments in which contractile dysfunction was evoked by replacing TnC in myofibers from healthy control subjects with recombinant, mutant TnC. Conversely, the contractile dysfunction of myofibers from patients was repaired by replacing endogenous, mutant TnC with recombinant, wild-type TnC. Finally, we tested the therapeutic potential of the fast skeletal muscle troponin activator tirasemtiv in patients' myofibers and showed that the contractile dysfunction was repaired. Thus, our data reveal that pathogenic variants in TNNC2 cause congenital muscle disease, and they provide therapeutic angles to repair muscle contractility.

Clinical and genomic analysis of a large Chinese family with familial cortical myoclonic tremor with epilepsy and SAMD12 intronic repeat expansion

Objective

Our goal was to perform detailed clinical and genomic analysis of a large multigenerational Chinese family with 21 individuals showing symptoms of Familial Cortical Myoclonic Tremor with Epilepsy (FCMTE) that we have followed for over 20 years.

Methods

Patients were subjected to clinical evaluation, routine EEG, and structural magnetic resonance imaging. Whole exome sequencing, repeat-primed PCR, long-range PCR, and PacBio sequencing were performed to characterize the disease-causing mutation in this family.

Results

All evaluated patients manifested adult-onset seizures and presented with progressive myoclonic postural tremors starting after the third or fourth decade of life. Seizures typically diminished markedly in frequency with implementation of antiseizure medications but did not completely cease. The electroencephalogram of affected individuals showed generalized or multifocal spikes and slow wave complexes. An expansion of TTTTA motifs with addition of TTTCA motifs in intron 4 of SAMD12 was identified to segregate with the disease phenotype in this family. Furthermore, we found that the mutant allele is unstable and can undergo both contraction and expansion by changes in the number of repeat motifs each time it is passed to the next generation. The size of mutant allele varied from 5 to 5.5 kb with 549-603 copies of TTTTA and 287-343 copies of TTTCA repeat motifs in this family.

Significance

Our study provides a detailed description of clinical progression of FCMTE symptoms and its management with antiseizure medications. Our method of repeat analysis by PacBio sequencing of long-range PCR products does not require high-quality DNA and hence can be easily applied to other families to elucidate any correlation between the repeat size and phenotypic variables, such as, age of onset, and severity of symptoms.

A negative storage model for precise but compact storage of genetic variation data

Falling sequencing costs and large initiatives are resulting in increasing amounts of data available for investigator use. However, there are informatics challenges in being able to access genomic data. Performance and storage are well-appreciated issues, but precision is critical for meaningful analysis and interpretation of genomic data. There is an inherent accuracy vs. performance trade-off with existing solutions. The most common approach (Variant-only Storage Model, VOSM) stores only variant data. Systems must therefore assume that everything not variant is reference, sacrificing precision and potentially accuracy. A more complete model (Full Storage Model, FSM) would store the state of every base (variant, reference and missing) in the genome thereby sacrificing performance. A compressed variation of the FSM can store the state of contiguous regions of the genome as blocks (Block Storage Model, BLSM), much like the file-based gVCF model. We propose a novel approach by which this state is encoded such that both performance and accuracy are maintained. The Negative Storage Model (NSM) can store and retrieve precise genomic state from different sequencing sources, including clinical and whole exome sequencing panels. Reduced storage requirements are achieved by storing only the variant and missing states and inferring the reference state. We evaluate the performance characteristics of FSM, BLSM and NSM and demonstrate dramatic improvements in storage and performance using the NSM approach.

The Association of MUC16 Mutation with Tumor Mutation Burden and Its Prognostic Implications in Cutaneous Melanoma

BACKGROUND

MUC16 is a mucin marker that is frequently mutated in melanoma, but whether MUC16 mutations could be useful as a surrogate biomarker for tumor mutation burden (TMB) remains unclear.

METHODS

This study rigorously evaluates the MUC16 mutation as a clinical biomarker in cutaneous melanoma by utilizing genomic and clinical data from patient samples from The Cancer Genome Atlas (TCGA) and two independent validation cohorts. We further extended the analysis to studies with patients treated with immunotherapies.

RESULTS

Analysis results showed that samples with MUC16 mutations had a higher TMB than the samples of wild-type, with strong statistical significance (P < 0.001) in all melanoma cohorts tested. Associations between MUC16 mutations and TMB remained statistically significant after adjusting for potential confounding factors in the TCGA cohort [OR, 9.28 (95% confidence interval (CI), 5.18-17.39); P < 0.001], Moffitt cohort [OR, 31.95 (95% CI, 8.71-163.90); P < 0.001], and Yale cohort [OR, 8.09 (95% CI, 3.12-23.79); P < 0.01]. MUC16 mutations were also found to be associated with overall survival in the TCGA [HR, 0.62; (95% CI, 0.45-0.85); P < 0.01] and Moffitt cohorts [HR, 0.49 (95% CI, 0.28-0.87); P = 0.014]. Strikingly, MUC16 is the only top frequently mutated gene for which prognostic significance was observed. MUC16 mutations were also found valuable in predicting anti-CTLA-4 and anti-PD-1 therapy responses.

CONCLUSIONS

MUC16 mutation appears to be a useful predictive marker of global TMB and patient survival in melanoma.

IMPACT

This is, to the best of our knowledge, the first systematic evaluation of MUC16 mutation as a clinical biomarker and a predictive biomarker for immunotherapy in melanoma.

Comparative clinical and genomic analysis of neurofibromatosis type 2-associated cranial and spinal meningiomas

Neurofibromatosis type 2 (NF2) is an autosomal dominant Mendelian tumor predisposition disorder caused by germline pathogenic variants in the tumor suppressor NF2. Meningiomas are the second most common neoplasm in NF2, often occurring in multiple intracranial and spinal locations within the same patient. In this prospective longitudinal study, we assessed volumes and growth rates of ten spinal and ten cranial benign meningiomas in seven NF2 patients that concluded with surgical resection and performed whole-exome sequencing and copy-number variant (CNV) analysis of the tumors. Our comparison of the volume and the growth rate of NF2-associated spinal and cranial meningiomas point to the differences in timing of tumor initiation and/or to the differences in tumor progression (e.g., non-linear, saltatory growth) at these two anatomical locations. Genomic investigation of these tumors revealed that somatic inactivation of NF2 is the principal and perhaps the only driver of tumor initiation; and that tumor progression likely occurs via accumulation of CNVs, rather than point mutations. Results of this study contribute to a better understanding of NF2-associated meningiomas clinical behavior and their genetic underpinnings.

Mutations in COMP cause familial carpal tunnel syndrome

Carpal tunnel syndrome (CTS) is the most common peripheral nerve entrapment syndrome, affecting a large proportion of the general population. Genetic susceptibility has been implicated in CTS, but the causative genes remain elusive. Here, we report the identification of two mutations in cartilage oligomeric matrix protein (COMP) that segregate with CTS in two large families with or without multiple epiphyseal dysplasia (MED). Both mutations impair the secretion of COMP by tenocytes, but the mutation associated with MED also perturbs its secretion in chondrocytes. Further functional characterization of the CTS-specific mutation reveals similar histological and molecular changes of tendons/ligaments in patients’ biopsies and the mouse models. The mutant COMP fails to oligomerize properly and is trapped in the ER, resulting in ER stress-induced unfolded protein response and cell death, leading to inflammation, progressive fibrosis and cell composition change in tendons/ligaments. The extracellular matrix (ECM) organization is also altered. Our studies uncover a previously unrecognized mechanism in CTS pathogenesis.

Familial carpal tunnel syndrome (CTS) is common, but causal genes are not characterized. Here the authors report two CTS-related mutations in two large families that impair secretion of COMP in tenocytes, leading to ER stress-induced unfolded protein response, inflammation and fibrosis in patients and mouse models.

Revisiting the complex architecture of ALS in Turkey: Expanding genotypes, shared phenotypes, molecular networks, and a public variant database

The last decade has proven that amyotrophic lateral sclerosis (ALS) is clinically and genetically heterogeneous, and that the genetic component in sporadic cases might be stronger than expected. This study investigates 1,200 patients to revisit ALS in the ethnically heterogeneous yet inbred Turkish population. Familial ALS (fALS) accounts for 20% of our cases. The rates of consanguinity are 30% in fALS and 23% in sporadic ALS (sALS). Major ALS genes explained the disease cause in only 35% of fALS, as compared with ~70% in Europe and North America. Whole exome sequencing resulted in a discovery rate of 42% (53/127). Whole genome analyses in 623 sALS cases and 142 population controls, sequenced within Project MinE, revealed well-established fALS gene variants, solidifying the concept of incomplete penetrance in ALS. Genome-wide association studies (GWAS) with whole genome sequencing data did not indicate a new risk locus. Coupling GWAS with a coexpression network of disease-associated candidates, points to a significant enrichment for cell cycle- and division-related genes. Within this network, literature text-mining highlights DECR1, ATL1, HDAC2, GEMIN4, and HNRNPA3 as important genes. Finally, information on ALS-related gene variants in the Turkish cohort sequenced within Project MinE was compiled in the GeNDAL variant browser (www.gendal.org).

Consecutive medical exome analysis at a tertiary center: Diagnostic and health-economic outcomes

The utility of whole exome analysis has been extensively demonstrated in research settings, but its clinical utility as a first-tier genetic test has not been well documented from diagnostic and health economic standpoints in real-life clinical settings. We performed medical exome analyses focusing on a clinically interpretable portion of the genome (4,813 genes) as a first-tier genetic test for 360 consecutive patients visiting a genetics clinic at a tertiary children's hospital in Japan, over a 3-year period. Bioinformatics analyses were conducted using standard software. A molecular diagnosis was made in 171 patients involving a total of 107 causative genes. Among these 107 causative genes, 57 genes were classified as genes with potential organ-specific interventions and management strategies. Clinically relevant results were obtained in 26% of the total cohort and 54% of the patients with a definitive molecular diagnosis. Performing the medical exome analysis at the time of the initial visit to the tertiary center, rather than after visits to pertinent specialists, brain MRI examination, and G-banded chromosome testing, would have reduced the financial cost by 197 euros according to retrospective calculation under multiple assumption. The present study demonstrated a high diagnostic yield (47.5%) for singleton medical exome analysis as a first-tier test in a real-life setting. Medical exome analysis yielded clinically relevant information in a quarter of the total patient cohort. The application of genomic testing during the initial visit to a tertiary medical center could be a rational approach to the diagnosis of patients with suspected genetic disorders.

webGQT: A Shiny Server for Genotype Query Tools for Model-Based Variant Filtering

Summary

Genotype Query Tools (GQT) were developed to discover disease-causing variations from billions of genotypes and millions of genomes, processes data at substantially higher speed over other existing methods. While GQT has been available to a wide audience as command-line software, the difficulty of constructing queries among non-IT or non-bioinformatics researchers has limited its applicability. To overcome this limitation, we developed webGQT, an easy-to-use tool with a graphical user interface. With pre-built queries across three modules, webGQT allows for pedigree analysis, case-control studies, and population frequency studies. As a package, webGQT allows researchers with less or no applied bioinformatics/IT experience to mine potential disease-causing variants from billions.

Results

webGQT offers a flexible and easy-to-use interface for model-based candidate variant filtering for Mendelian diseases from thousands to millions of genomes at a reduced computation time. Additionally, webGQT provides adjustable parameters to reduce false positives and rescue missing genotypes across all modules. Using a case study, we demonstrate the applicability of webGQT to query non-human genomes. In addition, we demonstrate the scalability of webGQT on large data sets by implementing complex population-specific queries on the 1000 Genomes Project Phase 3 data set, which includes 8.4 billion variants from 2504 individuals across 26 different populations. Furthermore, webGQT supports filtering single-nucleotide variants, short insertions/deletions, copy number or any other variant genotypes supported by the VCF specification. Our results show that webGQT can be used as an online web service, or deployed on personal computers or local servers within research groups.

Availability

webGQT is made available to the users in three forms: 1) as a webserver available at https://vm1138.kaj.pouta.csc.fi/webgqt/, 2) as an R package to install on personal computers, and 3) as part of the same R package to configure on the user's own servers. The application is available for installation at https://github.com/arumds/webgqt.

Genomic diversity affects the accuracy of bacterial single-nucleotide polymorphism-calling pipelines

BACKGROUND

Accurately identifying single-nucleotide polymorphisms (SNPs) from bacterial sequencing data is an essential requirement for using genomics to track transmission and predict important phenotypes such as antimicrobial resistance. However, most previous performance evaluations of SNP calling have been restricted to eukaryotic (human) data. Additionally, bacterial SNP calling requires choosing an appropriate reference genome to align reads to, which, together with the bioinformatic pipeline, affects the accuracy and completeness of a set of SNP calls obtained. This study evaluates the performance of 209 SNP-calling pipelines using a combination of simulated data from 254 strains of 10 clinically common bacteria and real data from environmentally sourced and genomically diverse isolates within the genera Citrobacter, Enterobacter, Escherichia, and Klebsiella.

RESULTS

We evaluated the performance of 209 SNP-calling pipelines, aligning reads to genomes of the same or a divergent strain. Irrespective of pipeline, a principal determinant of reliable SNP calling was reference genome selection. Across multiple taxa, there was a strong inverse relationship between pipeline sensitivity and precision, and the Mash distance (a proxy for average nucleotide divergence) between reads and reference genome. The effect was especially pronounced for diverse, recombinogenic bacteria such as Escherichia coli but less dominant for clonal species such as Mycobacterium tuberculosis.

CONCLUSIONS

The accuracy of SNP calling for a given species is compromised by increasing intra-species diversity. When reads were aligned to the same genome from which they were sequenced, among the highest-performing pipelines was Novoalign/GATK. By contrast, when reads were aligned to particularly divergent genomes, the highest-performing pipelines often used the aligners NextGenMap or SMALT, and/or the variant callers LoFreq, mpileup, or Strelka.

High-throughput custom capture sequencing identifies novel mutations in coloboma-associated genes: Mutation in DNA-binding domain of retinoic acid receptor beta affects nuclear localization causing ocular coloboma

Uveal coloboma is a potentially blinding congenital ocular malformation caused by the failure of optic fissure closure during the fifth week of human gestation. We performed custom capture high‐throughput screening of 38 known coloboma‐associated genes in 66 families. Suspected causative novel variants were identified in TFAP2A and CHD7, as well as two previously reported variants of uncertain significance in RARB and BMP7. The variant in RARB, unlike previously reported disease mutations in the ligand‐binding domain, was a missense change in the highly conserved DNA‐binding domain predicted to affect the protein's DNA‐binding ability. In vitro studies revealed lower steady‐state protein levels, reduced transcriptional activity, and incomplete nuclear localization of the mutant RARB protein compared with wild‐type. Zebrafish studies showed that human RARB messenger RNA partially reduced the ocular phenotype caused by morpholino knockdown of rarga gene, a zebrafish homolog of human RARB. Our study indicates that sequence alterations in known coloboma genes account for a small percentage of coloboma cases and that mutations in the RARB DNA‐binding domain could result in human disease.

Germline variants in cancer genes in high-risk non-BRCA patients from Puerto Rico

Inherited pathogenic variants in genes that confer moderate to high risk of breast cancer may explain up to 50% of familial breast cancer. This study aimed at identifying inherited pathogenic variants in breast cancer cases from Puerto Rico that were not linked to BRCA1 or BRCA2. Forty-eight breast cancer patients that met the clinical criteria for BRCA testing but had received a negative BRCA1/2 result were recruited. Fifty-three genes previously implicated in hereditary cancer predisposition were captured using the BROCA Agilent cancer risk panel followed by massively parallel sequencing. Missense variants of uncertain clinical significance in CHEK2 were evaluated using an in vitro kinase assays to determine their impact on function. Pathogenic variants were identified in CHEK2, MUTYH, and RAD51B in four breast cancer patients, which represented 8.3% of the cohort. We identified three rare missense variants of uncertain significance in CHEK2 and two variants (p.Pro484Leu and p.Glu239Lys) showed markedly decreased kinase activity in vitro comparable to a known pathogenic variant. Interestingly, the local ancestry at the RAD51B locus in the carrier of p.Arg47* was predicted to be of African origin. In this cohort, 12.5% of the BRCA-negative breast cancer patients were found to carry a known pathogenic variant or a variant affecting protein activity. This study reveals an unmet clinical need of genetic testing that could benefit a significant proportion of at-risk Latinas. It also highlights the complexity of Hispanic populations as pathogenic factors may originate from any of the ancestral populations that make up their genetic backgrounds.

Proteogenomic landscape of squamous cell lung cancer

How genomic and transcriptomic alterations affect the functional proteome in lung cancer is not fully understood. Here, we integrate DNA copy number, somatic mutations, RNA-sequencing, and expression proteomics in a cohort of 108 squamous cell lung cancer (SCC) patients. We identify three proteomic subtypes, two of which (Inflamed, Redox) comprise 87% of tumors. The Inflamed subtype is enriched with neutrophils, B-cells, and monocytes and expresses more PD-1. Redox tumours are enriched for oxidation-reduction and glutathione pathways and harbor more NFE2L2/KEAP1 alterations and copy gain in the 3q2 locus. Proteomic subtypes are not associated with patient survival. However, B-cell-rich tertiary lymph node structures, more common in Inflamed, are associated with better survival. We identify metabolic vulnerabilities (TP63, PSAT1, and TFRC) in Redox. Our work provides a powerful resource for lung SCC biology and suggests therapeutic opportunities based on redox metabolism and immune cell infiltrates.

Squamous cell lung cancer has dismal prognosis due to the dearth of effective treatments. Here, the authors perform an integrated proteogenomic analysis of the disease, revealing three proteomics-based subtypes and suggesting potential therapeutic opportunities.

eDiVA-Classification and prioritization of pathogenic variants for clinical diagnostics

Mendelian diseases have shown to be an and efficient model for connecting genotypes to phenotypes and for elucidating the function of genes. Whole‐exome sequencing (WES) accelerated the study of rare Mendelian diseases in families, allowing for directly pinpointing rare causal mutations in genic regions without the need for linkage analysis. However, the low diagnostic rates of 20–30% reported for multiple WES disease studies point to the need for improved variant pathogenicity classification and causal variant prioritization methods. Here, we present the exome Disease Variant Analysis (eDiVA; http://ediva.crg.eu), an automated computational framework for identification of causal genetic variants (coding/splicing single‐nucleotide variants and small insertions and deletions) for rare diseases using WES of families or parent–child trios. eDiVA combines next‐generation sequencing data analysis, comprehensive functional annotation, and causal variant prioritization optimized for familial genetic disease studies. eDiVA features a machine learning‐based variant pathogenicity predictor combining various genomic and evolutionary signatures. Clinical information, such as disease phenotype or mode of inheritance, is incorporated to improve the precision of the prioritization algorithm. Benchmarking against state‐of‐the‐art competitors demonstrates that eDiVA consistently performed as a good or better than existing approach in terms of detection rate and precision. Moreover, we applied eDiVA to several familial disease cases to demonstrate its clinical applicability.

Female genitourinary tract melanoma: mutation analysis with clinicopathologic correlation: a single-institution experience

Female genitourinary tract melanoma (FGTM) is a rare and often-fatal form of mucosal melanoma. We describe our institutional experience with 55 cases of FGTM, 16 of which were evaluated with next-generation sequencing targeting 151 cancer-associated genes. Tumors tended to be thicker than conventional melanoma at presentation (median: 3.2 mm), were frequently ulcerated (50%), and characterized by incomplete initial resections. Regional lymph nodes showed tumor involvement at presentation in 28% of cases. With a median follow-up of 23.6 months, the median recurrence free survival was 14.5 months and the median overall survival was 29.6 months. Genomic analysis revealed mutually exclusive mutations in TP53 and KIT in 25%, while 19% of cases showed BRAF mutation. NRAS mutation was found in 13% of cases. Mutation in ATRX, previously undescribed in mucosal melanoma, was seen in three (10%) of 16 patients. Only invasive melanoma cases were included in statistical analyses. Patients with three or more mutations had marginally worse overall survival rates than those with two or less (P=0.07). Further studies are required for potential adjuvant treatment modalities to improve survival outcomes of FGTM.

Biallelic HEPHL1 variants impair ferroxidase activity and cause an abnormal hair phenotype

Maintenance of the correct redox status of iron is functionally important for critical biological processes. Multicopper ferroxidases play an important role in oxidizing ferrous iron, released from the cells, into ferric iron, which is subsequently distributed by transferrin. Two well-characterized ferroxidases, ceruloplasmin (CP) and hephaestin (HEPH) facilitate this reaction in different tissues. Recently, a novel ferroxidase, Hephaestin like 1 (HEPHL1), also known as zyklopen, was identified. Here we report a child with compound heterozygous mutations in HEPHL1 (NM_001098672) who presented with abnormal hair (pili torti and trichorrhexis nodosa) and cognitive dysfunction. The maternal missense mutation affected mRNA splicing, leading to skipping of exon 5 and causing an in-frame deletion of 85 amino acids (c.809_1063del; p.Leu271_ala355del). The paternal mutation (c.3176T>C; p.Met1059Thr) changed a highly conserved methionine that is part of a typical type I copper binding site in HEPHL1. We demonstrated that HEPHL1 has ferroxidase activity and that the patient's two mutations exhibited loss of this ferroxidase activity. Consistent with these findings, the patient's fibroblasts accumulated intracellular iron and exhibited reduced activity of the copper-dependent enzyme, lysyl oxidase. These results suggest that the patient's biallelic variants are loss-of-function mutations. Hence, we generated a Hephl1 knockout mouse model that was viable and had curly whiskers, consistent with the hair phenotype in our patient. These results enhance our understanding of the function of HEPHL1 and implicate altered ferroxidase activity in hair growth and hair disorders.

Genetic Variant Screening of DNA Repair Genes in Myelodysplastic Syndrome Identifies a Novel Mutation in the XRCC2 Gene

BACKGROUND

We aimed to detect single nucleotide polymorphisms (SNPs) and mutations in DNA repair genes and their possible association with myelodysplastic syndrome (MDS).

METHODS

Targeted enrichment resequencing of 84 DNA repair genes was initially performed on a screening cohort of MDS patients. Real-time polymerase chain reaction was used for genotyping selected SNPs in the validation cohort of patients.

RESULTS

A heterozygous frameshift mutation in the XRCC2 gene was identified. It leads to the formation of a truncated non-functional protein and decreased XRCC2 expression level. Decreased expression levels of all DNA repair genes functionally connected with mutated XRCC2 were also present. Moreover, a synonymous substitution in the PRKDC gene and 2 missense mutations in the SMUG1 and XRCC1 genes were also found. In the screening cohort, 6 candidate SNPs were associated with the tendency to develop MDS: rs4135113 (TDG, p = 0.03), rs12917 (MGMT, p = 0.003), rs2230641 (CCNH, p = 0.01), rs2228529 and rs2228526 (ERCC6, p = 0.04 and p = 0.03), and rs1799977 (MLH1, p = 0.04). In the validation cohort, only a polymorphism in MLH1 was significantly associated with development of MDS in patients with poor cytogenetics (p = 0.0004).

CONCLUSION

Our study demonstrates that genetic variants are present in DNA repair genes of MDS patients and may be associated with susceptibility to MDS.

Cysteinyl-tRNA Synthetase Mutations Cause a Multi-System, Recessive Disease That Includes Microcephaly, Developmental Delay, and Brittle Hair and Nails

Aminoacyl-tRNA synthetases (ARSs) are essential enzymes responsible for charging tRNA molecules with cognate amino acids. Consistent with the essential function and ubiquitous expression of ARSs, mutations in 32 of the 37 ARS-encoding loci cause severe, early-onset recessive phenotypes. Previous genetic and functional data suggest a loss-of-function mechanism; however, our understanding of the allelic and locus heterogeneity of ARS-related disease is incomplete. Cysteinyl-tRNA synthetase (CARS) encodes the enzyme that charges tRNA^Cys with cysteine in the cytoplasm. To date, CARS variants have not been implicated in any human disease phenotype. Here, we report on four subjects from three families with complex syndromes that include microcephaly, developmental delay, and brittle hair and nails. Each affected person carries bi-allelic CARS variants: one individual is compound heterozygous for c.1138C>T (p.Gln380^∗) and c.1022G>A (p.Arg341His), two related individuals are compound heterozygous for c.1076C>T (p.Ser359Leu) and c.1199T>A (p.Leu400Gln), and one individual is homozygous for c.2061dup (p.Ser688Glnfs^∗2). Measurement of protein abundance, yeast complementation assays, and assessments of tRNA charging indicate that each CARS variant causes a loss-of-function effect. Compared to subjects with previously reported ARS-related diseases, individuals with bi-allelic CARS variants are unique in presenting with a brittle-hair-and-nail phenotype, which most likely reflects the high cysteine content in human keratins. In sum, our efforts implicate CARS variants in human inherited disease, expand the locus and clinical heterogeneity of ARS-related clinical phenotypes, and further support impaired tRNA charging as the primary mechanism of recessive ARS-related disease.

Expanding the phenotype of phospholipid remodelling disease due to MBOAT7 gene defect

MBOAT7 gene codes O-acyltransferase domain containing seven proteins which is one of four enzymes involved in remodeling of phosphoinositol phosphate (PIP) in LANDs cycle. We present clinical, neuroimaging, and genetic findings of 12 patients from 7 families with MBOAT7 gene defect, a recently defined novel phospholipid remodelling disease. To the best of our knowledge, our case series is the second report on patients with MBOAT7 gene defect. The patients present with global developmental delay particularly in speech and language skills, intellectual disability, stereotypical behavior, ataxic gait, early onset epilepsy with well response to medical treatment, strabismus and similar facial features. Common neuroimaging findings of the patients were folium dysgenesis of the cerebellum with a particular appearance, mild-to-moderate cerebellar atrophy, T2 hyperintensity of bilateral globus pallidius and dentate nuclei, enlarged perivascular areas, and mild thinning of the corpus callosum. Genome-wide genotyping and exome sequencing identified five different types of homozygous mutations in the MBOAT7 gene in all seven families which are p.Arg87*, p.Leu227ProfsX65, p.Gln376Lys, p.Trp426*, and chr19:54.666.173-54.677.766/11594 bp del. We conclude that clinical and neuroimaging findings of MBOAT7 gene defect may suggest the diagnosis and guide genetic tests.

A suite of automated sequence analyses reduces the number of candidate deleterious variants and reveals a difference between probands and unaffected siblings

PURPOSE

Develop an automated exome analysis workflow that can produce a very small number of candidate variants yet still detect different numbers of deleterious variants between probands and unaffected siblings.

METHODS

Ninety-seven outbred nuclear families from the Undiagnosed Diseases Program/Network included single probands and the corresponding unaffected sibling(s). Single-nucleotide polymorphism (SNP) chip and exome analyses were performed on all, with proband and unaffected sibling considered independently as the target. The total burden of candidate genetic variants was summed for probands and siblings over all considered disease models.

RESULTS

Exome analysis workflow include automated programs for ethnicity-matched genotype calling, salvage pathway for Mendelian inconsistency, compound heterozygous recessive detection, BAM file regional curation, population frequency filtering, pedigree-aware BAM file noise evaluation, and exon deletion filtration. This workflow relied heavily on BAM file analysis. A greater average pathogenic variant number was found compared with unaffected siblings. This was significant (p < 0.05) when using published recommended thresholds, and implies that causal variants are retained in many probands' lists.

CONCLUSION

Using Mendelian and non-Mendelian models, this agnostic exome analysis shows a difference between a small group of probands and their unaffected siblings. This workflow produces candidate lists small enough to pursue with laboratory validation.

Biallelic sequence variants in INTS1 in patients with developmental delays, cataracts, and craniofacial anomalies

The Integrator complex subunit 1 (INTS1) is a component of the integrator complex that comprises 14 subunits and associates with RPB1 to catalyze endonucleolytic cleavage of nascent snRNAs and assist RNA polymerase II in promoter-proximal pause-release on protein-coding genes. We present five patients, including two sib pairs, with biallelic sequence variants in INTS1. The patients manifested absent or severely limited speech, an abnormal gait, hypotonia and cataracts. Exome sequencing revealed biallelic variants in INTS1 in all patients. One sib pair demonstrated a missense variant, p.(Arg77Cys), and a frameshift variant, p.(Arg1800Profs*20), another sib pair had a homozygous missense variant, p.(Pro1874Leu), and the fifth patient had a frameshift variant, p.(Leu1764Cysfs*16) and a missense variant, p.(Leu2164Pro). We also report additional clinical data on three previously described individuals with a homozygous, loss of function variant, p.(Ser1784*) in INTS1 that shared cognitive delays, cataracts and dysmorphic features with these patients. Several of the variants affected the protein C-terminus and preliminary modeling showed that the p.(Pro1874Leu) and p.(Leu2164Pro) variants may interfere with INTS1 helix folding. In view of the cataracts observed, we performed in-situ hybridization and demonstrated expression of ints1 in the zebrafish eye. We used Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)/Cas9 to make larvae with biallelic insertion/deletion (indel) variants in ints1. The mutant larvae developed typically through gastrulation, but sections of the eye showed abnormal lens development. The distinctive phenotype associated with biallelic variants in INTS1 points to dysfunction of the integrator complex as a mechanism for intellectual disability, eye defects and craniofacial anomalies.

Novel SPEG mutations in congenital myopathies: Genotype-phenotype correlations

INTRODUCTION

Centronuclear myopathies (CNMs) are a subtype of congenital myopathies (CMs) characterized by muscle weakness, predominant type 1 fibers, and increased central nuclei. SPEG (striated preferentially expressed protein kinase) mutations have recently been identified in 7 CM patients (6 with CNMs). We report 2 additional patients with SPEG mutations expanding the phenotype and evaluate genotype-phenotype correlations associated with SPEG mutations.

METHODS

Using whole exome/genome sequencing in CM families, we identified novel recessive SPEG mutations in 2 patients.

RESULTS

Patient 1, with severe muscle weakness requiring respiratory support, dilated cardiomyopathy, ophthalmoplegia, and findings of nonspecific CM on muscle biopsy carried a homozygous SPEG mutation (p.Val3062del). Patient 2, with milder muscle weakness, ophthalmoplegia, and CNM carried compound heterozygous mutations (p.Leu728Argfs*82) and (p.Val2997Glyfs*52).

CONCLUSIONS

The 2 patients add insight into genotype-phenotype correlations of SPEG-associated CMs. Clinicians should consider evaluating a CM patient for SPEG mutations even in the absence of CNM features. Muscle Nerve 59:357-362, 2019.

Joubert Syndrome: Ophthalmological Findings in Correlation with Genotype and Hepatorenal Disease in 99 Patients Prospectively Evaluated at a Single Center

PURPOSE

Joubert syndrome (JS) is caused by mutations in >34 genes that encode proteins involved with primary (nonmotile) cilia and the cilium basal body. This study describes the varying ocular phenotypes in JS patients, with correlation to systemic findings and genotype.

DESIGN

Patients were systematically and prospectively examined at the National Institutes of Health (NIH) Clinical Center in the setting of a dedicated natural history clinical trial.

PARTICIPANTS

Ninety-nine patients with JS examined at a single center.

METHODS

All patients underwent genotyping for JS, followed by complete age-appropriate ophthalmic examinations at the NIH Clinical Center, including visual acuity (VA), fixation behavior, lid position, motility assessment, slit-lamp biomicroscopy, dilated fundus examination with an indirect ophthalmoscope, and retinoscopy. Color and fundus autofluorescence imaging, Optos wide-field photography (Dunfermline, Scotland, UK), and electroretinography (ERG) were performed when possible.

MAIN OUTCOME MEASURES

The VA (with longitudinal follow-up where possible), ptosis, extraocular muscle function, retinal and optic nerve status, and retinal function as measured by ERG.

RESULTS

Among patients with JS with quantifiable VA (68/99), values ranged from 0 logarithm of the minimum angle of resolution (logMAR) (Snellen 20/20) to 1.5 logMAR (Snellen 20/632). Strabismus (71/98), nystagmus (66/99), oculomotor apraxia (60/77), ptosis (30/98), coloboma (28/99), retinal degeneration (20/83), and optic nerve atrophy (8/86) were identified.

CONCLUSIONS

We recommend regular monitoring for ophthalmological manifestations of JS beginning soon after birth or diagnosis. We demonstrate delayed visual development and note that the amblyogenic time frame may last significantly longer in JS than is typical. In general, patients with coloboma were less likely to display retinal degeneration, and those with retinal degeneration did not have coloboma. Severe retinal degeneration that is early and aggressive is seen in disease caused by specific genes, such as CEP290- and AHI1-associated JS. Retinal degeneration in INPP5E-, MKS1-, and NPHP1-associated JS was generally milder. Finally, ptosis surgery can be helpful in a subset of patients with JS; decisions as to timing and benefit/risk ratio need to be made on an individual basis according to expert consultation.

Variants in EXOSC9 Disrupt the RNA Exosome and Result in Cerebellar Atrophy with Spinal Motor Neuronopathy

The exosome is a conserved multi-protein complex that is essential for correct RNA processing. Recessive variants in exosome components EXOSC3, EXOSC8, and RBM7 cause various constellations of pontocerebellar hypoplasia (PCH), spinal muscular atrophy (SMA), and central nervous system demyelination. Here, we report on four unrelated affected individuals with recessive variants in EXOSC9 and the effect of the variants on the function of the RNA exosome in vitro in affected individuals' fibroblasts and skeletal muscle and in vivo in zebrafish. The clinical presentation was severe, early-onset, progressive SMA-like motor neuronopathy, cerebellar atrophy, and in one affected individual, congenital fractures of the long bones. Three affected individuals of different ethnicity carried the homozygous c.41T>C (p.Leu14Pro) variant, whereas one affected individual was compound heterozygous for c.41T>C (p.Leu14Pro) and c.481C>T (p.Arg161^∗). We detected reduced EXOSC9 in fibroblasts and skeletal muscle and observed a reduction of the whole multi-subunit exosome complex on blue-native polyacrylamide gel electrophoresis. RNA sequencing of fibroblasts and skeletal muscle detected significant >2-fold changes in genes involved in neuronal development and cerebellar and motor neuron degeneration, demonstrating the widespread effect of the variants. Morpholino oligonucleotide knockdown and CRISPR/Cas9-mediated mutagenesis of exosc9 in zebrafish recapitulated aspects of the human phenotype, as they have in other zebrafish models of exosomal disease. Specifically, portions of the cerebellum and hindbrain were absent, and motor neurons failed to develop and migrate properly. In summary, we show that variants in EXOSC9 result in a neurological syndrome combining cerebellar atrophy and spinal motoneuronopathy, thus expanding the list of human exosomopathies.

A Screening Approach to Identify Clinically Actionable Variants Causing Congenital Heart Disease in Exome Data

BACKGROUND

Congenital heart disease (CHD)-structural abnormalities of the heart that arise during embryonic development-is the most common inborn malformation, affecting ≤1% of the population. However, currently, only a minority of cases can be explained by genetic abnormalities. The goal of this study was to identify disease-causal genetic variants in 30 families affected by CHD.

METHODS

Whole-exome sequencing was performed with the DNA of multiple family members. We utilized a 2-tiered whole-exome variant screening and interpretation procedure. First, we manually curated a high-confidence list of 90 genes known to cause CHD in humans, identified predicted damaging variants in genes on this list, and rated their pathogenicity using American College of Medical Genetics and Genomics-Association for Molecular Pathology guidelines.

RESULTS

In 3 families (10%), we found pathogenic variants in known CHD genes TBX5, TFAP2B, and PTPN11, explaining the cardiac lesions. Second, exomes were comprehensively analyzed to identify additional predicted damaging variants that segregate with disease in CHD candidate genes. In 10 additional families (33%), likely disease-causal variants were uncovered in PBX1, CNOT1, ZFP36L2, TEK, USP34, UPF2, KDM5A, KMT2C, TIE1, TEAD2, and FLT4.

CONCLUSIONS

The pathogenesis of CHD could be explained using our high-confidence CHD gene list for variant filtering in a subset of cases. Furthermore, our unbiased screening procedure of family exomes implicates additional genes and variants in the pathogenesis of CHD, which suggest themselves for functional validation. This 2-tiered approach provides a means of (1) identifying clinically actionable variants and (2) identifying additional disease-causal genes, both of which are essential for improving the molecular diagnosis of CHD.

Characteristics of Liver Disease in 100 Individuals With Joubert Syndrome Prospectively Evaluated at a Single Center

BACKGROUND AND AIMS

Joubert Syndrome (JS) is a rare, inherited, ciliopathy defined by cerebellar and brainstem malformations and is variably associated with liver, kidney, and ocular dysfunction. This study characterizes the hepatic findings in JS and identifies factors associated with probable portal hypertension.

METHODS

Hundred individuals with JS were prospectively evaluated at the National Institutes of Health Clinical Center. Laboratory tests, imaging, and DNA sequencing were performed. Patients were stratified based on the spleen length/patient height ratio as a marker of splenomegaly, used as a surrogate for probable portal hypertension.

RESULTS

Forty-three patients (43%) had liver involvement based on elevated liver enzymes and/or liver hyperechogenicity and/or splenomegaly. None of the patients had macroscopic liver cysts or bile duct dilatation. Based on the spleen length/patient height ratio, 13 patients were stratified into a probable portal hypertension group. We observed significant elevations in alkaline phosphatase (269 vs 169 U/L, P ≤ 0.001), alanine aminotransferase (92 vs 42 U/L, P = 0.004), aspartate aminotransferase (77 vs 40 U/L, P = 0.002), and gamma-glutamyl transferase (226 vs 51 U/L, P ≤ 0.001) in the probable portal hypertension group. Platelets were lower in the probable portal hypertension cohort (229 vs 299 × 10 cells/μL, P = 0.008), whereas synthetic function was intact in both groups. Probable portal hypertension was also more prevalent in patients with kidney disease (P = 0.001) and colobomas (P = 0.02), as well as mutations in the TMEM67 gene (P = 0.001).

CONCLUSIONS

In JS, probable portal hypertension is associated with abnormal hepatic enzymes, as well as presence of kidney disease, coloboma, and/or mutation in TMEM67. These findings may allow early identification of JS patients who have or are more likely to develop liver disease.

High-Throughput Tabular Data Processor - Platform independent graphical tool for processing large data sets

High-throughput technologies generate considerable amount of data which often requires bioinformatic expertise to analyze. Here we present High-Throughput Tabular Data Processor (HTDP), a platform independent Java program. HTDP works on any character-delimited column data (e.g. BED, GFF, GTF, PSL, WIG, VCF) from multiple text files and supports merging, filtering and converting of data that is produced in the course of high-throughput experiments. HTDP can also utilize itemized sets of conditions from external files for complex or repetitive filtering/merging tasks. The program is intended to aid global, real-time processing of large data sets using a graphical user interface (GUI). Therefore, no prior expertise in programming, regular expression, or command line usage is required of the user. Additionally, no a priori assumptions are imposed on the internal file composition. We demonstrate the flexibility and potential of HTDP in real-life research tasks including microarray and massively parallel sequencing, i.e. identification of disease predisposing variants in the next generation sequencing data as well as comprehensive concurrent analysis of microarray and sequencing results. We also show the utility of HTDP in technical tasks including data merge, reduction and filtering with external criteria files. HTDP was developed to address functionality that is missing or rudimentary in other GUI software for processing character-delimited column data from high-throughput technologies. Flexibility, in terms of input file handling, provides long term potential functionality in high-throughput analysis pipelines, as the program is not limited by the currently existing applications and data formats. HTDP is available as the Open Source software (https://github.com/pmadanecki/htdp).

SIX3 deletions and incomplete penetrance in families affected by holoprosencephaly

Holoprosencephaly (HPE) is failure of the forebrain to divide completely during embryogenesis. Incomplete penetrance has not been reported previously in SIX3 whole gene deletions, which are known to cause HPE. Both chromosomal microarray and whole exome sequencing (WES) were used to evaluate families with inherited HPE. Two families showed inherited deletions that contain SIX3 and were incompletely penetrant for HPE. Using WES, we ruled out parental mosaicism, a SIX3 hypomorph, and clinically significant variants in genes that are known to interact with SIX3 as causes of incomplete penetrance. We demonstrate the importance of molecular cascade testing in families with HPE and we answer important questions about incomplete penetrance.

A novel inborn error of the coenzyme Q10 biosynthesis pathway: cerebellar ataxia and static encephalomyopathy due to COQ5 C-methyltransferase deficiency

Primary coenzyme Q10 (CoQ10 ; MIM# 607426) deficiencies are an emerging group of inherited mitochondrial disorders with heterogonous clinical phenotypes. Over a dozen genes are involved in the biosynthesis of CoQ10 , and mutations in several of these are associated with human disease. However, mutations in COQ5 (MIM# 616359), catalyzing the only C-methylation in the CoQ10 synthetic pathway, have not been implicated in human disease. Here, we report three female siblings of Iraqi-Jewish descent, who had varying degrees of cerebellar ataxia, encephalopathy, generalized tonic-clonic seizures, and cognitive disability. Whole-exome and subsequent whole-genome sequencing identified biallelic duplications in the COQ5 gene, leading to reduced levels of CoQ10 in peripheral white blood cells of all affected individuals and reduced CoQ10 levels in the only muscle tissue available from one affected proband. CoQ10 supplementation led to clinical improvement and increased the concentrations of CoQ10 in blood. This is the first report of primary CoQ10 deficiency caused by loss of function of COQ5, with delineation of the clinical, laboratory, histological, and molecular features, and insights regarding targeted treatment with CoQ10 supplementation.

Tumor exome sequencing and copy number alterations reveal potential predictors of intrinsic resistance to multi-targeted tyrosine kinase inhibitors

Multi-targeted tyrosine kinase inhibitors (TKIs) have broad efficacy and similar FDA-approved indications, suggesting shared molecular drug targets across cancer types. Irrespective of tumor type, 20-30% of patients treated with multi-targeted TKIs demonstrate intrinsic resistance, with progressive disease as a best response. We conducted a retrospective cohort study to identify tumor (somatic) point mutations, insertion/deletions, and copy number alterations (CNA) associated with intrinsic resistance to multi-targeted TKIs. Using a candidate gene approach (n=243), tumor next-generation sequencing and CNA data was associated with resistant and non-resistant outcomes. Resistant individuals (n=11) more commonly harbored somatic point mutations in NTRK1, KDR, TGFBR2, and PTPN11 and CNA in CDK4, CDKN2B, and ERBB2 compared to non-resistant (n=26, p<0.01). Using a random forest classification model for variable reduction and a decision tree classification model, we were able to differentiate intrinsically resistant from non-resistant patients. CNA in CDK4 and CDKN2B were the most important analytical features, implicating the cyclin D pathway as a potentially important factor in resistance to multi-targeted TKIs. Replication of these results in a larger, independent patient cohort has potential to inform personalized prescribing of these widely utilized agents.

De novo, deleterious sequence variants that alter the transcriptional activity of the homeoprotein PBX1 are associated with intellectual disability and pleiotropic developmental defects

We present eight patients with de novo, deleterious sequence variants in the PBX1 gene. PBX1 encodes a three amino acid loop extension (TALE) homeodomain transcription factor that forms multimeric complexes with TALE and HOX proteins to regulate target gene transcription during development. As previously reported, Pbx1 homozygous mutant mice (Pbx1-/-) develop malformations and hypoplasia or aplasia of multiple organs, including the craniofacial skeleton, ear, branchial arches, heart, lungs, diaphragm, gut, kidneys, and gonads. Clinical findings similar to those in Pbx mutant mice were observed in all patients with varying expressivity and severity, including external ear anomalies, abnormal branchial arch derivatives, heart malformations, diaphragmatic hernia, renal hypoplasia and ambiguous genitalia. All patients but one had developmental delays. Previously reported patients with congenital anomalies affecting the kidney and urinary tract exhibited deletions and loss of function variants in PBX1. The sequence variants in our cases included missense substitutions adjacent to the PBX1 homeodomain (p.Arg184Pro, p.Met224Lys, and p.Arg227Pro) or within the homeodomain (p.Arg234Pro, and p.Arg235Gln), whereas p.Ser262Glnfs*2, and p.Arg288* yielded truncated PBX1 proteins. Functional studies on five PBX1 sequence variants revealed perturbation of intrinsic, PBX-dependent transactivation ability and altered nuclear translocation, suggesting abnormal interactions between mutant PBX1 proteins and wild-type TALE or HOX cofactors. It is likely that the mutations directly affect the transcription of PBX1 target genes to impact embryonic development. We conclude that deleterious sequence variants in PBX1 cause intellectual disability and pleiotropic malformations resembling those in Pbx1 mutant mice, arguing for strong conservation of gene function between these two species.

Plasticity of the MFS1 Promoter Leads to Multidrug Resistance in the Wheat Pathogen Zymoseptoria tritici

The ascomycete Zymoseptoria tritici is the causal agent of Septoria leaf blotch on wheat. Disease control relies mainly on resistant wheat cultivars and on fungicide applications. The fungus displays a high potential to circumvent both methods. Resistance against all unisite fungicides has been observed over decades. A different type of resistance has emerged among wild populations with multidrug-resistant (MDR) strains. Active fungicide efflux through overexpression of the major facilitator gene MFS1 explains this emerging resistance mechanism. Applying a bulk-progeny sequencing approach, we identified in this study a 519-bp long terminal repeat (LTR) insert in the MFS1 promoter, a relic of a retrotransposon cosegregating with the MDR phenotype. Through gene replacement, we show the insert as a mutation responsible for MFS1 overexpression and the MDR phenotype. Besides this type I insert, we found two different types of promoter inserts in more recent MDR strains. Type I and type II inserts harbor potential transcription factor binding sites, but not the type III insert. Interestingly, all three inserts correspond to repeated elements present at different genomic locations in either IPO323 or other Z. tritici strains. These results underline the plasticity of repeated elements leading to fungicide resistance in Z. tritici and which contribute to its adaptive potential. IMPORTANCE Disease control through fungicides remains an important means to protect crops from fungal diseases and to secure the harvest. Plant-pathogenic fungi, especially Zymoseptoria tritici, have developed resistance against most currently used active ingredients, reducing or abolishing their efficacy. While target site modification is the most common resistance mechanism against single modes of action, active efflux of multiple drugs is an emerging phenomenon in fungal populations reducing additionally fungicides' efficacy in multidrug-resistant strains. We have investigated the mutations responsible for increased drug efflux in Z. tritici field strains. Our study reveals that three different insertions of repeated elements in the same promoter lead to multidrug resistance in Z. tritici. The target gene encodes the membrane transporter MFS1 responsible for drug efflux, with the promoter inserts inducing its overexpression. These results underline the plasticity of repeated elements leading to fungicide resistance in Z. tritici.

The FOXA2 transcription factor is frequently somatically mutated in uterine carcinosarcomas and carcinomas

BACKGROUND

Uterine carcinosarcomas (UCSs) are a rare but clinically aggressive form of cancer. They are biphasic tumors consisting of both epithelial and sarcomatous components. The majority of uterine carcinosarcomas are clonal, with the carcinomatous cells undergoing metaplasia to give rise to the sarcomatous component. The objective of the current study was to identify novel somatically mutated genes in UCSs.

METHODS

We whole exome sequenced paired tumor and nontumor DNAs from 14 UCSs and orthogonally validated 464 somatic variants using Sanger sequencing. Fifteen genes that were somatically mutated in at least 2 tumor exomes were Sanger sequenced in another 39 primary UCSs.

RESULTS

Overall, among 53 UCSs in the current study, the most frequently mutated of these 15 genes were tumor protein p53 (TP53) (75.5%), phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha (PIK3CA) (34.0%), protein phosphatase 2, regulatory subunit A, alpha (PPP2R1A) (18.9%), F-box and WD repeat domain containing 7 (FBXW7) (18.9%), chromodomain helicase DNA binding protein 4 (CHD4) (17.0%), and forkhead box A2 (FOXA2) (15.1%). FOXA2 has not previously been implicated in UCSs and was predominated by frameshift and nonsense mutations. One UCS with a FOXA2 frameshift mutation expressed truncated FOXA2 protein by immunoblotting. Sequencing of FOXA2 in 160 primary endometrial carcinomas revealed somatic mutations in 5.7% of serous, 22.7% of clear cell, 9% of endometrioid, and 11.1% of mixed endometrial carcinomas, the majority of which were frameshift mutations.

CONCLUSIONS

Collectively, the findings of the current study provide compelling genetic evidence that FOXA2 is a pathogenic driver gene in the etiology of primary uterine cancers, including UCSs. Cancer 2018;124:65-73. © 2017 American Cancer Society.

Compound heterozygosity for loss-of-function GARS variants results in a multisystem developmental syndrome that includes severe growth retardation

Aminoacyl-tRNA synthetases (ARSs) are ubiquitously expressed enzymes that ligate amino acids onto tRNA molecules. Genes encoding ARSs have been implicated in myriad dominant and recessive disease phenotypes. Glycyl-tRNA synthetase (GARS) is a bifunctional ARS that charges tRNA^Gly in the cytoplasm and mitochondria. GARS variants have been associated with dominant Charcot-Marie-Tooth disease but have not been convincingly implicated in recessive phenotypes. Here, we describe a patient from the NIH Undiagnosed Diseases Program with a multisystem, developmental phenotype. Whole-exome sequence analysis revealed that the patient is compound heterozygous for one frameshift (p.Glu83Ilefs*6) and one missense (p.Arg310Gln) GARS variant. Using in vitro and in vivo functional studies, we show that both GARS variants cause a loss-of-function effect: the frameshift variant results in depleted protein levels and the missense variant reduces GARS tRNA charging activity. In support of GARS variant pathogenicity, our patient shows striking phenotypic overlap with other patients having ARS-related recessive diseases, including features associated with variants in both cytoplasmic and mitochondrial ARSs; this observation is consistent with the essential function of GARS in both cellular locations. In summary, our clinical, genetic, and functional analyses expand the phenotypic spectrum associated with GARS variants.

Compound heterozygous alterations in intraflagellar transport protein CLUAP1 in a child with a novel Joubert and oral-facial-digital overlap syndrome

Disruption of normal ciliary function results in a range of diseases collectively referred to as ciliopathies. Here we report a child with a phenotype that overlapped with Joubert, oral–facial–digital, and Pallister–Hall syndromes including brain, limb, and craniofacial anomalies. We performed exome-sequence analysis on a proband and both parents, filtered for putative causative variants, and Sanger-verified variants of interest. Identified variants in CLUAP1 were functionally analyzed in a Xenopus system to determine their effect on ciliary function. Two variants in CLUAP1 were identified through exome-sequence analysis, Chr16:g.3558407T>G, c.338T>G, p.(Met113Arg) and Chr16:g.3570011C>T, c.688C>T, p.(Arg230Ter). These variants were rare in the Exome Aggregation Consortium (ExAC) data set of 65,000 individuals (one and two occurrences, respectively). Transfection of mutant CLUAP1 constructs into Xenopus embryos showed reduced protein levels p.(Arg230Ter) and reduced intraflagellar transport p.(Met113Arg). The genetic data show that these variants are present in an affected child, are rare in the population, and result in reduced, but not absent, intraflagellar transport. We conclude that biallelic mutations in CLUAP1 resulted in this novel ciliopathy syndrome in the proband.

Myophosphorylase (PYGM) mutations determined by next generation sequencing in a cohort from Turkey with McArdle disease

This study aimed to identify PYGM mutations in patients with McArdle disease from Turkey by next generation sequencing (NGS). Genomic DNA was extracted from the blood of the McArdle patients (n = 67) and unrelated healthy volunteers (n = 53). The PYGM gene was sequenced with NGS and the observed mutations were validated by direct Sanger sequencing. A diagnostic algorithm was developed for patients with suspected McArdle disease. A total of 16 deleterious PYGM mutations were identified, of which 5 were novel, including 1 splice-site donor, 1 frame-shift, and 3 non-synonymous variants. The p.Met1Val (27-patients/11-families) was the most common PYGM mutation, followed by p.Arg576* (6/4), c.1827+7A>G (5/4), c.772+2_3delTG (5/3), p.Phe710del (4/2), p.Lys754Asnfs (2/1), and p.Arg50* (1/1). A molecular diagnostic flowchart is proposed for the McArdle patients in Turkey, covering the 6 most common PYGM mutations found in Turkey as well as the most common mutation in Europe. The diagnostic algorithm may alleviate the need for muscle biopsies in 77.6% of future patients. A prevalence of any of the mutations to a geographical region in Turkey was not identified. Furthermore, the NGS approach to sequence the entire PYGM gene was successful in detecting a common missense mutation and discovering novel mutations in this population study.

P4HA1 mutations cause a unique congenital disorder of connective tissue involving tendon, bone, muscle and the eye

Collagen prolyl 4-hydroxylases (C-P4Hs) play a central role in the formation and stabilization of the triple helical domain of collagens. P4HA1 encodes the catalytic α(I) subunit of the main C-P4H isoenzyme (C-P4H-I). We now report human bi-allelic P4HA1 mutations in a family with a congenital-onset disorder of connective tissue, manifesting as early-onset joint hypermobility, joint contractures, muscle weakness and bone dysplasia as well as high myopia, with evidence of clinical improvement of motor function over time in the surviving patient. Similar to P4ha1 null mice, which die prenatally, the muscle tissue from P1 and P2 was found to have reduced collagen IV immunoreactivity at the muscle basement membrane. Patients were compound heterozygous for frameshift and splice site mutations leading to reduced, but not absent, P4HA1 protein level and C-P4H activity in dermal fibroblasts compared to age-matched control samples. Differential scanning calorimetry revealed reduced thermal stability of collagen in patient-derived dermal fibroblasts versus age-matched control samples. Mutations affecting the family of C-P4Hs, and in particular C-P4H-I, should be considered in patients presenting with congenital connective tissue/myopathy overlap disorders with joint hypermobility, contractures, mild skeletal dysplasia and high myopia.

Cytoplasmic body pathology in severe ACTA1-related myopathy in the absence of typical nemaline rods

Mutations in ACTA1 cause a group of myopathies with expanding clinical and histopathological heterogeneity. We describe three patients with severe ACTA1-related myopathy who have muscle fiber cytoplasmic bodies but no classic nemaline rods. Patient 1 is a five-year-old boy who presented at birth with severe weakness and respiratory failure, requiring mechanical ventilation. Whole exome sequencing identified a heterozygous c.282C>A (p.Asn94Lys) ACTA1 mutation. Patients 2 and 3 were twin boys with hypotonia, severe weakness, and respiratory insufficiency at birth requiring mechanical ventilation. Both died at 6 months of age. The same heterozygous c.282C>A (p.Asn94Lys) ACTA1 mutation was identified by whole exome sequencing. We conclude that clinically severe ACTA1-related myopathy can present with muscle morphological findings suggestive of cytoplasmic body myopathy in the absence of definite nemaline rods. The Asn94Lys mutation in skeletal muscle sarcomeric α-actin may be linked to this histological appearance. These novel ACTA1 cases also illustrate the successful application of whole exome sequencing in directly arriving at a candidate genetic diagnosis in patients with unexpected phenotypic and histologic features for a known neuromuscular gene.

HSP and deafness: Neurocristopathy caused by a novel mosaic SOX10 mutation

OBJECTIVE

To identify the underlying genetic cause in 2 sisters affected with progressive lower extremity spasticity, neuropathy, and early-onset deafness.

METHODS

Whole-exome sequencing was performed, and segregation testing of variants was investigated using targeted Sanger sequencing. An inherited paternal mosaic mutation was further evaluated through quantitative analysis of the ratio of mutant vs wild-type allele in genomic DNA from various tissues, including blood, dermal fibroblasts, and saliva.

RESULTS

A novel heterozygous nonsense mutation (c.1140C>A; p.Y380X) in SOX10 was identified in the affected sisters. Paternal mosaicism was suspected based on a small chromatogram peak, which was less than the heterozygous peak of the mutated allele. Consistent with mosaicism, the mosaic paternal samples had notable variability in the ratio of mutant vs wild-type allele in various tissues (compared with the fully heterozygous daughter), with the highest paternal mutant levels in saliva (32.7%) and lowest in dermal fibroblasts (13.9%). Targeted clinical re-examination of the father revealed a sensorimotor neuropathy that was previously clinically unrecognized.

CONCLUSIONS

These findings expand the phenotypic spectrum of SOX10-related neurocristopathy. Mutations in SOX10 should be considered in patients presenting with a complicated form of hereditary spastic paraplegia that includes neuropathy and deafness. Diagnostic workup may be complicated, as SOX10 mutations can present in a mosaic state, with a mild clinical manifestation.

Somatic mutation profiles of clear cell endometrial tumors revealed by whole exome and targeted gene sequencing

BACKGROUND

The molecular pathogenesis of clear cell endometrial cancer (CCEC), a tumor type with a relatively unfavorable prognosis, is not well defined. We searched exome-wide for novel somatically mutated genes in CCEC and assessed the mutational spectrum of known and candidate driver genes in a large cohort of cases.

METHODS

We conducted whole exome sequencing of paired tumor-normal DNAs from 16 cases of CCEC (12 CCECs and the CCEC components of 4 mixed histology tumors). Twenty-two genes-of-interest were Sanger-sequenced from another 47 cases of CCEC. Microsatellite instability (MSI) and microsatellite stability (MSS) were determined by genotyping 5 mononucleotide repeats.

RESULTS

Two tumor exomes had relatively high mutational loads and MSI. The other 14 tumor exomes were MSS and had 236 validated nonsynonymous or splice junction somatic mutations among 222 protein-encoding genes. Among the 63 cases of CCEC in this study, we identified frequent somatic mutations in TP53 (39.7%), PIK3CA (23.8%), PIK3R1 (15.9%), ARID1A (15.9%), PPP2R1A (15.9%), SPOP (14.3%), and TAF1 (9.5%), as well as MSI (11.3%). Five of 8 mutations in TAF1, a gene with no known role in CCEC, localized to the putative histone acetyltransferase domain and included 2 recurrently mutated residues. Based on patterns of MSI and mutations in 7 genes, CCEC subsets molecularly resembled serous endometrial cancer (SEC) or endometrioid endometrial cancer (EEC).

CONCLUSION

Our findings demonstrate molecular similarities between CCEC and SEC and EEC and implicate TAF1 as a novel candidate CCEC driver gene. Cancer 2017;123:3261-8. © 2017 American Cancer Society.

Clinical and molecular phenotyping of a child with Hermansky-Pudlak syndrome-7, an uncommon genetic type of HPS

PURPOSE

Hermansky-Pudlak syndrome (HPS) is a rare inherited disorder with ten reported genetic types; each type has defects in subunits of either Adaptor Protein-3 complex or Biogenesis of Lysosome-related Organelles Complex (BLOC)-1, -2, or -3. Very few patients with BLOC-1 deficiency (HPS-7, -8, and -9 types) have been diagnosed. We report results of comprehensive clinical testing and molecular analyses of primary fibroblasts from a new case of HPS-7.

RESULTS

A 6-year old Paraguayan male presented with hypopigmentation, ocular albinism, nystagmus, reduced visual acuity, and easy bruising. He also experienced delayed motor and language development as a very young child; head and chest trauma resulted in intracranial hemorrhage with subsequent right hemiparesis and lung scarring. There was no clinical evidence of immunodeficiency or colitis. Whole mount transmission electron microscopy revealed absent platelet delta granules; platelet aggregation testing was abnormal. Exome sequencing revealed a homozygous nonsense mutation in the Dystrobrevin binding protein 1 (DTNBP1) gene [NM_032122.4: c.307C>T; p.Gln103*], previously reported in a Portuguese adult. The gene encodes the dysbindin subunit of BLOC-1. Dysbindin protein expression was negligible in our patient's dermal fibroblasts, while his DTNBP1 mRNA level was similar to that of a normal control.

CONCLUSIONS

Comprehensive clinical evaluation of the first pediatric case reported with HPS-7 reveals oculocutaneous albinism and platelet storage pool deficiency; his phenotype is consistent with findings in other patients with BLOC-1 disorders. This patient's markedly reduced Dysbindin protein expression in HPS-7 resulted from a mechanism other than nonsense mediated decay.

Recessive coding and regulatory mutations in FBLIM1 underlie the pathogenesis of chronic recurrent multifocal osteomyelitis (CRMO)

Chronic recurrent multifocal osteomyelitis (CRMO) is a rare, pediatric, autoinflammatory disease characterized by bone pain due to sterile osteomyelitis, and is often accompanied by psoriasis or inflammatory bowel disease. There are two syndromic forms of CRMO, Majeed syndrome and DIRA, for which the genetic cause is known. However, for the majority of cases of CRMO, the genetic basis is unknown. Via whole-exome sequencing, we detected a homozygous mutation in the filamin-binding domain of FBLIM1 in an affected child with consanguineous parents. Microarray analysis of bone marrow macrophages from the CRMO murine model (cmo) determined that the Fblim1 ortholog is the most differentially expressed gene, downregulated over 20-fold in the cmo mouse. We sequenced FBLIM1 in 96 CRMO subjects and found a second proband with a novel frameshift mutation in exon 6 and a rare regulatory variant. In SaOS2 cells, overexpressing the regulatory mutation showed the flanking region acts as an enhancer, and the mutation ablates enhancer activity. Our data implicate FBLIM1 in the pathogenesis of sterile bone inflammation and our findings suggest CRMO is a disorder of chronic inflammation and imbalanced bone remodeling.