The UCSC Genome Browser database: extensions and updates 2013

Site-specific DNA demethylation during spermatogenesis presets the sites of nucleosome retention in mouse sperm

DNA methylation patterns are inherited from the parental germline to the embryo. In mature sperm, the sites of unmethylated DNA are tightly coupled to sites of histone retention at gene regulatory elements that are implicated in paternal epigenetic inheritance. The timing and mechanism of site-specific DNA demethylation in the male germline currently remains unknown. Here, we perform genome-wide profiling of DNA methylation during spermatogenesis by capturing methylated DNA through interaction with a methyl-DNA binding protein domain (MBD). Our data demonstrate that there is a site-specific change in DNA methylation during the mitosis-to-meiosis transition. Importantly, the genomic sites that are demethylated during this transition predetermine nucleosome retention sites in spermatozoa. These results suggest that site-specific DNA demethylation during the mitosis-to-meiosis transition of spermatogenesis prepares embryonic gene expression after fertilization. We therefore propose DNA demethylation during spermatogenesis as a novel phase of epigenetic reprogramming that contributes to embryonic gene regulation.

GGNBP2 regulates MDA5 sensing triggered by self double-stranded RNA following loss of ADAR1 editing

Adenosine-to-inosine (A-to-I) editing of double-stranded RNA (dsRNA) by ADAR1 is an essential modifier of the immunogenicity of cellular dsRNA. The role of MDA5 in sensing unedited cellular dsRNA and the downstream activation of type I interferon (IFN) signaling are well established. However, we have an incomplete understanding of pathways that modify the response to unedited dsRNA. We performed a genome-wide CRISPR screen and showed that GGNBP2, CNOT10, and CNOT11 interact and regulate sensing of unedited cellular dsRNA. We found that GGNBP2 acts between dsRNA transcription and its cytoplasmic sensing by MDA5. GGNBP2 loss prevented induction of type I IFN and autoinflammation after the loss of ADAR1 editing activity by modifying the subcellular distribution of endogenous A-to-I editing substrates and reducing cytoplasmic dsRNA load. These findings reveal previously undescribed pathways to modify diseases associated with ADAR mutations and may be determinants of response or resistance to small-molecule ADAR1 inhibitors.

CRB2 Depletion Induces YAP Signaling and Disrupts Mechanosensing in Podocytes

Focal Segmental Glomerulosclerosis (FSGS) is a histologic lesion caused by a variety of injurious stimuli that lead to dysfunction/loss of glomerular visceral epithelial cells (i.e. podocytes). Pathogenic mutations in CRB2, encoding the type 1 transmembrane protein Crumb 2 Homolog Protein, have been shown to cause early-onset corticosteroid-resistant nephrotic syndrome (SRNS)/FSGS. Here, we identified a 2-generation East Asian kindred (DUK40595) with biopsy-proven SRNS/FSGS caused by a compound heterozygous mutation in CRB2 comprised of the previously described truncating mutation p.Gly1036_Alafs*43 and a rare 9-bp deletion mutation p.Leu1074_Asp1076del. Because compound heterozygous mutations involving the truncating p.Gly1036_Alafs*43 variant have been associated with reduced CRB2 expression in podocytes and autosomal recessive SRNS/FSGS, we sought to define the pathogenic effects of CRB2 deficiency in podocytes. We show that CRB2 knockdown induces YAP activity and target gene expression in podocytes. It upregulates YAP-mediated mechanosignaling and increases the density of focal adhesion and F-actin. Using Elastic Resonator Interference Stress Microscopy (ERISM), we demonstrate that CRB2 knockdown also enhances podocyte contractility in a substrate stiffness-dependent manner. The knockdown effect decreases with increasing substrate stiffness, indicating impaired mechanosensing in CRB2 knockdown cells at low substrate stiffness. While the mechanical activation of CRB2 knockdown cells is associated with increased YAP activity, the enhanced cell contractility is not significantly reduced by the selective YAP inhibitors K-975 and verteporfin, suggesting that multiple pathways may be involved in mechanosignaling downstream of CRB2. Taken together, these studies provide the first evidence that CRB2 deficiency may impair podocyte mechanotransduction via disruption of YAP signaling in podocytes.

APC mutations dysregulate alternative polyadenylation in cancer

BACKGROUND

Alternative polyadenylation (APA) affects most human genes and is recurrently dysregulated in all studied cancers. However, the mechanistic origins of this dysregulation are incompletely understood.

RESULTS

We describe an unbiased analysis of molecular regulators of poly(A) site selection across The Cancer Genome Atlas and identify that colorectal adenocarcinoma is an outlier relative to all other cancer subtypes. This distinction arises from the frequent presence of loss-of-function APC mutations in colorectal adenocarcinoma, which are strongly associated with long 3' UTR expression relative to tumors lacking APC mutations. APC knockout similarly dysregulates APA in human colon organoids. By mining previously published APC eCLIP data, we show that APC preferentially binds G- and C-rich motifs just upstream of proximal poly(A) sites. Lastly, we find that reduced APC expression is associated with APA dysregulation in tumor types lacking recurrent APC mutations.

CONCLUSIONS

As APC has been previously identified as an RNA-binding protein that preferentially binds 3' UTRs during mouse neurogenesis, our results suggest that APC promotes proximal poly(A) site use and that APC loss and altered expression contribute to pervasive APA dysregulation in cancers.

An autoregulatory poison exon in Smndc1 is conserved across kingdoms and influences organism growth

Many of the most highly conserved elements in the human genome are "poison exons," alternatively spliced exons that contain premature termination codons and permit post-transcriptional regulation of mRNA abundance through induction of nonsense-mediated mRNA decay (NMD). Poison exons are widely assumed to be highly conserved due to their presumed importance for organismal fitness, but this functional importance has never been tested in the context of a whole organism. Here, we report that a poison exon in Smndc1 is conserved across mammals and plants and plays a molecular autoregulatory function in both kingdoms. We generated mouse and A. thaliana models lacking this poison exon to find its loss leads to deregulation of SMNDC1 protein levels, pervasive alterations in mRNA processing, and organismal size restriction. Together, these models demonstrate the importance of poison exons for both molecular and organismal phenotypes that likely explain their extraordinary conservation.

RBM10 loss induces aberrant splicing of cytoskeletal and extracellular matrix mRNAs and promotes metastatic fitness

RBM10 modulates transcriptome-wide cassette exon splicing. Loss-of-function RBM10 mutations are enriched in thyroid cancers with distant metastases. Analysis of transcriptomes and genes mis-spliced by RBM10 loss showed pro-migratory and RHO/RAC signaling signatures. RBM10 loss increases cell velocity. Cytoskeletal and ECM transcripts subject to exon-inclusion events included vinculin (VCL), tenascin C (TNC) and CD44. Knockdown of the VCL exon inclusion transcript in RBM10-null cells reduced cell velocity, whereas knockdown of TNC and CD44 exon-inclusion isoforms reduced invasiveness. RAC1-GTP levels were increased in RBM10-null cells. Mouse Hras G12V /Rbm1O KO thyrocytes develop metastases that are reversed by RBM10 or by combined knockdown of VCL, CD44 and TNC inclusion isoforms. Thus, RBM10 loss generates exon inclusions in transcripts regulating ECM-cytoskeletal interactions, leading to RAC1 activation and metastatic competency. Moreover, a CRISPR-Cas9 screen for synthetic lethality with RBM10 loss identified NFkB effectors as central to viability, providing a therapeutic target for these lethal thyroid cancers.

DUX4 is a common driver of immune evasion and immunotherapy failure in metastatic cancers

Cancer immune evasion contributes to checkpoint immunotherapy failure in many patients with metastatic cancers. The embryonic transcription factor DUX4 was recently characterized as a suppressor of interferon-γ signaling and antigen presentation that is aberrantly expressed in a small subset of primary tumors. Here, we report that DUX4 expression is a common feature of metastatic tumors, with ~10-50% of advanced bladder, breast, kidney, prostate, and skin cancers expressing DUX4. DUX4 expression is significantly associated with immune cell exclusion and decreased objective response to PD-L1 blockade in a large cohort of urothelial carcinoma patients. DUX4 expression is a significant predictor of survival even after accounting for tumor mutational burden and other molecular and clinical features in this cohort, with DUX4 expression associated with a median reduction in survival of over 1 year. Our data motivate future attempts to develop DUX4 as a biomarker and therapeutic target for checkpoint immunotherapy resistance.

Sex-biased gene regulation varies across human populations as a result of adaptive evolution

OBJECTIVES

Studies of human sexual dimorphism and gender disparities in health focus on ostensibly universal molecular sex differences, such as sex chromosomes and circulating hormone levels, while ignoring the extraordinary diversity in biology, behavior, and culture acquired by different human populations over their unique evolutionary histories.

MATERIALS AND METHODS

Using RNA-Seq data and whole genome sequences from 1000G and HGDP, we investigate variation in sex-biased gene expression across 11 human populations and test whether population-level variation in sex-biased expression may have resulted from adaptive evolution in regions containing sex-specific regulatory variants.

RESULTS

We find that sex-biased gene expression in humans is highly variable, mostly population-specific, and demonstrates between population reversals. Expression quantitative trait locus mapping reveals sex-specific regulatory regions with evidence of recent positive natural selection, suggesting that variation in sex-biased expression may have evolved as an adaptive response to ancestral environments experienced by human populations.

DISCUSSION

These results indicate that sex-biased gene expression is more flexible than previously thought and is not generally shared among human populations. Instead, molecular phenotypes associated with sex depend on complex interactions between population-specific molecular evolution and physiological responses to contemporary socioecologies.

Multiplexed screening reveals how cancer-specific alternative polyadenylation shapes tumor growth in vivo

Alternative polyadenylation (APA) is strikingly dysregulated in many cancers. Although global APA dysregulation is frequently associated with poor prognosis, the importance of most individual APA events is controversial simply because few have been functionally studied. Here, we address this gap by developing a CRISPR-Cas9-based screen to manipulate endogenous polyadenylation and systematically quantify how APA events contribute to tumor growth in vivo. Our screen reveals individual APA events that control mouse melanoma growth in an immunocompetent host, with concordant associations in clinical human cancer. For example, forced Atg7 3' UTR lengthening in mouse melanoma suppresses ATG7 protein levels, slows tumor growth, and improves host survival; similarly, in clinical human melanoma, a long ATG7 3' UTR is associated with significantly prolonged patient survival. Overall, our study provides an easily adaptable means to functionally dissect APA in physiological systems and directly quantifies the contributions of recurrent APA events to tumorigenic phenotypes.

scTEA-db: a comprehensive database of novel terminal exon isoforms identified from human single cell transcriptomes

The usage of alternative terminal exons results in messenger RNA (mRNA) isoforms that differ in their 3' untranslated regions (3' UTRs) and often also in their protein-coding sequences. Alternative 3' UTRs contain different sets of cis-regulatory elements known to regulate mRNA stability, translation and localization, all of which are vital to cell identity and function. In previous work, we revealed that ∼25 percent of the experimentally observed RNA 3' ends are located within regions currently annotated as intronic, indicating that many 3' end isoforms remain to be uncovered. Also, the inclusion of not yet annotated terminal exons is more tissue specific compared to the already annotated ones. Here, we present the single cell-based Terminal Exon Annotation database (scTEA-db, www.scTEA-db.org) that provides the community with 12 063 so far not yet annotated terminal exons and associated transcript isoforms identified by analysing 53 069 publicly available single cell transcriptomes. Our scTEA-db web portal offers an array of features to find and explore novel terminal exons belonging to 5538 human genes, 110 of which are known cancer drivers. In summary, scTEA-db provides the foundation for studying the biological role of large numbers of so far not annotated terminal exon isoforms in cell identity and function.

Multigene testing panels reveal pathogenic variants in sporadic breast cancer patients in northern China

Background: Breast cancer, the most prevalent malignancy in women worldwide, presents diverse onset patterns and genetic backgrounds. This study aims to examine the genetic landscape and clinical implications of rare mutations in Chinese breast cancer patients. Methods: Clinical data from 253 patients, including sporadic and familial cases, were analyzed. Comprehensive genomic profiling was performed, categorizing identified rare variants according to the American College of Medical Genetics (ACMG) guidelines. In silico protein modeling was used to analyze potentially pathogenic variants' impact on protein structure and function. Results: We detected 421 rare variants across patients. The most frequently mutated genes were ALK (22.2%), BARD1 (15.6%), and BRCA2 (15.0%). ACMG classification identified 7% of patients harboring Pathogenic/Likely Pathogenic (P/LP) variants, with one case displaying a pathogenic BRCA1 mutation linked to triple-negative breast cancer (TNBC). Also identified were two pathogenic MUTYH variants, previously associated with colon cancer but increasingly implicated in breast cancer. Variants of uncertain significance (VUS) were identified in 112 patients, with PTEN c.C804A showing the highest frequency. The role of these variants in sporadic breast cancer oncogenesis was suggested. In-depth exploration of previously unreported variants led to the identification of three potential pathogenic variants: ATM c.C8573T, MSH3 c.A2723T, and CDKN1C c.C221T. Their predicted impact on protein structure and stability suggests a functional role in cancer development. Conclusion: This study reveals a comprehensive overview of the genetic variants landscape in Chinese breast cancer patients, highlighting the prevalence and potential implications of rare variants. We emphasize the value of comprehensive genomic profiling in breast cancer management and the necessity of continuous research into understanding the functional impacts of these variants.

Genome-wide methylation patterns from canine nanopore assemblies

Recent advances in long-read sequencing have enabled the creation of reference-quality genome assemblies for multiple individuals within a species. In particular, 8 long-read genome assemblies have recently been published for the canine model (dogs and wolves). These assemblies were created using a range of sequencing and computational approaches, with only limited comparisons described among subsets of the assemblies. Here we present 3 high-quality de novo reference assemblies based upon Oxford Nanopore long-read sequencing: 2 Bernese Mountain Dogs (BD & OD) and a Cairn terrier (CA611). These breeds are of particular interest due to the enrichment of unresolved genetic disorders. Leveraging advancement in software technologies, we utilized published data of Labrador Retriever (Yella) to generate a new assembly, resulting in a ∼280-fold increase in continuity (N50 size of 91 kbp vs 25.75 Mbp). In conjunction with these 4 new assemblies, we uniformly assessed 8 existing assemblies for generalized quality metrics, sequence divergence, and a detailed BUSCO assessment. We identified a set of ∼400 conserved genes during the BUSCO analysis missing in all assemblies. Genome-wide methylation profiles were generated from the nanopore sequencing, resulting in broad concordance with existing whole-genome and reduced-representation bisulfite sequencing, while highlighting superior overage of mobile elements. These analyses demonstrate the ability of Nanopore sequencing to resolve the sequence and epigenetic profile of canine genomes.

Establishment and characterization of CSCRi006-A: an induced pluripotent stem cell line generated from a patient with Diamond-Blackfan Anemia (DBA) carrying ribosomal protein S19 (RPS19) mutation

Diamond-Blackfan anemia (DBA) is a congenital hypoplastic anemia characterized by ineffective erythropoiesis. DBA is majorly caused by mutations in the ribosomal protein (RP) genes (Gadhiya and Wills in Diamond-Blackfan Anemia, https://www.statpearls.com/ ; 2023). A suitable disease model that yields a continuous supply of erythroid cells is required to study disease pathogenesis and drug discovery. Toward this, we reprogrammed dermal fibroblasts from a DBA patient with a heterozygous mutation c.22-23delAG in the RPS19 gene identified through exome sequencing. To generate induced pluripotent stem cells (iPSCs), we induced episomal expression of the reprogramming factors OTC3/4, L-MYC, LIN28, SOX2, and KLF4, and a p53 shRNA2. The DBA-iPSC line CSCRi006-A generated during this study was extensively characterized for its pluripotency and genome stability. The clone retained normal karyotype and showed high expression levels of pluripotency markers, OCT4, NANOG, SOX2, TRA-I-60, TRA-I-81, and SSEA4. It could differentiate into cells originating from all three germ cell layers, as identified by immunostaining for SOX17 (endoderm), Brachyury (mesoderm), and PAX6 (ectoderm). IPSCs provide a renewable source of cells for in vitro disease modeling. CSCRi006-A, a thoroughly characterized iPSC line carrying heterozygous RPS19 c.22-23delAG mutation, is a valuable cell line for the disease modeling of DBA. This iPSC line can be differentiated into different blood cell types to study the mechanisms of disease development and identify potential treatments.

Nutrition impact on ILC3 maintenance and function centers on a cell-intrinsic CD71-iron axis

Iron metabolism is pivotal for cell fitness in the mammalian host; however, its role in group 3 innate lymphoid cells (ILC3s) is unknown. Here we show that transferrin receptor CD71 (encoded by Tfrc)-mediated iron metabolism cell-intrinsically controls ILC3 proliferation and host protection against Citrobacter rodentium infection and metabolically affects mitochondrial respiration by switching of oxidative phosphorylation toward glycolysis. Iron deprivation or Tfrc ablation in ILC3s reduces the expression and/or activity of the aryl hydrocarbon receptor (Ahr), a key ILC3 regulator. Genetic ablation or activation of Ahr in ILC3s leads to CD71 upregulation or downregulation, respectively, suggesting Ahr-mediated suppression of CD71. Mechanistically, Ahr directly binds to the Tfrc promoter to inhibit transcription. Iron overload partially restores the defective ILC3 compartment in the small intestine of Ahr-deficient mice, consistent with the compensatory upregulation of CD71. These data collectively demonstrate an under-appreciated role of the Ahr-CD71-iron axis in the regulation of ILC3 maintenance and function.

OpenPBTA: The Open Pediatric Brain Tumor Atlas

Pediatric brain and spinal cancers are collectively the leading disease-related cause of death in children; thus, we urgently need curative therapeutic strategies for these tumors. To accelerate such discoveries, the Children's Brain Tumor Network (CBTN) and Pacific Pediatric Neuro-Oncology Consortium (PNOC) created a systematic process for tumor biobanking, model generation, and sequencing with immediate access to harmonized data. We leverage these data to establish OpenPBTA, an open collaborative project with over 40 scalable analysis modules that genomically characterize 1,074 pediatric brain tumors. Transcriptomic classification reveals universal TP53 dysregulation in mismatch repair-deficient hypermutant high-grade gliomas and TP53 loss as a significant marker for poor overall survival in ependymomas and H3 K28-mutant diffuse midline gliomas. Already being actively applied to other pediatric cancers and PNOC molecular tumor board decision-making, OpenPBTA is an invaluable resource to the pediatric oncology community.

Recursive splicing discovery using lariats in total RNA sequencing

Recursive splicing is a non-canonical splicing mechanism in which an intron is removed in segments via multiple splicing reactions. Relatively few recursive splice sites have been identified with high confidence in human introns, and more comprehensive analyses are needed to better characterize where recursive splicing happens and whether or not it has a regulatory function. In this study, we use an unbiased approach using intron lariats to search for recursive splice sites in constitutive introns and alternative exons in the human transcriptome. We find evidence for recursive splicing in a broader range of intron sizes than previously reported and detail a new location for recursive splicing at the distal ends of cassette exons. In addition, we identify evidence for the conservation of these recursive splice sites among higher vertebrates and the use of these sites to influence alternative exon exclusion. Together, our data demonstrate the prevalence of recursive splicing and its potential influence on gene expression through alternatively spliced isoforms.

The proto-oncogene SRC phosphorylates cGAS to inhibit an antitumor immune response

Cyclic GMP-AMP synthase (cGAS) is a DNA sensor and responsible for inducing an antitumor immune response. Recent studies reveal that cGAS is frequently inhibited in cancer, and therapeutic targets to promote antitumor cGAS function remain elusive. SRC is a proto-oncogene tyrosine kinase and is expressed at elevated levels in numerous cancers. Here, we demonstrate that SRC expression in primary and metastatic bladder cancer negatively correlates with innate immune gene expression and immune cell infiltration. We determine that SRC restricts cGAS signaling in human cell lines through SRC small molecule inhibitors, depletion, and overexpression. cGAS and SRC interact in cells and in vitro, while SRC directly inhibits cGAS enzymatic activity and DNA binding in a kinase-dependent manner. SRC phosphorylates cGAS, and inhibition of cGAS Y248 phosphorylation partially reduces SRC inhibition. Collectively, our study demonstrates that cGAS antitumor signaling is hindered by the proto-oncogene SRC and describes how cancer-associated proteins can regulate the innate immune system.

Over-expression of ADAR1 in mice does not initiate or accelerate cancer formation in vivo

Adenosine to inosine editing (A-to-I) in regions of double stranded RNA (dsRNA) is mediated by adenosine deaminase acting on RNA 1 (ADAR1) or ADAR2. ADAR1 and A-to-I editing levels are increased in many human cancers. Inhibition of ADAR1 has emerged as a high priority oncology target, however, whether ADAR1 overexpression enables cancer initiation or progression has not been directly tested. We established a series of in vivo models to allow overexpression of full-length ADAR1, or its individual isoforms, to test if increased ADAR1 expression was oncogenic. Widespread over-expression of ADAR1 or the p110 or p150 isoforms individually as sole lesions was well tolerated and did not result in cancer initiation. Therefore, ADAR1 overexpression alone is not sufficient to initiate cancer. We demonstrate that endogenous ADAR1 and A-to-I editing increased upon immortalization in murine cells, consistent with the observations from human cancers. We tested if ADAR1 over-expression could co-operate with cancer initiated by loss of tumour suppressors using a model of osteosarcoma. We did not see a disease potentiating or modifying effect of overexpressing ADAR1 or its isoforms in the models assessed. We conclude that increased ADAR1 expression and A-to-I editing in cancers is most likely a consequence of tumor formation.

The phenotype of the most common human ADAR1p150 Zα mutation P193A in mice is partially penetrant

ADAR1 -mediated A-to-I RNA editing is a self-/non-self-discrimination mechanism for cellular double-stranded RNAs. ADAR mutations are one cause of Aicardi-Goutières Syndrome, an inherited paediatric encephalopathy, classed as a "Type I interferonopathy." The most common ADAR1 mutation is a proline 193 alanine (p.P193A) mutation, mapping to the ADAR1p150 isoform-specific Zα domain. Here, we report the development of an independent murine P195A knock-in mouse, homologous to human P193A. The Adar1P195A/P195A mice are largely normal and the mutation is well tolerated. When the P195A mutation is compounded with an Adar1 null allele (Adar1P195A/- ), approximately half the animals are runted with a shortened lifespan while the remaining Adar1P195A/- animals are normal, contrasting with previous reports. The phenotype of the Adar1P195A/- animals is both associated with the parental genotype and partly non-genetic/environmental. Complementation with an editing-deficient ADAR1 (Adar1P195A/E861A ), or the loss of MDA5, rescues phenotypes in the Adar1P195A/- mice.

The targeted next-generation sequence revealed SMAD4, AKT1, and TP53 mutations from circulating cell-free DNA of breast cancer and its effect on protein structure - A computational approach

Breast cancer biomarkers that detect marginally advanced stages are still challenging. The detection of specific abnormalities, targeted therapy selection, prognosis, and monitoring of treatment effectiveness over time are all made possible by circulating free DNA (cfDNA) analysis. The proposed study will detect specific genetic abnormalities from the plasma cfDNA of a female breast cancer patient by sequencing a cancer-related gene panel (MGM455 - Oncotrack Ultima), including 56 theranostic genes (SNVs and small INDELs). Initially, we determined the pathogenicity of the observed mutations using PredictSNP, iStable, Align-GVGD, and ConSurf servers. As a next step, molecular dynamics (MD) was implemented to determine the functional significance of SMAD4 mutation (V465M). Lastly, the mutant gene relationships were examined using the Cytoscape plug-in GeneMANIA. Using ClueGO, we determined the gene's functional enrichment and integrative analysis. The structural characteristics of SMAD4 V465M protein by MD simulation analysis further demonstrated that the mutation was deleterious. The simulation showed that the native structure was more significantly altered by the SMAD4 (V465M) mutation. Our findings suggest that SMAD4 V465M mutation might be significantly associated with breast cancer, and other patient-found mutations (AKT1-E17K and TP53-R175H) are synergistically involved in the process of SMAD4 translocate to nuclease, which affects the target gene translation. Therefore, this combination of gene mutations could alter the TGF-β signaling pathway in BC. We further proposed that the SMAD4 protein loss may contribute to an aggressive phenotype by inhibiting the TGF-β signaling pathway. Thus, breast cancer's SMAD4 (V465M) mutation might increase their invasive and metastatic capabilities.Communicated by Ramaswamy H. Sarma.

Cost-effective genotyping for classical congenital adrenal hyperplasia (CAH) due to 21-hydroxylase deficiency (21-OHD) in resource-poor settings: multiplex ligation probe amplification (MLPA) with/without sequential next-generation sequencing (NGS)

PURPOSE

Genotyping of classic congenital adrenal hyperplasia (CAH) due to 21-hydroxylase deficiency (21-OHD) is becoming increasingly significant beyond prenatal counseling in the current era of emerging gene therapy/editing technologies. While the knowledge of common variants helps in designing cost-effective genotyping strategies, limited data are currently available from the Indian subcontinent, especially South India, mainly due to financial constraints. The aim of this study is to assess the genotype of individuals with classic CAH from a South Indian cohort in a cost-effective manner.

METHODS

The genotypes of 46 unrelated subjects with classic CAH were studied through initial multiplex ligation-dependent probe amplification (MLPA) using the SALSA MLPA Probe-mix P050 CAH (MRC Holland). Next-generation sequencing (NGS) was done in 10 subjects, as their MLPA was either negative or showed heterozygous variants.

RESULTS

The common variants observed in our study population of 46 subjects were large deletions (35.8%), intron 2 variant [c.293-13A/C > G] (35.8%), 8 bp del [c.332_339del p.(Gly111Valfs*21)] (7.7%), and R356W [c.1069 C > T p.(Arg357Trp)] (6.6%). MLPA alone detected pathogenic variants in 78.2% of the initial study samples (36/46). Sequential NGS resulted in a 100% detection rate in our study population.

CONCLUSION

MLPA appears to be an effective first genotyping modality for this South Indian cohort due to the high prevalence of large deletions and common variants. MLPA as a first initial screening genotyping test with sequential NGS when required may be a cost-effective and highly sensitive approach to CYP21A2 genotyping in our part of the world and in resource-poor settings.

An integrated investigation of structural and pathway alteration caused by PIK3CA and TP53 mutations identified in cfDNA of metastatic breast cancer

In peripheral blood, cell-free DNA (cfDNA) contains circulating tumor DNA (ctDNA), which indicates molecular abnormalities in metastatic breast tumor tissue. The sequencing of cfDNA of Metastatic Breast Cancer (MBC) patients allows assessment of therapy response and noninvasive treatment. In the proposed study, clinically significant alterations in PIK3CA and TP53 genes associated with MBC resulting in a missense substitution of His1047Arg and Arg282Trp from an next-generation sequencing-based multi-gene panel were reported in a cfDNA of a patient with MBC. To investigate the impact of the reported mutation, we used molecular docking, molecular dynamics simulation, network analysis, and pathway analysis. Molecular Docking analysis determined the distinct binding pattern revealing H1047R-ATP complex has a higher number of Hydrogen bonds (H-bonds) and binding affinity with a slight difference compared to the PIK3CA-ATP complex. Following, molecular dynamics simulation for 200 ns, of which H1047R-ATP complex resulted in the instability of PIK3CA. Similarly, for TP53 mutant R282W, the zinc-free state (apo) and zinc-bounded (holo) complexes were investigated for conformational change between apo and holo complexes, of which the holo complex mutant R282W was unstable. To validate the conformational change of PIK3CA and TP53, 80% mutation of H1047R in the kinase domain of p110α expressed ubiquitously in PIK3CA protein that alters PI3K pathway, while R282W mutation in DNA binding helix (H2) region of P53 protein inhibits the transcription factor in P53 pathway causing MBC. According to our findings, the extrinsic (hypoxia, oxidative stress, and acidosis); intrinsic factors (MYC amplification) in PIK3CA and TP53 mutations will provide potential insights for developing novel therapeutic methods for MBC therapy.

Modulation of RNA splicing enhances response to BCL2 inhibition in leukemia

Therapy resistance is a major challenge in the treatment of cancer. Here, we performed CRISPR-Cas9 screens across a broad range of therapies used in acute myeloid leukemia to identify genomic determinants of drug response. Our screens uncover a selective dependency on RNA splicing factors whose loss preferentially enhances response to the BCL2 inhibitor venetoclax. Loss of the splicing factor RBM10 augments response to venetoclax in leukemia yet is completely dispensable for normal hematopoiesis. Combined RBM10 and BCL2 inhibition leads to mis-splicing and inactivation of the inhibitor of apoptosis XIAP and downregulation of BCL2A1, an anti-apoptotic protein implicated in venetoclax resistance. Inhibition of splicing kinase families CLKs (CDC-like kinases) and DYRKs (dual-specificity tyrosine-regulated kinases) leads to aberrant splicing of key splicing and apoptotic factors that synergize with venetoclax, and overcomes resistance to BCL2 inhibition. Our findings underscore the importance of splicing in modulating response to therapies and provide a strategy to improve venetoclax-based treatments.

How to survive in the world's third poplar: Insights from the genome of the highest altitude woody plant, Hippophae tibetana (Elaeagnaceae)

Hippophae tibetana (Tibetan sea-buckthorn) is one of the highest distributed woody plants in the world (3,000-5,200 meters a.s.l.). It is characterized by adaptation to extreme environment and important economic values. Here, we combined PacBio Hifi platform and Hi-C technology to assemble a 1,452.75 Mb genome encoding 33,367 genes with a Contig N50 of 74.31 Mb, and inferred its sexual chromosome. Two Hippophae-specific whole-genome duplication events (18.7-21.2 million years ago, Ma; 28.6-32.4 Ma) and long terminal repeats retroelements (LTR-RTs) amplifications were detected. Comparing with related species at lower altitude, Ziziphus jujuba (<1, 700 meters a.s.l.), H. tibetana had some significantly rapid evolving genes involved in adaptation to high altitude habitats. However, comparing with Hippophae rhamnoides (<3, 700 meters a.s.l.), no rapid evolving genes were found except microtubule and microtubule-based process genes, H. tibetana has a larger genome, with extra 2, 503 genes (7.5%) and extra 680.46 Mb transposable elements (TEs) (46.84%). These results suggest that the changes in the copy number and regulatory pattern of genes play a more important role for H. tibetana adapting to more extreme and variable environments at higher altitude by more TEs and more genes increasing genome variability and expression plasticity. This suggestion was supported by two findings: nitrogen-fixing genes of H. tibetana having more copies, and intact TEs being significantly closer genes than fragmentary TEs. This study provided new insights into the evolution of alpine plants.

Molecular evaluation of gene mutation profiles and copy number variations in pediatric acute myeloid leukemia

BACKGROUND

The objectives of this study were to investigate the mutation profiles of targeted genes and copy number variations (CNVs) in normal cytogenetics (CN) pediatric acute myeloid leukemia (AML).

METHODS

This prospective study was conducted from October 2018 to December 2020. The next-generation sequencing (NGS) and chromosomal microarray analyses (CMA) were performed in pediatric CN-AML patients.

RESULTS

Out of 94 children (aged ≤18 years), 70 patients with AML (24 excluded) underwent conventional karyotyping/cytogenetic analyses. Forty-five (64.3%) of patients had abnormal/ recurrent cytogenetic abnormalities and 25 (35.7%) had normal cytogenetics. Twenty-three out of 25 CN-AML were further processed for gene mutation profile and CNVs using NGS and CMA, respectively. Twenty-two out of 23 (95.7%) patients were detected to have mutations in various genes. The common mutations were: NRAS, NPM1, CEBPA, KRAS, KIT, RUNX1, NOTCH1, WT1, GATA1, GATA2, FLT3, KMT2D, FLT3-TKD, and PHF6. Copy number variations (CNVs) were detected in nine patients (39%), and eight (34.8%) had a long contiguous stretch of homozygosity (LCSH) /loss of heterozygosity (LOH). An LCSH was detected on chromosomes 5, 7, 11, and 19. The gains were more common than losses (8 vs 2). The gains were observed on chromosomes 8, 9, 14, 19, 21, and 22, and the losses were detected on chromosomes 7 and 10. Monosomy was observed in three patients. Three patients (monosomy7, n = 2, and FLT-ITD, n = 1) were reclassified into the high-risk category. Post-induction, complete remission was achieved in all evaluable patients.

CONCLUSION

CN-AML patients have genetic abnormalities that can be detected by more advanced techniques like NGS and CMA. These genetic abnormalities play a role in risk stratification that may remain hidden in otherwise CN-AML.

Next-Generation Sequencing Highlights of Diffuse Large B-cell Lymphoma in a Tertiary Care Hospital in North India

Introduction: Next-generation sequencing (NGS) elucidates the diffuse large B-cell lymphoma (DLBCL) genetic characteristics by finding recurrent and novel somatic mutations. This observational study attempted to create an NGS panel with a focus on identifying novel somatic mutations which could have potential clinical and therapeutic implications. This panel was created to look for mutations in 133 genes chosen on basis of a literature review and it was used to sequence the tumor DNA of 20 DLBCL patients after a centralized histopathologic review.

Methods: The study included 20 patients having DLBCL. The quality and quantity of tumor cells were accessed by H&E staining and correlated with histopathology and Immunohistochemistry (IHC) status. Patients were grouped as ABC (activated B-cell), PMBL (primary mediastinal large B-cell lymphoma), and other or unclassified subtypes. The lymphoma panel of 133 was designed on targeted sequencing of multiple genes for the coding regions through NGS. The libraries were prepared and sequenced using the Illumina platform. The alignment of obtained sequences was performed using Burrows-Wheeler Aligner and identification of somatic mutations was done using LoFreq (version 2) variant caller. The mutations were annotated using an annotation pipeline (VariMAT). Previously published literature and databases were used for the annotation of clinically relevant mutations. The common variants were filtered for reporting based on the presence in various population databases (1000G, ExAC, EVS, 1000Japanese, dbSNP, UK10K, MedVarDb). A custom read-depth-based algorithm was used to determine CNV (Copy Number Variants) from targeted sequencing experiments. Rare CNVs were detected using a comparison of the test data read-depths with the matched reference dataset. Reportable mutations were prioritized and prepared based on AMP-ASCO-CAP (Association for Molecular Pathology-American Society of Clinical Oncology-College of American Pathologists), WHO guidelines, and also based on annotation metrics from OncoMD (a knowledge base of genomic alterations).

Results: The informativity of the panel was 95 percent. NOTCH 1 was the most frequently mutated gene in 16.1% of patients followed by 12.9% who had ARID1A mutations. MYD88 and TP53 mutations were detected in 9.6% of the patient while 6.4% of patients had CSF3R mutations. NOTCH 1 and TP 53 are the most frequently reported gene in the middle age group (40-60). Mutation in MYD88 is reported in every age group. MYD88 (51%) is the most common mutation in ABC subtypes of DLBCL, followed by NOTCH 1 (44%) and SOCS 1 (33%) according to our findings. NOTCH 1 mutations are frequent in ABC and PMBL subtypes. Closer investigation reveals missense mutation is the most frequent mutation observed in the total cohort targeting 68.4% followed by frameshift deletion reported in 26.3%. Six novel variants have been discovered in this study.

Conclusions: This study demonstrates the high yield of information in DLBCL using the NGS Lymphoma panel. Results also highlight the molecular heterogeneity of DLBCL subtypes which indicates the need for further studies to make the results of the NGS more clinically relevant.

Integrated DNA Copy Number and Expression Profiling Identifies IGF1R as a Prognostic Biomarker in Pediatric Osteosarcoma

Osteosarcoma is a primary malignant bone tumor arising from bone-forming mesenchymal cells in children and adolescents. Despite efforts to understand the biology of the disease and identify novel therapeutics, the survival of osteosarcoma patients remains dismal. We have concurrently profiled the copy number and gene expression of 226 osteosarcoma samples as part of the Strategic Partnering to Evaluate Cancer Signatures (SPECS) initiative. Our results demonstrate the heterogeneous landscape of osteosarcoma in younger populations by showing the presence of genome-wide copy number abnormalities occurring both recurrently among samples and in a high frequency. Insulin growth factor receptor 1 (IGF1R) is a receptor tyrosine kinase which binds IGF1 and IGF2 to activate downstream pathways involved in cell apoptosis and proliferation. We identify prevalent amplification of IGF1R corresponding with increased gene expression in patients with poor survival outcomes. Our results substantiate previously tenuously associated copy number abnormalities identified in smaller datasets (13q34+, 20p13+, 4q35-, 20q13.33-), and indicate the significance of high fibroblast growth factor receptor 2 (FGFR2) expression in distinguishing patients with poor prognosis. FGFR2 is involved in cellular proliferation processes such as division, growth and angiogenesis. In summary, our findings demonstrate the prognostic significance of several genes associated with osteosarcoma pathogenesis.

Synthetic introns enable splicing factor mutation-dependent targeting of cancer cells

Many cancers carry recurrent, change-of-function mutations affecting RNA splicing factors. Here, we describe a method to harness this abnormal splicing activity to drive splicing factor mutation-dependent gene expression to selectively eliminate tumor cells. We engineered synthetic introns that were efficiently spliced in cancer cells bearing SF3B1 mutations, but unspliced in otherwise isogenic wild-type cells, to yield mutation-dependent protein production. A massively parallel screen of 8,878 introns delineated ideal intronic size and mapped elements underlying mutation-dependent splicing. Synthetic introns enabled mutation-dependent expression of herpes simplex virus-thymidine kinase (HSV-TK) and subsequent ganciclovir (GCV)-mediated killing of SF3B1-mutant leukemia, breast cancer, uveal melanoma and pancreatic cancer cells in vitro, while leaving wild-type cells unaffected. Delivery of synthetic intron-containing HSV-TK constructs to leukemia, breast cancer and uveal melanoma cells and GCV treatment in vivo significantly suppressed the growth of these otherwise lethal xenografts and improved mouse host survival. Synthetic introns provide a means to exploit tumor-specific changes in RNA splicing for cancer gene therapy.

Comprehensive analysis of genetic factors predicting overall survival in Myelodysplastic syndromes

Myelodysplastic syndromes (MDS) are a group of clonal hematological disease with high risk of progression to AML. Accurate risk stratification is of importance for the proper management of MDS. Genetic lesions (Cytogenetic and Molecular mutations) are known to help in prognosticating the MDS patients. We have studied 152 MDS patients using cytogenetics and next generation sequencing (NGS). These patients were evaluated and as per cytogenetic prognostic group, majority (92.1%) of the patients classified as good (81.6%) and intermediate (10.5%) group. The NGS identified 38 different gene mutations in our cohort. Among 111 MDS patients with mutations, the most frequent mutated genes were SF3B1 (25.2%), SRSF2 (19%) U2AF1 (14.4%) ASXL1 (9.9%) RUNX1 (9.9%) TET2 (9%), TP53 (9%), ATM (6.3%), NRAS (5.4%) and JAK2/3 (5.4%). The survival analysis revealed that the mutations in TP53, JAK2/3, KRAS, NRAS and ASXL1 were significantly (P < 0.05) associated with poor survival of the patients. The univariate cox and multivariate cox analysis of our study suggested that the age, marrow morphology, cytogenetic and gene mutations with IPSS-R should be considered for prognosticating the MDS patients. We have proposed M-IPSS-R which changed the risk stratification i.e. 66.3% patients had decreased risk whereas 33.75% showed increased risk compared to IPSS-R. The survival analysis also showed that the M-IPSS-R were more significant in separating the patients as per their risk than the IPSS-R alone. The change in risk stratification could help in proper strategy for the treatment planning.

Mycobacterium leprae Transcriptome During In Vivo Growth and Ex Vivo Stationary Phases

Mycobacterium leprae, the causative agent of leprosy, is an obligate intracellular pathogen primarily residing within host macrophages and Schwann cells. Whole genome sequencing predicts a highly degraded genome with approximately one third of the coding capacity resulting in the loss of many catabolic pathways. Therefore, it can be assumed that M. leprae obtains many of the necessary metabolites for intracellular survival and growth from the host cells. In this study, global transcriptomic analyses were done on freshly harvested M. leprae growing in athymic mouse footpads for five months (MFP5) and compared to those held in axenic medium for 48 (ML48) and 96 (ML96) hours. Results show that all of the genes and pseudogenes were transcribed under both in vivo and in vitro conditions. 24% and 33% of gene transcript levels were significantly altered in ML48 and ML96 respectively, compared to MFP5. Approximately 45% (39/86) of lipid metabolism genes were significantly downregulated in ML96 compared to MFP5, majority of which are in the β-oxidation pathway. Cholesterol oxidase, acyl-CoA dehydrogenase, and coenzyme F420-dependent oxidoreductase, were significantly upregulated in both ML48 and ML96 compared to MFP5. 30% of cell wall and cell processes functional category genes had altered gene transcription at 96hr compared to MFP5. 40% of 57 genes associated with mycobacterial virulence showed significantly altered transcript levels with 52% significantly downregulated in ML96, including most of the Pro-Glu/Pro-Pro-Glu genes. All 111 hypothetical protein genes with unknown function were expressed. Adenosine triphosphate (ATP) synthesis in M. leprae appears to be significantly downregulated under ex vivo conditions. This is the first study comparing M. leprae global gene expression during in vivo growth and ex vivo stationery phase in axenic medium confirming that during the growth phase in the footpads of experimentally infected mice, M. leprae is metabolically active and its primary source of energy production is probably lipids.

Ribosomal leaky scanning through a translated uORF requires eIF4G2

eIF4G2 (DAP5 or Nat1) is a homologue of the canonical translation initiation factor eIF4G1 in higher eukaryotes but its function remains poorly understood. Unlike eIF4G1, eIF4G2 does not interact with the cap-binding protein eIF4E and is believed to drive translation under stress when eIF4E activity is impaired. Here, we show that eIF4G2 operates under normal conditions as well and promotes scanning downstream of the eIF4G1-mediated 40S recruitment and cap-proximal scanning. Specifically, eIF4G2 facilitates leaky scanning for a subset of mRNAs. Apparently, eIF4G2 replaces eIF4G1 during scanning of 5′ UTR and the necessity for eIF4G2 only arises when eIF4G1 dissociates from the scanning complex. In particular, this event can occur when the leaky scanning complexes interfere with initiating or elongating 80S ribosomes within a translated uORF. This mechanism is therefore crucial for higher eukaryotes which are known to have long 5′ UTRs with highly frequent uORFs. We suggest that uORFs are not the only obstacle on the way of scanning complexes towards the main start codon, because certain eIF4G2 mRNA targets lack uORF(s). Thus, higher eukaryotes possess two distinct scanning complexes: the principal one that binds mRNA and initiates scanning, and the accessory one that rescues scanning when the former fails.

Targeted Exome Sequencing in South Indian patients with Familial Hypercholesterolemia

BACKGROUND

Familial hypercholesterolemia (FH) is an autosomal dominant genetic disorder with elevated LDL-C levels which can ultimately lead to premature Coronary Artery Disease (CAD).

OBJECTIVES

In presence of limited genetic data on FH in India, the present study was aimed to determine the mutation spectrum in Indian FH patients using a targeted exome sequencing.

METHODS

54 FH cases (31 index cases + 23 extended family members) were categorized according to Dutch Lipid Clinic Network Criteria (DLCNC). Targeted exome sequencing was performed using 23 gene panel associated with lipid metabolism.

RESULTS

All subjects showed the presence of family history of CAD, 38(70%) patients had corneal arcus whereas only 06(11%) subjects had xanthomas. As per the DLCNC, definite, probable, possible and unlikely FH were 48%, 30%, 11% and 11% respectively. Mutations were observed in 12 of the 23 gene panel with CETP, APOA5, EPHX2 and SREBP2 genes were identified for the first time in Indian FH patients. All 19 mutations including a novel frame-shift mutation in LDLR gene were reported for the first time in Indian FH patients. These mutations were identified in 28(52%) subjects and interestingly ∼73% of the clinically identified FH patients didn't harbour mutations in FH classical genes (LDLR, ApoB, PCSK9).

CONCLUSION

This is the first study in the South Indian FH patients to perform targeted exome sequencing. Absence of mutations in the FH classical genes strongly indicates the polygenic nature of FH, further underscoring the importance of targeted exome sequencing for identifying mutations in genetically diverse Indian population.

Super-resolution microscopy reveals stochastic initiation of replication in Drosophila polytene chromosomes

Studying the probability distribution of replication initiation along a chromosome is a huge challenge. Drosophila polytene chromosomes in combination with super-resolution microscopy provide a unique opportunity for analyzing the probabilistic nature of replication initiation at the ultrastructural level. Here, we developed a method for synchronizing S-phase induction among salivary gland cells. An analysis of the replication label distribution in the first minutes of S phase and in the following hours after the induction revealed the dynamics of replication initiation. Spatial super-resolution structured illumination microscopy allowed identifying multiple discrete replication signals and to investigate the behavior of replication signals in the first minutes of the S phase at the ultrastructural level. We identified replication initiation zones where initiation occurs stochastically. These zones differ significantly in the probability of replication initiation per time unit. There are zones in which initiation occurs on most strands of the polytene chromosome in a few minutes. In other zones, the initiation on all strands takes several hours. Compact bands are free of replication initiation events, and the replication runs from outer edges to the middle, where band shapes may alter.

Parity-induced changes to mammary epithelial cells control NKT cell expansion and mammary oncogenesis

Pregnancy reprograms mammary epithelial cells (MECs) to control their responses to pregnancy hormone re-exposure and carcinoma progression. However, the influence of pregnancy on the mammary microenvironment is less clear. Here, we used single-cell RNA sequencing to profile the composition of epithelial and non-epithelial cells in mammary tissue from nulliparous and parous female mice. Our analysis indicates an expansion of γδ natural killer T-like immune cells (NKTs) following pregnancy and upregulation of immune signaling molecules in post-pregnancy MECs. We show that expansion of NKTs following pregnancy is due to elevated expression of the antigen-presenting molecule CD1d on MECs. Loss of CD1d expression on post-pregnancy MECs, or overall lack of activated NKTs, results in mammary oncogenesis. Collectively, our findings illustrate how pregnancy-induced changes modulate the communication between MECs and the immune microenvironment and establish a causal link between pregnancy, the immune microenvironment, and mammary oncogenesis.

Comparison of RNA content from hydrophobic interaction chromatography-isolated seminal plasma exosomes from intrauterine insemination (IUI) pregnancies

Male factors account for roughly half of infertility cases, with most male infertility diagnosed as idiopathic. Researchers predicting intrauterine insemination success rates have identified multiple prognostic factors, including semen parameters and seminal fluid composition. Seminal plasma contains extracellular exosomes, which contain RNAs and proteins involved in spermatogenesis. The contents of seminal plasma exosomes may be an indicator of overall sperm quality or fertility potential; therefore, analysis of exosomes may provide a measure for sperm viability and fertilization potential. In our study, exosomes were isolated and purified from seminal plasma obtained from IUI treatments with known pregnancy outcomes. We used a unique method to isolate the exosomes which made use of the hydrophobic interaction chromatography method. RNASeq was performed on RNAs from the purified exosomes. This analysis revealed holistic trends, including increased expression associated with RNA originating from chromosomes 1, 10, 12, 16 and 21. Overall, total RNA was significantly decreased and rRNA was significantly increased in successful IUI attempts. Furthermore, we found specific mRNAs and lincRNAs associated with positive versus negative pregnancy outcomes. Our study isolated and purified seminal plasma exosomes without ultracentrifugation, and it provides further evidence for differences in seminal plasma exosome molecular contents associated with pregnancy status.

The landscape of coding RNA editing events in pediatric cancer

BACKGROUND

RNA editing leads to post-transcriptional variation in protein sequences and has important biological implications. We sought to elucidate the landscape of RNA editing events across pediatric cancers.

METHODS

Using RNA-Seq data mapped by a pipeline designed to minimize mapping ambiguity, we investigated RNA editing in 711 pediatric cancers from the St. Jude/Washington University Pediatric Cancer Genome Project focusing on coding variants which can potentially increase protein sequence diversity. We combined de novo detection using paired tumor DNA-RNA data with analysis of known RNA editing sites.

RESULTS

We identified 722 unique RNA editing sites in coding regions across pediatric cancers, 70% of which were nonsynonymous recoding variants. Nearly all editing sites represented the canonical A-to-I (n = 706) or C-to-U sites (n = 14). RNA editing was enriched in brain tumors compared to other cancers, including editing of glutamate receptors and ion channels involved in neurotransmitter signaling. RNA editing profiles of each pediatric cancer subtype resembled those of the corresponding normal tissue profiled by the Genotype-Tissue Expression (GTEx) project.

CONCLUSIONS

In this first comprehensive analysis of RNA editing events in pediatric cancer, we found that the RNA editing profile of each cancer subtype is similar to its normal tissue of origin. Tumor-specific RNA editing events were not identified indicating that successful immunotherapeutic targeting of RNA-edited peptides in pediatric cancer should rely on increased antigen presentation on tumor cells compared to normal but not on tumor-specific RNA editing per se.

Utilization of Whole Exome Sequencing in Non-Syndromic Premature Ovarian Failure: Ficolin-3 Gene Mutation in an Iranian Family

Background

Premature ovarian failure is a heterogeneous disorder, leading to early menopause. Several genes have been identified as the cause of non-syndromic premature ovarian failure (POF). Our aim was to explore the genetic defects in Iranian patients with POF.

Methods

We studied a family with three females exhibiting non-syndromic POF. WES was performed for one of the affected individuals after ruling out the presence of CGG repeat expansion at fragile X mental retardation 1 gene in the family. Sanger sequencing was used to confirm the candidate sequence variants in the proband, and screening of the detected mutation was performed for the other affected and unaffected members of the family.

Results

A homozygous frameshift mutation, c.349delC, was identified in ficolin-3 (FCN3) gene in the proband and two other patients. The parents and two healthy brothers were heterozygous for the mutation, and an unaffected sister was homozygous for wild type.

Conclusion

This is the first report of a mutation in FCN3 gene in a family with POF. Our findings can lead to the enhancement of genetic databases of patients with POF, specifically for families with high-risk background.

ZBTB33 is mutated in clonal hematopoiesis and myelodysplastic syndromes and impacts RNA splicing

Clonal hematopoiesis results from somatic mutations in cancer driver genes in hematopoietic stem cells. We sought to identify novel drivers of clonal expansion using an unbiased analysis of sequencing data from 84,683 persons and identified common mutations in the 5-methylcytosine reader, ZBTB33, as well as in YLPM1, SRCAP, and ZNF318. We also identified these mutations at low frequency in myelodysplastic syndrome patients. Zbtb33 edited mouse hematopoietic stem and progenitor cells exhibited a competitive advantage in vivo and increased genome-wide intron retention. ZBTB33 mutations potentially link DNA methylation and RNA splicing, the two most commonly mutated pathways in clonal hematopoiesis and MDS.

Retinoblastoma genetics screening and clinical management

BACKGROUND

India accounts for 20% of the global retinoblastoma (RB) burden. However, the existing data on RB1 gene germline mutations and its influence on clinical decisions is minimally explored.

METHODS

Fifty children with RB underwent complete clinical examination and appropriate multidisciplinary management. Screening of germline RB1 gene mutations was performed through next-generation sequencing and Multiplex Ligation-dependent Probe Amplification (MLPA) analysis. The mutation and non-mutation groups were compared for clinical parameters especially severity, progression and recurrence.

RESULTS

Twenty-nine patients had bilateral RB (BLRB) and 21 had unilateral RB (ULRB). The genetic analysis revealed 20 RB1 variations in 29 probands, inclusive of 3 novel mutations, known 16 mutations and heterozygous whole gene deletions. The mutation detection rate (MDR) was 86.2% in BLRB and 19% in ULRB. Associations of disease recurrence (p = 0.021), progression (p = 0.000) and higher percentage of optic nerve invasion, subretinal seeds and high-risk pathological factors were observed in the mutation group. Clinical management was influenced by the presence of germline mutations, particularly while deciding on enucleation, frequency of periodic follow up and radiotherapy.

CONCLUSIONS

We identified novel RB1 mutations, and our mutation detection rate was on par with the previous global studies. In our study, genetic results influenced clinical management and we suggest that it should be an essential and integral component of RB-care in India and elsewhere.

Pharmacologic modulation of RNA splicing enhances anti-tumor immunity

Although mutations in DNA are the best-studied source of neoantigens that determine response to immune checkpoint blockade, alterations in RNA splicing within cancer cells could similarly result in neoepitope production. However, the endogenous antigenicity and clinical potential of such splicing-derived epitopes have not been tested. Here, we demonstrate that pharmacologic modulation of splicing via specific drug classes generates bona fide neoantigens and elicits anti-tumor immunity, augmenting checkpoint immunotherapy. Splicing modulation inhibited tumor growth and enhanced checkpoint blockade in a manner dependent on host T cells and peptides presented on tumor MHC class I. Splicing modulation induced stereotyped splicing changes across tumor types, altering the MHC I-bound immunopeptidome to yield splicing-derived neoepitopes that trigger an anti-tumor T cell response in vivo. These data definitively identify splicing modulation as an untapped source of immunogenic peptides and provide a means to enhance response to checkpoint blockade that is readily translatable to the clinic.

Defective folate metabolism causes germline epigenetic instability and distinguishes Hira as a phenotype inheritance biomarker

The mechanism behind transgenerational epigenetic inheritance is unclear, particularly through the maternal grandparental line. We previously showed that disruption of folate metabolism in mice by the Mtrr hypomorphic mutation results in transgenerational epigenetic inheritance of congenital malformations. Either maternal grandparent can initiate this phenomenon, which persists for at least four wildtype generations. Here, we use genome-wide approaches to reveal genetic stability in the Mtrr model and genome-wide differential DNA methylation in the germline of Mtrr mutant maternal grandfathers. We observe that, while epigenetic reprogramming occurs, wildtype grandprogeny and great grandprogeny exhibit transcriptional changes that correlate with germline methylation defects. One region encompasses the Hira gene, which is misexpressed in embryos for at least three wildtype generations in a manner that distinguishes Hira transcript expression as a biomarker of maternal phenotypic inheritance.

Utility of clinical exome sequencing in the evaluation of neonates with suspected genetic condition - An observational study from tertiary neonatal care unit in South India

OBJECTIVES

To study the utility of clinical exome sequencing (CES) using next generation sequencing (NGS) in evaluating neonates with suspected genetic conditions.

METHODS

This is an observational study conducted in a tertiary care neonatal unit. We included neonates with suspected genetic conditions, for whom CES were done either by direct sampling or from stored DNA. Data was collected from the Sri Ramachandra centre of excellence in perinatal health (SCOPE) case records of 2016-2019. Yield of CES, percentage of pathogenic, non-pathogenic and variant of uncertain significance (VUS) and associated disorders were studied.

RESULTS

CES was done in 36 neonates. Variants were detected in 78% (28/36). However, significant variants with clinical correlation were present in 20 (56%) babies. Test was carried out from the stored sample in 10 (28%) babies. Mean turn-around time was 39 ± 7 days. Specialist was involved in 1 and treatment changes were done in 5 neonates. Five out of 8 VUS were clinically correlating. Inborn errors of metabolism were the commonest (60%). Two VUS were ascertained as likely pathogenic after parental segregation analysis.

CONCLUSION

CES has a definite role in evaluation of suspected genetic conditions for diagnosis and prognostication. It also helps scientific society to build in additional evidence so that the "VUS" could be asserted as "likely pathogenic" . Our experience reiterates the importance of storing and archiving DNA of the affected child.

Minor intron retention drives clonal hematopoietic disorders and diverse cancer predisposition

Most eukaryotes harbor two distinct pre-mRNA splicing machineries: the major spliceosome, which removes >99% of introns, and the minor spliceosome, which removes rare, evolutionarily conserved introns. Although hypothesized to serve important regulatory functions, physiologic roles of the minor spliceosome are not well understood. For example, the minor spliceosome component ZRSR2 is subject to recurrent, leukemia-associated mutations, yet functional connections among minor introns, hematopoiesis and cancers are unclear. Here, we identify that impaired minor intron excision via ZRSR2 loss enhances hematopoietic stem cell self-renewal. CRISPR screens mimicking nonsense-mediated decay of minor intron-containing mRNA species converged on LZTR1, a regulator of RAS-related GTPases. LZTR1 minor intron retention was also discovered in the RASopathy Noonan syndrome, due to intronic mutations disrupting splicing and diverse solid tumors. These data uncover minor intron recognition as a regulator of hematopoiesis, noncoding mutations within minor introns as potential cancer drivers and links among ZRSR2 mutations, LZTR1 regulation and leukemias.

N, Ramprasad VL, Gupta R, Seshagiri S, Nongthomba U, Phalke S. Matrilineal analysis of mutations in the DMD gene in a multigenerational South Indian cohort using DMD gene panel sequencing

BACKGROUND

Duchenne muscular dystrophy (DMD) is an X-linked recessive neuromuscular disorder characterised by progressive irreversible muscle weakness, primarily of the skeletal and the cardiac muscles. DMD is characterised by mutations in the dystrophin gene, resulting in the absence or sparse quantities of dystrophin protein. A precise and timely molecular detection of DMD mutations encourages interventions such as carrier genetic counselling and in undertaking therapeutic measures for the DMD patients.

RESULTS

In this study, we developed a 2.1 Mb custom DMD gene panel that spans the entire DMD gene, including the exons and introns. The panel also includes the probes against 80 additional genes known to be mutated in other muscular dystrophies. This custom DMD gene panel was used to identify single nucleotide variants (SNVs) and large deletions with precise breakpoints in 77 samples that included 24 DMD patients and their matrilineage across four generations. We used this panel to evaluate the inheritance pattern of DMD mutations in maternal subjects representing 24 DMD patients.

CONCLUSION

Here we report our observations on the inheritance pattern of DMD gene mutations in matrilineage samples across four generations. Additionally, our data suggest that the DMD gene panel designed by us can be routinely used as a single genetic test to identify all DMD gene variants in DMD patients and the carrier mothers.

Genomic insights into population history and biological adaptation in Oceania

The Pacific region is of major importance for addressing questions regarding human dispersals, interactions with archaic hominins and natural selection processes1. However, the demographic and adaptive history of Oceanian populations remains largely uncharacterized. Here we report high-coverage genomes of 317 individuals from 20 populations from the Pacific region. We find that the ancestors of Papuan-related ('Near Oceanian') groups underwent a strong bottleneck before the settlement of the region, and separated around 20,000-40,000 years ago. We infer that the East Asian ancestors of Pacific populations may have diverged from Taiwanese Indigenous peoples before the Neolithic expansion, which is thought to have started from Taiwan around 5,000 years ago2-4. Additionally, this dispersal was not followed by an immediate, single admixture event with Near Oceanian populations, but involved recurrent episodes of genetic interactions. Our analyses reveal marked differences in the proportion and nature of Denisovan heritage among Pacific groups, suggesting that independent interbreeding with highly structured archaic populations occurred. Furthermore, whereas introgression of Neanderthal genetic information facilitated the adaptation of modern humans related to multiple phenotypes (for example, metabolism, pigmentation and neuronal development), Denisovan introgression was primarily beneficial for immune-related functions. Finally, we report evidence of selective sweeps and polygenic adaptation associated with pathogen exposure and lipid metabolism in the Pacific region, increasing our understanding of the mechanisms of biological adaptation to island environments.

New susceptibility alleles associated with severe coronary artery stenosis in the Lebanese population

BACKGROUND

Coronary Artery Disease (CAD) is the narrowing or blockage of the coronary arteries. It is closely associated with numerous genetics and environmental factors that have been extensively evaluated in various populations. In recent studies, severe phenotypes have been strongly linked to genetic risk factors.

METHODS

This study investigated the association of clinical, demographic, and genetic factors with severe coronary artery stenosis phenotypes in our population composed of 1734 individuals with severe coronary stenosis (≥ 50% in coronary vessels) and comparing them to 757 controls with no evidence of stenosis on angiography. We performed generalized linear model (GLM) genome-wide association studies to evaluate three stratification models and their associations to characteristics of the clinical disease. In model 1, patients were not stratified. In model 2, patients were stratified based on presence or absence of CAD family history (FxCAD). In model 3, patients were stratified by young age of CAD onset.

RESULTS

Eight SNPs (single nucleotide polymorphism) were significantly associated with severe CAD phenotypes in the various models [Formula: see text], four of these SNPs were associated with severe CAD and the four others were specifically significant for young CAD patients. While these SNPs were not previously reported for association with CAD, six of them are present in genes that have already been linked to coronary disease.

CONCLUSION

In conclusion, this study presents new genetic factors associated with severe stenosis and highlights different risk factors associated with a young age at diagnosis of CAD.

Homozygous Missense Mutation on Exon 22 of PKHD1 Gene Causing Fatal Autosomal Recessive Polycystic Kidney Disease

Autosomal recessive polycystic kidney disease, described as a congenital hepatorenal fibrocystic syndrome, is a significant inherited cause of end stage renal failure in children with reported incidence of 1 in 20,000 live births. The clinical spectrum is wide. Antenatal findings of echogenic reniform enlarged kidneys associated with evidence of intrauterine renal failure in the form of severe oligoamnios are pathognomonic. Postnatal illness ranges from fatal respiratory failure due to pulmonary hypoplasia in neonates to chronic kidney disease in children, or later presentation of ductal plate malformation and portal hypertension. Advances in genetic diagnostic techniques have allowed recognition of genotypes. We report a novel homozygous missense variant on exon 22 of PKHD1 gene (chr6:51915067G > A; c.2167C > T) that results in the amino acid substitution of cysteine for arginine at codon 723 (p.Arg723Cys). The affected neonate presented with antenatal anhydramnios, classical radiological features, and severe hypoxic respiratory failure likely due to pulmonary hypoplasia and succumbed. The parents were found to be heterozygous carriers. Detection of the specific variant in the proband facilitated prenatal investigation in the next pregnancy.

Genome-wide ancestry and introgression in a Zambian baboon hybrid zone

Hybridization in nature offers unique insights into the process of natural selection in incipient species and their hybrids. In order to evaluate the patterns and targets of selection, we examine a recently discovered baboon hybrid zone in the Kafue River Valley of Zambia, where Kinda baboons (Papio kindae) and grey-footed chacma baboons (P. ursinus griseipes) coexist with hybridization. We genotyped baboons at 14,962 variable genome-wide autosomal markers using double-digest RADseq. We compared ancestry patterns from this genome-wide data set to previously reported ancestry from mitochondrial-DNA and Y-chromosome sources. We also fit a Bayesian genomic cline model to scan for genes with extreme patterns of introgression. We show that the Kinda baboon Y chromosome has penetrated the species boundary to a greater extent than either mitochondrial DNA or the autosomal chromosomes. We also find evidence for overall restricted introgression in the JAK/STAT signalling pathway. Echoing results in other species including humans, we find evidence for enhanced and/or directional introgression of immune-related genes or pathways including the toll-like receptor pathway, the blood coagulation pathway, and the LY96 gene. Finally we show enhanced introgression and excess chacma baboon ancestry in the sperm tail gene ODF2. Together, our results elucidate the dynamics of introgressive hybridization in a primate system while identifying genes and pathways possibly under selection.

Fluorimetric detection of methylated DNA of Sept9 promoter by silver nanoclusters at intrastrand 6C-loop

Methylation of DNA at carbon 5 of cytosines is the most common epigenetic modification of human genome. Due to its critical role in many normal cell processes such as growth and development, any aberrant methylation pattern in a particular locus may lead to abnormal functions and diseases such as cancer. Development of methods to detect methylation state of DNA which may eliminate labor-intensive chemical or enzymatic treatments has received considerable attention in recent years. Herein, we report a DNA methylation detection procedure based on fluorescence turn-on strategy. Target sequence was selected from Sept9 promoter region that has been reported as one of the most frequently methylated sites in colorectal cancer. Probe DNA was designed to be complementary to this sequence with an additional six cytosines in the middle to form an internal loop to host silver nanoclusters. The fluorescence intensity of the synthesized silver nanoclusters with the duplexes of probe-non-methylated target was significantly different from that of probe-methylated target. The fluorescence enhanced with increasing the methylated DNA concentration with a linear relation in the range of 1.0 × 10^-8 M to 5.0 × 10^-7 M with the detection limit of 8.2 × 10^-9 M, and quenched with non-methylated ones. The method was very specific in the presence of non-complementary sequences with maximum similarity of 40%. Circular dichroism spectra indicated that silver ions significantly affected the structure of methylated and non-methylated DNA into different extents which could further influence the nanocluster fluorescence. Finally, a method was introduced to meet the concerns in the applicability of the proposed method in real situation.