Regulatory Regions Research Articles

Generation of humanized mouse models to support therapeutic development for SYNGAP1 and STXBP1 disorders.

Heterozygous variants in SYNGAP1 and STXBP1 lead to distinct neurodevelopmental disorders caused by haploinsufficient levels of post-synaptic SYNGAP1 and pre-synaptic STXBP1, which are critical for normal synaptic function. While several gene-targeted therapeutic approaches have proven efficacious in vitro, these often target regions of the human gene that are not conserved in rodents, hindering the pre-clinical development of these compounds and their transition to the clinic. To overcome this limitation, here we generate and characterize Syngap1 and Stxbp1 humanized mouse models in which we replaced the mouse Syngap1 and Stxbp1 gene, respectively, with the human counterpart, including regulatory and non-coding regions. Fully humanized Syngap1 mice present normal viability and can be successfully crossed with currently available Syngap1 haploinsufficiency mouse models to generate Syngap1 humanized haploinsufficient mice. Stxbp1 mice were successfully humanized, yet exhibit impaired viability (particularly males) and reduced STXBP1 protein abundance. Mouse viability could be improved by outcrossing this model to other mouse strains, while Stxbp1 humanized females and hybrid mice can be used to evaluate target engagement of human-specific therapeutics. Overall, these humanized mouse models represent a broadly available tool to further pre-clinical therapeutic development for SYNGAP1 and STXBP1 disorders.

Variant-to-function mapping of late-onset Alzheimer's disease GWAS signals in human microglial cell models implicates RTFDC1 at the CASS4 locus.

Late-onset Alzheimer's disease (LOAD) research has principally focused on neurons over the years due to their known role in the production of amyloid beta plaques and neurofibrillary tangles. In contrast, recent genomic studies of LOAD have implicated microglia as culprits of the prolonged inflammation exacerbating the neurodegeneration observed in patient brains. Indeed, recent LOAD genome-wide association studies (GWAS) have reported multiple loci near genes related to microglial function, including TREM2, ABI3, and CR1. However, GWAS alone cannot pinpoint underlying causal variants or effector genes at such loci, as most signals reside in non-coding regions of the genome and could presumably confer their influence frequently via long-range regulatory interactions. We elected to carry out a combination of ATAC-seq and high-resolution promoter-focused Capture-C in two human microglial cell models (iPSC-derived microglia and HMC3) in order to physically map interactions between LOAD GWAS-implicated candidate causal variants and their corresponding putative effector genes. Notably, we observed consistent evidence that rs6024870 at the GWAS CASS4 locus contacted the promoter of nearby gene, RTFDC1. We subsequently observed a directionallly consistent decrease in RTFDC1 expression with the the protective minor A allele of rs6024870 via both luciferase assays in HMC3 cells and expression studies in primary human microglia. Through CRISPR-Cas9-mediated deletion of the putative regulatory region harboring rs6024870 in HMC3 cells, we observed increased pro-inflammatory cytokine secretion and decreased DNA double strand break repair related, at least in part, to RTFDC1 expression levels. Our variant-to-function approach therefore reveals that the rs6024870-harboring regulatory element at the LOAD 'CASS4' GWAS locus influences both microglial inflammatory capacity and DNA damage resolution, along with cumulative evidence implicating RTFDC1 as a novel candidate effector gene.

A conserved germline-specific Dsn1 alternative splice isoform supports oocyte and embryo development

The chromosome segregation and cell division programs associated with somatic mitosis and germline meiosis display dramatic differences such as kinetochore orientation, cohesin removal, or the presence of a gap phase.1,2,3,4,5,6 These changes in chromosome segregation require alterations to the established cell division machinery.5,6 It remains unclear what aspects of kinetochore function and its regulatory control differ between the mitotic and meiotic cell divisions to rewire these core processes. Alternative RNA splicing can generate distinct protein isoforms to allow for the differential control of cell processes across cell types. However, alternative splice isoforms that differentially modulate distinct cell division programs have remained elusive. Here, we demonstrate that mammalian germ cells express an alternative mRNA splice isoform for the kinetochore component, DSN1, a subunit of the MIS12 complex that links the centromeres to spindle microtubules during chromosome segregation. This germline DSN1 isoform bypasses the requirement for Aurora kinase phosphorylation for its centromere localization due to the absence of a key regulatory region allowing DSN1 to display persistent centromere localization. Expression of the germline DSN1 isoform in somatic cells results in constitutive kinetochore localization, chromosome segregation errors, and growth defects, providing an explanation for its tight cell-type-specific expression. Reciprocally, precisely eliminating expression of the germline-specific DSN1 splice isoform in mouse models disrupts oocyte maturation and early embryonic divisions coupled with a reduction in fertility. Together, this work identifies a germline-specific splice isoform for a chromosome segregation component and implicates its role in mammalian fertility.

Unveiling Genetic Potential for Equine Meat Production: A Bioinformatics Approach.

In view of the predicted significant increase in global meat production, alternative sources such as horsemeat are becoming increasingly important due to their lower environmental impact and high nutritional value. This study aimed to identify SNP markers on the GeneSeek® Genomic Profiler™ Equine (Neogen, Lansing, MI, USA) that are important for horsemeat production traits. First, orthologous genes related to meat yield in cattle and common genes between horses and cattle within QTLs for body size and weight were identified. Markers for these genes were then evaluated based on predicted variant consequences, GERP scores, and positions within constrained elements and orthologous regulatory regions in pigs. A total of 268 markers in 57 genes related to meat production were analyzed. This resulted in 27 prioritized SNP markers in 22 genes, including notable markers in LCORL, LASP1, IGF1R, and MSTN. These results will benefit smallholder farmers by providing genetic insights for selective breeding that could improve meat yield. This study also supports future large-scale genetic analyses such as GWAS and Genomic Best Linear Unbiased Prediction (GBLUP). The results of this study may be helpful in improving the accuracy of genomic breeding values. However, limitations include reliance on bioinformatics without experimental validation. Future research can validate these markers and consider a wider range of traits to ensure accuracy in equine breeding.

Using Functional Genomic Data in Monocytes/Macrophages and Genotyping to Nominate Disease-Driving Single Nucleotide Polymorphisms and Target Genes in Juvenile Idiopathic Arthritis.

GWAS have identified multiple regions that confer risk for juvenile idiopathic arthritis (JIA). However, identifying the single nucleotide polymorphisms (SNPs) that drive disease risk is impeded by the SNPs' that identify risk loci being in linkage disequilibrium (LD) with hundreds of other SNPs. Since the causal SNPs remain unknown, it is difficult to identify target genes and use genetic information to inform patient care. We used genotyping and functional data in primary human monocytes/macrophages to nominate disease-driving SNPs on JIA risk haplotypes and identify their likely target genes. We identified JIA risk haplotypes using Immunochip data from Hinks et al (Nature Gen 2013) and the meta-analysis from McIntosh et al (Arthritis Rheum 2017). We used genotyping data from 3,939 children with JIA and 14,412 healthy controls to identify SNPs that: (1) were situated within open chromatin in multiple immune cell types and (2) were more common in children with JIA than the controls (p< 0.05). We intersected the chosen SNPs (n=846) with regions of bi-directional transcription initiation characteristic of non-coding regulatory regions detected using dREG to analyze GRO-seq data. Finally, we used MicroC data to identify gene promoters interacting with the regulatory regions harboring the candidate causal SNPs. We identified 190 SNPs that overlap with dREG peaks in monocytes and126 SNPs that overlap with dREG peaks in macrophages. Of these SNPs, 101 were situated within dREG peaks in both monocytes and macrophages, suggesting that these SNPs exert their effects independent of the cellular activation state. MicroC data in monocytes identified 20 genes/transcripts whose promoters interact with the enhancers harboring the SNPs of interest. SNPs in JIA risk regions that are candidate causal variants can be further screened using functional data such as GRO-seq. This process identifies a finite number of candidate causal SNPs, the majority of which are likely to exert their biological effects independent of cellular activation state in monocytes. Three-dimensional chromatin data generated with MicroC identifies genes likely to be influenced by these SNPs. These studies demonstrate the importance of investigations into the role of innate immunity in JIA.

Mutations of CYP1B1 and FOXC1 genes for childhood glaucoma in Japanese individuals.

To explore the frequency and positions of genetic mutations in CYP1B1 and FOXC1 in a Japanese population. Molecular genetic analysis. Genomic DNA was extracted from 31 Japanese patients with childhood glaucoma (CG) from 29 families. We examined the CYP1B, FOXC1, and MYOC genes using Sanger sequencing and whole-exome sequencing (WES). For CYP1B1, we identified 9 families that harbored novel mutations, p.A202T, p.D274E, p.Q340*, and p.V420G; the remaining mutations had been previously reported. When mapped to the CYP1B1 protein structure, all mutations appeared to influence the enzymatic activity of CYP1B1 by provoking structural deformity. Five patients were homozygotes or compound heterozygotes, supporting the recessive inheritance of the CYP1B1 mutations in CG. In contrast, four patients were heterozygous for the CYP1B1 mutation, suggesting the presence of regulatory region mutations or strong modifiers. For the FOXC1 gene, we identified 3 novel mutations, p.Q23fs, p.Q70R, and p.E163*, all of which were identified in a heterozygous state. No mutation was found in the MYOC gene in these CG patients. All individuals with CYP1B1 and FOXC1 mutations were severely affected by early-onset CG. In the CYP1B1-, FOXC1-, and MYOC-negative families, we also searched for variants in the other candidate genes reported for CG through WES, but could not find any mutations in these genes. Our analyses of 29 CG families revealed 9 families with point mutations in the CYP1B1 gene, and four of those patients appeared to be heterozygotes, suggesting the presence of complex pathogenic mechanisms. FOXC1 appears to be another major causal gene of CG, indicating that panel sequencing of CYP1B1 and FOXC1 will be useful for diagnosis of CG in Japanese individuals.

Functional characterization of variants in the regulatory regions of PCSK9 gene

Functional characterization of variants in the regulatory regions of PCSK9 gene

A multi-layered integrative analysis reveals a cholesterol metabolic program in outer radial glia with implications for human brain evolution.

The definition of molecular and cellular mechanisms contributing to brain ontogenetic trajectories is essential to investigate the evolution of our species. Yet their functional dissection at an appropriate level of granularity remains challenging. Capitalizing on recent efforts that have extensively profiled neural stem cells from the developing human cortex, we develop an integrative computational framework to perform trajectory inference and gene regulatory network reconstruction, (pseudo)time-informed non-negative matrix factorization for learning the dynamics of gene expression programs, and paleogenomic analysis for a higher-resolution mapping of derived regulatory variants in our species in comparison with our closest relatives. We provide evidence for cell type-specific regulation of gene expression programs during indirect neurogenesis. In particular, our analysis uncovers a key role for a cholesterol program in outer radial glia, regulated by zinc-finger transcription factor KLF6. A cartography of the regulatory landscape impacted by Homo sapiens-derived variants reveals signals of selection clustering around regulatory regions associated with GLI3, a well-known regulator of radial glial cell cycle, and impacting KLF6 regulation. Our study contributes to the evidence of significant changes in metabolic pathways in recent human brain evolution.

3D genomic features across >50 diverse cell types reveal insights into the genomic architecture of childhood obesity.

The prevalence of childhood obesity is increasing worldwide, along with the associated common comorbidities of type 2 diabetes and cardiovascular disease in later life. Motivated by evidence for a strong genetic component, our prior genome-wide association study (GWAS) efforts for childhood obesity revealed 19 independent signals for the trait; however, the mechanism of action of these loci remains to be elucidated. To molecularly characterize these childhood obesity loci we sought to determine the underlying causal variants and the corresponding effector genes within diverse cellular contexts. Integrating childhood obesity GWAS summary statistics with our existing 3D genomic datasets for 57 human cell types, consisting of high-resolution promoter-focused Capture-C/Hi-C, ATAC-seq, and RNA-seq, we applied stratified LD score regression and calculated the proportion of genome-wide SNP heritability attributable to cell type-specific features, revealing pancreatic alpha cell enrichment as the most statistically significant. Subsequent chromatin contact-based fine-mapping was carried out for genome-wide significant childhood obesity loci and their linkage disequilibrium proxies to implicate effector genes, yielded the most abundant number of candidate variants and target genes at the BDNF, ADCY3, TMEM18 and FTO loci in skeletal muscle myotubes and the pancreatic beta-cell line, EndoC-BH1. One novel implicated effector gene, ALKAL2 - an inflammation-responsive gene in nerve nociceptors - was observed at the key TMEM18 locus across multiple immune cell types. Interestingly, this observation was also supported through colocalization analysis using expression quantitative trait loci (eQTL) derived from the Genotype-Tissue Expression (GTEx) dataset, supporting an inflammatory and neurologic component to the pathogenesis of childhood obesity. Our comprehensive appraisal of 3D genomic datasets generated in a myriad of different cell types provides genomic insights into pediatric obesity pathogenesis.

Hsa-miR-874-3p Reduces Endogenous Expression of RGS4-1 Isoform In Vitro.

The level of the regulator of G-protein signaling 4-1 (RGS4-1) isoform, the longest RGS4 isoform, is significantly reduced in the dorsolateral prefrontal cortex (DLPFC) of people with schizophrenia. However, the mechanism behind this has not been clarified. The 3'untranslated regions (3'UTRs) are known to regulate the levels of their mRNA splice variants. We constructed recombinant pmir-GLO vectors with a truncated 3' regulatory region of the RGS4 gene (3R1, 3R2, 3R3, 3R4, 3R5, and 3R6). The dual-luciferase reporter assay was conducted to find functional regions in HEK-293, SK-N-SH, and U87cells and then predicted miRNA binding to these regions. We performed a dual-luciferase reporter assay and a Western blot analysis after transiently transfecting the predicted miRNAs. The dual-luciferase reporter assay found that regions +401-+789, +789-+1152, and +1562-+1990 (with the last base of the termination codon being +1) might be functional regions. Hsa-miR-874-3p, associated with many psychiatric disorders, might target the +789-+1152 region in the 3'UTR of the RGS4 gene. In the dual-luciferase reporter assay, the hsa-miR-874-3p mimic, co-transfected with 3R1, down-regulated the relative fluorescence intensities. However, this was reversed when the hsa-miR-874-3p mimic was co-transfected with m3R1 (deletion of +853-+859). The hsa-miR-874-3p mimic significantly decreased the endogenous expression of the RGS4-1 isoform in HEK-293 cells. Hsa-miR-874-3p inhibits the expression of the RGS4-1 isoform by targeting +853-+859.

Multimodal analysis of methylation and fragmentomic profiles in plasma cell free DNA for differentiation of benign and malignant breast tumors.

7 Background: Breast cancer ranks as the second leading cause of cancer-related mortality among women globally. Early detection of breast cancer is crucial for improving patient outcomes and reducing mortality rates. Liquid biopsy, which relies on circulating tumor DNA (ctDNA) shed by breast tumors into the bloodstream, presents a promising non-invasive approach for early breast cancer detection. However, accurately distinguishing between benign breast abnormalities and malignant tumors remains a significant clinical challenge, as misdiagnosis can lead to unnecessary invasive procedures. Methods: Herein, we employed a multimodal analysis approach, namely SPOT-MAS (Screen for the Presence of Tumor by DNA Methylation and Size), to profile alterations in methylation and fragment length patterns of cell-free DNA (cfDNA) from 169 breast cancer-confirmed patients and 99 patients diagnosed with benign breast lumps including cysts, fibroadenomas, and fibrocystic changes. A robust machine learning model was constructed using these signatures to differentiate between breast cancer patients and individuals with benign lesions. Results: Our genome-wide analyses identified distinct profiles of methylation changes, copy number alterations, and end motifs in cfDNA, enabling discrimination between breast cancer patients and individuals with benign conditions. Notably, we observed a significant enrichment of Thymine and Adenine at cfDNA cleavage sites in breast cancer patients. Moreover, our target sequencing analyses uncovered distinct methylation patterns in the regulatory regions of multiple genes, including hypermethylation in GPR126 or hypomethylation in TOP1 or MAFB in breast cancer cfDNA. Our multi-featured model achieved an AUC of 0.92 (95% CI: 0.88–0.97), a specificity of 97.44% and sensitivities of 57.14% and 61.70% for stage I and stage II–III patients, respectively. Furthermore, our multimodal assay effectively differentiated multiple molecular subtypes of breast tumors from benign lesions, achieving the highest sensitivity of 71.43% with Luminal A, followed by Luminal B, Triple Negative Breast Cancer (TNBC), Luminal B-HER2, and HER2 with sensitivities of 66.67%, 60%, 58.33% and 55.56%, respectively. Conclusions: Our findings demonstrate the potential of cancer-specific methylation and fragmentomic patterns in plasma cfDNA as novel biomarkers for accurately discriminating between breast cancer and benign lesions. This capability reduces the false-positive rate and helps avoid unnecessary biopsies in current clinical practice.

Analysis of Bacterial Community Characteristics in Maize Root Zones Under Maize-soybean Compound Planting Mode

Maize-soybean compound intercropping has the potential to increase yield and is being tested for spreading in Huang-Huai-hai Plain. However, the main regulatory regions of this cropping pattern on soil microbial communities have not been clarified. In the present study, the tested samples were collected from three maize root zones of bulk soil, rhizosphere soil, and roots under mono- and intercropping planting modes, respectively. The non-rhizosphere soil chemical properties and enzyme activities were determined, and bacterial communities were characterized using high-throughput sequencing of the 16S rRNA gene V3-V4 region. Compared with monocropping, the maize bulk soil electric conductivity （EC）, soil organic matter （SOM）, available potassium （AK）, available phosphorus （AP）, total nitrogen （TN）, and enzyme activities of intercropping were significantly increased. The α diversities and β diversity of the bacterial community in rhizosphere soil were significantly different between the two planting modes. There were 11 bacteria genera with significantly higher abundance in the rhizosphere soil of compound planting than that of monoculture, and TN, AP, and catalase were the three most important factors contributing to their distribution. The abundances of 8 genera among the 11 genera mentioned above, unclassified Vicinamibacterales, unclassified Geminicoccaceae, MND1, unclassified Gemmatimonadaceae, Acidibacter, unclassified Vicinamibacteraceae, Sphingomonas, and unclassified Comamonadaceae were significantly positively correlated with TN. As for the bacteria distribution in maize root, AK contributed the most and had a significantly negative correlation with unclassified Rhizobiaceae and unclassified Microscillaceae and a positive correlation with Haliangium. Maize-soybean compound intercropping affected mainly the bacterial community of maize rhizosphere and had an evident effect on soil fertilizer cultivation and microbial diversity regulation, which provides a theoretical basis and practical guidance for rational intercropping to maintain agroecosystem biodiversity.

Transcriptional and epigenetic characterization of a new in vitro platform to model the formation of human pharyngeal endoderm

BackgroundThe Pharyngeal Endoderm (PE) is an extremely relevant developmental tissue, serving as the progenitor for the esophagus, parathyroids, thyroids, lungs, and thymus. While several studies have highlighted the importance of PE cells, a detailed transcriptional and epigenetic characterization of this important developmental stage is still missing, especially in humans, due to technical and ethical constraints pertaining to its early formation.ResultsHere we fill this knowledge gap by developing an in vitro protocol for the derivation of PE-like cells from human Embryonic Stem Cells (hESCs) and by providing an integrated multi-omics characterization. Our PE-like cells robustly express PE markers and are transcriptionally homogenous and similar to in vivo mouse PE cells. In addition, we define their epigenetic landscape and dynamic changes in response to Retinoic Acid by combining ATAC-Seq and ChIP-Seq of histone modifications. The integration of multiple high-throughput datasets leads to the identification of new putative regulatory regions and to the inference of a Retinoic Acid-centered transcription factor network orchestrating the development of PE-like cells.ConclusionsBy combining hESCs differentiation with computational genomics, our work reveals the epigenetic dynamics that occur during human PE differentiation, providing a solid resource and foundation for research focused on the development of PE derivatives and the modeling of their developmental defects in genetic syndromes.

DNA Methylation Effects on Van der Woude Syndrome Phenotypic Variability.

Van der Woude Syndrome (VWS) presents with combinations of lip pits (LP) and cleft lip and/or cleft palate (CL/P, CPO). VWS phenotypic heterogeneity even amongst relatives, suggests that epigenetic factors may act as modifiers. IRF6, causal for 70% of VWS cases, and TP63 interact in a regulatory loop coordinating epithelial proliferation and differentiation in palatogenesis. We hypothesize that differential DNA methylation within IRF6 and TP63 regulatory regions underlie VWS phenotypic discordance. DNA methylation of CpG sites in IRF6 and TP63 promoters and in an IRF6 enhancer element was compared amongst blood or saliva DNA samples of 78 unrelated cases. Analyses were done separately for blood and saliva, within each sex and in combination, and to address cleft type (CL/P ± LP vs. CPO ± LP) and phenotypic severity (any cleft + LP vs. any cleft only). For cleft type, blood samples showed higher IRF6 and TP63 promoter methylation on males with CPO ± LP compared to CL/P ± LP and on individuals with CPO ± LP compared to those with CL/P ± LP, respectively. Saliva samples showed higher IRF6 enhancer methylation on individuals with CPO ± LP compared to CL/P ± LP and contrary to above, lower TP63 promoter methylation on CPO ± LP compared to CL/P ± LP. For phenotypic severity, blood samples showed no differences; however, saliva samples showed higher IRF6 promoter methylation in individuals with any cleft + LP compared to those without lip pits. We observed differential methylation in IRF6 and TP63 regulatory regions associated with cleft type and phenotypic severity, indicating that epigenetic changes in IRF6 and TP63 can contribute to phenotypic heterogeneity in VWS.

DNA hypo-methylation and expression of GBP4 induces T cell exhaustion in pancreatic cancer

Immunotherapy for pancreatic ductal carcinoma (PDAC) remains disappointing due to the repressive tumor microenvironment and T cell exhaustion, in which the roles of interferon-stimulated genes were largely unknown. Here, we focused on a typical interferon-stimulated gene, GBP4, and investigated its potential diagnostic and therapeutic value in pancreatic cancer. Expression analysis on both local samples and public databases indicated that GBP4 was one of the most dominant GBP family members present in the PDAC microenvironment, and the expression level of GBP4 was negatively associated with patient survival. We then identified DNA hypo-methylation in regulatory regions of GBP4 in PDAC, and validated its regulatory role on GBP4 expression via performing targeted methylation using dCas9-SunTag-DNMAT3A-sgRNA-targeted methylation system on selected DNA locus. After that, we investigated the downstream functions of GBP4, and chemotaxis assays indicated that GBP4 overexpression significantly improved the infiltration of CD8+T cells, but also induced upregulation of immune checkpoint genes and T cell exhaustion. Lastly, in vitro T cell killing assays using primary organoids suggested that the PDAC samples with high level of GBP4 expression displayed significantly higher sensitivity to anti-PD-1 treatment. Taken together, our studies revealed the expression patterns and epigenetic regulatory mechanisms of GBP4 in pancreatic cancer and clarified the effects of GBP4 on T cell exhaustion and antitumor immunology.

Whole-exome and whole-genome sequencing of 1064 individuals with type 1 diabetes reveals novel genes for diabetic kidney disease

Aims/hypothesisDiabetic kidney disease (DKD) is a severe diabetic complication that affects one third of individuals with type 1 diabetes. Although several genes and common variants have been shown to be associated with DKD, much of the predicted inheritance remains unexplained. Here, we performed next-generation sequencing to assess whether low-frequency variants, extending to a minor allele frequency (MAF) ≤10% (single or aggregated) contribute to the missing heritability in DKD.MethodsWe performed whole-exome sequencing (WES) of 498 individuals and whole-genome sequencing (WGS) of 599 individuals with type 1 diabetes. After quality control, next-generation sequencing data were available for a total of 1064 individuals, of whom 541 had developed either severe albuminuria or end-stage kidney disease, and 523 had retained normal albumin excretion despite a long duration of type 1 diabetes. Single-variant and gene-aggregate tests for protein-altering variants (PAV) and protein-truncating variants (PTV) were performed separately for WES and WGS data and combined in a meta-analysis. We also performed genome-wide aggregate analyses on genomic windows (sliding window), promoters and enhancers using the WGS dataset.ResultsIn the single-variant meta-analysis, no variant reached genome-wide significance, but a suggestively associated common THAP7 rs369250 variant (p=1.50 × 10−5, MAF=49%) was replicated in the FinnGen general population genome-wide association study (GWAS) data for chronic kidney disease and DKD phenotypes. The gene-aggregate meta-analysis provided suggestive evidence (p<4.0 × 10−4) at four genes for DKD, of which NAT16 (MAFPAV≤10%) and LTA (also known as TNFβ, MAFPAV≤5%) are replicated in the FinnGen general population GWAS data. The LTA rs2229092 C allele was associated with significantly lower TNFR1, TNFR2 and TNFR3 serum levels in a subset of FinnDiane participants. Of the intergenic regions suggestively associated with DKD, the enhancer on chromosome 18q12.3 (p=3.94 × 10−5, MAFvariants≤5%) showed interaction with the METTL4 gene; the lead variant was replicated, and predicted to alter binding of the MafB transcription factor.Conclusions/interpretationOur sequencing-based meta-analysis revealed multiple genes, variants and regulatory regions that were suggestively associated with DKD. However, as no variant or gene reached genome-wide significance, further studies are needed to validate the findings.Graphical

Regulation of Electron Delocalization Region in 2D Heteroligand-Based Copper-Organic Framework to Enhance Charge Storage.

The rechargeable aqueous ammonium ion battery shows great potential in low-cost energy storage system because of its long life and environmental friendliness. However, most inorganic host materials used in ammonium ion batteries are still limited by slow diffusion kinetics. Herein, it is identified that a 2D heteroligand-based copper-organic framework featuring numerous ammonium ion adsorption site in the π-conjugated periodic skeleton supplies multiple accessible redox-active sites for high-performance ammonium storage. Benefitting from the effective regulation of electron delocalization by heteroligand and the inherent hydrogen bond cage mechanism between ammonium ions, the resultant full battery delivers a large specific energy density of 211.84Whkg-1, and it can be stably operated for 12000 cycles at 5Ag-1 for over 80 days. This explanatory understanding provides a new idea for the rational design of high-performance MOF-based ammonium ion battery cathode materials for efficient energy storage and conversion in the future.

Ki67 deficiency impedes chromatin accessibility and BCR gene rearrangement.

The proliferation marker Ki67 has been attributed critical functions in maintaining mitotic chromosome morphology and heterochromatin organization during the cell cycle, indicating a potential role in developmental processes requiring rigid cell-cycle control. Here, we discovered that despite normal fecundity and organogenesis, germline deficiency in Ki67 resulted in substantial defects specifically in peripheral B and T lymphocytes. This was not due to impaired cell proliferation but rather to early lymphopoiesis at specific stages where antigen-receptor gene rearrangements occurred. We identified that Ki67 was required for normal global chromatin accessibility involving regulatory regions of genes critical for checkpoint stages in B cell lymphopoiesis. In line with this, mRNA expression of Rag1 was diminished and gene rearrangement was less efficient in the absence of Ki67. Transgenes encoding productively rearranged immunoglobulin heavy and light chains complemented Ki67 deficiency, completely rescuing early B cell development. Collectively, these results identify a unique contribution from Ki67 to somatic antigen-receptor gene rearrangement during lymphopoiesis.

MSI2 regulates NLK-mediated EMT and PI3K/AKT/mTOR pathway to promote pancreatic cancer progression

BackgroundThe incidence of pancreatic cancer is increasing by years, and the 5-year survival rate is very low. Our team have revealed that Musashi2 (MSI2) could promote aggressive behaviors in pancreatic cancer by downregulating Numb and p53. MSI2 also facilitates EMT in pancreatic cancer induced by EGF through the ZEB1-ERK/MAPK signaling pathway. This study aims to further explore the molecular mechanisms of MSI2-regulated downstream pathways in pancreatic cancer.MethodsIn vitro and in vivo experiments were conducted to investigate the role and mechanism of MSI2 in promoting malignant behaviors of pancreatic cancer through regulation of NLK.ResultsGenes closely related to MSI2 were screened from the GEPIA and TCGA databases. We found that NLK showed the most significant changes in mRNA levels with consistent changes following MSI2 interference and overexpression. The high correlation between MSI2 and NLK was also observed at the protein level. Multivariate analysis revealed that both MSI2 and NLK were independent adverse indicators of survival in pancreatic cancer patients, as well as join together. In vitro, silencing or overexpressing NLK altered cell invasion and migration, by regulating EMT and the PI3K-AKT-mTOR pathway. Silencing MSI2 reduced protein expression in the EMT and PI3K-AKT-mTOR pathways, leading to decreased cell invasion and migration abilities, while these effects could be reversed by overexpression of NLK. In vivo, MSI2 silencing inhibited liver metastasis, which could be reversed by overexpressing NLK. Mechanistically, MSI2 directly binds to the translation regulatory region of NLK mRNA at positions 79–87 nt, enhancing its transcriptional activity and exerting post-transcriptional regulatory roles. The analysis of molecular docking showed the close relationship between MSI2 and NLK in pancreatic cancer patients.ConclusionsOur findings elucidate the regulatory mechanisms of the MSI2-NLK axis in modulating aggressive behaviors of pancreatic cancer cells, which providing new evidence for therapeutic strategies in pancreatic cancer.

Assembly and annotation of the first complete mitochondrial genome of Phenacoccus solenopsis Tinsley (Hemiptera: Pseudococcidae)

The cotton mealybug, Phenacoccus solenopsis Tinsley (Hemiptera: Pseudococcidae), is a significant economic pest with a global distribution. Despite the increasing volume of literature on mitochondrial genome sequencing, complete mitochondrial genome sequences have not been reported for P. solenopsis yet. Here, we assembled the complete mitochondrial genome of P. solenopsis using high throughput DNA sequencing technology. The genome is 14,831 bp in length and is comprised of 37 genes: 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, and two ribosomal RNA (rRNA) genes. The total length of all PCGs is 9,678 bp, accounting for 65.2% of the overall mitogenome. Among the PCGs, atp8 is the smallest gene (99 bp) and nad5 is the largest (1,593 bp). Of the 13 PCGs, seven (nad3, cox3, atp6, atp8, cox2, cox1, nad2) are encoded on the majority strand (J-strand), while six (nad1, nad6, cob, nad4L, nad4, nad5) are on the minority strand (N-strand). The analysis revealed a predominant A+T base composition, making up 90.3% of the total genome. However, the non-coding regulatory control region (CR) is missing due to the non-overlapping endpoints of the linear assembly. This assembly provides comprehensive information for investigating the evolutionary relationships between scale insects and for the precise identification of insects.