ggmsa: a visual exploration tool for multiple sequence alignment and associated data

The Myotubularin Related Proteins and the Untapped Interaction Potential of Their Disordered C-Terminal Regions

Intrinsically disordered regions (IDRs) of proteins remain understudied with enigmatic sequence features relevant to their functions. Members of the myotubularin-related protein (MTMR) family contain uncharacterized IDRs. After decades of research on their phosphatase activity, recent work on the C-terminal IDRs of MTMR7 revealed new interactions and important new functions beyond the phosphatase function. Here we take a broader look at the C-terminal domains (CTDs) of 14 human MTMRs and use bioinformatic tools and biophysical methods to ask which other functions may be probable in this protein family. The predictions show that the CTDs are disordered and carry short linear motifs (SLiMs) important for targeting of MTMRs to defined subcellular compartments and implicating them in signaling, phase separation, interaction with diverse proteins, including transcription factors and are of relevance for cancer research and neuroscience. We also present experimental methods to study the CTDs and use them to characterize the coiled coil (CC) domains of MTMR7 and MTMR9. We show homo- and hetero-oligomerization with preference for MTMR7-CC to form dimers, while MTMR9-CC forms trimers. We relate the results to sequence features and make predictions for the structural landscape of other MTMRs. Our work gives a broad insight into the so far unrecognized features and SLiMs in MTMR-CTDs, and provides the basis for more in-depth experimental research on this diverse protein family and understudied IDRs in proteins in general.

Virulence factor discovery identifies associations between the Fic gene family and Fap2+ fusobacteria in colorectal cancer microbiomes

Fusobacterium is a bacterium associated with colorectal cancer (CRC) tumorigenesis, progression, and metastasis. Fap2 is a fusobacteria-specific outer membrane galactose-binding lectin that mediates Fusobacterium adherence to and invasion of CRC tumors. Advances in omics analyses provide an opportunity to profile and identify microbial genomic features that correlate with the cancer-associated bacterial virulence factor Fap2. Here, we analyze genomes of Fusobacterium colon tumor isolates and find that a family of post-translational modification enzymes containing Fic domains is associated with Fap2 positivity in these strains. We demonstrate that Fic family genes expand with the presence of Fap2 in the fusobacterial pangenome. Through comparative genomic analysis, we find that Fap2+ Fusobacteriota are highly enriched with Fic gene families compared to other cancer-associated and human gut microbiome bacterial taxa. Using a global data set of CRC shotgun metagenomes, we show that fusobacterial Fic and Fap2 genes frequently co-occur in the fecal microbiomes of individuals with late-stage CRC. We further characterize specific Fic gene families harbored by Fap2+ Fusobacterium animalis genomes and detect recombination events and elements of horizontal gene transfer via synteny analysis of Fic gene loci. Exposure of a F. animalis strain to a colon adenocarcinoma cell line increases gene expression of fusobacterial Fic and virulence-associated adhesins. Finally, we demonstrate that Fic proteins are synthesized by F. animalis as Fic peptides are detectable in F. animalis monoculture supernatants. Taken together, our study uncovers Fic genes as potential virulence factors in Fap2+ fusobacterial genomes.IMPORTANCEAccumulating data support that bacterial members of the intra-tumoral microbiota critically influence colorectal cancer progression. Yet, relatively little is known about non-adhesin fusobacterial virulence factors that may influence carcinogenesis. Our genomic analysis and expression assays in fusobacteria identify Fic domain-containing genes, well-studied virulence factors in pathogenic bacteria, as potential fusobacterial virulence features. The Fic family proteins that we find are encoded by fusobacteria and expressed by Fusobacterium animalis merit future investigation to assess their roles in colorectal cancer development and progression.

Structural and biophysical characterization of the cytoplasmic domains of HprS kinase and its interactions with the cognate regulator HprR

The HprSR constitutes the bacterial two-component regulatory system engaged by Escherichia coli to reduce the damaging effects of reactive chlorine and oxygen species present in its cytosol. Hypochlorous acid (HOCl) has been shown to be the molecule capable of activating of the HprSR system. HOCl is produced upon pathogen invasion by phagocytic cells of the human innate immune system, particularly neutrophils, to take advantage of its powerful antimicrobial attributes. Therefore, comprehensive studies concerning bacterial sensing and regulatory HprSR system are indispensable in understanding and effectively eliminating pathogens. Here we present the first crystal structure, solved at 1.7 Å resolution, of the HprS cytoplasmic domains arranged as a homodimer. In both protomers, the catalytic ATP-binding domain contains a non-hydrolysable ATP analog coordinated by a magnesium ion. This structure allowed us to provide a detailed characterization of kinase-substrate interaction. Furthermore, the structural data are supported by biophysical studies of kinase interaction with cognate response regulator HprR and substrate ATP. The kinase activity is also assessed in the presence or absence of HprR.

Acaricide resistance mechanisms and host plant responses in the tomato specialist Aculops lycopersici

BACKGROUND

The mite Aculops lycopersici is a major tomato pest with extremely reduced gene families involved in chemoreception and detoxification. How this limited detoxification toolbox affects the evolution of resistance to acaricides in tomato russet mite(s) (TRM) remains enigmatic. Moreover, although a tomato specialist, TRM has been observed on other Solanaceae and Convolvulaceae plant species, raising questions about transcriptional plasticity underlying host exchange.

RESULTS

We identified a field strain with strongly decreased susceptibility to both abamectin and spiromesifen. We detected target-site resistance caused by mutations at conserved positions in two glutamate-gated chloride channels (GluCl), as well as four overexpressed detoxification genes. We then examined transcriptional responses after host shift from tomato to two Solanaceae (potato and black nightshade) and two Convolvulaceae (sweet potato and hedge bindweed) species, as more challenging host plants. Transcriptional responses varied significantly between host plant families, with key differentially expressed genes (DEGs) related to proteolytic, metabolic and detoxification processes. Last, we also identified DEGs encoding for secreted proteins potentially involved in TRM-host plant interactions.

CONCLUSIONS

Despite a limited detoxification toolbox, A. lycopersici might quickly evolve target-site resistance, probably facilitated by strong selection pressure on the genetic variation associated with enormous population size in field settings. Responses to host plant changes include plasticity in genes related to digestion and detoxification, while most responsive genes are of unknown function and remain to be investigated. © 2024 Society of Chemical Industry.

Loss of cped1 does not affect bone and lean tissue in zebrafish

Human genetic studies have nominated cadherin-like and PC-esterase domain-containing 1 (CPED1) as a candidate target gene mediating bone mineral density (BMD) and fracture risk heritability. Recent efforts to define the role of CPED1 in bone in mouse and human models have revealed complex alternative splicing and inconsistent results arising from gene targeting, making its function in bone difficult to interpret. To better understand the role of CPED1 in adult bone mass and morphology, we conducted a comprehensive genetic and phenotypic analysis of cped1 in zebrafish, an emerging model for bone and mineral research. We analyzed two different cped1 mutant lines and performed deep phenotyping to characterize more than 200 measures of adult vertebral, craniofacial, and lean tissue morphology. We also examined alternative splicing of zebrafish cped1 and gene expression in various cell/tissue types. Our studies fail to support an essential role of cped1 in adult zebrafish bone. Specifically, homozygous mutants for both cped1 mutant alleles, which are expected to result in loss-of-function and impact all cped1 isoforms, exhibited no significant differences in the measures examined when compared to their respective wildtype controls, suggesting that cped1 does not significantly contribute to these traits. We identified sequence differences in critical residues of the catalytic triad between the zebrafish and mouse orthologs of CPED1, suggesting that differences in key residues, as well as distinct alternative splicing, could underlie different functions of CPED1 orthologs in the two species. Our studies fail to support a requirement of cped1 in zebrafish bone and lean tissue, adding to evidence that variants at 7q31.31 can act independently of CPED1 to influence BMD and fracture risk.

Combining Directed Evolution with Machine Learning Enables Accurate Genotype-to-Phenotype Predictions

Linking sequence variation to phenotypic effects is critical for efficient exploitation of large genomic datasets. Here we present a novel approach combining directed evolution with protein language modeling to characterize naturally-evolved variants of a rice immune receptor. Using high-throughput directed evolution, we engineered the rice immune receptor Pik-1 to bind and recognize the fungal proteins Avr-PikC and Avr-PikF, which evade detection by currently characterized Pik-1 alleles. A protein language model was fine-tuned on this data to correlate sequence variation with ligand binding behavior. This modeling was then used to characterize Pik-1 variants found in the 3,000 Rice Genomes Project dataset. Two variants scored highly for binding against Avr-PikC, and in vitro analyses confirmed their improved ligand binding over the wild-type Pik-1 receptor. Overall, this machine learning approach identified promising sources of disease resistance in rice and shows potential utility for exploring the phenotypic variation of other proteins of interest.

Characterizing Circulating microRNA Signatures of Type 2 Diabetes Subtypes

Type 2 diabetes (T2D) is a heterogeneous disease influenced by both genetic and environmental factors. Recent studies suggest that T2D subtypes may exhibit distinct gene expression profiles. In this study, we aimed to identify T2D cluster-specific miRNA expression signatures for the previously reported five clinical subtypes that characterize the underlying pathophysiology of long-standing T2D: severe insulin-resistant diabetes (SIRD), severe insulin-deficient diabetes (SIDD), mild age-related diabetes (MARD), mild obesity-related diabetes (MOD), and mild early-onset diabetes (MEOD). We analyzed the circulating microRNAs (miRNAs) in 45 subjects representing the five T2D clusters and 7 non-T2D healthy controls by single-end small RNA sequencing. Bioinformatic analyses identified a total of 430 known circulating miRNAs and 13 previously unreported novel miRNAs. Of these, 71 were upregulated and 37 were downregulated in either controls or individual clusters. Each T2D subtype was associated with a specific dysregulated miRNA profile, distinct from that of healthy controls. Specifically, 3 upregulated miRNAs were unique to SIRD, 1 to MARD, 9 to MOD, and 18 to MEOD. Among the downregulated miRNAs, 11 were specific to SIRD, 9 to SIDD, 2 to MARD, and 1 to MEOD. Our study confirms the heterogeneity of T2D, represented by distinguishable subtypes both clinically and epigenetically and highlights the potential of miRNAs as markers for distinguishing the pathophysiology of T2D subtypes.

Marsilea: an intuitive generalized paradigm for composable visualizations

Biological data visualization is challenged by the growing complexity of datasets. Traditional single-data plots or simple juxtapositions often fail to fully capture dataset intricacies and interrelations. To address this, we introduce "cross-layout," a novel visualization paradigm that integrates multiple plot types in a cross-like structure, with a central main plot surrounded by secondary plots for enhanced contextualization and interrelation insights. We also introduce "Marsilea," a Python-based implementation of cross-layout visualizations, available in both programmatic and web-based interfaces to support users of all experience levels. This paradigm and its implementation offer a customizable, intuitive approach to advance biological data visualization.

An orphan viral genome with unclear evolutionary status sheds light on a distinct lineage of flavi-like viruses infecting plants

Advancements in high-throughput sequencing and associated bioinformatics methods have significantly expanded the RNA virus repertoire, including novel viruses with highly divergent genomes encoding "orphan" proteins that apparently lack homologous sequences. This absence of homologs in routine sequence similarity search complicates their taxonomic classification and raises a fundamental question: Do these orphan viral genomes represent bona ide viruses? In 2022, an orphan viral genome encoding a large polyprotein was identified in alfalfa (Medicago sativa) and thrips (Frankliniella occidentalis), and named Snake River alfalfa virus (SRAV). SRAV was initially proposed as an uncommon flavi-like virus identified in a plant host distantly related to family Flaviviridae. Subsequently, another research group showed its common occurrence in alfalfa but challenged its taxonomic position, suggesting it belongs to the family Endornaviridae. In this study, a large-scale analysis of 77 publicly available small RNA datasets indicates that SRAV could be detected across various tissues and cultivars of alfalfa, and has a broad geographical distribution. Moreover, profiles of the SRAV-derived small interfering RNAs (vsiRNAs) exhibited typical characteristics of viruses in plant hosts. The evolutionary analysis suggests that SRAV represents a unique class of plant-hosted flavi-like viruses with an unusual genome organization and evolutionary status, distinct from previously identified flavi-like viruses documented to infect plants. The latter shows a close evolutionary relationship to flavi-like viruses primarily found in plant-feeding invertebrates and lacks evidence of triggering host RNA interference (RNAi) responses so far. Moreover, mining the transcriptome shotgun assembly (TSA) database identified two novel viral sequences with a similar genome organization and evolutionary status to SRAV. In summary, our study resolves the disagreement regarding the taxonomic status of SRAV and suggests the potential existence of two distinct clades of plant-hosted flavi-like viruses with independent evolutionary origins. Furthermore, our research provides the first evidence of plant-hosted flavi-like viruses triggering the host's RNAi antiviral response. The widespread occurrence of SRAV underscores its potential ecological significance in alfalfa, a crop of substantial economic importance.

Identification of clonally expanded γδ T-cell populations during CAR-T cell therapy

Anti-CD19 Chimeric Antigen Receptor (CAR)-T cell therapies have shown promise for treating B cell malignancies, but the clinical outcome is influenced by both the CAR-T product and the patient's immune system. The role of γδ T cells in the context of CAR-T cell therapy remains poorly understood. This study investigates the transcriptional heterogeneity, clonal expansion and dynamics of γδ T cells in patients undergoing anti-CD19 CAR-T cell therapy. Longitudinal single cell multi-omics analysis was performed on γδ T cells from four patients receiving anti-CD19 CAR-T cell therapy. Single cell RNA-seq, antibody-based protein profiling (AbSeq) and full-length TCRγδ sequences revealed clonally expanded populations displaying plasticity in T cell differentiation, and temporal dynamics of large clones, suggesting ongoing expansion and differentiation. Clonally expanded γδ T cells had heterogeneous gene expression profiles, occupying seven transcriptionally distinct clusters. Analysis of chemokine markers indicated cluster-specific homing tendencies of circulating γδ T cells to peripheral tissues. We found unexpectedly high frequencies of Vδ1 and Vδ3 cells in the blood with distinct gene and protein expression profiles. This analysis provides insights into the dynamic and heterogeneous nature of γδ T cells following anti-CD19 CAR-T cell therapy, contributing valuable information for optimizing CAR-T cell therapies in B cell malignancies.

Granulins rescue inflammation, lysosome dysfunction, lipofuscin, and neuropathology in a mouse model of progranulin deficiency

Progranulin (PGRN) deficiency is linked to neurodegenerative diseases, including frontotemporal dementia (FTD), Alzheimer's disease, and Parkinson's disease. Proper PGRN levels are critical for brain health; however, the function of PGRN is unclear. PGRN is composed of 7.5 repeat domains, called granulins, and processed into granulins inside the lysosome. PGRN is beneficial for neuronal health, but the role of individual granulins is controversial and unclear. We find that the expression of single granulins broadly rescues disease pathology in Grn-/- mice. Adeno-associated virus (AAV)-mediated expression of human granulin-2/F, granulin-4/A, or PGRN in Grn-/- mouse brain ameliorates dysregulated lysosomal proteins and lipids, microgliosis, and lipofuscinosis. Mechanistically, granulins localize to lysosomes in Grn-/- mouse brains or fibroblasts. These data support the hypothesis that PGRN is a precursor to granulins, which share a beneficial function inside the lysosome to maintain lipid and protein homeostasis to prevent neurodegeneration. Thus, granulins are potential therapeutics to treat FTD-GRN and related diseases.

Genetic Variants Affecting FADS2 Enzyme Dynamics and Gene Expression in Cogenetic Oysters with Different PUFA Levels Provide New Tools to Improve Unsaturated Fatty Acids

Long-chain polyunsaturated fatty acids (LC-PUFAs) are crucial for human health and cannot be produced internally. Bivalves, such as oysters, serve as valuable sources of high-quality PUFAs. The enzyme fatty acid desaturase (FADS) plays a key role in the metabolism of LC-PUFAs. In this study, we conducted a thorough genome-wide analysis of the genes belong to the FADS family in Crassostrea gigas and Crassostrea angulata, with the objective of elucidating the function of the FADS2 and investigating the genetic variations that affect PUFA biosynthesis. We identified six FADS genes distributed across four chromosomes, categorized into three subfamilies. The coding region of FADS2 revealed five non-synonymous mutations that were shown to influence protein structure and stability through molecular dynamics simulations. The promoter region of FADS2 contains ten SNPs and three indels significantly correlated with PUFA content. These genetic variations may explain the differences in PUFA levels observed between the two oyster species and could have potential applications in enhancing PUFA content. This study improves the molecular understanding of PUFA metabolism in oysters and presents a potential strategy for selecting oysters with high PUFA levels.

Identification and analysis of major latex protein (MLP) family genes in Rosa chinensis responsive to Botrytis cinerea infection by RNA-seq approaches

Roses (Rosa chinensis) are among the most cherished ornamental plants globally, yet they are highly susceptible to infections by Botrytis cinerea, the causative agent of gray mold disease. Here we inoculated the resistant rose variety 'Yellow Leisure Liness' with B. cinerea to investigate its resistance mechanisms against gray mold disease. Through transcriptome sequencing, we identified 578 differentially expressed genes (DEGs) that were significantly upregulated at 24, 48, and 72 hours post-inoculation, with these genes significantly enriched for three defense response-related GO terms. Further domain analysis of the genes in these GO terms reveal that 21 DEGs contain the Bet v 1 family domain, belonging to the major latex protein (MLP) gene family, suggesting their potential key role in rose disease resistance. Furthermore, we systematically identified 46 RcMLP genes in roses and phylogenetically categorized them into two distinct subfamilies: group I and II. Genomic duplication analysis indicates that tandem duplication is the main driver for the expansion of the RcMLP family, and these genes have undergone by purifying selection. Additionally, detailed analyses of gene structure, motif composition, and promoter regions reveal that RcMLP genes contain numerous stress-responsive elements, with 32 RcMLP genes harboring fungal elicitor/wound-responsive elements. The constructed potential transcription factor regulatory network showed significant enrichment of the ERF transcription factor family in the regulation of RcMLP genes. Gene expression analysis reveal that DEGs are mainly distributed in subfamily II, where four highly expressed genes (RcMLP13, RcMLP28, RcMLP14, and RcMLP27) are identified in a small branch, with their fold change exceeding ten folds and verified by qRT-PCR. In summary, our research results underscore the potential importance of the RcMLP gene family in response to B. cinerea infection and provide comprehensive basis for further function exploration of the MLP gene family in rose resistance to fungal infections.

Dominance of recombinant DWV genomes with changing viral landscapes as revealed in national US honey bee and varroa mite survey

Honey bees are essential pollinators for global agriculture. The viromes of US commercial apiaries and their ectoparasitic mites are poorly characterized at a strain level and there is a need to integrate genomics into pathogen surveillance. We sequenced RNA viromes from 383 adult bees and 173 mites pooled samples from 11 major US beekeeping hubs in 2021, assembling 45 complete and 1702 partial genomes. Protein sequence similarity networks and recombinant genome identification revealed a new viral landscape. Sinaivirus (n = 312), Iflavirus sacbroodi (n = 280), and Iflavirus aladeformis (DWV, n = 135) genomes were common. Recombinant DWV genomes with high nucleotide identity were widespread, and DWV type A master variants were rare, with an indication that RT-PCR surveillance may over-detect type A due to the prevalence of recombinant DWV genomes. Future work should use genomic strategies to avoid misidentification of common honey bee virus genomes and their impact on colony health.

Arabidopsis thaliana argininosuccinate lyase structure uncovers the role of serine as the catalytic base

Arginine is an important amino acid in plants, as it not only plays a structural role and serves as nitrogen storage but is also a precursor for various molecules, including polyamines and proline. Arginine is produced by argininosuccinate lyase (ASL) which catalyzes the cleavage of argininosuccinate to arginine and fumarate. ASL belongs to the fumarate lyase family and while many members of this family were well-characterized, little is known about plant ASLs. Here we present the first crystal structures of ASL from the model plant, Arabidopsis thaliana (AtASL). One of the structures represents the unliganded form of the AtASL homotetramer. The other structure, obtained from a crystal soaked in argininosuccinate, accommodates the substrate or the reaction products in one of four active sites of the AtASL tetramer. Each active site is located at the interface of three neighboring protomers. The AtASL structure with ligands allowed us to analyze the enzyme-substrate and the enzyme-product interactions in detail. Furthermore, based on our analyses, we describe residues of AtASL crucial for catalysis. The structure of AtASL gives the rationale for the open-to-close transition of the GSS mobile loop and indicates the importance of serine 333 from this loop for the enzymatic action of the enzyme. Finally, we supplemented the structural data with the identification of sequence motifs characteristic for ASLs.

Early Detection of Both Pyrenophora teres f. teres and f. maculata in Asymptomatic Barley Leaves Using Digital Droplet PCR (ddPCR)

Efficient early pathogen detection, before symptom apparition, is crucial for optimizing disease management. In barley, the fungal pathogen Pyrenophora teres is the causative agent of net blotch disease, which exists in two forms: P. teres f. sp. teres (Ptt), causing net-form of net blotch (NTNB), and P. teres f. sp. maculata (Ptm), responsible for spot-form of net blotch (STNB). In this study, we developed primers and a TaqMan probe to detect both Ptt and Ptm. A comprehensive k-mer based analysis was performed across a collection of P. teres genomes to identify the conserved regions that had potential as universal genetic markers. These regions were then analyzed for their prevalence and copy number across diverse Moroccan P. teres strains, using both a k-mer analysis for sequence identification and a phylogenetic assessment to establish genetic relatedness. The designed primer-probe set was successfully validated through qPCR, and early disease detection, prior to symptom development, was achieved using ddPCR. The k-mer analysis performed across the available P. teres genomes suggests the potential for these sequences to serve as universal markers for P. teres, transcending environmental variations.

Releasing a sugar brake generates sweeter tomato without yield penalty

In tomato, sugar content is highly correlated with consumer preferences, with most consumers preferring sweeter fruit1-4. However, the sugar content of commercial varieties is generally low, as it is inversely correlated with fruit size, and growers prioritize yield over flavour quality5-7. Here we identified two genes, tomato (Solanum lycopersicum) calcium-dependent protein kinase 27 (SlCDPK27; also known as SlCPK27) and its paralogue SlCDPK26, that control fruit sugar content. They act as sugar brakes by phosphorylating a sucrose synthase, which promotes degradation of the sucrose synthase. Gene-edited SlCDPK27 and SlCDPK26 knockouts increased glucose and fructose contents by up to 30%, enhancing perceived sweetness without fruit weight or yield penalty. Although there are fewer, lighter seeds in the mutants, they exhibit normal germination. Together, these findings provide insight into the regulatory mechanisms controlling fruit sugar accumulation in tomato and offer opportunities to increase sugar content in large-fruited cultivars without sacrificing size and yield.

CLN3 transcript complexity revealed by long-read RNA sequencing analysis

BACKGROUND

Batten disease is a group of rare inherited neurodegenerative diseases. Juvenile CLN3 disease is the most prevalent type, and the most common pathogenic variant shared by most patients is the "1-kb" deletion which removes two internal coding exons (7 and 8) in CLN3. Previously, we identified two transcripts in patient fibroblasts homozygous for the 1-kb deletion: the 'major' and 'minor' transcripts. To understand the full variety of disease transcripts and their role in disease pathogenesis, it is necessary to first investigate CLN3 transcription in "healthy" samples without juvenile CLN3 disease.

METHODS

We leveraged PacBio long-read RNA sequencing datasets from ENCODE to investigate the full range of CLN3 transcripts across various tissues and cell types in human control samples. Then we sought to validate their existence using data from different sources.

RESULTS

We found that a readthrough gene affects the quantification and annotation of CLN3. After taking this into account, we detected over 100 novel CLN3 transcripts, with no dominantly expressed CLN3 transcript. The most abundant transcript has median usage of 42.9%. Surprisingly, the known disease-associated 'major' transcripts are detected. Together, they have median usage of 1.5% across 22 samples. Furthermore, we identified 48 CLN3 ORFs, of which 26 are novel. The predominant ORF that encodes the canonical CLN3 protein isoform has median usage of 66.7%, meaning around one-third of CLN3 transcripts encode protein isoforms with different stretches of amino acids. The same ORFs could be found with alternative UTRs. Moreover, we were able to validate the translational potential of certain transcripts using public mass spectrometry data.

CONCLUSION

Overall, these findings provide valuable insights into the complexity of CLN3 transcription, highlighting the importance of studying both canonical and non-canonical CLN3 protein isoforms as well as the regulatory role of UTRs to fully comprehend the regulation and function(s) of CLN3. This knowledge is essential for investigating the impact of the 1-kb deletion and rare pathogenic variants on CLN3 transcription and disease pathogenesis.

Tracing the Evolutionary History of the Temperature-Sensing Prion-like Domain in EARLY FLOWERING 3 Highlights the Uniqueness of AtELF3

Plants have evolved mechanisms to anticipate and adjust their growth and development in response to environmental changes. Understanding the key regulators of plant performance is crucial to mitigate the negative influence of global climate change on crop production. EARLY FLOWERING 3 (ELF3) is one such regulator playing a critical role in the circadian clock and thermomorphogenesis. In Arabidopsis thaliana, ELF3 contains a prion-like domain (PrLD) that acts as a thermosensor, facilitating liquid-liquid phase separation at high ambient temperatures. To assess the conservation of this function across the plant kingdom, we traced the evolutionary emergence of ELF3, with a focus on the presence of PrLDs. We found that the PrLD, primarily influenced by the length of polyglutamine (polyQ) repeats, is most prominent in Brassicales. Analyzing 319 natural A. thaliana accessions, we confirmed the previously described wide range of polyQ length variation in AtELF3, but found it to be only weakly associated with geographic origin, climate conditions, and classic temperature-responsive phenotypes. Interestingly, similar polyQ length variation was not observed in several other investigated Bassicaceae species. Based on these findings, available prediction tools and limited experimental evidence, we conclude that the emergence of PrLD, and particularly polyQ length variation, is unlikely to be a key driver of environmental adaptation. Instead, it likely adds an additional layer to ELF3's role in thermomorphogenesis in A. thaliana, with its relevance in other species yet to be confirmed.

A master regulator of central carbon metabolism directly activates virulence gene expression in attaching and effacing pathogens

The ability of the attaching and effacing pathogens enterohaemorrhagic Escherichia coli (EHEC) and Citrobacter rodentium to overcome colonisation resistance is reliant on a type 3 secretion system used to intimately attach to the colonic epithelium. This crucial virulence factor is encoded on a pathogenicity island known as the Locus of Enterocyte Effacement (LEE) but its expression is regulated by several core-genome encoded transcription factors. Here, we unveil that the core transcription factor PdhR, traditionally known as a regulator of central metabolism in response to cellular pyruvate levels, is a key activator of the LEE. Through genetic and molecular analyses, we demonstrate that PdhR directly binds to a specific motif within the LEE master regulatory region, thus activating type 3 secretion directly and enhancing host cell adhesion. Deletion of pdhR in EHEC significantly impacted the transcription of hundreds of genes, with pathogenesis and protein secretion emerging as the most affected functional categories. Furthermore, in vivo studies using C. rodentium, a murine model for EHEC infection, revealed that PdhR is essential for effective host colonization and maximal LEE expression within the host. Our findings provide new insights into the complex regulatory networks governing bacterial pathogenesis. This research highlights the intricate relationship between virulence and metabolic processes in attaching and effacing pathogens, demonstrating how core transcriptional regulators can be co-opted to control virulence factor expression in tandem with the cell's essential metabolic circuitry.

Interaction of chikungunya virus glycoproteins with macrophage factors controls virion production

Despite their role as innate sentinels, macrophages can serve as cellular reservoirs of chikungunya virus (CHIKV), a highly-pathogenic arthropod-borne alphavirus that has caused large outbreaks among human populations. Here, with the use of viral chimeras and evolutionary selection analysis, we define CHIKV glycoproteins E1 and E2 as critical for virion production in THP-1 derived human macrophages. Through proteomic analysis and functional validation, we further identify signal peptidase complex subunit 3 (SPCS3) and eukaryotic translation initiation factor 3 subunit K (eIF3k) as E1-binding host proteins with anti-CHIKV activities. We find that E1 residue V220, which has undergone positive selection, is indispensable for CHIKV production in macrophages, as its mutation attenuates E1 interaction with the host restriction factors SPCS3 and eIF3k. Finally, we show that the antiviral activity of eIF3k is translation-independent, and that CHIKV infection promotes eIF3k translocation from the nucleus to the cytoplasm, where it associates with SPCS3. These functions of CHIKV glycoproteins late in the viral life cycle provide a new example of an intracellular evolutionary arms race with host restriction factors, as well as potential targets for therapeutic intervention.

Understanding species-specific and conserved RNA-protein interactions in vivo and in vitro

While evolution is often considered from a DNA- and protein-centric view, RNA-based regulation can also impact gene expression and protein sequences. Here we examine interspecies differences in RNA-protein interactions using the conserved neuronal RNA-binding protein, Unkempt (UNK) as model. We find that roughly half of mRNAs bound in human are also bound in mouse. Unexpectedly, even when transcript-level binding was conserved across species differential motif usage was prevalent. To understand the biochemical basis of UNK-RNA interactions, we reconstitute the human and mouse UNK-RNA interactomes using a high-throughput biochemical assay. We uncover detailed features driving binding, show that in vivo patterns are captured in vitro, find that highly conserved sites are the strongest bound, and associate binding strength with downstream regulation. Furthermore, subtle sequence differences surrounding motifs are key determinants of species-specific binding. We highlight the complex features driving protein-RNA interactions and how these evolve to confer species-specific regulation.

Multiplex PCR Approach for Rapid African Swine Fever Virus Genotyping

African swine fever virus (ASFV) has been spreading through Europe, Asia, and the Caribbean after its introduction in Georgia in 2007 and, due to its particularly high mortality rate, poses a continuous threat to the pig industry. The golden standard to trace back the ASFV is whole genome sequencing, but it is a cost and time-intensive methodology. A more efficient way of tracing the virus is to amplify only specific genomic regions relevant for genotyping. This is mainly accomplished by amplifying single amplicons by PCR followed by Sanger sequencing. To reduce costs and processivity time, we evaluated a multiplex PCR based on the four primer sets routinely used for ASFV genotyping (B646L, E183L, B602L, and intergenic I73R-I329L), which was followed by Nanopore ligation-based amplicon sequencing. We show that with this protocol, we can genotype ASFV DNA originating from different biological matrices and correctly classify multiple genotypes and strains using a single PCR reaction. Further optimization of this method can be accomplished by adding or swapping the primer sets used for amplification based on the needs of a specific country or region, making it a versatile tool that can speed up the processing time and lower the costs of genotyping during ASFV outbreaks.

Comparative Genomics Uncovers the Evolutionary Dynamics of Detoxification and Insecticide Target Genes Across 11 Phlebotomine Sand Flies

Sand flies infect more than 1 million people annually with Leishmania parasites and other bacterial and viral pathogens. Progress in understanding sand fly adaptations to xenobiotics has been hampered by the limited availability of genomic resources. To address this gap, we sequenced, assembled, and annotated the transcriptomes of 11 phlebotomine sand fly species. Subsequently, we leveraged these genomic resources to generate novel evolutionary insights pertaining to their adaptations to xenobiotics, including those contributing to insecticide resistance. Specifically, we annotated over 2,700 sand fly detoxification genes and conducted large-scale phylogenetic comparisons to uncover the evolutionary dynamics of the five major detoxification gene families: cytochrome P450s (CYPs), glutathione-S-transferases (GSTs), UDP-glycosyltransferases (UGTs), carboxyl/cholinesterases (CCEs), and ATP-binding cassette (ABC) transporters. Using this comparative approach, we show that sand flies have evolved diverse CYP and GST gene repertoires, with notable lineage-specific expansions in gene groups evolutionarily related to known xenobiotic metabolizers. Furthermore, we show that sand flies have conserved orthologs of (i) CYP4G genes involved in cuticular hydrocarbon biosynthesis, (ii) ABCB genes involved in xenobiotic toxicity, and (iii) two primary insecticide targets, acetylcholinesterase-1 (Ace1) and voltage gated sodium channel (VGSC). The biological insights and genomic resources produced in this study provide a foundation for generating and testing hypotheses regarding the molecular mechanisms underlying sand fly adaptations to xenobiotics.

Codon Usage Analysis Reveals Distinct Evolutionary Patterns and Host Adaptation Strategies in Duck Hepatitis Virus 1 (DHV-1) Phylogroups

Duck hepatitis virus 1 (DHV-1) is a major threat to the global poultry industry, causing significant economic losses due to high mortality rates in young ducklings. To better understand the evolution and host adaptation strategies of DHV-1, we conducted a comprehensive codon usage analysis of DHV-1 genomes. Our phylogenetic analysis revealed three well-supported DHV-1 phylogroups (Ia, Ib, and II) with distinct genetic diversity patterns. Comparative analyses of the codon usage bias and dinucleotide abundance uncovered a strong preference for A/U-ended codons and a biased pattern of dinucleotide usage in the DHV-1 genome, with CG dinucleotides being extremely underrepresented. Effective number of codons (ENC) analysis indicated a low codon usage bias in the DHV-1 ORF sequences, suggesting adaptation to host codon usage preferences. PR2 bias, ENC plot, and neutrality analyses revealed that both mutation pressure and natural selection influence the codon usage patterns of DHV-1. Notably, the three DHV-1 phylogroups exhibited distinct evolutionary trends, with phylogroups Ia and Ib showing evidence of neutral evolution accompanied by selective pressure, while the phylogroup II evolution was primarily driven by random genetic drift. Comparative analysis of the codon usage indices (CAI, RCDI, and SiD) among the phylogroups highlighted significant differences between subgroups Ia and Ib, suggesting distinct evolutionary pressures or adaptations influencing their codon usage. These findings contribute to our understanding of DHV-1 evolution and host adaptation, with potential implications for the development of effective control measures and vaccines.

Diverse RNA viruses discovered in multiple seagrass species

Seagrasses are marine angiosperms that form highly productive and diverse ecosystems. These ecosystems, however, are declining worldwide. Plant-associated microbes affect critical functions like nutrient uptake and pathogen resistance, which has led to an interest in the seagrass microbiome. However, despite their significant role in plant ecology, viruses have only recently garnered attention in seagrass species. In this study, we produced original data and mined publicly available transcriptomes to advance our understanding of RNA viral diversity in Zostera marina, Zostera muelleri, Zostera japonica, and Cymodocea nodosa. In Z. marina, we present evidence for additional Zostera marina amalgavirus 1 and 2 genotypes, and a complete genome for an alphaendornavirus previously evidenced by an RNA-dependent RNA polymerase gene fragment. In Z. muelleri, we present evidence for a second complete alphaendornavirus and near complete furovirus. Both are novel, and, to the best of our knowledge, this marks the first report of a furovirus infection naturally occurring outside of cereal grasses. In Z. japonica, we discovered genome fragments that belong to a novel strain of cucumber mosaic virus, a prolific pathogen that depends largely on aphid vectoring for host-to-host transmission. Lastly, in C. nodosa, we discovered two contigs that belong to a novel virus in the family Betaflexiviridae. These findings expand our knowledge of viral diversity in seagrasses and provide insight into seagrass viral ecology.

Loss of cped1 does not affect bone and lean tissue in zebrafish

Human genetic studies have nominated Cadherin-like and PC-esterase Domain-containing 1 (CPED1) as a candidate target gene mediating bone mineral density (BMD) and fracture risk heritability. Recent efforts to define the role of CPED1 in bone in mouse and human models have revealed complex alternative splicing and inconsistent results arising from gene targeting, making its function in bone difficult to interpret. To better understand the role of CPED1 in adult bone mass and morphology, we conducted a comprehensive genetic and phenotypic analysis of cped1 in zebrafish, an emerging model for bone and mineral research. We analyzed two different cped1 mutant lines and performed deep phenotyping to characterize more than 200 measures of adult vertebral, craniofacial, and lean tissue morphology. We also examined alternative splicing of zebrafish cped1 and gene expression in various cell/tissue types. Our studies fail to support an essential role of cped1 in adult zebrafish bone. Specifically, homozygous mutants for both cped1 mutant alleles, which are expected to result in loss-of-function and impact all cped1 isoforms, exhibited no significant differences in the measures examined when compared to their respective wildtype controls, suggesting that cped1 does not significantly contribute to these traits. We identified sequence differences in critical residues of the catalytic triad between the zebrafish and mouse orthologs of CPED1, suggesting that differences in key residues, as well as distinct alternative splicing, could underlie different functions of CPED1 orthologs in the two species. Our studies fail to support a requirement of cped1 in zebrafish bone and lean tissue, adding to evidence that variants at 7q31.31 can act independently of CPED1 to influence BMD and fracture risk.

Viral piracy of host RNA phosphatase DUSP11 by avipoxviruses

Proper recognition of viral pathogens is an essential part of the innate immune response. A common viral replicative intermediate and chemical signal that cells use to identify pathogens is the presence of a triphosphorylated 5' end (5'ppp) RNA, which activates the cytosolic RNA sensor RIG-I and initiates downstream antiviral signaling. While 5'pppRNA generated by viral RNA-dependent RNA polymerases (RdRps) can be a potent activator of the immune response, endogenous RNA polymerase III (RNAPIII) transcripts can retain the 5'pppRNA generated during transcription and induce a RIG-I-mediated immune response. We have previously shown that host RNA triphosphatase dual-specificity phosphatase 11 (DUSP11) can act on both host and viral RNAs, altering their levels and reducing their ability to induce RIG-I activation. Our previous work explored how artificially altered DUSP11 can impact immune activation, prompting further exploration into natural contexts of altered DUSP11. Here, we have identified viral DUSP11 homologs (vDUSP11s) present in some avipoxviruses. Consistent with the known functions of endogenous DUSP11, we have shown that expression of vDUSP11s: 1) reduces levels of endogenous RNAPIII transcripts, 2) reduces a cell's sensitivity to 5'pppRNA-mediated immune activation, and 3) restores virus infection defects seen in the absence of DUSP11. Our results identify a virus-relevant context where DUSP11 activity has been co-opted to alter RNA metabolism and influence the outcome of infection.

The New SH3b_T Domain Increases the Structural and Functional Variability Among SH3b-Like CBDs from Staphylococcal Phage Endolysins

Endolysins, proteins encoded by phages to lyse their hosts and release their progeny, have evolved to adapt to the structural features of each host. The endolysins from Staphylococcus-infecting phages typically feature complex architectures with two enzymatically active domains (EADs) and one cell wall-binding domain (CBD) belonging to the bacterial SH3 (SH3b) superfamily. This study focuses on three SH3b-like CBDs from representative staphylococcal phage endolysins (LysRODI, LysC1C and LysIPLA5) that were structurally and functionally characterized. While RODI_CBD and C1C_CBD were assigned to the well-known SH3_5 family, a new family, SH3b_T (PF24246), was identified using the CBD from LysIPLA5 as a model. GFP-fused CBDs were created to assess their differential binding to a collection of staphylococcal strains. IPLA5_CBD showed enhanced binding to Staphylococcus epidermidis, while RODI_CBD and C1C_CBD exhibited distinct binding profiles, with RODI_CBD targeting Staphylococcus aureus specifically and C1C_CBD displaying broad binding. Sequence comparisons suggested that a few differences in key amino acids could be responsible for the latter binding difference. The CBDs modulated the activity spectrum of synthetic EAD-CBD combinations in accordance with the previous binding profiles, but in a manner that was also dependent on the EAD present in the fusion protein. These results serve as a context for the diversity and versatility of SH3b domains in staphylococcal endolysins, providing insights on how (i) the CBDs from this superfamily have diverged to adapt to diverse bacterial ligands in spite of sharing a common fold; and (ii) the evolution of specificity relies on the EAD-CBD combination rather than solely the CBD.

Exploring the role of candidalysin in the pathogenicity of Candida albicans by gene set enrichment analysis and evolutionary dynamics

AIMS

To explore the pathogenic mechanisms of Candida albicans (C. albicans), focusing on its impact on human health, particularly through invasive infections in the gastrointestinal and respiratory tracts.

METHODS

In this study, we evaluated the demographic and clinical profiles of 7 pneumonia patients. Meanwhile, we used Gene Set Enrichment Analysis (GSEA) and Evolutionary Dynamics method to analyze the role of candidalysin in C. albicans pathogenicity.

RESULTS

By analyzing genomic data and conducting biomedical text mining, we identified novel mutation sites in the candidalysin coding gene ECE1-III, shedding light into the genetic diversity within C. albicans strains and their potential implications for antifungal resistance. Our results revealed significant associations between C. albicans and respiratory as well as gastrointestinal diseases, emphasizing the fungus's role in the pathogenesis of these diseases. Additionally, we identified a new mutation site in the C. albicans strain YF2-5, isolated from patients with pneumonia. This mutation may be associated with its heightened pathogenicity.

CONCLUSION

Our research advances the understanding of C. albicans pathogenicity and opens new avenues for developing targeted antifungal therapies. By focusing on the molecular basis of fungal virulence, we aim to contribute to the development of more effective treatment strategies, addressing the challenge of multidrug resistance in invasive fungal infections.

Preparation and Epitope Identification of Monoclonal Antibodies against the NS6 Protein of Porcine Deltacoronavirus (PDCoV)

Porcine deltacoronavirus (PDCoV) is an emerging enteric pathogen that causes substantial economic losses in the swine industry worldwide. The PDCoV NS6 protein is an accessory protein that plays a pivotal role in the viral life cycle and immune evasion. However, the functions of NS6 and its role in PDCoV pathogenesis remain largely unknown. In this study, we prepared a monoclonal antibody (mAb) 5-A11 that specifically recognizes the PDCoV NS6 protein. The mAb 5-A11 exhibited high specificity for PDCoV, with no cross-reactivity with several major porcine pathogenic viruses. Furthermore, the epitope recognized by mAb 5-A11 was precisely mapped to residues 70EYGSIYGKDFI80 of the NS6 protein using Western blot analysis. Notably, this epitope is highly conserved among different PDCoV isolates. Substantial variations were observed when comparing this epitope with the corresponding regions in the NS6 proteins of other δ coronaviruses, suggesting potential differences in the structure, function, and antigenicity of their NS6 proteins. Our findings provide valuable tools and insights for further elucidating the functions of the NS6 protein and its role in PDCoV pathogenesis, as well as for developing diagnostic and therapeutic strategies against PDCoV infection.

CTCF mutation at R567 causes developmental disorders via 3D genome rearrangement and abnormal neurodevelopment

The three-dimensional genome structure organized by CTCF is required for development. Clinically identified mutations in CTCF have been linked to adverse developmental outcomes. Nevertheless, the underlying mechanism remains elusive. In this investigation, we explore the regulatory roles of a clinically relevant R567W point mutation, located within the 11th zinc finger of CTCF, by introducing this mutation into both murine models and human embryonic stem cell-derived cortical organoid models. Mice with homozygous CTCFR567W mutation exhibit growth impediments, resulting in postnatal mortality, and deviations in brain, heart, and lung development at the pathological and single-cell transcriptome levels. This mutation induces premature stem-like cell exhaustion, accelerates the maturation of GABAergic neurons, and disrupts neurodevelopmental and synaptic pathways. Additionally, it specifically hinders CTCF binding to peripheral motifs upstream to the core consensus site, causing alterations in local chromatin structure and gene expression, particularly at the clustered protocadherin locus. Comparative analysis using human cortical organoids mirrors the consequences induced by this mutation. In summary, this study elucidates the influence of the CTCFR567W mutation on human neurodevelopmental disorders, paving the way for potential therapeutic interventions.

Protocol for fast clonal family inference and analysis from large-scale B cell receptor repertoire sequencing data

The expeditious identification and comprehensive analysis of clonal families from extensive B cell receptor (BCR) repertoire sequencing data are imperative for elucidating the intricacies of B cell immune responses. Here, we introduce a computational pipeline designed to swiftly deduce clonal families from bulk BCR heavy-chain sequencing data, accompanied by a suite of functional modules tailored to streamline post-clustering analysis. The outlined methodology encompasses guidelines for software installation, meticulous data preparation, and the systematic inference and analysis of clonal families. For complete details on the use and execution of this protocol, please refer to Wang et al.1.

The substrate quality of CK2 target sites has a determinant role on their function and evolution

Most biological processes are regulated by signaling modules that bind to short linear motifs. For protein kinases, substrates may have full or only partial matches to the kinase recognition motif, a property known as "substrate quality." However, it is not clear whether differences in substrate quality represent neutral variation or if they have functional consequences. We examine this question for the kinase CK2, which has many fundamental functions. We show that optimal CK2 sites are phosphorylated at maximal stoichiometries and found in many conditions, whereas minimal substrates are more weakly phosphorylated and have regulatory functions. Optimal CK2 sites tend to be more conserved, and substrate quality is often tuned by selection. For intermediate sites, increases or decreases in substrate quality may be deleterious, as we demonstrate for a CK2 substrate at the kinetochore. The results together suggest a strong role for substrate quality in phosphosite function and evolution. A record of this paper's transparent peer review process is included in the supplemental information.

BGCFlow: systematic pangenome workflow for the analysis of biosynthetic gene clusters across large genomic datasets

Genome mining is revolutionizing natural products discovery efforts. The rapid increase in available genomes demands comprehensive computational platforms to effectively extract biosynthetic knowledge encoded across bacterial pangenomes. Here, we present BGCFlow, a novel systematic workflow integrating analytics for large-scale genome mining of bacterial pangenomes. BGCFlow incorporates several genome analytics and mining tools grouped into five common stages of analysis such as: (i) data selection, (ii) functional annotation, (iii) phylogenetic analysis, (iv) genome mining, and (v) comparative analysis. Furthermore, BGCFlow provides easy configuration of different projects, parallel distribution, scheduled job monitoring, an interactive database to visualize tables, exploratory Jupyter Notebooks, and customized reports. Here, we demonstrate the application of BGCFlow by investigating the phylogenetic distribution of various biosynthetic gene clusters detected across 42 genomes of the Saccharopolyspora genus, known to produce industrially important secondary/specialized metabolites. The BGCFlow-guided analysis predicted more accurate dereplication of BGCs and guided the targeted comparative analysis of selected RiPPs. The scalable, interoperable, adaptable, re-entrant, and reproducible nature of the BGCFlow will provide an effective novel way to extract the biosynthetic knowledge from the ever-growing genomic datasets of biotechnologically relevant bacterial species.

Development and application of a whole genome amplicon sequencing method for infectious salmon anemia virus (ISAV)

Infectious salmon anemia (ISA) is an infectious disease primarily affecting farmed Atlantic salmon, Salmo salar, which is caused by the ISA virus (ISAV). ISAV belongs to the Orthomyxoviridae family. The disease is a serious condition resulting in reduced fish welfare and high mortality. In this study, we designed an amplicon-based sequencing protocol for whole genome sequencing of ISAV. The method consists of 80 ISAV-specific primers that cover 92% of the virus genome and was designed to be used on an Illumina MiSeq platform. The sequencing accuracy was investigated by comparing sequences with previously published Sanger sequences. The sequences obtained were nearly identical to those obtained by Sanger sequencing, thus demonstrating that sequences produced by this amplicon sequencing protocol had an acceptable accuracy. The amplicon-based sequencing method was used to obtain the whole genome sequence of 12 different ISAV isolates from a small local epidemic in the northern part of Norway. Analysis of the whole genome sequences revealed that segment reassortment took place between some of the isolates and could identify which segments that had been reassorted.

A parasite odyssey: An RNA virus concealed in Toxoplasma gondii

We are entering a 'Platinum Age of Virus Discovery', an era marked by exponential growth in the discovery of virus biodiversity, and driven by advances in metagenomics and computational analysis. In the ecosystem of a human (or any animal) there are more species of viruses than simply those directly infecting the animal cells. Viruses can infect all organisms constituting the microbiome, including bacteria, fungi, and unicellular parasites. Thus the complexity of possible interactions between host, microbe, and viruses is unfathomable. To understand this interaction network we must employ computationally assisted virology as a means of analyzing and interpreting the millions of available samples to make inferences about the ways in which viruses may intersect human health. From a computational viral screen of human neuronal datasets, we identified a novel narnavirus Apocryptovirus odysseus (Ao) which likely infects the neurotropic parasite Toxoplasma gondii. Previously, several parasitic protozoan viruses (PPVs) have been mechanistically established as triggers of host innate responses, and here we present in silico evidence that Ao is a plausible pro-inflammatory factor in human and mouse cells infected by T. gondii. T. gondii infects billions of people worldwide, yet the prognosis of toxoplasmosis disease is highly variable, and PPVs like Ao could function as a hitherto undescribed hypervirulence factor. In a broader screen of over 7.6 million samples, we explored phylogenetically proximal viruses to Ao and discovered nineteen Apocryptovirus species, all found in libraries annotated as vertebrate transcriptome or metatranscriptomes. While samples containing this genus of narnaviruses are derived from sheep, goat, bat, rabbit, chicken, and pigeon samples, the presence of virus is strongly predictive of parasitic Apicomplexa nucleic acid co-occurrence, supporting the fact that Apocryptovirus is a genus of parasite-infecting viruses. This is a computational proof-of-concept study in which we rapidly analyze millions of datasets from which we distilled a mechanistically, ecologically, and phylogenetically refined hypothesis. We predict that this highly diverged Ao RNA virus is biologically a T. gondii infection, and that Ao, and other viruses like it, will modulate this disease which afflicts billions worldwide.

Nudix hydrolase 23 post-translationally regulates carotenoid biosynthesis in plants

Carotenoids are essential for photosynthesis and photoprotection. Plants must evolve multifaceted regulatory mechanisms to control carotenoid biosynthesis. However, the regulatory mechanisms and the regulators conserved among plant species remain elusive. Phytoene synthase (PSY) catalyzes the highly regulated step of carotenogenesis and geranylgeranyl diphosphate synthase (GGPPS) acts as a hub to interact with GGPP-utilizing enzymes for the synthesis of specific downstream isoprenoids. Here, we report a function of Nudix hydrolase 23 (NUDX23), a Nudix domain-containing protein, in post-translational regulation of PSY and GGPPS for carotenoid biosynthesis. NUDX23 expresses highly in Arabidopsis (Arabidopsis thaliana) leaves. Overexpression of NUDX23 significantly increases PSY and GGPPS protein levels and carotenoid production, whereas knockout of NUDX23 dramatically reduces their abundances and carotenoid accumulation in Arabidopsis. NUDX23 regulates carotenoid biosynthesis via direct interactions with PSY and GGPPS in chloroplasts, which enhances PSY and GGPPS protein stability in a large PSY-GGPPS enzyme complex. NUDX23 was found to co-migrate with PSY and GGPPS proteins and to be required for the enzyme complex assembly. Our findings uncover a regulatory mechanism underlying carotenoid biosynthesis in plants and offer promising genetic tools for developing carotenoid-enriched food crops.

Gene expression, evolution, and the genetics of electrosensing in the smalltooth sawfish, Pristis pectinata

Sawfishes (Pristidae) are large, highly threatened rays named for their tooth-studded rostrum, which is used for prey sensing and capture. Of all five species, the smalltooth sawfish, Pristis pectinata, has experienced the greatest decline in range, currently found in only ~20% of its historic range. To better understand the genetic underpinnings of these taxonomically and morphologically unique animals, we collected transcriptomic data from several tissue types, mapped them to the recently completed reference genome, and contrasted the patterns observed with comparable data from other elasmobranchs. Evidence of positive selection was detected in 79 genes in P. pectinata, several of which are involved in growth factor/receptor tyrosine kinase signaling and body symmetry and may be related to the unique morphology of sawfishes. Changes in these genes may impact cellular responses to environmental conditions such as temperature, dissolved oxygen, and salinity. Data acquired also allow for examination of the molecular components of P. pectinata electrosensory systems, which are highly developed in sawfishes and have likely been influential in their evolutionary success.

Long-read sequencing of extrachromosomal circular DNA and genome assembly of a Solanum lycopersicum breeding line revealed active LTR retrotransposons originating from S. Peruvianum L. introgressions

Transposable elements (TEs) are a major force in the evolution of plant genomes. Differences in the transposition activities and landscapes of TEs can vary substantially, even in closely related species. Interspecific hybridization, a widely employed technique in tomato breeding, results in the creation of novel combinations of TEs from distinct species. The implications of this process for TE transposition activity have not been studied in modern cultivars. In this study, we used nanopore sequencing of extrachromosomal circular DNA (eccDNA) and identified two highly active Ty1/Copia LTR retrotransposon families of tomato (Solanum lycopersicum), called Salsa and Ketchup. Elements of these families produce thousands of eccDNAs under controlled conditions and epigenetic stress. EccDNA sequence analysis revealed that the major parts of eccDNA produced by Ketchup and Salsa exhibited low similarity to the S. lycopersicum genomic sequence. To trace the origin of these TEs, whole-genome nanopore sequencing and de novo genome assembly were performed. We found that these TEs occurred in a tomato breeding line via interspecific introgression from S. peruvianum. Our findings collectively show that interspecific introgressions can contribute to both genetic and phenotypic diversity not only by introducing novel genetic variants, but also by importing active transposable elements from other species.

Phylogeny and structural modeling of the transcription factor CsqR (YihW) from Escherichia coli

CsqR (YihW) is a local transcription factor that controls expression of yih genes involved in degradation of sulfoquinovose in Escherichia coli. We recently showed that expression of the respective gene cassette might be regulated by lactose. Here, we explore the phylogenetic and functional traits of CsqR. Phylogenetic analysis revealed that CsqR had a conserved Met25. Western blot demonstrated that CsqR was synthesized in the bacterial cell as two protein forms, 28.5 (CsqR-l) and 26 kDa (CsqR-s), the latter corresponding to start of translation at Met25. CsqR-s was dramatically activated during growth with sulfoquinovose as a sole carbon source, and displaced CsqR-l in the stationary phase during growth on rich medium. Molecular dynamic simulations revealed two possible states of the CsqR-s structure, with the interdomain linker being represented by either a disordered loop or an ɑ-helix. This helix allowed the hinge-like motion of the N-terminal domain resulting in a switch of CsqR-s between two conformational states, "open" and "compact". We then modeled the interaction of both CsqR forms with putative effectors sulfoquinovose, sulforhamnose, sulfoquinovosyl glycerol, and lactose, and revealed that they all preferred the same pocket in CsqR-l, while in CsqR-s there were two possible options dependent on the linker structure.

Streptomyces-based whole-cell biosensors for detecting diverse cell envelope-targeting antibiotics

Cell envelope-targeting antibiotics are potent therapeutic agents against various bacterial infections. The emergence of multiple antibiotic-resistant strains underscores the significance of identifying potent antimicrobials specifically targeting the cell envelope. However, current drug screening approaches are tedious and lack sufficient specificity and sensitivity, warranting the development of more efficient methods. Genetic circuit-based whole-cell biosensors hold great promise for targeted drug discovery from natural products. Here, we performed comparative transcriptomic analysis of Streptomyces coelicolor M1146 exposed to diverse cell envelope-targeting antibiotics, aiming to identify regulatory elements involved in perceiving and responding to these compounds. Differential gene expression analysis revealed significant activation of VanS/R two-component system in response to the glycopeptide class of cell envelope-acting antibiotics. Therefore, we engineered a pair of VanS/R-based biosensors that exhibit functional complementarity and possess exceptional sensitivity and specificity for glycopeptides detection. Additionally, through promoter screening and characterization, we expanded the biosensor's detection range to include various cell envelope-acting antibiotics beyond glycopeptides. Our genetically engineered biosensor exhibits superior performance, including a dynamic range of up to 887-fold for detecting subtle antibiotic concentration changes in a rapid 2-h response time, enabling high-throughput screening of natural product libraries for antimicrobial agents targeting the bacterial cell envelope.

The complex polyploid genome architecture of sugarcane

Sugarcane, the world's most harvested crop by tonnage, has shaped global history, trade and geopolitics, and is currently responsible for 80% of sugar production worldwide1. While traditional sugarcane breeding methods have effectively generated cultivars adapted to new environments and pathogens, sugar yield improvements have recently plateaued2. The cessation of yield gains may be due to limited genetic diversity within breeding populations, long breeding cycles and the complexity of its genome, the latter preventing breeders from taking advantage of the recent explosion of whole-genome sequencing that has benefited many other crops. Thus, modern sugarcane hybrids are the last remaining major crop without a reference-quality genome. Here we take a major step towards advancing sugarcane biotechnology by generating a polyploid reference genome for R570, a typical modern cultivar derived from interspecific hybridization between the domesticated species (Saccharum officinarum) and the wild species (Saccharum spontaneum). In contrast to the existing single haplotype ('monoploid') representation of R570, our 8.7 billion base assembly contains a complete representation of unique DNA sequences across the approximately 12 chromosome copies in this polyploid genome. Using this highly contiguous genome assembly, we filled a previously unsized gap within an R570 physical genetic map to describe the likely causal genes underlying the single-copy Bru1 brown rust resistance locus. This polyploid genome assembly with fine-grain descriptions of genome architecture and molecular targets for biotechnology will help accelerate molecular and transgenic breeding and adaptation of sugarcane to future environmental conditions.

Cystic echinococcosis in cattle and sheep caused by Echinococcus granulosus sensu stricto genotypes G1 and G3 in the USA

BACKGROUND

Endemic domestic dog-ruminant cycles and human cystic echinococcosis caused by Echinococcus granulosus have been sporadically reported in the United States. However, there is a paucity of molecular data describing the genotypes and haplotypes of this important cestode in domestic ruminant hosts.

METHODS

Ninety-four cysts from the lungs and/or livers of slaughtered beef cattle (76 samples), dairy cows (five samples) and sheep (13 samples) were collected from abattoirs in four states of the USA. Samples were genotyped at two mitochondrial loci, cox1 and nad5. Sequences were used to determine species, genotypes and haplotypes using median joining networks and Bayesian phylogenetic analyses. Cyst fertility was assessed in hematoxylin and eosin-stained sections. Additionally, previously reported autochthonous E. granulosus infections in the USA in various hosts were mapped.

RESULTS

Based on cox1 sequences obtained from 94 cysts, 89 (94.7%) were identified as E. granulosus G1/G3, while five (5.3%) were Taenia hydatigena. Taenia hydatigena were only isolated from sheep. Based on nad5 sequences obtained from 89 hydatid cysts, 96.6% and 3.4% belonged to E. granulosus sensu stricto genotypes G1 and G3 respectively. Two haplotypes were found among E. granulosus cox1 sequences, neither of which was geographically unique. Six haplotypes were found among nad5 sequences in genotype G1, of which five were novel, while one haplotype was found in genotype G3. In the concatenated cox1-nad5 dataset, seven haplotypes were identified, of which six were geographically unique. All cysts from cattle were non-fertile. Four cysts from sheep were fertile.

CONCLUSIONS

All genotyped samples belonged to E. granulosus s.s. This is the first study to our knowledge to confirm the presence of genotypes G1 and G3 in domestic cattle and sheep intermediate hosts in the USA and provide data for future diagnostic and epidemiological studies. Sequences have been deposited in GenBank (cox1 sequences: OR398494-OR398496, nad5 sequences: OR400695-OR400702).

Common mitochondrial deletions in RNA-Seq: evaluation of bulk, single-cell, and spatial transcriptomic datasets

Common mitochondrial DNA (mtDNA) deletions are large structural variants in the mitochondrial genome that accumulate in metabolically active tissues with age and have been investigated in various diseases. We applied the Splice-Break2 pipeline (designed for high-throughput quantification of mtDNA deletions) to human RNA-Seq datasets and describe the methodological considerations for evaluating common deletions in bulk, single-cell, and spatial transcriptomics datasets. A robust evaluation of 1570 samples from 14 RNA-Seq studies showed: (i) the abundance of some common deletions detected in PCR-amplified mtDNA correlates with levels observed in RNA-Seq data; (ii) RNA-Seq library preparation method has a strong effect on deletion detection; (iii) deletions had a significant, positive correlation with age in brain and muscle; (iv) deletions were enriched in cortical grey matter, specifically in layers 3 and 5; and (v) brain regions with dopaminergic neurons (i.e., substantia nigra, ventral tegmental area, and caudate nucleus) had remarkable enrichment of common mtDNA deletions.

Comparative analysis of jujube and sour jujube gave insight into their difference in genetic diversity and suitable habitat

Jujube (Ziziphus jujuba var. jujuba Mill.) and sour jujube (Z. jujuba var. spinosa (Bunge) Hu ex H.F.Chow.) are economically, nutritionally, and ecologically significant members of the Rhamnaceae family. Despite their importance, insufficient research on their genetics and habitats has impeded effective conservation and utilization. To address this knowledge gap, we conducted plastome sequencing, integrated distribution data from China, and assessed genetic diversity and suitable habitat. The plastomes of both species exhibited high conservation and low genetic diversity. A new-found 23 bp species-specific Indel in the petL-petG enabled us to develop a rapid Indel-based identification marker for species discrimination. Phylogenetic analysis and dating illuminated their genetic relationship, showing speciation occurred 6.9 million years ago, in a period of dramatic global temperature fluctuations. Substantial variations in suitable climatic conditions were observed, with the mean temperature of the coldest quarter as the primary factor influencing distributions (-3.16°C-12.73°C for jujube and -5.79°C to 4.11°C for sour jujube, suitability exceeding 0.6). Consequently, distinct conservation strategies are warranted due to differences in suitable habitats, with jujube having a broader distribution and sour jujube concentrated in Northern China. In conclusion, disparate habitats and climatic factors necessitate tailored conservation approaches. Comparing genetic diversity and developing rapid species-specific primers will further enhance the sustainable utilization of these valuable species.

In silico characterization of five novel disease-resistance proteins in Oryza sativa sp. japonica against bacterial leaf blight and rice blast diseases

In the current study, gene network analysis revealed five novel disease-resistance proteins against bacterial leaf blight (BB) and rice blast (RB) diseases caused by Xanthomonas oryzae pv. oryzae (Xoo) and Magnaporthe oryzae (M. oryzae), respectively. In silico modeling, refinement, and model quality assessment were performed to predict the best structures of these five proteins and submitted to ModelArchive for future use. An in-silico annotation indicated that the five proteins functioned in signal transduction pathways as kinases, phospholipases, transcription factors, and DNA-modifying enzymes. The proteins were localized in the nucleus and plasma membrane. Phylogenetic analysis showed the evolutionary relation of the five proteins with disease-resistance proteins (XA21, OsTRX1, PLD, and HKD-motif-containing proteins). This indicates similar disease-resistant properties between five unknown proteins and their evolutionary-related proteins. Furthermore, gene expression profiling of these proteins using public microarray data showed their differential expression under Xoo and M. oryzae infection. This study provides an insight into developing disease-resistant rice varieties by predicting novel candidate resistance proteins, which will assist rice breeders in improving crop yield to address future food security through molecular breeding and biotechnology.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03893-5.

Understanding species-specific and conserved RNA-protein interactions in vivo and in vitro

While evolution is often considered from a DNA- and protein-centric view, RNA-based regulation can also impact gene expression and protein sequences. Here we examined interspecies differences in RNA-protein interactions using the conserved neuronal RNA binding protein, Unkempt (UNK) as model. We find that roughly half of mRNAs bound in human are also bound in mouse. Unexpectedly, even when transcript-level binding was conserved across species differential motif usage was prevalent. To understand the biochemical basis of UNK-RNA interactions, we reconstituted the human and mouse UNK-RNA interactomes using a high-throughput biochemical assay. We uncover detailed features driving binding, show that in vivo patterns are captured in vitro, find that highly conserved sites are the strongest bound, and associate binding strength with downstream regulation. Furthermore, subtle sequence differences surrounding motifs are key determinants of species-specific binding. We highlight the complex features driving protein-RNA interactions and how these evolve to confer species-specific regulation.

Genomic and panproteomic analysis of the development of infant immune responses to antigenically-diverse pneumococci

Streptococcus pneumoniae (pneumococcus) is a nasopharyngeal commensal and respiratory pathogen. This study characterises the immunoglobulin G (IgG) repertoire recognising pneumococci from birth to 24 months old (mo) in a prospectively-sampled cohort of 63 children using a panproteome array. IgG levels are highest at birth, due to transplacental transmission of maternal antibodies. The subsequent emergence of responses to individual antigens exhibit distinct kinetics across the cohort. Stable differences in the strength of individuals' responses, correlating with maternal IgG concentrations, are established by 6 mo. By 12 mo, children develop unique antibody profiles that are boosted by re-exposure. However, some proteins only stimulate substantial responses in adults. Integrating genomic data on nasopharyngeal colonisation demonstrates rare pneumococcal antigens can elicit strong IgG levels post-exposure. Quantifying such responses to the diverse core loci (DCL) proteins is complicated by cross-immunity between variants. In particular, the conserved N terminus of DCL protein zinc metalloprotease B provokes the strongest early IgG responses. DCL proteins' ability to inhibit mucosal immunity likely explains continued pneumococcal carriage despite hosts' polyvalent antibody repertoire. Yet higher IgG levels are associated with reduced incidence, and severity, of pneumonia, demonstrating the importance of the heterogeneity in response strength and kinetics across antigens and individuals.

Dissection of the role of a Src homology 3 domain in the evolution of binding preference of paralogous proteins

Protein-protein interactions (PPIs) drive many cellular processes. Some interactions are directed by Src homology 3 (SH3) domains that bind proline-rich motifs on other proteins. The evolution of the binding specificity of SH3 domains is not completely understood, particularly following gene duplication. Paralogous genes accumulate mutations that can modify protein functions and, for SH3 domains, their binding preferences. Here, we examined how the binding of the SH3 domains of 2 paralogous yeast type I myosins, Myo3 and Myo5, evolved following duplication. We found that the paralogs have subtly different SH3-dependent interaction profiles. However, by swapping SH3 domains between the paralogs and characterizing the SH3 domains freed from their protein context, we find that very few of the differences in interactions, if any, depend on the SH3 domains themselves. We used ancestral sequence reconstruction to resurrect the preduplication SH3 domains and examined, moving back in time, how the binding preference changed. Although the most recent ancestor of the 2 domains had a very similar binding preference as the extant ones, older ancestral domains displayed a gradual loss of interaction with the modern interaction partners when inserted in the extant paralogs. Molecular docking and experimental characterization of the free ancestral domains showed that their affinity with the proline motifs is likely not the cause for this loss of binding. Taken together, our results suggest that a SH3 and its host protein could create intramolecular or allosteric interactions essential for the SH3-dependent PPIs, making domains not functionally equivalent even when they have the same binding specificity.