Genome-wide Target Research Articles

HASCH - A high-throughput amplicon-based SNP-platform for medicinal cannabis and industrial hemp genotyping applications

BackgroundCannabis sativa is seeing a global resurgence as a food, fiber and medicinal crop for industrial hemp and medicinal Cannabis industries respectively. However, a widespread moratorium on the use and research of C. sativa throughout most of the 20th century has seen the development of improved cultivars for specific end uses lag behind that of conventional crops. While C. sativa research and development has seen significant investments in the recent past, resulting in a suite of publicly available genomic resources and tools, a versatile and cost-effective mid-density genotyping platform for applied purposes in breeding and pre-breeding is lacking. Here we report on a first mid-density fixed-target SNP platform for C. sativa.ResultsThe High-throughput Amplicon-based SNP-platform for medicinal Cannabis and industrial Hemp (HASCH) was designed using a combination of filtering and Integer Linear Programming on publicly available whole-genome sequencing and RNA sequencing data, supplemented with in-house generated genotyping-by-sequencing (GBS) data. HASCH contains 1,504 genome-wide targets of high call rate (97% mean) and even distribution across the genome, designed to be highly informative (> 0.3 minor allele frequency) across both medicinal cannabis and industrial hemp gene pools. Average numbers of mismatch SNP between any two accessions were 251 for medicinal cannabis (N = 116) and 272 for industrial hemp (N = 87). Comparing HASCH data with corresponding GBS data on a collection of diverse C. sativa accessions demonstrated high concordance and resulted in comparable phylogenies and genetic distance matrices. Using HASCH on a segregating F2 population derived from a cross between a tetrahydrocannabinol (THC)-dominant and a cannabidiol (CBD)-dominant accession resulted in a genetic map consisting of 310 markers, comprising 10 linkage groups and a total size of 582.7 cM. Quantitative Trait Locus (QTL) mapping identified a major QTL for CBD content on chromosome 7, consistent with previous findings.ConclusionHASCH constitutes a versatile, easy to use and cost-effective genotyping solution for the rapidly growing Cannabis research community. It provides consistent genetic fingerprints of 1504 SNPs with wide applicability genetic resource management, quantitative genetics and breeding.

CRISPR-Cas9 and Cas12a target site richness reflects genomic diversity in natural populations of Anopheles gambiae and Aedes aegypti mosquitoes

Due to limitations in conventional disease vector control strategies including the rise of insecticide resistance in natural populations of mosquitoes, genetic control strategies using CRISPR gene drive systems have been under serious consideration. The identification of CRISPR target sites in mosquito populations is a key aspect for developing efficient genetic vector control strategies. While genome-wide Cas9 target sites have been explored in mosquitoes, a precise evaluation of target sites focused on coding sequence (CDS) is lacking. Additionally, target site polymorphisms have not been characterized for other nucleases such as Cas12a, which require a different DNA recognition site (PAM) and would expand the accessibility of mosquito genomes for genetic engineering. We undertook a comprehensive analysis of potential target sites for both Cas9 and Cas12a nucleases within the genomes of natural populations of Anopheles gambiae and Aedes aegypti from multiple continents. We demonstrate that using two nucleases increases the number of targets per gene. Also, we identified differences in nucleotide diversity between North American and African Aedes populations, impacting the abundance of good target sites with a minimal degree of polymorphisms that can affect the binding of gRNA. Lastly, we screened for gRNAs targeting sex-determination genes that could be widely applicable for developing field genetic control strategies. Overall, this work highlights the utility of employing both Cas9 and Cas12a nucleases and underscores the importance of designing universal genetic strategies adaptable to diverse mosquito populations.

Polygenic risk score of metabolic dysfunction-associated steatotic liver disease amplifies the health impact on severe liver disease and metabolism-related outcomes

BackgroundAlthough the inherited risk factors associated with fatty liver disease are well understood, little is known about the genetic background of metabolic dysfunction-associated steatotic liver disease (MASLD) and its related health impacts. Compared to non-alcoholic fatty liver disease (NAFLD), MASLD presents significantly distinct diagnostic criteria, and epidemiological and clinical features, but the related genetic variants are yet to be investigated. Therefore, we conducted this study to assess the genetic background of MASLD and interactions between MASLD-related genetic variants and metabolism-related outcomes.MethodsParticipants from the UK Biobank were grouped into discovery and replication cohorts for an MASLD genome-wide association study (GWAS), and base and target cohorts for polygenic risk score (PRS) analysis. Autosomal genetic variants associated with NAFLD were compared with the MASLD GWAS results. Kaplan–Meier and Cox regression analyses were used to assess associations between MASLD and metabolism-related outcomes.ResultsSixteen single-nucleotide polymorphisms (SNPs) were identified at genome-wide significance levels for MASLD and duplicated in the replication cohort. Differences were found after comparing these SNPs with the results of NAFLD-related genetic variants. MASLD cases with high PRS had a multivariate-adjusted hazard ratio of 3.15 (95% confidence interval, 2.54–3.90) for severe liver disease (SLD), and 2.81 (2.60–3.03) for type 2 diabetes mellitus. The high PRS amplified the impact of MASLD on SLD and extrahepatic outcomes.ConclusionsHigh PRS of MASLD GWAS amplified the impact of MASLD on SLD and metabolism-related outcomes, thereby refining the process of identification of individuals at high risk of MASLD. Supplementation of this process with relevant genetic backgrounds may lead to more effective MASLD prevention and management.

Nuclear receptor Nur77 and Yin-Yang 1 synergistically increase mitochondrial abundance and activity in macrophages.

Mitochondrial biogenesis requires precise regulation of both mitochondrial-encoded and nuclear-encoded genes. Nuclear receptor Nur77 is known to regulate mitochondrial metabolism in macrophages and skeletal muscle. Here, we compared genome-wide Nur77 binding site and target gene expression in these two cell types, which revealed conserved regulation of mitochondrial genes and enrichment of motifs for the transcription factor Yin-Yang 1 (YY1). We show that Nur77 and YY1 interact, that YY1 increases Nur77 activity, and that their binding sites are co-enriched at mitochondrial ribosomal protein gene loci in macrophages. Nur77 and YY1 co-expression synergistically increases Mrpl1 expression as well as mitochondrial abundance and activity in macrophages but not skeletal muscle. As such, we identify a macrophage-specific Nur77-YY1 interaction that enhances mitochondrial metabolism.

Abstract 5636: CUX1 transcriptionally regulates intestinal epithelial homeostasis

Abstract The small intestine digests and absorbs nutrients, a role made possible by specialized secretory and absorptive cells that constantly replenish from the intestinal stem cell niche. Disruption in intestinal epithelium homeostasis has been linked to pathologies including inflammatory bowel disease (IBD) and cancer. However, the processes regulating differentiation in healthy and disease setting are incompletely understood. CUX1 is a highly conserved homeodomain-containing transcription factor that is expressed in small intestinal crypts and has been previously shown to be essential in mice and Drosophila. CUX1 is upregulated upon inflammatory stress conditions, including intestinal damage and in IBD. In addition, CUX1 is predicted to be protective against IBD from genome-wide association studies. The aim of the current work is to determine the mechanistic role of CUX1 in secretory intestinal cells differentiation and epithelial homeostasis. To this end, we have generated transgenic mouse models that modulate CUX1 levels in an inducible manner. Mice with near-null CUX1 levels exhibit rapid weight loss requiring euthanasia while mice with ~50% residual CUX1 halt normal weight gain. Histologic examination of tissues from mice with CUX1 knock-down reveals a lack of mature secretory cells, including goblet and Paneth cells, the latter playing a supportive role to intestinal stem cells. Moreover, mouse intestinal tissues from our knock-down models demonstrated decreased proliferation with increased apoptosis. In transcriptome profiling, we observed downregulation of proliferation, intestinal stem cell, and digestive system developmental gene signatures and decreased levels of key secretory cell regulators, pointing towards a role for CUX1 in stem cell maintenance and differentiation. Indeed, reducing CUX1 levels decreased β-catenin protein levels, while Wnt/β-catenin signaling is required for the organization of the crypt-villus axis. Finally, CUX1 knock-down compromises the ability to form budding organoids, demonstrating the cell intrinsic role for CUX1 in stem cell differentiation. In ongoing studies, we are identifying CUX1 genome-wide binding targets and the epigenetic impact of CUX1 in intestinal epithelium cells. Our findings indicate that CUX1 is a critical regulator of intestinal epithelium homeostasis that has a role in differentiation and potential role in injury-induced de-differentiation in healthy as well as IBD patients. We propose that CUX1 regulates stem cell maintenance and crypt composition by regulating Wnt/β-catenin signaling and inducing fetal-like transcriptomic repair program shown to be coopted by colon cancers. Resolving an intricate network of transcriptional regulation of differentiation in the intestinal crypt holds the potential to successfully target inflammation-related damage in human intestines. Citation Format: Katarzyna Zawieracz, Ningfei An, Megan McNerney. CUX1 transcriptionally regulates intestinal epithelial homeostasis [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 5636.

In silico identification of papaya genome-encoded microRNAs to target begomovirus genes in papaya leaf curl disease.

Papaya leaf curl disease (PaLCuD) is widespread and classified in the genus begomovirus (Geminiviridae), disseminated by the vector whitefly Bemisia tabaci. RNA interference (RNAi)-based antiviral innate immunity stands as a pivotal defense mechanism and biological process in limiting viral genomes to manage plant diseases. The current study aims to identify and analyze Carica Papaya locus-derived capa-microRNAs with predicted potential for targeting divergent begomovirus species-encoded mRNAs using a 'four integrative in silico algorithms' approach. This research aims to experimentally activate the RNAi catalytic pathway using in silico-predicted endogenous capa-miRNAs and create papaya varieties capable of assessing potential resistance against begomovirus species and monitoring antiviral capabilities. This study identified 48 predicted papaya locus-derived candidates from 23 miRNA families, which were further investigated for targeting begomovirus genes. Premised all the four algorithms combined, capa-miR5021 was the most anticipated miRNA followed by capa-miR482, capa-miR5658, capa-miR530b, capa-miR3441.2, and capa-miR414 'effective' papaya locus-derived candidate capa-miRNA and respected putative binding sites for targets at the consensus nucleotide position. It was predicted to bind and target mostly to AC1 gene of the complementary strand and the AV1 gene of the virion strand of different begomovirus isolates, which were associated with replication-associated protein and encapsidation, respectively, during PaLCuD. These miRNAs were also found targeting betaC1 gene of betasatellite which were associated with retardation in leaf growth and developmental abnormalities with severe symptoms during begomovirus infection. To validate target prediction accuracy, we created an integrated Circos plot for comprehensive visualization of host-virus interaction. In silico-predicted papaya genome-wide miRNA-mediated begomovirus target gene regulatory network corroborated interactions that permit in vivo analysis, which could provide biological material and valuable evidence, leading to the development of begomovirus-resistant papaya plants. The integrative nature of our research positions it at the forefront of efforts to ensure the sustainable cultivation of papaya, particularly in the face of evolving pathogenic threats. As we move forward, the knowledge gained from this study provides a solid foundation for continued exploration and innovation in the field of papaya virology, and to the best of our knowledge, this study represents a groundbreaking endeavor, undertaken for the first time in the context of PaLCuD research.

Rewiring of the epigenome and chromatin architecture by exogenously induced retinoic acid signaling during zebrafish embryonic development.

Retinoic acid (RA) is the ligand of RA receptors (RARs), transcription factors that bind to RA response elements. RA signaling is required for multiple processes during embryonic development, including body axis extension, hindbrain antero-posterior patterning and forelimb bud initiation. Although some RA target genes have been identified, little is known about the genome-wide effects of RA signaling during in vivo embryonic development. Here, we stimulate the RA pathway by treating zebrafish embryos with all-trans-RA (atRA) and use a combination of RNA-seq, ATAC-seq, ChIP-seq and HiChIP to gain insight into the molecular mechanisms by which exogenously induced RA signaling controls gene expression. We find that RA signaling is involved in anterior/posterior patterning, central nervous system development, and the transition from pluripotency to differentiation. AtRA treatment also alters chromatin accessibility during early development and promotes chromatin binding of RARαa and the RA targets Hoxb1b, Meis2b and Sox3, which cooperate in central nervous system development. Finally, we show that exogenous RA induces a rewiring of chromatin architecture, with alterations in chromatin 3D interactions involving target genes. Altogether, our findings identify genome-wide targets of RA signaling and provide a molecular mechanism by which developmental signaling pathways regulate target gene expression by altering chromatin topology.

Genomic overview of INA-induced NPR1 targeting and transcriptional cascades in Arabidopsis.

The phytohormone salicylic acid (SA) triggers transcriptional reprogramming that leads to SA-induced immunity in plants. NPR1 is an SA receptor and master transcriptional regulator in SA-triggered transcriptional reprogramming. Despite the indispensable role of NPR1, genome-wide direct targets of NPR1 specific to SA signaling have not been identified. Here, we report INA (functional SA analog)-specific genome-wide targets of Arabidopsis NPR1 in plants expressing GFP-fused NPR1 under its native promoter. Analyses of NPR1-dependently expressed direct NPR1 targets revealed that NPR1 primarily activates genes encoding transcription factors upon INA treatment, triggering transcriptional cascades required for INA-induced transcriptional reprogramming and immunity. We identified genome-wide targets of a histone acetyltransferase, HAC1, including hundreds of co-targets shared with NPR1, and showed that NPR1 and HAC1 regulate INA-induced histone acetylation and expression of a subset of the co-targets. Genomic NPR1 targeting was principally mediated by TGACG-motif binding protein (TGA) transcription factors. Furthermore, a group of NPR1 targets mostly encoding transcriptional regulators was already bound to NPR1 in the basal state and showed more rapid and robust induction than other NPR1 targets upon SA signaling. Thus, our study unveils genome-wide NPR1 targeting, its role in transcriptional reprogramming, and the cooperativity between NPR1, HAC1, and TGAs in INA-induced immunity.

Genome-Wide Identification of PAP1 Direct Targets in Regulating Seed Anthocyanin Biosynthesis in Arabidopsis.

Anthocyanins are widespread water-soluble pigments in the plant kingdom. Anthocyanin accumulation is activated by the MYB-bHLH-WD40 (MBW) protein complex. In Arabidopsis, the R2R3-MYB transcription factor PAP1 activates anthocyanin biosynthesis. While prior research primarily focused on seedlings, seeds received limited attention. This study explores PAP1's genome-wide target genes in anthocyanin biosynthesis in seeds. Our findings confirm that PAP1 is a positive regulator of anthocyanin biosynthesis in Arabidopsis seeds. PAP1 significantly increased anthocyanin content in developing and mature seeds in Arabidopsis. Transcriptome analysis at 12 days after pollination reveals the upregulation of numerous genes involved in anthocyanin accumulation in 35S:PAP1 developing seeds. Chromatin immunoprecipitation and dual luciferase reporter assays demonstrate PAP1's direct promotion of ten key genes and indirect upregulation of TT8, TTG1, and eight key genes during seed maturation, thus enhancing seed anthocyanin accumulation. These findings enhance our understanding of PAP1's novel role in regulating anthocyanin accumulation in Arabidopsis seeds.

Resveratrol and 3,3'-Diindolylmethane Differentially Regulate Aryl Hydrocarbon Receptor and Estrogen Receptor Alpha Activity through Multiple Transcriptomic Targets in MCF-7 Human Breast Cancer Cells.

Inhibitory crosstalk between estrogen receptor alpha (ERα) and aryl hydrocarbon receptor (AHR) regulates 17β-estradiol (E2)-dependent breast cancer cell signaling. ERα and AHR are transcription factors activated by E2 and 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), respectively. Dietary ligands resveratrol (RES) and 3,3'diindolylmethane (DIM) also activate ERα while only DIM activates AHR and RES represses it. DIM and RES are reported to have anti-cancer and anti-inflammatory properties. Studies with genome-wide targets and AHR- and ERα-regulated genes after DIM and RES are unknown. We used chromatin immunoprecipitation with high-throughput sequencing and transcriptomics to study ERα as well as AHR coregulation in MCF-7 human breast cancer cells treated with DIM, RES, E2, or TCDD alone or E2+TCDD for 1 and 6 h, respectively. ERα bound sites after being DIM enriched for the AHR motif but not after E2 or RES while AHR bound sites after being DIM and E2+TCDD enriched for the ERE motif but not after TCDD. More than 90% of the differentially expressed genes closest to an AHR binding site after DIM or E2+TCDD also had an ERα site, and 60% of the coregulated genes between DIM and E2+TCDD were common. Collectively, our data show that RES and DIM differentially regulate multiple transcriptomic targets via ERα and ERα/AHR coactivity, respectively, which need to be considered to properly interpret their cellular and biological responses. These novel data also suggest that, when both receptors are activated, ERα dominates with preferential recruitment of AHR to ERα target genes.

Elucidation of genome-wide understudied proteins targeted by PROTAC-induced degradation using interpretable machine learning.

Proteolysis-targeting chimeras (PROTACs) are hetero-bifunctional molecules that induce the degradation of target proteins by recruiting an E3 ligase. PROTACs have the potential to inactivate disease-related genes that are considered undruggable by small molecules, making them a promising therapy for the treatment of incurable diseases. However, only a few hundred proteins have been experimentally tested for their amenability to PROTACs, and it remains unclear which other proteins in the entire human genome can be targeted by PROTACs. In this study, we have developed PrePROTAC, an interpretable machine learning model based on a transformer-based protein sequence descriptor and random forest classification. PrePROTAC predicts genome-wide targets that can be degraded by CRBN, one of the E3 ligases. In the benchmark studies, PrePROTAC achieved a ROC-AUC of 0.81, an average precision of 0.84, and over 40% sensitivity at a false positive rate of 0.05. When evaluated by an external test set which comprised proteins from different structural folds than those in the training set, the performance of PrePROTAC did not drop significantly, indicating its generalizability. Furthermore, we developed an embedding SHapley Additive exPlanations (eSHAP) method, which extends conventional SHAP analysis for original features to an embedding space through in silico mutagenesis. This method allowed us to identify key residues in the protein structure that play critical roles in PROTAC activity. The identified key residues were consistent with existing knowledge. Using PrePROTAC, we identified over 600 novel understudied proteins that are potentially degradable by CRBN and proposed PROTAC compounds for three novel drug targets associated with Alzheimer's disease.

SIX1 and EWS/FLI1 co-regulate an anti-metastatic gene network in Ewing Sarcoma

Ewing sarcoma (ES), which is characterized by the presence of oncogenic fusion proteins such as EWS/FLI1, is an aggressive pediatric malignancy with a high rate of early dissemination and poor outcome after distant spread. Here we demonstrate that the SIX1 homeoprotein, which enhances metastasis in most tumor types, suppresses ES metastasis by co-regulating EWS/FLI1 target genes. Like EWS/FLI1, SIX1 promotes cell growth/transformation, yet dramatically inhibits migration and invasion, as well as metastasis in vivo. We show that EWS/FLI1 promotes SIX1 protein expression, and that the two proteins share genome-wide binding profiles and transcriptional regulatory targets, including many metastasis-associated genes such as integrins, which they co-regulate. We further show that SIX1 downregulation of integrins is critical to its ability to inhibit invasion, a key characteristic of metastatic cells. These data demonstrate an unexpected anti-metastatic function for SIX1, through coordinate gene regulation with the key oncoprotein in ES, EWS/FLI1.

Regulation of acetate tolerance by small ORF-encoded polypeptides modulating efflux pump specificity in Methylomonas sp. DH-1

BackgroundMethanotrophs have emerged as promising hosts for the biological conversion of methane into value-added chemicals, including various organic acids. Understanding the mechanisms of acid tolerance is essential for improving organic acid production. WatR, a LysR-type transcriptional regulator, was initially identified as involved in lactate tolerance in a methanotrophic bacterium Methylomonas sp. DH-1. In this study, we investigated the role of WatR as a regulator of cellular defense against weak organic acids and identified novel target genes of WatR.ResultsBy conducting an investigation into the genome-wide binding targets of WatR and its role in transcriptional regulation, we identified genes encoding an RND-type efflux pump (WatABO pump) and previously unannotated small open reading frames (smORFs), watS1 to watS5, as WatR target genes activated in response to acetate. The watS1 to watS5 genes encode polypeptides of approximately 50 amino acids, and WatS1 to WatS4 are highly homologous with one predicted transmembrane domain. Deletion of the WatABO pump genes resulted in decreased tolerance against formate, acetate, lactate, and propionate, suggesting its role as an efflux pump for a wide range of weak organic acids. WatR repressed the basal expression of watS genes but activated watS and WatABO pump genes in response to acetate stress. Overexpression of watS1 increased tolerance to acetate but not to other acids, only in the presence of the WatABO pump. Therefore, WatS1 may increase WatABO pump specificity toward acetate, switching the general weak acid efflux pump to an acetate-specific efflux pump for efficient cellular defense against acetate stress.ConclusionsOur study has elucidated the role of WatR as a key transcription factor in the cellular defense against weak organic acids, particularly acetate, in Methylomonas sp. DH-1. We identified the genes encoding WatABO efflux pump and small polypeptides (WatS1 to WatS5), as the target genes regulated by WatR for this specific function. These findings offer valuable insights into the mechanisms underlying weak acid tolerance in methanotrophic bacteria, thereby contributing to the development of bioprocesses aimed at converting methane into value-added chemicals.

Insights into the Diversity of Transcription Activator-Like Effectors in Indian Pathotype Strains of Xanthomonas oryzae pv. oryzae.

Xanthomonas oryzae pv. oryzae (Xoo) is a major rice pathogen, and its genome harbors extensive inter-strain and inter-lineage variations. The emergence of highly virulent pathotypes of Xoo that can overcome major resistance (R) genes deployed in rice breeding programs is a grave threat to rice cultivation. The present study reports on a long-read Oxford nanopore-based complete genomic investigation of Xoo isolates from 11 pathotypes that are reported based on their reaction toward 10 R genes. The investigation revealed remarkable variation in the genome structure in the strains belonging to different pathotypes. Furthermore, transcription activator-like effector (TALE) proteins secreted by the type III secretion system display marked variation in content, genomic location, classes, and DNA-binding domain. We also found the association of tal genes in the vicinity of regions with genome structural variations. Furthermore, in silico analysis of the genome-wide rice targets of TALEs allowed us to understand the emergence of pathotypes compatible with major R genes. Long-read, cost-effective sequencing technologies such as nanopore can be a game changer in the surveillance of major and emerging pathotypes. The resource and findings will be invaluable in the management of Xoo and in appropriate deployment of R genes in rice breeding programs.

HNF1A binds and regulates the expression of SLC51B to facilitate the uptake of estrone sulfate in human renal proximal tubule epithelial cells

Renal defects in maturity onset diabetes of the young 3 (MODY3) patients and Hnf1a-/- mice suggest an involvement of HNF1A in kidney development and/or its function. Although numerous studies have leveraged on Hnf1α-/- mice to infer some transcriptional targets and function of HNF1A in mouse kidneys, species-specific differences obviate a straightforward extrapolation of findings to the human kidney. Additionally, genome-wide targets of HNF1A in human kidney cells have yet to be identified. Here, we leveraged on human in vitro kidney cell models to characterize the expression profile of HNF1A during renal differentiation and in adult kidney cells. We found HNF1A to be increasingly expressed during renal differentiation, with peak expression on day 28 in the proximal tubule cells. HNF1A ChIP-Sequencing (ChIP-Seq) performed on human pluripotent stem cell (hPSC)-derived kidney organoids identified its genome-wide putative targets. Together with a qPCR screen, we found HNF1A to activate the expression of SLC51B, CD24, and RNF186 genes. Importantly, HNF1A-depleted human renal proximal tubule epithelial cells (RPTECs) and MODY3 human induced pluripotent stem cell (hiPSC)-derived kidney organoids expressed lower levels of SLC51B. SLC51B-mediated estrone sulfate (E1S) uptake in proximal tubule cells was abrogated in these HNF1A-deficient cells. MODY3 patients also exhibit significantly higher excretion of urinary E1S. Overall, we report that SLC51B is a target of HNF1A responsible for E1S uptake in human proximal tubule cells. As E1S serves as the main storage form of nephroprotective estradiol in the human body, lowered E1S uptake and increased E1S excretion may reduce the availability of nephroprotective estradiol in the kidneys, contributing to the development of renal disease in MODY3 patients.

Novel mechanisms of androgen receptor-centered transcriptional regulatory network in regulating CD8 +T cell exhaustion and sex bias in cancer

Abstract Sex bias in cancers arising from nonreproductive organs is known but poorly understood. We recently reported on a T cell-intrinsic role of androgen receptor (AR) in driving CD8 +T cell exhaustion which underlies the poorer tumor control in males (Kwon et al., Sci Immunol, 2022). This study aims to understand how AR regulates CD8 +tumor-infiltrating lymphocytes (TILs) by identifying its genome-wide targets. We created a CD8 +T cell-specific AR knockout (KO) mouse model by crossing E8I-cre mice with Ar-floxed mice, which we challenged with syngeneic bladder tumor MB49. We then monitored tumor growth and performed spectral flow cytometry using a T cell exhaustion panel. Further, we used Cleavage Under Targets and Tagmentation – sequencing (CUT&Tag-seq) to map the entire AR targets in CD8 +TILs. Loss of AR in CD8 +T cells significantly slowed the growth of MB49 in male but not female mice. CD8 +TILs from male CD8 AR KO mice showed reduced TOX +TCF1 −terminally exhausted subset (69%; P≤ 0.05) and TOX expression (72%; P = 0.0583) compared to controls. Similarly, TOX expression decreased by 43% in AR-deleted CD8 +T cells following chronic TCR stimulation in vitro (P ≤ 0.0001). Finally, using CUT&Tag-seq, we found that AR binds directly to promoters of multiple key transcriptional regulators of T cell exhaustion, including Tcf7 and Tox. Through identifying the genome-wide targets of AR in CD8 +TILs, we characterized the cell-intrinsic mechanism by which AR regulates CD8 +T cell exhaustion. Our work suggests AR is a master transcriptional regulator that negatively impacts CD8 +T cell immunity, which may contribute to poorer survival in males and supports combining androgen deprivation therapy with immune checkpoint blockage for cancer treatment. Special Institutional Funding

DeepCGP: A Deep Learning Method to Compress Genome-Wide Polymorphisms for Predicting Phenotype of Rice.

Genomic selection (GS) is expected to accelerate plant and animal breeding. During the last decade, genome-wide polymorphism data have increased, which has raised concerns about storage cost and computational time. Several individual studies have attempted to compress the genome data and predict phenotypes. However, compression models lack adequate quality of data after compression, and prediction models are time consuming and use original data to predict the phenotype. Therefore, a combined application of compression and genomic prediction modeling using deep learning could resolve these limitations. A Deep Learning Compression-based Genomic Prediction (DeepCGP) model that can compress genome-wide polymorphism data and predict phenotypes of a target trait from compressed information was proposed. The DeepCGP model contained two parts: (i) an autoencoder model based on deep neural networks to compress genome-wide polymorphism data, and (ii) regression models based on random forests (RF), genomic best linear unbiased prediction (GBLUP), and Bayesian variable selection (BayesB) to predict phenotypes from compressed information. Two datasets with genome-wide marker genotypes and target trait phenotypes in rice were applied. The DeepCGP model obtained up to 99% prediction accuracy to the maximum for a trait after 98% compression. BayesB required extensive computational time among the three methods, and showed the highest accuracy; however, BayesB could only be used with compressed data. Overall, DeepCGP outperformed state-of-the-art methods in terms of both compression and prediction. Our code and data are available at https://github.com/tanzilamohita/DeepCGP.

Unique and overlapping functions for the transcriptional regulators KANADI1 and ULTRAPETALA1 in Arabidopsis gynoecium and stamen gene regulation.

Plants generate their reproductive organs, the stamens and the carpels, de novo within the flowers that form when the plant reaches maturity. The carpels comprise the female reproductive organ, the gynoecium, a complex organ that develops along several axes of polarity and is crucial for plant reproduction, fruit formation, and seed dispersal. The epigenetic trithorax group (trxG) protein ULTRAPETALA1 (ULT1) and the GARP domain transcription factor KANADI1 (KAN1) act cooperatively to regulate Arabidopsis thaliana gynoecium patterning along the apical-basal polarity axis; however, the molecular pathways through which this patterning activity is achieved remain to be explored. In this study, we used transcriptomics to identify genome-wide ULT1 and KAN1 target genes during reproductive development. We discovered 278 genes in developing flowers that are regulated by ULT1, KAN1, or both factors together. Genes involved in developmental and reproductive processes are overrepresented among ULT1 and/or KAN1 target genes, along with genes involved in biotic or abiotic stress responses. Consistent with their function in regulating gynoecium patterning, a number of the downstream target genes are expressed in the developing gynoecium, including a unique subset restricted to the stigmatic tissue. Further, we also uncovered a number of KAN1- and ULT1-induced genes that are transcribed predominantly or exclusively in developing stamens. These findings reveal a potential cooperative role for ULT1 and KAN1 in male as well as female reproductive development that can be investigated with future genetic and molecular experiments.

Biotinylated Tn5 transposase-mediated CUT&Tag efficiently profiles transcription factor-DNA interactions in plants.

In contrast to CUT&Tag approaches for profiling bulk histone modifications, current CUT&Tag methods for analysing specific transcription factor (TF)-DNA interactions remain technically challenging due to TFs having relatively low abundance. Moreover, an efficient CUT&Tag strategy for plant TFs is not yet available. Here, we first applied biotinylated Tn5 transposase-mediated CUT&Tag (B-CUT&Tag) to produce high-quality libraries for interrogating TF-DNA interactions. B-CUT&Tag combines streptavidin-biotin-based DNA purification with routine CUT&Tag, optimizing the removal of large amounts of intact chromatin not targeted by specific TFs. The biotinylated chromatin fragments are then purified for construction of deep sequencing libraries or qPCR analysis. We applied B-CUT&Tag to probe genome-wide DNA targets of Squamosa promoter-binding-like protein 9 (SPL9), a well-established TF in Arabidopsis; the resulting profiles were efficient and consistent in demonstrating its well-established target genes in juvenile-adult transition/flowering, trichome development, flavonoid biosynthesis, wax synthesis and branching. Interestingly, our results indicate functions of AtSPL9 in modulating growth-defence trade-offs. In addition, we established a method for applying qPCR after CUT&Tag (B-CUT&Tag-qPCR) and successfully validated the binding of SPL9 in Arabidopsis and PHR2 in rice. Our study thus provides a convenient and highly efficient CUT&Tag strategy for profiling TF-chromatin interactions that is widely applicable to the annotation of cis-regulatory elements for crop improvement.

Genomic regulation of Krüppel-like-factor family members by corticosteroid receptors in the rat brain

Hippocampal mineralocorticoid receptors (MRs) and glucocorticoid receptors (GRs) mediate glucocorticoid hormone (GC) action in the hippocampus. These receptors bind to glucocorticoid responsive elements (GREs) within target genes, eliciting transcriptional effects in response to stress and circadian variation. Until recently, little was known about the genome-wide targets of hippocampal MRs and GRs under physiological conditions. Following on from our genome-wide MR and GR ChIP-seq and Ribo-Zero RNA-seq studies on rat hippocampus, we investigated the Krüppel-like factors (KLFs) as targets of MRs and GRs throughout the brain under circadian variation and after acute stress. In particular, Klf2, Klf9 and Klf15 are known to be stress and/or GC responsive and play a role in neurobiological processes including synaptic plasticity and neuronal differentiation. We found increased binding of MR and GR to GREs within Klf2, Klf9 and Klf15 in the hippocampus, amygdala, prefrontal cortex, and neocortex after acute stress and resulting from circadian variation, which was accompanied by upregulation of corresponding hnRNA and mRNA levels. Adrenalectomy abolished transcriptional upregulation of specific Klf genes. These results show that MRs and GRs regulate Klf gene expression throughout the brain following exposure to acute stress or in response to circadian variation, likely alongside other transcription factors.