Search

Abstract Background Affymetrix 3' GeneChip microarrays are widely used to profile the expression of thousands of genes simultaneously. They differ from many other microarray types in that GeneChips are hybridised using a single labelled extract and because they contain multiple 'match' and 'mismatch' sequences for each transcript. Most algorithms extract the signal from GeneChip experiments in a sequence of separate steps, including background correction and normalisation, which inhibits the simultaneous use of all available information. They principally provide a point estimate of gene expression and, in contrast to BGX, do not fully integrate the uncertainty arising from potentially heterogeneous responses of the probes. Results BGX is a new Bioconductor R package that implements an integrated Bayesian approach to the analysis of 3' GeneChip data. The software takes into account additive and multiplicative error, non-specific hybridisation and replicate summarisation in the spirit of the model outlined in 1. It also provides a posterior distribution for the expression of each gene. Moreover, BGX can take into account probe affinity effects from probe sequence information where available. The package employs a novel adaptive Markov chain Monte Carlo (MCMC) algorithm that raises considerably the efficiency with which the posterior distributions are sampled from. Finally, BGX incorporates various ways to analyse the results, such as ranking genes by expression level as well as statistically based methods for estimating the amount of up and down regulated genes between two conditions. Conclusion BGX performs well relative to other widely used methods at estimating expression levels and fold changes. It has the advantage that it provides a statistically sound measure of uncertainty for its estimates. BGX includes various analysis functions to visualise and exploit the rich output that is produced by the Bayesian model.

Genome-wide platelet transcriptomics is increasingly used to uncover new aspects of platelet biology and as a diagnostic and prognostic tool. Nevertheless, platelet isolation methods for transcriptomic studies are not standardized, introducing challenges for cross-study comparisons, data integration, and replication. In this prospective multicenter study, called "Standardizing Platelet Transcriptomics for Discovery, Diagnostics, and Therapeutics in the Thrombosis and Hemostasis Community (STRIDE)" by the International Society on Thrombosis and Haemostasis Scientific and Standardization Committees, we assessed how 3 of the most commonly used platelet isolation protocols influence metrics from next-generation bulk RNA sequencing and functional assays. Compared with washing alone, more stringent removal of leukocytes by anti-CD45 beads or PALL filters resulted in a sufficient quantity of RNA for next-generation sequencing and similar quality of RNA sequencing metrics. Importantly, stringent removal of leukocytes resulted in the lower relative expression of known leukocyte-specific genes and the higher relative expression of known platelet-specific genes. The results were consistent across enrolling sites, suggesting that the techniques are transferrable and reproducible. Moreover, all 3 isolation techniques did not influence basal platelet reactivity, but agonist-induced integrin αIIbβ 3 activation is reduced by anti-CD45 bead isolation compared with washing alone. In conclusion, the isolation technique chosen influences genome-wide transcriptional and functional assays in platelets. These results should help the research community make informed choices about platelet isolation techniques in their own platelet studies., Competing Interests: Declaration of competing interests M.T.R. reports patents pending or issued on using platelet transcriptomics. All other authors have no conflicts to declare., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Members of the leucine rich repeat (LRR) and PDZ domain (LAP) protein family are essential for animal development and histogenesis. Densin-180, encoded by LRRC7, is the only LAP protein selectively expressed in neurons. Densin-180 is a postsynaptic scaffold at glutamatergic synapses, linking cytoskeletal elements with signalling proteins such as the α-subunit of Ca 2+ /calmodulin-dependent protein kinase II. We have previously observed an association between high impact variants in LRRC7 and Intellectual Disability; also three individual cases with variants in LRRC7 had been described. We identify here 33 individuals (one of them previously described) with a dominant neurodevelopmental disorder due to heterozygous missense or loss-of-function variants in LRRC7. The clinical spectrum involves intellectual disability, autism, ADHD, aggression and, in several cases, hyperphagia-associated obesity. A PDZ domain variant interferes with synaptic targeting of Densin-180 in primary cultured neurons. Using in vitro systems (two hybrid, BioID, coimmunoprecipitation of tagged proteins from 293T cells) we identified new candidate interaction partners for the LRR domain, including protein phosphatase 1 (PP1), and observed that variants in the LRR reduced binding to these proteins. We conclude that LRRC7 encodes a major determinant of intellectual development and behaviour., (© 2024. The Author(s).)

The major spliceosome comprises the five snRNAs U1, U2, U4, U5 and U6. We recently showed that mutations in RNU4- 2, which encodes U4 snRNA, cause one of the most prevalent monogenic neurodevelopmental disorders. Here, we report that recurrent germline mutations in RNU2-2P , a 191bp gene encoding U2 snRNA, are responsible for a related disorder. By genetic association, we implicated recurrent de novo single nucleotide mutations at nucleotide positions 4 and 35 of RNU2-2P among nine cases. We replicated this finding in six additional cases, bringing the total to 15. The disorder is characterized by intellectual disability, neurodevelopmental delay, autistic behavior, microcephaly, hypotonia, epilepsy and hyperventilation. All cases display a severe and complex seizure phenotype. Our findings cement the role of major spliceosomal snRNAs in the etiologies of neurodevelopmental disorders.

Most people with intellectual disability (ID) do not receive a molecular diagnosis following genetic testing. To identify new etiologies of ID, we performed a genetic association analysis comparing the burden of rare variants in 41,132 noncoding genes between 5,529 unrelated cases and 46,401 unrelated controls. RNU4-2, which encodes U4 small nuclear RNA, a critical component of the spliceosome, was the most strongly associated gene. We implicated de novo variants among 47 cases in two regions of RNU4-2 in the etiology of a syndrome characterized by ID, microcephaly, short stature, hypotonia, seizures and motor delay. We replicated this finding in three collections, bringing the number of unrelated cases to 73. Analysis of national genomic diagnostic data showed RNU4-2 to be a more common etiological gene for neurodevelopmental abnormality than any previously reported autosomal gene. Our findings add to the growing evidence of spliceosome dysfunction in the etiologies of neurological disorders., (© 2024. The Author(s).)

Individuals with ultrarare disorders pose a structural challenge for healthcare systems since expert clinical knowledge is required to establish diagnoses. In TRANSLATE NAMSE, a 3-year prospective study, we evaluated a novel diagnostic concept based on multidisciplinary expertise in Germany. Here we present the systematic investigation of the phenotypic and molecular genetic data of 1,577 patients who had undergone exome sequencing and were partially analyzed with next-generation phenotyping approaches. Molecular genetic diagnoses were established in 32% of the patients totaling 370 distinct molecular genetic causes, most with prevalence below 1:50,000. During the diagnostic process, 34 novel and 23 candidate genotype-phenotype associations were identified, mainly in individuals with neurodevelopmental disorders. Sequencing data of the subcohort that consented to computer-assisted analysis of their facial images with GestaltMatcher could be prioritized more efficiently compared with approaches based solely on clinical features and molecular scores. Our study demonstrates the synergy of using next-generation sequencing and phenotyping for diagnosing ultrarare diseases in routine healthcare and discovering novel etiologies by multidisciplinary teams., (© 2024. The Author(s).)

The etiologies of newborn deaths in neonatal intensive care units usually remain unknown, even after genetic testing. Whole-genome sequencing, combined with artificial intelligence-based methods for predicting the effects of non-coding variants, provide an avenue for resolving these deaths. Using one such method, SpliceAI, we identified a maternally inherited deep intronic PKHD1 splice variant (chr6:52030169T>C), in trans with a pathogenic missense variant (p.Thr36Met), in a newborn who died of autosomal recessive polycystic kidney disease at age 2 days. We validated the deep intronic variant's impact in maternal urine-derived cells expressing PKHD1. Reverse transcription polymerase chain reaction followed by Sanger sequencing showed that the variant causes inclusion of 147bp of the canonical intron between exons 29 and 30 of PKHD1 into the mRNA, including a premature stop codon. Allele-specific expression analysis at a heterozygous site in the mother showed that the mutant allele completely suppresses canonical splicing. In an unrelated healthy control, there was no evidence of transcripts including the novel splice junction. We returned a diagnostic report to the parents, who underwent in vitro embryo selection., (Copyright © 2024 National Kidney Foundation, Inc. Published by Elsevier Inc. All rights reserved.)

Purpose: Pathogenic variants of FIG4 generate enlarged lysosomes and neurological and developmental disorders. To identify additional genes regulating lysosomal volume, we carried out a genome-wide activation screen to detect suppression of enlarged lysosomes in FIG4 -/- cells., Methods: The CRISPR-a gene activation screen utilized sgRNAs from the promoters of protein-coding genes. Fluorescence-activated cell sorting separated cells with correction of the enlarged lysosomes from uncorrected cells. Patient variants of SLC12A9 were identified by exome or genome sequencing and studied by segregation analysis and clinical characterization., Results: Overexpression of SLC12A9, a solute co-transporter, corrected lysosomal swelling in FIG4 -/- cells. SLC12A9 (NP_064631.2) colocalized with LAMP2 at the lysosome membrane. Biallelic variants of SLC12A9 were identified in 3 unrelated probands with neurodevelopmental disorders. Common features included intellectual disability, skeletal and brain structural abnormalities, congenital heart defects, and hypopigmented hair. Patient 1 was homozygous for nonsense variant p.(Arg615∗), patient 2 was compound heterozygous for p.(Ser109Lysfs∗20) and a large deletion, and proband 3 was compound heterozygous for p.(Glu290Glyfs∗36) and p.(Asn552Lys). Fibroblasts from proband 1 contained enlarged lysosomes that were corrected by wild-type SLC12A9 cDNA. Patient variant p.(Asn552Lys) failed to correct the lysosomal defect., Conclusion: Impaired function of SLC12A9 results in enlarged lysosomes and a recessive disorder with a recognizable neurodevelopmental phenotype., Competing Interests: Conflict of Interest The authors declare no conflicts of interest., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Genetic studies of platelet reactivity (PR) phenotypes may identify novel antiplatelet drug targets. However, such studies have been limited by small sample sizes (n < 5000) because of the complexity of measuring PR. We trained a model to predict PR from complete blood count (CBC) scattergrams. A genome-wide association study of this phenotype in 29 806 blood donors identified 21 distinct associations implicating 20 genes, of which 6 have been identified previously. The effect size estimates were significantly correlated with estimates from a study of flow cytometry-measured PR and a study of a phenotype of in vitro thrombus formation. A genetic score of PR built from the 21 variants was associated with the incidence rates of myocardial infarction and pulmonary embolism. Mendelian randomization analyses showed that PR was causally associated with the risks of coronary artery disease, stroke, and venous thromboembolism. Our approach provides a blueprint for using phenotype imputation to study the determinants of hard-to-measure but biologically important hematological traits., (© 2023 by The American Society of Hematology. Licensed under Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International (CC BY-NC-ND 4.0), permitting only noncommercial, nonderivative use with attribution. All other rights reserved.)

The genetic etiologies of more than half of rare diseases remain unknown. Standardized genome sequencing and phenotyping of large patient cohorts provide an opportunity for discovering the unknown etiologies, but this depends on efficient and powerful analytical methods. We built a compact database, the 'Rareservoir', containing the rare variant genotypes and phenotypes of 77,539 participants sequenced by the 100,000 Genomes Project. We then used the Bayesian genetic association method BeviMed to infer associations between genes and each of 269 rare disease classes assigned by clinicians to the participants. We identified 241 known and 19 previously unidentified associations. We validated associations with ERG, PMEPA1 and GPR156 by searching for pedigrees in other cohorts and using bioinformatic and experimental approaches. We provide evidence that (1) loss-of-function variants in the Erythroblast Transformation Specific (ETS)-family transcription factor encoding gene ERG lead to primary lymphoedema, (2) truncating variants in the last exon of transforming growth factor-β regulator PMEPA1 result in Loeys-Dietz syndrome and (3) loss-of-function variants in GPR156 give rise to recessive congenital hearing impairment. The Rareservoir provides a lightweight, flexible and portable system for synthesizing the genetic and phenotypic data required to study rare disease cohorts with tens of thousands of participants., (© 2023. The Author(s).)

Transcriptional profiling of hematopoietic cell subpopulations has helped to characterize the developmental stages of the hematopoietic system and the molecular bases of malignant and non-malignant blood diseases. Previously, only the genes targeted by expression microarrays could be profiled genome-wide. High-throughput RNA sequencing, however, encompasses a broader repertoire of RNA molecules, without restriction to previously annotated genes. We analyzed the BLUEPRINT consortium RNA-sequencing data for mature hematopoietic cell types. The data comprised 90 total RNA-sequencing samples, each composed of one of 27 cell types, and 32 small RNA-sequencing samples, each composed of one of 11 cell types. We estimated gene and isoform expression levels for each cell type using existing annotations from Ensembl. We then used guided transcriptome assembly to discover unannotated transcripts. We identified hundreds of novel non-coding RNA genes and showed that the majority have cell type-dependent expression. We also characterized the expression of circular RNA and found that these are also cell type-specific. These analyses refine the active transcriptional landscape of mature hematopoietic cells, highlight abundant genes and transcriptional isoforms for each blood cell type, and provide a valuable resource for researchers of hematologic development and diseases. Finally, we made the data accessible via a web-based interface: https://blueprint.haem.cam.ac.uk/bloodatlas/.

Diagnosing mitochondrial disorders remains challenging. This is partly because the clinical phenotypes of patients overlap with those of other sporadic and inherited disorders. Although the widespread availability of genetic testing has increased the rate of diagnosis, the combination of phenotypic and genetic heterogeneity still makes it difficult to reach a timely molecular diagnosis with confidence. An objective, systematic method for describing the phenotypic spectra for each variant provides a potential solution to this problem. We curated the clinical phenotypes of 6688 published individuals with 89 pathogenic mitochondrial DNA (mtDNA) mutations, collating 26 348 human phenotype ontology (HPO) terms to establish the MitoPhen database. This enabled a hypothesis-free definition of mtDNA clinical syndromes, an overview of heteroplasmy-phenotype relationships, the identification of under-recognized phenotypes, and provides a publicly available reference dataset for objective clinical comparison with new patients using the HPO. Studying 77 patients with independently confirmed positive mtDNA diagnoses and 1083 confirmed rare disease cases with a non-mitochondrial nuclear genetic diagnosis, we show that HPO-based phenotype similarity scores can distinguish these two classes of rare disease patients with a false discovery rate <10% at a sensitivity of 80%. Enriching the MitoPhen database with more patients will improve predictions for increasingly rare variants., (© The Author(s) 2021. Published by Oxford University Press on behalf of Nucleic Acids Research.)

We have identified a rare missense variant on chromosome 9, position 125145990 (GRCh37), in exon 8 in PTGS1 (the gene encoding cyclo-oxygenase 1, COX-1, the target of anti-thrombotic aspirin therapy). We report that in the homozygous state within a large consanguineous family this variant is associated with a bleeding phenotype and alterations in platelet reactivity and eicosanoid production. Western blotting and confocal imaging demonstrated that COX-1 was absent in the platelets of three family members homozygous for the PTGS1 variant but present in their leukocytes. Platelet reactivity, as assessed by aggregometry, lumi-aggregometry and flow cytometry, was impaired in homozygous family members, as were platelet adhesion and spreading. The productions of COX-derived eicosanoids by stimulated platelets were greatly reduced but there were no changes in the levels of urinary metabolites of COX-derived eicosanoids. The proband exhibited additional defects in platelet aggregation and spreading which may explain why her bleeding phenotype was slightly more severe than those of other homozygous affected relatives. This is the first demonstration in humans of the specific loss of platelet COX-1 activity and provides insight into its consequences for platelet function and eicosanoid metabolism. Notably despite the absence of thromboxane A2 (TXA2) formation by platelets, urinary TXA2 metabolites were in the normal range indicating these cannot be assumed as markers of in vivo platelet function. Results from this study are important benchmarks for the effects of aspirin upon platelet COX-1, platelet function and eicosanoid production as they define selective platelet COX-1 ablation within humans.

Background: NRG1 gene fusions may be clinically actionable, since cancers carrying the fusion transcripts can be sensitive to tyrosine kinase inhibitors. The NRG1 gene encodes ligands for the HER2(ERBB2)-ERBB3 heterodimeric receptor tyrosine kinase, and the gene fusions are thought to lead to autocrine stimulation of the receptor. The NRG1 fusion expressed in the breast cancer cell line MDA-MB-175 serves as a model example of such fusions, showing the proposed autocrine loop and exceptional drug sensitivity. However, its structure has not been properly characterised, its oncogenic activity has not been fully explained, and there is limited data on such fusions in breast cancer., Methods: We analysed genomic rearrangements and transcripts of NRG1 in MDA-MB-175 and a panel of 571 breast cancers., Results: We found that the MDA-MB-175 fusion-originally reported as a DOC4(TENM4)-NRG1 fusion, lacking the cytoplasmic tail of NRG1-is in reality a double fusion, PPP6R3-TENM4-NRG1, producing multiple transcripts, some of which include the cytoplasmic tail. We hypothesise that many NRG1 fusions may be oncogenic not for lacking the cytoplasmic domain but because they do not encode NRG1's nuclear-localised form. The fusion in MDA-MB-175 is the result of a very complex genomic rearrangement, which we partially characterised, that creates additional expressed gene fusions, RSF1-TENM4, TPCN2-RSF1, and MRPL48-GAB2. We searched for NRG1 rearrangements in 571 breast cancers subjected to genome sequencing and transcriptome sequencing and found four cases (0.7%) with fusions, WRN-NRG1, FAM91A1-NRG1, ARHGEF39-NRG1, and ZNF704-NRG1, all splicing into NRG1 at the same exon as in MDA-MB-175. However, the WRN-NRG1 and ARHGEF39-NRG1 fusions were out of frame. We identified rearrangements of NRG1 in many more (8% of) cases that seemed more likely to inactivate than to create activating fusions, or whose outcome could not be predicted because they were complex, or both. This is not surprising because NRG1 can be pro-apoptotic and is inactivated in some breast cancers., Conclusions: Our results highlight the complexity of rearrangements of NRG1 in breast cancers and confirm that some do not activate but inactivate. Careful interpretation of NRG1 rearrangements will therefore be necessary for appropriate patient management.

Gray platelet syndrome (GPS) is a rare recessive disorder caused by biallelic variants in NBEAL2 and characterized by bleeding symptoms, the absence of platelet α-granules, splenomegaly, and bone marrow (BM) fibrosis. Due to the rarity of GPS, it has been difficult to fully understand the pathogenic processes that lead to these clinical sequelae. To discern the spectrum of pathologic features, we performed a detailed clinical genotypic and phenotypic study of 47 patients with GPS and identified 32 new etiologic variants in NBEAL2. The GPS patient cohort exhibited known phenotypes, including macrothrombocytopenia, BM fibrosis, megakaryocyte emperipolesis of neutrophils, splenomegaly, and elevated serum vitamin B12 levels. Novel clinical phenotypes were also observed, including reduced leukocyte counts and increased presence of autoimmune disease and positive autoantibodies. There were widespread differences in the transcriptome and proteome of GPS platelets, neutrophils, monocytes, and CD4 lymphocytes. Proteins less abundant in these cells were enriched for constituents of granules, supporting a role for Nbeal2 in the function of these organelles across a wide range of blood cells. Proteomic analysis of GPS plasma showed increased levels of proteins associated with inflammation and immune response. One-quarter of plasma proteins increased in GPS are known to be synthesized outside of hematopoietic cells, predominantly in the liver. In summary, our data show that, in addition to the well-described platelet defects in GPS, there are immune defects. The abnormal immune cells may be the drivers of systemic abnormalities such as autoimmune disease., (© 2020 by The American Society of Hematology.)

Hereditary hemorrhagic telangiectasia (HHT) is an autosomal dominant vascular dysplasia. Care delivery for HHT patients is impeded by the need for laborious, repeated phenotyping and gaps in knowledge regarding the relationships between causal DNA variants in ENG, ACVRL1, SMAD4 and GDF2, and clinical manifestations. To address this, we analyzed DNA samples from 183 previously uncharacterized, unrelated HHT and suspected HHT cases using the ThromboGenomics high-throughput sequencing platform. We identified 127 rare variants across 168 heterozygous genotypes. Applying modified American College of Medical Genetics and Genomics Guidelines, 106 variants were classified as pathogenic/likely pathogenic and 21 as nonpathogenic (variant of uncertain significance/benign). Unlike the protein products of ACVRL1 and SMAD4, the extracellular ENG amino acids are not strongly conserved. Our inferences of the functional consequences of causal variants in ENG were therefore informed by the crystal structure of endoglin. We then compared the accuracy of predictions of the causal gene blinded to the genetic data using 2 approaches: subjective clinical predictions and statistical predictions based on 8 Human Phenotype Ontology terms. Both approaches had some predictive power, but they were insufficiently accurate to be used clinically, without genetic testing. The distributions of red cell indices differed by causal gene but not sufficiently for clinical use in isolation from genetic data. We conclude that parallel sequencing of the 4 known HHT genes, multidisciplinary team review of variant calls in the context of detailed clinical information, and statistical and structural modeling improve the prognostication and treatment of HHT., (© 2020 by The American Society of Hematology.)