The phosphorylation landscape of infection-related development by the rice blast fungus.

Many of the world's most devastating crop diseases are caused by fungal pathogens that elaborate specialized infection structures to invade plant tissue. Here, we present a quantitative mass-spectrometry-based phosphoproteomic analysis of infection-related development by the rice blast fungus Magnaporthe oryzae, which threatens global food security. We mapped 8,005 phosphosites on 2,062 fungal proteins following germination on a hydrophobic surface, revealing major re-wiring of phosphorylation-based signaling cascades during appressorium development. Comparing phosphosite conservation across 41 fungal species reveals phosphorylation signatures specifically associated with biotrophic and hemibiotrophic fungal infection. We then used parallel reaction monitoring (PRM) to identify phosphoproteins regulated by the fungal Pmk1 MAPK that controls plant infection by M. oryzae. We define 32 substrates of Pmk1 and show that Pmk1-dependent phosphorylation of regulator Vts1 is required for rice blast disease. Defining the phosphorylation landscape of infection therefore identifies potential therapeutic interventions for the control of plant diseases.

An atlas of human vector-borne microbe interactions reveals pathogenicity mechanisms.

Vector-borne diseases are a leading cause of death worldwide and pose a substantial unmet medical need. Pathogens binding to host extracellular proteins (the "exoproteome") represents a crucial interface in the etiology of vector-borne disease. Here, we used bacterial selection to elucidate host-microbe interactions in high throughput (BASEHIT)-a technique enabling interrogation of microbial interactions with 3,324 human exoproteins-to profile the interactomes of 82 human-pathogen samples, including 30 strains of arthropod-borne pathogens and 8 strains of related non-vector-borne pathogens. The resulting atlas revealed 1,303 putative interactions, including hundreds of pairings with potential roles in pathogenesis, including cell invasion, tissue colonization, immune evasion, and host sensing. Subsequent functional investigations uncovered that Lyme disease spirochetes recognize epidermal growth factor as an environmental cue of transcriptional regulation and that conserved interactions between intracellular pathogens and thioredoxins facilitate cell invasion. In summary, this interactome atlas provides molecular-level insights into microbial pathogenesis and reveals potential host-directed targets for next-generation therapeutics.

Dual pathogenicity island transfer by piggybacking lateral transduction.

Lateral transduction (LT) is the process by which temperate phages mobilize large sections of bacterial genomes. Despite its importance, LT has only been observed during prophage induction. Here, we report that superantigen-carrying staphylococcal pathogenicity islands (SaPIs) employ a related but more versatile and complex mechanism of gene transfer to drive chromosomal hypermobility while self-transferring with additional virulence genes from the host. We found that after phage infection or prophage induction, activated SaPIs form concatamers in the bacterial chromosome by switching between parallel genomic tracks in replication bubbles. This dynamic life cycle enables SaPIbov1 to piggyback its LT of staphylococcal pathogenicity island vSaα, which encodes an array of genes involved in host-pathogen interactions, allowing both islands to be mobilized intact and transferred in a single infective particle. Our findings highlight previously unknown roles of pathogenicity islands in bacterial virulence and show that their evolutionary impact extends beyond the genes they carry.

Virological characteristics of the SARS-CoV-2 Omicron BA.2 spike.

Soon after the emergence and global spread of the SARS-CoV-2 Omicron lineage BA.1, another Omicron lineage, BA.2, began outcompeting BA.1. The results of statistical analysis showed that the effective reproduction number of BA.2 is 1.4-fold higher than that of BA.1. Neutralization experiments revealed that immunity induced by COVID vaccines widely administered to human populations is not effective against BA.2, similar to BA.1, and that the antigenicity of BA.2 is notably different from that of BA.1. Cell culture experiments showed that the BA.2 spike confers higher replication efficacy in human nasal epithelial cells and is more efficient in mediating syncytia formation than the BA.1 spike. Furthermore, infection experiments using hamsters indicated that the BA.2 spike-bearing virus is more pathogenic than the BA.1 spike-bearing virus. Altogether, the results of our multiscale investigations suggest that the risk of BA.2 to global health is potentially higher than that of BA.1.

Virological characteristics of the SARS-CoV-2 Omicron BA.2 subvariants, including BA.4 and BA.5.

After the global spread of the SARS-CoV-2 Omicron BA.2, some BA.2 subvariants, including BA.2.9.1, BA.2.11, BA.2.12.1, BA.4, and BA.5, emerged in multiple countries. Our statistical analysis showed that the effective reproduction numbers of these BA.2 subvariants are greater than that of the original BA.2. Neutralization experiments revealed that the immunity induced by BA.1/2 infections is less effective against BA.4/5. Cell culture experiments showed that BA.2.12.1 and BA.4/5 replicate more efficiently in human alveolar epithelial cells than BA.2, and particularly, BA.4/5 is more fusogenic than BA.2. We further provided the structure of the BA.4/5 spike receptor-binding domain that binds to human ACE2 and considered how the substitutions in the BA.4/5 spike play roles in ACE2 binding and immune evasion. Moreover, experiments using hamsters suggested that BA.4/5 is more pathogenic than BA.2. Our multiscale investigations suggest that the risk of BA.2 subvariants, particularly BA.4/5, to global health is greater than that of original BA.2.

Pathogenic Germline Variants in 10,389 Adult Cancers.

We conducted the largest investigation of predisposition variants in cancer to date, discovering 853 pathogenic or likely pathogenic variants in 8% of 10,389 cases from 33 cancer types. Twenty-one genes showed single or cross-cancer associations, including novel associations of SDHA in melanoma and PALB2 in stomach adenocarcinoma. The 659 predisposition variants and 18 additional large deletions in tumor suppressors, including ATM, BRCA1, and NF1, showed low gene expression and frequent (43%) loss of heterozygosity or biallelic two-hit events. We also discovered 33 such variants in oncogenes, including missenses in MET, RET, and PTPN11 associated with high gene expression. We nominated 47 additional predisposition variants from prioritized VUSs supported by multiple evidences involving case-control frequency, loss of heterozygosity, expression effect, and co-localization with mutations and modified residues. Our integrative approach links rare predisposition variants to functional consequences, informing future guidelines of variant classification and germline genetic testing in cancer.

Structural Basis for Virulence Activation of Francisella tularensis.

The bacterium Francisella tularensis (Ft) is one of the most infectious agents known. Ft virulence is controlled by a unique combination of transcription regulators: the MglA-SspA heterodimer, PigR, and the stress signal, ppGpp. MglA-SspA assembles with the σ70-associated RNAP holoenzyme (RNAPσ70), forming a virulence-specialized polymerase. These factors activate Francisella pathogenicity island (FPI) gene expression, which is required for virulence, but the mechanism is unknown. Here we report FtRNAPσ70-promoter-DNA, FtRNAPσ70-(MglA-SspA)-promoter DNA, and FtRNAPσ70-(MglA-SspA)-ppGpp-PigR-promoter DNA cryo-EM structures. Structural and genetic analyses show MglA-SspA facilitates σ70 binding to DNA to regulate virulence and virulence-enhancing genes. Our Escherichia coli RNAPσ70-homodimeric EcSspA structure suggests this is a general SspA-transcription regulation mechanism. Strikingly, our FtRNAPσ70-(MglA-SspA)-ppGpp-PigR-DNA structure reveals ppGpp binding to MglA-SspA tethers PigR to promoters. PigR in turn recruits FtRNAP αCTDs to DNA UP elements. Thus, these studies unveil a unique mechanism for Ft pathogenesis involving a virulence-specialized RNAP that employs two (MglA-SspA)-based strategies to activate virulence genes.

The phc Quorum-Sensing System in Ralstonia solanacearum Species Complex.

Ralstonia solanacearum species complex (RSSC) strains are devastating plant pathogens distributed worldwide. The primary cell density-dependent gene expression system in RSSC strains is phc quorum sensing (QS). It regulates the expression of about 30% of all genes, including those related to cellular activity, primary and secondary metabolism, pathogenicity, and more. The phc regulatory elements encoded by the phcBSRQ operon and phcA gene play vital roles. RSSC strains use methyl 3-hydroxymyristate (3-OH MAME) or methyl 3-hydroxypalmitate (3-OH PAME) as the QS signal. Each type of RSSC strain has specificity in generating and receiving its QS signal, but their signaling pathways might not differ significantly. In this review, I describe the genetic and biochemical factors involved in QS signal input and the regulatory network and summarize control of the phc QS system, new cell-cell communications, and QS-dependent interactions with soil fungi.

Evolution of Tuberculosis Pathogenesis.

Mycobacterium tuberculosis is a globally distributed, lethal pathogen of humans. The virulence armamentarium of M. tuberculosis appears to have been developed on a scaffold of antiphagocytic defenses found among diverse, mostly free-living species of Mycobacterium. Pathoadaptation was further aided by the modularity, flexibility, and interactivity characterizing mycobacterial effectors and their regulators. During emergence of M. tuberculosis, novel genetic material was acquired, created, and integrated with existing tools. The major mutational mechanisms underlying these adaptations are discussed in this review, with examples. During its evolution, M. tuberculosis lost the ability and/or opportunity to engage in lateral gene transfer, but despite this it has retained the adaptability that characterizes mycobacteria. M. tuberculosis exemplifies the evolutionary genomic mechanisms underlying adoption of the pathogenic niche, and studies of its evolution have uncovered a rich array of discoveries about how new pathogens are made.

Malassezia: A Commensal, Pathogen, and Mutualist of Human and Animal Skin.

Identified in the late nineteenth century as a single species residing on human skin, Malassezia is now recognized as a diverse genus comprising 18 species inhabiting not only skin but human gut, hospital environments, and even deep-sea sponges. All cultivated Malassezia species are lipid dependent, having lost genes for lipid synthesis and carbohydrate metabolism. The surging interest in Malassezia results from development of tools to improve sampling, culture, identification, and genetic engineering, which has led to findings implicating it in numerous skin diseases, Crohn disease, and pancreatic cancer. However, it has become clear that Malassezia plays a multifaceted role in human health, with mutualistic activity in atopic dermatitis and a preventive effect against other skin infections due to its potential to compete with skin pathogens such as Candida auris. Improved understanding of complex microbe-microbe and host-microbe interactions will be required to define Malassezia's role in human and animal health and disease so as to design targeted interventions.

Hijacking the Hijackers: Escherichia coli Pathogenicity Islands Redirect Helper Phage Packaging for Their Own Benefit.

Phage-inducible chromosomal islands (PICIs) represent a novel and universal class of mobile genetic elements, which have broad impact on bacterial virulence. In spite of their relevance, how the Gram-negative PICIs hijack the phage machinery for their own specific packaging and how they block phage reproduction remains to be determined. Using genetic and structural analyses, we solve the mystery here by showing that the Gram-negative PICIs encode a protein that simultaneously performs these processes. This protein, which we have named Rpp (for redirecting phage packaging), interacts with the phage terminase small subunit, forming a heterocomplex. This complex is unable to recognize the phage DNA, blocking phage packaging, but specifically binds to the PICI genome, promoting PICI packaging. Our studies reveal the mechanism of action that allows PICI dissemination in nature, introducing a new paradigm in the understanding of the biology of pathogenicity islands and therefore of bacterial pathogen evolution.

tRNA epitranscriptome determines pathogenicity of the opportunistic pathogen Pseudomonas aeruginosa.

The success of bacterial pathogens depends on the coordinated expression of virulence determinants. Regulatory circuits that drive pathogenesis are complex, multilayered, and incompletely understood. Here, we reveal that alterations in tRNA modifications define pathogenic phenotypes in the opportunistic pathogen Pseudomonas aeruginosa. We demonstrate that the enzymatic activity of GidA leads to the introduction of a carboxymethylaminomethyl modification in selected tRNAs. Modifications at the wobble uridine base (cmnm5U34) of the anticodon drives translation of transcripts containing rare codons. Specifically, in P. aeruginosa the presence of GidA-dependent tRNA modifications modulates expression of genes encoding virulence regulators, leading to a cellular proteomic shift toward pathogenic and well-adapted physiological states. Our approach of profiling the consequences of chemical tRNA modifications is general in concept. It provides a paradigm of how environmentally driven tRNA modifications govern gene expression programs and regulate phenotypic outcomes responsible for bacterial adaption to challenging habitats prevailing in the host niche.

An allelic-series rare-variant association test for candidate-gene discovery.

Allelic series are of candidate therapeutic interest because of the existence of a dose-response relationship between the functionality of a gene and the degree or severity of a phenotype. We define an allelic series as a collection of variants in which increasingly deleterious mutations lead to increasingly large phenotypic effects, and we have developed a gene-based rare-variant association test specifically targeted to identifying genes containing allelic series. Building on the well-known burden test and sequence kernel association test (SKAT), we specify a variety of association models covering different genetic architectures and integrate these into a Coding-Variant Allelic-Series Test (COAST). Through extensive simulations, we confirm that COAST maintains the type I error and improves the power when the pattern of coding-variant effect sizes increases monotonically with mutational severity. We applied COAST to identify allelic-series genes for four circulating-lipid traits and five cell-count traits among 145,735 subjects with available whole-exome sequencing data from the UK Biobank. Compared with optimal SKAT (SKAT-O), COAST identified 29% more Bonferroni-significant associations with circulating-lipid traits, on average, and 82% more with cell-count traits. All of the gene-trait associations identified by COAST have corroborating evidence either from rare-variant associations in the full cohort (Genebass, n = 400,000) or from common-variant associations in the GWAS Catalog. In addition to detecting many gene-trait associations present in Genebass by using only a fraction (36.9%) of the sample, COAST detects associations, such as that between ANGPTL4 and triglycerides, that are absent from Genebass but that have clear common-variant support.

De novo mutation hotspots in homologous protein domains identify function-altering mutations in neurodevelopmental disorders.

Variant interpretation remains a major challenge in medical genetics. We developed Meta-Domain HotSpot (MDHS) to identify mutational hotspots across homologous protein domains. We applied MDHS to a dataset of 45,221 de novo mutations (DNMs) from 31,058 individuals with neurodevelopmental disorders (NDDs) and identified three significantly enriched missense DNM hotspots in the ion transport protein domain family (PF00520). The 37 unique missense DNMs that drive enrichment affect 25 genes, 19 of which were previously associated with NDDs. 3D protein structure modeling supports the hypothesis of function-altering effects of these mutations. Hotspot genes have a unique expression pattern in tissue, and we used this pattern alongside in silico predictors and population constraint information to identify candidate NDD-associated genes. We also propose a lenient version of our method, which identifies 32 hotspot positions across 16 different protein domains. These positions are enriched for likely pathogenic variation in clinical databases and DNMs in other genetic disorders.

Prediction of cancer driver genes and mutations: the potential of integrative computational frameworks.

The vast amount of available sequencing data allows the scientific community to explore different genetic alterations that may drive cancer or favor cancer progression. Software developers have proposed a myriad of predictive tools, allowing researchers and clinicians to compare and prioritize driver genes and mutations and their relative pathogenicity. However, there is little consensus on the computational approach or a golden standard for comparison. Hence, benchmarking the different tools depends highly on the input data, indicating that overfitting is still a massive problem. One of the solutions is to limit the scope and usage of specific tools. However, such limitations force researchers to walk on a tightrope between creating and using high-quality tools for a specific purpose and describing the complex alterations driving cancer. While the knowledge of cancer development increases daily, many bioinformatic pipelines rely on single nucleotide variants or alterations in a vacuum without accounting for cellular compartments, mutational burden or disease progression. Even within bioinformatics and computational cancer biology, the research fields work in silos, risking overlooking potential synergies or breakthroughs. Here, we provide an overview of databases and datasets for building or testing predictive cancer driver tools. Furthermore, we introduce predictive tools for driver genes, driver mutations, and the impact of these based on structural analysis. Additionally, we suggest and recommend directions in the field to avoid silo-research, moving towards integrative frameworks.

Two-speed genome evolution drives pathogenicity in fungal pathogens of animals.

The origins and evolution of virulence in amphibian-infecting chytrids Batrachochytrium dendrobatidis (Bd) and Batrachochytrium salamandrivorans (Bsal) are largely unknown. Here, we use deep nanopore sequencing of Bsal and comparative genomics against 21 high-quality genome assemblies that span the fungal Chytridiomycota. We discover that Bsal has the most repeat-rich genome of the Chytridiomycota, comprising 40.9% repetitive elements; this genome has expanded to more than 3× the length of its conspecific Bd, with autonomous and fully functional LTR/Gypsy elements contributing significantly to the expansion. The M36 metalloprotease virulence factors are highly expanded (n = 177) in Bsal, most of which (53%) are flanked by transposable elements, suggesting they have a repeat-associated expansion. We find enrichment upstream of M36 metalloprotease genes of three novel repeat families belonging to the repeat superfamily of LINEs that are implicated with gene copy number variations. Additionally, Bsal has a highly compartmentalized genome architecture, with virulence factors enriched in gene-sparse/repeat-rich compartments, while core conserved genes are enriched in gene-rich/repeat-poor compartments. Genes upregulated during infection are primarily found in the gene-sparse/repeat-rich compartment in both Bd and Bsal. Furthermore, genes with signatures of positive selection in Bd are enriched in repeat-rich regions, suggesting these regions are a cradle for the evolution of chytrid pathogenicity. These are the hallmarks of two-speed genome evolution, and this study provides evidence of two-speed genomes in an animal pathogen, shedding light on the evolution of fungal pathogens of vertebrates driving global declines and extinctions.

Enhanced pathogenicity of Th17 cells due to natalizumab treatment: Implications for MS disease rebound.

After natalizumab (NAT) cessation, some multiple sclerosis (MS) patients experience a severe disease rebound. The rebound pathophysiology is still unclear; however, it has been linked to interleukin-17-producing T-helper (Th17) cells. We demonstrate that during NAT treatment, MCAM+CCR6+Th17 cells gradually acquire a pathogenic profile, including proinflammatory cytokine production, pathogenic transcriptional signatures, brain endothelial barrier impairment, and oligodendrocyte damage via induction of apoptotic pathways. This is accompanied by an increase in Th17 cell frequencies in the cerebrospinal fluid of NAT-treated patients. Notably, Th17 cells derived from NAT-treated patients, who later developed a disease rebound upon treatment cessation, displayed a distinct transcriptional pathogenicity profile associated with altered migratory properties. Accordingly, increased brain infiltration of patient Th17 cells was illustrated in a humanized mouse model and brain histology from a rebound patient. Therefore, peripheral blood-accumulated MCAM+CCR6+Th17 cells might be involved in rebound pathophysiology, and monitoring of changes in Th17 cell pathogenicity in patients before/during NAT treatment cessation might enable rebound risk assessment in the future.

Dimerization of the Alzheimer's disease pathogenic receptor SORLA regulates its association with retromer.

SORL1, the gene encoding the large multidomain SORLA protein, has emerged as only the fourth gene that when mutated can by itself cause Alzheimer's disease (AD), and as a gene reliably linked to both the early- and late-onset forms of the disease. SORLA is known to interact with the endosomal trafficking regulatory complex called retromer in regulating the recycling of endosomal cargo, including the amyloid precursor protein (APP) and the glutamate receptor GluA1. Nevertheless, SORLA's precise structural-functional relationship in endosomal recycling tubules remains unknown. Here, we address these outstanding questions by relying on crystallographic and artificial-intelligence evidence to generate a structural model for how SORLA folds and fits into retromer-positive endosomal tubules, where it is found to dimerize via both SORLA's fibronectin-type-III (3Fn)- and VPS10p-domains. Moreover, we identify a SORLA fragment comprising the 3Fn-, transmembrane, and cytoplasmic domains that has the capacity to form a dimer, and to enhance retromer-dependent recycling of APP by decreasing its amyloidogenic processing. Collectively, these observations generate a model for how SORLA dimer (and possibly polymer) formation can function in stabilizing and enhancing retromer function at endosome tubules. These findings can inform investigation of the many AD-associated SORL1 variants for evidence of pathogenicity and can guide discovery of novel drugs for the disease.

In vivo versus in silico assessment of potentially pathogenic missense variants in human reproductive genes.

Infertility is a heterogeneous condition, with genetic causes thought to underlie a substantial fraction of cases. Genome sequencing is becoming increasingly important for genetic diagnosis of diseases including idiopathic infertility; however, most rare or minor alleles identified in patients are variants of uncertain significance (VUS). Interpreting the functional impacts of VUS is challenging but profoundly important for clinical management and genetic counseling. To determine the consequences of these variants in key fertility genes, we functionally evaluated 11 missense variants in the genes ANKRD31, BRDT, DMC1, EXO1, FKBP6, MCM9, M1AP, MEI1, MSH4 and SEPT12 by generating genome-edited mouse models. Nine variants were classified as deleterious by most functional prediction algorithms, and two disrupted a protein-protein interaction (PPI) in the yeast two hybrid (Y2H) assay. Though these genes are essential for normal meiosis or spermiogenesis in mice, only one variant, observed in the MCM9 gene of a male infertility patient, compromised fertility or gametogenesis in the mouse models. To explore the disconnect between predictions and outcomes, we compared pathogenicity calls of missense variants made by ten widely used algorithms to 1) those annotated in ClinVar and 2) those evaluated in mice. All the algorithms performed poorly in terms of predicting the effects of human missense variants modeled in mice. These studies emphasize caution in the genetic diagnoses of infertile patients based primarily on pathogenicity prediction algorithms and emphasize the need for alternative and efficient in vitro or in vivo functional validation models for more effective and accurate VUS description to either pathogenic or benign categories.

Syntaxin 3 SPI-2 dependent crosstalk facilitates the division of Salmonella containing vacuole.

Intracellular membrane fusion is mediated by membrane-bridging complexes of soluble N-ethylmaleimide-sensitive factor attachment protein receptors (SNAREs). SNARE proteins are one of the key players in vesicular transport. Several reports shed light on intracellular bacteria modulating host SNARE machinery to establish infection successfully. The critical SNAREs in macrophages responsible for phagosome maturation are Syntaxin 3 (STX3) and Syntaxin 4 (STX4). Reports also suggest that Salmonella actively modulates its vacuole membrane composition to escape lysosomal fusion. Salmonella containing vacuole (SCV) harbours recycling endosomal SNARE Syntaxin 12 (STX12). However, the role of host SNAREs in SCV biogenesis and pathogenesis remains unclear. Upon knockdown of STX3, we observed a reduction in bacterial proliferation, which is concomitantly restored upon the overexpression of STX3. Live-cell imaging of Salmonella-infected cells showed that STX3 localises to the SCV membranes and thus might help in the fusion of SCV with intracellular vesicles to acquire membrane for its division. We also found the interaction STX3-SCV was abrogated when we infected with SPI-2 encoded Type 3 secretion system (T3SS) apparatus mutant (STM ∆ssaV) but not with SPI-1 encoded T3SS apparatus mutant (STM ∆invC). These observations were also consistent in the mice model of Salmonella infection. Together, these results shed light on the effector molecules secreted through T3SS encoded by SPI-2, possibly involved in interaction with host SNARE STX3, which is essential to maintain the division of Salmonella in SCV and help to maintain a single bacterium per vacuole.