Functional Genomics Research Articles

Creating and leveraging bespoke large-scale knowledge graphs for comparative genomics and multi-omics drug discovery with SocialGene.

The rapid expansion of multi-omics data has transformed biological research, offering unprecedented opportunities to explore complex genomic relationships across diverse organisms. However, the vast volume and heterogeneity of these datasets presents significant challenges for analyses. Here we introduce SocialGene, a comprehensive software suite designed to collect, analyze, and organize multi-omics data into structured knowledge graphs, with the ability to handle small projects to repository-scale analyses. Originally developed to enhance genome mining for natural product drug discovery, SocialGene has been effective across various applications, including functional genomics, evolutionary studies, and systems biology. SocialGene's concerted Python and Nextflow libraries streamline data ingestion, manipulation, aggregation, and analysis, culminating in a custom Neo4j database. The software not only facilitates the exploration of genomic synteny but also provides a foundational knowledge graph supporting the integration of additional diverse datasets and the development of advanced search engines and analyses. This manuscript introduces some of SocialGene's capabilities through brief case studies including targeted genome mining for drug discovery, accelerated searches for similar and distantly related biosynthetic gene clusters in biobank-available organisms, integration of chemical and analytical data, and more. SocialGene is free, open-source, MIT-licensed, designed for adaptability and extension, and available from github.com/socialgene.

New insights on β-glycan synthases using in vitro GT-array (i-GT-ray) platform

Cellulose and hemicellulose are the major structural β-glycan polysaccharides in cell walls of land plants. They are characterized by a backbone of β-(1,3)- and/or β-(1,4)-linked sugars such as glucose, mannose, or xylose. The backbones of these polymers are produced by processive glycosyltransferases (GTs) called synthases having multiple transmembrane domains anchoring them to the membrane. Thus, they are among the most difficult membrane proteins to test in vitro and to purify. Recently, we developed an in vitro GT-array (i-GTray) platform and showed that non-processive type II membrane GTs could be produced via cell-free system in a soluble and active form and tested in this platform. To determine whether i-GT-ray platform is adequate for the production and testing of β-glycan synthases, we tested five synthases involved in cellulose, xyloglucan, (gluco)mannan, and β-(1,3)(1,4)-mixed-linkage glucan synthesis. Our results revealed unsuspected features of these enzymes. For example, all these synthases could be produced in a soluble and active form and are active in the absence of detergent or membrane lipids, and none of them required a primer for initiation of synthesis. All synthases produced ethanol-insoluble products that were susceptible to the appropriate hydrolases (i.e., cellulase, lichenase, mannanase). Using this platform, we showed that AtCslC4 and AtXXT1 interact directly to form an active xyloglucan synthase that produced xylosylated cello-oligosaccharides (up to three xylosyl residues) when supplied with UDP-Glc and UDP-Xyl. i-GTray platform represents a simple and powerful functional genomics tool for discovery of new insights of synthase activities and can be adapted to other enzymes.

The genomic potential of photosynthesis in piconanoplankton is functionally redundant but taxonomically structured at a global scale.

Carbon fixation is a key metabolic function shaping marine life, but the underlying taxonomic and functional diversity involved is only partially understood. Using metagenomic resources targeted at marine piconanoplankton, we provide a reproducible machine learning framework to derive the potential biogeography of genomic functions through the multi-output regression of gene read counts on environmental climatologies. Leveraging the Marine Atlas of Tara Oceans Unigenes, we investigate the genomic potential of primary production in the global ocean. The latter is performed by ribulose-1,5-bisphosphate carboxylase/oxygenase (RUBISCO) and is often associated with carbon concentration mechanisms in piconanoplankton, major marine unicellular photosynthetic organisms. We show that the genomic potential supporting C4 enzymes and RUBISCO exhibits strong functional redundancy and important affinity toward tropical oligotrophic waters. This redundancy is taxonomically structured by the dominance of Mamiellophyceae and Prymnesiophyceae in mid and high latitudes. These findings enhance our understanding of the relationship between functional and taxonomic diversity of microorganisms and environmental drivers of key biogeochemical cycles.

Two paralogous PHD finger proteins participate in natural genome editing in Paramecium tetraurelia.

The unicellular eukaryote Paramecium tetraurelia contains functionally distinct nuclei: germline micronuclei (MICs) and a somatic macronucleus (MAC). During sex, the MIC genome is reorganized into a new MAC genome and the old MAC is lost. Almost 45,000 unique internal eliminated sequences (IESs) distributed throughout the genome require precise excision to guarantee a functional new MAC genome. Here, we characterize a pair of paralogous PHD finger proteins involved in DNA elimination. DevPF1, the early-expressed paralog, is present in only some of the gametic and post-zygotic nuclei during meiosis. Both DevPF1 and DevPF2 localize in the new developing MACs, where IES excision occurs. Upon DevPF2 knockdown (KD), long IESs are preferentially retained and late-expressed small RNAs decrease; no length preference for retained IESs was observed in DevPF1-KD and development-specific small RNAs were abolished. The expression of at least two genes from the new MAC with roles in genome reorganization seems to be influenced by DevPF1- and DevPF2-KD. Thus, both PHD fingers are crucial for new MAC genome development, with distinct functions, potentially via regulation of non-coding and coding transcription in the MICs and new MACs.

OpenNucleome for high-resolution nuclear structural and dynamical modeling.

The intricate structural organization of the human nucleus is fundamental to cellular function and gene regulation. Recent advancements in experimental techniques, including high-throughput sequencing and microscopy, have provided valuable insights into nuclear organization. Computational modeling has played significant roles in interpreting experimental observations by reconstructing high-resolution structural ensembles and uncovering organization principles. However, the absence of standardized modeling tools poses challenges for furthering nuclear investigations. We present OpenNucleome-an open-source software designed for conducting GPU-accelerated molecular dynamics simulations of the human nucleus. OpenNucleome offers particle-based representations of chromosomes at a resolution of 100 KB, encompassing nuclear lamina, nucleoli, and speckles. This software furnishes highly accurate structural models of nuclear architecture, affording the means for dynamic simulations of condensate formation, fusion, and exploration of non-equilibrium effects. We applied OpenNucleome to uncover the mechanisms driving the emergence of 'fixed points' within the nucleus-signifying genomic loci robustly anchored in proximity to specific nuclear bodies for functional purposes. This anchoring remains resilient even amidst significant fluctuations in chromosome radial positions and nuclear shapes within individual cells. Our findings lend support to a nuclear zoning model that elucidates genome functionality. We anticipate OpenNucleome to serve as a valuable tool for nuclear investigations, streamlining mechanistic explorations and enhancing the interpretation of experimental observations.

OpenNucleome for high-resolution nuclear structural and dynamical modeling

The intricate structural organization of the human nucleus is fundamental to cellular function and gene regulation. Recent advancements in experimental techniques, including high-throughput sequencing and microscopy, have provided valuable insights into nuclear organization. Computational modeling has played significant roles in interpreting experimental observations by reconstructing high-resolution structural ensembles and uncovering organization principles. However, the absence of standardized modeling tools poses challenges for furthering nuclear investigations. We present OpenNucleome—an open-source software designed for conducting GPU-accelerated molecular dynamics simulations of the human nucleus. OpenNucleome offers particle-based representations of chromosomes at a resolution of 100 KB, encompassing nuclear lamina, nucleoli, and speckles. This software furnishes highly accurate structural models of nuclear architecture, affording the means for dynamic simulations of condensate formation, fusion, and exploration of non-equilibrium effects. We applied OpenNucleome to uncover the mechanisms driving the emergence of ‘fixed points’ within the nucleus—signifying genomic loci robustly anchored in proximity to specific nuclear bodies for functional purposes. This anchoring remains resilient even amidst significant fluctuations in chromosome radial positions and nuclear shapes within individual cells. Our findings lend support to a nuclear zoning model that elucidates genome functionality. We anticipate OpenNucleome to serve as a valuable tool for nuclear investigations, streamlining mechanistic explorations and enhancing the interpretation of experimental observations.

A novel hierarchical network-based approach to unveil the complexity of functional microbial genome

Biological networks serve a crucial role in elucidating intricate biological processes. While interspecies environmental interactions have been extensively studied, the exploration of gene interactions within species, particularly among individual microorganisms, is less developed. The increasing amount of microbiome genomic data necessitates a more nuanced analysis of microbial genome structures and functions. In this context, we introduce a complex structure using higher-order network theory, “Solid Motif Structures (SMS)”, via a hierarchical biological network analysis of genomes within the same genus, effectively linking microbial genome structure with its function. Leveraging 162 high-quality genomes of Microcystis, a key freshwater cyanobacterium within microbial ecosystems, we established a genome structure network. Employing deep learning techniques, such as adaptive graph encoder, we uncovered 27 critical functional subnetworks and their associated SMSs. Incorporating metagenomic data from seven geographically distinct lakes, we conducted an investigation into Microcystis’ functional stability under varying environmental conditions, unveiling unique functional interaction models for each lake. Our work compiles these insights into an extensive resource repository, providing novel perspectives on the functional dynamics within Microcystis. This research offers a hierarchical network analysis framework for understanding interactions between microbial genome structures and functions within the same genus.

Chromosome-level genome assembly and annotation of Flueggea virosa (Phyllanthaceae)

Flueggea virosa (Roxb. ex Willd.) Royle, an evergreen shrub and small tree in the Phyllanthaceae family, holds significant potential in garden landscaping and pharmacological applications. However, the lack of genomic data has hindered further scientific understanding of its horticultural and medicinal values. In this study, we have assembled a haplotype-resolved genome of F. virosa for the first time. The two haploid genomes, named haplotype A genome and haplotype B genome, are 487.33 Mb and 477.53 Mb in size, respectively, with contig N50 lengths of 31.45 Mb and 32.81 Mb. More than 99% of the assembled sequences were anchored to 13 pairs of pseudo-chromosomes. Furthermore, 21,587 and 21,533 protein-coding genes were predicted in haplotype A and haplotype B genomes, respectively. The availability of this chromosome-level genome fills the gap in genomic data for F. virosa and provides valuable resources for molecular studies of this species, supporting future research on speciation, functional genomics, and comparative genomics within the Phyllanthaceae family.

Multifunctional histone variants in genome function.

Histones are integral components of eukaryotic chromatin that have a pivotal role in the organization and function of the genome. The dynamic regulation of chromatin involves the incorporation of histone variants, which can dramatically alter its structural and functional properties. Contrary to an earlier view that limited individual histone variants to specific genomic functions, new insights have revealed that histone variants exert multifaceted roles involving all aspects of genome function, from governing patterns of gene expression at precise genomic loci to participating in genome replication, repair and maintenance. This conceptual change has led to a new understanding of the intricate interplay between chromatin and DNA-dependent processes and how this connection translates into normal and abnormal cellular functions.

Functional analysis of the whole CYPome and Fdxome of Streptomyces venezuelae ATCC 15439

Cytochrome P450 enzymes (CYPs or P450s) and ferredoxins (Fdxs) are ubiquitously distributed in all domains of life. Bacterial P450s are capable of catalyzing various oxidative reactions with two electrons usually donated by Fdxs. Particularly in Streptomyces, there are abundant P450s that have exhibited outstanding biosynthetic capacity of bioactive metabolites and great potential for xenobiotic metabolisms. However, no systematic study has been conducted on physiological functions of the whole cytochrome P450 complement (CYPome) and ferredoxin complement (Fdxome) of any Streptomyces strain to date, leaving a significant knowledge gap in microbial functional genomics. Herein, we functionally analyze the whole CYPome and Fdxome of Streptomyces venezuelae ATCC 15439 by investigating groups of single and sequential P450 deletion mutants, single P450 overexpression mutants, and Fdx gene deletion or repression mutants. Construction of an unprecedented P450-null mutant strain indicates that none of P450 genes are essential for S. venezuelae in maintaining its survival and normal morphology. The non-housekeeping Fdx1 and housekeeping Fdx3 not only jointly support the cellular activity of the prototypic P450 enzyme PikC, but also play significant regulatory functions. These findings significantly advance the understandings of the native functionality of P450s and Fdxs as well as their cellular interactions.

Optimization of Genotyping‐by‐Sequencing (GBS) for Germplasm Fingerprinting and Trait Mapping in Faba Bean

ABSTRACTFaba bean (Vicia faba L.) is an important pulse crop with a wide range of agroecological adaptations. The development of genomic tools and a comprehensive catalog of extant genetic diversity are crucial for developing improved faba bean cultivars. The lack of a cost‐effective genotyping platform limits the characterization of large germplasm collections, understanding of genetic diversity across populations, and implementing breeder's tools like genomic selection. Genotyping‐by‐sequencing (GBS) offers high‐resolution genotyping for both model and crop plant species, even without a reference genome sequence. The genome fragments targeted by GBS depend substantially on the restriction enzyme (RE) used for the complexity reduction step. Species with complex genomic architecture require optimization of GBS with proper RE to realize the full potential of GBS. Here, we evaluated various REs in the GBS method and identified that the combination of ApeKI/MseI proved to be the most appropriate for faba bean based on the best library quality, a high number of genomic loci spread across chromosomes, and high enrichment loci associated with the gene space. With the new optimized protocol, we constructed a genetic map using a recombinant inbred line (RIL) population and identified a QTL for seed hilum color on Chromosome 1. In addition, we also genotyped a diversity panel and performed a genome‐wide association studies (GWAS) for important agronomic traits, including plant height (PH), flowering time (FT), and number of pods per plant (PPP). We identified six SNP markers significantly associated with these traits and listed potential candidate genes. The optimized faba bean‐specific GBS procedure will facilitate access to the untapped genetic diversity for genetic research and breeding and may facilitate functional genomics.

HMGA1 orchestrates chromatin compartmentalization and sequesters genes into 3D networks coordinating senescence heterogeneity

HMGA1 is an abundant non-histone chromatin protein that has been implicated in embryonic development, cancer, and cellular senescence, but its specific role remains elusive. Here, we combine functional genomics approaches with graph theory to investigate how HMGA1 genomic deposition controls high-order chromatin networks in an oncogene-induced senescence model. While the direct role of HMGA1 in gene activation has been described previously, we find little evidence to support this. Instead, we show that the heterogeneous linear distribution of HMGA1 drives a specific 3D chromatin organization. HMGA1-dense loci form highly interactive networks, similar to, but independent of, constitutive heterochromatic loci. This, coupled with the exclusion of HMGA1-poor chromatin regions, leads to coordinated gene regulation through the repositioning of genes. In the absence of HMGA1, the whole process is largely reversed, but many regulatory interactions also emerge, amplifying the inflammatory senescence-associated secretory phenotype. Such HMGA1-mediated fine-tuning of gene expression contributes to the heterogeneous nature of senescence at the single-cell level. A similar ‘buffer’ effect of HMGA1 on inflammatory signalling is also detected in lung cancer cells. Our study reveals a mechanism through which HMGA1 modulates chromatin compartmentalization and gene regulation in senescence and beyond.

DNA methylation in mammalian development and disease.

The DNA methylation field has matured from a phase of discovery and genomic characterization to one seeking deeper functional understanding of how this modification contributes to development, ageing and disease. In particular, the past decade has seen many exciting mechanistic discoveries that have substantially expanded our appreciation for how this generic, evolutionarily ancient modification can be incorporated into robust epigenetic codes. Here, we summarize the current understanding of the distinct DNA methylation landscapes that emerge over the mammalian lifespanand discuss how they interact with other regulatory layers to support diverse genomic functions. We thenreview the rising interest in alternative patterns foundduring senescence and the somatic transition to cancer. Alongside advancements in single-cell and long-read sequencing technologies, the collective insights made across these fields offer new opportunities to connect the biochemical and genetic features of DNA methylation to cell physiology, developmental potential and phenotype.

Hmo1: A versatile member of the high mobility group box family of chromosomal architecture proteins.

Eukaryotic chromatin consisting of nucleosomes connected by linker DNA is organized into higher order structures, which is facilitated by linker histone H1. Formation of chromatin compacts and protects the genome, but also hinders DNA transactions. Cells have evolved mechanisms to modify/remodel chromatin resulting in chromatin states suitable for genome functions. The high mobility group box (HMGB) proteins are non-histone chromatin architectural factors characterized by one or more HMGB motifs that bind DNA in a sequence nonspecific fashion. They play a major role in chromatin dynamics. The Saccharomyces cerevisiae (yeast hereafter) HMGB protein Hmo1 contains two HMGB motifs. However, unlike a canonical HMGB protein that has an acidic C-terminus, Hmo1 ends with a lysine rich, basic, C-terminus, resembling linker histone H1. Hmo1 exhibits characteristics of both HMGB proteins and linker histones in its multiple functions. For instance, Hmo1 promotes transcription by RNA polymerases I and II like canonical HMGB proteins but makes chromatin more compact/stable like linker histones. Recent studies have demonstrated that Hmo1 destabilizes/disrupts nucleosome similarly as other HMGB proteins in vitro and acts to maintain a common topological architecture of genes in yeast genome. This minireview reviews the functions of Hmo1 and the underlying mechanisms, highlighting recent discoveries.

Oncogenic transcription factors instruct promoter-enhancer hubs in individual triple negative breast cancer cells.

Sequencing-based mapping of ensemble pairwise interactions among regulatory elements support the existence of topological assemblies known as promoter-enhancer hubs or cliques in cancer. Yet, prevalence, regulators, and functions of promoter-enhancer hubs in individual cancer cells remain unclear. Here, we systematically integrated functional genomics, transcription factor screening, and optical mapping of promoter-enhancer interactions to identify key promoter-enhancer hubs, examine heterogeneity of their assembly, determine their regulators, and elucidate their role in gene expression control in individual triple negative breast cancer (TNBC) cells. Optical mapping of individual SOX9 and MYC alleles revealed the existence of frequent multiway interactions among promoters and enhancers within spatial hubs. Our single-allele studies further demonstrated that lineage-determining SOX9 and signaling-dependent NOTCH1 transcription factors compact MYC and SOX9 hubs. Together, our findings suggest that promoter-enhancer hubs are dynamic and heterogeneous topological assemblies, which are controlled by oncogenic transcription factors and facilitate subtype-restricted gene expression in cancer.

Glycosphingolipid synthesis mediates immune evasion in KRAS-driven cancer.

Cancer cells frequently alter their lipids to grow and adapt to their environment1-3. Despite the critical functions of lipid metabolism in membrane physiology, signalling and energy production, how specific lipids contribute to tumorigenesis remains incompletely understood. Here, using functional genomics and lipidomic approaches, we identified de novo sphingolipid synthesis as an essential pathway for cancer immune evasion. Synthesis of sphingolipids is surprisingly dispensable for cancer cell proliferation in culture or in immunodeficient mice but required for tumour growth in multiple syngeneic models. Blocking sphingolipid production in cancer cells enhances the anti-proliferative effects of natural killer and CD8+ T cells partly via interferon-γ (IFNγ) signalling. Mechanistically, depletion of glycosphingolipids increases surface levels of IFNγ receptor subunit 1 (IFNGR1), which mediates IFNγ-induced growth arrest and pro-inflammatory signalling. Finally, pharmacological inhibition of glycosphingolipid synthesis synergizes with checkpoint blockade therapy to enhance anti-tumour immune response. Altogether, our work identifies glycosphingolipids as necessary and limiting metabolites for cancer immune evasion.

Unveiling the Secrets of Oil Palm Genetics: A Look into Omics Research.

Oil palm is a versatile oil crop with numerous applications. Significant progress has been made in applying histological techniques in oil palm research in recent years. Whole genome sequencing of oil palm has been carried out to explain the function and structure of the order genome, facilitating the development of molecular markers and the construction of genetic maps, which are crucial for studying important traits and genetic resources in oil palm. Transcriptomics provides a powerful tool for studying various aspects of plant biology, including abiotic and biotic stresses, fatty acid composition and accumulation, and sexual reproduction, while proteomics and metabolomics provide opportunities to study lipid synthesis and stress responses, regulate fatty acid composition based on different gene and metabolite levels, elucidate the physiological mechanisms in response to abiotic stresses, and explain intriguing biological processes in oil palm. This paper summarizes the current status of oil palm research from a multi-omics perspective and hopes to provide a reference for further in-depth research on oil palm.

Barcoded overexpression screens in gut Bacteroidales identify genes with roles in carbon utilization and stress resistance

A mechanistic understanding of host-microbe interactions in the gut microbiome is hindered by poorly annotated bacterial genomes. While functional genomics can generate large gene-to-phenotype datasets to accelerate functional discovery, their applications to study gut anaerobes have been limited. For instance, most gain-of-function screens of gut-derived genes have been performed in Escherichia coli and assayed in a small number of conditions. To address these challenges, we develop Barcoded Overexpression BActerial shotgun library sequencing (Boba-seq). We demonstrate the power of this approach by assaying genes from diverse gut Bacteroidales overexpressed in Bacteroides thetaiotaomicron. From hundreds of experiments, we identify new functions and phenotypes for 29 genes important for carbohydrate metabolism or tolerance to antibiotics or bile salts. Highlights include the discovery of a d-glucosamine kinase, a raffinose transporter, and several routes that increase tolerance to ceftriaxone and bile salts through lipid biosynthesis. This approach can be readily applied to develop screens in other strains and additional phenotypic assays.

Functional genomics and taxonomic insights into heavymetaltolerant novel bacterium Brevibacterium metallidurans sp. nov. NCCP-602T isolated from tannery effluent in Pakistan.

The strain designated NCCP-602T was isolated from tannery effluent, and displayed aerobic, gram-positive, rod-shaped cells that were characterized by oxidase negative, catalase positive, and non-motile features. The most favourable growth conditions were observed at a temperature of 30°C, pH 7.0, and NaCl concentration of 1% (w/v). It tolerated heavy metals at high concentrations of chromium (3600 ppm), copper (3300 ppm), cadmium (3000 ppm), arsenic (1200 ppm) and lead (1500 ppm). The results of phylogenetic analysis, derived from sequences of the 16S rRNA gene, indicated the position of strain NCCP-602T within genus Brevibacterium and showed that it was closely related to Brevibacterium ammoniilyticum JCM 17537T. Strain NCCP-602T formed a robust branch that was clearly separate from closely related taxa. A comparison of 16S rRNA gene sequence similarity and dDDH values between the closely related type strains and strain NCCP-602T provided additional evidence supporting the classification of strain NCCP-602T as a distinct novel genospecies. The polar lipid profile included diphosphatidylglycerol, glycolipid, phospholipids and amino lipids. MK-7 and MK-8 were found as the respiratory quinones, while anteiso-C15:0, iso-C15:0, iso-C16:0, iso-C17:0, and anteiso-C17:0 were identified as the predominant cellular fatty acids (> 10%). Considering the convergence of phylogenetic, phenotypic, chemotaxonomic, and genotypic traits, it is suggested that strain NCCP-602T be classified as a distinct species Brevibacterium metallidurans sp. nov. within genus Brevibacterium with type strain NCCP-602T (JCM 18882T = CGMCC1.62055T).

Anticancer Plant Secondary Metabolites Induce Linker Histone Depletion from Chromatin.

Many plant secondary metabolites (PSMs) were shown to intercalate into DNA helix or interact with DNA grooves. This may influence histone-DNA interactions changeing chromatin structure and genome functioning. Nucleosome stability and linker histone H1.2, H1.4 and H1.5 localizations were studied in HeLa cells after the treatment with 15 PSMs, which are DNA-binders and possess anticancer activity according to published data. Chromatin remodeler CBL0137 was used as a control. Effects of PSMs were studied using fluorescent microscopy, flowcytometry, quantitative reverse transcriptase-polymerase chain reaction (RT-qPCR), western-blotting. We showed that 1-hour treatment with CBL0137 strongly inhibited DNA synthesis and caused intensive linker histone depletion consistent with nucleosome destabilization. None of PSMs caused nucleosome destabilization, while most of them demonstrated significant influence on linker histone localizations. In particular, cell treatment with 11 PSMs at non-toxic concentrations induced significant translocation of the histone H1.5 to nucleoli and most of PSMs caused depletion of the histones H1.2 and H1.4 from chromatin fraction. Curcumin, resveratrol, berberine, naringenin, and quercetin caused significant redistribution of all three variants of the studied linker histones showing some overlap of PSM effects on linker histone DNA-binding. We demonstrated that PSMs, which induced the most significant redistribution of the histone H1.5 (berberine, curcumin and naringenin), influence the proportion of cells synthesizing DNA, expressing or non-expressing cyclin B and influence cell cycle distribution. Berberine induction of H1.5 translocations to nucleoli was shown to occur independently on the phases of cell cycle (metaphase was not analyzed). For the first time we revealed PSM influence on linker histone location in cell nuclei that opens a new direction of PSM research as anticancer agents.