Qualimap: evaluating next-generation sequencing alignment data

Host Dependent-Transposon for a Plasmid Found in Aeromonas salmonicida subsp. salmonicida That Bears a catB3 Gene for Chloramphenicol Resistance

Plasmids that carry antibiotic resistance genes occur frequently in Aeromonas salmonicida subsp. salmonicida, an aquatic pathogen with severe consequences in salmonid farming. Here, we describe a 67 kb plasmid found in the A. salmonicida subsp. salmonicida Strain SHY15-2939 from Quebec, Canada. This new plasmid, named pAsa-2939 and identified by high throughput sequencing, displays features never found before in this bacterial species. It contains a transposon related to the Tn21 family, but with an unusual organization. This transposon bears a catB3 gene (chloramphenicol resistance) that has not been detected yet in A. salmonicida subsp. salmonicida. The plasmid is transferable by conjugation into Aeromonas hydrophila, but not into Escherichia coli. Based on PCR analysis and genomic sequencing (Illumina and PacBio), we determined that the transposon is unstable in A. salmonicida subsp. salmonicida Strain SHY15-2939, but it is stable in A. hydrophila trans-conjugants, which explains the chloramphenicol resistance variability observed in SHY15-2939. These results suggest that this bacterium is likely not the most appropriate host for this plasmid. The presence of pAsa-2939 in A. salmonicida subsp. salmonicida also strengthens the reservoir role of this bacterium for antibiotic resistance genes, even those that resist antibiotics not used in aquaculture in Québec, such as chloramphenicol.

The effect of in vitro consecutive passages and culture medium on the genetic variations in BCG Pasteur 1173P2 vaccine

Since the introduction of the Bacillus Calmette-Guérin (BCG) vaccine, the genomes of vaccine strains have undergone variations due to repeated passages in different laboratories and vaccine production facilities. Genetic variations have been considered as one of the effective factors in the BCG variable protective efficacy. Consecutive subcultures have been shown to play an essential role in causing genetic variations in several microorganisms, including Mycobacterium bovis BCG. Therefore, the world health organization (WHO) recommendation to limit the passages of master seed lot in the BCG vaccine production should be considered. Besides, the role of other external variables such as quality of the raw ingredients of the culture media, the type of the culture medium and the cultivation methods in the vaccine production has been poorly studied. Here, the effect of passages and culture medium on genetic variations in a BCG seed lot was investigated during a year. The findings of this study revealed a total of 19 variants compared to seed lot while the passages were more than the number recommended by WHO. The first culture of seed lot in the Sauton broth and Middlebrook 7H9 media, and the last subculture in Sauton broth had the least and the most variants, respectively. The observation of the higher number of variants in the last cultures on Sauton broth and Middlebrook 7H9 in comparison to the first and the middle cultures may indicate the effect of passages on the genetic variations in BCG. Additionally, more variants in BCG grown in the Sauton broth do not necessarily represent the greater ability of this medium to cause genetic mutations. For a better conclusion, it is required to examine the medium components as independent variables.

Distinct Cellular Origins and Differentiation Process Account for Distinct Oncogenic and Clinical Behaviors of Leiomyosarcomas

In leiomyosarcoma (LMS), a very aggressive disease, a relatively transcriptionally uniform subgroup of well-differentiated tumors has been described and is associated with poor survival. The question raised how differentiation and tumor progression, two apparently antagonist processes, coexist and allow tumor malignancy. We first identified the most transcriptionally homogeneous LMS subgroup in three independent cohorts, which we named 'hLMS'. The integration of multi-omics data and functional analysis suggests that hLMS originate from vascular smooth muscle cells and show that hLMS transcriptional program reflects both modulations of smooth muscle contraction activity controlled by MYOCD/SRF regulatory network and activation of the cell cycle activity controlled by E2F/RB1 pathway. We propose that the phenotypic plasticity of vascular smooth muscle cells coupled with MYOCD/SRF pathway amplification, essential for hLMS survival, concomitant with PTEN absence and RB1 alteration, could explain how hLMS balance this uncommon interplay between differentiation and aggressiveness.

PIK3CA is recurrently mutated in canine mammary tumors, similarly to in human mammary neoplasia

Biological features of neoplastic disease affecting mammary gland tissue are shared between canines and humans. Research performed in either species has translational value and early phase clinical trials performed in canines with spontaneous disease could be informative for human trials. The purpose of this study was to investigate the somatic genetic aberrations occurring in canine mammary neoplasia by exome capture and next generation sequencing. Based on 55 tumor-normal pairs we identified the PIK3CA gene as the most commonly mutated gene in canine mammary tumors, with 25% of samples carrying mutations in this gene. A recurrent missense mutation was identified, p.H1047R, which is homologous to the human PIK3CA hotspot mutation found in different types of breast neoplasia. Mutations homologous to other known human mutation hotspots such as the PIK3CA p.E545K and the KRAS p.G12V/D were also identified. We identified copy number aberrations affecting important tumor suppressor and oncogenic pathways including deletions affecting the PTEN tumor suppressor gene. We suggest that activation of the KRAS or PIK3CA oncogenes or loss of the PTEN suppressor gene may be important for mammary tumor development in dogs. This data endorses the conservation of cancer across species and the validity of studying cancer in non-human species.

Circulating Cell-Free RNA in Blood as a Host Response Biomarker for the Detection of Tuberculosis

Tuberculosis (TB) remains a leading cause of death from an infectious disease worldwide. This is partly due to a lack of tools to effectively screen and triage individuals with potential TB. Whole blood RNA signatures have been extensively studied as potential biomarkers for TB, but they have failed to meet the World Health Organization's (WHOs) target product profiles (TPPs) for a non-sputum triage or diagnostic test. In this study, we investigated the utility of plasma cell-free RNA (cfRNA) as a host response biomarker for TB. We used RNA profiling by sequencing to analyze plasma samples from 182 individuals with a cough lasting at least two weeks, who were seen at outpatient clinics in Uganda, Vietnam, and the Philippines. Of these individuals, 100 were diagnosed with microbiologically-confirmed TB. Our analysis of the plasma cfRNA transcriptome revealed 541 differentially abundant genes, the top 150 of which were used to train 15 machine learning models. The highest performing model led to a 9-gene signature that had a diagnostic accuracy of 89.1% (95% CI: 83.6-93.4%) and an area under the curve of 0.934 (95% CI: 0.8674-1) for microbiologically-confirmed TB. This 9-gene signature exceeds the optimal WHO TPPs for a TB triage test (sensitivity: 96.2% [95% CI: 80.9-100%], specificity: 89.7% [95% CI: 72.4-100%]) and was robust to differences in sample collection, geographic location, and HIV status. Overall, our results demonstrate the utility of plasma cfRNA for the detection of TB and suggest the potential for a point-of-care, gene expression-based assay to aid in early detection of TB.

Copy Number Variation and Allele Ratio Analysis in Candida albicans Using Whole Genome Sequencing Data

Whole genome sequencing of human fungal pathogens has revolutionized the speed and accuracy in which sequence variants that cause antifungal resistance can be identified. Genome rearrangements resulting in copy number variation (CNV) are a significant source of acquired antifungal drug resistance across diverse fungal species. Some CNVs are transient in nature, while other CNVs are stable and well tolerated even in the absence of antifungal drugs. By visualizing whole genome sequencing read depth as a function of genomic location, CNVs and CNV breakpoints (genomic positions where the copy number changes occur relative to the rest of the genome) are rapidly identified. A similar analysis can be used to visualize allele ratio changes that occur across the genomes of heterozygous fungal species, both in the presence and absence of CNVs. This protocol walks through the bioinformatic analysis of CNVs and allele ratios utilizing free, open-source visualization tools. We provide code to use with an example dataset (matched antifungal drug-sensitive and drug-resistant Candida albicans isolates) and notes on how to expand this protocol to other fungal genomes.

Genome Sequence Data Reveal at Least Two Distinct Incursions of the Tropical Race 4 Variant of Fusarium Wilt into South America

The global banana industry is threatened by one of the most devastating diseases: Fusarium wilt of banana. Fusarium wilt of banana is caused by the soilborne fungus Fusarium oxysporum f. sp. cubense (Foc), which almost annihilated the banana production in the late 1950s. A new strain of Foc, known as tropical race 4 (TR4), attacks a wide range of banana varieties, including Cavendish clones, which are the source of 99% of banana exports. In 2019, Foc TR4 was reported in Colombia, and more recently (2021) in Peru. In this study, we sequenced three fungal isolates identified as Foc TR4 from La Guajira (Colombia) and compared them against 19 whole-genome sequences of Foc TR4 publicly available, including four genome sequences recently released from Peru. To understand the genetic relatedness of the Colombian Foc TR4 isolates and those from Peru, we conducted a phylogenetic analysis based on a genome-wide set of single nucleotide polymorphisms (SNPs). Additionally, we compared the genomes of the 22 available Foc TR4 isolates, looking for the presence-absence of gene polymorphisms and genomic regions. Our results reveal that (i) the Colombian and Peruvian isolates are genetically distant, which could be better explained by independent incursions of the pathogen to the continent, and (ii) there is a high correspondence between the genetic relatedness and geographic origin of Foc TR4. The profile of present/absent genes and the distribution of missing genomic regions showed a high correspondence to the clades recovered in the phylogenetic analysis, supporting the results obtained by SNP-based phylogeny.

Detection of Circulating Tumor DNA in Plasma Using Targeted Sequencing

Cell-free DNA (cfDNA) is the degradation product of extracellular DNA. Circulating tumor DNA (ctDNA), as a fraction of cfDNA, comes from tumor cells and contains variations, including mutation, deletion, insertion, rearrangement, copy number variation, and methylation. Therefore, biomarkers identified in ctDNA show promising clinical applications in early diagnosis, recurrence monitoring, and conducting individualized treatment. In this chapter, we introduce experimental workflow and bioinformatic pipeline of targeted sequencing of cfDNA.

Morphological and transcriptomic analysis revealing morphological variations and genetic clues in one Lentinula edodes abnormal browning strain

Strain abnormal browning is a common problem during cultivation of Lentinula edodes. In this study, the L. edodes strain mycelia isolated from Le-WB and cultured on MYG (Le-WP) isolated from an abnormal browning bag was compared with its normal control mycelia isolated from Le-BB and cultured on MYG (Le-BP). The aerial hyphae of Le-WP were white, and the hyphal growth was significantly reduced. Morphological observation of Le-WP under scanning electron microscope (SEM) and transmission electron microscopy (TEM) revealed abnormal organelle structures. Through transcriptomic analysis, more differentially expressed genes (DEGs) were expressed in the metabolic process and catalytic activity in Le-WP than Le-BP. Two Kyoto encyclopedia of genes and genomes (KEGG) pathways named pentose and glucorunate interconversions, and starch and sucrose metabolism were found to be enriched in Le-WP. The gene expression profiles involved in these two pathways were further analyzed and 12 key genes were selected to be verified by quantitative real-time PCR (qRT-PCR), and the results showed that most of these genes were upregulated in Le-WP. Additionally, the content of 1,3-beta-glucan in Le-WP was also significantly higher than in other samples. This research suggests that abnormal strains may be related to the abnormal synthesis of 1,3-beta-glucan, and it needs further research. This research exhibits possible morphological and genetic clues of Le-WP and lays the foundation for understanding the degeneration of L. edodes strains.

RNA-seq analysis of synchronized developing pollen isolated from a single anther

Pollen development, from unicellular microspores to anthesis, is a complex process involving the coordinated specification, differentiation and functions of different cell types. Key to understanding this development is identifying the genes expressed at precise stages of development. However, transcriptomic studies on pollen prior to anthesis are complicated by the inaccessible nature of pollen developing in the anther and the resistant pollen wall. To assist with understanding gene expression during pollen development we have developed a protocol to perform RNA-Seq on pollen isolated from a single anther (SA RNA-Seq). The protocol involves removing pollen from a single anther for analysis and viewing the remaining pollen to determine the developmental stage. The isolated pollen is chemically lysed and mRNA isolated from the lysate using an oligo-dT column before library preparation. Here, we report on the development and testing of our method and the generation of a transcriptome for three stages of pollen development from Arabidopsis (Arabidopsis thaliana) and two stages from male kiwifruit (Actinidia chinensis). This protocol enables the transcriptome of precise developmental stages of pollen to be analyzed, and uses a small number of plants, potentially facilitating studies that require a range of treatments or the analysis of the first generation of transgenic plants.

CanB is a metabolic mediator of antibiotic resistance in Neisseria gonorrhoeae

The evolution of the obligate human pathogen Neisseria gonorrhoeae has been shaped by selective pressures from diverse host niche environments and antibiotics. The varying prevalence of antibiotic resistance across N. gonorrhoeae lineages suggests that underlying metabolic differences may influence the likelihood of acquisition of specific resistance mutations. We hypothesized that the requirement for supplemental CO₂, present in approximately half of isolates, reflects one such example of metabolic variation. Here, using a genome-wide association study and experimental investigations, we show that CO₂ dependence is attributable to a single substitution in a β-carbonic anhydrase, CanB. CanB^19E is necessary and sufficient for growth in the absence of CO₂, and the hypomorphic CanB^19G variant confers CO₂ dependence. Furthermore, ciprofloxacin resistance is correlated with CanB^19G in clinical isolates, and the presence of CanB^19G increases the likelihood of acquisition of ciprofloxacin resistance. Together, our results suggest that metabolic variation has affected the acquisition of fluoroquinolone resistance.

A Multi-Year, Multi-Cultivar Approach to Differential Expression Analysis of High- and Low-Protein Soybean (Glycine max)

Soybean (Glycine max (L.) Merr.) is among the most valuable crops based on its nutritious seed protein and oil. Protein quality, evaluated as the ratio of glycinin (11S) to β-conglycinin (7S), can play a role in food and feed quality. To help uncover the underlying differences between high and low protein soybean varieties, we performed differential expression analysis on high and low total protein soybean varieties and high and low 11S soybean varieties grown in four locations across Eastern and Western Canada over three years (2018-2020). Simultaneously, ten individual differential expression datasets for high vs. low total protein soybeans and ten individual differential expression datasets for high vs. low 11S soybeans were assessed, for a total of 20 datasets. The top 15 most upregulated and the 15 most downregulated genes were extracted from each differential expression dataset and cross-examination was conducted to create shortlists of the most consistently differentially expressed genes. Shortlisted genes were assessed for gene ontology to gain a global appreciation of the commonly differentially expressed genes. Genes with roles in the lipid metabolic pathway and carbohydrate metabolic pathway were differentially expressed in high total protein and high 11S soybeans in comparison to their low total protein and low 11S counterparts. Expression differences were consistent between East and West locations with the exception of one, Glyma.03G054100. These data are important for uncovering the genes and biological pathways responsible for the difference in seed protein between high and low total protein or 11S cultivars.

Aberrant Methylation of Somatostatin Receptor 2 Gene Is Initiated in Aged Gastric Mucosa Infected with Helicobacter pylori and Consequential Gene Silencing Is Associated with Establishment of Inflammatory Microenvironment In Vitro Study

The loss-of-function variants are thought to be associated with inflammation in the stomach. We here aimed to evaluate the extent and role of methylation at the SSTR2 promoter in inflammation and gastric tumor formation. A whole-genome bisulfite sequencing analysis revealed that the SSTR2 promoter was significantly hypermethylated in gastric tumors, dysplasia, and intestinal metaplasia compared to non-tumor tissues from patients with gastric cancer. Using public data, we confirmed SSTR2 promoter methylation in primary gastric tumors and intestinal metaplasia, and even aged gastric mucosae infected with Helicobacter pylori, suggesting that aberrant methylation is initiated in normal gastric mucosa. The loss-of-function of SSTR2 in SNU638 cell-induced cell proliferation in vitro, while stable transfection of SSTR2 in AGS and MKN74 cells inhibited cell proliferation and tumorigenesis in vitro and in vivo. As revealed by a comparison of target genes differentially expressed in these cells with hallmark molecular signatures, inflammation-related pathways were distinctly induced in SSTR2-KO SNU638 cell. By contrast, inflammation-related pathways were inhibited in AGS and MKN74 cells ectopically expressing SSTR2. Collectively, we propose that SSTR2 silencing upon promoter methylation is initiated in aged gastric mucosae infected with H. pylori and promotes the establishment of an inflammatory microenvironment via the intrinsic pathway. These findings provide novel insights into the initiation of gastric carcinogenesis.

Pantothenate and L-Carnitine Supplementation Improves Pathological Alterations in Cellular Models of KAT6A Syndrome

Mutations in several genes involved in the epigenetic regulation of gene expression have been considered risk alterations to different intellectual disability (ID) syndromes associated with features of autism spectrum disorder (ASD). Among them are the pathogenic variants of the lysine-acetyltransferase 6A (KAT6A) gene, which causes KAT6A syndrome. The KAT6A enzyme participates in a wide range of critical cellular functions, such as chromatin remodeling, gene expression, protein synthesis, cell metabolism, and replication. In this manuscript, we examined the pathophysiological alterations in fibroblasts derived from three patients harboring KAT6A mutations. We addressed survival in a stress medium, histone acetylation, protein expression patterns, and transcriptome analysis, as well as cell bioenergetics. In addition, we evaluated the therapeutic effectiveness of epigenetic modulators and mitochondrial boosting agents, such as pantothenate and L-carnitine, in correcting the mutant phenotype. Pantothenate and L-carnitine treatment increased histone acetylation and partially corrected protein and transcriptomic expression patterns in mutant KAT6A cells. Furthermore, the cell bioenergetics of mutant cells was significantly improved. Our results suggest that pantothenate and L-carnitine can significantly improve the mutant phenotype in cellular models of KAT6A syndrome.

Late Pleistocene-dated divergence between South Hemisphere populations of the non-conventional yeast L. cidri

Most organisms belonging to the Saccharomycotina subphylum have high genetic diversity and a vast repertoire of metabolisms and lifestyles. Lachancea cidri is an ideal yeast model for exploring the interplay between genetics, ecological function and evolution. Lachancea cidri diverged from the Saccharomyces lineage before the whole-genome duplication and is distributed across the South Hemisphere, displaying an important ecological success. We applied phylogenomics to investigate the genetic variation of L. cidri isolates obtained from Australia and South America. Our approach revealed the presence of two main lineages according to their geographic distribution (Aus and SoAm). Estimation of the divergence time suggests that SoAm and Aus lineages diverged near the last glacial maximum event during the Pleistocene (64-8 KYA). Interestingly, we found that the French reference strain is closely related to the Australian strains, with a recent divergence (405-51 YA), likely associated to human movements. Additionally, we identified different lineages within the South American population, revealing that Patagonia contains a similar genetic diversity comparable to that of other lineages in S. cerevisiae. These findings support the idea of a Pleistocene-dated divergence between South Hemisphere lineages, where the Nothofagus and Araucaria ecological niches likely favoured the extensive distribution of L. cidri in Patagonia.

Ataxin-2, Twenty-four, and Dicer-2 are components of a noncanonical cytoplasmic polyadenylation complex

Identification of components of a noncanonical cytoplasmic polyadenylation machinery in Drosophila expands the diversity of RNA-binding proteins involved in poly(A) tail length control.

Cytoplasmic polyadenylation is a mechanism to promote mRNA translation in a wide variety of biological contexts. A canonical complex centered around the conserved RNA-binding protein family CPEB has been shown to be responsible for this process. We have previously reported evidence for an alternative noncanonical, CPEB-independent complex in Drosophila, of which the RNA-interference factor Dicer-2 is a component. Here, we investigate Dicer-2 mRNA targets and protein cofactors in cytoplasmic polyadenylation. Using RIP-Seq analysis, we identify hundreds of potential Dicer-2 target transcripts, ∼60% of which were previously found as targets of the cytoplasmic poly(A) polymerase Wispy, suggesting widespread roles of Dicer-2 in cytoplasmic polyadenylation. Large-scale immunoprecipitation revealed Ataxin-2 and Twenty-four among the high-confidence interactors of Dicer-2. Complex analyses indicated that both factors form an RNA-independent complex with Dicer-2 and mediate interactions of Dicer-2 with Wispy. Functional poly(A)-test analyses showed that Twenty-four and Ataxin-2 are required for cytoplasmic polyadenylation of a subset of Dicer-2 targets. Our results reveal components of a novel cytoplasmic polyadenylation complex that operates during Drosophila early embryogenesis.

Campylobacter majalis sp. nov. and Campylobacter suis sp. nov., novel Campylobacter species isolated from porcine gastrointestinal mucosa

Strains LMG 7974^T and LMG 8286^T represent single, novel Campylobacter lineages with Campylobacter pinnipediorum and Campylobacter mucosalis as nearest phylogenomic neighbours, respectively. The results of average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) analyses of LMG 7974^T, LMG 8286^T and type strains of species of the genus Campylobacter confirmed that these strains represent novel species of the genus Campylobacter. The 16S rRNA gene sequences of both strains showed highest identity towards C. mucosalis (97.84 and 98.74 %, respectively). Strains LMG 7974^T and LMG 8286^T shared 72.5 and 73.7% ANI, respectively, with their nearest phylogenomic neighbours and less than 21 % dDDH. The draft genome sizes of LMG 7974^T and LMG 8286^T are 1 945429 bp and 1 708214 bp in length with percentage DNA G+C contents of 33.8 and 37.2 %, respectively. Anomalous biochemical characteristics and low MALDI-TOF mass spectrometry log scores supported their designation as representing novel species of the genus Campylobacte. We therefore propose to classify strain LMG 7974^T (=CCUG 20705^T) as the type strain of the novel species Campylobacter majalis sp. nov. and strain LMG 8286^T (=CCUG 24193^T, NCTC 11879^T) as the type strain of the novel species Campylobacter suis sp. nov.

Genome stability during serial subculturing in hyperepidemic multidrug-resistant Klebsiella pneumoniae and Escherichia coli

OBJECTIVES

Core-genome single nucleotide variant (cgSNV) analysis represents a powerful tool for epidemiological investigations of multidrug-resistant (MDR) bacteria. However, cgSNV thresholds to confirm whether isolates are the same clone are not formally defined.

METHODS

We implemented hybrid whole-genome sequencing to study the genomic changes of four MDR isolates belonging to hyperepidemic sequence types (STs) during 20 propagation steps (T20) on MacConkey and CHROMID^(R) ESBL plates. The following strains were analyzed: Klebsiella pneumoniae AE-2247421 (OXA-48/NDM-1-producing, ST101), K. pneumoniae MCL-2017-2 (CTX-M-15-producing, ST307), Escherichia coli Ec-042 (OXA-181-producing, ST410), and E. coli Ec-050 (NDM-5-producing, ST167). The genome assembly at T5 and T20 was compared to that at time point zero (T0) and to two reference genomes.

RESULTS

At T20, AE-2247421 lost the IncL bla_OXA-48-carrying plasmid when grown on CHROMID^(R) ESBL plates, while a large fragment encompassing bla_NDM-1 was lost from its IncC plasmid when grown on both plates. In contrast, no structural changes were noted for the other three strains. Regarding the cgSNVs, the following results were obtained at T5 and T20 (ranges considering the different agar plates and reference genomes): AE-2247421 (1-8 and 2-12 cgSNVs), MCL-2017-2 (both 1-2 cgSNVs), Ec-042 (both 0 cgSNVs), and Ec-050 (0-6 and 0-9 cgSNVs).

CONCLUSION

We showed that structural changes and accumulation of cgSNVs can occur in few propagation steps under laboratory conditions. These changes might also arise in the clinical context in a short time, especially under antibiotics treatment. This phenomenon should be carefully considered because it might affect the final interpretation of epidemiological genomic analyses.

Kazak mitochondrial genomes provide insights into the human population history of Central Eurasia

As a historical nomadic group in Central Asia, Kazaks have mainly inhabited the steppe zone from the Altay Mountains in the East to the Caspian Sea in the West. Fine scale characterization of the genetic profile and population structure of Kazaks would be invaluable for understanding their population history and modeling prehistoric human expansions across the Eurasian steppes. With this mind, we characterized the maternal lineages of 200 Kazaks from Jetisuu at mitochondrial genome level. Our results reveal that Jetisuu Kazaks have unique mtDNA haplotypes including those belonging to the basal branches of both West Eurasian (R0, H, HV) and East Eurasian (A, B, C, D) lineages. The great diversity observed in their maternal lineages may reflect pivotal geographic location of Kazaks in Eurasia and implies a complex history for this population. Comparative analyses of mitochondrial genomes of human populations in Central Eurasia reveal a common maternal genetic ancestry for Turko-Mongolian speakers and their expansion being responsible for the presence of East Eurasian maternal lineages in Central Eurasia. Our analyses further indicate maternal genetic affinity between the Sherpas from the Tibetan Plateau with the Turko-Mongolian speakers.

ONT-Based Alternative Assemblies Impact on the Annotations of Unique versus Repetitive Features in the Genome of a Romanian Strain of Drosophila melanogaster

To date, different strategies of whole-genome sequencing (WGS) have been developed in order to understand the genome structure and functions. However, the analysis of genomic sequences obtained from natural populations is challenging and the biological interpretation of sequencing data remains the main issue. The MinION device developed by Oxford Nanopore Technologies (ONT) is able to generate long reads with minimal costs and time requirements. These valuable assets qualify it as a suitable method for performing WGS, especially in small laboratories. The long reads resulted using this sequencing approach can cover large structural variants and repetitive sequences commonly present in the genomes of eukaryotes. Using MinION, we performed two WGS assessments of a Romanian local strain of Drosophila melanogaster, referred to as Horezu_LaPeri (Horezu). In total, 1,317,857 reads with a size of 8.9 gigabytes (Gb) were generated. Canu and Flye de novo assembly tools were employed to obtain four distinct assemblies with both unfiltered and filtered reads, achieving maximum reference genome coverages of 94.8% (Canu) and 91.4% (Flye). In order to test the quality of these assemblies, we performed a two-step evaluation. Firstly, we considered the BUSCO scores and inquired for a supplemental set of genes using BLAST. Subsequently, we appraised the total content of natural transposons (NTs) relative to the reference genome (ISO1 strain) and mapped the mdg1 retroelement as a resolution assayer. Our results reveal that filtered data provide only slightly enhanced results when considering genes identification, but the use of unfiltered data had a consistent positive impact on the global evaluation of the NTs content. Our comparative studies also revealed differences between Flye and Canu assemblies regarding the annotation of unique versus repetitive genomic features. In our hands, Flye proved to be moderately better for gene identification, while Canu clearly outperformed Flye for NTs analysis. Data concerning the NTs content were compared to those obtained with ONT for the D. melanogaster ISO1 strain, revealing that our strategy conducted to better results. Additionally, the parameters of our ONT reads and assemblies are similar to those reported for ONT experiments performed on various model organisms, revealing that our assembly data are appropriate for a proficient annotation of the Horezu genome.

Norepinephrine transporter defects lead to sympathetic hyperactivity in Familial Dysautonomia models

Familial dysautonomia (FD), a rare neurodevelopmental and neurodegenerative disorder affects the sympathetic and sensory nervous system. Although almost all patients harbor a mutation in ELP1, it remains unresolved exactly how function of sympathetic neurons (symNs) is affected; knowledge critical for understanding debilitating disease hallmarks, including cardiovascular instability or dysautonomic crises, that result from dysregulated sympathetic activity. Here, we employ the human pluripotent stem cell (hPSC) system to understand symN disease mechanisms and test candidate drugs. FD symNs are intrinsically hyperactive in vitro, in cardiomyocyte co-cultures, and in animal models. We report reduced norepinephrine transporter expression, decreased intracellular norepinephrine (NE), decreased NE re-uptake, and excessive extracellular NE in FD symNs. SymN hyperactivity is not a direct ELP1 mutation result, but may connect to NET via RAB proteins. We found that candidate drugs lowered hyperactivity independent of ELP1 modulation. Our findings may have implications for other symN disorders and may allow future drug testing and discovery.

Complete Genome Sequence of the Type Strain Aliivibrio fischeri DSM 507

The complete genome sequence of the type strain Aliivibrio fischeri comb. nov. DSM 507 is presented. The genome consists of two circular chromosomes comprising 2,970,859 bp (chromosome 1) and 1,532,347 bp (chromosome 2) with GC contents of 38.94 and 37.13%, respectively.

De novo genome assembly and analysis of Zalaria sp. Him3, a novel fructooligosaccharides producing yeast

Background

Zalaria sp. Him3 was reported as a novel fructooligosaccharides (FOS) producing yeast. However, Zalaria spp. have not been widely known and have been erroneously classified as a different black yeast, Aureobasidium pullulans. In this study, de novo genome assembly and analysis of Zalaria sp. Him3 was demonstrated to confirm the existence of a potential enzyme that facilitates FOS production and to compare with the genome of A. pullulans.

Results

The genome of Zalaria sp. Him3 was analyzed; the total read bases and total number of reads were 6.38 Gbp and 42,452,134 reads, respectively. The assembled genome sequence was calculated to be 22.38 Mbp, with 207 contigs, N50 of 885,387, L50 of 10, GC content of 53.8%, and 7,496 genes. g2419, g3120, and g3700 among the predicted genes were annotated as cellulase, xylanase, and β-fructofuranosidase (FFase), respectively. When the read sequences were mapped to A. pullulans EXF-150 genome as a reference, a small amount of reads (3.89%) corresponded to the reference genome. Phylogenetic tree analysis, which was based on the conserved sequence set consisting of 2,362 orthologs in the genome, indicated genetic differences between Zalaria sp. Him3 and Aureobasidium spp.

Conclusion

The differences between Zalaria and Aureobasidium spp. were evident at the genome level. g3700 identified in the Zalaria sp. Him3 likely does not encode a highly transfructosyl FFase because the motif sequences were unlike those in other FFases involved in FOS production. Therefore, strain Him3 may produce another FFase. Furthermore, several genes with promising functions were identified and might elicit further interest in Zalaria yeast.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-022-01094-2.

Transcriptome analyses unveiled differential regulation of AGO and DCL genes by pepino mosaic virus strains

Pepino mosaic virus (PepMV) is a single-stranded (ss), positive-sense (+) RNA potexvirus that affects tomato crops worldwide. We have described an in planta antagonistic interaction between PepMV isolates of two strains in which the EU isolate represses the accumulation of the CH2 isolate during mixed infections. Reports describing transcriptomic responses to mixed infections are scant. We carried out transcriptomic analyses of tomato plants singly and mixed-infected with two PepMV isolates of both strains. Comparison of the transcriptomes of singly infected plants showed that deeper transcriptomic alterations occurred at early infection times, and also that each of the viral strains modulated the host transcriptome differentially. Mixed infections caused transcriptomic alterations similar to those for the sum of single infections at early infection times, but clearly differing at later times postinfection. We next tested the hypothesis that PepMV-EU, in either single or mixed infections, deregulates host gene expression differentially so that virus accumulation of both strains gets repressed. That seemed to be the case for the genes AGO1a, DCL2d, AGO2a, and DCL2b, which are involved in the antiviral silencing pathway and were upregulated by PepMV-EU but not by PepMV-CH2 at early times postinfection. The pattern of AGO2a expression was validated by reverse transcription-quantitative PCR in tomato and Nicotiana benthamiana plants. Using an N. benthamiana ago2 mutant line, we showed that AGO2 indeed plays an important role in the antiviral defence against PepMV, but it is not the primary determinant of the outcome of the antagonistic interaction between the two PepMV strains.

Selective sweeps on different pigmentation genes mediate convergent evolution of island melanism in two incipient bird species

Insular organisms often evolve predictable phenotypes, like flightlessness, extreme body sizes, or increased melanin deposition. The evolutionary forces and molecular targets mediating these patterns remain mostly unknown. Here we study the Chestnut-bellied Monarch (Monarcha castaneiventris) from the Solomon Islands, a complex of closely related subspecies in the early stages of speciation. On the large island of Makira M. c. megarhynchus has a chestnut belly, whereas on the small satellite islands of Ugi, and Santa Ana and Santa Catalina (SA/SC) M. c. ugiensis is entirely iridescent blue-black (i.e., melanic). Melanism has likely evolved twice, as the Ugi and SA/SC populations were established independently. To investigate the genetic basis of melanism on each island we generated whole genome sequence data from all three populations. Non-synonymous mutations at the MC1R pigmentation gene are associated with melanism on SA/SC, while ASIP, an antagonistic ligand of MC1R, is associated with melanism on Ugi. Both genes show evidence of selective sweeps in traditional summary statistics and statistics derived from the ancestral recombination graph (ARG). Using the ARG in combination with machine learning, we inferred selection strength, timing of onset and allele frequency trajectories. MC1R shows evidence of a recent, strong, soft selective sweep. The region including ASIP shows more complex signatures; however, we find evidence for sweeps in mutations near ASIP, which are comparatively older than those on MC1R and have been under relatively strong selection. Overall, our study shows convergent melanism results from selective sweeps at independent molecular targets, evolving in taxa where coloration likely mediates reproductive isolation with the neighboring chestnut-bellied subspecies.

Chestnut-bellied Monarchs (Monarcha castaneiventris ugiensis) from two archipelagos in the Solomon Islands have evolved entirely black plumage from a chestnut ancestor (Monarcha castaneiventris megarhynchus), a phenomenon known as island melanism. We obtain and analyze whole genome sequences using traditional summary statistics and new methods that combine inference of the ancestral recombination graph with machine learning. We find multiple lines of evidence for independent selective sweeps on the MC1R and ASIP genes, a receptor/ligand pair which regulates the production of melanin. Melanism on each archipelago is mediated by mutations in one of these two genes. Mutations in and around MC1R underwent a recent soft sweep experiencing strong selection on the islands of Santa Ana and Santa Catalina, whereas selection was also strong but comparatively older for ASIP on the island of Ugi. We show how melanism originated under positive selection on independent molecular targets, evolving convergently in taxa where coloration mediates reproductive isolation.

Overcoming variant mutation-related impacts on viral sequencing and detection methodologies

Prompt and accurate pathogen identification, by diagnostics and sequencing, is an effective tool for tracking and potentially curbing pathogen spread. Targeted detection and amplification of viral genomes depends on annealing complementary oligonucleotides to genomic DNA or cDNA. However, genomic mutations that occur during viral evolution may perturb annealing, which can result in incomplete sequence coverage of the genome and/or false negative diagnostic test results. Herein, we demonstrate how to assess, test, and optimize sequencing and detection methodologies to attenuate the negative impact of mutations on genome targeting efficiency. This evaluation was conducted using in vitro-transcribed (IVT) RNA as well as RNA extracted from clinical SARS-CoV-2 variant samples, including the heavily mutated Omicron variant. Using SARS-CoV-2 as a current example, these results demonstrate how to maintain reliable targeted pathogen sequencing and how to evaluate detection methodologies as new variants emerge.

Two type I topoisomerases maintain DNA topology in human mitochondria

Genetic processes require the activity of multiple topoisomerases, essential enzymes that remove topological tension and intermolecular linkages in DNA. We have investigated the subcellular localisation and activity of the six human topoisomerases with a view to understanding the topological maintenance of human mitochondrial DNA. Our results indicate that mitochondria contain two topoisomerases, TOP1MT and TOP3A. Using molecular, genomic and biochemical methods we find that both proteins contribute to mtDNA replication, in addition to the decatenation role of TOP3A, and that TOP1MT is stimulated by mtSSB. Loss of TOP3A or TOP1MT also dysregulates mitochondrial gene expression, and both proteins promote transcription elongation in vitro. We find no evidence for TOP2 localisation to mitochondria, and TOP2B knockout does not affect mtDNA maintenance or expression. Our results suggest a division of labour between TOP3A and TOP1MT in mtDNA topology control that is required for the proper maintenance and expression of human mtDNA.

A Screening of Epigenetic Therapeutic Targets for Non-Small Cell Lung Cancer Reveals PADI4 and KDM6B as Promising Candidates

Despite advances in diagnostic and therapeutic approaches for lung cancer, new therapies targeting metastasis by the specific regulation of cancer genes are needed. In this study, we screened a small library of epigenetic inhibitors in non-small-cell lung cancer (NSCLC) cell lines and evaluated 38 epigenetic targets for their potential role in metastatic NSCLC. The potential candidates were ranked by a streamlined approach using in silico and in vitro experiments based on publicly available databases and evaluated by real-time qPCR target gene expression, cell viability and invasion assays, and transcriptomic analysis. The survival rate of patients with lung adenocarcinoma is inversely correlated with the gene expression of eight epigenetic targets, and a systematic review of the literature confirmed that four of them have already been identified as targets for the treatment of NSCLC. Using nontoxic doses of the remaining inhibitors, KDM6B and PADI4 were identified as potential targets affecting the invasion and migration of metastatic lung cancer cell lines. Transcriptomic analysis of KDM6B and PADI4 treated cells showed altered expression of important genes related to the metastatic process. In conclusion, we showed that KDM6B and PADI4 are promising targets for inhibiting the metastasis of lung adenocarcinoma cancer cells.

Editing melon eIF4E associates with virus resistance and male sterility

The cap-binding protein eIF4E, through its interaction with eIF4G, constitutes the core of the eIF4F complex, which plays a key role in the circularization of mRNAs and their subsequent cap-dependent translation. In addition to its fundamental role in mRNA translation initiation, other functions have been described or suggested for eIF4E, including acting as a proviral factor and participating in sexual development. We used CRISPR/Cas9 genome editing to generate melon eif4e knockout mutant lines. Editing worked efficiently in melon, as we obtained transformed plants with a single-nucleotide deletion in homozygosis in the first eIF4E exon already in a T0 generation. Edited and non-transgenic plants of a segregating F2 generation were inoculated with Moroccan watermelon mosaic virus (MWMV); homozygous mutant plants showed virus resistance, while heterozygous and non-mutant plants were infected, in agreement with our previous results with plants silenced in eIF4E. Interestingly, all homozygous edited plants of the T0 and F2 generations showed a male sterility phenotype, while crossing with wild-type plants restored fertility, displaying a perfect correlation between the segregation of the male sterility phenotype and the segregation of the eif4e mutation. Morphological comparative analysis of melon male flowers along consecutive developmental stages showed postmeiotic abnormal development for both microsporocytes and tapetum, with clear differences in the timing of tapetum degradation in the mutant versus wild-type. An RNA-Seq analysis identified critical genes in pollen development that were down-regulated in flowers of eif4e/eif4e plants, and suggested that eIF4E-specific mRNA translation initiation is a limiting factor for male gametes formation in melon.

Rare genetic variant burden in DPYD predicts severe fluoropyrimidine-related toxicity risk

Preemptive targeted pharmacogenetic testing of candidate variations in DPYD is currently being used to limit toxicity associated with fluoropyrimidines. The use of innovative next generation sequencing (NGS) approaches could unveil additional rare (minor allele frequency <1%) genetic risk variants. However, their predictive value and management in clinical practice are still controversial, at least partly due to the challenges associated with functional analyses of rare variants. The aim of this study was to define the predictive power of rare DPYD variants burden on the risk of severe fluoropyrimidine-related toxicity.

The DPYD coding sequence and untranslated regions were analyzed by NGS in 120 patients developing grade 3–5 (NCI-CTC vs3.0) fluoropyrimidine-related toxicity and 104 matched controls (no-toxicity). The functional impact of rare variants was assessed using two different in silico predictive tools (i.e., Predict2SNP and ADME Prediction Framework) and structural modeling. Plasma concentrations of uracil (U) and dihydrouracil (UH2) were quantified in carriers of the novel variants.

Here, we demonstrate that the burden of rare variants was significantly higher in patients with toxicity compared to controls (p = 0.007, Mann-Whitney test). Carriers of at least one rare missense DPYD variant had a 16-fold increased risk in the first cycle and an 11-fold increased risk during the entire course of chemotherapy of developing a severe adverse event compared to controls (p = 0.013 and p = 0.0250, respectively by multinomial regression model). Quantification of plasmatic U/UH2 metabolites and in silico visualization of the encoded protein were consistent with the predicted functional effect for the novel variations.

Analysis and consideration of rare variants by DPYD-sequencing could improve prevention of severe toxicity of fluoropyrimidines and improve patients’ quality of life.

•

DPYD genotype-guided dosing reduces fluoropyrimidine (FP) toxicity risk.

•

Rare DPYD variants associate with severe FP toxicities.

•

Carriers of rare DPYD variants have 11-fold increased risk of toxicity.

•

DPYD sequencing and in silico functional prediction could prevent FP toxicity events.

Uncovering anthocyanin diversity in potato landraces (Solanum tuberosum L. Phureja) using RNA-seq

Potato (Solanum tuberosum L.) is the third largest source of antioxidants in the human diet, after maize and tomato. Potato landraces have particularly diverse contents of antioxidant compounds such as anthocyanins. We used this diversity to study the evolutionary and genetic basis of anthocyanin pigmentation. Specifically, we analyzed the transcriptomes and anthocyanin content of tubers from 37 landraces with different colorations. We conducted analyses of differential expression between potatoes with different colorations and used weighted correlation network analysis to identify genes whose expression is correlated to anthocyanin content across landraces. A very significant fraction of the genes identified in these two analyses had annotations related to the flavonoid-anthocyanin biosynthetic pathway, including 18 enzymes and 5 transcription factors. Importantly, the causal genes at the D, P and R loci governing anthocyanin accumulation in potato cultivars also showed correlations to anthocyanin production in the landraces studied here. Furthermore, we found that 60% of the genes identified in our study were located within anthocyanin QTLs. Finally, we identified new candidate enzymes and transcription factors that could have driven the diversification of anthocyanins. Our results indicate that many anthocyanins biosynthetic genes were manipulated in ancestral potato breeding and can be used in future breeding programs.

A catalog of the genetic causes of hereditary angioedema in the Canary Islands (Spain)

Hereditary angioedema (HAE) is a rare disease where known causes involve C1 inhibitor dysfunction or dysregulation of the kinin cascade. The updated HAE management guidelines recommend performing genetic tests to reach a precise diagnosis. Unfortunately, genetic tests are still uncommon in the diagnosis routine. Here, we characterized for the first time the genetic causes of HAE in affected families from the Canary Islands (Spain). Whole-exome sequencing data was obtained from 41 affected patients and unaffected relatives from 29 unrelated families identified in the archipelago. The Hereditary Angioedema Database Annotation (HADA) tool was used for pathogenicity classification and causal variant prioritization among the genes known to cause HAE. Manual reclassification of prioritized variants was used in those families lacking known causal variants. We detected a total of eight different variants causing HAE in this patient series, affecting essentially SERPING1 and F12 genes, one of them being a novel SERPING1 variant (c.686-12A>G) with a predicted splicing effect which was reclassified as likely pathogenic in one family. Altogether, the diagnostic yield by assessing previously reported causal genes and considering variant reclassifications according to the American College of Medical Genetics guidelines reached 66.7% (95% Confidence Interval [CI]: 30.1-91.0) in families with more than one affected member and 10.0% (95% CI: 1.8-33.1) among cases without family information for the disease. Despite the genetic causes of many patients remain to be identified, our results reinforce the need of genetic tests as first-tier diagnostic tool in this disease, as recommended by the international WAO/EAACI guidelines for the management of HAE.

Epigenetic machinery is functionally conserved in cephalopods

BACKGROUND

Epigenetic regulatory mechanisms are divergent across the animal kingdom, yet these mechanisms are not well studied in non-model organisms. Unique features of cephalopods make them attractive for investigating behavioral, sensory, developmental, and regenerative processes, and recent studies have elucidated novel features of genome organization and gene and transposon regulation in these animals. However, it is not known how epigenetics regulates these interesting cephalopod features. We combined bioinformatic and molecular analysis of Octopus bimaculoides to investigate the presence and pattern of DNA methylation and examined the presence of DNA methylation and 3 histone post-translational modifications across tissues of three cephalopod species.

RESULTS

We report a dynamic expression profile of the genes encoding conserved epigenetic regulators, including DNA methylation maintenance factors in octopus tissues. Levels of 5-methyl-cytosine in multiple tissues of octopus, squid, and bobtail squid were lower compared to vertebrates. Whole genome bisulfite sequencing of two regions of the brain and reduced representation bisulfite sequencing from a hatchling of O. bimaculoides revealed that less than 10% of CpGs are methylated in all samples, with a distinct pattern of 5-methyl-cytosine genome distribution characterized by enrichment in the bodies of a subset of 14,000 genes and absence from transposons. Hypermethylated genes have distinct functions and, strikingly, many showed similar expression levels across tissues while hypomethylated genes were silenced or expressed at low levels. Histone marks H3K27me3, H3K9me3, and H3K4me3 were detected at different levels across tissues of all species.

CONCLUSIONS

Our results show that the DNA methylation and histone modification epigenetic machinery is conserved in cephalopods, and that, in octopus, 5-methyl-cytosine does not decorate transposable elements, but is enriched on the gene bodies of highly expressed genes and could cooperate with the histone code to regulate tissue-specific gene expression.

Myxoma virus lacking the host range determinant M062 stimulates cGAS-dependent type 1 interferon response and unique transcriptomic changes in human monocytes/macrophages

The evolutionarily successful poxviruses possess effective and diverse strategies to circumvent or overcome host defense mechanisms. Poxviruses encode many immunoregulatory proteins to evade host immunity to establish a productive infection and have unique means of inhibiting DNA sensing-dependent type 1 interferon (IFN-I) responses, a necessity given their dsDNA genome and exclusively cytoplasmic life cycle. We found that the key DNA sensing inhibition by poxvirus infection was dominant during the early stage of poxvirus infection before DNA replication. In an effort to identify the poxvirus gene products which subdue the antiviral proinflammatory responses (e.g., IFN-I response), we investigated the function of one early gene that is the known host range determinant from the highly conserved poxvirus host range C7L superfamily, myxoma virus (MYXV) M062. Host range factors are unique features of poxviruses that determine the species and cell type tropism. Almost all sequenced mammalian poxviruses retain at least one homologue of the poxvirus host range C7L superfamily. In MYXV, a rabbit-specific poxvirus, the dominant and broad-spectrum host range determinant of the C7L superfamily is the M062R gene. The M062R gene product is essential for MYXV infection in almost all cells tested from different mammalian species and specifically inhibits the function of host Sterile αMotif Domain-containing 9 (SAMD9), as M062R-null (ΔM062R) MYXV causes abortive infection in a SAMD9-dependent manner. In this study we investigated the immunostimulatory property of the ΔM062R. We found that the replication-defective ΔM062R activated host DNA sensing pathway during infection in a cGAS-dependent fashion and that knocking down SAMD9 expression attenuated proinflammatory responses. Moreover, transcriptomic analyses showed a unique feature of the host gene expression landscape that is different from the dsDNA alone-stimulated inflammatory state. This study establishes a link between the anti-neoplastic function of SAMD9 and the regulation of innate immune responses.

Poxviruses encode a group of genes called host range determinants to maintain or expand their host tropism. The mechanism by which many viral host range factors function remains elusive. Some host range factors possess immunoregulatory functions responsible for evading or subduing host immune defense mechanisms. Most known immunoregulatory proteins encoded by poxviruses are dispensable for viral replication in vitro. The uniqueness of MYXV M062R is that it is essential for viral infection in vitro and belongs to one of the most conserved poxvirus host range families, the C7L superfamily. There is one known host target of the MYXV M062 protein, SAMD9. SAMD9 is constitutively expressed in mammalian cells and exclusively present in the cytoplasm with an anti-neoplastic function. Humans with deleterious mutations in SAMD9 present disease that ranges from lethality at a young age to a predisposition to myelodysplastic syndromes (MDS) that often require bone marrow transplantation. More importantly, SAMD9 serves as an important antiviral intrinsic molecule to many viruses. The cellular function of SAMD9 remains unclear mostly due to the difficulty of studying this protein, i.e., its large size, long half-life, and its constitutive expression in most cells. In this study we used M062R-null MYXV as a tool to study SAMD9 function and report a functional link between SAMD9 and the regulation of the proinflammatory responses triggered by cGAS-dependent DNA sensing.

Molecular convergence between Down syndrome and fragile X syndrome identified using human pluripotent stem cell models

Down syndrome (DS), driven by an extra copy of chromosome 21 (HSA21), and fragile X syndrome (FXS), driven by loss of the RNA-binding protein FMRP, are two common genetic causes of intellectual disability and autism. Based upon the number of DS-implicated transcripts bound by FMRP, we hypothesize that DS and FXS may share underlying mechanisms. Comparing DS and FXS human pluripotent stem cell (hPSC) and glutamatergic neuron models, we identify increased protein expression of select targets and overlapping transcriptional perturbations. Moreover, acute upregulation of endogenous FMRP in DS patient cells using CRISPRa is sufficient to significantly reduce expression levels of candidate proteins and reverse 40% of global transcriptional perturbations. These results pinpoint specific molecular perturbations shared between DS and FXS that can be leveraged as a strategy for target prioritization; they also provide evidence for the functional relevance of previous associations between FMRP targets and disease-implicated genes.

3D microenvironment attenuates simulated microgravity-mediated changes in T cell transcriptome

Human space travel and exploration are of interest to both the industrial and scientific community. However, there are many adverse effects of spaceflight on human physiology. In particular, there is a lack of understanding of the extent to which microgravity affects the immune system. T cells, key players of the adaptive immune system and long-term immunity, are present not only in blood circulation but also reside within the tissue. As of yet, studies investigating the effects of microgravity on T cells are limited to peripheral blood or traditional 2D cell culture that recapitulates circulating blood. To better mimic interstitial tissue, 3D cell culture has been well established for physiologically and pathologically relevant models. In this work, we utilize 2D cell culture and 3D collagen matrices to gain an understanding of how simulated microgravity, using a random positioning machine, affects both circulating and tissue-resident T cells. T cells were studied in both resting and activated stages. We found that 3D cell culture attenuates the effects of simulated microgravity on the T cells transcriptome and nuclear irregularities compared to 2D cell culture. Interestingly, simulated microgravity appears to have less effect on activated T cells compared to those in the resting stage. Overall, our work provides novel insights into the effects of simulated microgravity on circulating and tissue-resident T cells which could provide benefits for the health of space travellers.

Whole-exome Sequencing of Nigerian Prostate Tumors from the Prostate Cancer Transatlantic Consortium (CaPTC) Reveals DNA Repair Genes Associated with African Ancestry

In this study, we used whole-exome sequencing of a cohort of 45 advanced-stage, treatment-naïve Nigerian (NG) primary prostate cancer tumors and 11 unmatched nontumor tissues to compare genomic mutations with African American (AA) and European American (EA) The Cancer Genome Atlas (TCGA) prostate cancer. NG samples were collected from six sites in central and southwest Nigeria. After whole-exome sequencing, samples were processed using GATK best practices. BRCA1 (100%), BARD1 (45%), BRCA2 (27%), and PMS2(18%) had germline alterations in at least two NG nontumor samples. Across 111 germline variants, the AA cohort reflected a pattern [BRCA1 (68%), BARD1 (34%), BRCA2 (28%), and PMS2 (16%)] similar to NG samples. Of the most frequently mutated genes, BRCA1 showed a statistically (P ≤ 0.05) higher germline mutation frequency in men of African ancestry (MAA) and increasing variant frequency with increased African ancestry. Disaggregating gene-level mutation frequencies by variants revealed both ancestry-linked and NG-specific germline variant patterns. Driven by rs799917 (T>C), BRCA1 showed an increasing mutation frequency as African ancestry increased. BRCA2_rs11571831 was present only in MAA, and BRCA2_rs766173 was elevated in NG men. A total of 133 somatic variants were present in 26 prostate cancer-associated genes within the NG tumor cohort. BRCA2 (27%), APC (20%), ATM (20%), BRCA1 (13%), DNAJC6 (13%), EGFR (13%), MAD1L1 (13%), MLH1 (11%), and PMS2 (11%) showed mutation frequencies >10%. Compared with TCGA cohorts, NG tumors showed statistically significant elevated frequencies of BRCA2, APC, and BRCA1. The NG cohort variant pattern shared similarities (cosign similarities ≥0.734) with Catalogue of Somatic Mutations in Cancer signatures 5 and 6, and mutated genes showed significant (q < 0.001) gene ontology (GO) and functional enrichment in mismatch repair and non-homologous repair deficiency pathways. Here, we showed that mutations in DNA damage response genes were higher in NG prostate cancer samples and that a portion of those mutations correlate with African ancestry. Moreover, we identified variants of unknown significance that may contribute to population-specific routes of tumorigenesis and treatment. These results present the most comprehensive characterization of the NG prostate cancer exome to date and highlight the need to increase diversity of study populations.

Significance

MAA have higher rates of prostate cancer incidence and mortality, however, are severely underrepresented in genomic studies. This is the first study utilizing whole-exome sequencing in NG men to identify West African ancestry-linked variant patterns that impact DNA damage repair pathways.

Genomic characteristics of two most widely used BCG vaccine strains: Danish 1331 and Pasteur 1173P2

Background

Bacillus Calmette–Guérin (BCG) refers to a group of vaccine strains with unique genetic characteristics. BCG is the only available vaccine for preventing tuberculosis (TB). Genetic and biochemical variations among the BCG vaccine strains have been considered as one of the significant parameters affecting the variable protective efficacy of the vaccine against pulmonary tuberculosis. To track genetic variations, here two vaccine strains (Danish 1331 and Pasteur 1173P2) popularly used according to the BCG World Atlas were subjected to a comparative analysis against the Mycobacterium tuberculosis H37Rv, Mycobacterium bovis AF2122/97, and Mycobacterium tuberculosis variant bovis BCG str. Pasteur 1173P2 reference genomes. Besides, the presence or absence of the experimentally verified human T cell epitopes was examined.

Results

Only two variants were identified in BCG Danish 1331 that have not been reported previously in any BCG strains with the complete submitted genome yet. Furthermore, we identified a DU1-like 14,577 bp region in BCG Danish 1331; The duplication which was previously seemed to be exclusive to the BCG Pasteur. We also found that 35% of the T cell epitopes are absent from both strains, and epitope sequences are more conserved than the rest of the genome.

Conclusions

We provided a comprehensive catalog of single nucleotide polymorphisms (SNPs) and short insertions and deletions (indels) in BCG Danish 1331 and BCG Pasteur 1173P2. These findings may help determine the effect of genetic variations on the variable protective efficacy of BCG vaccine strains.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08826-9.

Expansion of the pRAS3 Plasmid Family in Aeromonas salmonicida subsp. salmonicida and Growing Evidence of Interspecies Connections for These Plasmids

Aeromonas salmonicida subsp. salmonicida is a pathogenic bacterium responsible for furunculosis in salmonids. Following an outbreak of furunculosis, the infection can be treated with antibiotics, but it is common to observe ineffective treatment due to antibiotic resistance. This bacterium has a wide variety of plasmids responsible for this resistance. Among them, pRAS3 carries a tetracycline resistance gene. Several variants of this plasmid have been discovered over the years (pRAS3-3432 and pRAS3.1 to 3.4). During the present study, two new variants of the plasmid pRAS3 were identified (pRAS3.5 and pRAS3-3759) in strains of A. salmonicida subsp. salmonicida. Plasmid pRAS3-3759, which has been found in many strains from the same region over the past three years, has an additional genetic element identical to one found in pRAS3-3432. This genetic element was also found in Chlamydia suis, a swine pathogen. In this study, we analyzed the bacteria’s resistance to tetracycline, the number of copies of the plasmids, and the growth of the strains that carry five of the pRAS3 variants (pRAS3.3 to 3.5, pRAS3-3432, and pRAS3-3759). The results show no particular trend despite the differences between the plasmids, except for the resistance to tetracycline when analyzed in an isogenic background. Blast analysis also revealed the presence of pRAS3 plasmids in other bacterial species, which suggests that this plasmid family has widely spread. This study once again highlights the ability of A. salmonicida subsp. salmonicida to adapt to furunculosis antibiotic treatments, and the still-growing family of pRAS3 plasmids.

Genetic confirmation of a hybrid between two highly divergent cardinalid species: A rose-breasted grosbeak (Pheucticus ludovicianus) and a scarlet tanager (Piranga olivacea)

Using low-coverage whole-genome sequencing, analysis of vocalizations, and inferences from natural history, we document a first-generation hybrid between a rose-breasted grosbeak (Pheucticus ludovicianus) and a scarlet tanager (Piranga olivacea). These two species occur sympatrically throughout much of eastern North America, although were not previously known to interbreed. Following the field identification of a putative hybrid, we use genetic and bioacoustic data to show that a rose-breasted grosbeak was the maternal parent and a scarlet tanager was the paternal parent of the hybrid, whose song was similar to the latter species. These two species diverged >10 million years ago, and thus it is surprising to find a hybrid formed under natural conditions in the wild. Notably, the hybrid has an exceptionally heterozygous genome, with a conservative estimate of a heterozygous base every 100 bp. The observation that this hybrid of such highly divergent parental taxa has survived until adulthood serves as another example of the capacity for hybrid birds to survive with an exceptionally divergent genomic composition.

Robertsonian fusion and centromere repositioning contributed to the formation of satellite-free centromeres during the evolution of zebras

Centromeres are epigenetically specified by the histone H3 variant CENP-A and typically associated to highly repetitive satellite DNA. We previously discovered natural satellite-free neocentromeres in Equus caballus and E. asinus. Here, through ChIP-seq with an anti-CENP-A antibody, we found an extraordinarily high number of centromeres lacking satellite DNA in the zebras E. burchelli (15 of 22) and E. grevyi (13 of 23), demonstrating that the absence of satellite DNA at the majority of centromeres is compatible with genome stability and species survival and challenging the role of satellite DNA in centromere function. Nine satellite-free centromeres are shared between the two species in agreement with their recent separation. We assembled all centromeric regions and improved the reference genome of E. burchelli. Sequence analysis of the CENP-A binding domains revealed that they are LINE-1 and AT-rich with four of them showing DNA amplification. In the two zebras, satellite-free centromeres emerged from centromere repositioning or following Robertsonian fusion. In five chromosomes, the centromeric function arose near the fusion points, which are located within regions marked by traces of ancestral pericentromeric sequences. Therefore, besides centromere repositioning, Robertsonian fusions are an important source of satellite-free centromeres during evolution. Finally, in one case, a satellite-free centromere was seeded on an inversion breakpoint. At eleven chromosomes, whose primary constrictions seemed to be associated to satellite repeats by cytogenetic analysis, satellite-free neocentromeres were instead located near the ancestral inactivated satellite-based centromeres, therefore, the centromeric function has shifted away from a satellite repeat containing locus to a satellite-free new position.

Techniques for Profiling the Cellular Immune Response and Their Implications for Interventional Oncology

In recent years there has been increased interest in using the immune contexture of the primary tumors to predict the patient's prognosis. The tumor microenvironment of patients with cancers consists of different types of lymphocytes, tumor-infiltrating leukocytes, dendritic cells, and others. Different technologies can be used for the evaluation of the tumor microenvironment, all of which require a tissue or cell sample. Image-guided tissue sampling is a cornerstone in the diagnosis, stratification, and longitudinal evaluation of therapeutic efficacy for cancer patients receiving immunotherapies. Therefore, interventional radiologists (IRs) play an essential role in the evaluation of patients treated with systemically administered immunotherapies. This review provides a detailed description of different technologies used for immune assessment and analysis of the data collected from the use of these technologies. The detailed approach provided herein is intended to provide the reader with the knowledge necessary to not only interpret studies containing such data but also design and apply these tools for clinical practice and future research studies.

Leishmania parasites exchange drug-resistance genes through extracellular vesicles

Leishmania are eukaryotic parasites that have retained the ability to produce extracellular vesicles (EVs) through evolution. To date, it has been unclear if different DNA entities could be associated with Leishmania EVs and whether these could constitute a mechanism of horizontal gene transfer (HGT). Herein, we investigate the DNA content of EVs derived from drug-resistant parasites, as well as the EVs' potential to act as shuttles for DNA transfer. Next-generation sequencing and PCR assays confirm the enrichment of amplicons carrying drug-resistance genes associated with EVs. Transfer assays of drug-resistant EVs highlight a significant impact on the phenotype of recipient parasites induced by the expression of the transferred DNA. Recipient parasites display an enhanced growth and better control of oxidative stress. We provide evidence that eukaryotic EVs function as efficient mediators in HGT, thereby facilitating the transmission of drug-resistance genes and increasing the fitness of cells when encountering stressful environments.

Plasma contains ultrashort single-stranded DNA in addition to nucleosomal cell-free DNA

Plasma cell-free DNA is being widely explored as a biomarker for clinical screening. Currently, methods are optimized for the extraction and detection of double-stranded mononucleosomal cell-free DNA of ∼160bp in length. We introduce uscfDNA-seq, a single-stranded cell-free DNA next-generation sequencing pipeline, which bypasses previous limitations to reveal a population of ultrashort single-stranded cell-free DNA in human plasma. This species has a modal size of 50nt and is distinctly separate from mononucleosomal cell-free DNA. Treatment with single-stranded and double-stranded specific nucleases suggests that ultrashort cell-free DNA is primarily single-stranded. It is distributed evenly across chromosomes and has a similar distribution profile over functional elements as the genome, albeit with an enrichment over promoters, exons, and introns, which may be suggestive of a terminal state of genome degradation. The examination of this cfDNA species could reveal new features of cell death pathways or it can be used for cell-free DNA biomarker discovery.

A previously uncharacterized gene, PA2146, contributes to biofilm formation and drug tolerance across the ɣ-Proteobacteria

Transcriptomic studies have revealed a large number of uncharacterized genes that are differentially expressed in biofilms, which may be important in regulating biofilm phenotypes such as resistance to antimicrobial agents. To identify biofilm genes of unknown function in P. aeruginosa, we made use of RNA-seq and selected 27 uncharacterized genes that were induced upon biofilm growth. Biofilms by respective mutants were subsequently analyzed for two biofilm characteristics, the biofilm architecture and drug susceptibility. The screen revealed 12 out of 27 genes to contribute to biofilm formation and 13 drug susceptibility, with 8 genes affecting both biofilm phenotypes. Amongst the genes affecting both biofilm phenotypes was PA2146, encoding a small hypothetical protein that exhibited some of the most substantial increases in transcript abundance during biofilm growth by P. aeruginosa PAO1 and clinical isolates. PA2146 is highly conserved in ɣ-proteobacteria. Inactivation of PA2146 affected both biofilm phenotypes in P. aeruginosa PAO1, with inactivation of homologs in Klebsiella pneumoniae and Escherichia coli having similar effects. Heterologous expression of PA2146 homologs complemented the P. aeruginosa ∆PA2146, suggesting that PA2146 homologs substitute for and play a similar role as PA2146 in P. aeruginosa.

Enhancing Mesenchymal Stromal Cell Potency: Inflammatory Licensing via Mechanotransduction

Mesenchymal stromal cells (MSC) undergo functional maturation upon their migration from bone marrow and introduction to a site of injury. This inflammatory licensing leads to heightened immune regulation via cell-to-cell interaction and the secretion of immunomodulatory molecules, such as anti-inflammatory mediators and antioxidants. Pro-inflammatory cytokines are a recognized catalyst of inflammatory licensing; however, biomechanical forces, such as fluid shear stress, are a second, distinct class of stimuli that incite functional maturation. Here we show mechanotransduction, achieved by exposing MSC to various grades of wall shear stress (WSS) within a scalable conditioning platform, enhances the immunomodulatory potential of MSC independent of classical pro-inflammatory cytokines. A dose-dependent effect of WSS on potency is evidenced by production of prostaglandin E2 (PGE₂) and indoleamine 2,3 dioxygenase 1 (IDO1), as well as suppression of tumor necrosis factor-α (TNF- α) and interferon-γ (IFN-γ) production by activated immune cells. Consistent, reproducible licensing is demonstrated in adipose tissue and bone marrow human derived MSC without significant impact on cell viability, cellular yield, or identity. Transcriptome analysis of WSS-conditioned BM-MSC elucidates the broader phenotypic implications on the differential expression of immunomodulatory factors. These results suggest mechanotransduction as a viable, scalable pre-conditioning alternative to pro-inflammatory cytokines. Enhancing the immunomodulatory capacity of MSC via biomechanical conditioning represents a novel cell therapy manufacturing approach.

PARP inhibition is a modulator of anti-tumor immune response in BRCA-deficient tumors

PARP inhibitors are synthetically lethal with BRCA1/2 mutations, and in this setting, accumulation of DNA damage leads to cell death. Because increased DNA damage and subsequent immune activation can prime an anti-tumor immune response, we studied the impact of olaparib ± immune checkpoint blockade (ICB) on anti-tumor activity and the immune microenvironment. Concurrent combination of olaparib, at clinically relevant exposures, with ICB gave durable and deeper anti-tumor activity in the Brca1m BR5 model vs. monotherapies. Olaparib and combination treatment modulated the immune microenvironment, including increases in CD8+ T cells and NK cells, and upregulation of immune pathways, including type I IFN and STING signaling. Olaparib also induced a dose-dependent upregulation of immune pathways, including JAK/STAT, STING and type I IFN, in the tumor cell compartment of a BRCA1m (HBCx-10) but not a BRCA WT (HBCx-9) breast PDX model. In vitro, olaparib induced BRCAm tumor cell–specific dendritic cell transactivation. Relevance to human disease was assessed using patient samples from the MEDIOLA (NCT02734004) trial, which showed increased type I IFN, STING, and JAK/STAT pathway expression following olaparib treatment, in line with preclinical findings. These data together provide evidence for a mechanism and schedule underpinning potential benefit of ICB combination with olaparib.

Experimental evaluation of ecological principles to understand and modulate the outcome of bacterial strain competition in gut microbiomes

It is unclear if coexistence theory can be applied to gut microbiomes to understand their characteristics and modulate their composition. Through experiments in gnotobiotic mice with complex microbiomes, we demonstrated that strains of Akkermansia muciniphila and Bacteroides vulgatus could only be established if microbiomes were devoid of these species. Strains of A. muciniphila showed strict competitive exclusion, while B. vulgatus strains coexisted but populations were still influenced by competitive interactions. These differences in competitive behavior were reflective of genomic variation within the two species, indicating considerable niche overlap for A. muciniphila strains and a broader niche space for B. vulgatus strains. Priority effects were detected for both species as strains’ competitive fitness increased when colonizing first, which resulted in stable persistence of the A. muciniphila strain colonizing first and competitive exclusion of the strain arriving second. Based on these observations, we devised a subtractive strategy for A. muciniphila using antibiotics and showed that a strain from an assembled community can be stably replaced by another strain. By demonstrating that competitive outcomes in gut ecosystems depend on niche differences and are historically contingent, our study provides novel information to explain the ecological characteristics of gut microbiomes and a basis for their modulation.

Genetic basis of the lobedness degree in tomato fruit morphology

Fruit shape is a key trait in tomato (Solanum lycopersicum L.). Since most studies focused on proximo-distal fruit morphology, we hypothesized that unknown QTLs for medio-lateral direction ones could be found analysing segregating populations where major shape genes are fixed. We examined the diversity of fruit morphology in medio-lateral direction; defined divergent traits in cultivars carrying identical genetic constitution at LC and FAS genes; and identified QTLs for lobedness degree (LD) by a QTL-seq approach. We found that LC and FAS genes were not enough to explain LD variability in a large tomato collection. Then, we derived F₂ populations crossing cultivars divergent for LD where LC and FAS were fixed (Yellow Stuffer x Heinz 1439 [F₂YSxH] and Voyage x Old Brooks [F₂VxOB]). By QTL-seq we identified a QTL for LD on chromosome 8 in both F₂, which was validated in F₂YSxH by interval mapping accounting for ~ 17% of the variability. Other two QTLs located on chromosomes 6 and 11 with epistasis explained ~ 61% of the variability in the F₂VxOB. In conclusion, three novel QTLs with major effect for LD (ld6, ld8, and ld11) were identified through the study of diversity and genetic segregation in intraspecific tomato crosses.

Complete Genome Sequence of Bacillus frigoritolerans JHS1

Bacillus frigoritolerans JHS1 was isolated from the soil of a tomato plant (Solanum lycopersicum). The genome consists of one circular chromosome (5,552,463 bp) and a plasmid (16,118 bp) with an overall GC content of 40.57%. Using TYGS for taxonomic classification, strain JHS1 was assigned to the species Bacillus frigoritolerans.