Genomic Loci Research Articles

Multiomic characterization of RNA microenvironments by oligonucleotide-mediated proximity-interactome mapping.

RNA molecules form complex networks of molecular interactions that are central to their function and to cellular architecture. But these interaction networks are difficult to probe in situ. Here, we introduce Oligonucleotide-mediated proximity-interactome MAPping (O-MAP), a method for elucidating the biomolecules near an RNA of interest, within its native context. O-MAP uses RNA-fluorescence in situ hybridization-like oligonucleotide probes to deliver proximity-biotinylating enzymes to a target RNA in situ, enabling nearby molecules to be enriched by streptavidin pulldown. This induces exceptionally precise biotinylation that can be easily optimized and ported to new targets or sample types. Using the noncoding RNAs 47S, 7SK and Xist as models, we develop O-MAP workflows for discovering RNA-proximal proteins, transcripts and genomic loci, yielding a multiomic characterization of these RNAs' subcellular compartments and new regulatory interactions. O-MAP requires no genetic manipulation, uses exclusively off-the-shelf parts and requires orders of magnitude fewer cells than established methods, making it accessible to most laboratories.

Micro-capture-C in vessel wall cell types identifies mechanisms underlying coronary artery disease susceptibility loci

Abstract Genetics is estimated to account for 50% of Coronary Artery Disease (CAD) risk. Genome-Wide Association Studies (GWAS) have revealed 100s of genomic loci which are associated with CAD (1). However, moving from GWAS signals to causal genes has proved challenging, slowing progress in developing interventions for CAD and other diseases. The greatest challenge in interpreting GWAS data is that only 10% of variants are in protein-coding regions and 90% are in "non-coding" genomic regions. Many of these non-coding variants are found in regulatory elements such as enhancers. Enhancers are distal regulatory elements that bind transcription factors and regulate genes up to 1 million base pairs away. This is achieved by DNA looping bringing enhancers into physical proximity with target genes. However, it is not possible to predict which gene(s) a given enhancer will regulate based on genomic location alone. We aimed to develop a generalisable workflow to solve CAD GWAS signals accurately and sensitively, by identifying causal enhancer SNPs and linking them mechanistically to output genes in a cell-type specific manner. We combined cell-specific epigenomic mapping with machine learning to prioritise 21,959 SNPs from 861 candidate CAD loci, in endothelial cells (HUVECs) and smooth muscle cells (HCASMCs). Micro-Capture-C (MCC) was used to resolve enhancer-promoter interactions at base-pair resolution (2). MNase footprinting determined differential transcription factor binding at specific enhancer SNPs. In HCASMCs, 20 putative enhancers containing CAD SNPs were found to interact with 32 genes, and in HUVECs 25 enhancers interacted with 48 genes. These genes were enriched for pathways including TGFb signalling, RUNX3 activity, and MAPK6/MAPK4 signalling. To prove causality at specific loci, heterozygous SNPs were identified in the individual primary cell lines. The interaction frequency at these elements was statistically compared between risk and non-risk alleles. For example, at the 9p21.3 locus, an enhancer containing a non-risk T allele was found to interact more strongly with the CDKN2B gene, compared to the risk G allele. Interestingly, this imbalance between allele interaction was potentiated by stimulation with TGFb, a key mediator in atherosclerosis. Furthermore, the transcription factor footprint at the SNP location was altered by the risk allele, indicating the mechanistic cause of the CAD-associated change in CDKN2B expression. We demonstrate a combined computational and experimental approach to identify causal non-coding variants in CAD risk loci, at scale. Among the key outputs of this study is the identification of the risk variant (rs1537373), output gene (CDKN2B), and specific cell type (HCASMC) underlying the 9p21.3 locus association with CAD risk. Using this approach to translate GWAS signals into cell-type specific outputs genes, will lead to a better understanding of the pathogenesis of CAD, and new therapeutic opportunities.

An epigenetic editor targeting human PCSK9 efficiently and durably lowers Low Density Lipoprotein Cholesterol (LDL-C) in non-human primates

Abstract Background Reducing the risk of atherosclerotic cardiovascular disease is dependent on both the magnitude and cumulative duration of LDL-C lowering. However, in clinical practice achieving treatment goals is often hampered by low adherence to standard of care. Although PCSK9 inhibitors represent an effective option for LDL-C lowering, they require chronic life-long treatment. Epigenetic editing is a new therapeutic approach designed to durably silence genes by leveraging nature’s endogenous cellular mechanism for gene regulation via CpG methylation without the inherent genotoxic risks of genome editing approaches that nick or cut the DNA. We are developing a PCSK9 epigenetic editor (PCSK9-EE) to durably silence PCSK9, with the promise of lifelong reduction in LDL-C. Methods Here, we describe the development of an epigenetic editor targeting human PCSK9. This epigenetic editor consists of a DNA targeting component fused to a transcriptional repressor domain and a DNA methyltransferase domain. PCSK9-EE was first screened and evaluated for specificity and efficacy in primary human hepatocytes (PHH). In vivo efficacy and durability were then evaluated in a transgenic mouse carrying the human PCSK9 genomic locus (hPCSK9-Tg mouse) and non-human primates (NHP). Results Our initial screen of PCSK9-EEs identified several hits that efficiently suppressed secreted PCSK9 levels in immortalized liver cells. We confirmed that our top PCSK9-EE candidates robustly decreased PCSK9 secretion in PHH and were found to be highly specific at targeting PCSK9 as assessed by RNA-seq, methylation array, and whole-genome methylation sequencing. To examine in vivo activity, we treated hPCSK9-Tg mice with a single administration of a lipid nanoparticle delivering mRNA encoding our PCSK9-EE and observed >97% reduction in liver PCSK9 mRNA, and >98% reduction in plasma PCSK9 for the duration of the study (one year). We then delivered PCSK9-EE in NHP and observed a robust reduction in plasma PCSK9 with a concomitant reduction in LDL-C. These effects were found to be durable for at least 6 months following a single administration of PCSK9-EE in NHP. To establish a mechanistic relationship between PCSK9-EE’s molecular action and durable PCSK9 silencing, we performed serial liver biopsies in each NHP approximately 2 months apart. Robust CpG methylation at the PCSK9 genomic locus was observed in the first liver biopsy and was comparable to that observed in the second liver biopsy suggesting that PCSK9-EE induced stable CpG methylation in NHP livers. Conclusions We believe PCSK9-EE may offer an effective, safe, and durable approach to suppress hepatic PCSK9 expression and achieve long-lasting reduction in LDL-C. Furthermore, epigenetic editing may hold promise for one-and-done approaches for the treatment of atherosclerotic cardiovascular disease without cutting, nicking, or altering the DNA sequence.

Combined Analysis of Second- and Third-Generation Transcriptome Sequencing for Gene Characteristics and Identification of Key Splicing Variants in Wound Healing of Ganxi Goat Skin

Ganxi goat is a local breed of goat unique to Jiangxi Province, China, named for its primary distribution in the Ganxi region. Ganxi goats are primarily raised through grazing, showing good adaptation to the local humid and hot climate, strong disease resistance, and occupying an important position in the local livestock industry. The skin, as the main barrier of the body, plays an indispensable role in resisting the invasion of external pathogenic factors and has received increasing attention in the medical and scientific fields. In this study, Ganxi goat skin was used as the research subject. Full-length transcriptome sequencing of Ganxi goat skin was performed using PacBio third-generation sequencing technology to supplement and improve the annotation information of the Ganxi goat genome. A combined analysis of second- and third-generation transcriptome sequencing was used to analyze the splicing variant events of hub genes (CDC20, MMP2, TIMP1, and EDN1) and the expression changes in each splicing variant in skin samples on day 0 and day 5 after surgical wounding. The regulatory role of related hub gene splicing variants in wound healing was analyzed. A total of 926,667 full-length non-chimeric sequences were obtained, optimizing the annotation information of 3794 genomic gene loci and identifying 2834 new genes, 256 new LncRNAs, 12,283 alternative splicing events, 549 genes with polyadenylation, and 112 fusion genes. Three splicing variant forms were identified in both the CDC20 and EDN1 genes, seven in MMP2, and two in TIMP1. The expression levels of most splicing variants showed significant changes in the skin samples on days 0 and 5 after wounding, potentially participating in the regulation of wound healing. This study provides fundamental data for the annotation of the goat genome and offers a reference for studying the regulatory mechanisms of wound healing.

Genome-wide computational analysis of the dirigent gene family in Solanum lycopersicum

BackgroundDirigent (DIR) genes play a key role in the development of organic products in plants. They confer conformational influence on processes that lack stereoselectivity and regioselectivity through processes that are mostly understood. They are required to produce lignans, which are a unique and widely distributed family of plant secondary metabolites with intriguing pharmacological characteristics and potential role in plant development. DIR genes are implicated in the process of lignification and protect plants from environmental stresses, including biotic and abiotic stresses. Nevertheless, no research has been performed on the DIR gene family in Solanum lycopersicum. This study provides detailed information on the DIR gene family in S. lycopersicum.Methods and resultsThe conserved domain analysis, phylogenetic analysis, evolutionary adaptation, cis-acting elements, proteomic analysis, signal peptide detection, transmembrane potential analysis, sequence identity and similarity analysis, gene assembly, genomic localization, duplication of gene analysis, and evolutionary linkage of 31 potential DIR genes were studied. All these analyses provide a deep understanding of DIR genes in the S. lycopersicum genome that will provide a useful reference for further functional analysis of the DIR genes in S. lycopersicum.ConclusionThis research provides an in-depth and comprehensive explanation of the detailed process and structural characterization of DIR genes in the genome of S. lycopersicum, laying the groundwork for future plant genetic engineering and crop development exploration. This work will provide valuable information for identifying DIR genes in higher plants and support future research on the DIR gene family.

Genomic localization bias of secondary metabolite gene clusters and association with histone modifications in Aspergillus

Abstract Fungi are well-known producers of bioactive secondary metabolites (SMs), which have been exploited for decades by humankind for various medical applications like therapeutics and antibiotics. SMs are synthesized by biosynthetic gene clusters (BGCs) – physically co-localized and co-regulated genes. Because BGCs are often regulated by histone post-translational modifications (PTMs), it was suggested that their chromosomal location is important for their expression. Studies in a few fungal species indicated an enrichment of BGCs in sub-telomeric regions; however, there is no evidence that BGCs with distinct genomic localization are regulated by different histone PTMs. Here, we used 174 Aspergillus species covering 22 sections to determine the correlation between BGC genomic localization, gene expression and histone PTMs. We found a high abundance and diversity of SM backbone genes across the Aspergillus genus, with notable unique genes within sections. Being unique or conserved in many species, BGCs showed a strong bias for being localized in low-synteny regions, regardless of their position in chromosomes. Using chromosome-level assemblies, we also confirmed a significantly biased localization in sub-telomeric regions. Notably, SM backbone genes in sub-telomeric regions and about half of those in low-synteny regions exhibit higher gene expression variability, likely due to the similar higher variability in H3K4me3 and H3K36me3 histone PTMs; while variations in histone H3 acetylation and H3K9me3 are not correlated to genomic localization and expression variation, as analyzed in two Aspergillus species. Expression variability across four Aspergillus species further supports that BGCs tend to be located in low-synteny regions and that regulation of expression in those regions likely involves different histone PTMs than the most commonly studied modifications.

Codon optimality influences homeostatic gene expression in zebrafish.

The ribosome plays a crucial role in translating mRNA into protein; however, the genetic code extends beyond merely specifying amino acids. Upon translation, codons, the three-nucleotide sequences interpreted by ribosomes, have regulatory properties affecting mRNA stability, a phenomenon known as codon optimality. Codon optimality has been previously observed in vertebrates during embryogenesis, where specific codons can influence the stability and degradation rates of mRNA transcripts. In our previous work, we demonstrated that codon optimality impacts mRNA stability in human cell lines. However, the extent to which codon content influences vertebrate gene expression in vivo remained unclear. In this study, we expand on our previous findings by demonstrating that codon optimality has a robust effect on homeostatic mRNA and protein levels in whole zebrafish during normal physiological conditions. Using reporters with nearly identical nucleotide sequences but different codon compositions, all expressed from the same genomic locus, we show that codon composition can significantly influence gene expression. This study provides new insights into the regulatory roles of codon usage in vertebrate gene expression and underscores the importance of considering codon optimality in genetic and translational research. These findings have broad implications for understanding the complexities of gene regulation and could inform the design of synthetic genes and therapeutic strategies targeting mRNA stability.

Predicting protein synergistic effect in Arabidopsis using epigenome profiling

Histone modifications can regulate transcription epigenetically by marking specific genomic loci, which can be mapped using chromatin immunoprecipitation sequencing (ChIP-seq). Here we present QHistone, a predictive database of 1534 ChIP-seqs from 27 histone modifications in Arabidopsis, offering three key functionalities. Firstly, QHistone employs machine learning to predict the epigenomic profile of a query protein, characterized by its most associated histone modifications, and uses these modifications to infer the protein’s role in transcriptional regulation. Secondly, it predicts synergistic regulatory activities between two proteins by comparing their profiles. Lastly, it detects previously unexplored co-regulating protein pairs by screening all known proteins. QHistone accurately identifies histone modifications associated with specific known proteins, and allows users to computationally validate their results using gene expression data from various plant tissues. These functions demonstrate an useful approach to utilizing epigenome data for gene regulation analysis, making QHistone a valuable resource for the scientific community (https://qhistone.paoyang.ipmb.sinica.edu.tw).

Lytic promoter activity during herpes simplex virus latency is dependent on genome location.

Herpes simplex virus 1 (HSV-1) is a significant pathogen that establishes lifelong latent infections with intermittent episodes of resumed disease. In mouse models of HSV infection, sporadic low-level lytic gene expression has been detected during latency in the absence of reactivation events that lead to production of new viruses. This viral activity during latency has been reported using a sensitive Cre-marking model for several lytic gene promoters placed in one location in the HSV-1 genome. Here, we extend these findings in the same model by examining first, the activity of an ectopic lytic gene promoter in several places in the genome and second, whether any promoters might be active in their natural context. We found that Cre expression was detected during latency from ectopic and native promoters, but only in locations near the ends of the unique long genome segment. This location is significant because it is in close proximity to the region from which latency-associated transcripts (LATs) are derived. These results show that native HSV-1 lytic gene promoters can produce protein products during latency, but that this activity is only detectable when they are located close to the LAT locus.IMPORTANCEHSV is a significant human pathogen and the best studied model of mammalian virus latency. Traditionally, the active (lytic) and inactive (latent) phases of infection were considered to be distinct, but the notion of latency being entirely quiescent is evolving due to the detection of some lytic gene expression during latency. Here, we add to this literature by finding that the activity can be found for native lytic gene promoters as well as for constructs placed ectopically in the HSV genome. However, this activity was only detectable when these promoters were located close by a region known to be transcriptionally active during latency. These data have implications for our understanding of HSV gene regulation during latency and the extent to which transcriptionally active regions are insulated from adjacent parts of the viral genome.

Assessment of the Degree Of DNA Methylation in Lymphocytes аfter а Single Blood Irradiation in vitro

DNA methylation is one of the processes of epigenetic regulation of the genome, which is sensitive to the influence of endogenous and exogenous factors. The effect of ionizing radiation on the genome is accompanied by a change in the degree of DNA methylation, which can be dose-dependent and persist for a long time after radiation exposure. The objective of the study was to assess the degree of DNA methylation of blood lymphocytes after a single exposure to gamma radiation at a dose of 1.5 Gy using wide-genome bisulfite sequencing. The study included 10 conditionally healthy male employees of the ionizing radiation facility who were not exposed to radiation and did not suffer from chronic diseases. The material was whole blood: 0 Gy (control samples) and 1.5 Gy (experimental samples irradiated with gamma radiation). After irradiation with subsequent cultivation of whole blood, DNA isolation and bisulfite sequencing of limited sets of genomic loci (Reduced representation bisulfite sequencing) was performed using XmaI restriction enzyme (XmaI-RRBS). 41 genes were identified, including 26 genes (HOXD4, PADI2, FOXK1, FTCD, PRDM16, TOM1, PPP1R14A, FLNB, OR1F1, RARA, CRTAC1, AP5B1, ARL5C, NOC2L, MAMDC4, FGFRL1, PPFIA3, CUX2, ANKRD20A19P, FAM83H-AS1, CBFA2T3, POLN, MIR4458HG, FNBP1, SPIRE2, and ZSCAN10) have a tendency to hypomethylation DNA, and another 15 genes (CHRNA4, SEPTIN9, ZNF174, ELK3, NFAM1, ALG10, SOX8, KLHL30, URI1, HBZ, KLF14, MYO16, MYEOV, DMKN, and PAX7) tend to have hypermethylated state detected in at least 50٪ of the experimental samples. Thus, the genes identified in this study can be promising markers of radiation exposure and, in the future, be used to develop a new type of biological dosimetry – epigenomic dosimetry of personnel in contact with ionizing radiation sources in the course of their professional activities.

Mulea: An R package for enrichment analysis using multiple ontologies and empirical false discovery rate

Traditional gene set enrichment analyses are typically limited to a few ontologies and do not account for the interdependence of gene sets or terms, resulting in overcorrected p-values. To address these challenges, we introduce mulea, an R package offering comprehensive overrepresentation and functional enrichment analysis. mulea employs a progressive empirical false discovery rate (eFDR) method, specifically designed for interconnected biological data, to accurately identify significant terms within diverse ontologies. mulea expands beyond traditional tools by incorporating a wide range of ontologies, encompassing Gene Ontology, pathways, regulatory elements, genomic locations, and protein domains. This flexibility enables researchers to tailor enrichment analysis to their specific questions, such as identifying enriched transcriptional regulators in gene expression data or overrepresented protein domains in protein sets. To facilitate seamless analysis, mulea provides gene sets (in standardised GMT format) for 27 model organisms, covering 22 ontology types from 16 databases and various identifiers resulting in almost 900 files. Additionally, the muleaData ExperimentData Bioconductor package simplifies access to these pre-defined ontologies. Finally, mulea's architecture allows for easy integration of user-defined ontologies, or GMT files from external sources (e.g., MSigDB or Enrichr), expanding its applicability across diverse research areas. mulea is distributed as a CRAN R package downloadable from https://cran.r-project.org/web/packages/mulea/ and https://github.com/ELTEbioinformatics/mulea. It offers researchers a powerful and flexible toolkit for functional enrichment analysis, addressing limitations of traditional tools with its progressive eFDR and by supporting a variety of ontologies. Overall, mulea fosters the exploration of diverse biological questions across various model organisms.

Prophages in Vibrio.

Although tailed bacteriophages (phages) of the class Caudoviricetes are thought to constitute the most abundant and ecologically relevant group of phages that can integrate their genome into the host chromosome, it is becoming increasingly clear that other prophages are widespread. Here, we show that prophages derived from filamentous and tailless phages with genome sizes below 16 kb make up the majority of prophages in marine bacteria of the genus Vibrio. To estimate prophage prevalence unaffected by database biases, we combined comparative genomics and chemical induction of 58 diverse Vibrio cyclitrophicus isolates, resulting in 107 well-curated prophages. Complemented with computationally predicted prophages, we obtained 1,158 prophages from 931 naturally co-existing strains of the family Vibrionaceae. Prophages resembling tailless and filamentous phages predominated, accounting for 80% of all prophages in V. cyclitrophicus and 60% across the Vibrionaceae. In our experimental model, prophages of all three viral realms actively replicated upon induction indicating their ability to transfer to new hosts. Indeed, prophages were rapidly gained and lost, as suggested by variable prophage content between closely related V. cyclitrophicus. Prophages related to filamentous and tailless phages were integrated into only three genomic locations and restored the function of their integration site. Despite their small size, they contained highly diverse accessory genes that may contribute to host fitness, such as phage defense systems. We propose that, like their well-studied tailed equivalent, tailless and filamentous temperate phages are active and highly abundant drivers of host ecology and evolution in marine Vibrio, which have been largely overlooked.

Sirtuin 1 regulates the phenotype and functions of dendritic cells through Ido1 pathway in obesity

Sirtuin 1 (SIRT1) is a class III histone deacetylase (HDAC3) that plays a crucial role in regulating the activation and differentiation of dendritic cells (DCs) as well as controlling the polarization and activation of T cells. Obesity, a chronic inflammatory condition, is characterized by the activation of immune cells in various tissues. We hypothesized that SIRT1 might influence the phenotype and functions of DCs through the Ido1 pathway, ultimately leading to the polarization towards pro-inflammatory T cells in obesity. In our study, we observed that SIRT1 activity was reduced in bone marrow-derived DCs (BMDCs) from obese animals. These BMDCs exhibited elevated oxidative phosphorylation (OXPHOS) and increased extracellular acidification rates (ECAR), along with enhanced expression of class II MHC, CD86, and CD40, and elevated secretion of IL-12p40, while the production of TGF-β was reduced. The kynurenine pathway activity was decreased in BMDCs from obese animals, particularly under SIRT1 inhibition. SIRT1 positively regulated the expression of Ido1 in DCs in a PPARγ-dependent manner. To support these findings, ATAC-seq analysis revealed that BMDCs from obese mice had differentially regulated open chromatin regions compared to those from lean mice, with reduced chromatin accessibility at the Sirt1 genomic locus in BMDCs from obese WT mice. Gene Ontology (GO) enrichment analysis indicated that BMDCs from obese animals had disrupted metabolic pathways, including those related to GTPase activity and insulin response. Differential expression analysis showed reduced levels of Pparg and Sirt1 in BMDCs from obese mice, which was challenged and confirmed using BMDCs from mice with conditional knockout of Sirt1 in dendritic cells (SIRT1∆). This study highlights that SIRT1 controls the metabolism and functions of DCs through modulation of the kynurenine pathway, with significant implications for obesity-related inflammation.

Template switching during DNA replication is a prevalent source of adaptive gene amplification.

Copy number variants (CNVs)-gains and losses of genomic sequences-are an important source of genetic variation underlying rapid adaptation and genome evolution. However, despite their central role in evolution little is known about the factors that contribute to the structure, size, formation rate, and fitness effects of adaptive CNVs. Local genomic sequences are likely to be an important determinant of these properties. Whereas it is known that point mutation rates vary with genomic location and local DNA sequence features, the role of genome architecture in the formation, selection, and the resulting evolutionary dynamics of CNVs is poorly understood. Previously, we have found that the GAP1 gene in Saccharomyces cerevisiae undergoes frequent and repeated amplification and selection under long-term experimental evolution in glutamine-limiting conditions. The GAP1 gene has a unique genomic architecture consisting of two flanking long terminal repeats (LTRs) and a proximate origin of DNA replication (autonomously replicating sequence, ARS), which are likely to promote rapid GAP1 CNV formation. To test the role of these genomic elements on CNV-mediated adaptive evolution, we performed experimental evolution in glutamine-limited chemostats using engineered strains lacking either the adjacent LTRs, ARS, or all elements. Using a CNV reporter system and neural network simulation-based inference (nnSBI) we quantified the formation rate and fitness effect of CNVs for each strain. We find that although GAP1 CNVs repeatedly form and sweep to high frequency in strains with modified genome architecture, removal of local DNA elements significantly impacts the rate and fitness effect of CNVs and the rate of adaptation. We performed genome sequence analysis to define the molecular mechanisms of CNV formation for 177 CNV lineages. We find that across all four strain backgrounds, between 26% and 80% of all GAP1 CNVs are mediated by Origin Dependent Inverted Repeat Amplification (ODIRA) which results from template switching between the leading and lagging strand during DNA synthesis. In the absence of the local ARS, a distal ARS can mediate CNV formation via ODIRA. In the absence of local LTRs, homologous recombination mechanisms still mediate gene amplification following de novo insertion of retrotransposon elements at the locus. Our study demonstrates the remarkable plasticity of the genome and reveals that template switching during DNA replication is a frequent source of adaptive CNVs.

Integrating Multi-Omics Data to Uncover Prostate Tissue DNA Methylation Biomarkers and Target Genes for Prostate Cancer Risk.

Previous studies have indicated that specific CpG sites may be linked to the risk of prostate cancer (PCa) by regulating the expression of PCa target genes. However, most existing studies aim to identify DNA methylation (DNAm) biomarkers through blood tissue genetic instruments, which impedes the identification of relevant biomarkers in prostate tissue. To identify PCa risk-associated CpG sites in prostate tissue, we established genetic prediction models of DNAm levels using data from normal prostate samples in the GTEx (N = 108) and assessed associations between genetically predicted DNAm in prostate and PCa risk by studying 122,188 cases and 604,640 controls. We observed significant associations for 3879 CpG sites, including 926 at novel genomic loci. Among them, DNAm levels of 80 CpG sites located at novel loci are significantly associated with expression levels of 45 neighboring genes in normal prostate tissue. Of these genes, 11 further exhibit significant associations with PCa risk for their predicted expression levels in prostate tissue. Intriguingly, a total of 31 CpG sites demonstrate consistent association patterns across the methylation-gene expression-PCa risk pathway. Our findings suggest that specific CpG sites may be related to PCa risk by modulating the expression of nearby target genes.

Long range segmentation of prokaryotic genomes by gene age and functionality.

Bacterial and archaeal genomes encompass numerous operons that typically consist of two to five genes. On larger scales, however, gene order is poorly conserved through the evolution of prokaryotes. Nevertheless, non-random localization of different classes of genes on prokaryotic chromosomes could reflect important functional and evolutionary constraints. We explored the patterns of genomic localization of evolutionarily conserved (ancient) and variable (young) genes across the diversity of bacteria and archaea. Nearly all bacterial and archaeal chromosomes were found to encompass large segments of 100-300 kb that were significantly enriched in either ancient or young genes. Similar clustering of genes with lethal knockout phenotype (essential genes) was observed as well. Mathematical modeling of genome evolution suggests that this long-range gene clustering in prokaryotic chromosomes reflects perpetual genome rearrangement driven by a combination of selective and neutral processes rather than evolutionary conservation.

Unveiling the shared genetic architecture between testosterone and polycystic ovary syndrome

Testosterone (T) is a critical predictor of polycystic ovary syndrome (PCOS) but the genetic overlap between T and PCOS has not been established. Here by leveraging genetic datasets from large-scale genome-wide association studies, we assessed the genetic correlation and polygenic overlap between PCOS and three T-related traits using linkage disequilibrium score regression and the bivariate causal mixture model methods. The conjunctional false discovery rate (conjFDR) method was employed to identify shared causal variants. Functional annotation of variants was conducted using FUMA. Total T and bioavailable T exhibited positive correlations with PCOS, while sex hormone-binding globulin (SHBG) showed a negative correlation. All three traits demonstrated extensive genetic overlap with PCOS, with a minimum of 68% of T-related variants influencing PCOS. The conjFDR revealed 4 to 6 causal variants within joint genomic loci shared between PCOS and T-related traits. Functional annotations suggested that these variants might impact PCOS by modulating nearby genes, such as FSHB. Our findings support the hypothesis that PCOS is significantly influenced by androgen abnormalities. Additionally, this study identified several causal variants potentially involved in shared biological mechanisms between PCOS and T regulation.

Reporter Alleles in hiPSCs: Visual Cues on Development and Disease.

Reporter alleles are essential for advancing research with human induced pluripotent stem cells (hiPSCs), notably in developmental biology and disease modeling. This study investigates the state-of-the-art gene-editing techniques tailored for generating reporter alleles in hiPSCs, emphasizing their effectiveness in investigating cellular dynamics and disease mechanisms. Various methodologies, including the application of CRISPR/Cas9 technology, are discussed for accurately integrating reporter genes into the specific genomic loci. The synthesis of findings from the studies utilizing these reporter alleles reveals insights into developmental processes, genetic disorder modeling, and therapeutic screening, consolidating the existing knowledge. These hiPSC-derived models demonstrate remarkable versatility in replicating human diseases and evaluating drug efficacy, thereby accelerating translational research. Furthermore, this review addresses challenges and future directions in refining the reporter allele design and application to bolster their reliability and relevance in biomedical research. Overall, this investigation offers a comprehensive perspective on the methodologies, applications, and implications of reporter alleles in hiPSC-based studies, underscoring their essential role in advancing both fundamental scientific understanding and clinical practice.

Dnajc3 (HSP40) Modulates Axon Regeneration in the Mouse Optic Nerve.

The present study is designed to identify the genes modulating optic nerve regeneration in the mouse. Using the BXD mouse strains as a genetic mapping panel, we examined differential responses to axon regeneration in order to map genomic loci modulating axonal regeneration. To study regeneration in the optic nerve, Pten was knocked down in the retinal ganglion cells using adeno-associated virus (AAV) delivery of an shRNA, followed by the induction of a mild inflammatory response by an intravitreal injection of Zymosan with CPT-cAMP. The axons of the retinal ganglion cells were damaged by optic nerve crush (ONC). Following a 12-day survival period, regenerating axons were labeled by Cholera Toxin B. Two days later, the regenerating axons within the optic nerve were examined to determine the number of regenerating axons and the distance traveled down the optic nerve. An integral genomic map was made using the regenerative response. Candidate genes were tested by knocking down expression using shRNA or by overexpressing the gene in AAV vectors. The analysis revealed a considerable amount of differential axonal regeneration across all 33 BXD strains, demonstrated by the number of axons regenerating and the length of the regenerating axons. Some strains (BXD99, BXD90, and BXD29) demonstrated significant axonal regeneration; while other strains (BXD13, BXD18, and BXD34) had very little axon regrowth. Within the regenerative data, there was a 4-fold increase in distance regenerated and a 7.5-fold difference in the number of regenerating axons. These data were used to map a quantitative trait locus modulating axonal regeneration to Chromosome 14 (115 to 119 Mb). Within this locus were 16 annotated genes. Subsequent testing revealed that one candidate gene, Dnajc3 , modulates axonal regeneration. Knocking down of Dnajc3 led to a decreased regeneration response in the high regenerative strains (BXD90), while overexpression of Dnajc3 resulted in an increased regeneration response in C57BL/6J and a low regenerative strain (BXD34). In this study, Dnajc3 (encodes Heat Shock Protein 40, HSP40, a molecular chaperone) was identified as a modulator of axon regeneration in mice. This is the first report defining the role of Dnajc3 (HSP40) in axon regeneration.

Regulatory genome annotation of 33 insect species

Annotation of newly sequenced genomes frequently includes genes, but rarely covers important non-coding genomic features such as the cis-regulatory modules—e.g., enhancers and silencers—that regulate gene expression. Here, we begin to remedy this situation by developing a workflow for rapid initial annotation of insect regulatory sequences, and provide a searchable database resource with enhancer predictions for 33 genomes. Using our previously developed SCRMshaw computational enhancer prediction method, we predict over 2.8 million regulatory sequences along with the tissues where they are expected to be active, in a set of insect species ranging over 360 million years of evolution. Extensive analysis and validation of the data provides several lines of evidence suggesting that we achieve a high true-positive rate for enhancer prediction. One, we show that our predictions target specific loci, rather than random genomic locations. Two, we predict enhancers in orthologous loci across a diverged set of species to a significantly higher degree than random expectation would allow. Three, we demonstrate that our predictions are highly enriched for regions of accessible chromatin. Four, we achieve a validation rate in excess of 70% using in vivo reporter gene assays. As we continue to annotate both new tissues and new species, our regulatory annotation resource will provide a rich source of data for the research community and will have utility for both small-scale (single gene, single species) and large-scale (many genes, many species) studies of gene regulation. In particular, the ability to search for functionally related regulatory elements in orthologous loci should greatly facilitate studies of enhancer evolution even among distantly related species.