Bat coronaviruses related to SARS-CoV-2 and infectious for human cells.

The animal reservoir of SARS-CoV-2 is unknown despite reports of SARS-CoV-2-related viruses in Asian Rhinolophus bats1-4, including the closest virus from R. affinis, RaTG13 (refs. 5,6), and pangolins7-9. SARS-CoV-2 has a mosaic genome, to which different progenitors contribute. The spike sequence determines the binding affinity and accessibility of its receptor-binding domain to the cellular angiotensin-converting enzyme 2 (ACE2) receptor and is responsible for host range10-12. SARS-CoV-2 progenitor bat viruses genetically close to SARS-CoV-2 and able to enter human cells through a human ACE2 (hACE2) pathway have not yet been identified, although they would be key in understanding the origin of the epidemic. Here we show that such viruses circulate in cave bats living in the limestone karstic terrain in northern Laos, in the Indochinese peninsula. We found that the receptor-binding domains of these viruses differ from that of SARS-CoV-2 by only one or two residues at the interface with ACE2, bind more efficiently to the hACE2 protein than that of the SARS-CoV-2 strain isolated in Wuhan from early human cases, and mediate hACE2-dependent entry and replication in human cells, which is inhibited by antibodies that neutralize SARS-CoV-2. None of these bat viruses contains a furin cleavage site in the spike protein. Our findings therefore indicate that bat-borne SARS-CoV-2-like viruses that are potentially infectious for humans circulate in Rhinolophus spp. in the Indochinese peninsula.

SARS-CoV-2-related bat virus behavior in human-relevant models sheds light on the origin of COVID-19.

Bat sarbecovirus BANAL-236 is highly related to SARS-CoV-2 and infects human cells, albeit lacking the furin cleavage site in its spike protein. BANAL-236 replicates efficiently and pauci-symptomatically in humanized mice and in macaques, where its tropism is enteric, strongly differing from that of SARS-CoV-2. BANAL-236 infection leads to protection against superinfection by a virulent strain. We find no evidence of antibodies recognizing bat sarbecoviruses in populations in close contact with bats in which the virus was identified, indicating that such spillover infections, if they occur, are rare. Six passages in humanized mice or in human intestinal cells, mimicking putative early spillover events, select adaptive mutations without appearance of a furin cleavage site and no change in virulence. Therefore, acquisition of a furin site in the spike protein is likely a pre-spillover event that did not occur upon replication of a SARS-CoV-2-like bat virus in humans or other animals. Other hypotheses regarding the origin of the SARS-CoV-2 should therefore be evaluated, including the presence of sarbecoviruses carrying a spike with a furin cleavage site in bats.

Circovirus Hepatitis in Immunocompromised Patient, Switzerland.

We identified a novel human circovirus in an immunocompromised 66-year-old woman with sudden onset of self-limiting hepatitis. We detected human circovirus 1 (HCirV-1) transcripts in hepatocytes and the HCirV-1 genome long-term in the patient's blood, stool, and urine. HCirV-1 is an emerging human pathogen that persists in susceptible patients.

Late-onset enteric virus infection associated with hepatitis (EVAH) in transplanted SCID patients.

BACKGROUND: Allogenic hematopoietic stem cell transplantation (HSCT) and gene therapy (GT) are potentially curative treatments for severe combined immunodeficiency (SCID). Late-onset posttreatment manifestations (such as persistent hepatitis) are not uncommon. OBJECTIVE: We sought to characterize the prevalence and pathophysiology of persistent hepatitis in transplanted SCID patients (SCIDH+) and to evaluate risk factors and treatments. METHODS: We used various techniques (including pathology assessments, metagenomics, single-cell transcriptomics, and cytometry by time of flight) to perform an in-depth study of different tissues from patients in the SCIDH+ group and corresponding asymptomatic similarly transplanted SCID patients without hepatitis (SCIDH-). RESULTS: Eleven patients developed persistent hepatitis (median of 6 years after HSCT or GT). This condition was associated with the chronic detection of enteric viruses (human Aichi virus, norovirus, and sapovirus) in liver and/or stools, which were not found in stools from the SCIDH- group (n = 12). Multiomics analysis identified an expansion of effector memory CD8+ T cells with high type I and II interferon signatures. Hepatitis was associated with absence of myeloablation during conditioning, split chimerism, and defective B-cell function, representing 25% of the 44 patients with SCID having these characteristics. Partially myeloablative retransplantation or GT of patients with this condition (which we have named as "enteric virus infection associated with hepatitis") led to the reconstitution of T- and B-cell immunity and remission of hepatitis in 5 patients, concomitantly with viral clearance. CONCLUSIONS: Enteric virus infection associated with hepatitis is related to chronic enteric viral infection and immune dysregulation and is an important risk for transplanted SCID patients with defective B-cell function.

Chronic Aichi Virus Infection As a Cause of Long-Lasting Multiorgan Involvement in Patients With Primary Immune Deficiencies.

BACKGROUND: Metagenomic next-generation sequencing (mNGS) was used to assess patients with primary or secondary immune deficiencies (PIDs and SIDs) who presented with immunopathological conditions related to immunodysregulation. METHODS: Thirty patients with PIDs or SIDs who presented with symptoms related to immunodysregulation and 59 asymptomatic patients with similar PIDs or SIDs were enrolled. mNGS was performed on organ biopsy. Specific Aichi virus (AiV) reverse-transcription polymerase chain reaction (RT-PCR) was used to confirm AiV infection and screen the other patients. In situ hybridization (ISH) assay was done on AiV-infected organs to identify infected cells. Virus genotype was determined by phylogenetic analysis. RESULTS: AiV sequences were detected using mNGS in tissue samples of 5 patients and by RT-PCR in peripheral samples of another patient, all of whom presented with PID and long-lasting multiorgan involvement, including hepatitis, splenomegaly, and nephritis in 4 patients. CD8+ T-cell infiltration was a hallmark of the disease. RT-PCR detected intermittent low viral loads in urine and plasma from infected patients but not from uninfected patients. Viral detection stopped after immune reconstitution obtained by hematopoietic stem cell transplantation. ISH demonstrated the presence of AiV RNA in hepatocytes (n = 1) and spleen tissue (n = 2). AiV belonged to genotype A (n = 2) or B (n = 3). CONCLUSIONS: The similarity of the clinical presentation, the detection of AiV in a subgroup of patients suffering from immunodysregulation, the absence of AiV in asymptomatic patients, the detection of viral genome in infected organs by ISH, and the reversibility of symptoms after treatment argue for AiV causality.

Circovirus Hepatitis Infection in Heart-Lung Transplant Patient, France.

In March 2022, a 61-year-old woman in France who had received a heart-lung transplant sought treatment with chronic hepatitis mainly characterized by increased liver enzymes. After ruling out common etiologies, we used metagenomic next-generation sequencing to analyze a liver biopsy sample and identified an unknown species of circovirus, tentatively named human circovirus 1 (HCirV-1). We found no other viral or bacterial sequences. HCirV-1 shared 70% amino acid identity with the closest known viral sequences. The viral genome was undetectable in blood samples from 2017-2019, then became detectable at low levels in September 2020 and peaked at very high titers (1010 genome copies/mL) in January 2022. In March 2022, we found >108 genome copies/g or mL in the liver and blood, concomitant with hepatic cytolysis. We detected HCirV-1 transcripts in 2% of hepatocytes, demonstrating viral replication and supporting the role of HCirV-1 in liver damage.

First Case of Lethal Encephalitis in Western Europe Due to European Bat Lyssavirus Type 1.

BACKGROUND: Inaccurate diagnosis of encephalitis is a major issue as immunosuppressive treatments can be deleterious in case of viral infection. The European bat lyssavirus type 1 (EBLV-1), a virus related to rabies virus, is endemic in European bats. No human case has yet been reported in Western Europe. A 59-year-old patient without specific past medical history died from encephalitis. A colony of bats lived in an outbuilding of his house. No diagnosis was made using standard procedures. METHODS: We used a next generation sequencing (NGS) based transcriptomic protocol to search for pathogens in autopsy samples (meninges and brain frontal lobe). Results were confirmed by polymerase chain reaction (PCR) and by antibody testing in serum. Immunochemistry was used to characterize inflammatory cells and viral antigens in brain lesions. Cells and mice were inoculated with brain extracts for virus isolation. RESULTS: The patient's brain lesions were severe and diffuse in white and gray matter. Perivascular inflammatory infiltrates were abundant and rich in plasma cells. NGS identified European bat lyssavirus type 1a in brain, which was confirmed by PCR. A high titer of neutralizing antibodies was found in serum. No viral antigen was detected, and the virus could not be isolated by cell culture or by mouse inoculation. CONCLUSIONS: The patient died from European bat lyssavirus type 1a infection. NGS was key to identifying this unexpected viral etiology in an epidemiological context that did not suggest rabies. People exposed to bats should be strongly advised to be vaccinated with rabies vaccines, which are effective against EBLV-1.

Colorectal cancer specific conditions promote Streptococcus gallolyticus gut colonization.

Colonization by Streptococcus gallolyticus subsp. gallolyticus (SGG) is strongly associated with the occurrence of colorectal cancer (CRC). However, the factors leading to its successful colonization are unknown, and whether SGG influences the oncogenic process or benefits from the tumor-prone environment to prevail remains an open question. Here, we elucidate crucial steps that explain how CRC favors SGG colonization. By using mice genetically prone to CRC, we show that SGG colonization is 1,000-fold higher in tumor-bearing mice than in normal mice. This selective advantage occurs at the expense of resident intestinal enterococci. An SGG-specific locus encoding a bacteriocin ("gallocin") is shown to kill enterococci in vitro. Importantly, bile acids strongly enhance this bacteriocin activity in vivo, leading to greater SGG colonization. Constitutive activation of the Wnt pathway, one of the earliest signaling alterations in CRC, and the decreased expression of the bile acid apical transporter gene Slc10A2, as an effect of the Apc founding mutation, may thereby sustain intestinal colonization by SGG. We conclude that CRC-specific conditions promote SGG colonization of the gut by replacing commensal enterococci in their niche.

OSBPL10, RXRA and lipid metabolism confer African-ancestry protection against dengue haemorrhagic fever in admixed Cubans.

Ethnic groups can display differential genetic susceptibility to infectious diseases. The arthropod-born viral dengue disease is one such disease, with empirical and limited genetic evidence showing that African ancestry may be protective against the haemorrhagic phenotype. Global ancestry analysis based on high-throughput genotyping in admixed populations can be used to test this hypothesis, while admixture mapping can map candidate protective genes. A Cuban dengue fever cohort was genotyped using a 2.5 million SNP chip. Global ancestry was ascertained through ADMIXTURE and used in a fine-matched corrected association study, while local ancestry was inferred by the RFMix algorithm. The expression of candidate genes was evaluated by RT-PCR in a Cuban dengue patient cohort and gene set enrichment analysis was performed in a Thai dengue transcriptome. OSBPL10 and RXRA candidate genes were identified, with most significant SNPs placed in inferred weak enhancers, promoters and lncRNAs. OSBPL10 had significantly lower expression in Africans than Europeans, while for RXRA several SNPs may differentially regulate its transcription between Africans and Europeans. Their expression was confirmed to change through dengue disease progression in Cuban patients and to vary with disease severity in a Thai transcriptome dataset. These genes interact in the LXR/RXR activation pathway that integrates lipid metabolism and immune functions, being a key player in dengue virus entrance into cells, its replication therein and in cytokine production. Knockdown of OSBPL10 expression in THP-1 cells by two shRNAs followed by DENV2 infection tests led to a significant reduction in DENV replication, being a direct functional proof that the lower OSBPL10 expression profile in Africans protects this ancestry against dengue disease.

High-fat diet modifies the PPAR-γ pathway leading to disruption of microbial and physiological ecosystem in murine small intestine.

Diet is among the most important factors contributing to intestinal homeostasis, and basic functions performed by the small intestine need to be tightly preserved to maintain health. Little is known about the direct impact of high-fat (HF) diet on small-intestinal mucosal defenses and spatial distribution of the microbiota during the early phase of its administration. We observed that only 30 d after HF diet initiation, the intervillous zone of the ileum-which is usually described as free of bacteria-became occupied by a dense microbiota. In addition to affecting its spatial distribution, HF diet also drastically affected microbiota composition with a profile characterized by the expansion of Firmicutes (appearance of Erysipelotrichi), Proteobacteria (Desulfovibrionales) and Verrucomicrobia, and decrease of Bacteroidetes (family S24-7) and Candidatus arthromitus A decrease in antimicrobial peptide expression was predominantly observed in the ileum where bacterial density appeared highest. In addition, HF diet increased intestinal permeability and decreased cystic fibrosis transmembrane conductance regulator (Cftr) and the Na-K-2Cl cotransporter 1 (Nkcc1) gene and protein expressions, leading to a decrease in ileal secretion of chloride, likely responsible for massive alteration in mucus phenotype. This complex phenotype triggered by HF diet at the interface between the microbiota and the mucosal surface was reversed when the diet was switched back to standard composition or when mice were treated for 1 wk with rosiglitazone, a specific agonist of peroxisome proliferator-activated receptor-γ (PPAR-γ). Moreover, weaker expression of antimicrobial peptide-encoding genes and intervillous bacterial colonization were observed in Ppar-γ-deficient mice, highlighting the major role of lipids in modulation of mucosal immune defenses.

Bacteriocin from epidemic Listeria strains alters the host intestinal microbiota to favor infection.

Listeria monocytogenes is responsible for gastroenteritis in healthy individuals and for a severe invasive disease in immunocompromised patients. Among the three identified L. monocytogenes evolutionary lineages, lineage I strains are overrepresented in epidemic listeriosis outbreaks, but the mechanisms underlying the higher virulence potential of strains of this lineage remain elusive. Here, we demonstrate that Listeriolysin S (LLS), a virulence factor only present in a subset of lineage I strains, is a bacteriocin highly expressed in the intestine of orally infected mice that alters the host intestinal microbiota and promotes intestinal colonization by L. monocytogenes, as well as deeper organ infection. To our knowledge, these results therefore identify LLS as the first bacteriocin described in L. monocytogenes and associate modulation of host microbiota by L. monocytogenes epidemic strains to increased virulence.

Meta-analysis of SHANK Mutations in Autism Spectrum Disorders: a gradient of severity in cognitive impairments.

SHANK genes code for scaffold proteins located at the post-synaptic density of glutamatergic synapses. In neurons, SHANK2 and SHANK3 have a positive effect on the induction and maturation of dendritic spines, whereas SHANK1 induces the enlargement of spine heads. Mutations in SHANK genes have been associated with autism spectrum disorders (ASD), but their prevalence and clinical relevance remain to be determined. Here, we performed a new screen and a meta-analysis of SHANK copy-number and coding-sequence variants in ASD. Copy-number variants were analyzed in 5,657 patients and 19,163 controls, coding-sequence variants were ascertained in 760 to 2,147 patients and 492 to 1,090 controls (depending on the gene), and, individuals carrying de novo or truncating SHANK mutations underwent an extensive clinical investigation. Copy-number variants and truncating mutations in SHANK genes were present in ∼1% of patients with ASD: mutations in SHANK1 were rare (0.04%) and present in males with normal IQ and autism; mutations in SHANK2 were present in 0.17% of patients with ASD and mild intellectual disability; mutations in SHANK3 were present in 0.69% of patients with ASD and up to 2.12% of the cases with moderate to profound intellectual disability. In summary, mutations of the SHANK genes were detected in the whole spectrum of autism with a gradient of severity in cognitive impairment. Given the rare frequency of SHANK1 and SHANK2 deleterious mutations, the clinical relevance of these genes remains to be ascertained. In contrast, the frequency and the penetrance of SHANK3 mutations in individuals with ASD and intellectual disability-more than 1 in 50-warrant its consideration for mutation screening in clinical practice.

Estimation of Recent and Ancient Inbreeding in a Small Endogamous Tunisian Community Through Genomic Runs of Homozygosity.

Runs of homozygosity (ROHs) are extended genomic regions of homozygous genotypes that record populations' mating patterns in the past. We performed microarray genotyping on 15 individuals from a small isolated Tunisian community. We estimated the individual and population genome-wide level of homozygosity from data on ROH above 0.5 Mb in length. We found a high average number of ROH per individual (48.2). The smallest ROH category (0.5-1.49 Mb) represents 0.93% of the whole genome, while medium-size (1.5-4.99 Mb) and long-size ROH (≥5 Mb) cover 1.18% and 0.95%, respectively. We found that genealogical individual inbreeding coefficients (Fped ) based on three- to four-generation pedigrees are not reliable indicators of the current proportion of genome-wide homozygosity inferred from ROH (FROH ) either for 0.5 or 1.5 Mb ROH length thresholds, while identity-by-descent sharing is a function of shared coancestry. This study emphasizes the effect of reproductive isolation and a prolonged practice of consanguinity that limits the genetic heterogeneity. It also provides evidence of both recent and ancient parental relatedness contribution to the current level of genome-wide homozygosity in the studied population. These findings may be useful for evaluation of long-term effects of inbreeding on human health and for future applications of ROHs in identifying recessive susceptibility genes.

The Listeria transcriptional landscape from saprophytism to virulence.

The bacterium Listeria monocytogenes is ubiquitous in the environment and can lead to severe food-borne infections. It has recently emerged as a multifaceted model in pathogenesis. However, how this bacterium switches from a saprophyte to a pathogen is largely unknown. Here, using tiling arrays and RNAs from wild-type and mutant bacteria grown in vitro, ex vivo and in vivo, we have analysed the transcription of its entire genome. We provide the complete Listeria operon map and have uncovered far more diverse types of RNAs than expected: in addition to 50 small RNAs (<500 nucleotides), at least two of which are involved in virulence in mice, we have identified antisense RNAs covering several open-reading frames and long overlapping 5' and 3' untranslated regions. We discovered that riboswitches can act as terminators for upstream genes. When Listeria reaches the host intestinal lumen, an extensive transcriptional reshaping occurs with a SigB-mediated activation of virulence genes. In contrast, in the blood, PrfA controls transcription of virulence genes. Remarkably, several non-coding RNAs absent in the non-pathogenic species Listeria innocua exhibit the same expression patterns as the virulence genes. Together, our data unravel successive and coordinated global transcriptional changes during infection and point to previously unknown regulatory mechanisms in bacteria.

Genetic and functional analyses of SHANK2 mutations suggest a multiple hit model of autism spectrum disorders.

Autism spectrum disorders (ASD) are a heterogeneous group of neurodevelopmental disorders with a complex inheritance pattern. While many rare variants in synaptic proteins have been identified in patients with ASD, little is known about their effects at the synapse and their interactions with other genetic variations. Here, following the discovery of two de novo SHANK2 deletions by the Autism Genome Project, we identified a novel 421 kb de novo SHANK2 deletion in a patient with autism. We then sequenced SHANK2 in 455 patients with ASD and 431 controls and integrated these results with those reported by Berkel et al. 2010 (n = 396 patients and n = 659 controls). We observed a significant enrichment of variants affecting conserved amino acids in 29 of 851 (3.4%) patients and in 16 of 1,090 (1.5%) controls (P = 0.004, OR = 2.37, 95% CI = 1.23-4.70). In neuronal cell cultures, the variants identified in patients were associated with a reduced synaptic density at dendrites compared to the variants only detected in controls (P = 0.0013). Interestingly, the three patients with de novo SHANK2 deletions also carried inherited CNVs at 15q11-q13 previously associated with neuropsychiatric disorders. In two cases, the nicotinic receptor CHRNA7 was duplicated and in one case the synaptic translation repressor CYFIP1 was deleted. These results strengthen the role of synaptic gene dysfunction in ASD but also highlight the presence of putative modifier genes, which is in keeping with the "multiple hit model" for ASD. A better knowledge of these genetic interactions will be necessary to understand the complex inheritance pattern of ASD.

Genome-wide analysis of heteroduplex DNA in mismatch repair-deficient yeast cells reveals novel properties of meiotic recombination pathways.

Meiotic DNA double-strand breaks (DSBs) initiate crossover (CO) recombination, which is necessary for accurate chromosome segregation, but DSBs may also repair as non-crossovers (NCOs). Multiple recombination pathways with specific intermediates are expected to lead to COs and NCOs. We revisited the mechanisms of meiotic DSB repair and the regulation of CO formation, by conducting a genome-wide analysis of strand-transfer intermediates associated with recombination events. We performed this analysis in a SK1 × S288C Saccharomyces cerevisiae hybrid lacking the mismatch repair (MMR) protein Msh2, to allow efficient detection of heteroduplex DNAs (hDNAs). First, we observed that the anti-recombinogenic activity of MMR is responsible for a 20% drop in CO number, suggesting that in MMR-proficient cells some DSBs are repaired using the sister chromatid as a template when polymorphisms are present. Second, we observed that a large fraction of NCOs were associated with trans-hDNA tracts constrained to a single chromatid. This unexpected finding is compatible with dissolution of double Holliday junctions (dHJs) during repair, and it suggests the existence of a novel control point for CO formation at the level of the dHJ intermediate, in addition to the previously described control point before the dHJ formation step. Finally, we observed that COs are associated with complex hDNA patterns, confirming that the canonical double-strand break repair model is not sufficient to explain the formation of most COs. We propose that multiple factors contribute to the complexity of recombination intermediates. These factors include repair of nicks and double-stranded gaps, template switches between non-sister and sister chromatids, and HJ branch migration. Finally, the good correlation between the strand transfer properties observed in the absence of and in the presence of Msh2 suggests that the intermediates detected in the absence of Msh2 reflect normal intermediates.

Performance of clinical metagenomics in France: a prospective observational study.

BACKGROUND: Metagenomic next-generation sequencing (mNGS) allows untargeted identification of a broad range of pathogens, including rare or novel microorganisms. Despite the recognition of mNGS as a valuable diagnostic tool for infections, the most relevant indications for this innovative strategy remain poorly defined. We aimed to assess the determinants of positivity and clinical utility of mNGS. METHODS: In this observational study, we prospectively performed short-read shotgun metagenomics analysis as a second-line test (in cases of negative first-line test or when the symptoms were not fully explained by initial positive results) or as a first-line test in life-threatening situations requiring urgent non-targeted pathogen identification at the Necker-Enfants Malades Hospital (Paris, France). All sample types, clinical indications, and patient populations were included. Samples were accompanied by a mandatory form completed by the senior clinician or pathologist, on which the clinical level of suspected infection (defined as high or low) was indicated. We assessed the variables (gender, age, immune status, initial suspicion of infection, indication, and sample type) associated with mNGS pathogen detection using odds ratios (ORs) from multivariate logistic regression. Additional investigations were carried out using specific PCR or culture techniques, to confirm positive mNGS results, or when infectious suspicion was particularly high despite a negative mNGS result. FINDINGS: Between Oct 29, 2019, and Nov 7, 2022, we analysed 742 samples collected from 523 patients. The initial suspicion of infection was either high (n=470, 63%) or low (n=272, 37%). Causative or possibly causative pathogens were detected in 117 (25%) samples from patients with high initial suspicion of infection, versus nine (3%) samples analysed to rule out infection (OR 9·1, 95% CI 4·6-20·4; p<0·0001). We showed that mNGS had higher odds of detecting a causative or possibly causative pathogenic virus on CNS biopsies than CSF samples (4·1, 1·7-10·7; p=0·0025) and in samples from immunodeficient compared with immunocompetent individuals (2·4, 1·4-4·1; p=0·0013). Concordance with conventional confirmatory tests results was 103 (97%) of 106, when mNGS detected causative or possibly causative pathogens. Altogether, among 231 samples investigated by both mNGS and subsequent specific tests, discordant results were found in 69 (30%) samples, of which 58 (84%) were mNGS positive and specific tests negative, and 11 (16%) mNGS negative and specific tests positive. INTERPRETATION: Major determinants of pathogen detection by mNGS are immune status and initial level of suspicion of infection. These findings will contribute, along with future studies, to refining the positioning of mNGS in diagnostic and treatment decision-making algorithms. FUNDING: Necker-Enfants Malades Hospital and Institut Pasteur. TRANSLATION: For the French translation of the abstract see Supplementary Materials section.

Comparative analysis of resistant and susceptible macrophage gene expression response to Leishmania major parasite.

BACKGROUND: Leishmania are obligated intracellular pathogens that replicate almost exclusively in macrophages. The outcome of infection depends largely on parasite pathogenicity and virulence but also on the activation status and genetic background of macrophages. Animal models are essential for a better understanding of pathogenesis of different microbes including Leishmania. RESULTS: Here we compared the transcriptional signatures of resistant (C57BL/6) and susceptible (BALB/c) mouse bone marrow-derived macrophages in response to Leishmania major (L. major) promastigotes infection.Microarray results were first analyzed for significant pathways using the Kyoto Encylopedia of Genes and Genomes (KEGG) database. The analysis revealed that a large set of the shared genes is involved in the immune response and that difference in the expression level of some chemokines and chemokine receptors could partially explain differences in resistance. We next focused on up-regulated genes unique to either BALB/c or C57BL/6 derived macrophages and identified, using KEGG database, signal transduction pathways among the most relevant pathways unique to both susceptible and resistant derived macrophages. Indeed, genes unique to C57BL/6 BMdMs were associated with target of rapamycin (mTOR) signaling pathway while a range of genes unique to BALB/c BMdMs, belong to p53 signaling pathway. We next investigated whether, in a given mice strain derived macrophages, the different up-regulated unique genes could be coordinately regulated. Using GeneMapp Cytoscape, we showed that the induced genes unique to BALB/c or C57BL/6 BMdMs are interconnected. Finally, we examined whether the induced pathways unique to BALB/c derived macrophages interfere with the ones unique to C57BL/6 derived macrophages. Protein-protein interaction analysis using String database highlights the existence of a cross-talk between p53 and mTOR signaling pathways respectively specific to susceptible and resistant BMdMs. CONCLUSIONS: Taken together our results suggest that strains specific pathogenesis may be due to a difference in the magnitude of the same pathways and/or to differentially expressed pathways in the two mouse strains derived macrophages. We identify signal transduction pathways among the most relevant pathways modulated by L. major infection, unique to BALB/c and C57BL/6 BMdM and postulate that the interplay between these potentially interconnected pathways could direct the macrophage response toward a given phenotype.

Virus Pop-Expanding Viral Databases by Protein Sequence Simulation.

The improvement of our knowledge of the virosphere, which includes unknown viruses, is a key area in virology. Metagenomics tools, which perform taxonomic assignation from high throughput sequencing datasets, are generally evaluated with datasets derived from biological samples or in silico spiked samples containing known viral sequences present in public databases, resulting in the inability to evaluate the capacity of these tools to detect novel or distant viruses. Simulating realistic evolutionary directions is therefore key to benchmark and improve these tools. Additionally, expanding current databases with realistic simulated sequences can improve the capacity of alignment-based searching strategies for finding distant viruses, which could lead to a better characterization of the "dark matter" of metagenomics data. Here, we present Virus Pop, a novel pipeline for simulating realistic protein sequences and adding new branches to a protein phylogenetic tree. The tool generates simulated sequences with substitution rate variations that are dependent on protein domains and inferred from the input dataset, allowing for a realistic representation of protein evolution. The pipeline also infers ancestral sequences corresponding to multiple internal nodes of the input data phylogenetic tree, enabling new sequences to be inserted at various points of interest in the group studied. We demonstrated that Virus Pop produces simulated sequences that closely match the structural and functional characteristics of real protein sequences, taking as an example the spike protein of sarbecoviruses. Virus Pop also succeeded at creating sequences that resemble real sequences not included in the databases, which facilitated the identification of a novel pathogenic human circovirus not included in the input database. In conclusion, Virus Pop is helpful for challenging taxonomic assignation tools and could help improve databases to better detect distant viruses.