The amino acid transporter SLC7A11 expression in breast cancer

Breast cancer (BC), characterized by its diverse molecular profiles and clinical outcomes, presents a significant challenge in the development of effective therapeutic strategies. Metabolic reprogramming, a defining characteristic of cancer, has emerged as a promising target for novel therapies. SLC7A11, an amino acid transporter that facilitates cysteine uptake in exchange for glutamate, plays a crucial role in sustaining the altered metabolism of cancer cells. This study delves into the comprehensive analysis of SLC7A11 at the genomic, transcriptomic, and protein levels in extensive BC datasets to elucidate its potential role in different BC subtypes. SLC7A11 gene copy number and mRNA expression were evaluated using the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) cohort (n = 1,980) and Breast Cancer Gene Expression Miner (n = 4,712). SLC7A11 protein was assessed using immunohistochemistry in a large BC cohort (n = 1,981). Additionally, The Cancer Genome Atlas (TCGA) dataset was used to explore SLC7A11 DNA methylation patterns using MethSurv (n = 782) and association of SLC7A11 mRNA expression with immune infiltrates using TIMER (n = 1,100). High SLC7A11 mRNA and SLC7A11 protein expression were significantly associated with high tumor grade (p ≤ .02), indicating a potential role in cancer progression. Interestingly, SLC7A11 copy number gain was observed in HER2+ tumors (p = .01), suggesting a subtype-specific association. In contrast, SLC7A11 mRNA expression was higher in the basal-like/triple-negative (TN; p < .001) and luminal B tumors (p = .02), highlighting its differential expression across BC subtypes. Notably, high SLC7A11 protein expression was predominantly observed in Estrogen Receptor (ER)-negative and Triple Negative (TN) BC, suggesting a role in these aggressive subtypes. Further analysis revealed that SLC7A11 was positively correlated with other amino acid transporters and enzymes associated with glutamine metabolism, implying a coordinated role in metabolic regulation. Additionally, SLC7A11 gene expression was positively associated with neutrophil and macrophage infiltration, suggesting a potential link between SLC7A11 and tumor immunity. Our findings suggest that SLC7A11 plays a significant role in BC metabolism, demonstrating differential expression across subtypes and associations with poor patient outcomes. Further functional studies are warranted to elucidate the precise mechanisms by which SLC7A11 contributes to BC progression and to explore its potential as a therapeutic target.

Genomic epidemiology and ceftazidime-avibactam high-level resistance mechanisms of Pseudomonas aeruginosa in China from 2010 to 2022

Ceftazidime-avibactam (CZA) resistance is a huge threat in the clinic; however, the underlying mechanism responsible for high-level CZA resistance in Pseudomonas aeruginosa (PA) isolates remains unknown. In this study, a total of 5,763 P. aeruginosa isolates were collected from 2010 to 2022 to investigate the ceftazidime-avibactam (CZA) high-level resistance mechanisms of Pseudomonas aeruginosa (PA) isolates in China. Fifty-six PER-producing isolates were identified, including 50 isolates carrying blaPER-1 in PA, and 6 isolates carrying blaPER-4. Of these, 82.1% (46/56) were classified as DTR-PA isolates, and 76.79% (43/56) were resistant to CZA. Importantly, blaPER-1 and blaPER-4 overexpression led to 16-fold and >1024-fold increases in the MICs of CZA, respectively. WGS revealed that the blaPER-1 gene was located in two different transferable IncP-2-type plasmids and chromosomes, whereas blaPER-4 was found only on chromosomes and was carried by a class 1 integron embedded in a Tn6485-like transposon. Overexpression of efflux pumps may be associated with high-level CZA resistance in blaPER-1-positive strains. Kinetic parameter analysis revealed that PER-4 exhibited a similar kcat/Km with ceftazidime and a high (∼3359-fold) IC50 value with avibactam compared to PER-1. Our study found that overexpression of PER-1 combined with enhanced efflux pump expression and the low affinity of PER-4 for avibactam contributes to high-level resistance to CZA. Additionally, the Tn6485-like transposon plays a significant role in disseminating blaPER. Urgent active surveillance is required to prevent the further spread of high-level CZA resistance in DTR-PA isolates.

Genomic and eco-epidemiological investigations in Uruguay reveal local Chikungunya virus transmission dynamics during its expansion across the Americas in 2023

Uruguay experienced its first Chikungunya virus outbreak in 2023, resulting in a significant burden to its healthcare system. We conducted analysis based on real-time genomic surveillance (30 novel whole genomes) to offer timely insights into recent local transmission dynamics and eco-epidemiological factors behind its emergence and spread in the country.

Genomic analysis of lineage-specific transmission of multidrug resistance tuberculosis in China

OBJECTIVES

We investigated the genetic diversities and lineage-specific transmission dynamics of multidrug-resistant tuberculosis (MDR-TB), with the goal of determining the potential factors driving the MDR epidemics in China.

METHODS

We curated a large nationwide Mycobacterium tuberculosis (M. tuberculosis) whole genome sequence data set, including 1313 MDR strains. We reconstructed the phylogeny and mapped the transmission networks of MDR-TB across China using Bayesian inference. To identify drug-resistance variants linked to enhanced transmissibility, we employed ordinary least-squares (OLS) regression analysis.

RESULT

The majority of MDR-TB strains in China belong to lineage 2.2.1. Transmission chain analysis has indicated that the repeated and frequent transmission of L2.2.1 plays a central role in the establishment of MDR epidemic in China, but no occurrence of a large predominant MDR outbreak was detected. Using OLS regression, the most common single nucleotide polymorphisms (SNPs) associated with resistance to isoniazid (katG_p.Ser315Thr and katG_p.Ser315Asn) and rifampicin (rpoB_p.Ser450Leu, rpoB_p.His445Tyr, rpoB_p.His445Arg, rpoB_p.His445Asp, and rpoB_p.His445Asn) were more likely to be found in L2 clustered strains. Several putative compensatory mutations in rpoA, rpoC, and katG were significantly associated with clustering. The eastern, central, and southern regions of China had a high level of connectivity for the migration of L2 MDR strains throughout the country. The skyline plot showed distinct population size expansion dynamics for MDR-TB lineages in China.

CONCLUSION

MDR-TB epidemic in China is predominantly driven by the spread of highly transmissible Beijing strains. A range of drug-resistance mutations of L2 MDR-TB strains displayed minimal fitness costs and may facilitate their transmission.

Insight into the complexity of male infertility: a multi-omics review

Male infertility is a reproductive disorder, accounting for 40-50% of infertility. Currently, in about 70% of infertile men, the cause remains unknown. With the introduction of novel omics and advancement in high-throughput technology, potential biomarkers are emerging. The main purpose of our work was to overview different aspects of omics approaches in association with idiopathic male infertility and highlight potential genes, transcripts, non-coding RNA, proteins, and metabolites worth further exploring. Using the Gene Ontology (GO) analysis, we aimed to compare enriched GO terms from each omics approach and determine their overlapping. A PubMed database screening for the literature published between February 2014 and June 2022 was performed using the keywords: male infertility in association with different omics approaches: genomics, epigenomics, transcriptomics, ncRNAomics, proteomics, and metabolomics. A GO enrichment analysis was performed using the Enrichr tool. We retrieved 281 global studies: 171 genomics (DNA level), 21 epigenomics (19 of methylation and two histone residue modifications), 15 transcriptomics, 31 non-coding RNA, 29 proteomics, two protein posttranslational modification, and 19 metabolomics studies. Gene ontology comparison showed that different omics approaches lead to the identification of different molecular factors and that the corresponding GO terms, obtained from different omics approaches, do not overlap to a larger extent. With the integration of novel omics levels into the research of idiopathic causes of male infertility, using multi-omic systems biology approaches, we will be closer to finding the potential biomarkers and consequently becoming aware of the entire spectrum of male infertility, their cause, prognosis, and potential treatment.

Within-host genomic evolution of methicillin-resistant Staphylococcus aureus in long-term carriers

Assessing the genomic evolution of Staphylococcus aureus can help us understand how the bacteria adapt to its environment. In this study, we aimed to assess the mutation rate within 144 methicillin-resistant Staphylococcus aureus (MRSA) carriers with a carriage time from 4 to 11 years, including some carriers who belonged to the same households. We found that 23 of the 144 individuals had completely different MRSA types over time and were therefore not long-term carriers of the same MRSA. From the remaining 121 individuals, we performed whole-genome sequencing (WGS) on 424 isolates and then compared these pairwise using core genome multilocus sequence typing (cgMLST) and single-nucleotide polymorphism (SNP) analyses. We found a median within-host mutation rate in long-term MRSA carriers of 4.9 (3.4-6.9) SNPs/genome/year and 2.7 (1.8-4.2) allelic differences/genome/year, when excluding presumed recombination. Furthermore, we stratified the cohort into subgroups and found no significant difference between the median mutation rate of members of households, individuals with presumed continued exposure, e.g., from travel and persons without known continued exposure. Finally, we found that SNPs occurred at random within the genes in our cohort. KEY POINTS: • Median mutation rate within long-term MRSA carriers of 4.9 (3.4-6.9) SNPs/genome/year • Similar median mutation rates in subgroups (households, travelers) • No hotspots for SNPs within the genome.

Transmission of fluoroquinolones resistance among multidrug-resistant tuberculosis in Shanghai, China: a retrospective population-based genomic epidemiology study

Fluoroquinolones (FQ) are essential for the treatment of multidrug-resistant tuberculosis (MDR-TB). The FQ resistance (FQ-R) rate in MDR-TB in China and its risk factors remain poorly understood. We conducted a retrospective, population-based genomic epidemiology study of MDR-TB patients in Shanghai, China, from 2009 to 2018. A genomic cluster was defined as strains with genetic distances ≤ 12 single nucleotide polymorphisms. The transmitted FQ-R was defined as the same FQ resistance-conferring mutations shared by ≥ 2 strains in a genomic cluster. We used multivariable logistic regression analysis to identify the risk factors for drug resistance. Among the total 850 MDR-TB patients included in the study, 72.8% (619/850) were male, the median age was 39 (interquartile range 28, 55) years, 52.7% (448/850) were migrants, and 34.5% (293/850) were previously treated patients. Most of the MDR-TB strains belong to the Beijing lineage (91.7%, 779/850). Overall, the genotypic resistance rate of FQ was 34.7% (295/850), and 47.1% (139/295) FQ-R patients were in genomic clusters, of which 98 (33.2%, 98/295) were presumed as transmitted FQ-R. Patients with treatment-naïve (aOR = 1.84; 95% CI: 1.09, 3.16), diagnosed in a district-level hospital (aOR = 2.69; 95% CI: 1.56, 4.75), and streptomycin resistance (aOR = 3.69; 95% CI: 1.65, 9.42) were significantly associated with the transmission of FQ-R. In summary, the prevalence of FQ-R among MDR-TB patients was high in Shanghai, and at least one-third were transmitted. Enforced interventions including surveillance of FQ drug susceptibility testing and screening among MDR-TB before initiation of treatment were urgently needed.

Fetal congenital talipes equinovarus: genomic abnormalities and obstetric follow-up results

OBJECTIVE

The etiology of congenital talipes equinovarus (CTEV) is unknown, and the relationship between chromosome microdeletion/microduplication and fetal CTEV is rarely reported. In this study, we retrospectively analyzed fetal CTEV to explore the relationship among the CTEV phenotype, chromosome microdeletion/microduplication, and obstetric outcomes.

METHODS

Chromosome karyotype analysis and single nucleotide polymorphism (SNP) array were performed for the 68 fetuses with CTEV.

RESULTS

An SNP array was performed for 68 fetuses with CTEV; pathogenic copy number variations (CNVs) were detected in eight cases (11.8%, 8/68). In addition to one case consistent with karyotype analysis, the SNP array revealed seven additional pathogenic CNVs, including three with 22q11.21 microdeletions, two with 17p12p11.2 microduplications, one with 15q11.2 microdeletions, and one with 7q11.23 microduplications. Of the seven cases carrying pathogenic CNVs, three were tested for family genetics; of these, one was de novo, and two were inherited from either the father or mother. In total, 68 fetuses with CTEV were initially identified, of which 66 cases successfully followed-up. Of these, 9 were terminated, 2 died in utero, and 55 were live births. In 9 cases, no clinical manifestations of CTEV were found at birth; the false-positive rate of prenatal ultrasound CTEVdiagnosis was thus 13.6% (9/66).

CONCLUSION

CTEV was associated with chromosome microdeletion/microduplication, the most common of which was 22q11.21 microdeletion, followed by 17p12p11.2 microduplication. Thus, further genomic detection is recommended for fetuses with CTEV showing no abnormalities on conventional karyotype analysis.

RNA binding motif 4 inhibits the replication of ebolavirus by directly targeting 3'-leader region of genomic RNA

Ebola virus (EBOV) belongs to Filoviridae family possessing single-stranded negative-sense RNA genome, which is a serious threat to human health. Nowadays, no therapeutics have been proven to be successful in efficiently decreasing the mortality rate. RNA binding proteins (RBPs) are reported to participate in maintaining cell integrity and regulation of viral replication. However, little is known about whether and how RBPs participate in regulating the life cycle of EBOV. In our study, we found that RNA binding motif protein 4 (RBM4) inhibited the replication of EBOV in HEK293T and Huh-7 cells by suppressing viral mRNA production. Such inhibition resulted from the direct interaction between the RRM1 domain of RBM4 and the "CU" enrichment elements located in the PE1 and TSS of the 3'-leader region within the viral genome. Simultaneously, RBM4 could upregulate the expression of some cytokines involved in the host innate immune responses to synergistically exert its antiviral function. The findings therefore suggest that RBM4 might serve as a novel target of anti-EBOV strategy.

Tracing histoplasmosis genomic epidemiology and species occurrence across the USA

ABSTRACTHistoplasmosis is an endemic mycosis in North America frequently reported along the Ohio and Mississippi River Valleys, although autochthonous cases occur in non-endemic areas. In the United States, the disease is provoked by two genetically distinct clades of Histoplasma capsulatum sensu lato, Histoplasma mississippiense (Nam1) and H. ohiense (Nam2). To bridge the molecular epidemiological gap, we genotyped 93 Histoplasma isolates (62 novel genomes) including clinical, environmental, and veterinarian samples from a broader geographical range by whole-genome sequencing, followed by evolutionary and species niche modelling analyses. We show that histoplasmosis is caused by two major lineages, H. ohiense and H. mississippiense; with sporadic cases caused by H. suramericanum in California and Texas. While H. ohiense is prevalent in eastern states, H. mississipiense was found to be prevalent in the central and western portions of the United States, but also geographically overlapping in some areas suggesting that these species might co-occur. Species Niche Modelling revealed that H. ohiense thrives in places with warmer and drier conditions, while H. mississippiense is endemic to areas with cooler temperatures and more precipitation. In addition, we predicted multiple areas of secondary contact zones where the two species co-occur, potentially facilitating gene exchange and hybridization. This study provides the most comprehensive understanding of the genomic epidemiology of histoplasmosis in the USA and lays a blueprint for the study of invasive fungal diseases.

Investigation of genomic and pathogenicity characteristics of Streptococcus suis ST1 human strains from Guangxi Zhuang Autonomous Region (GX) between 2005 and 2020 in China

Streptococcus suis is a significant and emerging zoonotic pathogen. ST1 and ST7 strains are the primary agents responsible for S. suis human infections in China, including the Guangxi Zhuang Autonomous Region (GX). To enhance our understanding of S. suis ST1 population characteristics, we conducted an investigation into the phylogenetic structure, genomic features, and virulence levels of 73 S. suis ST1 human strains from GX between 2005 and 2020. The ST1 GX strains were categorized into three lineages in phylogenetic analysis. Sub-lineage 3-1a exhibited a closer phylogenetic relationship with the ST7 epidemic strain SC84. The strains from lineage 3 predominantly harboured 89K-like pathogenicity islands (PAIs) which were categorized into four clades based on sequence alignment. The acquirement of 89K-like PAIs increased the antibiotic resistance and pathogenicity of corresponding transconjugants. We observed significant diversity in virulence levels among the 37 representative ST1 GX strains, that were classified as follows: epidemic (E)/highly virulent (HV) (32.4%, 12/37), virulent plus (V+) (29.7%, 11/37), virulent (V) (18.9%, 7/37), and lowly virulent (LV) (18.9%, 7/37) strains based on survival curves and mortality rates at different time points in C57BL/6 mice following infection. The E/HV strains were characterized by the overproduction of tumour necrosis factor (TNF)-α in serum and promptly established infection at the early phase of infection. Our research offers novel insights into the population structure, evolution, genomic features, and pathogenicity of ST1 strains. Our data also indicates the importance of establishing a scheme for characterizing and subtyping the virulence levels of S. suis strains.

Genetic characteristics and potential pathogenic agents in Campylobacter upsaliensis based on genomic analysis

Campylobacter upsaliensis was the most common Campylobacter species in pets' gastrointestinal tracts and has been isolated from patients with bacteremia, hemolytic-uremic syndrome, spontaneous abortion, and Guillain-Barré syndrome. However, the genetic characteristics and the full extent of its significance as a human pathogen remain to be fully understood. This study involved an investigation for genomic analysis of 154 strains from different sources and additional antimicrobial resistance profiles of 26 strains for this species. The genomes contained 1,558-1,971 CDS and the genome sizes were estimated to vary from 1.53 Mb to 1.86 Mb, with an average GC content of  34.71%. The entire analyzed genomes could be divided into three clades (A, B, and C) based on ANI and phylogenomic analysis. Significantly, nearly all strains in Clade B were isolated from patient samples, and the virulence-related sequences FlgD, GmhA, and CdtC might serve as determining factors for the classification of Clade B. Half of the tested isolates had MIC values over 64 μg mL-1 for nalidixic acid, gentamicin, and streptomycin. Isolates from pets in China carried more resistant elements in the genomes. This study both provided a comprehensive profile of C. upsaliensis for its genomic features and suggested some pathogenic agents for human infection with this species.

Genomic characters of Anaplasma bovis and genetic diversity in China

The emergence of Anaplasma bovis or A. bovis-like infection in humans from China and the United States of America has raised concern about the public health importance of this pathogen. Although A. bovis has been detected in a wide range of ticks and mammals in the world, no genome of the pathogen is available up to now, which has prohibited us from better understanding the genetic basis for its pathogenicity. Here we describe an A. bovis genome from metagenomic sequencing of an infected goat in China. Anaplasma bovis had the smallest genome of the genus Anaplasma, and relatively lower GC content. Phylogenetic analysis of single-copy orthologue sequence showed that A. bovis was closely related to A. platys and A. phagocytophilum, but relatively far from intraerythrocytic Anaplasma species. Anaplasma bovis had 116 unique orthogroups and lacked 51 orthogroups in comparison to other Anaplasma species. The virulence factors of A. bovis were significantly less than those of A. phagocytophilum, suggesting less pathogenicity of A. bovis. When tested by specific PCR assays, A. bovis was detected in 23 of 29 goats, with an infection rate up to 79.3% (95% CI: 64.6% ∼94.1%). The phylogenetic analyses based on partial 16S rRNA, gltA and groEL genes indicated that A. bovis had high genetic diversity. The findings of this study lay a foundation for further understanding of the biological characteristics and genetic diversity of A. bovis, and will facilitate the formulation of prevention and control strategies.

Exome-wide comparative analyses revealed differentiating genomic regions for performance traits in Indian native buffaloes

In India, 20 breeds of buffalo have been identified and registered, yet limited studies have been conducted to explore the performance potential of these breeds, especially in the Indian native breeds. This study is a maiden attempt to delineate the important variants and unique genes through exome sequencing for milk yield, milk composition, fertility, and adaptation traits in Indian local breeds of buffalo. In the present study, whole exome sequencing was performed on Chhattisgarhi (n = 3), Chilika (n = 4), Gojri (n = 3), and Murrah (n = 4) buffalo breeds and after stringent quality control, 4333, 6829, 4130, and 4854 InDels were revealed, respectively. Exome-wide FST along 100-kb sliding windows detected 27, 98, 38, and 35 outlier windows in Chhattisgarhi, Chilika, Gojri, and Murrah, respectively. The comparative exome analysis of InDels and subsequent gene ontology revealed unique breed specific genes for milk yield (CAMSAP3), milk composition (CLCN1, NUDT3), fertility (PTGER3) and adaptation (KCNA3, TH) traits. Study provides insight into mechanism of how these breeds have evolved under natural selection, the impact of these events on their respective genomes, and their importance in maintaining purity of these breeds for the traits under study. Additionally, this result will underwrite to the genetic acquaintance of these breeds for breeding application, and in understanding of evolution of these Indian local breeds.

Weighted single step GWAS reveals genomic regions associated with economic traits in Murrah buffaloes

The objective of the present study was to identify genomic regions influencing economic traits in Murrah buffaloes using weighted single step Genome Wide Association Analysis (WssGWAS). Data on 2000 animals, out of which 120 were genotyped using a double digest Restriction site Associated DNA (ddRAD) sequencing approach. The phenotypic data were collected from NDRI, India, on growth traits, viz., body weight at 6M (month), 12M, 18M and 24M, production traits like 305D (day) milk yield, lactation length (LL) and dry period (DP) and reproduction traits like age at first calving (AFC), calving interval (CI) and first service period (FSP). The biallelic genotypic data consisted of 49353 markers post-quality check. The heritability estimates were moderate to high, low to moderate, low for growth, production, reproduction traits, respectively. Important genomic regions explaining more than 0.5% of the total additive genetic variance explained by 30 adjacent SNPs were selected for further analysis of candidate genes. In this study, 105 genomic regions were associated with growth, 35 genomic regions with production and 42 window regions with reproduction traits. Different candidate genes were identified in these genomic regions, of which important are OSBPL8, NAP1L1 for growth, CNTNAP2 for production and ILDR2, TADA1 and POGK for reproduction traits.

Whole genome resequencing revealed genomic variants and functional pathways related to adaptation in Indian yak populations

The present study aims to identify genomic variants through a whole genome sequencing (WGS) approach and uncover biological pathways associated with adaptation and fitness in Indian yak populations. A total of 30 samples (10 from each population) were included from Arunachali, Himachali and Ladakhi yak populations. WGS analysis revealed a total of 32171644, 27260825, and 32632460 SNPs and 4865254, 4429941, and 4847513 Indels in the Arunachali, Himachali, and Ladakhi yaks, respectively. Genes such as RYR2, SYNE2, BOLA, HF1, and the novel transcript ENSBGRG00000011079 were found to have the maximum number of high impact variants in all three yak populations, and might play a major role in local adaptation. Functional enrichment analysis of genes harboring high impact SNPs revealed overrepresented pathways related to response to stress, immune system regulation, and high-altitude adaptation. This study provides comprehensive information about genomic variants and their annotation in Indian yak populations, thus would serve as a data resource for researchers working on the yaks. Furthermore, it could be well exploited for better yak conservation strategies by estimating population genetics parameters viz., effective population size, inbreeding, and observed and expected heterozygosity.

Genomic characterization, phylogenetic and expression analysis of foraging gene in Apis mellifera

The genomic characterization of the foraging gene and its expression analysis are required to better understand the behavior of honey bees (Apis mellifera). The present study performed a genome-wide characterization of the foraging gene, analyzing its physicochemical properties, phylogenetic features, and expression. An in silico analysis was carried out to characterize the foraging gene and the motifs and conserved domains of the encoded protein to predict its physicochemical properties. Moreover, a phylogenetic analysis of the foraging gene was performed in different species using MEGAX. The relative expression of the foraging gene was determined using qRT-PCR in two groups of forager bee samples (incoming and outgoing bees) during two seasons (five times per day). In addition, the queen effect was evaluated in another experiment. The results revealed that foraging gene expression and bee traffic were influenced by the interaction of season and daytime. The daily foraging traffic and transcription level of the foraging gene were the same in both seasons. The traffic of bees and the transcription abundance of the foraging gene were the highest in the middle and at the end of the day in the first and second seasons, respectively. Furthermore, the mRNA expression of the foraging gene was relatively higher in incoming bees than in outgoing bees. The queen also had a significant effect on the outgoing bees. We conclude that gene-environment interactions affect the foraging behavior of bees through the modulation of the foraging gene transcription.

Bad to the bone? - Genomic analysis of Enterococcus isolates from diverse environments reveals that most are safe and display potential as food fermentation microorganisms

Enterococci comprise a group of lactic acid bacteria (LAB) with considerable potential to serve as food fermentation microorganisms. Unfortunately, enterococci have received a lot of negative attention, due to the occurrence of pathogenic and multidrug resistant strains. In this study, we used genomics to select safe candidates among the forty-four studied enterococcal isolates. The genomes of the forty-four strains were fully sequenced and assessed for presence of virulence and antibiotic resistance genes. Nineteen isolates belonging to the species Enterococcus lactis, Enterococcus faecium, Enterococcus durans, and Enterococcus thailandicus, were deemed safe from the genome analysis. The presence of secondary metabolite gene clusters for bacteriocins was assessed, and twelve candidates were found to secrete antimicrobial compounds effective against Listeria monocytogenes isolated from cheese and Staphylococcus aureus. Physiological characterization revealed nineteen industrial potentials; all strains grew well at 42 °C and acidified 1.5 hours faster than their mesophilic counterpart Lactococcus lactis, with which they share metabolism and flavor forming ability. We conclude that a large fraction of the examined enterococci were safe and could serve as excellent food fermentation microorganisms with inherent bioprotective abilities.

Diverse genetic causes of amenorrhea in an ethnically homogeneous cohort and an evolving approach to diagnosis

RESEARCH QUESTION

Premature ovarian insufficiency (POI) is characterised by amenorrhea associated with elevated follicle stimulating hormone (FSH) under the age of 40 years and affects 1-3.7% women. Genetic factors explain 20-30% of POI cases, but most causes remain unknown despite genomic advancements.

DESIGN

We used whole exome sequencing (WES) in four Iranian families, validated variants via Sanger sequencing, and conducted the Acyl-cLIP assay to measure HHAT enzyme activity.

RESULTS

Despite ethnic homogeneity, WES revealed diverse genetic causes, including a novel homozygous nonsense variant in SYCP2L, impacting synaptonemal complex (SC) assembly, in the first family. Interestingly, the second family had two independent causes for amenorrhea - the mother had POI due to a novel homozygous loss-of-function variant in FANCM (required for chromosomal stability) and her daughter had primary amenorrhea due to a novel homozygous GNRHR (required for gonadotropic signalling) frameshift variant. WES analysis also provided cytogenetic insights. WES revealed one individual was in fact 46, XY and had a novel homozygous missense variant of uncertain significance in HHAT, potentially responsible for complete sex reversal although functional assays did not support impaired HHAT activity. In the remaining individual, WES indicated likely mosaic Turners with the majority of X chromosome variants having an allelic balance of ∼85% or ∼15%. Microarray validated the individual had 90% 45,XO.

CONCLUSIONS

This study demonstrates the diverse causes of amenorrhea in a small, isolated ethnic cohort highlighting how a genetic cause in one individual may not clarify familial cases. We propose that, in time, genomic sequencing may become a single universal test required for the diagnosis of infertility conditions such as POI.

Comprehensive genomic and plasmid characterization of multidrug-resistant bacterial strains by R10.4.1 nanopore sequencing

The escalating prevalence of multidrug-resistant (MDR) bacteria pose a significant public health threat. Understanding the genomic features and deciphering the antibiotic resistance profiles of these pathogens is crucial for development of effective surveillance and treatment strategies. In this study, we employed the R10.4.1 nanopore sequencing technology, specifically through the use of the MinION platform, to analyze eight MDR bacterial strains originating from clinical, ecological and food sources. A single 72-hour sequencing run could yield approximately 12 million reads which covered a total of 34 gigabases (Gbp). The nanopore R10.4.1 data was processed using the Flye assembler, successfully assembling the genomes of eight bacterial strains and their 18 plasmids. Notably, the assemblies generated solely from R10.4.1 nanopore data closely matched those from next-generation sequencing data. Diverse antibiotic resistance patterns and specific resistance genes in the test strains were identified. Hospital strains that exhibited multidrug resistance were found to harbor various resistance genes that encode efflux pumps and extended-spectrum β-lactamases. Environmental and food sources were found to display resistance profiles in a species-specific manner. The composition of structurally complex plasmids in the test strains could also be revealed by analysis of nanopore long reads, which also suggested evidence of horizontal transfer of plasmids between different bacterial species. These findings provide valuable insights into the genetic characteristics of MDR bacteria and demonstrating the practicality of nanopore sequencing technology for detecting of resistance elements in bacterial pathogens.

PHDtools: A platform for pathogen detection and multi-dimensional genetic signatures decoding to realize pathogen genomics data analyses online

OBJECTIVES

Facing the emerging diseases, rapid identification of the pathogen and multi-dimensional characterization of the genomic features at both isolate-level and population-level through high-throughput sequencing data can provide invaluable information to guide the development of antiviral agents and strategies. However, a user-friendly program is in urgent need for clinical laboratories without bioinformatics background to decode the complex big genomics data.

METHODS

In this study, we developed an interactive online platform named PHDtools with a total of 15 functions to analyze metagenomics data to identify the potential pathogen and decode multi-dimensional genetic signatures including intra-/inter-host variations and lineage-level variations. The platform was applied to analyze the meta-genomic data of the samples collected from the 172 imported COVID-19 cases.

RESULTS

According to the analytical results of mNGS, 27 patients were found to have the co-infections of SARS-CoV-2 with either influenza virus (n = 9) or human picobirnavirus (n = 19). Enough coverages of all the assembled SARS-CoV-2 genomes provided the sub-lineages of Omicron variant, and the number of mutations in the non-structural genes and M gene was increased, as well as the intra-host variations occurred in E and M gene were under positive selection (Ka/Ks > 1). These findings of increased or changed mutations in the SARS-CoV-2 genome characterized the current adaptive evolution patterns of Omicron sub-lineages, and revealed the evolution speed of these sub-lineages might increase.

CONCLUSIONS

Consequently, the application of PHDtools has proved that this platform is accurate, user-friendly and convenient for clinical users who are deficient in bioinformatics, and the clinical monitor of SARS-CoV-2 genomes by PHDtools also highlighted the potential evolution features of current SARS-CoV-2 and indicated that the development of anti-SARS-CoV-2 agents and new-designed vaccines should incorporate the gene variations other than S gene.

Genomic insights into glume pubescence in durum wheat: GWAS and haplotype analysis implicates TdELD1-1A as a candidate gene

Glume pubescence is an important morphological trait for the characterization of wheat cultivars. It shows tolerance to biotic and abiotic stresses to some extent. Hg1 (formerly named Hg) locus on chromosome 1AS controls glume pubescence in wheat. Its genetic analysis, fine-mapping and candidate gene analysis have been widely studied recently, however, the cloning of Hg1 has not yet been reported. Here, we conducted a GWAS between a dense panel of 171,103 SNPs and glume pubescence (Gp) in a durum wheat population of 145 lines, and further analyzed the candidate genes of Hg1 combined with the gene expression, functional annotation, and haplotype analysis. As a results, TRITD0Uv1G104670 (TdELD1-1A), encoding glycosyltransferase-like ELD1/KOBITO 1, was detected as the most promising candidate gene of Hg1 for glume pubescence in durum wheat. Our findings not only contribute to a deeper understanding of its cloning and functional validation but also underscore the significance of accurate genome sequences and annotations. Additionally, our study highlights the relevance of unanchored sequences in chrUn and the application of bioinformatics analysis for gene discovery in durum wheat.

Comparative genomics and phylogenetic analysis of chloroplast genomes of Asian Caryodaphnopsis taxa (Lauraceae)

The genus Caryodaphnopsis, a member of the Lauraceae family, is characterized by seeds that are rich in oil, as well as highly exploitable fruits and wood. The Asian taxa within this genus exhibit complex morphological variations, posing challenges to their accurate classification and impeding their effective use and development as a resource. In this study, we sequenced the chloroplast genomes of 31 individuals representing nine Asian taxa within the Caryodaphnopsis genus. Our primary objectives were to reveal structural variations in these chloroplast genomes through comparative analyses and to infer the species' phylogenetic relationships. Our findings revealed that all chloroplast genomes had a tetrad structure, ranged in length from 148,828 to 154,946 bp, and harbored 128-131 genes. Notably, contraction of the IR region led to the absence of some genes in eight taxa. A comprehensive analysis identified 1267 long repetitive sequences and 2176 SSRs, 286 SNPs, and 135 indels across the 31 chloroplast genomes. The Ka/Ks ratio analysis indicated potential positive selection on the matK, rpl22, and rpoC2 genes. Furthermore, we identified six variable regions as promising barcode regions. Phylogenetic analysis grouped the nine Asian taxa into six branches, with C. henryi forming the basal group from which three distinct complexes emerged. This study contributes significantly to the current understanding of the evolutionary dynamics and phylogenetic relationships within the genus Caryodaphnopsis. Furthermore, the identified molecular markers hold potential for molecular barcoding applications in population genetics, providing valuable tools for future research and conservation efforts within this diverse genus.

Comparative genomics analysis and transposon mutagenesis provides new insights into high menaquinone-7 biosynthetic potential of Bacillus subtilis natto

This research combined Whole-Genome sequencing, intraspecific comparative genomics and transposon mutagenesis to investigate the menaquinone-7 (MK-7) synthesis potential in Bacillus subtilis natto. First, Whole-Genome sequencing showed that Bacillus subtilis natto BN-P15-11-1 contains one single circular chromosome in size of 3,982,436 bp with a GC content of 43.85 %, harboring 4,053 predicted coding genes. Next, the comparative genomics analysis among strain BN-P15-11-1 with model Bacillus subtilis 168 and four typical Bacillus subtilis natto strains proves that the closer evolutionary relationship Bacillus subtilis natto BN-P15-11-1 and Bacillus subtilis 168 both exhibit strong biosynthetic potential. To further dig for MK-7 biosynthesis latent capacity of BN-P15-11-1, we constructed a mutant library using transposons and a high throughput screening method using microplates. We obtained a YqgQ deficient high MK-7 yield strain F4 with a yield 3.02 times that of the parent strain. Experiments also showed that the high yield mutants had defects in different transcription and translation regulatory factor genes, indicating that regulatory factor defects may affect the biosynthesis and accumulation of MK-7 by altering the overall metabolic level. The findings of this study will provide more novel insights on the precise identification and rational utilization of the Bacillus subtilis subspecies for biosynthesis latent capacity.

Amplification-free quantitative detection of genomic DNA using lateral flow strips for milk authentication

With the globalization and complexity of the food supply chain, the market is becoming increasingly competitive and food fraudulent activities are intensifying. The current state of food detection faced two primary challenges. Firstly, existing testing methods were predominantly laboratory-based, requiring complex procedures and precision instruments. Secondly, there was a lack of accurate and efficient quantitative detection methods. Taking cow's milk as an example, this study introduced a novel method for nucleic acid quantification in dairy products, based on lateral flow strips (LFS). The core idea of this method is to design single-stranded DNA (ssDNA) probes to hybridize with mitochondrial genes, which are abundant, stable, and species-specific in dairy products, as detection targets. Drawing inspiration from the principles of nucleic acid amplification, this research innovatively established a new DNA hybridization method, named LAMP-Like Hybridization (HybLAMP-Like). Leveraging the denaturation and DNA polymerization functions of the bst enzyme, efficient binding of the probe and template strand was achieved. This method eliminated the need for nucleic acid amplification, simplifying the procedure and mitigating aerosol contamination, thereby ensuring the accuracy of the detection results. The method exhibited exceptional sensitivity, capable of detecting extremely low to 12.5 ng in visual inspection and 3.125 ng when using a reader. In terms of practicality, it could achieve visual detection of cow's milk content as low as 1% in adulterated dairy products. When combined with a portable LFS reader, it also enabled precise quantitative analysis of milk adulteration.

Genomic analysis of feed efficiency traits in beef cattle using random regression models

Feed efficiency plays a major role in the overall profitability and sustainability of the beef cattle industry, as it is directly related to the reduction of the animal demand for input and methane emissions. Traditionally, the average daily feed intake and weight gain are used to calculate feed efficiency traits. However, feed efficiency traits can be analysed longitudinally using random regression models (RRMs), which allow fitting random genetic and environmental effects over time by considering the covariance pattern between the daily records. Therefore, the objectives of this study were to: (1) propose genomic evaluations for dry matter intake (DMI), body weight gain (BWG), residual feed intake (RFI) and residual weight gain (RWG) data collected during an 84-day feedlot test period via RRMs; (2) compare the goodness-of-fit of RRM using Legendre polynomials (LP) and B-spline functions; (3) evaluate the genetic parameters behaviour for feed efficiency traits and their implication for new selection strategies. The datasets were provided by the EMBRAPA-GENEPLUS beef cattle breeding program and included 2920 records for DMI, 2696 records for BWG and 4675 genotyped animals. Genetic parameters and genomic breeding values (GEBVs) were estimated by RRMs under ssGBLUP for Nellore cattle using orthogonal LPs and B-spline. Models were compared based on the deviance information criterion (DIC). The ranking of the average GEBV of each test week and the overall GEBV average were compared by the percentage of individuals in common and the Spearman correlation coefficient (top 1%, 5%, 10% and 100%). The highest goodness-of-fit was obtained with linear B-Spline function considering heterogeneous residual variance. The heritability estimates across the test period for DMI, BWG, RFI and RWG ranged from 0.06 to 0.21, 0.11 to 0.30, 0.03 to 0.26 and 0.07 to 0.27, respectively. DMI and RFI presented within-trait genetic correlations ranging from low to high magnitude across different performance test-day. In contrast, BWG and RWG presented negative genetic correlations between the first 3 weeks and the other days of performance tests. DMI and RFI presented a high-ranking similarity between the GEBV average of week eight and the overall GEBV average, with Spearman correlations and percentages of individuals selected in common ranging from 0.95 to 1.00 and 93 to 100, respectively. Week 11 presented the highest Spearman correlations (ranging from 0.94 to 0.98) and percentages of individuals selected in common (ranging from 85 to 94) of BWG and RWG with the average GEBV of the entire period of the test. In conclusion, the RRM using linear B-splines is a feasible alternative for the genomic evaluation of feed efficiency. Heritability estimates of DMI, RFI, BWG and RWG indicate enough additive genetic variance to achieve a moderate response to selection. A new selection strategy can be adopted by reducing the performance test to 56 days for DMI and RFI selection and 77 days for BWG and RWG selection.

Comparative pangenome analysis of Aspergillus flavus and Aspergillus oryzae reveals their phylogenetic, genomic, and metabolic homogeneity

Aspergillus flavus and Aspergillus oryzae are closely related fungal species with contrasting roles in food safety and fermentation. To comprehensively investigate their phylogenetic, genomic, and metabolic characteristics, we conducted an extensive comparative pangenome analysis using complete, dereplicated genome sets for both species. Phylogenetic analyses, employing both the entirety of the identified single-copy orthologous genes and six housekeeping genes commonly used for fungal classification, did not reveal clear differentiation between A. flavus and A. oryzae genomes. Upon analyzing the aflatoxin biosynthesis gene clusters within the genomes, we observed that non-aflatoxin-producing strains were dispersed throughout the phylogenetic tree, encompassing both A. flavus and A. oryzae strains. This suggests that aflatoxin production is not a distinguishing trait between the two species. Furthermore, A. oryzae and A. flavus strains displayed remarkably similar genomic attributes, including genome sizes, gene contents, and G + C contents, as well as metabolic features and pathways. The profiles of CAZyme genes and secondary metabolite biosynthesis gene clusters within the genomes of both species further highlight their similarity. Collectively, these findings challenge the conventional differentiation of A. flavus and A. oryzae as distinct species and highlight their phylogenetic, genomic, and metabolic homogeneity, potentially indicating that they may indeed belong to the same species.

Clinicopathologic features, concurrent genomic alterations, and clinical outcomes of patients with KRAS G12D mutations in resected lung adenocarcinoma

BACKGROUND

In light of the ongoing clinical development of KRAS G12D-specific inhibitors, we sought to investigate the clinicopathologic, co-occurring genomic features and outcomes of patients with KRAS G12D-mutant lung adenocarcinoma.

METHODS

3828 patients with completely resected primary lung adenocarcinomas were examined for KRAS mutations between 2008 and 2020. The association between KRAS G12D and clinicopathologic features, molecular profiles, and outcomes was investigated.

RESULTS

65 patients (1.7%) with KRAS G12D-mutant lung adenocarcinoma were identified. KRAS G12D mutation was more frequent in males, former/current smokers, radiologic solid tumors, and invasive mucinous adenocarcinoma. TP53 and STK11 were the two most frequent concomitant mutations in the KRAS G12D group. KRAS G12D mutation did not appear to be a prognostic factor in resected stage I-III lung adenocarcinomas, while KRAS non-G12D mutation was related to worse survival, especially in stage I tumors. KRAS G12D mutations were associated with positive but low (1-49%) PD-L1 expression compared to negative (<1%), while KRAS non-G12D mutation was associated with high PD-L1 expression (≥50%). TP53 co-mutation indicated higher PD-L1 expression, while STK11 co-mutation had a negligible impact on PD-L1 expression. Furthermore, data mining of MSK datasets from cBioPortal revealed that KRAS G12D and SKT11 co-mutation were associated with a diminished response to immunotherapy.

CONCLUSIONS

KRAS G12D-mutant lung adenocarcinoma harbored unique clinicopathologic and genomic characteristics. Despite not being prognostic in resected lung adenocarcinoma, KRAS G12D might be a valuable biomarker in combination with certain co-mutations for identifying relevant subgroups of patients that could eventually influence treatment regimens.

Validation of dried blood spots for capturing hepatitis C virus diversity for genomic surveillance

Dried blood spots (DBS) have emerged as a promising alternative to traditional venous blood for hepatitis C virus (HCV) testing. However, their capacity to accurately reflect the genetic diversity of HCV remains poorly understood. We employed deep sequencing and advanced phylogenetic analyses on paired plasma and DBS samples from two common subtypes to evaluate the suitability of DBS for genomic surveillance. Results demonstrated that DBS captured equivalent viral diversity compared to plasma with no phylogenetic discordance observed. The ability of DBS to accurately reflect the profile of viral genetic diversity suggests it may be a promising avenue for future surveillance efforts to curb HCV outbreaks.

Comprehensive landscape of non-CODIS STRs in global populations provides new insights into challenging DNA profiles

The worldwide implementation of short tandem repeats (STR) profiles in forensic genetics necessitated establishing and expanding the CODIS core loci set to facilitated efficient data management and exchange. Currently, the mainstay CODIS STRs are adopted in most general-purpose forensic kits. However, relying solely on these loci failed to yield satisfactory results for challenging tasks, such as bio-geographical ancestry inference, complex DNA mixture profile interpretation, and distant kinship analysis. In this context, non-CODIS STRs are potent supplements to enhance the systematic discriminating power, particularly when combined with the high-throughput next-generation sequencing (NGS) technique. Nevertheless, comprehensive evaluation on non-CODIS STRs in diverse populations was scarce, hindering their further application in routine caseworks. To address this gap, we investigated genetic variations of 178 historically available non-CODIS STRs from ethnolinguistically different worldwide populations and studied their characteristics and forensic potentials via high-coverage whole genome sequencing (WGS) data. Initially, we delineated the genomic properties of these non-CODIS markers through sequence searching, repeat structure scanning, and manual inspection. Subsequent population genetics analysis suggested that these non-CODIS STRs had comparable polymorphism levels and forensic utility to CODIS STRs. Furthermore, we constructed a theoretical next-generation sequencing (NGS) panel comprising 108 STRs (20 CODIS STRs and 88 non-CODIS STRs), and evaluated its performance in inferring bio-geographical ancestry origins, deconvoluting complex DNA mixtures, and differentiating distant kinships using real and simulated datasets. Our findings demonstrated that incorporating supplementary non-CODIS STRs enabled the extrapolation of multidimensional information from a single STR profile, thereby facilitating the analysis of challenging forensic tasks. In conclusion, this study presents an extensive genomic landscape of forensic non-CODIS STRs among global populations, and emphasized the imperative inclusion of additional polymorphic non-CODIS STRs in future NGS-based forensic systems.

Comprehensive Genomic Analysis Identifies a Diverse Landscape of Sideroblastic and Nonsideroblastic Iron-Related Anemias with Novel and Pathogenic Variants in an Iron-Deficient Endemic Setting

Inherited iron metabolism defects are possibly missed or underdiagnosed in iron-deficient endemic settings because of a lack of awareness or a methodical screening approach. Hence, we systematically evaluated anemia cases (2019 to 2021) based on clinical phenotype, normal screening tests (high-performance liquid chromatography, α gene sequencing, erythrocyte sedimentation rate, C-reactive protein, and tissue transglutaminase), and abnormal iron profile by targeted next-generation sequencing (26-gene panel) supplemented with whole-exome sequencing, multiplex ligation probe amplification/mitochondrial DNA sequencing, and chromosomal microarray. Novel variants in ALAS2, STEAP3, and HSPA9 genes were functionally validated. A total of 290 anemia cases were screened, and 41 (14%) enrolled for genomic testing as per inclusion criteria. Comprehensive genomic testing revealed pathogenic variants in 23 of 41 cases (56%). Congenital sideroblastic anemia was the most common diagnosis (14/23; 61%), with pathogenic variations in ALAS2 (n = 6), SLC25A38 (n = 3), HSPA9 (n = 2) and HSCB, SLC19A2, and mitochondrial DNA deletion (n = 1 each). Nonsideroblastic iron defects included STEAP3-related microcytic anemia (2/23; 8.7%) and hypotransferrenemia (1/23; 4.3%). A total of 6 of 22 cases (27%) revealed a non-iron metabolism gene defect on whole-exome sequencing. Eleven novel variants (including variants of uncertain significance) were noted in 13 cases. Genotype-phenotype correlation revealed a significant association of frameshift/nonsense/splice variants with lower presentation age (0.8 months versus 9 years; P < 0.01) compared with missense variants. The systematic evaluation helped uncover an inherited iron defect in 41% (17/41) of cases, suggesting the need for active screening and awareness for these rare diseases in an iron-deficient endemic population.

Sequence introgression from exogenous lineages underlies genomic and biological differences among Cryptosporidium parvum IOWA lines

The IOWA strain of Cryptosporidium parvum is widely used in studies of the biology and detection of the waterborne pathogens Cryptosporidium spp. While several lines of the strain have been sequenced, IOWA-II, the only reference of the original subtype (IIaA15G2R1), exhibits significant assembly errors. Here we generated a fully assembled genome of IOWA-CDC of this subtype using PacBio and Illumina technologies. In comparative analyses of seven IOWA lines maintained in different laboratories (including two sequenced in this study) and 56 field isolates, IOWA lines (IIaA17G2R1) with less virulence had mixed genomes closely related to IOWA-CDC but with multiple sequence introgressions from IOWA-II and unknown lineages. In addition, the IOWA-IIaA17G2R1 lines showed unique nucleotide substitutions and loss of a gene associated with host infectivity, which were not observed in other isolates analyzed. These genomic differences among IOWA lines could be the genetic determinants of phenotypic traits in C. parvum. These data provide a new reference for comparative genomic analyses of Cryptosporidium spp. and rich targets for the development of advanced source tracking tools.

Listeria monocytogenes prevalence and genomic diversity along the pig and pork production chain

The facultative intracellular bacterium Listeria monocytogenes (L. monocytogenes) is the causative agent of listeriosis, a severe invasive illness. This ubiquitous species is widely distributed in the environment, but infection occurs almost exclusively through ingestion of contaminated food. The pork production sector has been heavily affected by a series of L. monocytogenes-related foodborne outbreaks in the past around the world. Ready-to-eat (RTE) pork products represent one of the main food sources for strong-evidence listeriosis outbreaks. This pathogen is known to be present throughout the entire pig and pork production chain. Some studies hypothesized that the main source of contamination in final pork products was either living pigs or the food-processing environment. A detailed genomic picture of L. monocytogenes can provide a renewed understanding of the routes of contamination from pig farms to the final products. This review provides an overview of the prevalence, the genomic diversity and the genetic background linked to virulence of L. monocytogenes along the entire pig and pork production chain, from farm to fork.

BioProtIS: Streamlining protein-ligand interaction pipeline for analysis in genomic and transcriptomic exploration

The identification of protein-ligand interactions plays a pivotal role in elucidating biological processes and discovering potential bioproducts. Harnessing the capabilities of computational methods in drug discovery, we introduce an innovative Inverted Virtual Screening (IVS) pipeline. This pipeline Integrated molecular dynamics and docking analyses to ensure that protein structures are not only energetically favorable but also representative of stable conformations. The primary objective of this pipeline is to automate and streamline the analysis of protein-ligand interactions at both genomic and transcriptomic scales. In the contemporary post-genomic era, high-throughput computational screening for bioproducts, biological systems, and therapeutic drugs has become a cornerstone practice. This approach offers the promise of cost-effectiveness, time efficiency, and optimization of laboratory work. Nevertheless, a notable deficiency persists in the availability of efficient pipelines capable of automating the virtual screening process, seamlessly integrating input and output, and leveraging the full potential of open-source tools. To bridge this critical gap, we have developed a versatile pipeline known as BioProtIS. This tool seamlessly integrates a suite of state-of-the-art tools, including Modeller, AlphaFold, Gromacs, FPOCKET, and AutoDock Vina, thus facilitating the streamlined docking of ligands with an expansive repertoire of proteins sourced from genomes and transcriptomes, and substrates. To assess the pipeline's performance, we employed the transcriptomes of Cereus jamacaru (a cactus species) and Aspisoma lineatum (firefly), along with the genome of Homo sapiens. This integration not only improves the accuracy of ligand-protein interactions by minimizing replicability deviations but also optimizes the discovery process by enabling the simultaneous evaluation of multiple substrates. Furthermore, our pipeline accommodates distinct testing scenarios, such as blind docking or site-specific targeting, which are invaluable in applications ranging from drug repositioning to the exploration of new allosteric binding sites and toxicity assessments. BioProtIS has been designed with modularity at its core. This inherent flexibility empowers users to make custom modifications directly within the source code, tailoring the pipeline to their specific research needs. Moreover, it lays the foundation for seamless integration of diverse docking algorithms in future iterations, promising ongoing advancements in the field of computational biology. This pipeline is available for free distribution and can be download at: https://github.com/BBMDO/BioProtIS.

Comparative genomic analysis and functional investigations for MCs catabolism mechanisms and evolutionary dynamics of MCs-degrading bacteria in ecology

Microcystins (MCs) significantly threaten the ecosystem and public health. Biodegradation has emerged as a promising technology for removing MCs. Many MCs-degrading bacteria have been identified, including an indigenous bacterium Sphingopyxis sp. YF1 that could degrade MC-LR and Adda completely. Herein, we gained insight into the MCs biodegradation mechanisms and evolutionary dynamics of MCs-degrading bacteria, and revealed the toxic risks of the MCs degradation products. The biochemical characteristics and genetic repertoires of strain YF1 were explored. A comparative genomic analysis was performed on strain YF1 and six other MCs-degrading bacteria to investigate their functions. The degradation products were investigated, and the toxicity of the intermediates was analyzed through rigorous theoretical calculation. Strain YF1 might be a novel species that exhibited versatile substrate utilization capabilities. Many common genes and metabolic pathways were identified, shedding light on shared functions and catabolism in the MCs-degrading bacteria. The crucial genes involved in MCs catabolism mechanisms, including mlr and paa gene clusters, were identified successfully. These functional genes might experience horizontal gene transfer events, suggesting the evolutionary dynamics of these MCs-degrading bacteria in ecology. Moreover, the degradation products for MCs and Adda were summarized, and we found most of the intermediates exhibited lower toxicity to different organisms than the parent compound. These findings systematically revealed the MCs catabolism mechanisms and evolutionary dynamics of MCs-degrading bacteria. Consequently, this research contributed to the advancement of green biodegradation technology in aquatic ecology, which might protect human health from MCs.

Phylogeography and demographic history of macaques, fascicularis species group, in East Asia: Inferred from multiple genomic markers

Climate changes at larger scales have influenced dispersal and range shifts of many taxa in East Asia. The fascicularis species group of macaques is composed of four species and is widely distributed in Southeast and East Asia. However, its phylogeography and demographic histories are currently poorly understood. Herein, we assembled autosomal, mitogenome, and Y-chromosome data for 106 individuals, and combined them with 174 mtDNA dloop haplotypes of this species group, with particular focus on the demographic histories and dispersal routes of Macaca fuscata, M. cyclopis, and M. mulatta. The results showed: (1) three monophyletic clades for M. fuscata, M. cyclopis, and M. mulatta based on the multiple genomics analyses; (2) the disparate demographic trajectories of the three species after their split ∼1.0 Ma revealed that M. cyclopis and M. fuscata were derived from an ancestral M. mulatta population; (3) the speciation time of M. cyclopis was later than that of M. fuscata, and their divergence time occurred at the beginning of "Ryukyu Coral Sea Stage" (1.0-0.2 Ma) when the East China Sea land bridge was completely submerged by the sea level rose; and (4) the three parallel rivers (Nujiang, Lancangjiang, and Jinshajiang) of Southwestern China divided M. mulatta into Indian and Chinese genetic populations ∼200 kya. These results shed light on understanding not only the evolutionary history of the fascicularis species group but also the formation mechanism of faunal diversity in East Asia during the Pleistocene.

Primary cardiac sarcoma: demographics, genomic study correlation, and survival benefits of surgery with adjuvant therapy in U.S. population

BACKGROUND

Cardiac sarcomas are rare and aggressive tumors with little known about the demographics, genetics, or treatment outcomes.

OBJECTIVES

The objectives of this study were to characterize the demographics, treatment modality, and survival associated with cardiac sarcomas and evaluate the potential for mutation-directed therapies.

METHODS

All cases from 2000 to 2018 of cardiac sarcoma were extracted from the SEER database. Genomic comparison utilized The Cancer Genome Atlas (TCGA) database, as well as reviews and re-analysis of past applicable genomic studies.

RESULTS

Cardiac sarcomas occurred most often in White patients, compared with national census data cardiac sarcomas occurred at a significantly higher rate in Asians. The majority of cases were undifferentiated (61.7%) and without distant metastases (71%). Surgery was the most common primary treatment modality and offered survival benefit (HR 0.391 (p = 0.001) that was most pronounced and sustained as compared to patients who received chemotherapy (HR 0.423 (p = 0.001) or radiation (HR 0.826 (p = 0.241) monotherapy. There was no difference in survival when stratified by race or sex; however, younger patients (< 50) had better survival. Genomics data on histologically undifferentiated cardiac sarcomas revealed a significant number were likely poorly differentiated pulmonary intimal sarcomas and angiosarcomas.

CONCLUSIONS

Cardiac sarcoma is a rare disease with surgery continuing to be a cornerstone of therapy followed by traditional chemotherapy. Case studies have indicated the potential for therapies directed to specific genetic aberrations to improve survival for these patients and utilization of next-generation sequencing (NGS) will help improve both classification and these therapies for cardiac sarcoma patients.

The mycorrhizal symbiosis: research frontiers in genomics, ecology, and agricultural application

Mycorrhizal symbioses between plants and fungi are vital for the soil structure, nutrient cycling, plant diversity, and ecosystem sustainability. More than 250 000 plant species are associated with mycorrhizal fungi. Recent advances in genomics and related approaches have revolutionized our understanding of the biology and ecology of mycorrhizal associations. The genomes of 250+ mycorrhizal fungi have been released and hundreds of genes that play pivotal roles in regulating symbiosis development and metabolism have been characterized. rDNA metabarcoding and metatranscriptomics provide novel insights into the ecological cues driving mycorrhizal communities and functions expressed by these associations, linking genes to ecological traits such as nutrient acquisition and soil organic matter decomposition. Here, we review genomic studies that have revealed genes involved in nutrient uptake and symbiosis development, and discuss adaptations that are fundamental to the evolution of mycorrhizal lifestyles. We also evaluated the ecosystem services provided by mycorrhizal networks and discuss how mycorrhizal symbioses hold promise for sustainable agriculture and forestry by enhancing nutrient acquisition and stress tolerance. Overall, unraveling the intricate dynamics of mycorrhizal symbioses is paramount for promoting ecological sustainability and addressing current pressing environmental concerns. This review ends with major frontiers for further research.

RIDGE, a tool tailored to detect gene flow barriers across species pairs

Characterizing the processes underlying reproductive isolation between diverging lineages is central to understanding speciation. Here, we present RIDGE-Reproductive Isolation Detection using Genomic polymorphisms-a tool tailored for quantifying gene flow barrier proportion and identifying the relevant genomic regions. RIDGE relies on an Approximate Bayesian Computation with a model-averaging approach to accommodate diverse scenarios of lineage divergence. It captures heterogeneity in effective migration rate along the genome while accounting for variation in linked selection and recombination. The barrier detection test relies on numerous summary statistics to compute a Bayes factor, offering a robust statistical framework that facilitates cross-species comparisons. Simulations revealed RIDGE's efficiency in capturing signals of ongoing migration. Model averaging proved particularly valuable in scenarios of high model uncertainty where no migration or migration homogeneity can be wrongly assumed, typically for recent divergence times <0.1 2Ne generations. Applying RIDGE to four published crow data sets, we first validated our tool by identifying a well-known large genomic region associated with mate choice patterns. Second, while we identified a significant overlap of outlier loci using RIDGE and traditional genomic scans, our results suggest that a substantial portion of previously identified outliers are likely false positives. Outlier detection relies on allele differentiation, relative measures of divergence and the count of shared polymorphisms and fixed differences. Our analyses also highlight the value of incorporating multiple summary statistics including our newly developed outlier ones that can be useful in challenging detection conditions.

Real-time genomic epidemiologic investigation of a multispecies plasmid-associated hospital outbreak of NDM-5-producing Enterobacterales infections

OBJECTIVES

New Delhi metallo-β-lactamase (NDM) is an emergent mechanism of carbapenem resistance associated with high mortality and limited treatment options. Because the blaNDM resistance gene is often carried on plasmids, traditional infection prevention and control (IP&C) surveillance methods and reactive whole genome sequencing (WGS) may not detect plasmid transfer in multispecies outbreaks.

METHODS

Initial outbreak detection of NDM-producing Enterobacterales identified at an acute care hospital occurred via traditional IP&C methods and was supplemented by real-time WGS surveillance performed weekly. To resolve NDM-encoding plasmids, we performed long-read sequencing and constructed hybrid assemblies. WGS data for suspected outbreaks was shared with the IP&C team for assessment and intervention.

RESULTS

We observed a multispecies outbreak of NDM-5-producing Enterobacterales isolated from 15 patients between February 2021 and February 2023. The 19 clinical and surveillance isolates sequenced included 7 bacterial species encoding the same NDM-5 plasmid. WGS surveillance and epidemiologic investigation characterized 10 horizontal plasmid transfer events and 6 bacterial transmission events between patients in varying hospital units.

CONCLUSIONS

Our investigation revealed a complex, multispecies outbreak of NDM involving multiple plasmid transfer and bacterial transmission events. We highlight the utility of combining traditional IP&C and prospective genomic methods in identifying and containing plasmid-associated outbreaks.

HER2 expression and genOmic characterization of rESected brain metastases from colorectal cancer: the HEROES study

BACKGROUND

Little is known about prognostic factors of brain metastases (BM) from colorectal cancer (CRC). HER2 amplification/overexpression (HER2+) was previously described; its impact on prognosis remains uncertain.

METHODS

In the translational study HEROES, extensive molecular analysis was performed on primary CRC (prCRC) and their matched resected BM by means of NGS comprehensive genomic profiling and HER2 status as assessed by immunohistochemical/ in situ hybridization. Count of tumour-infiltrating lymphocytes (TILs) was also performed.

PRIMARY OBJECTIVE

to describe the molecular landscape of paired BM/prCRC.

SECONDARY OBJECTIVES

to search for new prognostic biomarkers of outcome after BM resection: intracranial-only Progression-Free Survival (BM-iPFS), Progression-Free Survival (BM-PFS), and Overall Survival (BM-OS).

RESULTS

Out of 22 patients having paired samples of prCRC and BM, HER2+ was found on 4 (18%) BM, 3 (75%) of which also HER2+ in matched prCRC. Lower tumour mutation burden (HR 3.08; 95%CI 1.06-8.93; p = 0.0386) and HER2-negative BM (HER2neg) (HR 7.75;95%CI 1.97-30.40; p = 0.0033) were associated with longer BM-iPFS; HER2neg BM (HR 3.44; 95%CI 1.03-11.53; p = 0.0449) and KRASmut BM (HR 0.31; 95%CI 0.12-0.80; p = 0.0153) conferred longer BM-PFS. Longer BM-OS was found in pts with TILs-enriched (≥1.6/HPF) BM (HR 0.11; 95%CI0.01-0.91; p = 0.0403).

CONCLUSIONS

This study shows HER2+ enrichment in both BM and their prCRC. TILs-enriched BM conferred better BM-OS.

Genomic basis of seed colour in quinoa inferred from variant patterns using extreme gradient boosting

Quinoa is an agriculturally important crop species originally domesticated in the Andes of central South America. One of its most important phenotypic traits is seed colour. Seed colour variation is determined by contrasting abundance of betalains, a class of strong antioxidant and free radicals scavenging colour pigments only found in plants of the order Caryophyllales. However, the genetic basis for these pigments in seeds remains to be identified. Here we demonstrate the application of machine learning (extreme gradient boosting) to identify genetic variants predictive of seed colour. We show that extreme gradient boosting outperforms the classical genome-wide association approach. We provide re-sequencing and phenotypic data for 156 South American quinoa accessions and identify candidate genes potentially controlling betalain content in quinoa seeds. Genes identified include novel cytochrome P450 genes and known members of the betalain synthesis pathway, as well as genes annotated as being involved in seed development. Our work showcases the power of modern machine learning methods to extract biologically meaningful information from large sequencing data sets.

Genome-wide association studies for diarrhoea outcomes identified genomic regions affecting resistance to a severe enteropathy in suckling rabbits

Selection and breeding strategies to improve resistance to enteropathies are essential to reaching the sustainability of the rabbit production systems. However, disease heterogeneity (having only as major visible symptom diarrhoea) and low disease heritability are two barriers for the implementation of these strategies. Diarrhoea condition can affect rabbits at different life stages, starting from the suckling period, with large negative economic impacts. In this study, from a commercial population of suckling rabbits (derived from 133 litters) that experienced an outbreak of enteropathy, we first selected a few animals that died with severe symptoms of diarrhoea and characterized their microbiota, using 16S rRNA gene sequencing data. Clostridium genus was consistently present in all affected specimens. In addition, with the aim to identify genetic markers in the rabbit genome that could be used as selection tools, we performed genome-wide association studies for symptoms of diarrhoea in the same commercial rabbit population. These studies were also complemented with FST analyses between the same groups of rabbits. A total of 332 suckling rabbits (151 with severe symptoms of diarrhoea, 42 with mild symptoms and 129 without any symptoms till the weaning period), derived from 45 different litters (a subset of the 133 litters) were genotyped with the Affymetrix Axiom OrcunSNP Array. In both genomic approaches, rabbits within litters were paired to constitute two groups (susceptible and resistant, including the mildly affected in one or the other group) and run case and control genome-wide association analyses. Genomic heritability estimated in the designed experimental structure integrated in a commercial breeding scheme was 0.19-0.21 (s.e. 0.09-0.10). A total of eight genomic regions on rabbit chromosome 2 (OCU2), OCU3, OCU7, OCU12, OCU13, OCU16 and in an unassembled scaffold had significant single nucleotide polymorphisms (SNPs) and/or markers that trespassed the FST percentile distribution. Among these regions, three main peaks of SNPs were identified on OCU12, OCU13 and OCU16. The QTL region on OCU13 encompasses several genes that encode members of a family of immunoglobulin Fc receptors (FCER1G, FCRLA, FCRLB and FCGR2A) involved in the immune innate system, which might be important candidate genes for this pathogenic condition. The results obtained in this study demonstrated that resistance to an enteropathy occurring in suckling rabbits is in part genetically determined and can be dissected at the genomic level, providing DNA markers that could be used in breeding programmes to increase resistance to enteropathies in meat rabbits.

Genomic investigation on genes related to mercury metabolism in Amazonian indigenous populations

Studies have identified elevated levels of mercury in Amazonian Indigenous individuals, highlighting them as one of the most exposed to risks. In the unique context of the Brazilian Indigenous population, it is crucial to identify genetic variants with clinical significance to better understand vulnerability to mercury and its adverse effects. Currently, there is a lack of research on the broader genomic profile of Indigenous people, particularly those from the Amazon region, concerning mercury contamination. Therefore, the aim of this study was to assess the genomic profile related to the processes of mercury absorption, distribution, metabolism, and excretion in 64 Indigenous individuals from the Brazilian Amazon. We aimed to determine whether these individuals exhibit a higher susceptibility to mercury exposure. Our study identified three high-impact variants (GSTA1 rs1051775, GSTM1 rs1183423000, and rs1241704212), with the latter two showing a higher frequency in the study population compared to global populations. Additionally, we discovered seven new variants with modifier impact and a genomic profile different from the worldwide populations. These genetic variants may predispose the study population to more harmful mercury exposure compared to global populations. As the first study to analyze broader genomics of mercury metabolism pathways in Brazilian Amazonian Amerindians, we emphasize that our research aims to contribute to public policies by utilizing genomic investigation as a method to identify populations with a heightened susceptibility to mercury exposure.

Advances towards genome-based acute myeloid leukemia classification: A comparative analysis of WHO-HAEM4R, WHO-HAEM5, and International Consensus Classification

Two recent guidelines, the 5th edition of the World Health Organization Classification of Haematolymphoid Tumours (WHO-HAEM5) and the International Consensus Classification (ICC), were published to refine the diagnostic criteria of acute myeloid leukemia (AML). They both consider genomic features more extensively and expand molecularly defined AML subtypes. In this study, we compared the classifications of 1135 AML cases under both criteria. According to WHO-HAEM5 and ICC, the integration of whole transcriptome sequencing, targeted gene mutation screening, and conventional cytogenetic analysis identified defining genetic abnormalities in 89% and 90% of AML patients, respectively. The classifications displayed discrepancies in 16% of AML cases after being classified using the two guidelines, respectively. Both new criteria significantly reduce the number of cases defined by morphology and differentiation. However, their clinical implementation heavily relies on comprehensive and sophisticated genomic analysis, including genome and transcriptome levels, alongside the assessment of pathogenetic somatic and germline variations. Discrepancies between WHO-HAEM5 and ICC, such as the assignment of RUNX1 mutations, the rationality of designating AML with mutated TP53 as a unique entity, and the scope of rare genetic fusions, along with the priority of concurrent AML-defining genetic abnormalities, are still pending questions requiring further research for more elucidated insights.

Exploring genomic inbreeding and selection signatures in a commercial Brangus herd through functional annotation

Composite breeds, including Brangus, are widely utilized in subtropical and tropical regions to harness the advantages of both Bos t. taurus and Bos t. indicus breeds. The formation and subsequent selection of composite breeds may result in discernible signatures of selection and shifts in genomic population structure. The objectives of this study were to 1) assess genomic inbreeding, 2) identify signatures of selection, 3) assign functional roles to these signatures in a commercial Brangus herd, and 4) contrast signatures of selection between selected and non-selected cattle from the same year. A total of 4035 commercial Brangus cattle were genotyped using the GGP-F250K array. Runs of Homozygosity (ROH) were used to identify signatures of selection and calculate genomic inbreeding. Quantitative trait loci (QTL) enrichment analysis and literature search identified phenotypic traits linked to ROH islands. Genomic inbreeding averaged 5%, primarily stemming from ancestors five or more generations back. A total of nine ROH islands were identified, QTL enrichment analysis revealed traits related to growth, milk composition, carcass, reproductive, and meat quality traits. Notably, the ROH island on BTA14 encompasses the pleiomorphic adenoma (PLAG1) gene, which has been linked to growth, carcass, and reproductive traits. Moreover, ROH islands associated with milk yield and composition were more pronounced in selected replacement heifers of the population, underscoring the importance of milk traits in cow-calf production. In summary, our research sheds light on the changing genetic landscape of the Brangus breed due to selection pressures and reveals key genomic regions impacting production traits.

Precision Medicine in Veterinary Science

Precision medicine focuses on the clinical management of the individual patient, not on population-based findings. Successes from human precision medicine inform veterinary oncology. Early evidence of success for canines shows how precision medicine can be integrated into practice. Decreasing genomic profiling costs will allow increased utilization and subsequent improvement of knowledge base from which to make better informed decisions. Utility of precision medicine in canine oncology will only increase for improved cancer characterization, enhanced therapy selection, and overall more successful management of canine cancer. As such, practitioners are called to interpret and leverage precision medicine reports for their patients.

Early Onset Intrahepatic Cholangiocarcinoma: Clinical Characteristics, Oncological Outcomes, and Genomic/Transcriptomic Features

INTRODUCTION

Data on clinical characteristics and disease-specific prognosis among patients with early onset intrahepatic cholangiocarcinoma (ICC) are currently limited.

METHODS

Patients undergoing hepatectomy for ICC between 2000 and 2020 were identified by using a multi-institutional database. The association of early (≤50 years) versus typical onset (>50 years) ICC with recurrence-free (RFS) and disease-specific survival (DSS) was assessed in the multi-institutional database and validated in an external cohort. The genomic and transcriptomic profiles of early versus late onset ICC were analyzed by using the Total Cancer Genome Atlas (TCGA) and Memorial Sloan Kettering Cancer Center databases.

RESULTS

Among 971 patients undergoing resection for ICC, 22.7% (n = 220) had early-onset ICC. Patients with early-onset ICC had worse 5-year RFS (24.1% vs. 29.7%, p < 0.05) and DSS (36.5% vs. 48.9%, p = 0.03) compared with patients with typical onset ICC despite having earlier T-stage tumors and lower rates of microvascular invasion. In the validation cohort, patients with early-onset ICC had worse 5-year RFS (7.4% vs. 20.5%, p = 0.002) compared with individuals with typical onset ICC. Using the TCGA cohort, 652 and 266 genes were found to be upregulated (including ATP8A2) and downregulated (including UTY and KDM5D) in early versus typical onset ICC, respectively. Genes frequently implicated as oncogenic drivers, including CDKN2A, IDH1, BRAF, and FGFR2 were infrequently mutated in the early-onset ICC patients.

CONCLUSIONS

Early-onset ICC has distinct clinical and genomic/transcriptomic features. Morphologic and clinicopathologic characteristics were unable to fully explain differences in outcomes among early versus typical onset ICC patients. The current study offers a preliminary landscape of the molecular features of early-onset ICC.