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

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.

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.

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.

Innovative approach for high-throughput exploiting sex-specific markers in Japanese parrotfish Oplegnathus fasciatus

BACKGROUND

The use of sex-specific molecular markers has become a prominent method in enhancing fish production and economic value, as well as providing a foundation for understanding the complex molecular mechanisms involved in fish sex determination. Over the past decades, research on male and female sex identification has predominantly employed molecular biology methodologies such as restriction fragment length polymorphism, random amplification of polymorphic DNA, simple sequence repeat, and amplified fragment length polymorphism. The emergence of high-throughput sequencing technologies, particularly Illumina, has led to the utilization of single nucleotide polymorphism and insertion/deletion variants as significant molecular markers for investigating sex identification in fish. The advancement of sex-controlled breeding encounters numerous challenges, including the inefficiency of current methods, intricate experimental protocols, high costs of development, elevated rates of false positives, marker instability, and cumbersome field-testing procedures. Nevertheless, the emergence and swift progress of PacBio high-throughput sequencing technology, characterized by its long-read output capabilities, offers novel opportunities to overcome these obstacles.

FINDINGS

Utilizing male/female assembled genome information in conjunction with short-read sequencing data survey and long-read PacBio sequencing data, a catalog of large-segment (>100 bp) insertion/deletion genetic variants was generated through a genome-wide variant site-scanning approach with bidirectional comparisons. The sequence tagging sites were ranked based on the long-read depth of the insertion/deletion site, with markers exhibiting lower long-read depth being considered more effective for large-segment deletion variants. Subsequently, a catalog of bulk primers and simulated PCR for the male/female variant loci was developed, incorporating primer design for the target region and electronic PCR (e-PCR) technology. The Japanese parrotfish (Oplegnathus fasciatus), belonging to the Oplegnathidae family within the Centrarchiformes order, holds significant economic value as a rocky reef fish indigenous to East Asia. The criteria for rapid identification of male and female differences in Japanese parrotfish were established through agarose gel electrophoresis, which revealed 2 amplified bands for males and 1 amplified band for females. A high-throughput identification catalog of sex-specific markers was then constructed using this method, resulting in the identification of 3,639 (2,786 INS/853 DEL, ♀ as reference) and 3,672 (2,876 INS/833 DEL, ♂ as reference) markers in conjunction with 1,021 and 894 high-quality genetic sex identification markers, respectively. Sixteen differential loci were randomly chosen from the catalog for validation, with 11 of them meeting the criteria for male/female distinctions. The implementation of cost-effective and efficient technological processes would facilitate the rapid advancement of genetic breeding through expediting the high-throughput development of sex genetic markers for various species.

CONCLUSIONS

Our study utilized assembled genome information from male and female individuals obtained from PacBio, in addition to data from short-read sequencing data survey and long-read PacBio sequencing data. We extensively employed genome-wide variant site scanning and identification, high-throughput primer design of target regions, and e-PCR batch amplification, along with statistical analysis and ranking of the long-read depth of the variant sites. Through this integrated approach, we successfully compiled a catalog of large insertion/deletion sites (>100 bp) in both male and female Japanese parrotfish.

Antibody Design for the Quantification of Photosynthetic Proteins and Their Isoforms

Antibodies are a valuable research tool, with uses including detection and quantification of specific proteins. By using peptide fragments to raise antibodies, they can be designed to differentiate between structurally similar proteins, or to bind conserved motifs in divergent proteins. Peptide sequence selection and antibody validation are crucial to ensure reliable results from antibody-based experiments. This chapter describes the steps for the identification of peptide sequences to produce protein- or isoform-specific antibodies using recombinant technologies as well as the subsequent validation of such antibodies. The photosynthetic protein Rubisco activase is used as a case study to explain the various steps involved and key aspects to take into consideration.

Dairy phages escape CRISPR defence of Streptococcus thermophilus via the anti-CRISPR AcrIIA3

Bacterial community collapse due to phage infection is a major risk in cheese making processes. As virulent phages are ubiquitous and diverse in milk fermentation factories, the use of phage-resistant lactic acid bacteria (LAB) is essential to obtain high-quality fermented dairy products. The LAB species Streptococcus thermophilus contains two type II-A CRISPR-Cas systems (CRISPR1 and CRISPR3) that can effectively protect against phage infection. However, virulent streptococcal phages carrying anti-CRISPR proteins (ACR) that block the activity of CRISPR-Cas systems have emerged in yogurt and cheese environments. For example, phages carrying AcrIIA5 can impede both CRISPR1 and CRISPR3 systems, while AcrIIA6 stops only CRISPR1. Here, we explore the activity and diversity of a third streptococcal phage anti-CRISPR protein, namely AcrIIA3. We were able to demonstrate that AcrIIA3 is efficiently active against the CRISPR3-Cas system of S. thermophilus. We used AlphaFold2 to infer the structure of AcrIIA3 and we predicted that this new family of functional ACR in virulent streptococcal phages has a new α-helical fold, with no previously identified structural homologs. Because ACR proteins are being explored as modulators in genome editing applications, we also tested AcrIIA3 against SpCas9. We found that AcrIIA3 could block SpCas9 in bacteria but not in human cells. Understanding the diversity and functioning of anti-defence mechanisms will be of importance in the design of long-term stable starter cultures.

Analysis of geometric morphometrics and molecular phylogeny for Anopheles species in the Republic of Korea

Human malaria, transmitted by Anopheles mosquitoes, is the most predominant mosquito-borne disease that is responsible for hundreds of thousands of deaths worldwide each year. In the Republic of Korea (ROK), there are currently several hundred malaria cases annually, mostly near the demilitarized zone (DMZ). Eight species of Anopheles mosquitoes are currently known to be present in the ROK. Similar to other major malaria vectors in Africa and India, it is very challenging to morphologically differentiate Anopheles mosquitoes in the ROK due to their extremely similar morphology. In this study, wing geometric morphometrics (WGM) were used to differentiate the eight Anopheles species collected at six locations near the DMZ, Seoul and Pyeongtaek from April-October 2021. Phylogenetic analysis was also performed using cytochrome c oxidase subunit 1 (COI), internal transcribed spacer 2 (ITS2), and tyrosine hydroxylase (TH) genes for comparison with WGM analysis and to infer evolutionary relationships. The results of cross-validation (overall accuracy = 74.8%) demonstrated that species identification using WGM alone was not possible with a high accuracy for all eight species. While phylogenetic analyses based on the COI region could not clearly distinguish some species, the analysis based on ITS2 and TH was more useful for resolving the phylogenetic correlation of the eight species. Our results may improve Anopheles species identification strategies for effective identification and control of malaria vectors in the ROK.

Matching variants for functional characterization of genetic variants

Rapid and low-cost sequencing, as well as computer analysis, have facilitated the diagnosis of many genetic diseases, resulting in a substantial rise in the number of disease-associated genes. However, genetic diagnosis of many disorders remains problematic due to the lack of interpretation for many genetic variants, especially missenses, the infeasibility of high-throughput experiments on mammals, and the shortcomings of computational prediction technologies. Additionally, the available mutant databases are not well-utilized. Toward this end, we used Caenorhabditis elegans mutant resources to delineate the functions of eight missense variants (V444I, V517D, E610K, L732F, E817K, H873P, R1105K, and G1205E) and two stop codons (W937stop and Q1434stop), including several matching variants (MatchVar) with human in ciliopathy associated IFT-140 (also called CHE-11)//IFT140 (intraflagellar transport protein 140). Moreover, MatchVars carrying C. elegans mutants, including IFT-140(G680S) and IFT-140(P702A) for the human (G704S) (dbSNP: rs150745099) and P726A (dbSNP: rs1057518064 and a conflicting variation) were created using CRISPR/Cas9. IFT140 is a key component of IFT complex A (IFT-A), which is involved in the retrograde transport of IFT along cilia and the entrance of G protein-coupled receptors into cilia. Functional analysis of all 10 variants revealed that P702A and W937stop, but not others phenocopied the ciliary phenotypes (short cilia, IFT accumulations, mislocalization of membrane proteins, and cilia entry of nonciliary proteins) of the IFT-140 null mutant, indicating that both P702A and W937stop are phenotypic in C. elegans. Our functional data offered experimental support for interpreting human variants, by using ready-to-use mutants carrying MatchVars and generating MatchVars with CRISPR/Cas9.

Phylogenetic and ecophysiological novelty of subsurface mercury methylators in mangrove sediments

Mangrove sediment is a crucial component in the global mercury (Hg) cycling and acts as a hotspot for methylmercury (MeHg) production. Early evidence has documented the ubiquity of well-studied Hg methylators in mangrove superficial sediments; however, their diversity and metabolic adaptation in the more anoxic and highly reduced subsurface sediments are lacking. Through MeHg biogeochemical assay and metagenomic sequencing, we found that mangrove subsurface sediments (20-100 cm) showed a less hgcA gene abundance but higher diversity of Hg methylators than superficial sediments (0-20 cm). Regional-scale investigation of mangrove subsurface sediments spanning over 1500 km demonstrated a prevalence and family-level novelty of Hg-methylating microbial lineages (i.e., those affiliated to Anaerolineae, Phycisphaerae, and Desulfobacterales). We proposed the candidate phylum Zixibacteria lineage with sulfate-reducing capacity as a currently understudied Hg methylator across anoxic environments. Unlike other Hg methylators, the Zixibacteria lineage does not use the Wood-Ljungdahl pathway but has unique capabilities of performing methionine synthesis to donate methyl groups. The absence of cobalamin biosynthesis pathway suggests that this Hg-methylating lineage may depend on its syntrophic partners (i.e., Syntrophobacterales members) for energy in subsurface sediments. Our results expand the diversity of subsurface Hg methylators and uncover their unique ecophysiological adaptations in mangrove sediments.

Varying Modes of Selection Among Toxin Families in the Venoms of the Giant Desert Hairy Scorpions (Hadrurus)

Venoms are primarily believed to evolve under strong diversifying selection resulting from persistent coevolution between predator and prey. Recent research has challenged this hypothesis, proposing that venoms from younger venomous lineages (e.g., snakes and cone snails) are governed predominantly by diversifying selection, while venoms from older venomous lineages (e.g., centipedes, scorpions, and spiders) are under stronger purifying selection. However, most research in older lineages has tested selection at more diverse phylogenetic scales. Although these tests are important for evaluating broad macroevolutionary trends underlying venom evolution, they are less equipped to detect species-level evolutionary trends, which likely have large impacts on venom variation seen at more diverse phylogenetic scales. To test for selection among closely related species from an older venomous lineage, we generated high-throughput venom-gland transcriptomes and venom proteomes for four populations of Giant Desert Hairy Scorpions (Hadrurus), including three Hadrurus arizonensis populations and one Hadrurus spadix population. We detected significant episodic and pervasive diversifying selection across a highly abundant toxin family that likely has a major role in venom function ([Formula: see text]KTxs), providing a contrast to the stronger purifying selection identified from other studies on scorpion venoms. Conversely, we detected weak episodic diversifying and/or stronger purifying selection in four toxin families (non-disulfide bridged peptides, phospholipase A2s, scorpine-like antimicrobial peptides, and serine proteases), most of which were less abundant and likely have ancillary functional roles. Finally, although we detected several major toxin families at disproportionate transcriptomic and/or proteomic abundances, we did not identify significant sex-based variation in Hadrurus venoms.

Comparison of contemporary invasive and non-invasive Streptococcus pneumoniae isolates reveals new insights into circulating anti-microbial resistance determinants

Streptococcus pneumoniae is a major human pathogen with a high burden of disease. Non-invasive isolates (those found in non-sterile sites) are thought to be a key source of invasive isolates (those found in sterile sites) and a reservoir of anti-microbial resistance (AMR) determinants. Despite this, pneumococcal surveillance has almost exclusively focused on invasive isolates. We aimed to compare contemporaneous invasive and non-invasive isolate populations to understand how they interact and identify differences in AMR gene distribution. We used a combination of whole-genome sequencing and phenotypic anti-microbial susceptibility testing and a data set of invasive (n = 1,288) and non-invasive (n = 186) pneumococcal isolates, collected in Victoria, Australia, between 2018 and 2022. The non-invasive population had increased levels of antibiotic resistance to multiple classes of antibiotics including beta-lactam antibiotics penicillin and ceftriaxone. We identified genomic intersections between the invasive and non-invasive populations and no distinct phylogenetic clustering of the two populations. However, this analysis revealed sub-populations overrepresented in each population. The sub-populations that had high levels of AMR were overrepresented in the non-invasive population. We determined that WamR-Pneumo was the most accurate in silico tool for predicting resistance to the antibiotics tested. This tool was then used to assess the allelic diversity of the penicillin-binding protein genes, which acquire mutations leading to beta-lactam antibiotic resistance, and found that they were highly conserved (≥80% shared) between the two populations. These findings show the potential of non-invasive isolates to serve as reservoirs of AMR determinants.

The transcriptional regulator EarA and intergenic terminator sequences modulate archaellation in Pyrococcus furiosus

The regulation of archaellation, the formation of archaeal-specific cell appendages called archaella, is crucial for the motility, adhesion, and survival of archaeal organisms. Although the heavily archaellated and highly motile Pyrococcus furiosus is a key model organism for understanding the production and function of archaella in Euryarchaea, the transcriptional regulation of archaellum assembly is so far unknown. Here we show that the transcription factor EarA is the master regulator of the archaellum (arl) operon transcription, which is further modulated by intergenic transcription termination signals. EarA deletion or overexpression strains demonstrate that EarA is essential for archaellation in P. furiosus and governs the degree of archaellation. Providing a single-molecule update on the transcriptional landscape of the arl operon in P. furiosus, we identify sequence motifs for EarA binding upstream of the arl operon and intergenic terminator sequences as critical elements for fine-tuning the expression of the multicistronic arl cluster. Furthermore, transcriptome re-analysis across different Thermococcales species demonstrated a heterogeneous production of major archaellins, suggesting a more diverse composition of archaella than previously recognized. Overall, our study provides novel insights into the transcriptional regulation of archaellation and highlights the essential role of EarA in Pyrococcus furiosus. These findings advance our understanding of the mechanisms governing archaellation and have implications for the functional diversity of archaella.