Genomic Analysis Research Articles

Genomic Analysis of Aspergillus Section Terrei Reveals a High Potential in Secondary Metabolite Production and Plant Biomass Degradation

Aspergillus terreus has attracted interest due to its application in industrial biotechnology, particularly for the production of itaconic acid and bioactive secondary metabolites. As related species also seem to possess a prosperous secondary metabolism, they are of high interest for genome mining and exploitation. Here, we present draft genome sequences for six species from Aspergillus section Terrei and one species from Aspergillus section Nidulantes. Whole-genome phylogeny confirmed that section Terrei is monophyletic. Genome analyses identified between 70 and 108 key secondary metabolism genes in each of the genomes of section Terrei, the highest rate found in the genus Aspergillus so far. The respective enzymes fall into 167 distinct families with most of them corresponding to potentially unique compounds or compound families. Moreover, 53% of the families were only found in a single species, which supports the suitability of species from section Terrei for further genome mining. Intriguingly, this analysis, combined with heterologous gene expression and metabolite identification, suggested that species from section Terrei use a strategy for UV protection different to other species from the genus Aspergillus. Section Terrei contains a complete plant polysaccharide degrading potential and an even higher cellulolytic potential than other Aspergilli, possibly facilitating additional applications for these species in biotechnology.

Survival Strategy, Metabolic Potential and Taxonomic Reframe of Kocuria polaris

Antarctica is renowned as the most inhospitable environment where microorganisms are thriving in the frontiers of life. In the past few years, many novel bacterial species have been reported from the Antarctic environment. During taxonomic re-evaluation of novel bacterial species from Antarctica, it was noticed that Kocuria polaris shared high 16S rRNA gene sequence similarity with Kocuria rosea. In the present study, the taxonomic position, metabolic potentials, and stress survival strategy of K. polaris were evaluated through genome analysis. K. polaris encodes genes for glycolysis, citrate cycle, pentose phosphate pathway, dissimilatory nitrate reduction, assimilatory sulfate reduction, etc. In addition, K. polaris also encodes genes for cold and salt stress. The 16S rRNA gene sequence extracted from K. polaris and K. rosea genomes showed 99.7% similarity. In the phylogenomic tree, K. polaris and K. rosea clustered together. The average nucleotide identity and digital DNA–DNA hybridization values between K. polaris and K. rosea exceeded the threshold (95-96% for ANI and 70% for dDDH) value for distinguishing species, showing that they are similar species. The present study shed light on K. polaris survival strategy in extreme conditions. We further propose to reclassify Kocuria polaris as a later heterotypic synonym of Kocuria rosea.

AMRomics: a scalable workflow to analyze large microbial genome collections

Whole genome analysis for microbial genomics is critical to studying and monitoring antimicrobial resistance strains. The exponential growth of microbial sequencing data necessitates a fast and scalable computational pipeline to generate the desired outputs in a timely and cost-effective manner. Recent methods have been implemented to integrate individual genomes into large collections of specific bacterial populations and are widely employed for systematic genomic surveillance. However, they do not scale well when the population expands and turnaround time remains the main issue for this type of analysis. Here, we introduce AMRomics, an optimized microbial genomics pipeline that can work efficiently with big datasets. We use different bacterial data collections to compare AMRomics against competitive tools and show that our pipeline can generate similar results of interest but with better performance. The software is open source and is publicly available at https://github.com/amromics/amromics under an MIT license.

Epigenetic modification of PHLDA2 is associated with tumor microenvironment and unfavorable outcome of immune checkpoint inhibitor-based therapies in clear cell renal cell carcinoma

BackgroundA substantial proportion of patients with metastatic clear cell renal cell carcinoma (ccRCC) cannot derive benefit from immune checkpoint inhibitor (ICI) plus anti-angiogenic agent combination therapy, making identification of predictive biomarkers an urgent need. The members of pleckstrin homology-like domain family A (PHLDA) play critical roles in multiple cancers, whereas their roles in ccRCC remain unknown.MethodsTranscriptomic, clinical, genetic alteration and DNA methylation data were obtained for integrated analyses from TCGA database. RNA sequencing was performed on 117 primary tumors and 79 normal kidney tissues from our center. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis, gene set enrichment analysis were performed to explore transcriptomic features. Data from three randomized controlled trials (RCT), including CheckMate025, IMmotion151, JAVELIN101, were obtained for validation.ResultsMembers of PHLDA family were dysregulated in pan-cancer. Elevated PHLDA2 expression was associated with adverse clinicopathologic parameters and worse prognosis in ccRCC. Aberrant DNA hypomethylation contributed to up-regulation of PHLDA2. An immunosuppressive microenvironment featured by high infiltrates of Tregs and cancer-associated fibroblasts, was observed in ccRCC with higher PHLDA2 expression. Utilizing data from three RCTs, the association of elevated PHLDA2 expression with poor therapeutic efficacy of ICI plus anti-angiogenic combination therapy was confirmed.ConclusionsOur study revealed that elevated PHLDA2 expression regulated by DNA hypomethylation was correlated with poor prognosis and immunosuppressive microenvironment, and highlighted the role of PHLDA2 as a robust biomarker for predicting therapeutic efficacy of ICI plus anti-angiogenic agent combination therapy in ccRCC, which expand the dimension of precision medicine.

Safety evaluation of Bifidobacterium animalis subsp. lactis BLa80 under in vitro, and in vivo conditions

Bifidobacterium animalis subsp. lactis BLa80 is a new probiotic strain with extensive applications in food products both domestically and internationally. Given the rising consumption of this probiotic, its safety assessment is increasingly crucial in the food industry. This study evaluates the safety of strain BLa80 using a combination of in vitro and in vivo assays along with genomic analysis. Methods included exposing the strain to artificial gastric and intestinal fluids, as well as a medium containing bile salts, to stimulate human digestive conditions. The strain showed high tolerance to gastric fluid at pH of 2.5 and to 0.3 % bile salts. It maintained a 99.92 % survival rate in intestinal fluid. Additional tests assessed hemolytic activity, antibiotic susceptibility (revealing sensitivity to 7 antibiotics), and biogenic amine production using HPLC-ELSD, confirming the absence of histamine, and other harmful amines. Bile salt hydrolase activity was demonstrated qualitatively, and metabolic byproducts were quantitatively analyzed using a D-/l-lactic acid assay kit, showing that BLa80 produces 1.48 mg/mL of l-lactic acid and no harmful d-lactic acid. Genomic analysis confirmed the absence of virulence or pathogenicity genes, and a 90-day oral toxicity study in rats confirmed no toxic effects at various doses. Overall, these findings support the safety classification of the strain BLa80.

Therapeutic and diagnostic applications of exosomes in colorectal cancer.

Exosomes play a key role in colorectal cancer (CRC) related processes. This review explores the various functions of exosomes in CRC and their potential as diagnostic markers, therapeutic targets, and drug delivery vehicles. Exosomal long non-coding RNAs (lncRNAs) and microRNAs (miRNAs) significantly influence CRC progression. Specific exosomal lncRNAs are linked to drug resistance and tumor growth, respectively, highlighting their therapeutic potential. Similarly, miRNAs like miR-21, miR-10b, and miR-92a-3p, carried by exosomes, contribute to chemotherapy resistance by altering signaling pathways and gene expression in CRC cells. The review also discusses exosomes' utility in CRC diagnosis. Exosomes from cancer cells have distinct molecular signatures compared to healthy cells, making them reliable biomarkers. Specific exosomal lncRNAs (e.g., CRNDE-h) and miRNAs (e.g., miR-17-92a) have shown effectiveness in early CRC detection and monitoring of treatment responses. Furthermore, exosomes show promise as vehicles for targeted drug delivery. The potential of mesenchymal stem cell (MSC)-derived exosomes in CRC treatment is also noted, with their role varying from promoting to inhibiting tumor progression. The application of multi-omics approaches to exosome research is highlighted, emphasizing the potential for discovering novel CRC biomarkers through comprehensive genomic, transcriptomic, proteomic, and metabolomic analyses. The review also explores the emerging field of exosome-based vaccines, which utilize exosomes' natural properties to elicit strong immune responses. In conclusion, exosomes represent a promising frontier in CRC research, offering new avenues for diagnosis, treatment, and prevention. Their unique properties and versatile functions underscore the need for continued investigation into their clinical applications and underlying mechanisms.

Isolation and Genomic Analysis of a Case of Staphylococcus argenteus ST2250 Related to Sepsis in Italy

Staphylococcus argenteus, identified in 2006, represents a challenging case of bacterial taxonomic identification because of its high similarity to Staphylococcus aureus. In this context, neither mass spectrometry (MS) nor 16S gene analysis cannot precisely reveal the difference between the two species. In our study, the sensitivity to antibiotics of S. argenteus isolated from blood culture was tested, and the investigation of the bacterial genome was performed by Multi-Locus Sequence Typing (MLST) and Whole-Genome Next-Generation Sequencing (WG-NGS). The pathogen was identified as ST2250 and presented perfectly matched resistance genes, namely aph(3′)-III, mgrA, and sepA, whereas the virulence gene detected was scn. Two plasmids were found: the pSAS plasmid, belonging to the family of Inc18, and plasmid pN315, belonging to the Rep3 group. The epidemiological distribution and the spread of S. argenteus infection are scarcely documented, particularly when associated with sepsis. Therefore, a correct taxonomy identification, antibiogram, and resistance gene analysis may help in acquiring knowledge about this bacterium and implement its detection and treatment.

3D genomic analysis reveals novel enhancer-hijacking caused by complex structural alterations that drive oncogene overexpression

Cancer genomes are composed of many complex structural alterations on chromosomes and extrachromosomal DNA (ecDNA), making it difficult to identify non-coding enhancer regions that are hijacked to activate oncogene expression. Here, we describe a 3D genomics-based analysis called HAPI (Highly Active Promoter Interactions) to characterize enhancer hijacking. HAPI analysis of HiChIP data from 34 cancer cell lines identified enhancer hijacking events that activate both known and potentially novel oncogenes such as MYC, CCND1, ETV1, CRKL, and ID4. Furthermore, we found enhancer hijacking among multiple oncogenes from different chromosomes, often including MYC, on the same complex amplicons such as ecDNA. We characterized a MYC-ERBB2 chimeric ecDNA, in which ERBB2 heavily hijacks MYC’s enhancers. Notably, CRISPRi of the MYC promoter led to increased interaction of ERBB2 with MYC enhancers and elevated ERBB2 expression. Our HAPI analysis tool provides a robust strategy to detect enhancer hijacking and reveals novel insights into oncogene activation.

Whole genome sequence and comparative genomic analysis of novel Rickettsia koreansis strain CNH17-7 isolated from human.

To determine the genomic feature of novel spotted fever-causing Rickettsia koreansis strain CNH17-7, which is different from R. japonica that is a causative agent for Japanese spotted fever (JSF), and to perform its comparative genomic analysis. Whole genome sequencing (WGS) was performed on R. koreansis strain CNH17-7 by using the Illumina Miseq system. After WGS, assembly and annotation were done by SPAdes. Then, its genomic features were compared with 19 different Rickettsia species. Based on the average nucleotide identity (ANI) value, an unweighted pair group method with an arithmetic mean (UPGMA) dendrogram was generated. Following the dendrogram analysis, pan-and core-genome analysis was performed. Then additional comparative analyses with two genetically closest Rickettsia species were conducted based on gene repertoire. R. koreansis strain CNH17-7 has a chromosome consisting of 1,392,633bp with GC content of 32.4%. The ANI-derived UPGMA showed that R. koreansis strain CNH17-7 is genetically close to R. japonica YH and R. heilongjiangensis 054 but is distinctively differentiated. The ANI value of R. koreansis strain CNH17-7 to R. japonica YH and R. heilongjiangensis 054 are 98.14% and 98.04% respectively, indicating R. koreansis strain CNH17-7 is sufficient to be classified as a new species. Other than ANI, R. koreansis strain CNH17-7 also contains novel CDS and its COG functional category proportion which is distinct compared to R. japonica YH and R. heilongjiangensis 054. We have revealed genomic features of the novel R. koreansis strain CNH17-7. Hence, we propose R. koreansis strain CNH17-7 as new Rickettsia species.

Evaluation of genotype imputation using Glimpse tools on low coverage ancient DNA.

Ancient DNA provides a unique frame for directly studying human population genetics in time and space. Still, since most of the ancient genomic data is low coverage, analysis is confronted with a low number of SNPs, genotype uncertainties, and reference-bias. Here, we for the first time benchmark the two distinct versions of Glimpse tools on 120 ancient human genomes from Eurasia including those largely from previously under-evaluated regions and compare the performance of genotype imputation with de facto analysis approaches for low coverage genomic data analysis. We further investigate the impact of two distinct reference panels on imputation accuracy for low coverage genomic data. We compute accuracy statistics and perform PCA and f4-statistics to explore the behaviour of genotype imputation on low coverage data regarding (i)two versions of Glimpse, (ii)two reference panels, (iii)four post-imputation filters and coverages, as well as (iv)data type and geographical origin of the samples on the analyses. Our results reveal that even for 0.1X coverage ancient human genomes, genotype imputation using Glimpse-v2 is suitable. Additionally, using the 1000 Genomes merged with Human Genome Diversity Panel improves the accuracy of imputation for the rare variants with low MAF, which might be important not only for ancient genomics but also for modern human genomic studies based on low coverage data and for haplotype-based analysis. Most importantly, we reveal that genotype imputation of low coverage ancient human genomes reduces the genetic affinity of the samples towards human reference genome. Through solving one of the most challenging biases in data analysis, so-called reference bias, genotype imputation using Glimpse v2 is promising for low coverage ancient human genomic data analysis and for rare-variant-based and haplotype-based analysis.

Genomic Insights into High-Grade Infarct-Associated Bone Sarcomas

Sarcomas rarely develop in bones previously compromised by infarcts. These infarct-associated sarcomas often present as undifferentiated pleomorphic sarcomas (UPS), and their genetic characteristics are poorly understood. High-grade spindle cell/UPS of bone are typically treated with a combination of surgery and chemotherapy, similar to osteosarcoma. We conducted a detailed clinicopathologic and genomic analysis of six cases of intraosseous sarcomas arising from histologically and radiographically confirmed bone infarcts. We analyzed 523 genes for sequence-level mutations using next-generation sequencing with the TruSight Oncology 500 panel and utilized whole-genome SNP Microarray (OncoScan CNV) to detect copy number alterations and loss of heterozygosity (LOH). Genomic instability was assessed through Homologous Recombination Deficiency (HRD) metrics, incorporating LOH, telomeric allelic imbalance, and large-scale state transitions. FISH and immunohistochemistry validated the findings. The cohort included three men and three women, with a median age of 70, and tumors located in the femur and tibia. Five of the six patients developed distant metastases. Treatment involved surgery and chemotherapy or immune checkpoint inhibitors. Genomic analysis revealed significant complexity and high HRD scores, ranging from 32 to 57 (with a cut-off of 32). Chromosome 12 alterations, including segmental amplification or chromothripsis, were observed in four cases. Notably, MDM2 amplification, confirmed by FISH, was detected in two cases. Homozygous deletion of CDKN2A/B was observed in all six cases. Tumor mutational burden (TMB) levels ranged from 2.4 to 7.9 mutations per megabase. Notable pathogenic mutations included H3-3A mutations (p.G35R and p.G35W), and mutations in HRAS, DNMT3A, NF2, PIK3CA, POLE, and TP53, each in one case. These results suggest that high-grade infarct-associated sarcomas of bone, while sharing high levels of structural variations with osteosarcoma, may exhibit potentially less frequent TP53 mutations and more common CDKN2A/B deletions. This points to the possibility that the mutation spectrum and disrupted pathways could be distinct from conventional osteosarcoma.

Una comune caduta per una rara diagnosi: un caso di emicrania emiplegica

The Authors describe the case of a three-year-old girl with hemiplegic migraine triggered by mild trauma; the diagnosis was confirmed through genomic analysis (heterozygous missense variant c.1091>T in the ATP1A2 gene).

Diagnosis of Two Unrelated Syndromes of Prader-Willi and Calpainopathy: Insight from Trio Whole Genome Analysis and Isodisomy Mapping

Diagnosis of Two Unrelated Syndromes of Prader-Willi and Calpainopathy: Insight from Trio Whole Genome Analysis and Isodisomy Mapping

Synergic activity of FGFR2 and MEK inhibitors in the treatment of FGFR2-amplified cancers of unknown primary

Patients with cancer of unknown primary (CUP) carry the double burden of an aggressive disease and reduced access to therapies. Experimental models are pivotal for CUP biology investigation and drug testing. We derived two CUP cell lines (CUP#55 and #96), and corresponding patient-derived xenografts (PDXs), from ascites tumor cells. CUP cell lines and PDXs underwent histological, immune-phenotypical, molecular, and genomic characterization confirming the features of the original tumor. The tissue-of-origin prediction was obtained from the tumor microRNA expression profile and confirmed by single-cell transcriptomics. Genomic testing and FISH analysis identified FGFR2 gene amplification in both models, in the form of homogenously staining region (HSR) in CUP#55 and double minutes in CUP#96. FGFR2 was recognized as the main oncogenic driver and therapeutic target. FGFR2-targeting drug BGJ-398 (infigratinib) in combination with the MEK inhibitor trametinib proved to be synergic and exceptionally active, both in vitro and in vivo. The effects of the combined treatment by single-cell gene expression analysis revealed a remarkable plasticity of tumor cells and the greater sensitivity of cells with epithelial phenotype. This study brings personalized therapy closer to CUP patients and provides the rationale for FGFR2 and MEK targeting in metastatic tumors with FGFR2 pathway activation.

Accelerated mitochondrial evolution and asymmetric fitness of hybrids contribute to the persistence of Helix thessalica in the Helix pomatia range.

Interbreeding and introgression between recently diverged species is common. However, the processes that prevent these species from merging where they co-occur are not well understood. We studied the mechanisms that allowed an isolated group of populations of the snail Helix thessalica to persist within the range of the related Helix pomatia despite high gene flow. Using genomic cline analysis, we found that the nuclear gene flow between the two taxa across the mosaic hybrid zone was not different from that expected under neutral admixture, but that the exchange of mtDNA was asymmetric. Tests showed that there is relaxed selection in the mitochondrial genome of H. thessalica and that the substitution rate is elevated compared to that of H. pomatia. A lack of hybrids that combine the mtDNA of H. thessalica with a mainly (>46%) H. pomatia genomic background indicates that the nuclear-encoded mitochondrial proteins of H. pomatia are not well adapted to the more rapidly evolving proteins and RNAs encoded by the mitochondrion of H. thessalica. The presumed reduction of fitness of hybrids with the fast-evolving mtDNA of H. thessalica and a high H. pomatia ancestry, similar to 'Darwin's Corollary to Haldane's rule', resulted in a relative loss of H. pomatia nuclear ancestry compared to H. thessalica ancestry in the hybrid zone. This probably prevents the H. thessalica populations from merging quickly with the surrounding H. pomatia populations and supports the hypothesis that incompatibilities between rapidly evolving mitochondrial genes and nuclear genes contribute to speciation.

Allelic variations in the chpG effector gene within Clavibacter michiganensis populations determine pathogen host range.

Plant pathogenic bacteria often have a narrow host range, which can vary among different isolates within a population. Here, we investigated the host range of the tomato pathogen Clavibacter michiganensis (Cm). We determined the genome sequences of 40 tomato Cm isolates and screened them for pathogenicity on tomato and eggplant. Our screen revealed that out of the tested isolates, five were unable to cause disease on any of the hosts, 33 were exclusively pathogenic on tomato, and two were capable of infecting both tomato and eggplant. Through comparative genomic analyses, we identified that the five non-pathogenic isolates lacked the chp/tomA pathogenicity island, which has previously been associated with virulence in tomato. In addition, we found that the two eggplant-pathogenic isolates encode a unique allelic variant of the putative serine hydrolase chpG (chpGC), an effector that is recognized in eggplant. Introduction of chpGC into a chpG inactivation mutant in the eggplant-non-pathogenic strain Cm101, failed to complement the mutant, which retained its ability to cause disease in eggplant and failed to elicit hypersensitive response (HR). Conversely, introduction of the chpG variant from Cm101 into an eggplant pathogenic Cm isolate (C48), eliminated its pathogenicity on eggplant, and enabled C48 to elicit HR. Our study demonstrates that allelic variation in the chpG effector gene is a key determinant of host range plasticity within Cm populations.

A systematic analysis of GASA family members in the Salix matsudana genome: Characterization, expression profile, and putative function in antioxidation

The gibberellic acid stimulated Arabidopsis (GASA) gene family, specific to plants, plays a pivotal role in growth, development, and response to diverse environmental stresses. Despite their importance, GASA genes in willow trees have been little studied. Here, 35 SmGASAs were identified from the Salix matsudana genome, which exhibited similar physical and chemical proprieties to GASAs in other species. A phylogenetic analysis revealed that GASAs from Populus trichocarpa, Vitis vinifera, Nicotiana tabacum, Oryza sativa, Arabidopsis thaliana, and S. matsudana are categorized into three distinct subfamilies. It appeared that polyploidization and segmental duplication have contributed to the expansion and diversification of the SmGASA family. These genes were significantly enriched in response to GA and oxygen-containing compounds. The diverse expression patterns observed across willow varieties suggested functional variation of SmGASAs throughout their evolutionary history. The expression profiles of SmGASAs were further analyzed in response to different treatments (i.e., salt, submergence, and polyethylene glycol (PEG)), indicating a potential role for SmGASAs in adapting to these stressors in willows. Remarkably, SmGASAs were shown to detoxify hydrogen peroxide (H2O2), highlighting their potent antioxidant properties. Overall, the comprehensive genomic analysis of willow SmGASAs shed light on their role in stress responses, offering valuable insights for future research and supporting the hypothesis of their involvement in stress adaptation mechanisms.

Assembly and comparative genome analysis of four mitochondrial genomes from Saccharum complex species

Assembly and comparative genome analysis of four mitochondrial genomes from Saccharum complex species

Assembly and comparative analysis of the complete mitochondrial genome of Vaccinium carlesii Dunn

Vaccinium L. is an important fruit tree with nutritional, medicinal, and ornamental values. However, the mitochondrial (mt) genome of Vaccinium L. remains largely unexplored. Vaccinium carlesii Dunn is an endemic wild resource in China, which is crucial for blueberry breeding. The V. carlesii mt genomes were sequenced using Illumina and Nanopore, which total length was 636,904 bp with 37 protein coding genes, 20 tRNA genes, and three rRNA genes. We found four pairs of long repeat fragments homologous recombination mediated the generation of substructures in the V. carlesii mt genome. We predicted 383 RNA editing sites, all converting cytosine (C) to uracil (U). According to the phylogenetic analysis, V. carlesii and V. macrocarpon of the Ericaceae exhibited the closest genetic relationship. This study provides a theoretical basis for understanding the evolution of higher plants, species classification and identification, and will also be useful for further utilization of Vaccinium germplasm resources.

Co-expression of a pair of interdependent regulators coding genes ovmZ and ovmW awakens the production of angucyclinones antibiotics in Streptomyces neyagawaensis

BackgroundMicrobial genome sequencing and analysis revealed the presence of abundant silent secondary metabolites biosynthetic gene clusters (BGCs) in streptomycetes. Activating these BGCs has great significance for discovering new compounds and novel biosynthetic pathways.ResultsIn this study, we found that ovmZ and ovmW homologs, a pair of interdependent transcriptional regulators coding genes, are widespread in actinobacteria and closely associated with the biosynthesis of secondary metabolites. Through co-overexpression of native ovmZ and ovmW in Streptomyces neyagawaensis NRRL B-3092, a silent type II polyketide synthase (PKS) gene cluster was activated to produce gephyromycin A, tetrangomycin and fridamycin E with the yields of 22.3 ± 8.0 mg/L, 4.8 ± 0.5 mg/L and 20.3 ± 4.1 mg/L respectively in the recombinant strain of S.ne/pZnWn. However, expression of either ovmZ or ovmW failed to activate this gene cluster. Interestingly, overexpression of the heterologous ovmZ and ovmW pair from oviedomycin BGC of S. ansochromogenes 7100 also led to awakening of this silent angucyclinone BGC in S. neyagawaensis.ConclusionA silent angucyclinone BGC was activated by overexpressing both ovmZ and ovmW in S. neyagawaensis. Due to the wide distribution of ovmZ and ovmW in the BGCs of actinobacteria, co-overexpression of ovmZ and ovmW could be a strategy for activating silent BGCs, thus stimulating the biosynthesis of secondary metabolites.