Bioinformatics of germline variant discovery for rare disease diagnostics: current approaches and remaining challenges.

Next-generation sequencing (NGS) has revolutionized the field of rare disease diagnostics. Whole exome and whole genome sequencing are now routinely used for diagnostic purposes; however, the overall diagnosis rate remains lower than expected. In this work, we review current approaches used for calling and interpretation of germline genetic variants in the human genome, and discuss the most important challenges that persist in the bioinformatic analysis of NGS data in medical genetics. We describe and attempt to quantitatively assess the remaining problems, such as the quality of the reference genome sequence, reproducible coverage biases, or variant calling accuracy in complex regions of the genome. We also discuss the prospects of switching to the complete human genome assembly or the human pan-genome and important caveats associated with such a switch. We touch on arguably the hardest problem of NGS data analysis for medical genomics, namely, the annotation of genetic variants and their subsequent interpretation. We highlight the most challenging aspects of annotation and prioritization of both coding and non-coding variants. Finally, we demonstrate the persistent prevalence of pathogenic variants in the coding genome, and outline research directions that may enhance the efficiency of NGS-based disease diagnostics.

HOCOMOCO in 2024: a rebuild of the curated collection of binding models for human and mouse transcription factors.

We present a major update of the HOCOMOCO collection that provides DNA binding specificity patterns of 949 human transcription factors and 720 mouse orthologs. To make this release, we performed motif discovery in peak sets that originated from 14 183 ChIP-Seq experiments and reads from 2554 HT-SELEX experiments yielding more than 400 thousand candidate motifs. The candidate motifs were annotated according to their similarity to known motifs and the hierarchy of DNA-binding domains of the respective transcription factors. Next, the motifs underwent human expert curation to stratify distinct motif subtypes and remove non-informative patterns and common artifacts. Finally, the curated subset of 100 thousand motifs was supplied to the automated benchmarking to select the best-performing motifs for each transcription factor. The resulting HOCOMOCO v12 core collection contains 1443 verified position weight matrices, including distinct subtypes of DNA binding motifs for particular transcription factors. In addition to the core collection, HOCOMOCO v12 provides motif sets optimized for the recognition of binding sites in vivo and in vitro, and for annotation of regulatory sequence variants. HOCOMOCO is available at https://hocomoco12.autosome.org and https://hocomoco.autosome.org.

scATAC-seq preprocessing and imputation evaluation system for visualization, clustering and digital footprinting.

Single-cell ATAC-seq (scATAC-seq) is a recently developed approach that provides means to investigate open chromatin at single cell level, to assess epigenetic regulation and transcription factors binding landscapes. The sparsity of the scATAC-seq data calls for imputation. Similarly, preprocessing (filtering) may be required to reduce computational load due to the large number of open regions. However, optimal strategies for both imputation and preprocessing have not been yet evaluated together. We present SAPIEnS (scATAC-seq Preprocessing and Imputation Evaluation System), a benchmark for scATAC-seq imputation frameworks, a combination of state-of-the-art imputation methods with commonly used preprocessing techniques. We assess different types of scATAC-seq analysis, i.e. clustering, visualization and digital genomic footprinting, and attain optimal preprocessing-imputation strategies. We discuss the benefits of the imputation framework depending on the task and the number of the dataset features (peaks). We conclude that the preprocessing with the Boruta method is beneficial for the majority of tasks, while imputation is helpful mostly for small datasets. We also implement a SAPIEnS database with pre-computed transcription factor footprints based on imputed data with their activity scores in a specific cell type. SAPIEnS is published at: https://github.com/lab-medvedeva/SAPIEnS. SAPIEnS database is available at: https://sapiensdb.com.

p38γ is essential for cell cycle progression and liver tumorigenesis.

The cell cycle is a tightly regulated process that is controlled by the conserved cyclin-dependent kinase (CDK)-cyclin protein complex1. However, control of the G0-to-G1 transition is not completely understood. Here we demonstrate that p38 MAPK gamma (p38γ) acts as a CDK-like kinase and thus cooperates with CDKs, regulating entry into the cell cycle. p38γ shares high sequence homology, inhibition sensitivity and substrate specificity with CDK family members. In mouse hepatocytes, p38γ induces proliferation after partial hepatectomy by promoting the phosphorylation of retinoblastoma tumour suppressor protein at known CDK target residues. Lack of p38γ or treatment with the p38γ inhibitor pirfenidone protects against the chemically induced formation of liver tumours. Furthermore, biopsies of human hepatocellular carcinoma show high expression of p38γ, suggesting that p38γ could be a therapeutic target in the treatment of this disease.

EpiFactors 2022: expansion and enhancement of a curated database of human epigenetic factors and complexes.

We present an update of EpiFactors, a manually curated database providing information about epigenetic regulators, their complexes, targets, and products which is openly accessible at http://epifactors.autosome.org. An updated version of the EpiFactors contains information on 902 proteins, including 101 histones and protamines, and, as a main update, a newly curated collection of 124 lncRNAs involved in epigenetic regulation. The amount of publications concerning the role of lncRNA in epigenetics is rapidly growing. Yet, the resource that compiles, integrates, organizes, and presents curated information on lncRNAs in epigenetics is missing. EpiFactors fills this gap and provides data on epigenetic regulators in an accessible and user-friendly form. For 820 of the genes in EpiFactors, we include expression estimates across multiple cell types assessed by CAGE-Seq in the FANTOM5 project. In addition, the updated EpiFactors contains information on 73 protein complexes involved in epigenetic regulation. Our resource is practical for a wide range of users, including biologists, bioinformaticians and molecular/systems biologists.

The outcome of boosting mitochondrial activity in alcohol-associated liver disease is organ-dependent.

BACKGROUND AND AIMS: Alcohol-associated liver disease (ALD) accounts for 70% of liver-related deaths in Europe, with no effective approved therapies. Although mitochondrial dysfunction is one of the earliest manifestations of alcohol-induced injury, restoring mitochondrial activity remains a problematic strategy due to oxidative stress. Here, we identify methylation-controlled J protein (MCJ) as a mediator for ALD progression and hypothesize that targeting MCJ may help in recovering mitochondrial fitness without collateral oxidative damage. APPROACH AND RESULTS: C57BL/6 mice [wild-type (Wt)] Mcj knockout and Mcj liver-specific silencing (MCJ-LSS) underwent the NIAAA dietary protocol (Lieber-DeCarli diet containing 5% (vol/vol) ethanol for 10 days, plus a single binge ethanol feeding at day 11). To evaluate the impact of a restored mitochondrial activity in ALD, the liver, gut, and pancreas were characterized, focusing on lipid metabolism, glucose homeostasis, intestinal permeability, and microbiota composition. MCJ, a protein acting as an endogenous negative regulator of mitochondrial respiration, is downregulated in the early stages of ALD and increases with the severity of the disease. Whole-body deficiency of MCJ is detrimental during ALD because it exacerbates the systemic effects of alcohol abuse through altered intestinal permeability, increased endotoxemia, and dysregulation of pancreatic function, which overall worsens liver injury. On the other hand, liver-specific Mcj silencing prevents main ALD hallmarks, that is, mitochondrial dysfunction, steatosis, inflammation, and oxidative stress, as it restores the NAD + /NADH ratio and SIRT1 function, hence preventing de novo lipogenesis and improving lipid oxidation. CONCLUSIONS: Improving mitochondrial respiration by liver-specific Mcj silencing might become a novel therapeutic approach for treating ALD.

A Diterpenoid of the Medicinal Plant Andrographis paniculata Targets Cutaneous TRPV3 Channel and Relieves Itch.

Transient receptor potential vanilloid subtype 3 (TRPV3) is an ion channel implicated in skin physiology and itch. TRPV3 inhibitors can present a novel strategy for combating debilitating itch conditions, and medicinal plants are a natural pool of such compounds. Here, we report the isolation of a TRPV3-inhibiting compound from Andrographis paniculata, a medicinal plant with anti-inflammatory properties whose bioactive components are poorly characterized in terms of molecular targets. Using 1H and 13C NMR and high-resolution mass spectrometry, the compound was identified as a labdane-type diterpenoid, 14-deoxy-11,12-didehydroandrographolide (ddA). The activity of the compound was evaluated by fluorescent calcium assay and manual whole-cell patch-clamp technique. ddA inhibited human TRPV3 in stably expressing CHO and HaCaT keratinocytes, acting selectively among other TRP channels implicated in itch and inflammation and not showing toxicity to HaCaT cells. Antipruritic effects of the compound were evaluated in scratching behavior models on ICR mice. ddA suppressed itch induced by the TRPV3 activator carvacrol. Additionally, ddA potently suppressed histamine-induced itch with efficacy comparable to loratadine, a clinically used antihistamine drug. These results suggest the potential of ddA as a possible safe and efficacious alternative for antipruritic therapy.

Comparison of passive targeted delivery of inorganic and organic nanocarriers among different types of tumors.

In this study, we have considered four types of nanoparticles (NPs): polylactic acid (PLA), gold (Au), calcium carbonate (CaCO3), and silica (SiO2) with similar sizes (TEM: 50-110 nm and DLS: 110-140 nm) to examine their passive accumulation in three different tumors: colon (CT26), melanoma (B16-F10), and breast (4T1) cancers. Our results demonstrate that each tumor model showed a different accumulation of NPs, in the following order: CT26 > B16-F10 > 4T1. The Au and PLA NPs were evidently characterized by a higher delivery efficiency in case of CT26 tumors compared to CaCO3 and SiO2 NPs. The Au NPs demonstrated the highest accumulation in B16-F10 cells compared to other NPs. These results were verified using SPECT, ex vivo fluorescence bioimaging, direct radiometry and histological analysis. Thus, this work contributes to new knowledge in passive tumor targeting of NPs and can be used for the development of new strategies for delivery of bioactive compounds.

Genome Analyses Reveal the Secondary Metabolites that Potentially Influence the Geographical Distribution of Fusarium pseudograminearum Populations.

Fusarium crown rot (FCR), caused by Fusarium pseudograminearum, significantly impacts wheat yield and quality in China's Huanghuai region. The rapid F. pseudograminearum epidemic and FCR outbreak within a decade remain unexplained. In this study, two high-quality, chromosome-level genomes of F. pseudograminearum strains producing 3-acetyl-deoxynivalenol (3AcDON) and 15-acetyl-deoxynivalenol (15AcDON) toxins were assembled. Additionally, 38 related strains were resequenced. Genomic differences such as single nucleotide polymorphisms (SNPs), insertions/deletions (indels), and structural variations (SVs) among F. pseudograminearum strains were analyzed. The whole-genome SNP locus-based population classification mirrored the toxin chemotype (3AcDON and 15AcDON)-based classification, indicating the presence of genes associated with the trichothecene toxin gene cluster. Further analysis of differential SNP, indel, and SV loci between the 3AcDON and 15AcDON populations revealed a predominant connection to secondary metabolite synthesis genes. Notably, the majority of the secondary metabolite biosynthesis gene cluster loci were located in SNP-dense genomic regions, suggesting high mutability and a possible contribution to F. pseudograminearum population structure and environmental adaptability. This study provides insightful perspectives on the distribution and evolution of F. pseudograminearum and for forecasting the spread of wheat FCR, thereby aiding in the development of preventive measures and control strategies.

Evolutionary Changes in Primate Glutamate Dehydrogenases 1 and 2 Influence the Protein Regulation by Ligands, Targeting and Posttranslational Modifications.

There are two paralogs of glutamate dehydrogenase (GDH) in humans encoded by the GLUD1 and GLUD2 genes as a result of a recent retroposition during the evolution of primates. The two human GDHs possess significantly different regulation by allosteric ligands, which is not fully characterized at the structural level. Recent advances in identification of the GDH ligand binding sites provide a deeper perspective on the significance of the accumulated substitutions within the two GDH paralogs. In this review, we describe the evolution of GLUD1 and GLUD2 after the duplication event in primates using the accumulated sequencing and structural data. A new gibbon GLUD2 sequence questions the indispensability of ancestral R496S and G509A mutations for GLUD2 irresponsiveness to GTP, providing an alternative with potentially similar regulatory features. The data of both GLUD1 and GLUD2 evolution not only confirm substitutions enhancing GLUD2 mitochondrial targeting, but also reveal a conserved mutation in ape GLUD1 mitochondrial targeting sequence that likely reduces its transport to mitochondria. Moreover, the information of GDH interactors, posttranslational modification and subcellular localization are provided for better understanding of the GDH mutations. Medically significant point mutations causing deregulation of GDH are considered from the structural and regulatory point of view.

Chemical Constituents of Stinging Nettle (Urtica dioica L.): A Comprehensive Review on Phenolic and Polyphenolic Compounds and Their Bioactivity.

Polyphenolic compounds are of great interest in today's science. Naturally, they occur in plants and other sources in many different forms. Their wide range of biological activity has attracted the attention of the scientific community. One of the sources of phenolic compounds is stinging nettle (Urtica dioica L.), a common plant in almost all parts of the world. A long tradition of utilization and an interesting chemical profile make this plant a fascinating and extensive object of study. The chemical profile also allows this plant to be used as a food and a pigment source in the food, pharmaceutical, and cosmetic industries. Previously conducted studies found phenolic acids and polyphenolic compounds in root, stalk, and stinging nettle leaves. Different extraction techniques were usually used to isolate them from the leaves. Obtained extracts were used to investigate biological activity further or formulate different functional food products. This study aimed to collect all available knowledge about this plant, its chemical composition, and biological activity and to summarize this knowledge with particular attention to polyphenolic compounds and the activity and mechanisms of their actions.

Unique Self-Phosphorylating Polybenzimidazole of the 6F Family for HT-PEM Fuel Cell Application.

High-temperature polymer-electrolyte membrane fuel cells (HT-PEMFCs) are a very important type of fuel cells since they operate at 150-200 °C, making it possible to use hydrogen contaminated with CO. However, the need to improve the stability and other properties of gas-diffusion electrodes still impedes their distribution. Self-supporting anodes based on carbon nanofibers (CNF) are prepared using the electrospinning method from a polyacrylonitrile solution containing zirconium salt, followed by pyrolysis. After the deposition of Pt nanoparticles on the CNF surface, the composite anodes are obtained. A new self-phosphorylating polybenzimidazole of the 6F family is applied to the Pt/CNF surface to improve the triple-phase boundary, gas transport, and proton conductivity of the anode. This polymer coating ensures a continuous interface between the anode and proton-conducting membrane. The polymer is investigated using CO2 adsorption, TGA, DTA, FTIR, GPC, and gas permeability measurements. The anodes are studied using SEM, HAADF STEM, and CV. The operation of the membrane-electrode assembly in the H2/air HT-PEMFC shows that the application of the new PBI of the 6F family with good gas permeability as a coating for the CNF anodes results in an enhancement of HT-PEMFC performance, reaching 500 mW/cm2 at 1.3 A/cm2 (at 180 °C), compared with the previously studied PBI-O-PhT-P polymer.

Prediction of the Aggressive Clinical Course of Papillary Thyroid Carcinoma Based on Fine Needle Aspiration Biopsy Molecular Testing.

Molecular genetic events are among the numerous factors affecting the clinical course of papillary thyroid carcinoma (PTC). Recent studies have demonstrated that aberrant expression of miRNA, as well as different thyroid-related genes, correlate with the aggressive clinical course of PTC and unfavorable treatment outcomes, which opens up new avenues for using them in the personalization of the treatment strategy for patients with PTC. In the present work, our goal was to assess the applicability of molecular markers in the preoperative diagnosis of aggressive variants of papillary thyroid cancer. The molecular genetic profile (expression levels of 34 different markers and BRAF mutations) was studied for 108 cytology specimens collected by fine-needle aspiration biopsy in patients with PTC having different clinical manifestations. Statistically significant differences with adjustment for multiple comparisons (p < 0.0015) for clinically aggressive variants of PTC were obtained for four markers: miRNA-146b, miRNA-221, fibronectin 1 (FN1), and cyclin-dependent kinase inhibitor 2A (CDKN2A) genes. A weak statistical correlation (0.0015 < p < 0.05) was observed for miRNA-31, -375, -551b, -148b, -125b, mtDNA, CITED1, TPO, HMGA2, CLU, NIS, SERPINA1, TFF3, and TMPRSS4. The recurrence risk of papillary thyroid carcinoma can be preoperatively predicted using miRNA-221, FN1, and CDKN2A genes.

Evaluation of Ouabain's Tissue Distribution in C57/Black Mice Following Intraperitoneal Injection, Using Chromatography and Mass Spectrometry.

Cardiotonic steroids (CTSs), such as digoxin, are used for heart failure treatment. However, digoxin permeates the brain-blood barrier (BBB), affecting central nervous system (CNS) functions. Finding a CTS that does not pass through the BBB would increase CTSs' applicability in the clinic and decrease the risk of side effects on the CNS. This study aimed to investigate the tissue distribution of the CTS ouabain following intraperitoneal injection and whether ouabain passes through the BBB. After intraperitoneal injection (1.25 mg/kg), ouabain concentrations were measured at 5 min, 15 min, 30 min, 1 h, 3 h, 6 h, and 24 h using HPLC-MS in brain, heart, liver, and kidney tissues and blood plasma in C57/black mice. Ouabain was undetectable in the brain tissue. Plasma: Cmax = 882.88 ± 21.82 ng/g; Tmax = 0.08 ± 0.01 h; T1/2 = 0.15 ± 0.02 h; MRT = 0.26 ± 0.01. Cardiac tissue: Cmax = 145.24 ± 44.03 ng/g (undetectable at 60 min); Tmax = 0.08 ± 0.02 h; T1/2 = 0.23 ± 0.09 h; MRT = 0.38 ± 0.14 h. Kidney tissue: Cmax = 1072.3 ± 260.8 ng/g; Tmax = 0.35 ± 0.19 h; T1/2 = 1.32 ± 0.76 h; MRT = 1.41 ± 0.71 h. Liver tissue: Cmax = 2558.0 ± 382.4 ng/g; Tmax = 0.35 ± 0.13 h; T1/2 = 1.24 ± 0.7 h; MRT = 0.98 ± 0.33 h. Unlike digoxin, ouabain does not cross the BBB and is eliminated quicker from all the analyzed tissues, giving it a potential advantage over digoxin in systemic administration. However, the inability of ouabain to pass though the BBB necessitates intracerebral administration when used to investigate its effects on the CNS.

A Combination of Library Screening and Rational Mutagenesis Expands the Available Color Palette of the Smallest Fluorogen-Activating Protein Tag nanoFAST.

NanoFAST is the smallest fluorogen-activating protein, consisting of only 98 amino acids, used as a genetically encoded fluorescent tag. Previously, only a single fluorogen with an orange color was revealed for this protein. In the present paper, using rational mutagenesis and in vitro screening of fluorogens libraries, we expanded the color palette of this tag. We discovered that E46Q is one of the key substitutions enabling the range of possible fluorogens to be expanded. The introduction of this and several other substitutions has made it possible to use not only orange but also red and green fluorogens with the modified protein.

Knockout of the Tnfa Gene Decreases Influenza Virus-Induced Histological Reactions in Laboratory Mice.

Tumor necrosis factor alpha (TNF-α) is a cytokine that is responsible for many processes associated with immune response and inflammation. It is involved in the development of an antiviral response to many virus infections. This factor was shown to be activated in influenza A virus infection, which enhances production of other cytokines. The overexpression of these cytokines can lead to a cytokine storm. To study the role of TNF-α in the development of pathologies associated with viral infection, we generated a Tnfa knockout mouse strain. We demonstrated that these mice were characterized by a significant increase in the number of viral genomes compared to that in the parental strain, but the amount of live virus did not differ. A histopathology of the lungs in the genetically modified animals was significantly lower in terms of interalveolar septal infiltration. The generated model may be used to further study pathological processes in viral infections.

Expectations and perspectives of cognitive behavioural therapy for childhood anxiety and related disorders.

BACKGROUND: Anxiety disorders are the most frequently diagnosed psychiatric conditions in children and adolescents. Cognitive behavioural therapy (CBT) is a well-established and effective treatment for anxiety and related disorders across the lifespan. Expectations of psychotherapy have been demonstrated to affect outcomes, yet there is sparse existing literature on adolescent patient and parent perspectives of CBT prior to engagement with treatment. AIMS: This study aimed to qualitatively explore the expectations and perceptions of CBT for anxiety and related disorders among adolescent patients and parents. METHOD: Fourteen adolescent patients and 16 parents participated in semi-structured individual interviews or focus groups consisting of 2-3 participants. Interview transcripts were analysed using inductive analysis. RESULTS: Three themes were identified: worries about CBT, expectations and knowledge of the CBT process, and the role of parents and families. Overall, we found that adolescents and parents had generally positive views of CBT. The outset of CBT saw adolescents and parents express concern about stigma as well as the ambiguity of CBT. Parents continued to express a lack of understanding of what CBT entailed during their child's treatment course. CONCLUSION: These results suggest that both adolescents and parents would benefit from early discussion and reinforcement of expectations for CBT treatment. Further research efforts are warranted and should be directed towards determining appropriate expectations for parental involvement in a child's CBT course and effective communication of treatment expectations to both adolescents and parents.

The effects of chronic high-dose morphine on microgliosis and the microglial transcriptome in rat spinal cord.

Background: Opioids are efficacious and safe analgesic drugs in short-term use for acute pain but chronic use can lead to tolerance and dependence. Opioid-induced microglial activation may contribute to the development of tolerance and this process may differ between males and females. A link is suggested between this microglial activation and inflammation, disturbances of circadian rhythms, and neurotoxic effects. We set out to further delineate the effects of chronic morphine on pain behaviour, microglial and neuronal staining, and the transcriptome of spinal microglia, to better understand the role of microglia in the consequences of long-term high-dose opioid administration. Experimental Approach: In two experiments, we administered increasing subcutaneous doses of morphine hydrochloride or saline to male and female rats. Thermal nociception was assessed with the tail flick and hot plate tests. In Experiment I, spinal cord (SC) samples were prepared for immunohistochemical staining for microglial and neuronal markers. In Experiment II, the transcriptome of microglia from the lumbar SC was analysed. Key Results: Female and male rats had similar antinociceptive responses to morphine and developed similar antinociceptive tolerance to thermal stimuli following chronic increasing high doses of s.c. morphine. The area of microglial IBA1-staining in SC decreased after 2 weeks of morphine administration in both sexes. Following morphine treatment, the differentially expressed genes identified in the microglial transcriptome included ones related to the circadian rhythm, apoptosis, and immune system processes. Conclusions: Female and male rats showed similar pain behaviour following chronic high doses of morphine. This was associated with decreased staining of spinal microglia, suggesting either decreased activation or apoptosis. High-dose morphine administration also associated with several changes in gene expression in SC microglia, e.g., those related to the circadian rhythm (Per2, Per3, Dbp). These changes should be considered in the clinical consequences of long-term high-dose administration of opioids.

Drought Tolerance of Legumes: Physiology and the Role of the Microbiome.

Water scarcity and global warming make drought-tolerant plant species more in-demand than ever. The most drastic damage exerted by drought occurs during the critical growth stages of seed development and reproduction. In the course of their evolution, plants form a variety of drought-tolerance mechanisms, including recruiting beneficial microorganisms. Legumes (one of the three largest groups of higher plants) have unique features and the potential to adapt to abiotic stress. The available literature discusses the genetic (breeding) and physiological aspects of drought tolerance in legumes, neglecting the role of the microbiome. Our review aims to fill this gap: starting with the physiological mechanisms of legume drought adaptation, we describe the symbiotic relationship of the plant host with the microbial community and its role in facing drought. We consider two types of studies related to microbiomes in low-water conditions: comparisons and microbiome engineering (modulation). The first type of research includes diversity shifts and the isolation of microorganisms from the various plant niches to which they belong. The second type focuses on manipulating the plant holobiont through microbiome engineering-a promising biotech strategy to improve the yield and stress-resistance of legumes.

PCID2 Subunit of the Drosophila TREX-2 Complex Has Two RNA-Binding Regions.

Drosophila PCID2 is a subunit of the TREX-2 mRNA nuclear export complex. Although the complex has long been studied in eukaryotes, it is still unclear how TREX-2 interacts with mRNA in multicellular organisms. Here, the interaction between Drosophila PCID2 and the ras2 RNA was studied by EMSA. We show that the C-terminal region of the WH domain of PCID2 specifically binds the 3'-noncoding region of the ras2 RNA. While the same region of PCID2 interacts with the Xmas-2 subunit of the TREX-2 complex, PCID2 interacts with RNA independently of Xmas-2. An additional RNA-binding region (M region) was identified in the N-terminal part of the PCI domain and found to bind RNA nonspecifically. Point mutations of evolutionarily conserved amino acid residues in this region completely abolish the PCID2-RNA interaction, while a deletion of the C-terminal domain only partly decreases it. Thus, the specific interaction of PCID2 with RNA requires nonspecific PCID2-RNA binding.