Gene Pairs Research Articles

Decrease in purifying selection pressures on wheat homoeologous genes: tetraploidization versus hexaploidization.

A series of polyploidizations in higher-order polyploids is the main event affecting gene content in a genome. Each polyploidization event can lead to massive functional divergence because of the subsequent decrease in selection pressure on duplicated genes; however, the causal relationship between multiple rounds of polyploidization and the functional divergence of duplicated genes is poorly understood. We focused on the Triticum-Aegilops complex lineage and compared selection pressure before and after tetraploidization and hexaploidization events. Although both events led to decreased selection pressure on homoeologous gene pairs (compared with diploids and tetraploids), the initial tetraploidization had a greater impact on selection pressure on homoeologous gene pairs than did subsequent hexaploidization. Consistent with this, selection pressure on expression patterns for the initial event relaxed more than those for the subsequent event. Surprisingly, the decreased selection pressure on these homoeologous genes was independent of the existence of in-paralogs within the same subgenome. Wheat homoeologous pairs had different evolutionary consequences compared with orthologs related to other mechanisms (ancient allopolyploidization, ancient autopolyploidization, and small-scale duplication). Furthermore, tetraploidization and hexaploidization also seemed to have different evolutionary consequences. This suggests that homoeologous genes retain unique functions, including functions that are unlikely to be preserved in genes generated by the other duplication mechanisms. We found that their unique functions differed between tetraploidization and hexaploidization (e.g., reproductive and chromosome segregation processes). These findings imply that the substantial number of gene pairs resulting from multiple allopolyploidization events, especially initial tetraploidization, may have been a unique source of functional divergence.

A type II toxin-antitoxin system, ECs3274-ECs3275, in enterohemorrhagic Escherichia coli O157.

The toxin-antitoxin (TA) genetic module controls various bacterial events. Novel toxins with different functions are still being discovered. This study aimed to determine whether the ECs3274-ECs3275 gene pair encoded by enterohemorrhagic Escherichia coli O157 functions as a TA system. To characterize this putative TA system, we analyzed the growth of E. coli expressing ECs3274, ECs3275, or both; the interaction between ECs3274 and ECs3275 using bacterial adenylate cyclase two-hybrid assays; and the DNA-binding ability of ECs3274 using gel-mobility shift assays. We observed that the ECs3274 antitoxin interacted with the ECs3275 toxin, was destabilized by Lon protease, and repressed its promoter activity via its helix-turn-helix (HTH) motif. These properties are consistent with those of typical type II TA antitoxins. Interestingly, ECs3275 has an HTH motif not observed in other TA toxins and is necessary for ECs3275 toxicity, suggesting that ECs3275 may exert its toxicity by regulating the expression of specific genes.

Genome-wide exploration of the CONSTANS-like (COL) gene family and its potential role in regulating plant flowering time in foxtail millet (Setaria italica)

In photoperiod-sensitive plants, the CO-like gene (CONSTANS-like, COL) has a crucial function in regulating the timing of flowering. The blooming period in foxtail millet is greatly influenced by the duration of daylight; however, there is a scarcity of data regarding the molecular properties of the COL genes in the foxtail millet. In this study, we conducted a comprehensive analysis of the COL gene family in foxtail millet at the genome-wide level. We identified 11 SiCOL genes and performed gene structure analysis, which showed pronounced variation in gene length and intron number among the genes. The examination of COL proteins in foxtail millet and other plant species using phylogenetic analysis revealed that they could be clustered into three distinct groups. Cis-element analysis identified elements related to light-responsiveness, hormones, and abiotic stress in the promoter region of the SiCOL gene. Furthermore, tissue-specific expression analysis showed widespread expression of all 11 SiCOL genes in various foxtail millet tissues and organs, particularly in leaves and panicles. Collinearity analysis identified 14 syntenic gene pairs in both foxtail millet and rice. The results also revealed diurnal oscillations in the transcription levels of SiCOL genes under different light conditions. Moreover, among the 11 genes, SiCO, SiCOL1, and SiCOL6 expression levels were negatively correlated with flowering time variation in two foxtail cultivars. Additionally, upon constructing a network of predicted molecular interactions, FLOWER LOCUS-like (FTL) and Phytochromea A (PHY A) were suggested to potentially interact with SiCO, SiCOL1, and SiCOL6. SiCO, SiCOL1, and SiCOL6 have the potential for flowering and heading in foxtail millet. This research enhances our comprehension of the role and control of the SiCOL gene family constituents in foxtail millet, establishing a basis for future investigations.

Identification and Characterization of the DOF Gene Family in Phoebe bournei and Its Role in Abiotic Stress-Drought, Heat and Light Stress.

Phoebe bournei is a second-class endangered and protected species unique to China, and it holds significant ecological and economic value. DNA binding one zinc finger (Dof) transcription factors are plant-specific regulators. Numerous studies have demonstrated that Dof genes are involved in plant growth, development and responses to abiotic stress. In this study, we identified and analyzed 34 PbDof gene members at the whole-genome level. The results indicated that the 34 PbDof genes were unevenly distributed across 12 chromosomes. We utilized the Dof genes from Arabidopsis thaliana and P. bournei to construct a phylogenetic tree and categorized these genes into eight subgroups. In the collinearity analysis, there were 16 homologous gene pairs between AtDof and PbDof and nine homologous gene pairs between ZmDof and PbDof. We conducted a cis-acting element analysis and found that cis-acting elements involved in light response were the most abundant in PbDof genes. Through SSR site prediction, we analyzed that the evolution level of Dof genes is low. Additionally, we assessed the expression profiles of eight PbDof genes under high temperature, drought, and light stress using qRT-PCR. In particular, PbDof08 and PbDof16 are significantly upregulated under the three stresses. This study provides foundational information for PbDof genes and offers new insights for further research on the mechanism of Dof transcription factors responding to stress, as well as the adaptation of P. bournei to environmental changes.

Partial correlation network analysis identifies coordinated gene expression within a regional cluster of COPD genome-wide association signals.

Chronic obstructive pulmonary disease (COPD) is a complex disease influenced by well-established environmental exposures (most notably, cigarette smoking) and incompletely defined genetic factors. The chromosome 4q region harbors multiple genetic risk loci for COPD, including signals near HHIP, FAM13A, GSTCD, TET2, and BTC. Leveraging RNA-Seq data from lung tissue in COPD cases and controls, we estimated the co-expression network for genes in the 4q region bounded by HHIP and BTC (~70MB), through partial correlations informed by protein-protein interactions. We identified several co-expressed gene pairs based on partial correlations, including NPNT-HHIP, BTC-NPNT and FAM13A-TET2, which were replicated in independent lung tissue cohorts. Upon clustering the co-expression network, we observed that four genes previously associated to COPD: BTC, HHIP, NPNT and PPM1K appeared in the same network community. Finally, we discovered a sub-network of genes differentially co-expressed between COPD vs controls (including FAM13A, PPA2, PPM1K and TET2). Many of these genes were previously implicated in cell-based knock-out experiments, including the knocking out of SPP1 which belongs to the same genomic region and could be a potential local key regulatory gene. These analyses identify chromosome 4q as a region enriched for COPD genetic susceptibility and differential co-expression.

Antibody sequence determinants of viral antigen specificity.

Throughout life, humans experience repeated exposure to viral antigens through infection and vaccination, resulting in the generation of diverse, antigen-specific antibody repertoires. A paramount feature of antibodies that enables their critical contributions in counteracting recurrent and novel pathogens, and consequently fostering their utility as valuable targets for therapeutic and vaccine development, is the exquisite specificity displayed against their target antigens. Yet, there is still limited understanding of the determinants of antibody-antigen specificity, particularly as a function of antibody sequence. In recent years, experimental characterization of antibody repertoires has led to novel insights into fundamental properties of antibody sequences but has been largely decoupled from at-scale antigen specificity analysis. Here, using the LIBRA-seq technology, we generated a large data set mapping antibody sequence to antigen specificity for thousands of B cells, by screening the repertoires of a set of healthy individuals against 20 viral antigens representing diverse pathogens of biomedical significance. Analysis uncovered virus-specific patterns in variable gene usage, gene pairing, somatic hypermutation, as well as the presence of convergent antiviral signatures across multiple individuals, including the presence of public antibody clonotypes. Notably, our results showed that, for B-cell receptors originating from different individuals but leveraging an identical combination of heavy and light chain variable genes, there is a specific CDRH3 identity threshold above which B cells appear to exclusively share the same antigen specificity. This finding provides a quantifiable measure of the relationship between antibody sequence and antigen specificity and further defines experimentally grounded criteria for defining public antibody clonality.IMPORTANCEThe B-cell compartment of the humoral immune system plays a critical role in the generation of antibodies upon new and repeated pathogen exposure. This study provides an unprecedented level of detail on the molecular characteristics of antibody repertoires that are specific to each of the different target pathogens studied here and provides empirical evidence in support of a 70% CDRH3 amino acid identity threshold in pairs of B cells encoded by identical IGHV:IGL(K)V genes, as a means of defining public clonality and therefore predicting B-cell antigen specificity in different individuals. This is of exceptional importance when leveraging public clonality as a method to annotate B-cell receptor data otherwise lacking antigen specificity information. Understanding the fundamental rules of antibody-antigen interactions can lead to transformative new approaches for the development of antibody therapeutics and vaccines against current and emerging viruses.

Genome-wide identification of AP2/ERF gene family in Coptis Chinensis Franch reveals its role in tissue-specific accumulation of benzylisoquinoline alkaloids

BackgroundThe Plant-specific AP2/ERF gene family encodes proteins involved in various biological and physiological processes. Although the genome of Coptis chinensis Franch, a plant producing benzylisoquinoline alkaloids (BIAs), has been sequenced at the chromosome level, studies on the AP2/ERF gene family in C. chinensis are lacking. Thus, a genome-wide identification of AP2/ERF gene family in C. chinensis was conducted to explore its role in BIAs biosynthesis.ResultsA total of 96 CcAP2/ERF genes were identified and categorized into five subfamilies, including 43 ERFs, 32 DREBs, 17 AP2s, 3 RAVs, and 1 Soloist, based on their structural domains. These CcAP2/ERF genes were unevenly distributed across nine chromosomes. Analysis of gene duplication events identified 17 CcAP2/ERF gene pairs in the genome, with 7 involved in tandem duplication events and 10 involved in segmental duplicate events, indicating that both types of duplications contributed to the expansion of the AP2/ERF gene family. The Ka/Ks ratio analysis suggested that the CcAP2/ERF gene family underwent strong purifying selection. Two phytohormones, methyl jasmonate and abscisic acid, were identified as potential key inducers of BIAs biosynthesis due to the cis-acting element prediction. Analysis of the spatial transcriptomic data revealed that 28 differentially expressed AP2/ERF genes had the highest or relatively higher expression levels in the rhizome, 17 of which positively correlated with the tissue-specific accumulation of BIAs. Further real-time PCR verification and protein-protein interaction analysis indicated that DREB1B might be one of the central regulators in a highly complex BIAs biosynthesis network.ConclusionThese findings provide significant insight into the function of AP2/ERF genes in C. chinensis, particularly in the regulatory network of BIAs biosynthesis in C. chinensis. This study also identifies candidate genes for metabolic engineering to increase BIAs content in C. chinensis.

The puzzle of two tandem acyl-CoA ligases of Pseudomonas putida F1.

Pseudomonas bacteria demonstrate a great deal of metabolic diversity and consequently colonize a wide range of ecological niches. A characteristic of these bacteria is a pair of genes in tandem annotated as acyl-CoA ligases involved in fatty acid degradation. The Pseudomonas putida F1 genome is annotated as having at least nine genes encoding acyl-CoA ligases which are scattered around the chromosome excepting the tandem pair. Since similar tandem pairs are found in other pseudomonads, we have constructed and characterized deletion mutants of the tandem ligases. We report that the encoded proteins are authentic acyl-CoA ligases involved in fatty acid degradation.

Transcriptional Analysis of Different Cultivars and Genome-Wide Identification of FAD Gene Family in Tree Peony

The tree peony (Paeonia ostii), a newly recognized woody oil plant endemic to China, is noteworthy for its high content of unsaturated fatty acids (UFA), particularly alpha-linolenic acid (ALA). Fatty acid desaturases (FADs) are integral to plant development and defense mechanisms. Nonetheless, there is limited understanding of (i) the molecular mechanism underlying FA biosynthesis in various varieties during seed maturation and (ii) a genome-wide analysis of FAD family genes within the tree peony. We selected three distinct cultivars of tree peony for transcriptome sequencing and performed an extensive analysis of PoFAD genes. In total, 67,542 unigenes were acquired and annotated with six protein databases available to the public. Forty-one differentially expressed genes (DEGs) pertinent to FA biosynthesis and lipid metabolism were identified in this study. Notably, genes such as PoFAD2, PoFAD6, and PoSAD were found to be significantly upregulated, contributing to a differential linolenic acid and linoleic acid content across the three cultivars. Herein, 24 PoFADs from the P. ostii genome were recognized and categorized into four distinct clusters according to their conserved structural features. The distribution of PoFADs was found to be random and uneven across five chromosomes, indicating a complex genomic architecture. Six colinear gene pairs were found between P. ostii and V. vinifera, indicating a potential link due to their close relationship. Together, these findings significantly enhance our knowledge of the molecular processes governing fatty acid synthesis, elucidate the functional roles of the FAD gene family, and lay the groundwork for using genetic manipulation to boost lipid levels.

Genome-wide identification and analysis of anthocyanin synthesis-related R2R3-MYB genes in Fragaria pentaphylla

BackgroundMYB transcription factors regulate anthocyanin biosynthesis across numerous plant species. However, comprehensive genome-wide investigations regarding the R2R3-MYB gene family and its involvement in regulating anthocyanin biosynthesis in the red and white fruit color morphs of Fragaria pentaphylla remain scarce.ResultsA total of 101 FpR2R3-MYB genes were identified from the F. pentaphylla genome and were divided into 34 subgroups based on phylogenetic analysis. Gene structure (exon/intron) and protein motifs were particularly conserved among the FpR2R3-MYB genes, especially members within the same subgroup. The FpR2R3-MYB genes were distributed over seven F. pentaphylla chromosomes. Analysis of gene duplication events revealed five pairs of tandem duplication genes and 16 pairs of segmental duplication genes, suggesting that segmental duplications are the major pattern for expansion of the FpR2R3-MYB gene family expansion in F. pentaphylla. Cis-regulatory elements of the FpR2R3-MYB promoters were involved in cellular development, phytohormones, environmental stress and photoresponse. Based on the analysis of the FpR2R3-MYB gene family and transcriptome sequencing (RNA-seq) data, FpMYB9 was identified as a key transcription factor involved in the regulation of anthocyanin synthesis in F. pentaphylla fruits. The expression of FpMYB9 increases significantly during the ripening stage of red fruits, as confirmed by reverse transcription quantitative real-time PCR. In addition, subcellular localization experiments further confirmed the nuclear presence of FpMYB9, supporting its role as a transcription factor involved in anthocyanin biosynthesis.ConclusionOur results showed that the FpR2R3-MYB genes are highly conserved and play important roles in the anthocyanin biosynthesis in F. pentaphylla fruits. Our results also provide a compelling basis for further understanding of the regulatory mechanism underlying the role of FpMYB9 in anthocyanin formation in F. pentaphylla fruits.

ELiAH: the atlas of E3 ligases in human tissues for targeted protein degradation with reduced off-target effect.

The development of therapeutic agents has mainly focused on designing small molecules to modulate target proteins or genes which are conventionally druggable. Therefore, targeted protein degradation (TPD) for undruggable cases has emerged as promising pharmaceutical approach. TPD, often referred PROTACs (PROteolysis TArgeting Chimeras), uses a linker to degrade target proteins by hijacking the ubiquitination system. Therefore, unravel the relationship including reversal and co-expression between E3 ligands and other possible target genes in various human tissues is essential to mitigate off-target effects of TPD. Here, we developed the atlas of E3 ligases in human tissues (ELiAH), to prioritize E3 ligase-target gene pairs for TPD. Leveraging over 2900 of RNA-seq profiles consisting of 11 human tissues from the GTEx (genotype-tissue expression) consortium, users of ELiAH can identify tissue-specific genes and E3 ligases (FDR P-value of Mann-Whitney test < .05). ELiAH unravels 933 830 relationships consisting of 614 E3 ligases and 20 924 of expressed genes considering degree of tissue specificity, which are indispensable for ubiquitination based TPD development. In addition, docking properties of those relationships are also modeled using RosettaDock. Therefore, ELiAH presents comprehensive repertoire of E3 ligases for ubiquitination-based TPD drug development avoiding off-target effects. Database URL: https://eliahdb.org.

Nuclear compensatory evolution driven by mito-nuclear incompatibilities

Mitochondrial function relies on the coordinated expression of mitochondrial and nuclear genes, exhibiting remarkable resilience despite high mitochondrial mutation rates. The nuclear compensation mechanism suggests deleterious mitochondrial alleles drive compensatory nuclear mutations to preserve mito-nuclear compatibility. However, prevalence and factors conditioning this phenomenon remain debated due to its conflicting evidence. Here, we investigate how mito-nuclear incompatibilities impact substitutions in a model for species radiation. Mating success depends on genetic compatibility (nuclear DNA) and spatial proximity. Populations evolve from partially compatible mito-nuclear states, simulating mitochondrial DNA (mtDNA) introgression. Mutations do not confer advantages nor disadvantages, but individual fecundity declines with increasing incompatibilities, selecting for mito-nuclear coordination. We find that selection for mito-nuclear compatibility affects each genome differently based on their initial state. In compatible gene pairs, selection reduces substitutions in both genomes, while in incompatible nuclear genes, it consistently promotes compensation, facilitated by more mismatches. Interestingly, high mitochondrial mutation rates can reduce nuclear compensation by increasing mtDNA rectification, while substitutions in initially compatible nuclear gene are boosted. Finally, the presence of incompatibilities accelerates species radiation, but equilibrium richness is not directly correlated to substitution rates, revealing the complex dynamics triggered by mitochondrial introgression and mito-nuclear coevolution. Our study provides a perspective on nuclear compensation and the role of mito-nuclear incompatibilities in speciation by exploring extreme scenarios and identifying trends that empirical data alone cannot reveal. We emphasize the challenges in detecting these dynamics and propose analyzing specific genomic signatures could shed light on this evolutionary process.

Genome-wide identification and expression analysis of the ZRT, IRT-like protein (ZIP) family in Nicotiana tabacum.

Iron (Fe) and Zinc (Zn) are essential micronutrients for plant growth and development. ZIP (ZRT, IRT-like protein) transporters, known for their role in the regulation of Zinc and Iron uptake, are pivotal in facilitating the absorption, transport, and maintenance of Fe/Zn homeostasis in plants. Nicotiana tabacum has been widely used as a model plant for gene function analysis; however, the tobacco ZIP genes have not been identified systematically. In this study, we have identified a comprehensive set of 32 NtZIP genes, which were phylogenetically categorized into three distinct clades. The gene structures, characterized by their exon/intron organization, and the protein motifs are relatively conserved, particularly among genes within the same clade. These NtZIP genes exhibit an uneven distribution across 12 chromosomes. The gene localization analysis revealed the presence of 11 pairs of homeologous locus genes and 7 pairs of tandem duplication genes within the genome. To further explore the functionality of these genes, real-time quantitative reverse transcription PCR was employed to assess their expression levels in roots subjected to metal deficiency. The results indicated that certain NtZIP genes are specifically upregulated in response to either Fe or Zn deficiency. Additionally, the presence of specific cis-elements within their promoter regions, such as the E-box associated with Fe deficiency response and the ZDRE box linked to Zn deficiency response, was identified. This study lays a foundational groundwork for future research into the biological functions of NtZIP genes in tobacco in micronutrient regulation and homeostasis.

The Role of PLATZ6 in Raffinose Family Oligosaccharides Loading of Leaves via PLATZ Family Characterization in Cucumber

The plant AT protein and zinc-binding protein (PLATZ) genes, a novel cluster of plant-specific zinc-finger-dependent DNA-binding proteins, play a crucial role in regulating stress response and plant development. However, there has been little study focus on the role of the cucumber PLATZ family in assimilating loading in leaves. (1) In this study, a total of 12 PLATZ genes were identified from the cucumber genome. The cucumber PLATZ genes were clustered into five groups, and unevenly distributed on five chromosomes. A single pair of cucumber PLATZ genes underwent segmental duplication. (2) The results of genome-wide expression analysis suggested that the cucumber PLATZ genes were widely expressed in a wide range of cucumber tissues, with three PLATZ (PLATZ2, PLATZ6, and PLATZ12) genes exhibiting high expression in the vascular tissues of cucumber leaves. PLATZ2, PLATZ6, and PLATZ12 proteins were primarily located in cytomembrane and nucleus. (3) In VIGS-PLATZ6 plants, the expression of Galactinol synthase 1 (GolS1) and STACHYOSE SYNTHASE (STS), two genes involved in the synthesis of raffinose family oligosaccharides (RFOs) were observed to be decreased in cucumber leaves. In conclusion, the comprehensive analysis of the cucumber PLATZ family and the preliminary functional verification of PLATZ6 lay the foundation for the molecular and physiological functions of cucumber PLATZ genes.

Genome-wide analysis of the WRKY gene family and their response to low-temperature stress in elephant grass.

Elephant grass (Cenchrus purpureus) is a perennial forage grass characterized by tall plants, high biomass and wide adaptability. Low-temperature stress severely limits elephant grass biomass and geographic distribution. WRKY is one of the largest families of plant-specific transcription factors and plays important roles in plant resistance to low-temperature. However, the understanding of the WRKY family in grasses is limited. In this study, we conducted a genome-wide characterization of WRKY proteins in elephant grass, including gene structure, phylogeny, expression, conserved motif organization, and functional annotation, to identify key CpWRKY candidates involved in cold tolerance. In this study, a total of 176 WRKY genes were identified in elephant grass. It was found that 172 were unevenly distributed across its 14 chromosomes, while the remaining 4 genes were not anchored to any chromosome. The genes were classified into three groups based on their WRKY conserved domains and zinc finger motifs. There were 12, 8, 19, 27, 12, 18 and 80 CpWRKYs belonging to group I, group IIa, group IIb, group IIc, group IId, group IIe and group III, respectively. We hypothesized that the ancient subgroup IIc WRKY gene is the ancestor of all WRKY genes in elephant grass. Most CpWRKYs in the same group have similar structure and motif composition. A total of 169 duplicate gene pairs were identified, suggesting that segmental duplication might have contributed to the expansion of the CpWRKY gene family. Ka/Ks analysis revealed that most of the CpWRKYs were subjected to purifying selection during the evolution. It was also found that six genes (CpWRKY51, CpWRKY81, CpWRKY100, CpWRKY101, CpWRKY140 and CpWRKY143) exhibited higher expression in roots compare to leaves, and were significantly induced by low temperature stress. Among them, CpWRKY81 had the highest expression under low-temperature stress, and its over-expression significantly enhanced the cold tolerance in yeast. In this study, we characterized WRKY genes in elephant grass and further investigated their physicochemical properties, evolution, and expression patterns under low-temperature stress. This research provides valuable resources for identifying key CpWRKY genes that contribute to cold tolerance in elephant grass.

Genome-Wide Identification and Expression Analysis of the BTB Gene Superfamily Provides Insight into Sex Determination and Early Gonadal Development of Alligator sinensis.

The BTB gene superfamily is widely distributed among higher eukaryotes and plays a significant role in numerous biological processes. However, there is limited knowledge about the structure and function of BTB genes in the critically endangered species Alligator sinensis, which is endemic to China. A total of 170 BTB genes were identified from the A. sinensis genome, classified into 13 families, and unevenly distributed across 16 chromosomes. Analysis of gene duplication events yielded eight pairs of tandem duplication genes and six pairs of segmental duplication genes. Phylogenetics shows that the AsBTB genes are evolutionarily conserved. The cis-regulatory elements in the AsBTB family promoter region reveal their involvement in multiple biological processes. Protein interaction network analysis indicates that the protein interactions of the AsBTB genes are centered around CLU-3, mainly participating in the regulation of biological processes through the ubiquitination pathway. The expression profile and protein interaction network analysis of AsBTB genes during sex differentiation and early gonadal development indicate that AsBTB genes are widely expressed in this process and involves numerous genes and pathways for regulation. This study provides a basis for further investigation of the role of the BTB gene in sex differentiation and gonadal development in A. sinensis.

Harnessing the power of AI in precision medicine: NGS-based therapeutic insights for colorectal cancer cohort.

Developing innovative precision and personalized cancer therapeutics is essential to enhance cancer survivability, particularly for prevalent cancer types such as colorectal cancer. This study aims to demonstrate various approaches for discovering new targets for precision therapies using artificial intelligence (AI) on a Polish cohort of colorectal cancer patients. We analyzed 71 patients with histopathologically confirmed advanced resectional colorectal adenocarcinoma. Whole exome sequencing was performed on tumor and peripheral blood samples, while RNA sequencing (RNAseq) was conducted on tumor samples. We employed three approaches to identify potential targets for personalized and precision therapies. First, using our in-house neoantigen calling pipeline, ARDentify, combined with an AI-based model trained on immunopeptidomics mass spectrometry data (ARDisplay), we identified neoepitopes in the cohort. Second, based on recurrent mutations found in our patient cohort, we selected corresponding cancer cell lines and utilized knock-out gene dependency scores to identify synthetic lethality genes. Third, an AI-based model trained on cancer cell line data was employed to identify cell lines with genomic profiles similar to selected patients. Copy number variants and recurrent single nucleotide variants in these cell lines, along with gene dependency data, were used to find personalized synthetic lethality pairs. We identified approximately 8,700 unique neoepitopes, but none were shared by more than two patients, indicating limited potential for shared neoantigenic targets across our cohort. Additionally, we identified three synthetic lethality pairs: the well-known APC-CTNNB1 and BRAF-DUSP4 pairs, along with the recently described APC-TCF7L2 pair, which could be significant for patients with APC and BRAF variants. Furthermore, by leveraging the identification of similar cancer cell lines, we uncovered a potential gene pair, VPS4A and VPS4B, with therapeutic implications. Our study highlights three distinct approaches for identifying potential therapeutic targets in cancer patients. Each approach yielded valuable insights into our cohort, underscoring the relevance and utility of these methodologies in the development of precision and personalized cancer therapies. Importantly, we developed a novel AI model that aligns tumors with representative cell lines using RNAseq and methylation data. This model enables us to identify cell lines closely resembling patient tumors, facilitating accurate selection of models needed for in vitro validation.

9284 Leveraging Inter-Individual Transcriptional Correlation Structure To Infer Discrete Signaling Mechanisms Across Metabolic Tissues

Abstract Disclosure: M. Zhou: None. I. Tamburini: None. C. Van: None. J. Molendijk: None. L.M. Velez: None. C. Nguyen: None. R. Yeo: None. C. Filho: None. A.L. Hevener: None. J. Justice: None. L.M. Sparks: None. E.E. Kershaw: None. D. Nicholas: None. M. Seldin: None. Inter-organ communication is a vital process to maintain physiologic homeostasis, and its dysregulation contributes to many human diseases. Given that circulating bioactive factors are stable in serum, occur naturally, and are easily assayed from blood, they present obvious focal molecules for therapeutic intervention and biomarker development. A major obstacle in the characterization of such soluble factors is that defining their tissues and pathways of action requires extensive experimental testing in cells and animal models. Recently, studies have shown that secreted proteins mediating inter-tissue signaling could be identified by “brute-force” surveys of all genes within RNA-sequencing measures across tissues within a population. Expanding on this intuition, we reasoned those parallel strategies could be leveraged to understand how individual genes mediate signaling across metabolic tissues through correlative analyses of genetic gene variation between individuals. Thus, genetics comparison of quantitative levels of gene expression relationships between organs in a population could aid in understanding cross-organ signaling by adopting a gene-centric approach. Here, we surveyed gene-gene genetic correlation structure for ∼6.1x10^12 gene pairs across 18 metabolic tissues in 310 human individuals and 7 tissues in 103 diverse strains of mice fed a normal chow or HFHS diet. where the variation of genes such as FGF21, ADIPOQ, GCG, and IL6 showed enrichments that recapitulate experimental observations. Further, similar analyses were applied to explore both local within-tissue signaling mechanisms (liver PCSK9) as well as genes encoding enzymes producing metabolites (adipose PNPLA2), where genetic inter-individual correlation structure aligned with known roles for these critical metabolic pathways. Examination of sex hormone receptor correlations in mice highlighted the difference in tissue-specific variation in relationships with metabolic traits. Finally, we utilized this resource to suggest new functions for metabolic coordination between organs. For example, we prioritized key proteins for putative signaling between skeletal muscle and hippocampus, and further suggest colon as a central coordinator for systemic circadian clocks. We refer to this resource as Genetically-Derived Correlations Across Tissues (GD-CAT) where all tools and data are built into a web portal enabling users to perform these analyses without a single line of code (gdcat.org). This resource enables the querying of any gene in any tissue to find genetic coregulation correlated patterns of genes, cell types, pathways, and network architectures across metabolic organs. Presentation: 6/1/2024

7226 Leveraging Inter-Individual Transcriptional Correlation Structure To Infer Discrete Signaling Mechanisms Across Metabolic Tissues

Abstract Disclosure: M. Zhou: None. I. Tamburini: None. C. Van: None. J. Molendijk: None. L.M. Velez: None. C. Nguyen: None. R. Yeo: None. C. Filho: None. A.L. Hevener: None. J. Justice: None. L.M. Sparks: None. E.E. Kershaw: None. D. Nicholas: None. M. Seldin: None. Inter-organ communication is a vital process to maintain physiologic homeostasis, and its dysregulation contributes to many human diseases. Given that circulating bioactive factors are stable in serum, occur naturally, and are easily assayed from blood, they present obvious focal molecules for therapeutic intervention and biomarker development. A major obstacle in the characterization of such soluble factors is that defining their tissues and pathways of action requires extensive experimental testing in cells and animal models. Recently, studies have shown that secreted proteins mediating inter-tissue signaling could be identified by “brute-force” surveys of all genes within RNA-sequencing measures across tissues within a population. Expanding on this intuition, we reasoned those parallel strategies could be leveraged to understand how individual genes mediate signaling across metabolic tissues through correlative analyses of genetic gene variation between individuals. Thus, genetics comparison of quantitative levels of gene expression relationships between organs in a population could aid in understanding cross-organ signaling by adopting a gene-centric approach. Here, we surveyed gene-gene genetic correlation structure for ∼6.1x10^12 gene pairs across 18 metabolic tissues in 310 human individuals and 7 tissues in 103 diverse strains of mice fed a normal chow or HFHS diet. where the variation of genes such as FGF21, ADIPOQ, GCG, and IL6 showed enrichments that recapitulate experimental observations. Further, similar analyses were applied to explore both local within-tissue signaling mechanisms (liver PCSK9) as well as genes encoding enzymes producing metabolites (adipose PNPLA2), where genetic inter-individual correlation structure aligned with known roles for these critical metabolic pathways. Examination of sex hormone receptor correlations in mice highlighted the difference in tissue-specific variation in relationships with metabolic traits. Finally, we utilized this resource to suggest new functions for metabolic coordination between organs. For example, we prioritized key proteins for putative signaling between skeletal muscle and hippocampus, and further suggest colon as a central coordinator for systemic circadian clocks. We refer to this resource as Genetically-Derived Correlations Across Tissues (GD-CAT) where all tools and data are built into a web portal enabling users to perform these analyses without a single line of code (gdcat.org). This resource enables the querying of any gene in any tissue to find genetic coregulation correlated patterns of genes, cell types, pathways, and network architectures across metabolic organs. Presentation: 6/1/2024

Tumor draining lymph nodes connected to cold triple-negative breast cancers are characterized by Th2-associated microenvironment

Tumor draining lymph nodes (TDLN) represent a key component of the tumor-immunity cycle. There are few studies describing how TDLNs impact lymphocyte infiltration into tumors. Here we directly compare tumor-free TDLNs draining “cold” and “hot” human triple negative breast cancers (TDLNCold and TDLNHot). Using machine-learning-based self-correlation analysis of immune gene expression, we find unbalanced intranodal regulations within TDLNCold. Two gene pairs (TBX21/GATA3-CXCR1) with opposite correlations suggest preferential priming of T helper 2 (Th2) cells by mature dendritic cells (DC) within TDLNCold. This is validated by multiplex immunofluorescent staining, identifying more mature-DC-Th2 spatial clusters within TDLNCold versus TDLNHot. Associated with this Th2 priming preference, more IL4 producing mast cells (MC) are found within sinus regions of TDLNCold. Downstream, Th2-associated fibrotic TME is found in paired cold tumors with increased Th2/T-helper-1-cell (Th1) ratio, upregulated fibrosis growth factors, and stromal enrichment of cancer associated fibroblasts. These findings are further confirmed in a validation cohort and public genomic data. Our results reveal a potential role of IL4+ MCs within TDLNs, associated with Th2 polarization and reduced immune infiltration into tumors.