A unified analysis of evolutionary and population constraint in protein domains highlights structural features and pathogenic sites.

Protein evolution is constrained by structure and function, creating patterns in residue conservation that are routinely exploited to predict structure and other features. Similar constraints should affect variation across individuals, but it is only with the growth of human population sequencing that this has been tested at scale. Now, human population constraint has established applications in pathogenicity prediction, but it has not yet been explored for structural inference. Here, we map 2.4 million population variants to 5885 protein families and quantify residue-level constraint with a new Missense Enrichment Score (MES). Analysis of 61,214 structures from the PDB spanning 3661 families shows that missense depleted sites are enriched in buried residues or those involved in small-molecule or protein binding. MES is complementary to evolutionary conservation and a combined analysis allows a new classification of residues according to a conservation plane. This approach finds functional residues that are evolutionarily diverse, which can be related to specificity, as well as family-wide conserved sites that are critical for folding or function. We also find a possible contrast between lethal and non-lethal pathogenic sites, and a surprising clinical variant hot spot at a subset of missense enriched positions.

Fuzzy kernel evidence Random Forest for identifying pseudouridine sites.

Pseudouridine is an RNA modification that is widely distributed in both prokaryotes and eukaryotes, and plays a critical role in numerous biological activities. Despite its importance, the precise identification of pseudouridine sites through experimental approaches poses significant challenges, requiring substantial time and resources.Therefore, there is a growing need for computational techniques that can reliably and quickly identify pseudouridine sites from vast amounts of RNA sequencing data. In this study, we propose fuzzy kernel evidence Random Forest (FKeERF) to identify pseudouridine sites. This method is called PseU-FKeERF, which demonstrates high accuracy in identifying pseudouridine sites from RNA sequencing data. The PseU-FKeERF model selected four RNA feature coding schemes with relatively good performance for feature combination, and then input them into the newly proposed FKeERF method for category prediction. FKeERF not only uses fuzzy logic to expand the original feature space, but also combines kernel methods that are easy to interpret in general for category prediction. Both cross-validation tests and independent tests on benchmark datasets have shown that PseU-FKeERF has better predictive performance than several state-of-the-art methods. This new method not only improves the accuracy of pseudouridine site identification, but also provides a certain reference for disease control and related drug development in the future.

Many purported pseudogenes in bacterial genomes are bona fide genes.

BACKGROUND: Microbial genomes are largely comprised of protein coding sequences, yet some genomes contain many pseudogenes caused by frameshifts or internal stop codons. These pseudogenes are believed to result from gene degradation during evolution but could also be technical artifacts of genome sequencing or assembly. RESULTS: Using a combination of observational and experimental data, we show that many putative pseudogenes are attributable to errors that are incorporated into genomes during assembly. Within 126,564 publicly available genomes, we observed that nearly identical genomes often substantially differed in pseudogene counts. Causal inference implicated assembler, sequencing platform, and coverage as likely causative factors. Reassembly of genomes from raw reads confirmed that each variable affects the number of putative pseudogenes in an assembly. Furthermore, simulated sequencing reads corroborated our observations that the quality and quantity of raw data can significantly impact the number of pseudogenes in an assembler dependent fashion. The number of unexpected pseudogenes due to internal stops was highly correlated (R2 = 0.96) with average nucleotide identity to the ground truth genome, implying relative pseudogene counts can be used as a proxy for overall assembly correctness. Applying our method to assemblies in RefSeq resulted in rejection of 3.6% of assemblies due to significantly elevated pseudogene counts. Reassembly from real reads obtained from high coverage genomes showed considerable variability in spurious pseudogenes beyond that observed with simulated reads, reinforcing the finding that high coverage is necessary to mitigate assembly errors. CONCLUSIONS: Collectively, these results demonstrate that many pseudogenes in microbial genome assemblies are actually genes. Our results suggest that high read coverage is required for correct assembly and indicate an inflated number of pseudogenes due to internal stops is indicative of poor overall assembly quality.

SNORA69 is up-regulated in the lateral habenula of individuals with major depressive disorder.

Major depressive disorder (MDD) is a complex and potentially debilitating illness whose etiology and pathology remains unclear. Non-coding RNAs have been implicated in MDD, where they display differential expression in the brain and the periphery. In this study, we quantified small nucleolar RNA (snoRNA) expression by small RNA sequencing in the lateral habenula (LHb) of individuals with MDD (n = 15) and psychiatrically-healthy controls (n = 15). We uncovered five snoRNAs that exhibited differential expression between MDD and controls (FDR < 0.01). Specifically, SNORA69 showed increased expression in MDD and was technically validated via RT-qPCR. We further investigated the expression of Snora69 in the LHb and peripheral blood of an unpredicted chronic mild stress (UCMS) mouse model of depression. Snora69 was specifically up-regulated in mice that underwent the UCMS paradigm. SNORA69 is known to guide pseudouridylation onto 5.8S and 18S rRNAs. We quantified the relative abundance of pseudouridines on 5.8S and 18S rRNA in human post-mortem LHb samples and found increased abundance of pseudouridines in the MDD group. Overall, our findings indicate the importance of brain snoRNAs in the pathology of MDD. Future studies characterizing SNORA69's role in MDD pathology is warranted.

Self-supervised learning on millions of primary RNA sequences from 72 vertebrates improves sequence-based RNA splicing prediction.

Language models pretrained by self-supervised learning (SSL) have been widely utilized to study protein sequences, while few models were developed for genomic sequences and were limited to single species. Due to the lack of genomes from different species, these models cannot effectively leverage evolutionary information. In this study, we have developed SpliceBERT, a language model pretrained on primary ribonucleic acids (RNA) sequences from 72 vertebrates by masked language modeling, and applied it to sequence-based modeling of RNA splicing. Pretraining SpliceBERT on diverse species enables effective identification of evolutionarily conserved elements. Meanwhile, the learned hidden states and attention weights can characterize the biological properties of splice sites. As a result, SpliceBERT was shown effective on several downstream tasks: zero-shot prediction of variant effects on splicing, prediction of branchpoints in humans, and cross-species prediction of splice sites. Our study highlighted the importance of pretraining genomic language models on a diverse range of species and suggested that SSL is a promising approach to enhance our understanding of the regulatory logic underlying genomic sequences.

The novel HLA-DPB1*1500:01N allele, identified by Sanger dideoxy nucleotide sequencing in a Chinese individual.

HLA-DPB1*1500:01N differs from HLA-DPB1*05:01:01:01 by one nucleotide in exon 3.

PDP1 promotes the progression of breast cancer through STAT3 pathway.

This study aimed to investigate the expression pattern and mechanisms of Pyruvate Dehydrogenase Phosphatase Catalytic Subunit 1 (PDP1) in the progression of breast cancer (BC). PDP1, known for its involvement in cell energy metabolism, was found to be overexpressed in BC tissues. Notably, low PDP1 expression aligns with improved overall survival (OS) in BC patients. In this study, we found that PDP1 was overexpressed among BC tissues and low PDP1 expression showed a better prognosis for the patients with BC. PDP1 knockdown suppressed cell amplification and migration and triggered cell apoptosis in BC cells. In vivo assessments through a xenograft model unveiled the pivotal role and underlying mechanisms of PDP1 knockdown. RNA sequencing and kyoto encyclopedia of genes and genomes analysis of RNAs from PDP1 knockdown and normal MCF7 cells revealed 1440 differentially expressed genes, spotlighting the involvement of the JAK/STAT3 signaling pathway in BC progression. Western blot results implied that PDP1 knockdown led to a loss of p-STAT3, whereas overexpression of PDP1 induced the p-STAT3 expression. Cell counting kit-8 assay showed that PDP1 overexpression significantly raised MDA-MB-231 and MCF7 cell viability while STAT3 inhibitor S3I-201 recovered the cell growth to normal level. To summarize, PDP1 promotes the progression of BC through STAT3 pathway by regulating p-STAT3. The findings contribute to understanding the molecular mechanisms underlying BC progression, and opening avenues for targeted therapeutic approaches.

Specific isoforms of the ubiquitin ligase gene WWP2 are targets of osteoarthritis genetic risk via a differentially methylated DNA sequence.

BACKGROUND: Transitioning from a genetic association signal to an effector gene and a targetable molecular mechanism requires the application of functional fine-mapping tools such as reporter assays and genome editing. In this report, we undertook such studies on the osteoarthritis (OA) risk that is marked by single nucleotide polymorphism (SNP) rs34195470 (A > G). The OA risk-conferring G allele of this SNP associates with increased DNA methylation (DNAm) at two CpG dinucleotides within WWP2. This gene encodes a ubiquitin ligase and is the host gene of microRNA-140 (miR-140). WWP2 and miR-140 are both regulators of TGFß signaling. METHODS: Nucleic acids were extracted from adult OA (arthroplasty) and foetal cartilage. Samples were genotyped and DNAm quantified by pyrosequencing at the two CpGs plus 14 flanking CpGs. CpGs were tested for transcriptional regulatory effects using a chondrocyte cell line and reporter gene assay. DNAm was altered using epigenetic editing, with the impact on gene expression determined using RT-qPCR. In silico analysis complemented laboratory experiments. RESULTS: rs34195470 genotype associates with differential methylation at 14 of the 16 CpGs in OA cartilage, forming a methylation quantitative trait locus (mQTL). The mQTL is less pronounced in foetal cartilage (5/16 CpGs). The reporter assay revealed that the CpGs reside within a transcriptional regulator. Epigenetic editing to increase their DNAm resulted in altered expression of the full-length and N-terminal transcript isoforms of WWP2. No changes in expression were observed for the C-terminal isoform of WWP2 or for miR-140. CONCLUSIONS: As far as we are aware, this is the first experimental demonstration of an OA association signal targeting specific transcript isoforms of a gene. The WWP2 isoforms encode proteins with varying substrate specificities for the components of the TGFß signaling pathway. Future analysis should focus on the substrates regulated by the two WWP2 isoforms that are the targets of this genetic risk.

Complete organelle genomes of Korean fir, Abies koreana and phylogenomics of the gymnosperm genus Abies using nuclear and cytoplasmic DNA sequence data.

Abies koreana E.H.Wilson is an endangered evergreen coniferous tree that is native to high altitudes in South Korea and susceptible to the effects of climate change. Hybridization and reticulate evolution have been reported in the genus; therefore, multigene datasets from nuclear and cytoplasmic genomes are needed to better understand its evolutionary history. Using the Illumina NovaSeq 6000 and Oxford Nanopore Technologies (ONT) PromethION platforms, we generated complete mitochondrial (1,174,803 bp) and plastid (121,341 bp) genomes from A. koreana. The mitochondrial genome is highly dynamic, transitioning from cis- to trans-splicing and breaking conserved gene clusters. In the plastome, the ONT reads revealed two structural conformations of A. koreana. The short inverted repeats (1186 bp) of the A. koreana plastome are associated with different structural types. Transcriptomic sequencing revealed 1356 sites of C-to-U RNA editing in the 41 mitochondrial genes. Using A. koreana as a reference, we additionally produced nuclear and organelle genomic sequences from eight Abies species and generated multiple datasets for maximum likelihood and network analyses. Three sections (Balsamea, Momi, and Pseudopicea) were well grouped in the nuclear phylogeny, but the phylogenomic relationships showed conflicting signals in the mitochondrial and plastid genomes, indicating a complicated evolutionary history that may have included introgressive hybridization. The obtained data illustrate that phylogenomic analyses based on sequences from differently inherited organelle genomes have resulted in conflicting trees. Organelle capture, organelle genome recombination, and incomplete lineage sorting in an ancestral heteroplasmic individual can contribute to phylogenomic discordance. We provide strong support for the relationships within Abies and new insights into the phylogenomic complexity of this genus.

Integrated analysis of single-cell and bulk RNA-sequencing reveals a novel signature based on NK cell marker genes to predict prognosis and immunotherapy response in gastric cancer.

Natural killer (NK) cells play essential roles in the tumor development, diagnosis, and prognosis of tumors. In this study, we aimed to establish a reliable signature based on marker genes in NK cells, thus providing a new perspective for assessing immunotherapy and the prognosis of patients with gastric cancer (GC). We analyzed a total of 1560 samples retrieved from the public database. We performed a comprehensive analysis of single-cell RNA-sequencing (scRNA-seq) data of gastric cancer and identified 377 marker genes for NK cells. By performing Cox regression analysis, we established a 12-gene NK cell-associated signature (NKCAS) for the Cancer Genome Atlas (TCGA) cohort, that assigned GC patients into a low-risk group (LRG) or a high-risk group (HRG). In the TCGA cohort, the areas under curve (AUC) value were 0.73, 0.81, and 0.80 at 1, 3, and 5 years. External validation of the predictive ability for the signature was then validated in the Gene Expression Omnibus (GEO) cohorts (GSE84437). The expression levels of signature genes were measured and validated in GC cell lines by real-time PCR. Moreover, NKCAS was identified as an independent prognostic factor by multivariate analysis. We combined this with a variety of clinicopathological characteristics (age, M stage, and tumor grade) to construct a nomogram to predict the survival outcomes of patients. Moreover, the LRG showed higher immune cell infiltration, especially CD8+ T cells and NK cells. The risk score was negatively associated with inflammatory activities. Importantly, analysis of the independent immunotherapy cohort showed that the LRG had a better prognosis and immunotherapy response when compared with the HRG. The identification of NK cell marker genes in this study suggests potential therapeutic targets. Additionally, the developed predictive signatures and nomograms may aid in the clinical management of GC.

The chromosomal-scale genome sequencing and assembly of Athetis lepigone.

Athetis lepigone is an emerging highly polyphagous insect pest reported to cause crop damage in several European and Asian countries. However, our understanding of its genetic adaptation mechanisms has been limited due to lack of high-quality genetic resources. In this study, we present a chromosomal-level genome of A. lepigone, representing the first species in the genus of Athetis. We employed PacBio long-read sequencing and Hi-C technologies to generate 612.49 Mb genome assembly which contains 42.43% repeat sequences with a scaffold N50 of 20.9 Mb. The contigs were successfully clustered into 31 chromosomal-size scaffolds with 37% GC content. BUSCO assessment revealed a genome completeness of 97.4% with 96.3 identified as core Arthropoda single copy orthologs. Among the 17,322 genes that were predicted, 15,965 genes were functionally annotated, representing a coverage of 92.17%. Furthermore, we revealed 106 P450, 37 GST, 27 UGT, and 74 COE gene families in the genome of A. lepigone. This genome provides a significant and invaluable genomic resource for further research across the entire genus of Athetis.

Characterization of the novel HLA-DQB1*02:01:01:21Q allele by sequencing-based typing.

HLA-DQB1*02:01:01:21Q differs from HLA-DQB1*02:01:01:01 by one nucleotide substitution in the splice site in the beginning of intron 3.

The novel HLA-DPB1*1467:01 allele characterised by two different sequencing-based typing techniques.

The novel allele HLA-DPB1*1467:01 differs from HLA-DPB1*09:01:01:01 by one non-synonymous nucleotide substitution in exon 2.

Recognition of the HLA-C*07:359 allele in Taiwanese individuals.

One nucleotide substitution in codon 264 of HLA-C*07:02:01:01 results in a novel allele, HLA-C*07:359.

Induction of the human CDC45 gene promoter activity by natural compound trans­resveratrol.

GGAA motifs in the human TP53 and HELB gene promoters play a part in responding to trans­resveratrol (Rsv) in HeLa S3 cells. This sequence is also present in the 5'­upstream region of the human CDC45 gene, which encodes a component of CMG DNA helicase protein complex. The cells were treated with Rsv (20 µM), then transcripts and the translated protein were analyzed by quantitative RT­PCR and western blotting, respectively. The results showed that the CDC45 gene and protein expression levels were induced after the treatment. To examine whether they were due to the activation of transcription, a 5'­upstream 556­bp of the CDC45 gene was cloned and inserted into a multi­cloning site of the Luciferase (Luc) expression vector. In the present study, various deletion/point mutation­introduced Luc expression plasmids were constructed and they were used for the transient transfection assay. The results showed that the GGAA motif, which is included in a putative RELB protein recognizing sequence, plays a part in the promoter activity with response to Rsv in HeLa S3 cells.

Successful application of modified crude DNA extraction from muscle tissues for various types of PCR amplifications.

BACKGROUND: One of the most challenging aspects of nucleic acid amplification tests is the extraction of genomic DNA. However, achieving satisfactory quality and quantity of genomic DNA is not always easy, while the demand for rapid, low-cost and less laborious DNA isolation methods is ever-increasing. METHODS AND RESULTS: We have developed a rapid (⁓2 min) crude DNA extraction method leading to direct-PCR that requires minimum reagents and laboratory equipment. It was developed by eliminating the time-consuming purification steps of DNA extraction, by processing the sample in optimized amounts of Taq KCl PCR buffer and DNARelease Additive/Proteinase K in only two minutes and carrying out amplification using conventional Taq DNA polymerase. The DNA preparation method was validated on muscle tissue samples from 12 different species as well as 48 cooked meat samples. Its compatibility was also successfully tested with different types of PCR amplification platforms extensively used for genetic analysis, such as simplex PCR, PCR-RFLP (Restriction Fragment Length Polymorphism), multiplex PCR, isothermal amplification, real-time PCR and DNA sequencing. CONCLUSIONS: The developed protocol provides sufficient amount of crude DNA from muscle tissues of different species for PCR amplifications to identify species-of-origin via different techniques coupled with PCR. The simplicity and robustness of this protocol make nucleic acid amplification assays more accessible and affordable to researchers and authorities for both laboratory and point-of-care tests.

Usefulness and Limitations of PFGE Diagnosis and Nucleotide Sequencing Method in the Analysis of Food Poisoning Pathogens Found in Cooking Employees.

In the case of a food poisoning outbreak, it is essential to understand the relationship between cooking workers and food poisoning. Many biological diagnostic methods have recently been developed to detect food poisoning pathogens. Among these diagnostic tools, this study presents PCR-based pulsed-field gel electrophoresis and nucleotide sequencing diagnostic analysis results for diagnosing food poisoning outbreaks associated with cooking employees in Chungcheongnam-do, Republic of Korea. Pulsed-field gel electrophoresis was useful in identifying the food poisoning outbreaks caused by Staphylococcus aureus and Enteropathogenic Escherichia coli. In the case of Norovirus, nucleotide sequencing was used to identify the relationship between cooking workers and the food poisoning outbreak. However, it is difficult to determine whether cooking employees directly caused the food poisoning outbreaks based on these molecular biological diagnostic results alone. A system is needed to integrate epidemiological and diagnostic information to identify a direct correlation between the food poisoning outbreak and cooking employees.

Dynamics of cuticle-associated transcript profiles during moulting of the bed bug Cimexlectularius.

The bed bug Cimex lectularius is a worldwide human pest. The sequenced genome allows molecular analyses of all aspects of bed bug biology. The present work was conducted to contribute to bed bug cuticle biology. As in other insect species, the C. lectularius cuticle consists of the three horizontal layers procuticle, epicuticle and envelope. To analyse the genes needed for the establishment of the stratified cuticle, we studied the expression pattern of 42 key cuticle-related genes at the transition of the penultimate nymphal stage to adult animals when a new cuticle is formed. Based on gene expression dynamics, in simplified model, we distinguish two key events during cuticle renewal in C. lectularius. First, upon blood feeding, modulation of ecdysone signalling culminates in the transcriptional activation of the transcription factor Clec-Ftz-F1 that possibly controls the expression of 32 of the 42 genes tested. Second, timed expression of Clec-Ftz-F1 seems to depend also on the insulin signalling pathway as RNA interference against transcripts of the insulin receptor delays Clec-Ftz-F1 expression and stage transition. An important observation of our transcript survey is that genes needed for the construction of the three cuticle layers are largely expressed simultaneously. Based on these data, we hypothesise a considerable synchronous mechanism of layer formation rather than a strictly sequential one. Together, this work provides a basis for functional analyses of cuticle formation in C. lectularius.

Ultrasimple size encoded microfluidic chip for rapid simultaneous multiplex detection of DNA sequences.

Simultaneous multiplexed analysis can provide comprehensive information for disease diagnosis. However, the current multiplex methods rely on sophisticated barcode technology, which hinders its wider application. In this study, an ultrasimple size encoding method is proposed for multiplex detection using a wedge-shaped microfluidic chip. Driving by negative pressure, microparticles are naturally arranged in distinct stripes based on their sizes within the chip. This size encoding method demonstrates a high level of precision, allowing for accuracy in distinguishing 3-5 sizes of microparticles with a remarkable accuracy rate of up to 99%, even the microparticles with a size difference as small as 0.5 µm. The entire size encoding process is completed in less than 5 min, making it ultrasimple, reliable, and easy to operate. To evaluate the function of this size encoding microfluidic chip, three commonly co-infectious viruses' nucleic acid sequences (including complementary DNA sequences of HIV and HCV, and DNA sequence of HBV) are employed for multiplex detection. Results indicate that all three DNA sequences can be sensitively detected without any cross-interference. This size-encoding microfluidic chip-based multiplex detection method is simple, rapid, and high-resolution, its successful application in serum samples renders it highly promising for potential clinical promotion.