Proximal Gene Research Articles

Discovery of genes that positively affect biomass and stress associated traits in poplar

Woody biomass serves as a renewable resource for various industries, including pulp and paper production, construction, biofuels, and electricity generation. However, the molecular mechanisms behind biomass traits are poorly understood, which significantly curtails the speed and efficiency of their improvement. We used activation tagging to discover genes that can positively affect tree biomass-associated traits. We generated and screened under greenhouse conditions a population of 2,700 independent activation tagging lines. A total of 761 lines, which had significantly and positively affected at least one biomass-associated trait, were discovered. The tag was positioned in the genome for forty lines which were affected in multiple traits and activation of proximal genes validated for a subset. For two lines we fully recapitulated the phenotype of the original lines through overexpression. Moreover, the overexpression led to more pronounced and additional improvements, not observed in the original lines. Importantly, the overexpression of a Fasciclin-like gene (PtaFLA10) and a Patatin-like gene (PtaPAT) was found to substantially improve biomass, with a 40% increase in dry-stem weight, and enhance drought tolerance, respectively. Additionally, PtaPAT overexpression increased cellulose content, which is crucial for biofuel production. Our work shows that the activation tagging approach applied even on a non-genome saturation scale in a poplar tree can be successfully used for the discovery of genes positively modify biomass productivity. Such dominant forward genetics approaches can aid in biotechnological manipulation of woody biomass traits and help unravel the functions and mechanisms of individual genes, gene families, and regulatory modules.

Sustained alterations in proximal tubule gene expression in primary culture associate with HNF4A loss

Primary cultures of proximal tubule cells are widely used to model the behavior of kidney epithelial cells in vitro. However, de-differentiation of primary cells upon culture has been observed and appreciated for decades, yet the mechanisms driving this phenomenon remain poorly understood. This confounds the interpretation of experiments using primary kidney epithelial cells and prevents their use to engineer functional kidney tissue ex vivo. In this report, we measure the dynamics of cell-state transformations in early primary culture of mouse proximal tubules to identify key pathways and processes that correlate with and may drive de-differentiation. Our data show that the loss of proximal-tubule-specific genes is rapid, uniform, and sustained even after confluent, polarized epithelial monolayers develop. This de-differentiation occurs uniformly across many common culture condition variations. Changes in early culture were strongly associated with the loss of HNF4A. Exogenous re-expression of HNF4A can promote expression of a subset of proximal tubule genes in a de-differentiated proximal tubule cell line. Using genetically labeled proximal tubule cells, we show that selective pressures very early in culture influence which cells grow to confluence. Together, these data indicate that the loss of in vivo function in proximal tubule cultures occurs very early and suggest that the sustained loss of HNF4A is a key regulatory event mediating this change.

Squash vein yellowing virus from California emerged in the Middle East via intragenic and intergeneric recombination events in the hypervariable potyvirus P1 and ipomovirus P1a genes.

We present the complete sequence of the genomic RNA of an isolate of squash vein yellowing virus (Ipomovirus cucurbitavenaflavi) from California (SqVYV-CA) and show it is a recombinant virus with a highly divergent 5' UTR and proximal P1a gene. The evolution of SqVYV-CA involved an intrageneric event between unknown potyviruses, related to isolates of papaya ringspot virus (Potyvirus papayanuli) from the Old World, and an intergeneric event between this recombinant potyvirus (minor parent) and an isolate of SqVYV from Israel (SqVYV-IL) (major parent). These events occurred in mixed infections and in the potyvirus P1 and ipomovirus P1a recombination hotspots and resulted in SqVYV-CA having a potyvirus 5' UTR and chimeric P1-P1a gene/protein and the remainder of the genome from SqVYV-IL. The SqVYV-CA chimeric P1-P1a gene is under positive selection, and the protein is intrinsically disordered and may localize to the nucleus via nuclear localization signals in the P1 part. The C-terminal SqVYV-IL P1a part also diverged but retained the conserved serine protease motif. Furthermore, substantial divergence in SqVYV isolates from the Middle East was associated with genetic drift and a long evolutionary history in this region. The finding that the host range and symptomatology in cucurbits of SqVYV-CA is similar to those of SqVYV from Florida and SqVYV-IL, indicated that the recombinant part of the genome had no obvious effect on the virus-host interaction. A divergent part of the P1 sequence of the SqVYV-CA P1-P1a gene was used to develop a primer pair and RT-PCR test for specific detection of SqVYV-CA. This test was used to detect spread of SqVYV-CA to a new production area of California in 2021 and 2022. Together, these results demonstrate (i) a high level of genetic diversity exists among isolates of SqVYV and involved intra- and intergeneric recombination and genetic drift (mutation), (ii) evidence that SqVYV originated in the Middle East and that there were independent introductions into the New World and (iii) the remarkable genetic flexibility of the 5' proximal genes of these viruses.

Vegetable phylloplane microbiomes harbour class 1 integrons in novel bacterial hosts and drive the spread of chlorite resistance

Bacterial hosts in vegetable phylloplanes carry mobile genetic elements, such as plasmids and transposons that are associated with integrons. These mobile genetic elements and their cargo genes can enter human microbiomes via consumption of fresh agricultural produce, including uncooked vegetables. This presents a risk of acquiring antimicrobial resistance genes from uncooked vegetables. To better understand horizontal gene transfer of class 1 integrons in these compartments, we applied epicPCR, a single-cell fusion-PCR surveillance technique, to link the class 1 integron integrase (intI1) gene with phylogenetic markers of their bacterial hosts. Ready-to-eat salads carried class 1 integrons from the phyla Bacteroidota and Pseudomonadota, including four novel genera that were previously not known to be associated with intI1. We whole-genome sequenced Pseudomonas and Erwinia hosts of pre-clinical class 1 integrons that are embedded in Tn402-like transposons. The proximal gene cassette in these integrons was identified as a chlorite dismutase gene cassette, which we showed experimentally to confer chlorite resistance. Chlorine-derived compounds such as acidified sodium chlorite and chloride dioxide are used to disinfectant raw vegetables in food processing facilities, suggesting selection for chlorite resistance in phylloplane integrons. The spread of integrons conferring chlorite resistance has the potential to exacerbate integron-mediated antimicrobial resistance (AMR) via co-selection of chlorite resistance and AMR, thus highlighting the importance of monitoring chlorite residues in agricultural produce. These results demonstrate the strength of combining epicPCR and culture-based isolation approaches for identifying hosts and dissecting the molecular ecology of class 1 integrons.

Elucidating the Subcellular Localization of GLRaV-3 Proteins Encoded by the Unique Gene Block in N. benthamiana Suggests Implications on Plant Host Suppression.

Grapevine leafroll-associated virus 3 (GLRaV-3) is a formidable threat to the stability of the global grape and wine industries. It is the primary etiological agent of grapevine leafroll disease (GLD) and significantly impairs vine health, fruit quality, and yield. GLRaV-3 is a member of the genus Ampelovirus, Closteroviridae family. Viral genes within the 3' proximal unique gene blocks (UGB) remain highly variable and poorly understood. The UGBs of Closteroviridae viruses include diverse open reading frames (ORFs) that have been shown to contribute to viral functions such as the suppression of the host RNA silencing defense response and systemic viral spread. This study investigates the role of GLRaV-3 ORF8, ORF9, and ORF10, which encode the proteins p21, p20A, and p20B, respectively. These genes represent largely unexplored facets of the GLRaV-3 genome. Here, we visualize the subcellular localization of wildtype and mutagenized GLRaV-3 ORFs 8, 9, and 10, transiently expressed in Nicotiana benthamiana. Our results indicate that p21 localizes to the cytosol, p20A associates with microtubules, and p20B is trafficked into the nucleus to carry out the suppression of host RNA silencing. The findings presented herein provide a foundation for future research aimed at the characterization of the functions of these ORFs. In the long run, it would also facilitate the development of innovative strategies to understand GLRaV-3, mitigate its spread, and impacts on grapevines and the global wine industry.

Analyzing Medicago spp. seed morphology using GWAS and machine learning

Alfalfa is widely recognized as an important forage crop. To understand the morphological characteristics and genetic basis of seed morphology in alfalfa, we screened 318 Medicago spp., including 244 Medicago sativa subsp. sativa (alfalfa) and 23 other Medicago spp., for seed area size, length, width, length-to-width ratio, perimeter, circularity, the distance between the intersection of length & width (IS) and center of gravity (CG), and seed darkness & red–green–blue (RGB) intensities. The results revealed phenotypic diversity and correlations among the tested accessions. Based on the phenotypic data of M. sativa subsp. sativa, a genome-wide association study (GWAS) was conducted using single nucleotide polymorphisms (SNPs) called against the Medicago truncatula genome. Genes in proximity to associated markers were detected, including CPR1, MON1, a PPR protein, and Wun1(threshold of 1E−04). Machine learning models were utilized to validate GWAS, and identify additional marker-trait associations for potentially complex traits. Marker S7_33375673, upstream of Wun1, was the most important predictor variable for red color intensity and highly important for brightness. Fifty-two markers were identified in coding regions. Along with strong correlations observed between seed morphology traits, these genes will facilitate the process of understanding the genetic basis of seed morphology in Medicago spp.

A genetic-epigenetic interplay at 1q21.1 locus underlies CHD1L-mediated vulnerability to primary progressive multiple sclerosis

Multiple Sclerosis (MS) is a heterogeneous inflammatory and neurodegenerative disease with an unpredictable course towards progressive disability. Treating progressive MS is challenging due to limited insights into the underlying mechanisms. We examined the molecular changes associated with primary progressive MS (PPMS) using a cross-tissue (blood and post-mortem brain) and multilayered data (genetic, epigenetic, transcriptomic) from independent cohorts. In PPMS, we found hypermethylation of the 1q21.1 locus, controlled by PPMS-specific genetic variations and influencing the expression of proximal genes (CHD1L, PRKAB2) in the brain. Evidence from reporter assay and CRISPR/dCas9 experiments supports a causal link between methylation and expression and correlation network analysis further implicates these genes in PPMS brain processes. Knock-down of CHD1L in human iPSC-derived neurons and knock-out of chd1l in zebrafish led to developmental and functional deficits of neurons. Thus, several lines of evidence suggest a distinct genetic-epigenetic-transcriptional interplay in the 1q21.1 locus potentially contributing to PPMS pathogenesis.

Efficient activation of hundreds of LTR12C elements reveals cis-regulatory function determined by distinct epigenetic mechanisms.

Long terminal repeats (LTRs), which often contain promoter and enhancer sequences of intact endogenous retroviruses (ERVs), are known to be co-opted as cis-regulatory elements for fine-tuning host-coding gene expression. Since LTRs are mainly silenced by the deposition of repressive epigenetic marks, substantial activation of LTRs has been found in human cells after treatment with epigenetic inhibitors. Although the LTR12C family makes up the majority of ERVs activated by epigenetic inhibitors, how these epigenetically and transcriptionally activated LTR12C elements can regulate the host-coding gene expression remains unclear due to genome-wide alteration of transcriptional changes after epigenetic inhibitor treatments. Here, we specifically transactivated >600 LTR12C elements by using single guide RNA-based dCas9-SunTag-VP64, a site-specific targeting CRISPR activation (CRISPRa)system, with minimal off-target events. Interestingly, most of the transactivated LTR12C elements acquired the H3K27ac-marked enhancer feature, while only 20% were co-marked with promoter-associated H3K4me3 modifications. The enrichment of the H3K4me3 signal was intricately associated with downstream regions of LTR12C, such as internal regions of intact ERV9 or other types of retrotransposons. Here, we leverage an optimized CRISPRa system to identify two distinct epigenetic signatures that define LTR12C transcriptional activation, which modulate the expression of proximal protein-coding genes.

Genomic dissection and mutation-specific target discovery for breast cancer PIK3CA hotspot mutations

BackgroundRecent advancements in high-throughput genomics and targeted therapies have provided tremendous potential to identify and therapeutically target distinct mutations associated with cancers. However, to date the majority of targeted therapies are used to treat all functional mutations within the same gene, regardless of affected codon or phenotype.ResultsIn this study, we developed a functional genomic analysis workflow with a unique isogenic cell line panel bearing two distinct hotspot PIK3CA mutations, E545K and H1047R, to accurately identify targetable differences between mutations within the same gene. We performed RNA-seq and ATAC-seq and identified distinct transcriptomic and epigenomic differences associated with each PIK3CA hotspot mutation. We used this data to curate a select CRISPR knock out screen to identify mutation-specific gene pathway vulnerabilities. These data revealed AREG as a E545K-preferential target that was further validated through in vitro analysis and publicly available patient databases.ConclusionsUsing our multi-modal genomics framework, we discover distinct differences in genomic regulation between PIK3CA hotspot mutations, suggesting the PIK3CA mutations have different regulatory effects on the function and downstream signaling of the PI3K complex. Our results demonstrate the potential to rapidly uncover mutation specific molecular targets, specifically AREG and a proximal gene regulatory region, that may provide clinically relevant therapeutic targets. The methods outlined provide investigators with an integrative strategy to identify mutation-specific targets for the treatment of other oncogenic mutations in an isogenic system.

Genome-wide association analyses using multilocus models on bananas (Musa spp.) reveal candidate genes related to morphology, fruit quality, and yield.

Bananas (Musa spp.) are an essential fruit worldwide and rank as the fourth most significant food crop for addressing malnutrition due to their rich nutrients and starch content. The potential of their genetic diversity remains untapped due to limited molecular breeding tools. Our study examined a phenotypically diverse group of 124 accessions from the Colombian Musaceae Collection conserved in AGROSAVIA. We assessed 12 traits categorized into morphology, fruit quality, and yield, alongside sequence data. Our sequencing efforts provided valuable insights, with an average depth of about 7× per accession, resulting in 187,133 single-nucleotide polymorphisms (SNPs) against Musa acuminata (A genome) and 220,451 against Musa balbisiana (B genome). Population structure analysis grouped samples into four and five clusters based on the reference genome. By using different association models, we identified marker-trait associations (MTAs). The mixed linear model revealed four MTAs, while the Bayesian-information and linkage-disequilibrium iteratively nested keyway and fixed and random model for circulating probability unification models identified 82 and 70 MTAs, respectively. We identified 38 and 40 candidate genes in linkage proximity to significant MTAs for the A genome and B genome, respectively. Our findings provide insights into the genetic underpinnings of morphology, fruit quality, and yield. Once validated, the SNP markers and candidate genes can potentially drive advancements in genomic-guided breeding strategies to enhance banana crop improvement.

Chimera RNA transcribed from integrated HPV18 genome with adjacent host genomic region promotes oncogenic gene expression through condensate formation.

Most cervical cancers are caused by human papillomavirus (HPV) infection. In HeLa cells, the HPV18 viral genome is integrated at chromosome 8q24.21 and activates transcription of the proto-oncogene c-Myc. However, the mechanism of how the integrated HPV genome and its transcribed RNAs exhibit transcription activation function has not been fully elucidated. In this study, we found that HPV18 transcripts contain an enhancer RNA-like function to activate proximal genes including CCAT1-5L and c-Myc. We showed that the human genome-integrated HPV18 genes are activated by transcription coregulators including BRD4 and Mediator. The transcribed HPV18 RNAs form a liquid-like condensate at chromosome 8q24.21 locus, which in turn accumulates RNA polymerase II. Moreover, we focused on a relatively uncharacterized transcript from the upstream region of CCAT1, named URC. The URC RNA is transcribed as a chimera RNA with HPV18 and is composed of the 3'-untranslated region of the HPV18 transcript. We experimentally showed that the URC contributes to stabilization of HPV18 RNAs by supplying a polyadenylation site for the HPV18 transcript. Our findings suggest that integrated HPV18 at 8q24.21 locus produces HPV18-URC chimera RNA and promotes tumorigenesis through RNA-based condensate formation.

Comparative analysis of sipeimine content, metabolome and chloroplast genome in cultivated and wild varieties of Fritillaria taipaiensis.

The wild variety Fritillaria taipaiensis E.B (EB) is known for its superior therapeutic effects, but its limited production cannot meet demand. As a result, the cultivated variety F. taipaiensis P. Y. Li (PY) has been widely grown. In this study, we conducted a comprehensive analysis comparing EB and PY in terms of external features, sipeimine content, metabolome and chloroplast genome to differentiate these two varieties. Our research revealed that the petals and pods of EB are green, while those of PY have purple markings. The bulbs of EB contain significantly higher levels of sipeimine compared to those of PY. Metabolomic analysis identified 56 differentially expressed metabolites (DMs), with 23 upregulated and 33 downregulated in EB bulbs. Particularly, 3-hydroxycinnamic acid and secoxyloganin may serve as distinctive DMs. These DMs were associated with 17 KEGG pathways, including pyrimidine metabolism, alanine, aspartate and glutamate metabolism, and galactose metabolism. Differences in the length of the chloroplast genome were primarily observed in the large single-copy (LSC) region, with the largest variation in the trnH-GUC-psbA region. The placement of the trnH gene and the rps gene in proximity to the LSC/IRb boundary differs between EB and PY. The results of this study provide valuable insights for the introduction and comprehensive development of wild F. taipaiensis from a scientific perspective. © 2024 Society of Chemical Industry.

Abstract 1049: Long Noncoding RNA Regulation Of Inflammation In Atherosclerosis

Atherosclerosis is a chronic inflammatory disease driven by cholesterol deposition in the artery wall. Resolution of atherosclerotic inflammation and plaque regression are important clinical goals. As long noncoding RNAs (lncRNAs) are key regulators of gene expression, there is an intense interest in understanding their roles in atherosclerosis and their potential as therapeutic targets. To identify lncRNAs regulated during plaque progression and regression, we performed single cell RNA-sequencing of CD45 + immune cells isolated from aortic arch plaques of Ldlr –/– mice under conditions of plaque progression (16 wks Western diet; WD) or regression (16 wks WD + 4 wks apoB ASO). We screened for lncRNAs acting in cis to modify the expression of proximal protein-coding genes and identified Caslr ( C tsc a nti s ense l nc R NA) , a novel lncRNA transcribed from the opposite strand of Cathepsin C gene ( Ctsc ) encoding a pro-atherosclerotic lysosomal protease. Caslr–Ctsc pair is upregulated in macrophages of progressing plaques and in response to pro-inflammatory stimuli ( i.e . LPS, IL-6) in vitro , and conversely downregulated during atherosclerosis regression and in response to anti-inflammatory stimulus ( i.e. IL-4). LncRNA knockdown approaches combined with transcriptomic analyses and functional assays showed that Caslr regulates Ctsc expression and associated inflammatory responses, including NF-κB signaling and inflammasome-dependent IL-1β release. RNA purification strategies combined with mass-spectrometry revealed that Caslr interacts with spliceosome-associated proteins to coordinate splicing of Ctsc mRNA. Notably, Caslr shares functional homology with a human lncRNA, CASLR, proximal to the CTSC gene, which regulates CTSC expression and protease activity. Human CASLR and CTSC are highly expressed in the blood of hyperlipidemic patients and their expression is coordinately decreased upon aggressive cholesterol lowering by anti-PCSK9 treatment. Collectively, our data identify Caslr as an evolutionary conserved lncRNA acting in cis to drive the expression of cathepsin C with downstream consequences on inflammatory signaling in plaque macrophages, and suggest that targeting this lncRNA may be beneficial in treating atherosclerosis.

Abstract 151: Single Cell Transcriptional And Epigenetic Map Of Human Coronary Artery Disease In A Multi-ethnic Population

Identifying target genes that are causally linked to coronary artery disease (CAD) remains an elusive problem despite the plethora of CAD genome-wide association study (GWAS) risk loci. Here, we performed single nucleus multiomic (RNA + ATAC) sequencing in 44 human coronary arteries and bulk Hi-C in 16 coronaries from multi-ethnic participants. From sequencing 126,804 nuclei, we identified cell type specific transcriptional signatures, chromatin peaks, and gene regulatory networks through peak to gene linkage analysis. Single-cell quantitative trait loci (QTL) analysis identified cell-specific caQTL’s including cell-type specific caQTLs within CAD GWAS loci. The overlap of CAD GWAS loci with ATAC peaks demonstrated that vascular smooth muscle and endothelial cells harbor the greatest genetic risk. To identify proximal and distal putative target genes for CAD loci in a cell type specific manner, we built a enhancer-gene map by leveraging coronary artery ChiP-seq, cell specific RNA/ATAC-seq, and Hi-C based chromatin contact loops. Using this map, we identified distal target genes linked to CAD GWAS loci. Finally, given the effect of ancestry on CAD GWAS risk interpretation, we leveraged Hi-C data to compare the chromatin architecture in coronary tissue from 8 European Americans vs 8 African Americans. A meta-map, sample level, and local ancestry analysis identified ancestry-specific differences in chromatin looping and topologically associated domains. Collectively, our results provide a cell-type specific transcriptomic and epigenetic map which along with bulk chromatin contacts can be used to interpret and prioritize CAD GWAS candidates and nominate cell specific causal target genes.

GLI1 Coamplification in Well-Differentiated/Dedifferentiated Liposarcomas: Clinicopathologic and Molecular Analysis of 92 Cases

GLI1(12q13.3) amplification is identified in a subset of mesenchymal neoplasms with a distinct nested round cell/epithelioid phenotype. MDM2 and CDK4 genes are situated along the oncogenic 12q13-15 segment, amplification of which defines well-differentiated liposarcoma (WDLPS)/dedifferentiated liposarcoma (DDLPS). The 12q amplicon can occasionally include GLI1, a gene in close proximity to CDK4. We hereby describe the first cohort of GLI1/MDM2/CDK4 coamplified WD/DDLPS. The departmental database was queried retrospectively for all cases of WD/DDLPS having undergone next-generation (MSK-IMPACT) sequencing with confirmed MDM2, CDK4, and GLI1 coamplification. Clinicopathologic data was obtained from a review of the medical chart and available histologic material. Four hundred eighty-six WD/DDLPS cases underwent DNA sequencing, 92 (19%) of which harbored amplification of the GLI1 locus in addition to that of MDM2 and CDK4. These included primary tumors (n = 60), local recurrences (n = 29), and metastases (n = 3). Primary tumors were most frequently retroperitoneal (47/60, 78%), mediastinal (4/60, 7%), and paratesticular (3/60, 5%). Average age was 63 years, with a male:female ratio of 3:2. The cohort was comprised of DDLPS (86/92 [93%], 6 of which were WDLPS with early dedifferentiation) and WDLPS without any longitudinal evidence of dedifferentiation (6/92, 7%). One-fifth (13/86, 17%) of DDLPS cases showed no evidence of a well-differentiated component in any of the primary, recurrent, or metastatic specimens. Dedifferentiated areas mostly showed high-grade undifferentiated pleomorphic sarcoma-like (26/86,30%) and high-grade myxofibrosarcoma-like (13/86,16%) morphologies. A disproportionately increased incidence of meningothelial whorls with/without osseous metaplasia was observed as the predominant pattern in 16/86 (19%) cases, and GLI1-altered morphology as described was identified in a total of 10/86 (12%) tumors. JUN (1p32.1), also implicated in the pathogenesis of WD/DDLPS, was coamplified with all 3 of MDM2, CDK4, and GLI1 in 7/91 (8%) cases. Additional loci along chromosomal arms 1p and 6q, including TNFAIP3, LATS1, and ESR1, were also amplified in a subset of cases. In this large-scale cohort of GLI1 coamplified WD/DDLPS, we elucidate uniquely recurrent features including meningothelial whorl-like and GLI-altered morphology in dedifferentiated areas. Assessment of tumor location (retroperitoneal or mediastinal), identification of a well-differentiated liposarcoma component, and coamplification of other spatially discrete genomic segments (1p and 6q) might aid in distinction from tumors with true driver GLI1 alterations.

A computational approach for deciphering the interactions between proximal and distal gene regulators in GC B-cell response.

Delineating the intricate interplay between promoter-proximal and -distal regulators is crucial for understanding the function of transcriptional mediator complexes implicated in the regulation of gene expression. The present study aimed to develop a computational method for accurately modeling the spatial proximal and distal regulatory interactions. Our method combined regression-based models to identify key regulators through gene expression prediction and a graph-embedding approach to detect coregulated genes. This approach enabled a detailed investigation of the gene regulatory mechanisms for germinal center B cells, accompanied by dramatic rearrangements of the genome structure. We found that while the promoter-proximal regulatory elements were the principal regulators of gene expression, the distal regulators fine-tuned transcription. Moreover, our approach unveiled the presence of modular regulators, such as cofactors and proximal/distal transcription factors, which were co-expressed with their target genes. Some of these modules exhibited abnormal expression patterns in lymphoma. These findings suggest that the dysregulation of interactions between transcriptional and architectural factors is associated with chromatin reorganization failure, which may increase the risk of malignancy. Therefore,our computational approachhelps decipherthe transcriptional cis-regulatory code spatially interacting.

The Immunometabolic Gene N-Acetylglucosamine Kinase Is Uniquely Involved in the Heritability of Multiple Sclerosis Severity.

The clinical severity of multiple sclerosis (MS), an autoimmune disorder of the central nervous system, is thought to be determined by environmental and genetic factors that have not yet been identified. In a recent genome-wide association study (GWAS), a single nucleotide polymorphism (SNP), rs10191329, has been associated with MS severity in two large independent cohorts of patients. Different approaches were followed by the authors to prioritize the genes that are transcriptionally regulated by such an SNP. It was concluded that the identified SNP regulates a group of proximal genes involved in brain resilience and cognitive abilities rather than immunity. Here, by conducting an alternative strategy for gene prioritization, we reached the opposite conclusion. According to our re-analysis, the main target of rs10191329 is N-Acetylglucosamine Kinase (NAGK), a metabolic gene recently shown to exert major immune functions via the regulation of the nucleotide-binding oligomerization domain-containing protein 2 (NOD2) pathway. To gain more insights into the immunometabolic functions of NAGK, we analyzed the currently known list of NAGK protein partners. We observed that NAGK integrates a dense network of human proteins that are involved in glucose metabolism and are highly expressed by classical monocytes. Our findings hold potentially major implications for the understanding of MS pathophysiology.

A case report of a family with developmental arrest of human prokaryotic stage zygote.

To study the genetic variation leading to the arrest phenotype of pronuclear (PN) zygotes. We recruited a family characterized by recurrent PN arrest during in vitro fertilization (IVF) and intracytoplasmic sperm injection cycles (ICSI) and performed whole-exome sequencing for 2 individuals. The transcriptome profiles of PN-arrest zygotes were assessed by single-cell RNA sequencing analysis. The variants were then validated by PCR amplification and Sanger sequencing in the affected individuals and other family members. A family characterized by recurrent PN arrest during IVF and ICSI cycles were enrolled after giving written informed consent. Peripheral blood samples were taken for DNA extraction. Three PN-arrest zygotes from patient III-3 were used for single-cell RNA-seq as described. This phenotype was reproduced after multiple cycles of egg retrieval and after trying different fertilization methods and multiple ovulation regimens. The mutant genes of whole exon sequencing were screened and verified. The missense variant c. C1630T (p.R544W) in RGS12 was responsible for a phenotype characterized by paternal transmission. RGS12 controls Ca2+ oscillation, which is required for oocyte activation after fertilization. Single-cell transcriptome profiling of PN-arrest zygotes revealed defective established translation, RNA processing and cell cycle, which explained the failure of complete oocyte activation. Furthermore, we identified proximal genes involved in Ca2+ oscillation-cytostatic factor-anaphase-promoting complex (Ca2+ oscillation-CSF-APC) signaling, including upregulated CaMKII, ORAI1, CDC20, and CDH1 and downregulated EMI1 and BUB3. The findings indicate abnormal spontaneous Ca2+ oscillations leading to oocytes with prolonged low CSF level and high APC level, which resulted in defective nuclear envelope breakdown and DNA replication. We have identified an RGS12 variant as the potential cause of female infertility characterized by arrest at the PN stage during multiple IVF and ICSI.

Abstract 6944: Transcriptional programs underpinning reoccurring cancer cell states, plasticity and tumour phylogenies in clear cell kidney cancer

Abstract Intratumor heterogeneity (ITH) is a substrate for tumor evolution and confers amenability to stressors and therapies. Sources of ITH include subclonal mutations, microenvironmental factors, and plasticity whereby a cell has the ability to adopt different cell states. Clear cell kidney cancer (ccRCC) has a well-characterized genetic initiating event - VHL inactivation - resulting in activation of HIF-1 and HIF-2 transcription factors. Although VHL inactivation is consistently clonal in nature, ccRCC nonetheless exhibits extensive genetic and morphological ITH. This makes ccRCC an excellent model in which to study ITH. We performed single cell RNA-seq on multiregional samples from kidney tumors and macroscopically normal kidney tissue. We obtained 55 samples from 16 tumors and used the 10x Genomics platform to obtain gene expression data from ~150,000 single cells, including ~80,000 cancer and normal epithelial cells. We identified three cancer cell states, co-existing in all clear cell kidney cancer tumors analyzed and marked by different gene expression programs: 1) differentiated-like cells marked by proximal tubule cell genes characteristic of the cancer cell of origin; 2) epithelial-mesenchymal transition (EMT)-like cells and 3) injured-like cells. Differentiated markers and EMT markers were associated with good and poor prognosis, respectively, indicating that these cell subpopulations are clinically important. Notably, these cell states are also present in the normal proximal tubule although in different proportions to the tumor samples, which contain more cells in the EMT-like and injured-like states. This suggests that the core features of ITH are inherent to the cell of origin of the cancer, but that the cell state trajectories are dysregulated in the tumor. Furthermore, in tumors these states exist as a continuum, whereby transitional cells in between states appear to exist, suggesting that cancer cells are continually and dynamically transitioning between states. Lastly, we used chromosomal copy number variation (CNV) patterns to study genetic ITH. We examined the effect of 14q loss, which is associated with ccRCC metastasis, on gene expression. We found that 14q loss is associated with downregulation of HIF-1 transcriptional activity together with a reciprocal increase in HIF-2 activity, MYC activity and oxidative phosphorylation, highlighting potential mechanisms by which 14q loss may driver ccRCC metastasis. Citation Format: Olivia Lombardi, Ran Li, Lisa Browning, Faiz Jabbar, Hannah Evans, Peter J. Ratcliffe, David R. Mole. Transcriptional programs underpinning reoccurring cancer cell states, plasticity and tumour phylogenies in clear cell kidney cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 6944.

Functional identification of cis-regulatory long noncoding RNAs at controlled false discovery rates.

A key attribute of some long noncoding RNAs (lncRNAs) is their ability to regulate expression of neighbouring genes in cis. However, such 'cis-lncRNAs' are presently defined using ad hoc criteria that, we show, are prone to false-positive predictions. The resulting lack of cis-lncRNA catalogues hinders our understanding of their extent, characteristics and mechanisms. Here, we introduce TransCistor, a framework for defining and identifying cis-lncRNAs based on enrichment of targets amongst proximal genes. TransCistor's simple and conservative statistical models are compatible with functionally defined target gene maps generated by existing and future technologies. Using transcriptome-wide perturbation experiments for 268 human and 134 mouse lncRNAs, we provide the first large-scale survey of cis-lncRNAs. Known cis-lncRNAs are correctly identified, including XIST, LINC00240 and UMLILO, and predictions are consistent across analysis methods, perturbation types and independent experiments. We detectcis-activityin a minority of lncRNAs, primarily involving activators over repressors. Cis-lncRNAs are detected by both RNA interference and antisense oligonucleotide perturbations. Mechanistically, cis-lncRNA transcripts are observed to physically associate with their target genes and are weakly enriched with enhancer elements. In summary, TransCistor establishes a quantitative foundation for cis-lncRNAs, opening a path to elucidating their molecular mechanisms and biological significance.