Meloidogyne javanica causing root-knot nematode in soybean is an important problem in soybean areas, leading to several yield losses. Some accessions have been identified carrying resistance loci to this nematode specie. In this study, a set of 317 soybean accessions were characterized for resistance to M. javanica. Genome-wide association study (GWAS) was performed using SNPs from genotyping-by-sequencing (GBS), and a region of 29.2 Kbp on chromosome 13 was identified. The haplotype analysis showed that SNPs were able to discriminate susceptible and resistant accessions, leading to 25 accessions sharing the resistance locus. Furthermore, 5 accessions may be new M. javanica resistance sources. The screening of the SNPs in the USDA soybean germplasm showed that several accessions previous reported as resistance to other nematodes also showed the resistance haplotype on chromosome 13. High levels of concordance among the phenotypes of Brazilian cultivars and the SNPs in chromosome 13 were observed. A in silico analysis of the mapped region on soybean genome revealed a presence of 5 genes with structural similarity with major resistance genes. The expression levels of the candidate genes in the interval demonstrated a potential pseudogene, and other two model genes up-regulated in the resistance source after pathogen infection. The SNPs associated to the region conferring resistance is a important tool for introgression of the resistance by marker-assisted selection in soybean breeding programs.

Background: Studies employing next-generation sequencing (NGS) show that the oral fungal community (mycobiome) is far more complex than hitherto thought. However, the role of the oral mycobiome in health and disease, including oral carcinogenesis, has not been explored. Objective: To characterize the mycobiome associated with oral squamous cell carcinoma (OSCC). Methods: Tissue biopsies [cases: 25 OSCC; controls: 27 intra-oral fibro-epithelial polyp (FEP)] were collected from oral and maxillofacial units in Sri Lanka. Total DNA was extracted and subjected to sequencing of the fungal ITS2 region using Illumina’s 2x300 bp chemistry. High quality, non-chimeric merged reads were classified to species level using a BLASTN-algorithm with UNITE’s named species sequences as reference. Downstream analyses were performed using QIIME and LEfSe. Results: 364 species representing 160 genera and 2 phyla (Ascomycota and Basidiomycota) were identified, with Candida and Malassezia making up 48% and 11% of the average mycobiome, respectively. However, only 5 species and 4 genera were detected in ≥50% of the samples. The species richness and diversity were significantly lower in OSCC. At the genus level, Candida, Hannaella and Gibberella were overrepresented in OSCC while Alternaria and Trametes were more abundant in FEP. Species-wise, C. albicans, C. etchellsii and Hannaella luteola-like species were enriched in OSCC while Malassezia restricta, Aspergillus tamarii, Alternaria alternate, Cladosporium halotolerans, and Hanseniaspora uvarum-like species were the most significantly abundant in FEP. Conclusions: A dysbiotic mycobiome dominated by C. albicans was found in association with OSCC. Whether this dysbiosis plays a role in oral carcinogenesis warrants further investigation.

Heteromorphic sex chromosomes are rarely found in plants. They were observed only in 47 species from phylogenetically distant families, suggesting that the evolution of sex chromosomes was independent in these species. It was shown that DNA repeat sequences are one of the major factors driving sex chromosomes evolution, and an accumulation or elimination of the repetitive DNA elements are closely linked with the formation of differences in the sex chromosomes. The goal of this study was to characterize the transposon composition in male and female plants of Cannabis sativa L., Humulus lupulus L. and Humulus japonicus Siebold & Zucc. For the first time, the male and female genomes of H. japonicus as well as male genomes of H. lupulus and C. sativa have been sequenced (there were no open data about them). The analysis of genome-wide sequencing data with using Repeatexplorer2 and author’s scripts was carried out. It was shown that accumulation of Ty3-gypsy may be associated with speciation in Cannabaceae family which is the opposite of the theory of speciation throw whole-genome duplication. Moreover, the sex-specific DNA repeat clusters in C. sativa and H. japonicus were found. The analysis also revealed that the concentration of Tekay, Retand and Ikeros repeats in the Y chromosome of C. sativa is lower than in the X chromosome and the Angela concentration is higher in the Y chromosome.

: The primary objective of this study was to investigate the association of certain genetic alterations and intraoperative fluorescent activity of 5-aminolevulinic acid (ALA) in brain metastasis (BM) of lung adenocarcinoma. A retrospective cohort study was conducted among 72 patients who underwent surgical resection of BM from lung adenocarcinoma at our institute for 5 years. Cancer cell infiltration was estimated by the intraoperative fluorescent activity of 5-ALA, and genetic alterations were analyzed by next-generation sequencing (NGS). The sensitivity and specificity for detecting cancer cell infiltration using 5-ALA were 87.5% and 96.4%, respectively. Genes associated with cell cycle regulation (p = 0.003) and cell proliferation (p = 0.044) were significantly associated with positive fluorescence activity of 5-ALA in the adjacent brain tissue. Genetic alterations in cell cycle regulation and cell proliferation were also associated with shorter recurrence-free survival (p = 0.013 and p = 0.042, respectively) and overall survival (p = 0.026 and p = 0.042, respectively) in the multivariate analysis. The results suggest that genetic alterations in cell cycle regulation and cell proliferation are associated with positive fluorescence activity of 5-ALA in the adjacent infiltrative brain tissue and influence the clinical outcome of BM of lung adenocarcinoma.

Cutaneous squamous cell carcinoma is a prevalent malignancy with a rising incidence and notably high mutational load. Exploring the genetic nuances of cSCC and investigating molecular approaches stands as a potential avenue for improving outcome in high-risk patients. This retrospective case-control study involved two cohorts of 14 patients (“discovery cohort”) and 12 patients (“validation cohort”) with cSCC located in head-neck anatomical district and diagnosed at pT2 stage. Overall, cases developed early local relapse of the disease whereas controls never re-lapsed during the entire follow-up period. Next generation sequencing (NGS) approach conducted on histological samples revealed that TP53 and CDKN2A were the most frequently mutated genes in our series. No specific mutation was identified as potential prognostic or therapeutic targets. Controls exhibited a tendency toward a higher mutational rate compared to cases. It is possible that an increased number of mutations could prompt the cSCC to expose more antigens, becoming more immunogenic and facilitating the recognition by the immune system. This could enhance and sustain immunological response, potentially preventing future recurrences.

The molecular basis of Down syndrome (DS) predisposition to leukemia is not fully understood but involves various factors such as chromosomal abnormalities, oncogenic mutations, epigenetic alterations, and changes in selection dynamics. Myeloid leukemia associated with DS (ML-DS) is preceded by a preleukemic phase called transient abnormal myelopoiesis driven by GATA1 gene mutations and progresses to ML-DS through additional mutations in cohesin genes, CTCF, RAS, or JAK/STAT pathway genes. DS-related ALL (ALL-DS) differs from non-DS ALL in terms of cytogenetic subgroups and genetic driver events and aberrant expression of CRLF2, JAK2 mutations, and RAS pathway activating mutations are frequent in ALL-DS. Recent advancements in single-cell multi-omics technologies have provided unprecedented insights into the cellular and molecular heterogeneity of DS-associated hematologic neoplasms. Single-cell RNA sequencing and digital spatial profiling enable the identification of rare cell subpopulations, characterization of clonal evolution dynamics, and exploration of the tumor microenvironment's role. These approaches may help identify new druggable targets and tailor therapeutic interventions based on distinct molecular profiles, ultimately improving patient outcomes with the potential to guide personalized medicine approaches and the development of targeted therapies.

Numerous dynamic and complicated processes characterize development from the oocyte to the embryo. However, given the importance of functional transcriptome profiles, long non-coding RNAs, single-nucleotide polymorphisms, and alternative splicing during embryonic development, the effect that these features have on the blastomeres of 2-, 4-, 8-, 16-cell, and morula stages of development have not been studied. Here, we conducted a scRNA-seq survey of cells from sheep from the oocyte to the blastocyst developmental stages. We then carried out experiments to identify and functionally analyze the transcriptome profiles, long non-coding RNAs, single-nucleotide polymorphisms (SNPs), and alternative splicing (AS). We founded that between the oocyte and zygote groups significantly down-regulated genes and the second-largest change in gene expression occurred between the 8- and 16-cell stages. We used various methods to construct a profile to characterize cellular and molecular features and systematically analyze the related GO and KEGG profile of cells of all stages from the oocyte to the blastocyst. This large-scale, single-cell atlas provides key cellular information and will likely assist clinical studies in improving preimplantation genetic diagnosis.

Whole-genome sequencing (WGS) was used for genomic characterization of 110 strains of Listeria innocua isolated from 23 cattle farms, 8 beef abattoirs, and 48 retail outlets in Gauteng province, South Africa. In silico multilocus sequence typing (MLST) was used to identify the isolates' sequence types (STs). BLAST-based analyses were used to identify antimicrobial and virulence genes. The study also linked the detection of the genes to the origin (sources and types of samples) of the L. innocua isolates. The study detected 12 STs, 13 resistance genes, and 23 virulence genes. Of the 12 STs de-tected, ST637 (26.4%), ST448 (20%), 537 (13.6%), and 1085 (12.7%) were predominant, and the frequency varied significantly (P<0.05). All 110 isolates of L. innocua were car-riers of one or more antimicrobial resistance genes, with resistance genes lin (100%), fosX (100%), and tet(M) (30%) being the most frequently detected (P<0.05). Of the 23 virulence genes recognized, 13 (clpC, clpE, clpP, hbp1.svpA, hbp2, iap/cwhA, lap, lpeA, lplA1, lspA, oatA, pdgA, and prsA2) were found in all 110 isolates of L. innocua. Overall, diversity and significant differences were detected in the frequencies of STs, resistance, and virulence genes according to the origins (source and sample type) of the L. innocua isolates. This, being the first genomic characterization of L. innocua recovered from the three levels (farm, abattoir, and retail outlets) of the beef production system in South Africa, provides data on the organism's distribution and potential food safety implications..

Rapid advances in the screening, diagnosis, and treatment of genetic disorders has increased the number of conditions that could be detected through universal newborn screening (NBS). However, addition of conditions to the Recommended Uniform Screening Panel (RUSP) and implementation of nationwide screening has been a slow process taking several years to accomplish for individual conditions. Here we describe web-based tools and resources developed and implemented by the newborn screening translational research network (NBSTRN) to advance newborn screening research and support NBS stakeholders. NBSTRN’s tools include the Longitudinal Pediatric Data Resource (LPDR), the NBS Condition Resource (NBS-CR), the NBS Virtual Repository (NBS-VR), and the Ethical, Legal, and Social Issues (ELSI) Advantage. Research programs including the Inborn Errors of Metabolism Information System (IBEM-IS), BabySeq, EarlyCheck, and Family Narratives Use Cases have utilized NBSTRN’s tools, and in turn contributed research data to further expand and refine these this resource. Additionally, we discuss ongoing tool development to facilitate expansion of genetic disease screening in increasingly diverse populations. In conclusion, NBSTRN’s tools and resources provide a trusted platform to enable NBS stakeholders to advance NBS research and improve clinical care for patients and their families.

Large population passages of the SARS-CoV-2 in the past two and a half years have allowed the circulating virus to accumulate an increasing number of mutations in its genome. The most recently emerging Omicron subvariants have the highest number of mutations in the Spike (S) protein gene and these mutations mainly occur in the receptor-binding domain (RBD) and the N-terminal domain (NTD) of the S gene. The eCDC and the WHO recommend partial Sanger sequencing of the SARS-CoV-2 S gene RBD and NTD on the PCR-positive samples in diagnostic laboratories as a practical means of determining the variants of concern to monitor a possible increased transmissibility, increased virulence, or reduced effectiveness of vaccines against them. The author’s diagnostic laboratory has implemented the eCDC/WHO recommendation by sequencing a 398-base segment of the N gene for the definitive detection of SARS-CoV-2 in clinical samples, and sequencing a 445-base segment of the RBD and a 490-509-base segment of the NTD for variant determination. This paper presents 5 selective cases to illustrate the challenges of using Sanger sequencing to diagnose Omicron subvariant when the samples harbor a high level of co-existing minor subvariant sequences with multi-allelic single nucleotide polymorphisms (SNPs) or possible recombinant Omicron subvariants containing a BA.1 NTD and a BA.2 RBD, which can only be detected by using specially designed PCR primers. The current large-scale surveillance programs using next-generation sequencing (NGS) do not face similar problems because NGS focuses on deriving consensus sequence.

A 23-month-old neutered male dog of unknown ancestry presented with a history of progressive neurological signs that included anxiety, cognitive impairment, tremors, seizure activity, ataxia, and pronounced visual impairment. The clinical signs were accompanied by global brain atrophy. Due to progression in the severity of disease signs, the dog was euthanized at 26 months of age. Examination of tissues collected at necropsy revealed dramatic intracellular accumulations of autofluorescent inclusions in the brain, retina, and cardiac muscle. The inclusions were immunopositive for subunit c of mitochondrial ATP synthase, and their ultrastructural appearances were similar to those of lysosomal storage bodies that accumulate in some of the neuronal ceroid lipofuscinosis (NCL) diseases. The also dog exhibited widespread neuroinflammation. Based on these findings, the dog was deemed likely to have suffered from a form of NCL. A whole genome sequence analysis of the proband’s DNA revealed a homozygous C to T substitution that altered the intron 3 – exon 4 splice site of CLN6. Other mutations in CLN6 cause NCL diseases in humans and animals including dogs. CLN6 protein was undetectable with immunolabeling in the tissues of the proband. Based on the clinical history, fluorescence and electron-microscopy, immunohistochemistry, and molecular genetic findings, the disorder in this dog is classified as an NCL resulting from the absence of CLN6 protein. Screening the dog’s genome for a panel of breed-specific polymorphisms indicated that its ancestry included numerous breeds, with no single breed predominating. This suggests that the CLN6 disease variant is likely to be present in other mixed breed dogs and at least some of the ancestral breeds.

Colletotrichum lini is a flax fungal pathogen, comprising differently virulent strains. In this work, we sequenced the genomes of three C. lini strains of different virulence: high virulent #390-1, medium virulent #757, and low virulent #771. We obtained more than 100× genome coverage with Oxford Nanopore Technologies reads and more than 50× genome coverage with Illumina data. To construct high-quality genomes, several assembly strategies were tested. The final assemblies were obtained using the Canu – Racon ×2 – Medaka – Polca scheme and were 54.0-55.3 Mb in length, contained 26-32 contigs, had an N50 > 5 Mb and the BUSCO completeness > 96%. Comparative analysis of the constructed genomes showed high similarity between mitochondrial and nuclear genomes. However, a rearrangement event and the loss of a 0.9 Mb contig were revealed. The candidate effectors were predicted in sets of proteins received after genome annotation. The obtained genome sequences of C. lini strains extend the knowledge on the genetic diversity of Colletotrichum species and are the basis for establishing molecular mechanisms of pathogenicity of the flax pathogen and effective control of flax anthracnose.

Background: After its initial detection in Wuhan, China, in December 2019, SARS-CoV-2 has spread rapidly, causing successive epidemic waves worldwide. This study aims to provide a genomic epidemiology of SARS-CoV-2 in Burkina Faso. Methods: Three hundred and seventy-seven SARS-CoV-2 genomes obtained from PCR-positive nasopharyngeal samples (PCR cycle threshold score <35) collected between May 5, 2020, and January 31, 2022 were analysed. Genomic sequences were assigned to phylogenetic clades using NextClade and to Pango lineages using pangolin. Phylogenetic and phylogeographic analyses were performed to determine the geographical sources and time of virus introduction in Burkina Faso. Results: The analyzed SARS-CoV-2 genomes could be assigned to 10 phylogenetic clades and 27 Pango lineages already described worldwide. Our analyses revealed the important role of cross-border human mobility in the successive SARS-CoV-2 introductions in Burkina Faso from neighboring countries. Conclusion: This study provides additional insights into the genomic epidemiology of SARS-CoV-2 in West Africa. It highlights the importance of land travel in the spread of the virus and the need to rapidly implement preventive policies. Regional cross-border collaborations and the adherence of the general population to government policies are key to prevent new epidemic waves.

In this report, we describe the course and successful treatment of complicated infective endocarditis (IE). A patient presented with a high-grade irregular fever with chills lasting at least 2 months, dyspnoea, chest pain, fatigue, weight loss, and night sweats during the previous 3 months. Above cardiac congenital disorders, he was found to have Granulicatella adiacens infective aortic valve endocarditis, presumably transmitted from the oral cavity niche. Validated metagenomic 16S rDNA next generation sequencing was used to perform taxonomic identification, allowing for specific adequate antibiotic therapy instead of empiric therapy. This paper highlights the critical role of rapid taxonomic identification of nutritionally variant streptococci crucial and the benefit of proper IE treatment to avoid relapses or fatal complications.

Atrial fibrillation (AF) is a form of sustained cardiac arrhythmia and microRNAs (miRs) play crucial roles in pathophysiology of AF. To identify novel miR-mRNA pairs, we performed RNAseq from atrial biopsies of AF and non-AF patients. Differentially expressed miRs (11-down and 9-up) and mRNAs (95-up and 82-down) were identified and hierarchically clustered in a heat-map. Subsequently, GO, KEGG, and GSEA analyses were run to identify deregulated pathways. Then, miR targets were predicted in miRDB database, and a regulatory network of negatively correlated miR-mRNA pairs was constructed using Cytoscape. To select potential candidate genes from GSEA analysis, top-50 enriched genes in GSEA were overlaid with predicted targets of differentially deregulated miRs. Besides, protein-protein-interaction (PPI) network of enriched genes in GSEA was constructed, and subsequently GO and canonical pathway analyses were run for genes in PPI network. Our analyses showed that TNF-α, p53, EMT, and SYDECAN1 signaling were among the highly affected pathways in AF samples. SDC-1 (syndecan-1) was the top-enriched gene in p53, EMT, and SYDECAN1 signaling. Consistently, SDC-1 mRNA and protein levels were significantly higher in atrial samples of AF patients. Among negatively correlated miRs, miR-302b-3p was experimentally validated to suppress SDC-1 transcript levels. Overall, our results suggested that miR-302b-3p/SDC-1 axis may involve in pathogenesis of AF.

Cancer immunotherapy reaches a remarkable achievement in various cancer types and brings new possibilities to improve cancer patients’ long-term survival. However, outcomes vary from case to case, and the present protocol benefits a small fraction of patients. One notable factor is the tumor microenvironment, especially the immune cell components, that may reflect the immune response's status quo on site. Thus, understanding the content of infiltrating immune cells in tumors is not only for research interesting but also a crucial subject toward precision medicine. We implement an algorithm for resolving relative proportions of twenty-one immune cell subclasses from a human tissue profiled transcriptome by microarray technology to reach the goal above. By selecting gene features and then adopting ?-Support Vector Regression, we can construct a deconvolution model and resolve the immune cell context. The excellent consistency between the estimated values and the correct immune-cell composition further demonstrates this approach provides a more natural alternative to revealing samples' immune cell content and reliable results like recent single-cell technologies. Based on this algorithm, the web-based deconvolution tool implemented named mySORT provides a user-friendly interface for estimating the immune cell content by uploading gene expression profiling. We also present comprehensive visualization 2D/3D plots in mySORT so that users can easily make a comparison between different samples. Finally, we synthesized pseudo-bulk expression data from single-cell transcriptomic datasets of 17 melanoma and 16 head and neck cancer patients. The deconvolution results of microarray-based data in the previous study and synthetic pseudo-bulk data all proved the excellent performance of mySORT. We believe that mySORT can help researchers in all fields easily understand complex immune microenvironment. The website of mySORT is freely accessible on https://symbiosis.iis.sinica.edu.tw/mySORT/.

Congenital Pouch Colon (CPC) is a rare anorectal anomaly common to North Western India specifically Rajasthan. Despite efforts to understand the clinical genetic makeup of CPC, no attempt on identifying non-coding RNAs was done. We have earlier reported CPC's rare variants from whole exome sequencing across 18 affected samples in a total of 64 subjects. A Smith-Waterman algorithm was used to infer a couple of lncRNAs from WES samples of CPC with predictions from the Noncode database. Further screening and quantification using PCR, we ascertained interactions using Micro Scale Thermophoresis (MST). We report the role of lnc-EPB41-1-1 shown to be promiscuously interacting with KIF13A substantiating their role in regulation.

Hypervirulent Klebsiella pneumoniae (hvKp) is a variant that has been increasingly linked to severe, life threatening infections including pyogenic liver abscess and bloodstream infections. HvKp belonging to the capsular serotypes K1 and K2 have been reported worldwide, however very scarce studies are available on the their genomics and virulence. In the current study we report four hypermucoviscous hvKp ESBL clinical strains of capsular serotype K1 and K2 isolated from pus and urine of critically ill patients in tertiary care hospitals in Oman. These strains belong to diverse sequence types (STs) namely ST-23(K1), ST-231(K2), ST-881(K2), and ST-14(K2). To study their virulence, Galleria Mellonella model and resistance to human serum killing were used. The G. mellonella model revealed that the K1/ST-23 isolate is the most virulent as 50% of the larvae died in the first day, followed by isolate K2/ST-231 for which 75% of the larvae died in the second day and isolate K2/ST-14 for which 50% of the larvae died in the second day. Resistance to human serum killing showed there was complete inhibition of bacterial growth of all four isolates by the end of the first hour and up to the third hour. For ST-23(K1), ST-231(K2), ST-881(K2), and ST-14(K2), there was a gradual decrease in the number of colonies up to the second to third hour, which indicates neutralization of bacterial cells by the serum components. However, this was followed by sudden increase of bacterial growth, indicating possible resistance of bacteria against human serum bactericidal activity. Whole genome sequencing (WGS) revealed that hvKp strains display unique genetic arrangement of k-loci. Whole genome SNP-based phylogenetic analysis revealed that these hvKp isolates are phylogenetically distinct, belonging to diverse clades and belong to different STs in comparison to global isolates. This is the first report from Oman detailing the WGS of hvKp clinical isolates and assessing their resistance and virulence genomics, which reinforce our understanding of their epidemiology and dissemination in clinical settings.

The aim was to study alterations of bacterial communities in patients undergoing hip or knee arthroplasty to assess the impact of chlorhexidine gluconate soap decolonisation and systemic antibiotic prophylaxis. A Swedish multicentre, prospective collection of samples obtained from elective arthroplasty patients (n=83) by swabbing anterior nares, skin sites in the groin and the site of planned surgery, before and after arthroplasty surgery, was analysed by 16S rRNA (V3-V4) gene sequencing and a complementary targeted tuf gene sequencing approach to comprehensively characterise alterations in staphylococcal communities. Significant reductions in alpha diversity was detected for both bacterial (p=0.04) and staphylococcal (p=0.03) groin communities after arthroplasty surgery with significant reductions in relative Corynebacterium (p=0.001) abundance and S. hominis (p=0.01) relative staphylococcal abundance. In nares, significant reductions occurred for S. hominis (p=0.02), S. haemolyticus (p=0.02), and S. pasteuri (p=0.003) relative to other staphylococci. S. aureus colonised 35% of anterior nares before and 26% after arthroplasty surgery. S. epidermidis was the most abundant staphylococcal species at all sampling sites. No bacterial genus or staphylococcal species increased significantly after arthroplasty surgery. Application of a targeted tuf gene sequencing approach provided auxiliary staphylococcal community profiles and allowed species-level characterisation directly from low biomass clinical samples.

Antimalarial drug resistance has become a real public health problem despite WHO measures. New sequencing technologies make it possible to investigate genomic variations associated with resistant phenotypes at the genome-wide scale. Based on the use of hemisynthetic nanopores, the PromethION technology from Oxford Nanopore Technologies can produce long-read sequences, in contrast to previous short-read technologies used as gold-standard to sequence Plasmodium. Two clones of P.falciparum (Pf3D7 and PfW2) were sequenced in long-read using the PromethION sequencer from Oxford Nanopore Technologies without genomic amplification. This made it possible to create a processing analysis pipeline for human Plasmodium with ONT Fastq only. De novo assembly revealed N50 lengths of 18,488 kb and 17,502 kb for the Pf3D7 and PfW2, respectively. The genome size was estimated at 23,235,407 base pairs for the Pf3D7 clone and 21,712,038 base pairs for the PfW2 clone. The average genome coverage depth was estimated at 787X and 653X for the Pf3D7 and PfW2 clones, respectively. This study proposes an assembly processing pipeline for the human Plasmodium genome using software’s adapted to large ONT data and the high AT percentage of Plasmodium. This search provides all the parameters, which have been optimized, for use with the software selected in the pipeline.

Investigation of bacterial communities is on the rise both in human and veterinary medicine. Their role in health maintenance and pathogenic mechanisms is in the limelight of infectious, metabolic, and cancer research. Among the most considered, gut bacterial communities takes the cake. Their part in animals was assessed mainly to improve animal production, public health, and pet management. In this regard, canaries deserve attention, being a popular pet and source of economic income for bird-keepers, for whom breeding represents a pivotal point. Thus, the aim of the present work was to follow gut bacterial communities’ evolution along on whole reproductive cycle of 12 healthy female canaries. Feces were collected during parental care, molting, and resting phase, and submitted for 16S rRNA sequencing. Data analysis a substantial presence of Lactobacillus aviarius along all the phases, and a relevant shift of microbiota during molting and rest due to an abrupt decrease of Vermiphilaceae family. Although the meaning of such change is not clear, future research may highlight unforeseen scenarios. Moreover, Lactobacillus aviarius may be deemed for normal bacteria flora restoration in debilitated birds, perhaps improving their health and productivity.

Background: Carbapenem-resistant Klebsiella pneumoniae (CRKP) is a serious public health issue. Aim of the study was to identify the antimicrobial resistance and accessory genes, the clonal relatedness and the evolution dynamics of selected CRKP isolates recovered in an adult and a pediatric intensive care unit of a tertiary hospital in Greece. Methods: Twenty-four CRKP isolates recovered during 2018-2022 were included in the study. Next-generation sequencing was performed using the Ion Torrent PGM Platform. The identification of the plasmid content, MLST and antimicrobial resistance genes, as well as the comparison of multiple genome alignments and the identification of core genome single-nucleotide polymorphism sites, were performed applying various bioinformatics software. Results: The isolates belonged to eight sequence types 11, 15, 30, 35, 39, 307, 323, and 512. A variety of carbapenemases (KPC, VIM, NDM and OXA-48) and resistance genes were detected. CRKP strains shared visually common genomic regions with the reference strain (NTUH-K2044). ST15, ST323, ST39, and ST11 CRKP isolates presented on average 17, 6, 16 and 866 recombined SNPs, respectively. All isolates belonging to ST15, ST323 and ST39 were classified in distinct phylogenetic branches, while ST11 isolates were assigned to a two-subclade branch. For large CRKP sets, the phylogeny seems to change approximately every seven SNPs. Conclusion: The current study provides insight into the genetic characterization of CRKP isolates in the ICUs of a tertiary hospital. Our results indicate clonal dispersion of ST15, ST323, and ST39 and highly diverged ST11 isolates.

Most existing cancer therapies negatively affect normal tissue as well as cancerous tissue. A potentially effective strategy for treating cancer that precludes off-target damage and could be an option for most patients would involve targeting one or more mutations that are ubiquitous in the given patient’s tumor(s). To effect this strategy, one would employ multi-region sequencing of a patient’s primary tumor and metastases to seek out mutations that are shared between all or at least most regions. Once the target or targets are known, one would ideally rapidly generate a molecular switch for at least one of said ubiquitous mutations that can distinguish the mutated DNA, RNA, or protein from the wild-type version and subsequently trigger a therapeutic response. I propose that the therapeutic response involve the replication of an oncolytic virus or intracellular bacterium, as any mutation can theoretically be detected by a vector that enters the cell - and automatic propagation could be very helpful. Moreover, the mutation “signal” can be easily enhanced through transcriptional and translational (if the target is an intracellular protein) enhancement. Importantly, RNA may make the best target for the molecular switches in terms of amplification of the signal and ease of targeting.

Abstract Successful treatment against infectious agents depends on rapid and accurate detection of the causative organisms. Misdiagnosis can hamper such success while leading to improper advising of antibiotics. In Bangladesh, the majority of the diagnostic centers detect and identify pathogens through culture and biochemical test-based methods and suggest antibiotics based solely on the results of disk-diffusion methods. This pilot study tried to validate the identity of the isolates characterized by diagnostic facilities near Dhaka. One hundred and twenty pre-characterized clinical isolates were analyzed biochemically and genotypically. Random Amplification of Polymorphic DNA-PCR, rcsA, and phoA genes-based PCR, and Loop-Mediated Isothermal Amplification (LAMP)-based identification of Klebsiella pneumoniae and Escherichia coli, respectively, followed by 16S rRNA sequencing confirmed misidentification of some clinical pathogens of other genera as Klebsiella spp. and E. coli. According to the antibiotic susceptibility testing guidelines, antibiotic choice, sensitivity pattern, and breakpoint measurement are different for each group of organisms. The lack of adherence to proper standards results in misdiagnosis and may facilitate the development of antibiotic resistance. The pilot study observers misidentification of clinical pathogens identified by the diagnostic centers. Well-characterized rapid molecular techniques like LAMP are suggested in clinical diagnosis to avoid misdiagnosis and subsequently circumvent antibiotic resistance development.

The holistic approach of One Health, which sees human, animal, plant, and environmental health as a unit, rather than discrete parts, requires not only interdisciplinary cooperation, but standardized methods for communicating and archiving data, enabling participants to easily share what they have learned and allow others to build upon their findings.Ongoing work by NCBI and the GenomeTrakr project illustrates how open data platforms can help meet the needs of federal and state regulators, public health laboratories, departments of agriculture, and universities. Here we describe how microbial pathogen surveillance can be transformed by having an open access database along with Best Practices for contributors to follow. First, we describe the open pathogen surveillance framework, hosted on the NCBI platform. We cover the current community standards for WGS quality, provide an SOP for assessing your own sequence quality and recommend QC thresholds for all submitters to follow. We then provide an overview of NCBI data submission along with step by step details. And finally, we provide curation guidance and an SOP for keeping your public data current within the database. These Best Practices can be models for other open data projects, thereby advancing the One Health goals of Findable, Accessible, Interoperable and Re-usable (FAIR) data.

The transition to reproduction is a crucial step in the life cycle of any organism. In Arabidopsis thaliana the establishment of reproductive growth can be divided into two phases: Firstly, cauline leaves with axillary meristems are formed and internode elongation begins. Secondly, lateral meristems develop into flowers with defined organs. Floral shoots are usually determinate and suppress the development of lateral shoots. Here, we describe a transposon insertion mutant in the Nossen accession with defects in floral development and growth. Most strikingly is the outgrowth of stems from the axillary bracts of the primary flower carrying secondary flowers. Therefore, we named this mutant flower-in-flower (fif). However, the transposon insertion in the annotated gene is not the cause for the fif phenotype. By means of classical and genome sequencing-based mapping, the mutation responsible for the fif phenotype was found to be in the LEAFY gene. The mutation, a G-to-A exchange in the second exon of LEAFY, creates a novel lfy allele and results in a cysteine-to-tyrosine exchange in the α1-helix of LEAFY´s DNA-binding domain. This exchange abolishes target DNA-binding, whereas subcellular localization and homomerization are not affected. To explain the strong fif phenotype against this molecular findings, several hypotheses are discussed.

Studies have previously described the occurrence of multidrug-resistant (MDR) Escherichia coli in human and veterinary medical settings, livestock, and, to a lesser extent, in the environment and food. While they mostly analyzed foodborne E. coli regarding phenotypic and sometimes genotypic antibiotic resistance and basic phylogenetic classification, we have limited understanding of the in vitro and in vivo virulence characteristics and global phylogenetic contexts of these bacteria. Here, we investigated an E. coli strain (PBIO3502) isolated from a pork sausage in Germany in 2021 in-depth. Whole-genome sequence analysis revealed sequence type (ST)58, which is an emerging international high-risk clonal lineage. In addition to its MDR phenotype that mostly matched the genotype, PBIO3502 demonstrated pronounced virulence features, including in vitro biofilm formation, siderophore secretion, serum resilience, and in vivo mortality in Galleria mellonella larvae. Along with the genomic analysis indicating close phylogenetic relatedness of our strain with publicly available, clinically relevant representatives of the same ST, these results suggest the zoonotic and pathogenic character of PBIO3502 with the potential to cause infection in humans and animals. Also, our study highlights the necessity of the One Health approach while integrating human, animal, and environmental health, as well as the role of meat products and food chains in the putative transmission of MDR pathogens.

We investigated and quantified bacterial accumulation in facemasks used by the general population, using 16S rRNA (Sanger Sequencing), culture and biochemical analysis along with Rose Bengal staining. Additionally, a systematic overview of the literature on face mask contamination was undertaken. We found an average bacterial load of 4.24 × 10⁴ CFU recovered/mask, with a maximum load of 2.85 × 10⁵ CFU. This maximum is 310 times higher than the limit value for contamination of ventilation system outlet surfaces specified by the German standard VDI 6022. Biochemical and molecular identification predominantly found Staphylococcus species (80%), including Staphylococcus aureus, along with endospore-forming Bacillus spp. Literature reports also indicate contamination of masks by bacterial and fungal opportunists of the genera Acinetobacter, Aspergillus, Alternaria, Bacillus, Cadosporium, Candida, Escherichia, Enterobacter, Enterococcus, Klebsiella (including K. pneumoniae), Micrococcus, Microsporum, Mucor, Pseudomonas, Staphylococcus and Streptococcus. Bacterial counts increase linearly with wearing duration, and prolonged use may affect the skin and respiratory microbiomes, promoting consequential eye, skin, oral and airway conditions. These aspects underscore the urgent need for further research and a risk-benefit analysis in respect of mask use, particularly given their unproven efficacy in disrupting the transmission of respiratory viruses and their adverse social consequences.

Background: SARS-CoV-2 Omicron JN.1 has the highest number of amino acid mutations in its receptor binding domain (RBD) and has acquired a hallmark L455S mutation. The immune evasion capability of JN.1 is a subject of scientific investigation. The US CDC used RT-qPCR SGTF as the proxy indicator of JN.1 infections for evaluation of the effectiveness of updated monovalent XBB.1.5 COVID-19 vaccines against JN.1 and recommended that all persons aged ≥6 months should receive an updated COVID-19 vaccine dose. Objective: Recommend Sanger sequencing for diagnosis of JN.1 infections to generate the data based on which guidelines are made to direct vaccination policies. Methods: The RNA in nasopharyngeal swab specimens from patients with clinical respiratory infection was subjected to nested RT-PCR, targeting a 398-base segment of the N-gene and a 445-base segment of the RBD of SARS-CoV-2 for amplification. The DNA sequences were analyzed for amino acid mutations. Results: The N-gene sequence showed Q229K, usually associated with Omicron BA.2.86 (+JN.1). The RBD sequence showed 24 amino acid mutations, including the hallmark L455S mutation for JN.1 and the V483del for BA.2.86 lineage. Conclusions: The CDC should consider Sanger sequencing of the RBD to diagnose JN.1 infections for statistical analysis in making vaccination policies.

Background: With continuous advances in next-generation sequencing (NGS), novel hepadnaviruses have been discovered in many species over the last 10 years. Methods: In this study, the cloud search and analysis of NGS data published in Nature by Edgar et al. was used as a basis to re-mine and reanalyze public NGS data for new hepadnaviruses. Results: Ultimately, at least 41 new species of hepadnaviruses were identified, including hepadnaviruses from the model animals hamster and mouse, frog hepadnavirus with pan-species infectivity, and diverse African cichlid hepadnaviruses that circulate within populations. Conclusions: The discovery of the new species of hepadnaviruses not only provides new clues for the study of the origin and evolution of hepadnaviruses, but also can be used to construct new hepadnavirus animal infection models, which will be helpful for the research of eradicative drugs for hepatitis B virus (HBV).

eDNA extracted from the gut contents of filter feeders with unique feeding habits can be used to identify biodiversity in aquatic ecosystems. In this study, we used eDNA from the gut contents of the clam Corbicula fluminea to examine estuarine ecosystem biodiversity. The field survey was conducted at three study sites in the Nakdong River Estuary, which is characterised by closed estuarine features resulting from the presence of an estuarine barrage. The collected C. fluminea samples were dissected to separate the gut contents, and the extracted eDNA was amplified using the 18S V9 primer targeting all eukaryotes. The amplified DNA was sequenced using next generation sequencing (NGS) techniques, and BLASTn was performed based on the NCBI database. We obtained 21 unique operational taxonomic units (OTUs), including fish (approximately 9.52%), copepods (approximately 14.29%), and green algae (approximately 23.81%), which represented a wide range of habitats such as marine, brackish, freshwater, and terrestrial environments. These results suggest that various organisms living in aquatic ecosystems can be identified through eDNA from the gut contents of C. fluminea and support the potential application of bivalves as eDNA samplers in diverse aquatic ecosystems.

In the post antimicrobial era, increasing attention is paid towards using bacteriophage (phage in short) therapy to control antibiotic-resistant bacteria. The first step in phage therapy applications is isolating highly efficient lytic phages or phage cocktails from various sources. When a double-layer- agar with around 0.7% agar in top agar is employed, it results in a low number of phage isolation with a poor resolution, and in many cases, you miss the phage. To address this problem, a low concentration of agar in top agar is examined for better phage isolation. Here, our results proved the efficiency of isolating phage upon formulating a double-layer agar with 0.3% agar in top agar. A sewage sample was collected then phages were isolated, purified, and spotted on a top layer agar with 0.3% agar. The results showed the possibility of isolating a higher number of phages on 0.3% top agar than 0.7%. The finding advocates using 0.3% top agar for the double-layer agar, as it will provide fast, better, and easy phage screening and isolation.

Next Generation Sequencing based cancer risk screening with multigene panels has become the most successful method for programming cancer prevention strategies. ATM germ-line heterozygosity has been described as able to increase tumor susceptibility. In particular, families that carry heterozygous germ-line variants of ATM gene show a 5- to 9-fold risk of developing breast cancer. Recent studies identified ATM as the second most mutated gene after CHEK2 in BRCA-negative patients. Nowadays, more than 170 potential missense variants and several truncating mutations have been identified in ATM gene. Here we present the molecular characterization a new ATM deletion, identified thanks to the CNV algorithm implemented in the NGS analysis pipeline. An automated workflow implementing the SOPHiA Genetics’ Hereditary Cancer Solution (HCS) protocol was used to generate NGS libraries that were sequenced on Illumina MiSeq Platform. NGS data analysis allowed to identify a new inactivating deletion of exons 19-27 of ATM gene. DNA breakpoint was characterized both at DNA and RNA level

Directed evolution methods are becoming increasingly popular, as they are extremely powerful toward developing new biomolecules with altered/novel activities, e.g., proteins with new catalytic functions or substrate specificities, and nucleic acids that recognize an intended target. Especially useful are systems that have incorporated continuous evolution, where the protein to be evolved undergoes continuous mutagenesis to evolve a desired trait with little to no input from the researcher once the system is started. However, continuous evolution methods can be challenging to implement in the lab and daunting for researchers to invest time and resources. Our intent is to provide basic information and helpful suggestions that we have gained from our experience with bacterial phage-assisted continuous evolution (PACE). Specifically, we review factors to consider before adopting PACE for a given evolution scheme, different types of selection circuits that can be utilized with particular focus on the PACE-B1H selection system, what optimization of a PACE selection circuit may look like using directed evolution of ME47 as a case study, and additional techniques that may be incorporated into a PACE experiment. With this information, researchers will be better equipped to determine if PACE is a valid strategy to use to evolve their proteins and how to set up a valid selection circuit.

Liquid biopsies are an integral part of the diagnosis of cancer. Here, we have extended previous validation studies of a new targeted NGS panel to include the detection of copy number amplifications (CNAs), fusions, and exon skipping variants. Detection of these gene classes included specimens from clinical and healthy donors and cell lines (fusions: ROS1, EML4-ALK, NTRK1), (exon skipping: MET exon 14), and (CNAs: HER2, CDK6, EGFR, MYC and MET). The limit of detection (LOD) for fusion/skipping was 42 copies (QC threshold = 3 copies) and was verified using 3 additional fusion/skipping variants. LOD for CNAs was 1.40-fold-change (QC threshold = 1.15-fold change) and was verified with 3 additional CNAs. In precision studies, all fusion/skipping variants were detected in all runs and all days of testing (n=18/18; 100%); average CV for repeatability was 20.5% (range 8.7% - 34.8%), and for reproducibility was 20.8% (range 15.7% - 30.5%). For CNAs, 28/29 (96.6%) copy gains were detected. For CNAs, average CV of 1.85% (range 0% to 5.49%) for repeatability and 6.59% (range 1.65% to 9.22%) for reproducibility. The test panel meets the criteria for being highly sensitive and specific and extends its utility for the detection of clinically relevant variants in cancer.

Viral integration within the host genome plays a pivotal role in carcinogenesis. Various disruptive mechanisms are involved, leading to genomic instability, mutations, and DNA damage. With next-generation sequencing (NGS), we can now precisely identify viral and host genomic breakpoints and chimeric sequences, which are useful for integration site analysis. In this study, we evaluated a commercial hybrid-capture NGS panel specifically designed for detecting three key viruses: HPV, HBV, and HIV-1. We also tested workflows for Viral Hybrid Capture (VHC) and Viral Integration Site (VIS) analysis, leveraging customized viral databases in CLC Microbial Ge-nomics. By analyzing sequenced data from virally infected cancer cell lines (including SiHa, HeLa, CaSki, C-33A, DoTc2, 2A3, SCC154 for HPV; 3B2, SNU-182 for HBV; and ACH-2 for HIV-1), we precisely pinpointed viral integration sites. The workflow also highlighted disrupted and neigh-boring human genes that may play a crucial role in tumor development. Our results included in-formative virus-host read mappings, genomic breakpoints, and integration circular plots. These visual representations enhance our understanding of the integration process. In conclusion, our seamless end-to-end workflow bridges the gap in understanding viral contributions to cancer development, paving the way for improved diagnostics and treatment strategies.

The rapid expansion of transcriptomics from increased affordability of next-generation sequencing (NGS) technologies generates rocketing amounts of gene expression data across biology and medicine, and notably in cancer research. Concomitantly, many bioinformatics tools were developed to streamline gene expression analysis and quantification. We tested the concordance of NGS RNA sequencing (RNA-seq) analysis outcomes between the two predominant programs for reads alignment, HISAT2 and STAR, and the two most popular programs for quantifying gene expression in NGS experiments, edgeR and DESeq2, using RNA-seq data from a series of breast cancer progression specimens, which include histologically confirmed normal, early neoplasia, ductal carcinoma in situ and infiltrating ductal carcinoma samples microdissected from formalin fixed, paraffin embedded (FFPE) breast tissue blocks. We identified significant differences in aligners’ performance: HISAT2 was prone to misalign reads to retrogene genomic loci, STAR generated more precise alignments, especially for early neoplasia samples. edgeR and DESeq2 produced similar lists of differentially expressed genes in stage comparisons, with edgeR producing more conservative, though shorter, lists of genes. Albeit, Gene Ontology (GO) enrichment analysis revealed no skewness in significant GO categories identified among differentially expressed genes by edgeR vs DESeq2. As transcriptome analysis of archived FFPE samples becomes a vanguard of precision medicine, identification and fine-tuning of bioinformatics tools becomes critical for clinical research. Our results indicate that STAR and edgeR are well-suited tools for differential gene expression analysis from FFPE samples.

The study examines the integration of postbiotics in food products through the use of attenuated probiotics, specifically lactic acid bacteria (LAB) in bread. Postbiotics, non-viable microorganisms or their metabolites, offer health benefits similar to probiotics without the risks associated with live bacteria. This research evaluates the regulatory aspects and safety of LAB in sourdough bread production, highlighting their historical and significant use in Europe before 1997. The study includes microbial quantification and Next-Generation Sequencing (NGS) to identify LAB in traditional sourdough, comparing them with historical and current EFSA Qualified Presumption of Safety (QPS) lists. Findings show that the LAB present in sourdough have been extensively and safely used in bread making, supporting their classification as non-novel foods under EU regulations. The stability and consistency of LAB metabolites in sourdough bread are also confirmed, ensuring quality and safety in each batch. The study concludes that LAB in sourdough, when inactivated through bread-making processes, are not considered novel foods, aligning with historical, scientific, and regulatory evidence.

Disorder of S-adenosylhomocysteine hydrolase (AHCY) activity leads to the potentially lethal rare disease AHCY deficiency, first described in 2004 by Baric and co-workers [1]. In order to shed new light on molecular aspects of the disease, in particular changes at transcriptome level, we enabled knock-down of AHCY expression in model systems such as the colon cancer cell line SW480 to simulate the environment occurring in AHCY deficiency patients. Further, we per-formed deep sequencing of mRNA, followed by differential expression and molecular pathway analysis. Fifteen differentially expressed networks were identified, and interestingly, we found a predicted effect of AHCY down-regulation on the expression of the Lymphoidphoid enhanc-er-binding factor 1 (LEF1) gene, indicating changes in the TCF4/LEF1 complex. LEF1, a member of the T-cell Factor (TCF)/LEF1 family of high-mobility group transcription factors, is a down-stream mediator of the Wnt/β-catenin signaling pathway [2]. LEF1 is essential in stem cell maintenance, and especially in its role in epithelial-mesenchymal transition (EMT). Western blot analysis of LEF1 protein expression confirmed our transcriptomic data predictions and revealed significantly increased LEF1 protein in AHCY- deficient cells, providing a novel link between AHCY and cancer cell phenotype.

Hepatocellular carcinoma (HCC) is one of the most common cancers and the main cause of cancer-related death globally. Immune dysregulation of CD4+ T cells has been identified as a role in the development of HCC. Nevertheless, the underlying molecular pathways of CD4+ T cells in HCC are not completely known. Thus, a better understanding of the dysregulation of lncRNA-miRNA/mRNA network might yield novel insights into the etiology or progression of HCC. In this study, circulating CD4+ T cells were isolated from the whole blood of 10 healthy controls and 10 HCC patients for next-generation sequencing of the expression of lncRNAs, miRNAs, and mRNAs. Our data showed that there was different expression of 34 transcripts (two lncRNAs, XIST and MIR222HG, 29 mRNA, and 3 other types of RNA) and 13 miRNAs in the circulating CD4+ T cells of HCC patients. The expression of the lncRNA XIST-related miRNAs and their target mRNAs was confirmed using real-time quantitative polymerase chain reaction (qPCR) on samples from 100 healthy controls and 60 HCC patients. The lncRNA–miRNA/mRNA regulation network was created using interaction data generated from ENCORI and revealed there are positive correlations in the infiltration of total CD4+ T cells and negative correlations in the infiltration of Th1 CD4+ T cells.

The majority of Candida species are known as non-pathogenic yeasts and rarely involved in human diseases. However, recently case reports of human infections caused by non-albicans Candida species have increased, mostly in immunocompromised hosts. Our study aimed to describe and caracterise as thoroughly as possible, a new species of the Candida genus, named here Candida massiliensis (PMML0037), isolated from a clinical sample of human sputum. We compared genetic data based on the sequences of four genetic regions: "Internal Transcribed Spacers" of rRNA, D1/D2 domains (28S large subunit rRNA) and part of the genes encoding Translation Elongation Factor 1-α and β-tubulin2, to morphological characters, from scanning electron microscopy (TM 4000 Plus, SU5000), physiological, including the results of oxidation and assimilation tests of different carbon sources by the Biolog system, and chemical mapping by Energy-Dispersive X-ray Spectroscopy. Lastly, the in vitro antifungal susceptibility profile was performed using the E-testTM exponential gradient method. The multilocus analysis supported the genetic position of Candida massiliensis (PMML0037) as a new species of the genus Candida, and the phenotypic analysis highlighted its unique morphological and chemical profile when compared to other Candida species included in the study.

Polymyxin resistance, determined by mcr genes located on plasmid DNA, currently pose a high epidemiological threat. Non-typhoid Salmonella (NTS) are one of the key pathogens causing diarrheal diseases. Here, we report the isolation and whole genome sequencing of multidrug colistin-resistant/susceptible isolates of non-typhoid Salmonella enterica serovars carries mcr genes. Non-typhoid strains of Salmonella enterica subsp. enterica were isolated during microbiological monitoring of the environment, food, and diarrheal disease patients between 2018 and 2020 in Russia (n=586). mcr-1 genes were detected using a previously developed qPCR assay and whole genome sequencing of mcr positive isolates was performed by both short-read (Illumina) and long-read (Oxford Nanopore) approaches. Three colistin-resistant isolates including two isolates of S. Enteritidis and one isolate of S. Bovismorbificans carried the mcr-1.1 gene located on IncX4 and IncI2 conjugative plasmids, respectively. The phenotypically colistin-susceptible isolate of S. Typhimurium carried a mcr-9 gene on plasmid IncHI2. In conclusion, we present the first three cases of mcr gene carrying NTS isolates detected in Russia with both outbreak and sporadic epidemiological background.

A growing number of emerging SARS-CoV-2 variants is being identified worldwide, potentially impacting the effectiveness of current vaccines. We report the data obtained in several Italian regions involved in the SARS-CoV-2 variant monitoring from the beginning of the epidemic and spanning the period from October 2020 to March 2021.

One in eight women will be diagnosed with breast cancer (BC) in their lifetime, resulting in over 2 million cases annually. BC is the most common cancer among women. Unfortunately, the etiology of majority of cases remains unknown. Recently, evidence has shown that the human microbiota plays an important role in health and disease. Intriguingly, studies have revealed the presence of microorganisms in human breast tissue, which was previously presumed to be sterile. Next-generation sequencing technologies have paved way for the investigation of breast microbiota, uncovering bacterial signatures that are associated with BC. Some of the bacterial species were found to possess pro-carcinogenic and/or anti-carcinogenic properties, suggesting that the breast microbiota has potentially crucial roles in maintenance of breast health. In this review, we summarize the recent findings on breast tissue microbiota and its interplay with BC. Bacterial signatures identified via next-generation sequencing as well as their impact on breast carcinogenesis and cancer therapies are reviewed. Correlation of breast tissue microbiota and other factors, such as geographical and racial differences, in BC is discussed. Additionally, we discuss the future directions of research on breast microbiota as well as its potential role in prevention, diagnosis and treatment of BC.

A genome-wide association study (GWAS) was performed on a set of 260 lines which belong to three different bi-parental flax mapping populations. These lines were sequenced to an averaged genome coverage of 19× using the Illumina Hi-Seq platform. Phenotypic data for 11 seed yield and oil quality traits were collected in eight year/location environments. A total of 17,288 single nucleotide polymorphisms were identified, which explained more than 80% of the phenotypic variation for days to maturity (DTM), iodine value (IOD), palmitic (PAL), stearic, linoleic (LIO) and linolenic (LIN) acid contents. Twenty-three unique genomic regions associated with 33 QTL for the studied traits were detected, thereby validating four genomic regions previously identified. The 33 QTL explained 48-73% of the phenotypic variation for oil content, IOD, PAL, LIO and LIN but only 8-14% for plant height, DTM and seed yield. A genome-wide selective sweep scan for selection signatures detected 114 genomic regions that accounted for 7.82% of the flax pseudomolecule and overlapped with the 11 GWAS-detected genomic regions associated with 18 QTL for 11 traits. The results demonstrate the utility of GWAS combined with selection signatures for dissection of the genetic structure of traits and for pinpointing genomic regions for breeding improvement.

Transposable elements (TEs) are repetitive elements which make up around 45% of the human genome. A class of TEs known as SINE-VNTR-Alu (SVA) demonstrate the capacity to mobilise throughout the genome, resulting in SVA polymorphisms for presence or absence within the population. Although studies have previously highlighted the involvement of TEs within neurodegenerative diseases, such as Parkinson’s disease and amyotrophic lateral sclerosis (ALS), however the exact mechanism has yet to be identified. In this study we used whole genome sequencing and RNA sequencing data of ALS patients and healthy controls from the New York Genome Center ALS Consortium, to elucidate the influence of reference SVA elements on gene expression genome-wide within central nervous system (CNS) tissues. To investigate this, we applied matrix expression quantitative trait loci analysis and demonstrated that reference SVA insertion polymorphisms can significantly modulate the expression of numerous genes, preferentially in the trans position, and in a tissue-specific manner. We also highlight that SVAs significantly regulate mitochondrial genes as well as genes within the HLA and MAPT loci, previously associated within neurodegenerative disease. In conclusion, this study continues to bring to light the effects of polymorphic SVAs on gene regulation and further highlights the importance of TEs within disease pathology.

Zinc is an important nutrient for plant growth and development. Its availability is influenced by zinc solubilizing microbes (ZSM). Effects of commonly promoted agronomic practices on the abundance of ZSM are so far not well understood. In this study, we assessed the effects of farmyard manure (FYM) application, either sole, or in combination with residue and/or inorganic fertilizer inputs on ZSM community structure using 11 treatments in a long-term (17 years) integrated soil fertility management experiment located in western Kenya, The study was conducted in 2019. Bacterial and fungal community composition (at 0-15 cm depth) was evaluated by amplicon sequencing on an Illumina Miseq platform following the manufacturer’s guidelines. Microbial diversity indices (Shannon, Simpson and Chao1) and statistical analyses were computed using R project packages. Putative ZSM (i.e., the ZSM generally considered to possess the zinc solubilizing capabilities) were clustered in two major clades based on either application or no application of FYM. Sole application of FYM significantly (P<0.05) increased the abundance of several ZSM under maize-tephrosia rotation. In addition, systems with combined application of FYM with other inputs generally showed significantly increasing trends for some ZSM under maize-tephrosia rotation. Moreover, combined application of FYM and P rather than P only significantly increased abundance of some ZSM under maize-monocropping systems. Furthermore, besides affecting ZSM abundance, soil chemical variables involving soil organic carbon (SOC), total N and Olsen P significantly increased with FYM application. This study indicated that management practices like application of FYM that increase SOC, and other soil chemical parameters, also/concomitantly increase ZSM abundance; implying enhanced capacities for microbial-linked zinc availability.

Cyclic vomiting syndrome (CVS) is a frequently disabling condition defined by severe, discrete, stereotypical episodes of nausea and vomiting. In our previous report [1], a candidate gene list was evaluated in a group of 80 unrelated participants with CVS with extensive genetic sequencing (whole exome/genome). Twenty-two good candidate genes for CVS association were identified, from which a cellular model of aberrant cation transport and energy metabolism was proposed. Herein, we present additional clinical information on those 80 participants, including eight case reports to highlight general principles and clinical practices. Our participants demonstrate substantial phenotypic and genotypic heterogeneity in CVS, including a propensity for specific co-morbidities in the participants and relatives, evolution towards episodes in which vomiting is less apparent, labeling with “non-organic” diagnoses, multiple genes, polygenic inheritance, and association with variants not designated as Pathogenic. Our case reports illustrate how genetic information can guide clinical management and argue to the high clinical utility of appropriate genetic testing and analysis in CVS. Expression profiling of our 22 candidate genes suggests an anatomical model in a unit defined by vagal afferents and adjacent cells. Our models are consistent with multiple current hypotheses of CVS and the efficacy of commonly-recommended treatments.

Until the 1950s, developmental dyslexia was defined as a hereditary visual disability, selectively affecting reading without compromising oral or non-verbal reasoning skills. This changed radically after the development of the phonological theory of dyslexia; this not only ruled out any role for visual processing in its aetiology, but also cast doubt on the use of discrepancy between reading and reasoning skills as a criterion for diagnosing it. Here I argue that this theory is set at too high a cognitive level to be explanatory; we need to understand the pathophysiological visual and auditory mechanisms that cause children’s phonological problems. I discuss how the ‘magnocellular theory’ attempts to do this in terms of slowed and error prone temporal processing which leads to dyslexics’ defective visual and auditory sequencing when attempting to read. I attempt to deal with the criticisms of this theory and show how it leads to a number of successful ways of helping dyslexic children to overcome their reading difficulties.

Both SARS-CoV-2 and SARS-CoV initially appeared in China and spread to other parts of the world. SARS-CoV-2 has generated a COVID-19 pandemic causing more than 6 million human deaths worldwide while the SARS outbreak quickly ended in six months with a global total of 774 reported deaths. One of the factors contributing to this stunning difference in the outcome between these two outbreaks is the inaccuracy of the RT-qPCR tests for SARS-CoV-2, which generated a large number of false-negative and false-positive test results that have misled patient management and public health policy-makers. This article presented Sanger sequencing evidence to show that the RT-PCR diagnostic protocol established in 2003 for SARS-CoV can in fact detect SARS-CoV-2 accurately due to the well-known nonspecific PCR amplification of DNAs with similar sequences. Using nested RT-PCR followed by Sanger sequencing to retest 50 patient samples collected in January, 2022 and sold as RT-qPCR positive reference confirmed 21 (42%) were false-positive. Although the other 29 positive isolates were categorized as Omicron variant by partial sequencing of the N gene, and the RBD and the NTD of the S gene, 9 (31%) showed focal to complete sequencing failure in the S gene segments due to multi-allelic SNPs. During the course of the study, an Omicron variant isolate containing a BA.1 NTD and a BA.2 RBD in its S gene was also detected. Routine partial S gene sequencing of all PCR-positive samples can timely discover multi-allelic SNPs and viral recombination in the circulating variants for investigation of their impacts on vaccine efficacies, therapeutics and diagnostics.

The field of mitochondrial genomics has advanced rapidly and revolutionized disciplines from molecular anthropology, population genetics, to medical/oncogenetics. However, mtDNA next-generation sequencing (NGS) analysis for matrilineal haplotyping and phylogeographic inference remains hindered by the lack of a consolidated, mitogenome database and efficient bioinformatics pipeline. To address this, we developed a customized human mitogenome database (hMITO DB) embedded in a CLC Genomics workflow for read mapping, variant analysis, haplotyping, and geo-mapping. The database was constructed from 4,286 mitogenomes. The macro-haplogroup (A to Z) distribution and representative phylogenetic tree were found consistent with published literature. The hMITO DB automated workflow was tested using mtDNA-NGS sequences derived from Pap smears and cervical cancer cell lines. The auto-generated read mapping, variants track, and table of haplotypes and geo-origins were completed in 15 min for 47 samples. The mtDNA workflow proved to be a rapid, efficient and accurate means of sequence analysis for translational mitogenomics.

Primary varicella-zoster virus (VZV) infection causes varicella (chickenpox) and the establishment of a lifelong latent infection in ganglionic neurons. VZV reactivates in about one-third of infected individuals to cause herpes zoster, often accompanied by neurological complications. The restricted host range of VZV and, until recently, the lack of suitable in vitro models to study VZV latency have seriously hampered molecular studies of viral latency. Nevertheless, recent technological advances facilitated a series of exciting studies that resulted in the discovery of a VZV latency-associated transcript (VLT) and have redefined our understanding of VZV latency and factors that initiate reactivation. Together, these findings pave the way for a new era of research that may finally unravel the precise molecular mechanisms that govern latency. In this review, we will summarize the implications of recent discoveries in the VZV latency field from both a virus and host perspective and provide a roadmap for future studies.

As SARS-CoV-2 continues to spread among human populations, genetic changes occur and accumulate in the circulating virus. Some of these genetic changes have caused amino acid mutations, including deletions, which may have potential impact on critical SARS-CoV-2 countermeasures, including vaccines, therapeutics, and diagnostics. Considerable efforts have been made to categorize the amino acid mutations of the angiotensin-converting enzyme 2 (ACE2) receptor binding domain (RBD) of the spike (S) protein along with certain mutations in other regions within the S protein as specific variants in an attempt to study the relationship between these mutations and the biological behavior of the virus. However, the currently used whole genome sequencing surveillance technologies can test only a small fraction of the positive specimens with high viral loads and often generate uncertainties in nucleic acid sequencing that needs additional verification for precision determination of mutations. This article introduces a generic protocol to routinely sequence a 437-bp nested RT-PCR cDNA amplicon of the ACE2 RBD and a 490-bp nested RT-PCR cDNA amplicon of the N-terminal domain (NTD) of the S gene for detection of the amino acid mutations needed for accurate determination of all variants of concern and variants of interest according to the definitions published by the U.S. Centers for Disease Control and Prevention. This protocol was able to amplify both nucleic acid targets into cDNA amplicons to be used as templates for Sanger sequencing on all 16 clinical specimens that were positive for SARS-CoV-2.

The observation that certain therapeutic strategies for targeting inflammation benefit patients with distinct immune-mediated inflammatory diseases (IMIDs) is exemplified by the success of TNF blockade in conditions including rheumatoid arthritis, ulcerative colitis and skin psoriasis, albeit only for subsets of individuals with each condition. This suggests intersecting "nodes" in inflammatory networks at a molecular and cellular level may drive and/or maintain IMIDs, being “shared” between traditionally distinct diagnoses without mapping neatly to a single clinical phenotype. In line with this proposition, integrative tumour tissue analyses in oncology have highlighted novel cell states acting across diverse cancers, with important implications for precision medicine. Drawing upon advances in the oncology field, this narrative review will first summarise learnings from the Human Cell Atlas in health as a platform for interrogating IMID tissues. It will then review cross-disease studies to date that inform this endeavour, before considering future directions in the field.

The infections caused by various bacterial pathogens both in clinical and community settings represent a significant threat to public healthcare worldwide. The growing resistance to antimicrobial drugs acquired by bacterial species causing healthcare-associated infections has already become a life-threatening danger noticed by the World Health Organization. Several groups or lineages of bacterial isolates usually called ‘the clones of high risk’ often drive the spread of resistance within particular species. Thus, it is vitally important to reveal and track the spread of such clones and the mechanisms by which they acquire antibiotic resistance and enhance their survival skills. Currently, the analysis of whole genome sequences for bacterial isolates of interest is increasingly used for these purposes, including epidemiological surveillance and developing of spread prevention measures. However, the availability and uniformity of the data derived from the genomic sequences often represents a bottleneck for such investigations. In this dataset, we present the results of a comprehensive genomic epidemiology analysis of 17,546 genomes of a dangerous bacterial pathogen Acinetobacter baumannii. Important typing information including multilocus sequence typing (MLST)-based sequence types (STs), intrinsic blaOXA-51-like gene variants, capsular (KL) and oligosaccharide (OCL) types, CRISPR-Cas systems, and cgMLST profiles are presented, as well as the assignment of particular isolates to nine known international clones of high risk. The presence of antimicrobial resistance genes within the genomes is also reported. These data will be useful for researchers in the field of A. baumannii genomic epidemiology, resistance analysis and prevention measure development.

Newcastle disease virus (NDV) genotype VII is a highly pathogenic Orthoavulavirus that has caused multiple outbreaks among poultry in Egypt since 2011. This study aimed to investigate the genetic diversity of NDV prevailing in domestic and wild birds in Egyptian governorates. A total of 37 oropharyngeal swabs from wild birds and 101 swabs from domestic bird flocks including chickens, ducks, turkeys, and swans were collected from different geographic regions within 13 governorates during 2019-2020. Virus isolation and propagation via embryonated eggs revealed 91 swab samples produced allantoic fluid containing hemagglutination activity, suggestive of virus presence. The use of RT-PCR targeted to F gene successfully detected NDV in 85 samples. The geographical prevalence of NDV spread to 12 governorates in domestic birds, migratory and non-migratory wild birds. Following whole genome sequencing, we assembled six NDV genome sequences (70 - 99% of genome coverage), including five full F gene sequences. All NDV strains carried high virulence, based on the presence of polybasic amino acids (RRQRF) at the F gene cleavage site. Phylogenetic analysis revealed that the NDV strains belonged to class II within genotype VII.1.1. The presence of genetically similar virulent NDV in wild birds further highlights their role in the dissemination of NDV in poultry populations across Egypt. Continued genomic surveillance in both wild birds and poultry would be necessary for monitoring NDV incursions and genetic diversification.

Pasmo (Septoria linicola) is a fungal disease causing major losses in seed yield and quality, and stem fibre quality in flax. Pasmo resistance (PR) is quantitative and has low heritability. To improve PR breeding efficiency, the accuracy of genomic prediction (GP) was evaluated using a diverse worldwide core collection of 370 accessions. Four marker sets, including three defined by 500, 134, and 67 previously identified quantitative trait loci (QTL) and one of 52,347 PR-correlated genome-wide single nucleotide polymorphisms, were used to build ridge regression best linear unbiased prediction (RR-BLUP) models using pasmo severity (PS) data collected from field experiments performed during five consecutive years. With five-fold random cross-validation, GP accuracy as high as 0.92 was obtained from the models using the 500 QTL when the average PS was used as the training dataset. GP accuracy increased with training population size, reaching values >0.9 with training population size greater than 185. Linear regression of the observed PS with the number of positive-effect QTL in accessions provided an alternative GP approach with an accuracy of 0.86. The results demonstrate the GP models based on marker information from all identified QTL and the 5-year PS average is highly effective for PR prediction.

Recent evidence indicates that the risk of being diagnosed with cancer in a tissue is strongly correlated (0.80) with the number of stem cell divisions accumulated by the tissue. Since cell division can generate random mutations during DNA replication, this correlation has been used to propose that cancer is largely caused by the accumulation of unavoidable mutations in driver genes. However, no correlation between the number of gene mutations and cancer risk across tissues has been reported. Because many somatic mutations in cancers originate prior to tumor initiation and the number of cell divisions occurring during tumor growth is similar among tissues, here I use whole genome sequencing information from 22,086 cancer samples and incidence data from the largest cancer registry in each continent to study the relationship between the number of gene mutations and the risk of cancer across 33 tissue types. Results show a weak positive correlation (mean = 0.14) between these two parameters in each of the five cancer registries. The correlation became stronger (mean = 0.50) when gender-related cancers were excluded. Results also show that 1,003 samples from 29 cancer types have zero mutations in genes. These data suggest that cancer etiology can be better explained by the accumulation of stem cell divisions than by the accumulation of gene mutations. Possible mechanisms by which the accumulation of cell divisions in stem cells increases the risk of cancer are discussed.

Homologous recombination deficiency (HRD) has become an important prognostic and predictive biomarker for patients with high-grade serous ovarian cancer who may benefit from poly-ADP ribose polymerase inhibitors (PARPi) and platinum-based therapies. HRD testing provides relevant information to personalize patients’ treatment options, and it has been progressively incorporated in diagnostic laboratories. In this study, we evaluated the feasibility of in-house HRD testing deployable in a diagnostic clinical setting, comparing results from two commercially available next-generation sequencing (NGS)-based tumor tests (SOPHiA DDMTM HRD Solution and AmoyDx HRD Focus Panel) with the reference assay from Myriad MyChoice CDx. A total of 85 ovarian cancer samples were subjected to HRD testing. An overall strong correlation was observed across the three assays evaluated, regardless of the different underlying methods employed to assess genomic instability, with the highest pairwise correlation between Myriad and SOPHiA (R=0.87, p-value=3.39x10-19). The comparison of the assigned HRD status to the reference Myriad’s test revealed a positive predictive value (PPV) and negative predictive value (NPV) of 90.9% and 96.3% for SOPHiA’s test, while AmoyDx’s test achieved 75% PPV and 100% NPV. This is the largest HRD testing evaluation using different methodologies and provides a clear picture of the robustness of NGS-based tests currently offered in the market. Our data show that the implementation of in-house HRD testing in diagnostic laboratories is technically feasible and it can be reliably performed with commercial assays. Also, the turnaround time is compatible with clinical needs, being an ideal alternative to offer to a broader number of patients, while still maintaining high-quality standards at more accessible price tiers.

Blackcurrant leaf chlorosis associated virus (BCLCaV) was detected recently by next-generation sequencing (NGS) and proposed as a new and distinct species in the genus Idaeovirus. Genomic components of BCLCaV that were detected and confirmed include: 1) RNA-1 that is monocistronic and encodes the replicase complex; 2) a bicistronic RNA-2 that encodes a movement protein (MP) and the coat protein (CP) of the virus, with open reading frames (ORF) that overlap by a single adenine (A) nucleotide (nt) representing the third position of an opal stop codon of the MP ORF2a and the first position of the start codon of the CP ORF2b; 3) a subgenomic form of RNA-2 (RNA-3) that contains ORF2b; and 4) a concatenated form of RNA-2 that consists of a complementary and inverted RNA-3 conjoined to the full-length RNA-2. Analysis of NGS-derived paired-end reads revealed the existence of bridge reads encompassing the 3’-terminus and 5’-terminus of RNA-2 or RNA-3 of BCLCaV. The full RNA-2 or RNA-3 could be amplified using outward facing or abutting primers; also RNA-2/RNA-3 could be detected even after three consecutive RNase R enzyme treatments with denaturation at 95 ^oC preceding each digestion. Evidence was obtained indicating that there are circular forms of BCLCaV RNA-2 and RNA-3.

Escherichia coli is an indicator microorganism in One Health antibiotic resistance surveillance programmes. The purpose of the study was to describe and compare E. coli isolates obtained from pigs and human contacts from a commercial farm in South Africa using conventional methods and whole genome sequencing (WGS). Porcine E. coli isolates were proportionally more resistant phenotypically and harboured a richer diversity of antibiotic resistance genes, as compared to human E. coli isolates. Different pathovars, namely ExPEC (12.43%, 21/169), ETEC (4.14%; 7/169), EPEC (2.96%, 5/169), EAEC (2.96%, 5/169) and STEC (1.18%, 2/169) were detected at low frequencies. Sequence type complex (STc) 10 were the most prevalent (85.51%, 59/169) among human and porcine isolates. Six STcs (STc10, STc86, STc168, STc206, STc278 and STc469) were shared at the human-livestock interface according to multilocus sequence typing (MLST). Core genome MLST and hierarchical clustering (HC) showed that human and porcine isolates were overall genetically diverse, but some clustering at HC2 to HC200 were observed. In conclusion, even though the isolates shared a spatio-temporal relationship, there were still differences in the virulence potential, antibiotic resistance profiles, and cgMLST and HC according to the source of isolation.

Fear extinction requires coordinated neural activity within the amygdala and medial prefrontal cortex (mPFC). Any behavior has a transcriptomic signature that is modified by environmental experiences, and specific genes are involved in functional plasticity and synaptic wiring during fear extinction. Here, we investigated the effects of optogenetic manipulations of prelimbic (PrL) pyramidal neurons on amygdala gene expression to analyze the specific transcriptional pathways involved in adaptive and maladaptive fear extinction. To this aim, transgenic mice were (or not) fear-conditioned and during the extinction phase they received optogenetic (or sham) stimulations over PrL pyramidal neurons. At the end of behavioral testing, electrophysiological (neural cellular excitability and Excitatory Post-Synaptic Currents) and morphological (spinogenesis) correlates were evaluated in the PrL pyramidal neurons. Furthermore, transcriptomic cell-specific RNA-analyses (differential gene expression profiling and functional enrichment analyses) were performed in amygdala pyramidal neurons. Our results show that the optogenetic activation of PrL pyramidal neurons in fear-conditioned mice induces fear extinction deficits, reflected in an increase of cellular excitability, excitatory neurotransmission, and spinogenesis of PrL pyramidal neurons, and in strong modifications of the transcriptome of amygdala pyramidal neurons. Understanding the electrophysiological, morphological and transcriptomic architecture of fear extinction may facilitate the comprehension of fear-related disorders.

The vast majority of patients with pancreatic ductal adenocarcinoma (PDAC) suffer cachexia. Although cachexia results from concurrent loss of adipose and muscle tissue, most studies focus on muscle alone. Emerging data demonstrate the prognostic value of fat loss in cachexia. Here we sought to identify the muscle and adipose gene profiles and pathways regulated in cachexia. Matched rectus abdominis muscle and subcutaneous adipose tissue were obtained at surgery from patients with benign conditions (n=11) and patients with PDAC (n=24). Self-reported weight loss and body composition measurements defined cachexia status. Gene profiling was done using Ion proton sequencing. Results were queried against external datasets for validation. 961 DE genes were identified from muscle and 2000 from adipose tissue, demonstrating greater response of adipose than muscle. In addition to known cachexia genes such as FOXO1, novel genes from muscle, including PPP1R8 and AEN correlated with cancer weight loss. All the adipose correlated genes including SCGN and EDR17 are novel for PDAC cachexia. Pathway analysis demonstrated shared pathways but largely non-overlapping genes in both tissues. Age related muscle loss predominantly had a distinct gene profiles compared to cachexia. This analysis of matched, externally validate gene expression points to novel targets in cachexia.

Manganese porphyrins reportedly exhibit synergic effects when combined with irradiation. However, an in-depth understanding of intratumoral heterogeneity and immune pathways, as affected by Mn porphyrins, remains limited. Here, we explored the mechanisms underlying immunomodulation of a clinical candidate, MnTnBuOE-2-PyP5+ (BMX-001, MnBuOE), using single-cell analysis in murine carcinoma model. Mice bearing 4T1 tumors were divided into 4 groups: control, MnBuOE, radiotherapy (RT), combined MnBuOE and radiotherapy (MnBuOE/RT). In epithelial cells, epithelial-mesenchymal transition, TNF-α signaling via NF-кB, angiogenesis, and hypoxia-related genes were significantly downregulated in the MnBuOE/RT compared to the RT. All subtypes of cancer-associated fibroblasts (CAFs) were clearly reduced in MnBuOE and MnBuOE/RT. Inhibitory receptor-ligand interactions, in which epithelial cells and CAFs interacted with CD8+ T cells, were significantly lower in the MnBuOE/RT than in the RT. Trajectory analysis showed that DC maturation-associated markers were increased in MnBuOE/RT. M1 macrophages were significantly increased in the MnBuOE/RT compared to the RT, whereas myeloid-derived suppressor cells were decreased. CellChat analysis showed that the number of cell-cell communications was the lowest in the MnBuOE/RT. Our study is the first to provide evidence for the combined radiotherapy with a novel Mn porphyrin clinical candidate, BMX-001 from the perspective of each cell-type within the tumor microenvironment.

(1) Background: Transcription factors (TFs) are main regulators of eukaryotic gene expression. The cooperative binding of at least two TFs to genomic DNA is the widespread mechanism of transcription regulation. Massive analysis of co-occurrence of overrepresented pairs of motifs (composite elements, CEs) for different target TFs studied in ChIP-seq experiments can clarify the mechanisms of TF cooperation. (2) Methods: We focused on homotypic CEs representing the same motif models and considered only CEs with spacers. We improved the capability of the Motifs Co-Occurrence Tool (MCOT) to predict asymmetric homotypic CEs in which one of par-ticipating motifs has higher conservation than the other. We categorized the target TFs from the data of M. musculus ChIP-seq and A. thaliana ChIP-seq/DAP-seq according to the structure of DNA-binding domains (DBDs) into classes. (3) Results: For all TF classes the predominance of asymmetric CEs over both symmetric and those with no enrichment in both directions was re-vealed. Target TFs from classes Basic helix-loop-helix factors (bHLH) and Basic leucine zipper factors (bZIP) showed the highest fractions of datasets with asymmetric CEs. We showed that pioneer TFs, despite their DBD types, have a higher significance of asymmetry within homotyp-ic CEs compared to other TFs. (4) Conclusion: Asymmetry within homotypic CEs is a promising new feature decrypting the mechanisms of gene transcription regulation.

The study of cell-free DNA (cfDNA) and other peripheral blood components (known as “liquid biopsies”) is promising and has been investigated especially in solid tumors. Nevertheless, it is increasingly showing greater utility in the diagnosis, prognosis, and response to treatment of hematological malignancies; in the future, it could prevent invasive techniques, such as bone marrow (BM) biopsy. Most of the studies about this topic have been focused on B cell lymphoid malignancies; some of them have shown that cfDNA can be used as a novel way for diagnosis and minimal residual monitoring in B cell lymphomas, using techniques such as next-generation sequencing (NGS). In myelodysplastic syndromes, multiple myeloma, or chronic lymphocytic leukemia, liquid biopsies may allow for an interesting genomic representation of the tumor clones affecting different lesions (spatial heterogeneity). In acute leukemias, it can be helpful in the monitoring of early treatment response and the prediction of treatment failure. In chronic lymphocytic leukemia, the evaluation of cfDNA permits the definition of clonal evolution and drug resistance in real-time. However, there are limitations such as the difficulty in obtaining sufficient circulating tumor DNA for achieving a high sensitivity to assess minimal residual disease or the lack of standardization of the method and clinical studies to confirm its prognostic impact. This review focuses on clinical applications of cfDNA on minimal residual disease in hematological malignancies.

Pathogenicity is unknown for the majority of human gene variants. For prioritization of sequenced somatic and germline mutation variants, in silico approaches can be utilized. In this study, 84 million non-synonymous Single Nucleotide Variants (SNVs) in the human coding genome were annotated using consensus Variant Effect Prediction (cVEP) method. An algorithm, implemented as a stacked ensemble of supervised learners, performed combination of the 39 functional, conservation mutation impact scores from dbNSFP4.0. Adding gene indispensability score, accounting for differences in the pathogenicities of the variants in the essential and the mutation-tolerant genes, improved the predictions. For each SNV the consensus combination gives either a continuous-value pathogenicity score, or a categorical score in five classes: pathogenic, likely pathogenic, uncertain significance, likely benign, benign. The provided class database is aimed for direct use in clinical practice. The trained prediction models were 5-fold cross-validated on the evidence-based categorical annotations from the ClinVar database. The rankings of the scores based on their ability to predict pathogenicity were obtained. A two-step strategy using the rankings, scores and class annotations is suggested for filtering and prioritization of the human exome mutations in clinical and biological applications of NGS technology.

Progesterone receptor (PGR) expression level determines biological characteristics of the serous ovarian carcinoma, and low PGR expression appears to be associated with chemoresistance and worse outcome. In this study, we aimed to find relationships between tumor progesterone receptor level and RNA-profile (miRNAs, piwiRNAs, and mRNAs), determining its biological characteristics and behavior. For this purpose, we applied next generation sequencing of small noncoding RNAs, quantitative RT-PCR, and immunohistochemistry to analyze FFPE and frozen tumor samples as well as blood plasma from patients with benign cystadenoma (BSC), serous borderline tumour (SBT), low-grade and high-grade serous ovarian carcinoma (LGSOC and HGSOC, respectively). We found significant upregulation of MMP7 and MUC16 and downregulation of PGR in LGSOC and HGSOC in comparison with BSC. The tissue content of miR-199a-5p, miR-214-3p, miR-424-3p, miR-424-5p, miR-125b-5p significantly inversely correlated with the expression level of MUC16 and blood serum CA125 concentration, and significantly directly correlated with the PGR expression level in tumor tissue. On the contrary, the tissue content of miR-16-5p, miR-17-5p, miR-20a-5p, miR-93-5p, responsible for epithelial-mesenchymal transition (EMT) of the cell, significantly directly correlated with the blood serum CA125 concentration and significantly inversely correlated with the PGR expression level in the tumor tissue. Levels of the EMT-associated miRNAs significantly directly correlated with the content of hsa_piR_022437, hsa_piR_009295, hsa_piR_020813, hsa_piR_004307, hsa_piR_019914 in tumor tissues. Among them, expression level of hsa_piR_004307 significantly inversly correlated with the PGR expression level in the tumor. Two optimal logistic regression models were developed based on the quantitation of hsa_piR_020813, miR-16-5p, hsa_piR_022437 or hsa_piR_004307, hsa_piR_019914, miR-93-5p in the tumor tissue, both of which significantly diagnose PGR-negative tumor phenotype with 93% sensitivity. According to FunRich3.1.3 functional enrichment analysis tool, 72 gene-targets of miRNA and piRNA identified here as markers of PGR-negative ovarian tumor phenotype were proven to be mutated in different cancers such as ovarian, breast, colorectal, liver, stomach, lung, endometrial, thyroid cancer. Among them, the blood plasma levels of miR-16-5p and hsa_piR_022437 can be used to diagnose PGR-negative tumor phenotype with 86% sensitivity before surgery and chemotherapy to choose the treatment strategy for this most aggressive type of ovarian cancer (for instance, neoadjuvant chemotherapy followed by cytoreduction in combination with hyperthermic intraperitoneal chemotherapy) to increase the effectiveness of treatment and longevity of the patient.

Background: Transcription factors (TFs) are main regulators of eukaryotic gene expression. The cooperative binding to genomic DNA of at least two TFs is the widespread mechanism of transcription regulation. Cooperating TFs can be revealed through the analysis of co-occurrence of their motifs. Methods: We applied Motifs Co-Occurrence Tool (MCOT) that predicted pairs of spaced or overlapped motifs (composite elements, CEs) for a single ChIP-seq dataset. We improved MCOT capability for prediction of asymmetric CEs with one of participating motifs possessing higher conservation than another does. Results: Analysis of 119 ChIP-seq datasets for 45 human TFs revealed that almost for all families of TFs the co-occurrence with an overlap between motifs of target TFs and more conserved partner motifs was significantly higher than that for less conserved partner motifs. The asymmetry toward partner TFs was the most clear for partner motifs of TFs from ETS family. Conclusion: Co-occurrence with an overlap of less conserved motif of a target TF and more conserved motifs of partner TFs explained a substantial portion of ChIP-seq data lacking conserved motifs of target TFs. Among other TF families, conservative motifs of TFs from ETS family were the most prone to mediate interaction of target TFs with its weak motifs in ChIP-seq.