NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins

Transcriptomics of Diphyllatea (CRuMs) from South Pacific crater lakes confirm new cryptic clades

The Diphyllatea (CRuMs) are heterotrophic protists currently divided into three distinct clades (Diphy I-III). Diphy I are biflagellates in the genus Diphylleia, whereas Diphy II and III represent cryptic clades comprising Collodictyon-type quadriflagellates that were recently distinguished based on rRNA gene phylogenies. Here, we isolated Diphyllatea from freshwater crater lakes on two South Pacific islands and generated high-quality transcriptomes from species representing each clade, including the first transcriptomic data from Diphy III. Phylogenomic analyses support the separation of Diphy II and III, while transcriptome completeness highlights the utility of these data for future studies. Lastly, we discuss the biogeography and ecology of Diphyllatea on these remote islands.

Genomics 2 Proteins portal: a resource and discovery tool for linking genetic screening outputs to protein sequences and structures

Recent advances in AI-based methods have revolutionized the field of structural biology. Concomitantly, high-throughput sequencing and functional genomics have generated genetic variants at an unprecedented scale. However, efficient tools and resources are needed to link disparate data types-to 'map' variants onto protein structures, to better understand how the variation causes disease, and thereby design therapeutics. Here we present the Genomics 2 Proteins portal ( https://g2p.broadinstitute.org/ ): a human proteome-wide resource that maps 20,076,998 genetic variants onto 42,413 protein sequences and 77,923 structures, with a comprehensive set of structural and functional features. Additionally, the Genomics 2 Proteins portal allows users to interactively upload protein residue-wise annotations (for example, variants and scores) as well as the protein structure beyond databases to establish the connection between genomics to proteins. The portal serves as an easy-to-use discovery tool for researchers and scientists to hypothesize the structure-function relationship between natural or synthetic variations and their molecular phenotypes.

Analyzing the performance of short-read classification tools on metagenomic samples toward proper diagnosis of diseases

Accurate knowledge of the genome, virus and bacteria that have invaded our bodies is crucial for diagnosing many human diseases. The field of bioinformatics encompasses the complex computational methods required for this purpose. Metagenomics employs next-generation sequencing (NGS) technology to study and identify microbial communities in environmental samples. This technique allows for the measurement of the relative abundance of different microbes. Various tools are available for detecting bacterial species in sequenced metagenomic samples. In this study, we focus on well-known taxonomic classification tools such as MetaPhlAn4, Centrifuge, Kraken2, and Bracken, and evaluate their performance at the species level using synthetic and real datasets. The results indicate that MetaPhlAn4 exhibited high precision in identifying species in the simulated dataset, while Kraken2 had the best area under the precision-recall curve (AUPR) performance. Centrifuge, Kraken2, and Bracken showed accurate estimation of species abundances, unlike MetaPhlAn4, which had a higher L2 distance. In the real dataset analysis with samples from an inflammatory bowel disease (IBD) research, MetaPhlAn4, and Kraken2 had faster execution times, with differences in performance at family and species levels among the tools. Enterobacteriaceae and Pasteurellaceae were highlighted as the most abundant families by Centrifuge, Kraken2, and MetaPhlAn4, with variations in abundance among ulcerative colitis (UC), Crohn's disease (CD), and control non-IBD (CN) groups. Escherichia coli (E. coli) has the highest abundance among Enterobacteriaceae species in the CD and UC groups in comparison with the CN group. Bracken overestimated E. coli abundance, emphasizing result interpretation caution. The findings of this research can assist in selecting the appropriate short-read classifier, thereby aiding in the diagnosis of target diseases.

The robustness of porin-cytochrome gene clusters from Geobacter metallireducens in extracellular electron transfer

To investigate their roles in extracellular electron transfer (EET), the porin-cytochrome (pcc) gene clusters Gmet0825-0828, Gmet0908-0910, and Gmet0911-0913 of the Gram-negative bacterium Geobacter metallireducens were deleted. Failure to delete all pcc gene clusters at the same time suggested their essential roles in extracellular reduction of Fe(III)-citrate by G. metallireducens. Deletion of Gmet0825-0828 had no impact on bacterial reduction of Fe(III)-citrate but diminished bacterial reduction of ferrihydrite and abolished anode reduction and direct interspecies electron transfer (DIET) to Methanosarcina barkeri and Geobacter sulfurreducens. Although it had no impact on the bacterial reduction of Fe(III)-citrate, deletion of Gmet0908-0910 delayed ferrihydrite reduction, abolished anode reduction, and diminished DIET. Deletion of Gmet0911-0913 had little impact on DIET but diminished bacterial reductions of Fe(III)-citrate, ferrihydrite, and anodes. Most importantly, deletions of both Gmet0825-0828 and Gmet0908-0910 restored bacterial reduction of ferrihydrite and anodes and DIET. Enhanced expression of Gmet0911-0913 in this double mutant when grown in coculture with G. sulfurreducens ΔhybLΔfdnG suggested that this cluster might compensate for impaired EET functions of deleting Gmet0825-0828 and Gmet0908-0910. Thus, these pcc gene clusters played essential, distinct, overlapping, and compensatory roles in EET of G. metallireducens that are difficult to characterize as deletion of some clusters affected expression of others. The robustness of these pcc gene clusters enabled G. metallireducens to mediate EET to different acceptors for anaerobic growth even when two of its three pcc gene clusters were inactivated by mutation. The results from this investigation provide new insights into the roles of pcc gene clusters in bacterial EET.

IMPORTANCE

The Gram-negative bacterium Geobacter metallireducens is of environmental and biotechnological significance. Crucial to the unique physiology of G. metallireducens is its extracellular electron transfer (EET) capability. This investigation sheds new light on the robust roles of the three porin-cytochrome (pcc) gene clusters, which are directly involved in EET across the bacterial outer membrane, in the EET of G. metallireducens. In addition to their essential roles, these gene clusters also play distinct, overlapping, and compensatory roles in the EET of G. metallireducens. The distinct roles of the pcc gene clusters enable G. metallireducens to mediate EET to a diverse group of electron acceptors for anaerobic respirations. The overlapping and compensatory roles of the pcc gene clusters enable G. metallireducens to maintain and restore its EET capability for anaerobic growth when one or two of its three pcc gene clusters are deleted from the genome.

RNA-Seq and WGCNA Analyses Reveal Key Regulatory Modules and Genes for Salt Tolerance in Cotton

The problem of soil salinization has seriously hindered agricultural development. Cotton is a pioneering salinity-tolerant crop, so harvesting its key salinity-tolerant genes is important for improving crop salt tolerance. In this study, we analyzed changes in the transcriptome expression profiles of the salt-tolerant cultivar Lu Mian 28 (LM) and the salt-sensitive cultivar Zhong Mian Suo 12 (ZMS) after applying salt stress, and we constructed weighted gene co-expression networks (WGCNA). The results indicated that photosynthesis, amino acid biosynthesis, membrane lipid remodeling, autophagy, and ROS scavenging are key pathways in the salt stress response. Plant-pathogen interactions, plant hormone signal transduction, the mitogen-activated protein kinase (MAPK) signaling pathway, and carotenoid biosynthesis are the regulatory networks associated with these metabolic pathways that confer cotton salt tolerance. The gene-weighted co-expression network was used to screen four modules closely related to traits, identifying 114 transcription factors, including WRKYs, ERFs, NACs, bHLHs, bZIPs, and MYBs, and 11 hub genes. This study provides a reference for acquiring salt-tolerant cotton and abundant genetic resources for molecular breeding.

Temporal stability of the rumen microbiome and its longitudinal associations with performance traits in beef cattle

The rumen microbiome is the focus of a growing body of research, mostly based on investigation of rumen fluid samples collected once from each animal. Exploring the temporal stability of rumen microbiome profiles is imperative, as it enables evaluating the reliability of findings obtained through single-timepoint sampling. We explored the temporal stability of rumen microbiomes considering taxonomic and functional aspects across the 7-month growing-finishing phase spanning 6 timepoints. We identified a temporally stable core microbiome, encompassing 515 microbial genera (e.g., Methanobacterium) and 417 microbial KEGG genes (e.g., K00856-adenosine kinase). The temporally stable core microbiome profiles collected from all timepoints were strongly associated with production traits with substantial economic and environmental impact (e.g., average daily gain, daily feed intake, and methane emissions); 515 microbial genera explained 45-83%, and 417 microbial genes explained 44-83% of their phenotypic variation. Microbiome profiles influenced by the bovine genome explained 54-87% of the genetic variation of bovine traits. Overall, our results provide evidence that the temporally stable core microbiome identified can accurately predict host performance traits at phenotypic and genetic level based on a single timepoint sample taken as early as 7 months prior to slaughter.

Genotyping-by-sequencing-based high-resolution mapping reveals a single candidate gene for the grapevine veraison locus Ver1

Veraison marks the transition from berry growth to berry ripening and is a crucial phenological stage in grapevine (Vitis vinifera): the berries become soft and begin to accumulate sugars, aromatic substances, and, in red cultivars, anthocyanins for pigmentation, while the organic acid levels begin to decrease. These changes determine the potential quality of wine. However, rising global temperatures lead to earlier flowering and ripening, which strongly influence wine quality. Here, we combined genotyping-by-sequencing with a bioinformatics pipeline on ∼150 F1 genotypes derived from a cross between the early ripening variety "Calardis Musqué" and the late-ripening variety "Villard Blanc". Starting from 20,410 haplotype-based markers, we generated a high-density genetic map and performed a quantitative trait locus analysis based on phenotypic datasets evaluated over 20 yrs. Through locus-specific marker enrichment and recombinant screening of ∼1,000 additional genotypes, we refined the originally postulated 5-mb veraison locus, Ver1, on chromosome 16 to only 112 kb, allowing us to pinpoint the ethylene response factor VviERF027 (VCost.v3 gene ID: Vitvi16g00942, CRIBIv1 gene ID: VIT_16s0100g00400) as veraison candidate gene. Furthermore, the early veraison allele could be traced back to a clonal "Pinot" variant first mentioned in the seventeenth century. "Pinot Precoce Noir" passed this allele over "Madeleine Royale" to the maternal grandparent "Bacchus Weiss" and, ultimately, to the maternal parent "Calardis Musqué". Our findings are crucial for ripening time control, thereby improving wine quality, and for breeding grapevines adjusted to climate change scenarios that have a major impact on agro-ecosystems in altering crop plant phenology.

The genetic cause of neurodevelopmental disorders in 30 consanguineous families

Objective

This study aims to clinically and genetically assess 30 unrelated consanguineous Pakistani families from various ethnic backgrounds, all exhibiting features of neurodevelopmental disorders (NDDs).

Methods

We conducted clinical, genetic, biochemical, and molecular analyses on 30 consanguineous families with NDDs enrolled from various regions of Pakistan. The likely molecular causes of primary microcephaly and NDDs were identified. Detailed clinical investigations and molecular diagnoses were performed using whole exome sequencing (WES) of the proband, followed by Sanger sequencing for validation and segregation in the available family members of the affected families.

Results

WES identified likely disease-causing homozygous variants in 30 unrelated consanguineous families. Six families presented newly described variants in known NDD-related genes: ABAT (c.1439 T > G; p.Phe480Cys) [OMIM613163], SLC12A6 (c.2865_2865insT; p.Glu955Asnfs*5) [OMIM 218000], SHANK3 (c.1305-3_1,305-2delTT; p.Gln29-_Gly305del) [OMIM 606232], BCKDK (c.356_356insC; p.Gly119Alafs*24) [OMIM 614923], DDHD2 (c.2065G > T; p.Asp689Tyr) [OMIM 615033], ERCC2 (c.1255G > A; p.Glu419Lys) [OMIM 610756]. Additionally, 12 families had previously reported disease-causing variants associated with different types of NDDs: ATRX (c.109C > T; p.Arg37*) [OMIM 309580], GPR56 [ADGRG1] (c.1423C > T; p.Arg475*) [OMIM 606854], NAGLU (c.1694G > A; p.Arg565Gln) [OMIM 252920], DOLK (c.3G > A; p.Met1Ile) [OMIM 610768], GPT2 (c.815C > T; p.Ser272Leu) [OMIM 616281], DYNC1I2 (c.607 + 1G > A; p.?) [OMIM 618492], FBXL3 (c.885delT; p.Leu295Phefs25*) [OMIM 606220], LINGO1 (c.869G > A; p.Arg290His) [OMIM 618103], and ASPM (c.3978G > A; Trp1326*, c.9557C > G; p.Ser3186*, c.6994C > T; p.Arg2332*) [OMIM 608716]. All the identified variants showed segregation compatible with autosomal recessive inheritance.

Conclusion

In the present study, we observed a high frequency of ASPM variants in the genetic analysis of 30 consanguineous families exhibiting features of NDDs, particularly those associated with autosomal recessive primary microcephaly. These findings contribute to studies on genotype-phenotype correlation, genetic counseling for families, and a deeper understanding of human brain function and development.

Integration of chromosome locations and functional aspects of enhancers and topologically associating domains in knowledge graphs enables versatile queries about gene regulation

Knowledge about transcription factor binding and regulation, target genes, cis-regulatory modules and topologically associating domains is not only defined by functional associations like biological processes or diseases but also has a determinative genome location aspect. Here, we exploit these location and functional aspects together to develop new strategies to enable advanced data querying. Many databases have been developed to provide information about enhancers, but a schema that allows the standardized representation of data, securing interoperability between resources, has been lacking. In this work, we use knowledge graphs for the standardized representation of enhancers and topologically associating domains, together with data about their target genes, transcription factors, location on the human genome, and functional data about diseases and gene ontology annotations. We used this schema to integrate twenty-five enhancer datasets and two domain datasets, creating the most powerful integrative resource in this field to date. The knowledge graphs have been implemented using the Resource Description Framework and integrated within the open-access BioGateway knowledge network, generating a resource that contains an interoperable set of knowledge graphs (enhancers, TADs, genes, proteins, diseases, GO terms, and interactions between domains). We show how advanced queries, which combine functional and location restrictions, can be used to develop new hypotheses about functional aspects of gene expression regulation.

HPC-T-Annotator: an HPC tool for de novo transcriptome assembly annotation

BACKGROUND

The availability of transcriptomic data for species without a reference genome enables the construction of de novo transcriptome assemblies as alternative reference resources from RNA-Seq data. A transcriptome provides direct information about a species' protein-coding genes under specific experimental conditions. The de novo assembly process produces a unigenes file in FASTA format, subsequently targeted for the annotation. Homology-based annotation, a method to infer the function of sequences by estimating similarity with other sequences in a reference database, is a computationally demanding procedure.

RESULTS

To mitigate the computational burden, we introduce HPC-T-Annotator, a tool for de novo transcriptome homology annotation on high performance computing (HPC) infrastructures, designed for straightforward configuration via a Web interface. Once the configuration data are given, the entire parallel computing software for annotation is automatically generated and can be launched on a supercomputer using a simple command line. The output data can then be easily viewed using post-processing utilities in the form of Python notebooks integrated in the proposed software.

CONCLUSIONS

HPC-T-Annotator expedites homology-based annotation in de novo transcriptome assemblies. Its efficient parallelization strategy on HPC infrastructures significantly reduces computational load and execution times, enabling large-scale transcriptome analysis and comparison projects, while its intuitive graphical interface extends accessibility to users without IT skills.

The plasma membrane inner leaflet PI(4,5)P2 is essential for the activation of proton-activated chloride channels

Proton-activated chloride (PAC) channels, ubiquitously expressed in tissues, regulate intracellular Cl- levels and cell death following acidosis. However, molecular mechanisms and signaling pathways involved in PAC channel modulation are largely unknown. Herein, we determine that phosphatidylinositol 4,5-bisphosphate [PI(4,5)P2] of the plasma membrane inner leaflet is essential for the proton activation of PAC channels. PI(4,5)P2 depletion by activating phosphatidylinositol 5-phosphatases or Gq protein-coupled muscarinic receptors substantially inhibits human PAC currents. In excised inside-out patches, PI(4,5)P2 application to the cytoplasmic side increases the currents. Structural simulation reveals that the putative PI(4,5)P2-binding site is localized within the cytosol in resting state but shifts to the cell membrane's inner surface in an activated state and interacts with inner leaflet PI(4,5)P2. Alanine neutralization of basic residues near the membrane-cytosol interface of the transmembrane helice 2 significantly attenuates PAC currents. Overall, our study uncovers a modulatory mechanism of PAC channel through inner membrane PI(4,5)P2.

Unveiling the microbial realm with VEBA 2.0: a modular bioinformatics suite for end-to-end genome-resolved prokaryotic, (micro)eukaryotic and viral multi-omics from either short- or long-read sequencing

The microbiome is a complex community of microorganisms, encompassing prokaryotic (bacterial and archaeal), eukaryotic, and viral entities. This microbial ensemble plays a pivotal role in influencing the health and productivity of diverse ecosystems while shaping the web of life. However, many software suites developed to study microbiomes analyze only the prokaryotic community and provide limited to no support for viruses and microeukaryotes. Previously, we introduced the Viral Eukaryotic Bacterial Archaeal (VEBA) open-source software suite to address this critical gap in microbiome research by extending genome-resolved analysis beyond prokaryotes to encompass the understudied realms of eukaryotes and viruses. Here we present VEBA 2.0 with key updates including a comprehensive clustered microeukaryotic protein database, rapid genome/protein-level clustering, bioprospecting, non-coding/organelle gene modeling, genome-resolved taxonomic/pathway profiling, long-read support, and containerization. We demonstrate VEBA's versatile application through the analysis of diverse case studies including marine water, Siberian permafrost, and white-tailed deer lung tissues with the latter showcasing how to identify integrated viruses. VEBA represents a crucial advancement in microbiome research, offering a powerful and accessible software suite that bridges the gap between genomics and biotechnological solutions.

Transcriptome analysis unravels key pathways and hub genes related to immature fruit abscission in Camellia oleifera

Immature fruit abscission of Camellia oleifera (C. oleifera) is a common problem limiting yield increases. However, the regulatory mechanisms underlying immature fruit abscission in C. oleifera are unclear. In this study, we systematically investigated changes in the morphological, physiological, and gene expression of fruit abscission zones (FAZs) of soon-to-abscise fruits (M2). We found that fruit abscission before ripening mainly occurs during the August abscission stage of 'Huashuo'. At the beginning of this stage, the FAZs of M2 have a marked dent, and the separation layer structures are preliminarily formed. Phytohormone analysis showed that the contents of indole-3-acetic acid (IAA) and jasmonic acid (JA) in the FAZs of M2 were significantly decreased compared with the non-abscised fruits, while the content of trans-zeatin (TZR) was increased. Transcriptome analysis identified differentially expressed genes (DEGs) mainly involved in phytohormone metabolism, including ethylene, auxin, JA, and the cis-zeatin signal transduction pathway. There were also many DEGs involved in cell wall catabolism. Weighted gene co-expression network analysis (WGCNA) further suggested that the transcription factors NAC100 and ERF114 participate in the immature fruit abscission of C. oleifera. This study provides insights into the fruit abscission mechanism of C. oleifera.

In Silico Evaluation of Lawsonia intracellularis Genes Orthologous to Genes Associated with Pathogenesis in Other Intracellular Bacteria

Proliferative enteropathy is an enteric disease caused by the bacterium Lawsonia intracellularis, which affects several species of domestic and wild animals. The mechanisms underlying the mechanisms employed by L. intracellularis to cause host cell proliferation are poorly understood, mostly because this bacterium is extremely difficult to isolate and propagate in vitro. Comparative genomics methods for searching for genes orthologous to genes known to be associated with pathogenesis allow identification of genes potentially involved in pathogenesis by the pathogen of interest. The goal of this study was to carry out in silico research on L. intracellularis genes orthologous to genes required for intracellular invasion and survival present in other pathogenic bacteria, particularly Brucella abortus, B. melitensis, B. suis, Listeria monocytogenes, Mycobacterium tuberculosis, Mycobacterium avium subspecies paratuberculosis, Salmonella enterica, Yersinia pestis, Y. enterocolitica, and Y. pseudotuberculosis. A total of 127 genes associated with invasion and intracellular survival from five known intracellular bacteria were mapped against the predicted proteomes of all L. intracellularis strains publicly available on GenBank, using the OrthoFinder program. A total of 45 L. intracellularis genes were orthologous to genes associated with pathogenesis of other intracellular bacteria. Genes putatively associated with signal the transduction of chemotaxis and cell motility were identified. Genes related to DNA binding and repair were also identified, with some of them supporting a possible association of bacteria with macrophages or inducing pro-inflammatory responses. The homology-based identification of these genes suggests their potential involvement in the virulence and pathogenicity of L. intracellularis, opening avenues for future research and insights into the molecular mechanisms of Lawsonia-elicited proliferative enteropathy.

Genome analysis of Salinimicrobium sp. 3283s, a deep-sea bacterium isolated from the sediments of South China Sea, China

Salinimicrobium sp. 3283s is an aerobic, golden-yellow pigment-producing, Flavobacteriaceae bacterium isolated from the sediments at the depth of 1751 m in the South China Sea. In this study, we present the complete genome sequence of strain 3283s, which only have a single circular chromosome comprising 3,702,683 bp with 41.41% G + C content and no circular plasmid. In total, 3257 protein coding genes, 45 tRNA, 9 rRNA, and 13 sRNA genes were obtained. In terms of the function of gene annotation, strain 3283s was more different from Salinimicrobium oceani J15B91, which was isolated from the South China Sea at a similar depth, and more similar to a Mariana Trench-derived strain Salinimicrobium profundisediminis MT39, which was closer in phylogenetic taxonomic status, suggesting that strain 3283s possesses a stronger potential to adapt to the deep-sea environment. Furthermore, the high- pressure simulations also confirmed that strain 3283s can grow in both 30 MPa and 60 MPa hydrostatic pressure environments, and that it grows better in 30 MPa hydrostatic pressure environments than in 60 MPa hydrostatic pressure environments. In addition, we found a large number of genes in strain 3283s that can promote better adaptation of the bacteria to the low oxygen and high hydrostatic pressure (HHP) environment of the deep sea, such as biosynthetic enzymes of antioxidant pigments, genes encoding cytochromes with enhanced affinity for oxygen, proteins for adaptation to HHP, and genes encoding TonB-dependent transporters in the absence of flagella.

Comparative genome-wide association study on body weight in Chinese native ducks using four models

The Jinling White duck represents a newly developed breed characterized by a rapid growth rate and a superior meat quality, offering significant economic value and research potential; however, the genetic basis underlying their body weight traits remains less understood. Here, we performed whole-genome resequencing for 201 diverse Jinling White male ducks and conducted population genomic analyses, suggesting a rich genetic diversity within the Jinling White duck population. Equipped with our genomic resources, we applied genome-wide association analysis for body weight on birth (BWB), body weight on 1 wk (BW1), body weight on 3 wk (BW3), body weight on 5 wk (BW5) and body weight on 7 wk (BW7) using 4 statistical models. Comparative studies indicated that factored spectrally transformed linear mixed models (FaST-LMM) demonstrated the most superior efficiency, yielding more results with the minimal false positives. We discovered that PUS7, FBXO11, FOXN2, MSH6, and SLC4A4 were associated with BWB. RAG2, and TMEFF2 were candidate genes for BW1, and STARD13, Klotho, ZAR1L are likely candidates for BW3 and BW5. PLXNC1, ATP1A1, CD58, FRYL, OCIAD1, and OCIAD2 were linked to BW7. These findings provide a genetic reference for the selection and breeding of Jinling White ducks, while also deepened our understanding of Growth and development phenotypic in ducks.

Differences in alternative splicing events in the adaptive strategies of two Daphnia galeata genotypes induced by fish kairomones

BACKGROUND

Daphnia galeata is a suitable model organism for investigating predator-induced defense. Genes and pathways exhibiting differential expression between fish kairomone-treated and untreated groups in D. galeata have been identified. However, understanding of the significance of alternative splicing, a crucial process of the regulation of gene expression in eukaryotes, to this mechanism remains limited. This study measured life-history traits and conducted short-read RNA sequencing and long-read isoform sequencing of two Korean D. galeata genotypes (KB1 and KE2) to uncover the genetic mechanism underlying their phenotypic plasticity under predation stress.

RESULTS

KB1 exhibited strategies to enhance fertility and decrease body length when exposed to fish kairomones, while KE2 deployed an adaptive strategy to increase body length. Full-length transcriptomes from KB1 and KE2 yielded 65,736 and 57,437 transcripts, respectively, of which 32 differentially expressed transcripts (DETs) were shared under predation stress across both genotypes. Prominent DETs common to both genotypes were related to energy metabolism and the immune system. Additionally, differential alternative splicing (DAS) events were detected in both genotypes in response to fish kairomones. DAS genes shared between both genotypes may indicate their significant role in the post-transcriptional stress response to fish predation. Calpain-3, involved in digestion and nutrient absorption, was identified as a DAS gene in both genotypes when exposed to fish kairomones. In addition, the gene encoding thymosin beta, which is related to growth, was found to be a statistically significant DAS only in KB1, while that encoding ultraspiracle protein, also associated with growth, was only identified in KE2. Moreover, transcripts encoding proteins such as EGF-like domain-containing protein, vitellogenin fused with superoxide dismutase, and others were identified overlapping between DAS events and DETs and potentially elucidating their association with the observed phenotypic variation in each genotype.

CONCLUSIONS

Our findings highlight the crucial role of alternative splicing in modulating transcriptome landscape under predation stress in D. galeata, emphasizing the requirement for integrating gene expression and splicing analyses in evolutionary adaptation studies.

Genomic insights into the 2022-2023Vibrio cholerae outbreak in Malawi

Malawi experienced its deadliest Vibrio cholerae (Vc) outbreak following devastating cyclones, with >58,000 cases and >1700 deaths reported between March 2022 and May 2023. Here, we use population genomics to investigate the attributes and origin of the Malawi 2022-2023 Vc outbreak isolates. Our results demonstrate the predominance of ST69 clone, also known as the seventh cholera pandemic El Tor (7PET) lineage, expressing O1 Ogawa (~ 80%) serotype followed by Inaba (~ 16%) and sporadic non-O1/non-7PET serogroups (~ 4%). Phylogenetic reconstruction revealed that the Malawi outbreak strains correspond to a recent importation from Asia into Africa (sublineage AFR15). These isolates harboured known antimicrobial resistance and virulence elements, notably the ICEGEN/ICEVchHai1/ICEVchind5 SXT/R391-like integrative conjugative elements and a CTXφ prophage with the ctxB7 genotype compared to historical Malawian Vc isolates. These data suggest that the devastating cyclones coupled with the recent importation of 7PET serogroup O1 strains, may explain the magnitude of the 2022-2023 cholera outbreak in Malawi.

Time-series metagenomics reveals changing protistan ecology of a temperate dimictic lake

BACKGROUND

Protists, single-celled eukaryotic organisms, are critical to food web ecology, contributing to primary productivity and connecting small bacteria and archaea to higher trophic levels. Lake Mendota is a large, eutrophic natural lake that is a Long-Term Ecological Research site and among the world's best-studied freshwater systems. Metagenomic samples have been collected and shotgun sequenced from Lake Mendota for the last 20 years. Here, we analyze this comprehensive time series to infer changes to the structure and function of the protistan community and to hypothesize about their interactions with bacteria.

RESULTS

Based on small subunit rRNA genes extracted from the metagenomes and metagenome-assembled genomes of microeukaryotes, we identify shifts in the eukaryotic phytoplankton community over time, which we predict to be a consequence of reduced zooplankton grazing pressures after the invasion of a invasive predator (the spiny water flea) to the lake. The metagenomic data also reveal the presence of the spiny water flea and the zebra mussel, a second invasive species to Lake Mendota, prior to their visual identification during routine monitoring. Furthermore, we use species co-occurrence and co-abundance analysis to connect the protistan community with bacterial taxa. Correlation analysis suggests that protists and bacteria may interact or respond similarly to environmental conditions. Cryptophytes declined in the second decade of the timeseries, while many alveolate groups (e.g., ciliates and dinoflagellates) and diatoms increased in abundance, changes that have implications for food web efficiency in Lake Mendota.

CONCLUSIONS

We demonstrate that metagenomic sequence-based community analysis can complement existing efforts to monitor protists in Lake Mendota based on microscopy-based count surveys. We observed patterns of seasonal abundance in microeukaryotes in Lake Mendota that corroborated expectations from other systems, including high abundance of cryptophytes in winter and diatoms in fall and spring, but with much higher resolution than previous surveys. Our study identified long-term changes in the abundance of eukaryotic microbes and provided context for the known establishment of an invasive species that catalyzes a trophic cascade involving protists. Our findings are important for decoding potential long-term consequences of human interventions, including invasive species introduction. Video Abstract.

Localisation and function of key axonemal microtubule inner proteins and dynein docking complex members reveal extensive diversity among vertebrate motile cilia

Vertebrate motile cilia are classified as (9+2) or (9+0), based on the presence or absence of the central pair apparatus, respectively. Cryogenic electron microscopy analyses of (9+2) cilia have uncovered an elaborate axonemal protein composition. The extent to which these features are conserved in (9+0) cilia remains unclear. CFAP53, a key axonemal filamentous microtubule inner protein (fMIP) and a centriolar satellites component, is essential for motility of (9+0), but not (9+2) cilia. Here, we show that in (9+2) cilia, CFAP53 functions redundantly with a paralogous fMIP, MNS1. MNS1 localises to ciliary axonemes, and combined loss of both proteins in zebrafish and mice caused severe outer dynein arm loss from (9+2) cilia, significantly affecting their motility. Using immunoprecipitation, we demonstrate that, whereas MNS1 can associate with itself and CFAP53, CFAP53 is unable to self-associate. We also show that additional axonemal dynein-interacting proteins, two outer dynein arm docking (ODAD) complex members, show differential localisation between types of motile cilia. Together, our findings clarify how paralogous fMIPs, CFAP53 and MNS1, function in regulating (9+2) versus (9+0) cilia motility, and further emphasise extensive structural diversity among these organelles.

INDELpred: Improving the prediction and interpretation of indel pathogenicity within the clinical genome

Small insertions and deletions (indels) are critical yet challenging genetic variations with significant clinical implications. However, the identification of pathogenic indels from neutral variants in clinical contexts remains an understudied problem. Here, we developed INDELpred, a machine-learning-based predictive model for discerning pathogenic from benign indels. INDELpred was established based on key features, including allele frequency, indel length, function-based features, and gene-based features. A set of comprehensive evaluation analyses demonstrated that INDELpred exhibited superior performance over competing methods in terms of computational efficiency and prediction accuracy. Importantly, INDELpred highlighted the crucial role of function-based features in identifying pathogenic indels, with a clear interpretability of the features in understanding the disease-causing variants. We envisage INDELpred as a desirable tool for the detection of pathogenic indels within large-scale genomic datasets, thereby enhancing the precision of genetic diagnoses in clinical settings.

Genomic Insights into the Landfill Microbial Community: Denitrifying Activity Supporting One-Carbon Utilization

In spite of the developments in understanding of denitrifying methylotrophy in the recent years, challenges still exist in unravelling the overall biochemistry of nitrate-dependent methane oxidation in novel or poorly characterized/not-yet-cultured bacteria. In the present study, landfill site was mined for novel C1-carbon-metabolizing bacteria which can use nitrate/nitrite as an electron acceptor. A high-throughput rapid plate assay identified three bacterial isolates with eminent ability for nitrate-dependent methane metabolism under anaerobic conditions. Taxonomic identification by whole-genome sequence-based overall genome relatedness indices accurately assigned the isolates AAK_M13, AAK_M29, and AAK_M39 at the species level to Enterobacter cloacae, Bacillus subtilis, and Bacillus halotolerans, respectively. Several genes encoding sub-components involved in alcohol utilization and denitrification pathways, such as adh, fdh, fdo, nar, nir, and nor, were identified in all the genomes. Though no gene clusters encoding MMO/AMO were annotated, sequencing of PCR amplicons revealed similarity with pMMO/AMO gene using translated nucleotide sequence of strains AAK_M29 and AAK_M39, while strain AAK_M13 showed similarity with XRE family transcriptional regulator. This suggests the horizontal gene transfer and/or presence of a truncated version of a housekeeping enzyme encoded by genes exhibiting partial sequence similarity with pMMO genes that mimicked its function at greenhouse gas emission sites. Owing to lack of conclusive evidence for presence of methane metabolism genes in the selected isolates, further experiment was performed to validate their nitrate-dependent methane oxidation capacities. Bacillus subtilis AAK_M29, Bacillus halotolerans AAK_M39, and Enterobacter cloacae AAK_M13 could oxidize 60%, 75%, and 85% of the added methane respectively accompanied by high nitrate reduction (56-62%) thus supporting the correlation between these two activities. The remarkable ability of these isolates for nitrate-dependent methane metabolism has highlighted their role in ecological contribution and biotechnological potential to serve as methane and nitrate sinks in the landfill sites.

Combination of computational techniques and RNAi reveal targets in Anopheles gambiae for malaria vector control

Increasing reports of insecticide resistance continue to hamper the gains of vector control strategies in curbing malaria transmission. This makes identifying new insecticide targets or alternative vector control strategies necessary. CLassifier of Essentiality AcRoss EukaRyote (CLEARER), a leave-one-organism-out cross-validation machine learning classifier for essential genes, was used to predict essential genes in Anopheles gambiae and selected predicted genes experimentally validated. The CLEARER algorithm was trained on six model organisms: Caenorhabditis elegans, Drosophila melanogaster, Homo sapiens, Mus musculus, Saccharomyces cerevisiae and Schizosaccharomyces pombe, and employed to identify essential genes in An. gambiae. Of the 10,426 genes in An. gambiae, 1,946 genes (18.7%) were predicted to be Cellular Essential Genes (CEGs), 1716 (16.5%) to be Organism Essential Genes (OEGs), and 852 genes (8.2%) to be essential as both OEGs and CEGs. RNA interference (RNAi) was used to validate the top three highly expressed non-ribosomal predictions as probable vector control targets, by determining the effect of these genes on the survival of An. gambiae G3 mosquitoes. In addition, the effect of knockdown of arginase (AGAP008783) on Plasmodium berghei infection in mosquitoes was evaluated, an enzyme we computationally inferred earlier to be essential based on chokepoint analysis. Arginase and the top three genes, AGAP007406 (Elongation factor 1-alpha, Elf1), AGAP002076 (Heat shock 70kDa protein 1/8, HSP), AGAP009441 (Elongation factor 2, Elf2), had knockdown efficiencies of 91%, 75%, 63%, and 61%, respectively. While knockdown of HSP or Elf2 significantly reduced longevity of the mosquitoes (p<0.0001) compared to control groups, Elf1 or arginase knockdown had no effect on survival. However, arginase knockdown significantly reduced P. berghei oocytes counts in the midgut of mosquitoes when compared to LacZ-injected controls. The study reveals HSP and Elf2 as important contributors to mosquito survival and arginase as important for parasite development, hence placing them as possible targets for vector control.

Unraveling the genomic secrets of Tritonibacter mobilis AK171: a plant growth-promoting bacterium isolated from Avicennia marina

The scarcity of freshwater resources resulting in a significant yield loss presents a pressing challenge in agriculture. To address this issue, utilizing abundantly available saline water could offer a smart solution. In this study, we demonstrate that the genome sequence rhizosphere bacterium Tritonibacter mobilis AK171, a halophilic marine bacterium recognized for its ability to thrive in saline and waterlogged environments, isolated from mangroves, has the remarkable ability to enable plant growth using saline irrigation. AK171 is characterized as rod-shaped cells, displays agile movement in free-living conditions, and adopts a rosette arrangement in static media. Moreover, The qualitative evaluation of PGP traits showed that AK171 could produce siderophores and IAA but could not solubilize phosphate nor produce hydrolytic enzymes it exhibits a remarkable tolerance to high temperatures and salinity. In this study, we conducted a comprehensive genome sequence analysis of T. mobilis AK171 to unravel the genetic mechanisms underlying its plant growth-promoting abilities in such challenging conditions. Our analysis revealed diverse genes and pathways involved in the bacterium's adaptation to salinity and waterlogging stress. Notably, T. mobilis AK171 exhibited a high level of tolerance to salinity and waterlogging through the activation of stress-responsive genes and the production of specific enzymes and metabolites. Additionally, we identified genes associated with biofilm formation, indicating its potential role in establishing symbiotic relationships with host plants. Furthermore, our analysis unveiled the presence of genes responsible for synthesizing antimicrobial compounds, including tropodithietic acid (TDA), which can effectively control phytopathogens. This genomic insight into T. mobilis AK171 provides valuable information for understanding the molecular basis of plant-microbial interactions in saline and waterlogged environments. It offers potential applications for sustainable agriculture in challenging conditions.

The Scr and Csc pathways for sucrose utilization co-exist in E. coli, but only the Scr pathway is widespread in other Enterobacteriaceae

Most Escherichia coli isolates from humans do not utilize D-sucrose as a substrate for fermentation or growth. Previous work has shown that the Csc pathway allows some E. coli to utilize sucrose for slow growth, and this pathway has been engineered in E. coli W strains to enhance use of sucrose as a feedstock for industrial applications. An alternative sucrose utilization pathway, Scr, was first identified in Klebsiella pneumoniae and has been reported in some E. coli and Salmonella enterica isolates. We show here that the Scr pathway is native to an important subset of E. coli phylogroup B2 lineages that lack the Csc pathway but grow rapidly on sucrose. Laboratory E. coli strains derived from MG1655 (phylogroup A, ST10) are unable to utilize sucrose and lack the scr and csc genes, but a recombinant plasmid-borne scr locus enables rapid growth on and fermentation of sucrose. Genome analyses of Enterobacteriaceae indicate that the scr locus is widespread in other Enterobacteriaceae; including Enterobacter and Klebsiella species, and some Citrobacter and Proteus species. In contrast, the Csc pathway is limited mostly to E. coli, some Shigella species (in which csc loci are rendered non-functional by various mutations), and Citrobacter freundii. The more efficient Scr pathway likely has greater potential than the Csc pathway for bioindustrial applications of E. coli and other Enterobacteriaceae using sucrose as a feedstock.

Quantifying functional redundancy in polysaccharide-degrading prokaryotic communities

BACKGROUND

Functional redundancy (FR) is widely present, but there is no consensus on its formation process and influencing factors. Taxonomically distinct microorganisms possessing genes for the same function in a community lead to within-community FR, and distinct assemblies of microorganisms in different communities playing the same functional roles are termed between-community FR. We proposed two formulas to respectively quantify the degree of functional redundancy within and between communities and analyzed the FR degrees of carbohydrate degradation functions in global environment samples using the genetic information of glycoside hydrolases (GHs) encoded by prokaryotes.

RESULTS

Our results revealed that GHs are each encoded by multiple taxonomically distinct prokaryotes within a community, and the enzyme-encoding prokaryotes are further distinct between almost any community pairs. The within- and between-FR degrees are primarily affected by the alpha and beta community diversities, respectively, and are also affected by environmental factors (e.g., pH, temperature, and salinity). The FR degree of the prokaryotic community is determined by deterministic factors.

CONCLUSIONS

We conclude that the functional redundancy of GHs is a stabilized community characteristic. This study helps to determine the FR formation process and influencing factors and provides new insights into the relationships between prokaryotic community biodiversity and ecosystem functions. Video Abstract.

Whole genome sequencing of uropathogenic E. coli from Ireland reveals diverse resistance mechanisms and strong correlation with phenotypic (EUCAST) susceptibility testing

Urinary tract infections (UTI) caused by uropathogenic Escherichia coli (UPEC) pose a global health concern. Resistance mechanisms, including genetic mutations in antimicrobial target genes, efflux pumps, and drug deactivating enzymes, hinder clinical treatment. These resistance factors often spread through mobile genetic elements. Molecular techniques like whole genome sequencing (WGS), multilocus sequence typing (MLST), and phylotyping help decode bacterial genomes and categorise resistance genes. In this study, we analysed 57 UPEC isolates from different UTI patients following EUCAST guidelines. A selection of 17 representative strains underwent WGS, phylotyping, MLST, and comparative analysis to connect laboratory susceptibility data with predictive genomics based on key resistance genes and chromosomal mutations in antimicrobial targets. Trimethoprim resistance consistently correlated with dfr genes, with six different alleles detected among the isolates. These dfr genes often coexisted with class 1 integrons, with the most common gene cassette combining dfr and aadA. Furthermore, 52.9% of isolates harboured the blaTem-1 gene, rendering resistance to ampicillin and amoxicillin. Ciprofloxacin-resistant strains exhibited mutations in GyrA, GyrB and ParC, plasmid-mediated quinolone resistance genes (qnrb10), and aac(6')-Ib-cr5. Nitrofurantoin resistance in one isolate stemmed from a four amino acid deletion in NfsB. These findings illustrate the varied strategies employed by UPEC to resist antibiotics and the correlation between clinical susceptibility testing and molecular determinants. As molecular testing gains prominence in clinical applications, understanding key resistance determinants becomes crucial for accurate susceptibility testing and guiding effective antimicrobial therapy.

Description of Asgardarchaeum abyssi gen. nov. spec. nov., a novel species within the class Asgardarchaeia and phylum Asgardarchaeota in accordance with the SeqCode

Asgardarchaeota, commonly referred to as Asgard archaea, is a candidatus phylum-rank archaeal clade that includes the closest archaeal relatives of eukaryotes. Despite their prevalence in the scientific literature, the name Asgardarchaeota lacks nomenclatural validation. Here, we describe a novel high-quality metagenome-assembled genome (MAG), AB3033_2TS, proposed to serve as the nomenclatural type for the species Asgardarchaeum abyssiTS according to the rules of the SeqCode. Based on protein content and compositional features, we infer that A. abyssi AB3033_2TS is an acetogenic chemoheterotroph, possibly a facultative lithoautotroph, and is adapted to a thermophilic lifestyle. Utilizing genomes from Asgard archaea, TACK, and Euryarchaea, we perform phylogenomic reconstructions using the GTDB archaeal marker genes, the current reference set for taxonomic classification. Calibrating relative evolutionary divergence (RED) values for Asgardarchaeota using established Thermoproteota lineages in the GTDB r207 reference tree, we establish a robust classification and propose Asgardarchaeum as the type genus for the family Asgardarchaeaceae (fam. nov)., the order Asgardarchaeales (ord. nov.), the class Asgardarchaeia (class. nov.), and the phylum Asgardarchaeota (phyl. nov.). This effort aims to preserve taxonomic congruence in the scientific literature.

Integrated genomics provides insights into the evolution of the polyphosphate accumulation trait of Ca. Accumulibacter

Candidatus Accumulibacter, a prominent polyphosphate-accumulating organism (PAO) in wastewater treatment, plays a crucial role in enhanced biological phosphorus removal (EBPR). The genetic underpinnings of its polyphosphate accumulation capabilities, however, remain largely unknown. Here, we conducted a comprehensive genomic analysis of Ca. Accumulibacter-PAOs and their relatives within the Rhodocyclaceae family, identifying 124 core genes acquired via horizontal gene transfer (HGT) at its least common ancestor. Metatranscriptomic analysis of an enrichment culture of Ca. Accumulibacter revealed active transcription of 44 of these genes during an EBPR cycle, notably including the polyphosphate kinase 2 (PPK2) gene instead of the commonly recognized polyphosphate kinase 1 (PPK1) gene. Intriguingly, the phosphate regulon (Pho) genes showed minimal transcriptions, pointing to a distinctive fact of Pho dysregulation, where PhoU, the phosphate signaling complex protein, was not regulating the high-affinity phosphate transport (Pst) system, resulting in continuous phosphate uptake. To prevent phosphate toxicity, Ca. Accumulibacter utilized the laterally acquired PPK2 to condense phosphate into polyphosphate, resulting in the polyphosphate-accumulating feature. This study provides novel insights into the evolutionary emergence of the polyphosphate-accumulating trait in Ca. Accumulibacter, offering potential advancements in understanding the PAO phenotype in the EBPR process.

Comprehensive Metagenomic Analysis of Veterinary Probiotics in Broiler Chickens

Probiotics are widely used in broiler chickens to support the gut microbiome, gut health, and to reduce the amount of antibiotics used. Despite their benefits, there is concern over their ability to carry and spread antimicrobial resistance genes (ARGs), posing a significant public health risk. This study utilized next-generation sequencing to investigate ARGs in probiotics approved for poultry, focusing on their potential to be transferred via mobile genetic elements such as plasmids and phages. We examined the gut microbiome and resistome changes in 60 broiler chickens over their rearing period, correlating these changes with different probiotic treatments. Specific resistance mechanisms against critically important antibiotics were identified, including genes related to fluoroquinolone resistance and peptide antibiotic resistance. We also found genes with significant relevance to public health (aadK, AAC(6')-Ii) and multiple drug-resistance genes (vmlR, ykkC, ykkD, msrC, clbA, eatAv). Only one phage-encoded gene (dfrA43) was detected, with no evidence of plasmid or mobile genetic element transmission. Additionally, metagenomic analysis of fecal samples showed no significant changes corresponding to time or diet across groups. Our findings highlight the potential risks associated with the use of probiotics in poultry, particularly regarding the carriage of ARGs. It is crucial to conduct further research into the molecular genetics of probiotics to develop strategies that mitigate the risk of resistance gene transfer in agriculture, ensuring the safe and effective use of probiotics in animal husbandry.

Selection and horizontal gene transfer underlie microdiversity-level heterogeneity in resistance gene fate during wastewater treatment

Activated sludge is the centerpiece of biological wastewater treatment, as it facilitates removal of sewage-associated pollutants, fecal bacteria, and pathogens from wastewater through semi-controlled microbial ecology. It has been hypothesized that horizontal gene transfer facilitates the spread of antibiotic resistance genes within the wastewater treatment plant, in part because of the presence of residual antibiotics in sewage. However, there has been surprisingly little evidence to suggest that sewage-associated antibiotics select for resistance at wastewater treatment plants via horizontal gene transfer or otherwise. We addressed the role of sewage-associated antibiotics in promoting antibiotic resistance using lab-scale sequencing batch reactors fed field-collected wastewater, metagenomic sequencing, and our recently developed bioinformatic tool Kairos. Here, we found confirmatory evidence that fluctuating levels of antibiotics in sewage are associated with horizontal gene transfer of antibiotic resistance genes, microbial ecology, and microdiversity-level differences in resistance gene fate in activated sludge.

Minimal and hybrid hydrogenases are active from archaea

Microbial hydrogen (H2) cycling underpins the diversity and functionality of diverse anoxic ecosystems. Among the three evolutionarily distinct hydrogenase superfamilies responsible, [FeFe] hydrogenases were thought to be restricted to bacteria and eukaryotes. Here, we show that anaerobic archaea encode diverse, active, and ancient lineages of [FeFe] hydrogenases through combining analysis of existing and new genomes with extensive biochemical experiments. [FeFe] hydrogenases are encoded by genomes of nine archaeal phyla and expressed by H2-producing Asgard archaeon cultures. We report an ultraminimal hydrogenase in DPANN archaea that binds the catalytic H-cluster and produces H2. Moreover, we identify and characterize remarkable hybrid complexes formed through the fusion of [FeFe] and [NiFe] hydrogenases in ten other archaeal orders. Phylogenetic analysis and structural modeling suggest a deep evolutionary history of hybrid hydrogenases. These findings reveal new metabolic adaptations of archaea, streamlined H2 catalysts for biotechnological development, and a surprisingly intertwined evolutionary history between the two major H2-metabolizing enzymes.

From genome to evolution: investigating type II methylotrophs using a pangenomic analysis

A comprehensive pangenomic approach was employed to analyze the genomes of 75 type II methylotrophs spanning various genera. Our investigation revealed 256 exact core gene families shared by all 75 organisms, emphasizing their crucial role in the survival and adaptability of these organisms. Additionally, we predicted the functionality of 12 hypothetical proteins. The analysis unveiled a diverse array of genes associated with key metabolic pathways, including methane, serine, glyoxylate, and ethylmalonyl-CoA (EMC) metabolic pathways. While all selected organisms possessed essential genes for the serine pathway, Methylooceanibacter marginalis lacked serine hydroxymethyltransferase (SHMT), and Methylobacterium variabile exhibited both isozymes of SHMT, suggesting its potential to utilize a broader range of carbon sources. Notably, Methylobrevis sp. displayed a unique serine-glyoxylate transaminase isozyme not found in other organisms. Only nine organisms featured anaplerotic enzymes (isocitrate lyase and malate synthase) for the glyoxylate pathway, with the rest following the EMC pathway. Methylovirgula sp. 4MZ18 stood out by acquiring genes from both glyoxylate and EMC pathways, and Methylocapsa sp. S129 featured an A-form malate synthase, unlike the G-form found in the remaining organisms. Our findings also revealed distinct phylogenetic relationships and clustering patterns among type II methylotrophs, leading to the proposal of a separate genus for Methylovirgula sp. 4M-Z18 and Methylocapsa sp. S129. This pangenomic study unveils remarkable metabolic diversity, unique gene characteristics, and distinct clustering patterns of type II methylotrophs, providing valuable insights for future carbon sequestration and biotechnological applications.

IMPORTANCE

Methylotrophs have played a significant role in methane-based product production for many years. However, a comprehensive investigation into the diverse genetic architectures across different genera of methylotrophs has been lacking. This study fills this knowledge gap by enhancing our understanding of core hypothetical proteins and unique enzymes involved in methane oxidation, serine, glyoxylate, and ethylmalonyl-CoA pathways. These findings provide a valuable reference for researchers working with other methylotrophic species. Furthermore, this study not only unveils distinctive gene characteristics and phylogenetic relationships but also suggests a reclassification for Methylovirgula sp. 4M-Z18 and Methylocapsa sp. S129 into separate genera due to their unique attributes within their respective genus. Leveraging the synergies among various methylotrophic organisms, the scientific community can potentially optimize metabolite production, increasing the yield of desired end products and overall productivity.

Methylation-Based Characterization of a New IDH2 Mutation in Sinonasal Undifferentiated Carcinoma

Mutations affecting codon 172 of the isocitrate dehydrogenase 2 (IDH2) gene define a subgroup of sinonasal undifferentiated carcinomas (SNUCs) with a relatively favorable prognosis and a globally hypermethylated phenotype. They are also recurrent (along with IDH1 mutations) in gliomas, acute myeloid leukemia, and intrahepatic cholangiocarcinoma. Commonly reported mutations, all associated with aberrant IDH2 enzymatic activity, include R172K, R172S, R172T, R172G, and R172M. We present a case of SNUC with a never-before-described IDH2 mutation, R172A. Our report compares the methylation pattern of our sample to other cases from the Gene Expression Omnibus database. Hierarchical clustering suggests a strong association between our sample and other IDH-mutant SNUCs and a clear distinction between sinonasal normal tissues and tumors. Principal component analysis (PCA), using 100 principal components explaining 94.5% of the variance, showed the position of our sample to be within 1.02 standard deviation of the other IDH-mutant SNUCs. A molecular modeling analysis of the IDH2 R172A versus other R172 variants provides a structural explanation to how they affect the protein active site. Our findings thus suggest that the R172A mutation in IDH2 confers a gain of function similar to other R172 mutations in IDH2, resulting in a similar hypermethylated profile.

Structural and functional analysis of Cyanovirin-N homologs: Carbohydrate binding affinities and antiviral potential of cyanobacterial peptides

Cyanobacteria, a group of photosynthetic prokaryotes, can sinthesize several substances due to their secondary metabolism, with notable properties, such as Cyanovirin-N(CVN), a carbohydrate-binding lectin, that exhibits antiviral activity against several pathogens, due to its ability to bind viral surface carbohydrates such as mannose, thus interfering with the viral entry on the cell. CVN has been described in several cyanobacterial strains and shows biotechnological potential for the development of drugs of pharmaceutical interest. This study focuses on the genomic exploration and characterization of Cyanovirin-N homologs to assess the conservation of carbohydrate-binding affinity within the group. The analysis of their antiviral properties was carried out using bioinformatics tools to study protein models through an in silico pipeline, following the steps of genomic prospection on public databases, homology modeling, docking, molecular dynamics and energetic analysis. Mannose served as the reference ligand, and the lectins' binding affinity with mannose was assessed across Cyanovirin-N homologs. Genomic mining identified 33 cyanobacterial lectin sequences, which underwent structural and functional characterization. The results obtained from this work indicate strong carbohydrate affinity on several homologs, pointing to the conservation of antiviral properties alongside the group. However, this affinity was not uniformly distributed among sequences, exhibiting significant heterogeneity in binding site residues, suggesting potential multi-ligand binding capabilities on the Cyanovirin-N homologs group. Studies focused on the properties involved in these molecules and the investigation of the genetic diversity of Cyanovirin-N homologs could provide valuable insights into the discovery of new drug candidates, harvesting the potential of bioinformatics for large-scale functional and structural analysis.

Genome analysis of Rossellomorea sp. y25, a deep sea bacterium isolated from the sediments of South China Sea

Rossellomorea sp. y25, a putative new species of yellow pigment-producing, aerobic and chemoheterotrophic bacterium belonging to the family Bacillaceae, was isolated from the sediments at the depth of 1829 m in the South China Sea. In this study, we present the complete genome sequences of strain y25, which consisted of only one circular chromosome with 4,633,006 bp and the content of G + C was 41.76%. A total of 4466 CDSs, 106 tRNA, 33 rRNA, and 101 sRNA genes were obtained. Genomic analysis of strain y25 showed that it has the ability to produce antioxidant carotenoids and a large number of heavy metal resistance genes, such as arsenic, cadmium and zinc. In addition, strain y25 contains a prophage that may contribute to host protection against lysis by related Bacillus-like phages. This is the first report of genome-wide information on a bacterium of the genus Rossellomorea isolated from the deep sea, providing insights into how microorganisms of this genus adapt to deep-sea environments.

Mechanisms of acquisition of the vanA operon among vancomycin-resistant Staphylococcus aureus genomes: The tip of the iceberg?

BACKGROUND

Vancomycin is frequently used as a last line of defence against infections due to multidrug-resistant Staphylococcus aureus (S. aureus). A recent finding described the acquisition of vancomycin-resistant S. aureus strains by the integration of an enterococcal plasmid containing the vanA operon into the S. aureus chromosome via homologous recombination involving a specific integration site called locus L2.

METHODS

To characterise all mechanisms of acquisition of vanA, this study analysed the 15 706 S. aureus genomes to look for vanA and described its genetic environment.

RESULTS

A complete vanA operon was found in 25 S. aureus strains isolated from 12 patients, including nine co-isolated with vancomycin-resistant Enterococcus strains. VanA was found within transposon Tn1546-like elements on 17 plasmids and eight chromosomes. VanA might be acquired through conjugation of enterococcal and staphylococcal plasmids, transposition of Tn1546 carrying vanA and plasmid integration into the chromosome. Further, L2 was detected in 2087 genomes (13.3%) of S. aureus strains across different continents. Six potential chromosomal hotspots for integration of the entire vanA-containing enterococcal plasmid were identified by homologous recombination via L2.

CONCLUSIONS

These findings suggest that the recently described scenario in a New York patient could be reproduced anywhere. Surveillance of this possibility is mandatory, especially in patients with vancomycin-resistant Enterococcus infection or colonisation.

Immunoinformatic prediction to identify Staphylococcus aureus peptides that bind to CD8+ T-cells as potential vaccine candidates

Background and Aim

Staphylococcus aureus, with its diverse virulence factors and immune response evasion mechanisms, presents a formidable challenge as an opportunistic pathogen. Developing an effective vaccine against S. aureus has proven elusive despite extensive efforts. Autologous Staphylococcus lysate (ASL) treatment has proven effective in triggering an immune response against bovine mastitis. Peptides that stimulate the immune response can be the subject of further research. The study aimed to use immunoinformatics tools to identify epitopes on S. aureus surface and secretory proteins that can bind to major histocompatibility complex class I (MHC I) and CD8+ T-cells. This method aids in discovering prospective vaccine candidates and elucidating the rationale behind ASL therapy's efficacy.

Materials and Methods

Proteins were identified using both literature search and the National Center for Biotechnology Information search engine Entrez. Self and non-self peptides, allergenicity predictions, epitope locations, and physicochemical characteristics were determined using sequence alignment, AllerTOP, SVMTriP, and Protein-Sol tools. Hex was employed for simulating the docking interactions between S. aureus proteins and the MHC I + CD8+ T-cells complex. The binding sites of S. aureus proteins were assessed using Computer Atlas of Surface Topography of Proteins (CASTp) while docked with MHC I and CD8+ T-cells.

Results

Nine potential S. aureus peptides and their corresponding epitopes were identified in this study, stimulating cytotoxic T-cell mediated immunity. The peptides were analyzed for similarity with self-antigens and allergenicity. 1d20, 2noj, 1n67, 1nu7, 1amx, and 2b71, non-self and stable, are potential elicitors of the cytotoxic T-cell response. The energy values from docking simulations of peptide-MHC I complexes with the CD8+ and T-cell receptor (TCR) indicate the stability and strength of the formed complexes. These peptides - 2noj, 1d20, 1n67, 2b71, 1nu7, 1yn3, 1amx, 2gi9, and 1edk - demonstrated robust MHC I binding, as evidenced by their low binding energies. Peptide 2gi9 exhibited the lowest energy value, followed by 2noj, 1nu7, 1n67, and 1d20, when docked with MHC I and CD8 + TCR, suggesting a highly stable complex. CASTp analysis indicated substantial binding pockets in the docked complexes, with peptide 1d20 showing the highest values for area and volume, suggesting its potential as an effective elicitor of immunological responses. These peptides - 2noj, 2gi9, 1d20, and 1n67 - stand out for vaccine development and T-cell activation against S. aureus.

Conclusion

This study sheds light on the design and development of S. aureus vaccines, highlighting the significance of employing computational methods in conjunction with experimental verification. The significance of T-cell responses in combating S. aureus infections is emphasized by this study. More experiments are needed to confirm the effectiveness of these vaccine candidates and discover their possible medical uses.

Genome characterisation and comparative analysis of Schaalia dentiphila sp. nov. and its subspecies, S. dentiphila subsp. denticola subsp. nov., from the human oral cavity

BACKGROUND

Schaalia species are primarily found among the oral microbiota of humans and other animals. They have been associated with various infections through their involvement in biofilm formation, modulation of host responses, and interaction with other microorganisms. In this study, two strains previously indicated as Actinomyces spp. were found to be novel members of the genus Schaalia based on their whole genome sequences.

RESULTS

Whole-genome sequencing revealed both strains with a genome size of 2.3 Mbp and GC contents of 65.5%. Phylogenetics analysis for taxonomic placement revealed strains NCTC 9931 and C24 as distinct species within the genus Schaalia. Overall genome-relatedness indices including digital DNA-DNA hybridization (dDDH), and average nucleotide/amino acid identity (ANI/AAI) confirmed both strains as distinct species, with values below the species boundary thresholds (dDDH < 70%, and ANI and AAI < 95%) when compared to nearest type strain Schaalia odontolytica NCTC 9935 T. Pangenome and orthologous analyses highlighted their differences in gene properties and biological functions compared to existing type strains. Additionally, the identification of genomic islands (GIs) and virulence-associated factors indicated their genetic diversity and potential adaptive capabilities, as well as potential implications for human health. Notably, CRISPR-Cas systems in strain NCTC 9931 underscore its adaptive immune mechanisms compared to strain C24.

CONCLUSIONS

Based on these findings, strain NCTC 9931T (= ATCC 17982T = DSM 43331T = CIP 104728T = CCUG 18309T = NCTC 14978T = CGMCC 1.90328T) represents a novel species, for which the name Schaalia dentiphila subsp. dentiphila sp. nov. subsp. nov. is proposed, while strain C24T (= NCTC 14980T = CGMCC 1.90329T) represents a distinct novel subspecies, for which the name Schaalia dentiphila subsp. denticola. subsp. nov. is proposed. This study enriches our understanding of the genomic diversity of Schaalia species and paves the way for further investigations into their roles in oral health.

SIGNIFICANCE

This research reveals two Schaalia strains, NCTC 9931 T and C24T, as novel entities with distinct genomic features. Expanding the taxonomic framework of the genus Schaalia, this study offers a critical resource for probing the metabolic intricacies and resistance patterns of these bacteria. This work stands as a cornerstone for microbial taxonomy, paving the way for significant advances in clinical diagnostics.

scCDC: a computational method for gene-specific contamination detection and correction in single-cell and single-nucleus RNA-seq data

In droplet-based single-cell and single-nucleus RNA-seq assays, systematic contamination of ambient RNA molecules biases the quantification of gene expression levels. Existing methods correct the contamination for all genes globally. However, there lacks specific evaluation of correction efficacy for varying contamination levels. Here, we show that DecontX and CellBender under-correct highly contaminating genes, while SoupX and scAR over-correct lowly/non-contaminating genes. Here, we develop scCDC as the first method to detect the contamination-causing genes and only correct expression levels of these genes, some of which are cell-type markers. Compared with existing decontamination methods, scCDC excels in decontaminating highly contaminating genes while avoiding over-correction of other genes.

A new chromosome-scale genome of wild Brassica oleracea provides insights into the domestication of Brassica crops

The cultivated diploid Brassica oleracea is an important vegetable crop, but the genetic basis of its domestication remains largely unclear in the absence of high-quality reference genomes of wild B. oleracea. Here, we report the first chromosome-level assembly of the wild Brassica oleracea L. W03 genome (total genome size, 630.7 Mb; scaffold N50, 64.6 Mb). Using the newly assembled W03 genome, we constructed a gene-based B. oleracea pangenome and identified 29 744 core genes, 23 306 dispensable genes, and 1896 private genes. We re-sequenced 53 accessions, representing six potential wild B. oleracea progenitor species. The results of the population genomic analysis showed that the wild B. oleracea populations had the highest level of diversity and represents the most closely related population to modern-day horticultural B. oleracea. In addition, the WUSCHEL gene was found to play a decisive role in domestication and to be involved in cauliflower and broccoli curd formation. We also illustrate the loss of disease-resistance genes during selection for domestication. Our results provide new insights into the domestication of B. oleracea and will facilitate the future genetic improvement of Brassica crops.

Combining DNA and protein alignments to improve genome annotation with LiftOn

As the number and variety of assembled genomes continues to grow, the number of annotated genomes is falling behind, particularly for eukaryotes. DNA-based mapping tools help to address this challenge, but they are only able to transfer annotation between closely-related species. Here we introduce LiftOn, a homology-based software tool that integrates DNA and protein alignments to enhance the accuracy of genome-scale annotation and to allow mapping between relatively distant species. LiftOn's protein-centric algorithm considers both types of alignments, chooses optimal open reading frames, resolves overlapping gene loci, and finds additional gene copies where they exist. LiftOn can reliably transfer annotation between genomes representing members of the same species, as we demonstrate on human, mouse, honey bee, rice, and Arabidopsis thaliana. It can further map annotation effectively across species pairs as far apart as mouse and rat or Drosophila melanogaster and D. erecta.

Pre-Bleaching Coral Microbiome Is Enriched in Beneficial Taxa and Functions

Coral reef health is tightly connected to the coral holobiont, which is the association between the coral animal and a diverse microbiome functioning as a unit. The coral holobiont depends on key services such as nitrogen and sulfur cycling mediated by the associated bacteria. However, these microbial services may be impaired in response to environmental changes, such as thermal stress. A perturbed microbiome may lead to coral bleaching and disease outbreaks, which have caused an unprecedented loss in coral cover worldwide, particularly correlated to a warming ocean. The response mechanisms of the coral holobiont under high temperatures are not completely understood, but the associated microbial community is a potential source of acquired heat-tolerance. Here we investigate the effects of increased temperature on the taxonomic and functional profiles of coral surface mucous layer (SML) microbiomes in relationship to coral-algal physiology. We used shotgun metagenomics in an experimental setting to understand the dynamics of microbial taxa and genes in the SML microbiome of the coral Pseudodiploria strigosa under heat treatment. The metagenomes of corals exposed to heat showed high similarity at the level of bacterial genera and functional genes related to nitrogen and sulfur metabolism and stress response. The coral SML microbiome responded to heat with an increase in the relative abundance of taxa with probiotic potential, and functional genes for nitrogen and sulfur acquisition. Coral-algal physiology significantly explained the variation in the microbiome at taxonomic and functional levels. These consistent and specific microbial taxa and gene functions that significantly increased in proportional abundance in corals exposed to heat are potentially beneficial to coral health and thermal resistance.

In silico genome wide identification of long non-coding RNAs differentially expressed during Candida auris host pathogenesis

Candida auris is an invasive fungal pathogen of high concern due to acquired drug tolerance against antifungals used in clinics. The prolonged persistence on biotic and abiotic surfaces can result in onset of hospital outbreaks causing serious health threat. An in depth understanding of pathology of C. auris is highly desirable for development of efficient therapeutics. Non-coding RNAs play crucial role in fungal pathology. However, the information about ncRNAs is scanty to be utilized. Herein our aim is to identify long noncoding RNAs with potent role in pathobiology of C. auris. Thereby, we analyzed the transcriptomics data of C. auris infection in blood for identification of potential lncRNAs with regulatory role in determining invasion, survival or drug tolerance under infection conditions. Interestingly, we found 275 lncRNAs, out of which 253 matched with lncRNAs reported in Candidamine, corroborating for our accurate data analysis pipeline. Nevertheless, we obtained 23 novel lncRNAs not reported earlier. Three lncRNAs were found to be under expressed throughout the course of infection, in the transcriptomics data. 16 of potent lncRNAs were found to be coexpressed with coding genes, emphasizing for their functional role. Noteworthy, these ncRNAs are expressed from intergenic regions of the genes associated with transporters, metabolism, cell wall biogenesis. This study recommends for possible association between lncRNA expression and C. auris pathogenesis.

Immunoinformatics assisted profiling of West Nile virus proteome to determine immunodominant epitopes for the development of next-generation multi-peptide vaccine

Emerging infectious diseases represent a significant threat to global health, with West Nile virus (WNV) being a prominent example due to its potential to cause severe neurological disorders alongside mild feverish conditions. Particularly prevalent in the continental United States, WNV has emerged as a global concern, with outbreaks indicating the urgent need for effective prophylactic measures. The current problem is that the absence of a commercial vaccine against WNV highlights a critical gap in preventive strategies against WNV. This study aims to address this gap by proposing a novel, multivalent vaccine designed using immunoinformatics approaches to elicit comprehensive humoral and cellular immune responses against WNV. The objective of the study is to provide a theoretical framework for experimental scientists to formulate of vaccine against WNV and tackle the current problem by generating an immune response inside the host. The research employs reverse vaccinology and subtractive proteomics methodologies to identify NP_041724.2 polyprotein and YP_009164950.1 truncated flavivirus polyprotein NS1 as the prime antigens. The selection process for epitopes focused on B and T-cell reactivity, antigenicity, water solubility, and non-allergenic properties, prioritizing candidates with the potential for broad immunogenicity and safety. The designed vaccine construct integrates these epitopes, connected via GPGPG linkers, and supplemented with an adjuvant with the help of another linker EAAAK, to enhance immunogenicity. Preliminary computational analyses suggest that the proposed vaccine could achieve near-universal coverage, effectively targeting approximately 99.74% of the global population, with perfect coverage in specific regions such as Sweden and Finland. Molecular docking and immune simulation studies further validate the potential efficacy of the vaccine, indicating strong binding affinity with toll-like receptor 3 (TLR-3) and promising immune response profiles, including significant antibody-mediated and cellular responses. These findings present the vaccine construct as a viable candidate for further development and testing. While the theoretical and computational results are promising, advancing from in-silico predictions to a tangible vaccine requires comprehensive laboratory validation. This next step is essential to confirm the vaccine's efficacy and safety in eliciting an immune response against WNV. Through this study, we propose a novel approach to vaccine development against WNV and contribute to the broader field of immunoinformatics, showcasing the potential to accelerate the design of effective vaccines against emerging viral threats. The journey from hypothesis to practical solution embodies the interdisciplinary collaboration essential for modern infectious disease management and prevention strategies.

Evolutionary mining and functional characterization of TnpB nucleases identify efficient miniature genome editors

As the evolutionary ancestor of Cas12 nuclease, the transposon (IS200/IS605)-encoded TnpB proteins act as compact RNA-guided DNA endonucleases. To explore their evolutionary diversity and potential as genome editors, we screened TnpBs from 64 annotated IS605 members and identified 25 active in Escherichia coli, of which three are active in human cells. Further characterization of these 25 TnpBs enables prediction of the transposon-associated motif (TAM) and the right-end element RNA (reRNA) directly from genomic sequences. We established a framework for annotating TnpB systems in prokaryotic genomes and applied it to identify 14 additional candidates. Among these, ISAam1 (369 amino acids (aa)) and ISYmu1 (382 aa) TnpBs demonstrated robust editing activity across dozens of genomic loci in human cells. Both RNA-guided genome editors demonstrated similar editing efficiency as SaCas9 (1,053 aa) while being substantially smaller. The enormous diversity of TnpBs holds potential for the discovery of additional valuable genome editors.

Hypothetical chloroplast reading frame 51 encodes a photosystem I assembly factor in cyanobacteria

Hypothetical chloroplast open reading frames (ycfs) are putative genes in the plastid genomes of photosynthetic eukaryotes. Many ycfs are also conserved in the genomes of cyanobacteria, the presumptive ancestors of present-day chloroplasts. The functions of many ycfs are still unknown. Here, we generated knock-out mutants for ycf51 (sll1702) in the cyanobacterium Synechocystis sp. PCC 6803. The mutants showed reduced photoautotrophic growth due to impaired electron transport between photosystem II (PSII) and PSI. This phenotype results from greatly reduced PSI content in the ycf51 mutant. The ycf51 disruption had little effect on the transcription of genes encoding photosynthetic complex components and the stabilization of the PSI complex. In vitro and in vivo analyses demonstrated that Ycf51 cooperates with PSI assembly factor Ycf3 to mediate PSI assembly. Furthermore, Ycf51 interacts with the PSI subunit PsaC. Together with its specific localization in the thylakoid membrane and the stromal exposure of its hydrophilic region, our data suggest that Ycf51 is involved in PSI complex assembly. Ycf51 is conserved in all sequenced cyanobacteria, including the earliest branching cyanobacteria of the Gloeobacter genus, and is also present in the plastid genomes of glaucophytes. However, Ycf51 has been lost from other photosynthetic eukaryotic lineages. Thus, Ycf51 is a PSI assembly factor that has been functionally replaced during the evolution of oxygenic photosynthetic eukaryotes.

Transposable elements cause the loss of self-incompatibility in citrus

Self-incompatibility (SI) is a widespread prezygotic mechanism for flowering plants to avoid inbreeding depression and promote genetic diversity. Citrus has an S-RNase-based SI system, which was frequently lost during evolution. We previously identified a single nucleotide mutation in Sm-RNase, which is responsible for the loss of SI in mandarin and its hybrids. However, little is known about other mechanisms responsible for conversion of SI to self-compatibility (SC) and we identify a completely different mechanism widely utilized by citrus. Here, we found a 786-bp miniature inverted-repeat transposable element (MITE) insertion in the promoter region of the FhiS2-RNase in Fortunella hindsii Swingle (a model plant for citrus gene function), which does not contain the Sm-RNase allele but are still SC. We demonstrate that this MITE plays a pivotal role in the loss of SI in citrus, providing evidence that this MITE insertion prevents expression of the S-RNase; moreover, transgenic experiments show that deletion of this 786-bp MITE insertion recovers the expression of FhiS2-RNase and restores SI. This study identifies the first evidence for a role for MITEs at the S-locus affecting the SI phenotype. A family-wide survey of the S-locus revealed that MITE insertions occur frequently adjacent to S-RNase alleles in different citrus genera, but only certain MITEs appear to be responsible for the loss of SI. Our study provides evidence that insertion of MITEs into a promoter region can alter a breeding strategy and suggests that this phenomenon may be broadly responsible for SC in species with the S-RNase system.

Centrifuger: lossless compression of microbial genomes for efficient and accurate metagenomic sequence classification

Centrifuger is an efficient taxonomic classification method that compares sequencing reads against a microbial genome database. In Centrifuger, the Burrows-Wheeler transformed genome sequences are losslessly compressed using a novel scheme called run-block compression. Run-block compression achieves sublinear space complexity and is effective at compressing diverse microbial databases like RefSeq while supporting fast rank queries. Combining this compression method with other strategies for compacting the Ferragina-Manzini (FM) index, Centrifuger reduces the memory footprint by half compared to other FM-index-based approaches. Furthermore, the lossless compression and the unconstrained match length help Centrifuger achieve greater accuracy than competing methods at lower taxonomic levels.

Comparative Genomics Reveals Genetic Diversity and Variation in Metabolic Traits in Fructilactobacillus sanfranciscensis Strains

Fructilactobacillus sanfranciscensis is a significant and dominant bacterial species of sourdough microbiota from ecological and functional perspectives. Despite the remarkable prevalence of different strains of this species in sourdoughs worldwide, the drivers behind the genetic diversity of this species needed to be clarified. In this research, 14 F. sanfranciscensis strains were isolated from sourdough samples to evaluate the genetic diversity and variation in metabolic traits. These 14 and 31 other strains (obtained from the NCBI database) genomes were compared. The values for genome size and GC content, on average, turned out to 1.31 Mbp and 34.25%, respectively. In 45 F. sanfranciscensis strains, there were 162 core genes and 0 to 51 unique genes present in each strain. The primary functions of core genes were related to nucleotide, lipid transport, and amino acid, as well as carbohydrate metabolism. The size of core genes accounted for 41.18% of the pan-genome size in 14 F. sanfranciscensis strains, i.e., 0.70 Mbp of 1.70 Mbp. There were genetic variations among the 14 strains involved in carbohydrate utilization and antibiotic resistance. Moreover, exopolysaccharides biosynthesis-related genes were annotated, including epsABD, wxz, wzy. The Type IIA & IE CRISPR-Cas systems, pediocin PA-1 and Lacticin_3147_A1 bacteriocins operons were also discovered in F. sanfranciscensis. These findings can help to select desirable F. sanfranciscensis strains to develop standardized starter culture for sourdough fermentation, and expect to provide traditional fermented pasta with a higher quality and nutritional value for the consumers.