BACKGROUND: Cancer immunotherapies including immune checkpoint inhibitors and Chimeric Antigen Receptor (CAR) T-cell therapy have shown variable response rates in paediatric patients highlighting the need to establish robust biomarkers for patient selection. While the tumour microenvironment in adults has been widely studied to delineate determinants of immune response, the immune composition of paediatric solid tumours remains relatively uncharacterized calling for investigations to identify potential immune biomarkers. METHODS: To inform immunotherapy approaches in paediatric cancers with embryonal origin, we performed an immunogenomic analysis of RNA-seq data from 925 treatment-naïve paediatric nervous system tumours (pedNST) spanning 12 cancer types from three publicly available data sets. RESULTS: Within pedNST, we uncovered four broad immune clusters: Paediatric Inflamed (10%), Myeloid Predominant (30%), Immune Neutral (43%) and Immune Desert (17%). We validated these clusters using immunohistochemistry, methylation immune inference and segmentation analysis of tissue images. We report shared biology of these immune clusters within and across cancer types, and characterization of specific immune cell frequencies as well as T- and B-cell repertoires. We found no associations between immune infiltration levels and tumour mutational burden, although molecular cancer entities were enriched within specific immune clusters. CONCLUSIONS: Given the heterogeneity of immune infiltration within pedNST, our findings suggest personalized immunogenomic profiling is needed to guide selection of immunotherapeutic strategies.

BACKGROUND: Human bone marrow stromal cells (BMSCs) are an easily accessible and expandable progenitor population with the capacity to generate neural cell types in addition to mesoderm. Lineage tracing studies in transgenic animals have indicated Nestin + BMSCs to be descended from the truncal neural crest. Single-cell analysis provides a means to identify the developmental origin and identity of human BMSC-derived neural progenitors when lineage tracing remains infeasible. This is a prerequisite towards translational application. METHODS: We attained transcriptomic profiles of embryonic long bone, adult human bone marrow, cultured BMSCs and BMSC-derived neurospheres. Integrated scRNAseq analysis was supplemented by characterization of cells during culture expansion and following provision of growth factors and signalling agonists to bias lineage. RESULTS: Reconstructed pseudotime upon the integrated dataset indicated distinct neural and osteogenic differentiation trajectories. The starting state towards the neural differentiation trajectory consisted of Nestin + /MKI67 + BMSCs, which could also be diverted towards the osteogenic trajectory via a branch point. Nestin + /PDGFRA + BMSCs responded to neurosphere culture conditions to generate a subpopulation of cells with a neuronal phenotype according to marker expression and gene ontogeny analysis that occupied the end state along the neural differentiation trajectory. Reconstructed pseudotime also revealed an upregulation of BMP4 expression during culture of BMSC-neurospheres. This provided the rationale for culture supplementation with the BMP signalling agonist SB4, which directed progenitors to upregulate Pax6 and downregulate Nestin. CONCLUSIONS: This study suggested BMSCs originating from truncal neural crest to be the source of cells within long bone marrow possessing neural differentiation potential. Unravelling the transcriptomic dynamics of BMSC-derived neural progenitors promises to enhance differentiation efficiency and safety towards clinical application in cell therapy and disease modelling.

BACKGROUND: The immune response to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection in COVID-19 patients has been extensively investigated. However, much less is known about the long-term effects of infection in patients and how it could affect the immune system and its capacity to respond to future perturbations. METHODS: Using a targeted single-cell multiomics approach, we have recently identified a prolonged anti-inflammatory gene expression signature in T and NK cells in type 1 diabetes patients treated with low-dose IL-2. Here, we investigated the dynamics of this signature in three independent cohorts of COVID-19 patients: (i) the Oxford COVID-19 Multi-omics Blood Atlas (COMBAT) dataset, a cross-sectional cohort including 77 COVID-19 patients and ten healthy donors; (ii) the INCOV dataset, consisting of 525 samples taken from 209 COVID-19 patients during and after infection; and (iii) a longitudinal dataset consisting of 269 whole-blood samples taken from 139 COVID-19 patients followed for a period of up to 7 months after the onset of symptoms using a bulk transcriptomic approach. RESULTS: We discovered that SARS-CoV-2 infection leads to a prolonged alteration of the gene expression profile of circulating T, B and NK cells and monocytes. Some of the genes affected were the same as those present in the IL-2-induced anti-inflammatory gene expression signature but were regulated in the opposite direction, implying a pro-inflammatory status. The altered transcriptional profile was detected in COVID-19 patients for at least 2 months after the onset of the disease symptoms but was not observed in response to influenza infection or sepsis. Gene network analysis suggested a central role for the transcriptional factor NF-κB in the regulation of the observed transcriptional alterations. CONCLUSIONS: SARS-CoV-2 infection causes a prolonged increase in the pro-inflammatory transcriptional status that could predispose post-acute patients to the development of long-term health consequences, including autoimmune disease, reactivation of other viruses and disruption of the host immune system-microbiome ecosystem.

BACKGROUND: Breast cancer patients from the indigenous Arab population present much earlier than patients from Western countries and have traditionally been underrepresented in cancer genomics studies. The contribution of polygenic and Mendelian risk toward the earlier onset of breast cancer in the population remains elusive. METHODS: We performed low-pass whole genome sequencing (lpWGS) and whole-exome sequencing (WES) from 220 female breast cancer patients unselected for positive family history from the indigenous Arab population. Using publicly available resources, we imputed population-specific variants and calculated breast cancer burden-sensitive polygenic risk scores (PRS). Variant pathogenicity was also evaluated on exome variants with high coverage. RESULTS: Variants imputed from lpWGS showed high concordance with paired exome (median dosage correlation: 0.9459, Interquartile range: 0.9410-0.9490). After adjusting the PRS to the Arab population, we found significant associations between PRS performance in risk prediction and first-degree relative breast cancer history prediction (Spearman rho=0.43, p = 0.03), where breast cancer patients in the top PRS decile are 5.53 (95% CI 1.76-17.97, p = 0.003) times more likely also to have a first-degree relative diagnosed with breast cancer compared to those in the middle deciles. In addition, we found evidence for the genetic liability threshold model of breast cancer where among patients with a family history of breast cancer, pathogenic rare variant carriers had significantly lower PRS than non-carriers (p = 0.0205, Mann-Whitney U test) while for non-carriers every standard deviation increase in PRS corresponded to 4.52 years (95% CI 8.88-0.17, p = 0.042) earlier age of presentation. CONCLUSIONS: Overall, our study provides a framework to assess polygenic risk in an understudied population using lpWGS and identifies common variant risk as a factor independent of pathogenic variant carrier status for earlier age of onset of breast cancer among indigenous Arab breast cancer patients.

BACKGROUND: Long-read sequencing (LRS) techniques have been very successful in identifying structural variants (SVs). However, the high error rate of LRS made the detection of small variants (substitutions and short indels < 20 bp) more challenging. The introduction of PacBio HiFi sequencing makes LRS also suited for detecting small variation. Here we evaluate the ability of HiFi reads to detect de novo mutations (DNMs) of all types, which are technically challenging variant types and a major cause of sporadic, severe, early-onset disease. METHODS: We sequenced the genomes of eight parent-child trios using high coverage PacBio HiFi LRS (~ 30-fold coverage) and Illumina short-read sequencing (SRS) (~ 50-fold coverage). De novo substitutions, small indels, short tandem repeats (STRs) and SVs were called in both datasets and compared to each other to assess the accuracy of HiFi LRS. In addition, we determined the parent-of-origin of the small DNMs using phasing. RESULTS: We identified a total of 672 and 859 de novo substitutions/indels, 28 and 126 de novo STRs, and 24 and 1 de novo SVs in LRS and SRS respectively. For the small variants, there was a 92 and 85% concordance between the platforms. For the STRs and SVs, the concordance was 3.6 and 0.8%, and 4 and 100% respectively. We successfully validated 27/54 LRS-unique small variants, of which 11 (41%) were confirmed as true de novo events. For the SRS-unique small variants, we validated 42/133 DNMs and 8 (19%) were confirmed as true de novo event. Validation of 18 LRS-unique de novo STR calls confirmed none of the repeat expansions as true DNM. Confirmation of the 23 LRS-unique SVs was possible for 19 candidate SVs of which 10 (52.6%) were true de novo events. Furthermore, we were able to assign 96% of DNMs to their parental allele with LRS data, as opposed to just 20% with SRS data. CONCLUSIONS: HiFi LRS can now produce the most comprehensive variant dataset obtainable by a single technology in a single laboratory, allowing accurate calling of substitutions, indels, STRs and SVs. The accuracy even allows sensitive calling of DNMs on all variant levels, and also allows for phasing, which helps to distinguish true positive from false positive DNMs.

Epigenetic characterization of cell-free DNA (cfDNA) is an emerging approach for detecting and characterizing diseases such as cancer. We developed a strategy using nanopore-based single-molecule sequencing to measure cfDNA methylomes. This approach generated up to 200 million reads for a single cfDNA sample from cancer patients, an order of magnitude improvement over existing nanopore sequencing methods. We developed a single-molecule classifier to determine whether individual reads originated from a tumor or immune cells. Leveraging methylomes of matched tumors and immune cells, we characterized cfDNA methylomes of cancer patients for longitudinal monitoring during treatment.

BACKGROUND: Array-CGH is the first-tier genetic test both in pre- and postnatal developmental disorders worldwide. Variants of uncertain significance (VUS) represent around 10~15% of reported copy number variants (CNVs). Even though VUS reanalysis has become usual in practice, no long-term study regarding CNV reinterpretation has been reported. METHODS: This retrospective study examined 1641 CGH arrays performed over 8 years (2010-2017) to demonstrate the contribution of periodically re-analyzing CNVs of uncertain significance. CNVs were classified using AnnotSV on the one hand and manually curated on the other hand. The classification was based on the 2020 American College of Medical Genetics (ACMG) criteria. RESULTS: Of the 1641 array-CGH analyzed, 259 (15.7%) showed at least one CNV initially reported as of uncertain significance. After reinterpretation, 106 of the 259 patients (40.9%) changed categories, and 12 of 259 (4.6%) had a VUS reclassified to likely pathogenic or pathogenic. Six were predisposing factors for neurodevelopmental disorder/autism spectrum disorder (ASD). CNV type (gain or loss) does not seem to impact the reclassification rate, unlike the length of the CNV: 75% of CNVs downgraded to benign or likely benign are less than 500 kb in size. CONCLUSIONS: This study's high rate of reinterpretation suggests that CNV interpretation has rapidly evolved since 2010, thanks to the continuous enrichment of available databases. The reinterpreted CNV explained the phenotype for ten patients, leading to optimal genetic counseling. These findings suggest that CNVs should be reinterpreted at least every 2 years.

BACKGROUND: The association between microbes and cancer has been reported repeatedly; however, it is not clear if molecular tumour properties are connected to specific microbial colonisation patterns. This is due mainly to the current technical and analytical strategy limitations to characterise tumour-associated bacteria. METHODS: Here, we propose an approach to detect bacterial signals in human RNA sequencing data and associate them with the clinical and molecular properties of the tumours. The method was tested on public datasets from The Cancer Genome Atlas, and its accuracy was assessed on a new cohort of colorectal cancer patients. RESULTS: Our analysis shows that intratumoural microbiome composition is correlated with survival, anatomic location, microsatellite instability, consensus molecular subtype and immune cell infiltration in colon tumours. In particular, we find Faecalibacterium prausnitzii, Coprococcus comes, Bacteroides spp., Fusobacterium spp. and Clostridium spp. to be strongly associated with tumour properties. CONCLUSIONS: We implemented an approach to concurrently analyse clinical and molecular properties of the tumour as well as the composition of the associated microbiome. Our results may improve patient stratification and pave the path for mechanistic studies on microbiota-tumour crosstalk.

BACKGROUND: Fumarate hydratase-deficient renal cell carcinoma (FH-RCC) is a rare highly aggressive subtype of kidney cancer for which the distinct genomic, transcriptomic, and evolutionary relationships between metastatic and primary lesions are still unclear. METHODS: In this study, whole-exome, RNA-seq, and DNA methylation sequencing were performed on primary-metastatic paired specimens from 19 FH-RCC cases, including 23 primary and 35 matched metastatic lesions. Phylogenetic and clonal evolutionary analyses were used to investigate the evolutionary characteristics of FH-RCC. Transcriptomic analyses, immunohistochemistry, and multiple immunofluorescence experiments were performed to identify the tumor microenvironmental features of metastatic lesions. RESULTS: Paired primary and metastatic lesions generally showed similar characteristics of tumor mutation burden, tumor neoantigen burden, microsatellite instability score, CNV burden, and genome instability index. Notably, we identified an FH-mutated founding MRCA (the most recent common ancestor) clone that dominated the early evolutionary trajectories in FH-RCC. Although both primary and metastatic lesions manifested high immunogenicity, metastatic lesions exhibited higher enrichment of T effector cells and immune-related chemokines, together with upregulation of PD-L1, TIGIT, and BTLA. In addition, we found that concurrent NF2 mutation may be associated with bone metastasis and upregulation of cell cycle signature in metastatic lesions. Furthermore, although in FH-RCC metastatic lesions in general shared similar CpG island methylator phenotype with primary lesions, we found metastatic lesions displaying hypomethylated chemokine and immune checkpoints related genomic loci. CONCLUSIONS: Overall, our study demonstrated the genomic, epigenomic, and transcriptomic features of metastatic lesions in FH-RCC and revealed their early evolutionary trajectory. These results provided multi-omics evidence portraying the progression of FH-RCC.

BACKGROUND: Transcriptional classification has been used to stratify colorectal cancer (CRC) into molecular subtypes with distinct biological and clinical features. However, it is not clear whether such subtypes represent discrete, mutually exclusive entities or molecular/phenotypic states with potential overlap. Therefore, we focused on the CRC Intrinsic Subtype (CRIS) classifier and evaluated whether assigning multiple CRIS subtypes to the same sample provides additional clinically and biologically relevant information. METHODS: A multi-label version of the CRIS classifier (multiCRIS) was applied to newly generated RNA-seq profiles from 606 CRC patient-derived xenografts (PDXs), together with human CRC bulk and single-cell RNA-seq datasets. Biological and clinical associations of single- and multi-label CRIS were compared. Finally, a machine learning-based multi-label CRIS predictor (ML2CRIS) was developed for single-sample classification. RESULTS: Surprisingly, about half of the CRC cases could be significantly assigned to more than one CRIS subtype. Single-cell RNA-seq analysis revealed that multiple CRIS membership can be a consequence of the concomitant presence of cells of different CRIS class or, less frequently, of cells with hybrid phenotype. Multi-label assignments were found to improve prediction of CRC prognosis and response to treatment. Finally, the ML2CRIS classifier was validated for retaining the same biological and clinical associations also in the context of single-sample classification. CONCLUSIONS: These results show that CRIS subtypes retain their biological and clinical features even when concomitantly assigned to the same CRC sample. This approach could be potentially extended to other cancer types and classification systems.

BACKGROUND: Hearing problems (HP) in adults are common and are associated with several comorbid conditions. Its prevalence increases with age, reflecting the cumulative effect of environmental factors and genetic predisposition. Although several risk loci have been already identified, HP biology and epidemiology are still insufficiently investigated by large-scale genetic studies. METHODS: Leveraging the UK Biobank, the Nurses' Health Studies (I and II), the Health Professionals Follow-up Study, and the Million Veteran Program, we conducted a comprehensive genome-wide investigation of HP in 748,668 adult participants (discovery N = 501,825; replication N = 226,043; cross-ancestry replication N = 20,800). We leveraged the GWAS findings to characterize HP polygenic architecture, exploring sex differences, polygenic risk across ancestries, tissue-specific transcriptomic regulation, cause-effect relationships with genetically correlated traits, and gene interactions with HP environmental risk factors. RESULTS: We identified 54 risk loci and demonstrated that HP polygenic risk is shared across ancestry groups. Our transcriptomic regulation analysis highlighted the potential role of the central nervous system in HP pathogenesis. The sex-stratified analyses showed several additional associations related to peripheral hormonally regulated tissues reflecting a potential role of estrogen in hearing function. This evidence was supported by the multivariate interaction analysis that showed how genes involved in brain development interact with sex, noise pollution, and tobacco smoking in relation to their HP associations. Additionally, the genetically informed causal inference analysis showed that HP is linked to many physical and mental health outcomes. CONCLUSIONS: The results provide many novel insights into the biology and epidemiology of HP in adults. Our sex-specific analyses and transcriptomic associations highlighted molecular pathways that may be targeted for drug development or repurposing. Additionally, the potential causal relationships identified may support novel preventive screening programs to identify individuals at risk.

BACKGROUND: High sequence identity between segmental duplications (SDs) can facilitate copy number variants (CNVs) via non-allelic homologous recombination (NAHR). These CNVs are one of the fundamental causes of genomic disorders such as the 3q29 deletion syndrome (del3q29S). There are 21 protein-coding genes lost or gained as a result of such recurrent 1.6-Mbp deletions or duplications, respectively, in the 3q29 locus. While NAHR plays a role in CNV occurrence, the factors that increase the risk of NAHR at this particular locus are not well understood. METHODS: We employed an optical genome mapping technique to characterize the 3q29 locus in 161 unaffected individuals, 16 probands with del3q29S and their parents, and 2 probands with the 3q29 duplication syndrome (dup3q29S). Long-read sequencing-based haplotype resolved de novo assemblies from 44 unaffected individuals, and 1 trio was used for orthogonal validation of haplotypes and deletion breakpoints. RESULTS: In total, we discovered 34 haplotypes, of which 19 were novel haplotypes. Among these 19 novel haplotypes, 18 were detected in unaffected individuals, while 1 novel haplotype was detected on the parent-of-origin chromosome of a proband with the del3q29S. Phased assemblies from 44 unaffected individuals enabled the orthogonal validation of 20 haplotypes. In 89% (16/18) of the probands, breakpoints were confined to paralogous copies of a 20-kbp segment within the 3q29 SDs. In one del3q29S proband, the breakpoint was confined to a 374-bp region using long-read sequencing. Furthermore, we categorized del3q29S cases into three classes and dup3q29S cases into two classes based on breakpoints. Finally, we found no evidence of inversions in parent-of-origin chromosomes. CONCLUSIONS: We have generated the most comprehensive haplotype map for the 3q29 locus using unaffected individuals, probands with del3q29S or dup3q29S, and available parents, and also determined the deletion breakpoint to be within a 374-bp region in one proband with del3q29S. These results should provide a better understanding of the underlying genetic architecture that contributes to the etiology of del3q29S and dup3q29S.

BACKGROUND: Single-cell RNA sequencing (scRNA-seq) provides valuable insights into human islet cell types and their corresponding stable gene expression profiles. However, this approach requires cell dissociation that complicates its utility in vivo. On the other hand, single-nucleus RNA sequencing (snRNA-seq) has compatibility with frozen samples, elimination of dissociation-induced transcriptional stress responses, and affords enhanced information from intronic sequences that can be leveraged to identify pre-mRNA transcripts. METHODS: We obtained nuclear preparations from fresh human islet cells and generated snRNA-seq datasets. We compared these datasets to scRNA-seq output obtained from human islet cells from the same donor. We employed snRNA-seq to obtain the transcriptomic profile of human islets engrafted in immunodeficient mice. In both analyses, we included the intronic reads in the snRNA-seq data with the GRCh38-2020-A library. RESULTS: First, snRNA-seq analysis shows that the top four differentially and selectively expressed genes in human islet endocrine cells in vitro and in vivo are not the canonical genes but a new set of non-canonical gene markers including ZNF385D, TRPM3, LRFN2, PLUT (β-cells); PTPRT, FAP, PDK4, LOXL4 (α-cells); LRFN5, ADARB2, ERBB4, KCNT2 (δ-cells); and CACNA2D3, THSD7A, CNTNAP5, RBFOX3 (γ-cells). Second, by integrating information from scRNA-seq and snRNA-seq of human islet cells, we distinguish three β-cell sub-clusters: an INS pre-mRNA cluster (β3), an intermediate INS mRNA cluster (β2), and an INS mRNA-rich cluster (β1). These display distinct gene expression patterns representing different biological dynamic states both in vitro and in vivo. Interestingly, the INS mRNA-rich cluster (β1) becomes the predominant sub-cluster in vivo. CONCLUSIONS: In summary, snRNA-seq and pre-mRNA analysis of human islet cells can accurately identify human islet cell populations, subpopulations, and their dynamic transcriptome profile in vivo.

BACKGROUND: Medulloblastoma (MB) is a malignant tumour of the cerebellum which can be classified into four major subgroups based on gene expression and genomic features. Single-cell transcriptome studies have defined the cellular states underlying each MB subgroup; however, the spatial organisation of these diverse cell states and how this impacts response to therapy remains to be determined. METHODS: Here, we used spatially resolved transcriptomics to define the cellular diversity within a sonic hedgehog (SHH) patient-derived model of MB and show that cells specific to a transcriptional state or spatial location are pivotal for CDK4/6 inhibitor, Palbociclib, treatment response. We integrated spatial gene expression with histological annotation and single-cell gene expression data from MB, developing an analysis strategy to spatially map cell type responses within the hybrid system of human and mouse cells and their interface within an intact brain tumour section. RESULTS: We distinguish neoplastic and non-neoplastic cells within tumours and from the surrounding cerebellar tissue, further refining pathological annotation. We identify a regional response to Palbociclib, with reduced proliferation and induced neuronal differentiation in both treated tumours. Additionally, we resolve at a cellular resolution a distinct tumour interface where the tumour contacts neighbouring mouse brain tissue consisting of abundant astrocytes and microglia and continues to proliferate despite Palbociclib treatment. CONCLUSIONS: Our data highlight the power of using spatial transcriptomics to characterise the response of a tumour to a targeted therapy and provide further insights into the molecular and cellular basis underlying the response and resistance to CDK4/6 inhibitors in SHH MB.

BACKGROUND: There has been a growing appreciation recently that mutagenic processes can be studied through the lenses of mutational signatures, which represent characteristic mutation patterns attributed to individual mutagens. However, the causal links between mutagens and observed mutation patterns as well as other types of interactions between mutagenic processes and molecular pathways are not fully understood, limiting the utility of mutational signatures. METHODS: To gain insights into these relationships, we developed a network-based method, named GENESIGNET that constructs an influence network among genes and mutational signatures. The approach leverages sparse partial correlation among other statistical techniques to uncover dominant influence relations between the activities of network nodes. RESULTS: Applying GENESIGNET to cancer data sets, we uncovered important relations between mutational signatures and several cellular processes that can shed light on cancer-related processes. Our results are consistent with previous findings, such as the impact of homologous recombination deficiency on clustered APOBEC mutations in breast cancer. The network identified by GENESIGNET also suggest an interaction between APOBEC hypermutation and activation of regulatory T Cells (Tregs), as well as a relation between APOBEC mutations and changes in DNA conformation. GENESIGNET also exposed a possible link between the SBS8 signature of unknown etiology and the Nucleotide Excision Repair (NER) pathway. CONCLUSIONS: GENESIGNET provides a new and powerful method to reveal the relation between mutational signatures and gene expression. The GENESIGNET method was implemented in python, and installable package, source codes and the data sets used for and generated during this study are available at the Github site https://github.com/ncbi/GeneSigNet.

BACKGROUND: The incidence of early-onset colorectal cancer (EOCRC; patients < 50 years old) has been rising rapidly, whereas the EOCRC genetic susceptibility remains incompletely investigated. Here, we aimed to systematically identify specific susceptible genetic variants for EOCRC. METHODS: Two parallel GWASs were conducted in 17,789 CRC cases (including 1490 EOCRC cases) and 19,951 healthy controls. A polygenic risk score (PRS) model was built based on identified EOCRC-specific susceptibility variants by using the UK Biobank cohort. We also interpreted the potential biological mechanisms of the prioritized risk variant. RESULTS: We identified 49 independent susceptibility loci that were significantly associated with the susceptibility to EOCRC and the diagnosed age of CRC (both P < 5.0×10-4), replicating 3 previous CRC GWAS loci. There are 88 assigned susceptibility genes involved in chromatin assembly and DNA replication pathways, mainly associating with precancerous polyps. Additionally, we assessed the genetic effect of the identified variants by developing a PRS model. Compared to the individuals in the low genetic risk group, the individuals in the high genetic risk group have increased EOCRC risk, and these results were replicated in the UKB cohort with a 1.63-fold risk (95% CI: 1.32-2.02, P = 7.67×10-6). The addition of the identified EOCRC risk loci significantly increased the prediction accuracy of the PRS model, compared to the PRS model derived from the previous GWAS-identified loci. Mechanistically, we also elucidated that rs12794623 may contribute to the early stage of CRC carcinogenesis via allele-specific regulating the expression of POLA2. CONCLUSIONS: These findings will broaden the understanding of the etiology of EOCRC and may facilitate the early screening and individualized prevention.

BACKGROUND: Sex chromosome aneuploidies (SCAs) give rise to a broad range of phenotypic traits and diseases. Previous studies based on peripheral blood samples have suggested the presence of ripple effects, caused by altered X chromosome number, affecting the methylome and transcriptome. Whether these alterations can be connected to disease-specific tissues, and thereby having clinical implication for the phenotype, remains to be elucidated. METHODS: We performed a comprehensive analysis of X chromosome number on the transcriptome and methylome in blood, fat, and muscle tissue from individuals with 45,X, 46,XX, 46,XY, and 47,XXY. RESULTS: X chromosome number affected the transcriptome and methylome globally across all chromosomes in a tissue-specific manner. Furthermore, 45,X and 47,XXY demonstrated a divergent pattern of gene expression and methylation, with overall gene downregulation and hypomethylation in 45,X and gene upregulation and hypermethylation in 47,XXY. In fat and muscle, a pronounced effect of sex was observed. We identified X chromosomal genes with an expression pattern different from what would be expected based on the number of X and Y chromosomes. Our data also indicate a regulatory function of Y chromosomal genes on X chromosomal genes. Fourteen X chromosomal genes were downregulated in 45,X and upregulated in 47,XXY, respectively, in all three tissues (AKAP17A, CD99, DHRSX, EIF2S3, GTPBP6, JPX, KDM6A, PP2R3B, PUDP, SLC25A6, TSIX, XIST, ZBED1, ZFX). These genes may be central in the epigenetic and genomic regulation of sex chromosome aneuploidies. CONCLUSION: We highlight a tissue-specific and complex effect of X chromosome number on the transcriptome and methylome, elucidating both shared and non-shared gene-regulatory mechanism between SCAs.

Since many single-cell RNA-seq (scRNA-seq) data are obtained after cell sorting, such as when investigating immune cells, tracking cellular landscape by integrating single-cell data with spatial transcriptomic data is limited due to cell type and cell composition mismatch between the two datasets. We developed a method, spSeudoMap, which utilizes sorted scRNA-seq data to create virtual cell mixtures that closely mimic the gene expression of spatial data and trains a domain adaptation model for predicting spatial cell compositions. The method was applied in brain and breast cancer tissues and accurately predicted the topography of cell subpopulations. spSeudoMap may help clarify the roles of a few, but crucial cell types.

BACKGROUND: Rapidly and efficiently identifying critically ill infants for whole genome sequencing (WGS) is a costly and challenging task currently performed by scarce, highly trained experts and is a major bottleneck for application of WGS in the NICU. There is a dire need for automated means to prioritize patients for WGS. METHODS: Institutional databases of electronic health records (EHRs) are logical starting points for identifying patients with undiagnosed Mendelian diseases. We have developed automated means to prioritize patients for rapid and whole genome sequencing (rWGS and WGS) directly from clinical notes. Our approach combines a clinical natural language processing (CNLP) workflow with a machine learning-based prioritization tool named Mendelian Phenotype Search Engine (MPSE). RESULTS: MPSE accurately and robustly identified NICU patients selected for WGS by clinical experts from Rady Children's Hospital in San Diego (AUC 0.86) and the University of Utah (AUC 0.85). In addition to effectively identifying patients for WGS, MPSE scores also strongly prioritize diagnostic cases over non-diagnostic cases, with projected diagnostic yields exceeding 50% throughout the first and second quartiles of score-ranked patients. CONCLUSIONS: Our results indicate that an automated pipeline for selecting acutely ill infants in neonatal intensive care units (NICU) for WGS can meet or exceed diagnostic yields obtained through current selection procedures, which require time-consuming manual review of clinical notes and histories by specialized personnel.

BACKGROUND: Next-generation sequencing (NGS) based population screening holds great promise for disease prevention and earlier diagnosis, but the costs associated with screening millions of humans remain prohibitive. New methods for population genetic testing that lower the costs of NGS without compromising diagnostic power are needed. METHODS: We developed double batched sequencing where DNA samples are batch-sequenced twice - directly pinpointing individuals with rare variants. We sequenced batches of at-birth blood spot DNA using a commercial 113-gene panel in an explorative (n = 100) and a validation (n = 100) cohort of children who went on to develop pediatric cancers. All results were benchmarked against individual whole genome sequencing data. RESULTS: We demonstrated fully replicable detection of cancer-causing germline variants, with positive and negative predictive values of 100% (95% CI, 0.91-1.00 and 95% CI, 0.98-1.00, respectively). Pathogenic and clinically actionable variants were detected in RB1, TP53, BRCA2, APC, and 19 other genes. Analyses of larger batches indicated that our approach is highly scalable, yielding more than 95% cost reduction or less than 3 cents per gene screened for rare disease-causing mutations. We also show that double batched sequencing could cost-effectively prevent childhood cancer deaths through broad genomic testing. CONCLUSIONS: Our ultracheap genetic diagnostic method, which uses existing sequencing hardware and standard newborn blood spots, should readily open up opportunities for population-wide risk stratification using genetic screening across many fields of clinical genetics and genomics.