Improvements in viral gene annotation using large language models and soft alignments.

BACKGROUND: The annotation of protein sequences in public databases has long posed a challenge in molecular biology. This issue is particularly acute for viral proteins, which demonstrate limited homology to known proteins when using alignment, k-mer, or profile-based homology search approaches. A novel methodology employing Large Language Models (LLMs) addresses this methodological challenge by annotating protein sequences based on embeddings. RESULTS: Central to our contribution is the soft alignment algorithm, drawing from traditional protein alignment but leveraging embedding similarity at the amino acid level to bypass the need for conventional scoring matrices. This method not only surpasses pooled embedding-based models in efficiency but also in interpretability, enabling users to easily trace homologous amino acids and delve deeper into the alignments. Far from being a black box, our approach provides transparent, BLAST-like alignment visualizations, combining traditional biological research with AI advancements to elevate protein annotation through embedding-based analysis while ensuring interpretability. Tests using the Virus Orthologous Groups and ViralZone protein databases indicated that the novel soft alignment approach recognized and annotated sequences that both blastp and pooling-based methods, which are commonly used for sequence annotation, failed to detect. CONCLUSION: The embeddings approach shows the great potential of LLMs for enhancing protein sequence annotation, especially in viral genomics. These findings present a promising avenue for more efficient and accurate protein function inference in molecular biology.

A ridge-to-reef ecosystem microbial census reveals environmental reservoirs for animal and plant microbiomes.

Microbes are found in nearly every habitat and organism on the planet, where they are critical to host health, fitness, and metabolism. In most organisms, few microbes are inherited at birth; instead, acquiring microbiomes generally involves complicated interactions between the environment, hosts, and symbionts. Despite the criticality of microbiome acquisition, we know little about where hosts' microbes reside when not in or on hosts of interest. Because microbes span a continuum ranging from generalists associating with multiple hosts and habitats to specialists with narrower host ranges, identifying potential sources of microbial diversity that can contribute to the microbiomes of unrelated hosts is a gap in our understanding of microbiome assembly. Microbial dispersal attenuates with distance, so identifying sources and sinks requires data from microbiomes that are contemporary and near enough for potential microbial transmission. Here, we characterize microbiomes across adjacent terrestrial and aquatic hosts and habitats throughout an entire watershed, showing that the most species-poor microbiomes are partial subsets of the most species-rich and that microbiomes of plants and animals are nested within those of their environments. Furthermore, we show that the host and habitat range of a microbe within a single ecosystem predicts its global distribution, a relationship with implications for global microbial assembly processes. Thus, the tendency for microbes to occupy multiple habitats and unrelated hosts enables persistent microbiomes, even when host populations are disjunct. Our whole-watershed census demonstrates how a nested distribution of microbes, following the trophic hierarchies of hosts, can shape microbial acquisition.

Leveraging deep contrastive learning for semantic interaction.

The semantic interaction process seeks to elicit a user's mental model as they interact with and query visualizations during a sense-making activity. Semantic interaction enables the development of computational models that capture user intent and anticipate user actions. Deep learning is proving to be highly effective for learning complex functions and is, therefore, a compelling tool for encoding a user's mental model. In this paper, we show that deep contrastive learning significantly enhances semantic interaction in visual analytics systems. Our approach does so by allowing users to explore alternative arrangements of their data while simultaneously training a parametric algorithm to learn their evolving mental model. As an example of the efficacy of our approach, we deployed our model in Z-Explorer, a visual analytics extension to the widely used Zotero document management system. The user study demonstrates that this flexible approach effectively captures users' mental data models without explicit hyperparameter tuning or even requiring prior machine learning expertise.

CoCoNet: an efficient deep learning tool for viral metagenome binning.

MOTIVATION: Metagenomic approaches hold the potential to characterize microbial communities and unravel the intricate link between the microbiome and biological processes. Assembly is one of the most critical steps in metagenomics experiments. It consists of transforming overlapping DNA sequencing reads into sufficiently accurate representations of the community's genomes. This process is computationally difficult and commonly results in genomes fragmented across many contigs. Computational binning methods are used to mitigate fragmentation by partitioning contigs based on their sequence composition, abundance or chromosome organization into bins representing the community's genomes. Existing binning methods have been principally tuned for bacterial genomes and do not perform favorably on viral metagenomes. RESULTS: We propose Composition and Coverage Network (CoCoNet), a new binning method for viral metagenomes that leverages the flexibility and the effectiveness of deep learning to model the co-occurrence of contigs belonging to the same viral genome and provide a rigorous framework for binning viral contigs. Our results show that CoCoNet substantially outperforms existing binning methods on viral datasets. AVAILABILITY AND IMPLEMENTATION: CoCoNet was implemented in Python and is available for download on PyPi (https://pypi.org/). The source code is hosted on GitHub at https://github.com/Puumanamana/CoCoNet and the documentation is available at https://coconet.readthedocs.io/en/latest/index.html. CoCoNet does not require extensive resources to run. For example, binning 100k contigs took about 4 h on 10 Intel CPU Cores (2.4 GHz), with a memory peak at 27 GB (see Supplementary Fig. S9). To process a large dataset, CoCoNet may need to be run on a high RAM capacity server. Such servers are typically available in high-performance or cloud computing settings. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The amphibian microbiome exhibits poor resilience following pathogen-induced disturbance.

Infectious pathogens can disrupt the microbiome in addition to directly affecting the host. Impacts of disease may be dependent on the ability of the microbiome to recover from such disturbance, yet remarkably little is known about microbiome recovery after disease, particularly in nonhuman animals. We assessed the resilience of the amphibian skin microbial community after disturbance by the pathogen, Batrachochytrium dendrobatidis (Bd). Skin microbial communities of laboratory-reared mountain yellow-legged frogs were tracked through three experimental phases: prior to Bd infection, after Bd infection (disturbance), and after clearing Bd infection (recovery period). Bd infection disturbed microbiome composition and altered the relative abundances of several dominant bacterial taxa. After Bd infection, frogs were treated with an antifungal drug that cleared Bd infection, but this did not lead to recovery of microbiome composition (measured as Unifrac distance) or relative abundances of dominant bacterial groups. These results indicate that Bd infection can lead to an alternate stable state in the microbiome of sensitive amphibians, or that microbiome recovery is extremely slow-in either case resilience is low. Furthermore, antifungal treatment and clearance of Bd infection had the additional effect of reducing microbial community variability, which we hypothesize results from similarity across frogs in the taxa that colonize community vacancies resulting from the removal of Bd. Our results indicate that the skin microbiota of mountain yellow-legged frogs has low resilience following Bd-induced disturbance and is further altered by the process of clearing Bd infection, which may have implications for the conservation of this endangered amphibian.

Transcriptional regulation of prolactin in a euryhaline teleost: Characterisation of gene promoters through in silico and transcriptome analyses.

The sensitivity of prolactin (Prl) cells of the Mozambique tilapia (Oreochromis mossambicus) pituitary to variations in extracellular osmolality enables investigations into how osmoreception underlies patterns of hormone secretion. Through the actions of their main secretory products, Prl cells play a key role in supporting hydromineral balance of fishes by controlling the major osmoregulatory organs (ie, gill, intestine and kidney). The release of Prl from isolated cells of the rostral pars distalis (RPD) occurs in direct response to physiologically relevant reductions in extracellular osmolality. Although the particular signal transduction pathways that link osmotic conditions to Prl secretion have been identified, the processes that underlie hyposmotic induction of prl gene expression remain unknown. In this short review, we describe two distinct tilapia gene loci that encode Prl177 and Prl188 . From our in silico analyses of prl177 and prl188 promoter regions (approximately 1000 bp) and a transcriptome analysis of RPDs from fresh water (FW)- and seawater (SW)-acclimated tilapia, we propose a working model for how multiple transcription factors link osmoreceptive processes with adaptive patterns of prl177 and prl188 gene expression. We confirmed via RNA-sequencing and a quantitative polymerase chain reaction that multiple transcription factors emerging as predicted regulators of prl gene expression are expressed in the RPD of tilapia. In particular, gene transcripts encoding pou1f1, stat3, creb3l1, pbxip1a and stat1a were highly expressed; creb3l1, pbxip1a and stat1a were elevated in fish acclimated to SW vs FW. Combined, our in silico and transcriptome analyses set a path for resolving how adaptive patterns of Prl secretion are achieved via the integration of osmoreceptive processes with the control of prl gene transcription.

Transcriptome Profile of Nicotinic Receptor-Linked Sensitization of Beta Amyloid Neurotoxicity.

Understanding the specific gene changes underlying the prodromic stages of Alzheimer's disease pathogenesis will aid the development of new, targeted therapeutic strategies for this neurodegenerative disorder. Here, we employed RNA-sequencing to analyze global differential gene expression in a defined model nerve cell line expressing α4ß2 nicotinic receptors (nAChRs), high-affinity targets for beta amyloid (Aß). The nAChR-expressing neuronal cells were treated with nanomolar Aß1-42 to gain insights into the molecular mechanisms underlying Aß-induced neurotoxicity in the presence of this sensitizing target receptor. We identified 15 genes (out of 15,336) that were differentially expressed upon receptor-linked Aß treatment. Genes up-regulated with Aß treatment were associated with calcium signaling and axonal vesicle transport (including the α4 nAChR subunit, the calcineurin regulator RCAN3, and KIF1C of the kinesin family). Downregulated genes were associated with metabolic, apoptotic or DNA repair pathways (including APBA3, PARP1 and RAB11). Validation of the differential expression was performed via qRT-PCR and immunoblot analysis in the defined model nerve cell line and primary mouse neurons. Further verification was performed using immunocytochemistry. In conclusion, we identified apparent changes in gene expression on Aß treatment in the presence of the sensitizing nAChRs, linked to early-stage Aß-induced neurotoxicity, which may represent novel therapeutic targets.

Critical symbiont signals drive both local and systemic changes in diel and developmental host gene expression.

The colonization of an animal's tissues by its microbial partners creates networks of communication across the host's body. We used the natural binary light-organ symbiosis between the squid Euprymna scolopes and its luminous bacterial partner, Vibrio fischeri, to define the impact of colonization on transcriptomic networks in the host. A night-active predator, E. scolopes coordinates the bioluminescence of its symbiont with visual cues from the environment to camouflage against moon and starlight. Like mammals, this symbiosis has a complex developmental program and a strong day/night rhythm. We determined how symbiont colonization impacted gene expression in the light organ itself, as well as in two anatomically remote organs: the eye and gill. While the overall transcriptional signature of light organ and gill were more alike, the impact of symbiosis was most pronounced and similar in light organ and eye, both in juvenile and adult animals. Furthermore, the presence of a symbiosis drove daily rhythms of transcription within all three organs. Finally, a single mutation in V. fischeri-specifically, deletion of the lux operon, which abrogates symbiont luminescence-reduced the symbiosis-dependent transcriptome of the light organ by two-thirds. In addition, while the gills responded similarly to light-organ colonization by either the wild-type or mutant, luminescence was required for all of the colonization-associated transcriptional responses in the juvenile eye. This study defines not only the impact of symbiont colonization on the coordination of animal transcriptomes, but also provides insight into how such changes might impact the behavior and ecology of the host.

Symbiotic organs shaped by distinct modes of genome evolution in cephalopods.

Microbes have been critical drivers of evolutionary innovation in animals. To understand the processes that influence the origin of specialized symbiotic organs, we report the sequencing and analysis of the genome of Euprymna scolopes, a model cephalopod with richly characterized host-microbe interactions. We identified large-scale genomic reorganization shared between E. scolopes and Octopus bimaculoides and posit that this reorganization has contributed to the evolution of cephalopod complexity. To reveal genomic signatures of host-symbiont interactions, we focused on two specialized organs of E. scolopes: the light organ, which harbors a monoculture of Vibrio fischeri, and the accessory nidamental gland (ANG), a reproductive organ containing a bacterial consortium. Our findings suggest that the two symbiotic organs within E. scolopes originated by different evolutionary mechanisms. Transcripts expressed in these microbe-associated tissues displayed their own unique signatures in both coding sequences and the surrounding regulatory regions. Compared with other tissues, the light organ showed an abundance of genes associated with immunity and mediating light, whereas the ANG was enriched in orphan genes known only from E. scolopes Together, these analyses provide evidence for different patterns of genomic evolution of symbiotic organs within a single host.

Ultra-Deep Sequencing Analysis on HIV Drug-Resistance-Associated Mutations Among HIV-Infected Individuals: First Report from the Philippines.

A sharp increase in the number of people living with HIV has been documented in the Philippines. In response, the government has instituted antiretroviral therapy (ART) nationwide through HIV treatment hubs. However, no data presently exist on the status of ART drug-resistance-associated mutations (DRMs). In this study, we aim at analyzing DRM profiles in the Philippines and at providing comprehensive data on DRMs to guide treatment decisions and prevent viral failures. We conducted a cross-sectional study in 119 volunteers who tested positive for HIV from more than 8,000 participants screened for HIV across the nation through the 2013 Integrated HIV Behavioral and Serologic Surveillance (IHBSS) program. Amplicons were generated from plasma RNA by using primers designed to analyze diverse HIV-1 isolates targeting the reverse transcriptase region and sequenced on a 454 ultra-deep sequencing (UDS) platform to assess DRMs. DRMs were defined by using the Stanford HIV drug resistance database, and we found only 2 from 110 evaluable individuals with major HIV variants (>20% prevalence) that were highly resistant to the non-nucleoside reverse transcriptase inhibitor (NNRTI: efavirenz and nevirapine). However, a larger fraction of individuals harbored minority drug-resistant HIV variants (0.5%-20% prevalence) and they were highly resistant to NNRTI nevirapine (89/110), rilpivirine (5/110), and efavirenz (49/110). This study is the first report on the presence of HIV drug resistance in the Philippines and demonstrates the utility of UDS in assisting the detection of HIV minor variants. Monitoring for ART-DRMs will assist in improving HIV management strategies in curtailing the evolving epidemic in the Philippines.

Comparative genomics explains the evolutionary success of reef-forming corals.

Transcriptome and genome data from twenty stony coral species and a selection of reference bilaterians were studied to elucidate coral evolutionary history. We identified genes that encode the proteins responsible for the precipitation and aggregation of the aragonite skeleton on which the organisms live, and revealed a network of environmental sensors that coordinate responses of the host animals to temperature, light, and pH. Furthermore, we describe a variety of stress-related pathways, including apoptotic pathways that allow the host animals to detoxify reactive oxygen and nitrogen species that are generated by their intracellular photosynthetic symbionts, and determine the fate of corals under environmental stress. Some of these genes arose through horizontal gene transfer and comprise at least 0.2% of the animal gene inventory. Our analysis elucidates the evolutionary strategies that have allowed symbiotic corals to adapt and thrive for hundreds of millions of years.

Identification of host genes leading to West Nile virus encephalitis in mice brain using RNA-seq analysis.

Differential host responses may be critical determinants of distinct pathologies of West Nile virus (WNV) NY99 (pathogenic) and WNV Eg101 (non-pathogenic) strains. We employed RNA-seq technology to analyze global differential gene expression in WNV-infected mice brain and to identify the host cellular factors leading to lethal encephalitis. We identified 1,400 and 278 transcripts, which were differentially expressed after WNV NY99 and WNV Eg101 infections, respectively, and 147 genes were common to infection with both the viruses. Genes that were up-regulated in infection with both the viruses were mainly associated with interferon signaling. Genes associated with inflammation and cell death/apoptosis were only expressed after WNV NY99 infection. We demonstrate that differences in the activation of key pattern recognition receptors resulted in the induction of unique innate immune profiles, which corresponded with the induction of interferon and inflammatory responses. Pathway analysis of differentially expressed genes indicated that after WNV NY99 infection, TREM-1 mediated activation of toll-like receptors leads to the high inflammatory response. In conclusion, we have identified both common and specific responses to WNV NY99 and WNV Eg101 infections as well as genes linked to potential resistance to infection that may be targets for therapeutics.

Characterization of class II ß chain major histocompatibility complex genes in a family of Hawaiian honeycreepers: 'amakihi (Hemignathus virens).

Hawaiian honeycreepers (Drepanidinae) have evolved in the absence of mosquitoes for over five million years. Through human activity, mosquitoes were introduced to the Hawaiian archipelago less than 200 years ago. Mosquito-vectored diseases such as avian malaria caused by Plasmodium relictum and Avipoxviruses have greatly impacted these vulnerable species. Susceptibility to these diseases is variable among and within species. Due to their function in adaptive immunity, the role of major histocompatibility complex genes (Mhc) in disease susceptibility is under investigation. In this study, we evaluate gene organization and levels of diversity of Mhc class II ß chain genes (exon 2) in a captive-reared family of Hawaii 'amakihi (Hemignathus virens). A total of 233 sequences (173 bp) were obtained by PCR+1 amplification and cloning, and 5720 sequences were generated by Roche 454 pyrosequencing. We report a total of 17 alleles originating from a minimum of 14 distinct loci. We detected three linkage groups that appear to represent three distinct haplotypes. Phylogenetic analysis revealed one variable cluster resembling classical Mhc sequences (DAB) and one highly conserved, low variability cluster resembling non-classical Mhc sequences (DBB). High net evolutionary divergence values between DAB and DBB resemble that seen between chicken BLB system and YLB system genes. High amino acid identity among non-classical alleles from 12 species of passerines (DBB) and four species of Galliformes (YLB) was found, suggesting that these non-classical passerine sequences may be related to the Galliforme YLB sequences.

Identification and preliminary characterization of chemosensory perception-associated proteins in the melon fly Bactrocera cucurbitae using RNA-seq.

An investigation into proteins involved in chemosensory perception in the melon fly, Bactrocera cucurbitae (Diptera: Tephritidae) is described here using a newly generated transcriptome dataset. The melon fly is a major agricultural pest, widely distributed in the Asia-Pacific region and some parts of Africa. For this study, a transcriptome dataset was generated using RNA extracted from 4-day-old adult specimens of the melon fly. The dataset was assembled and annotated via Gene Ontology (GO) analysis. Based on this and similarity searches to data from other species, a number of protein sequences putatively involved in chemosensory reception were identified and characterized in the melon fly. This included the highly conserved "Orco" along with a number of other less conserved odorant binding protein sequences. In addition, several sequences representing putative ionotropic and gustatory receptors were also identified. This study provides a foundation for future functional studies of chemosensory proteins in the melon fly and for making more detailed comparisons to other species. In the long term, this will ultimately help in the development of improved tools for pest management.

Preclinical development of HIvax: Human survivin highly immunogenic vaccines.

Our previous work involved the development of a recombinant fowlpox virus encoding survivin (FP-surv) vaccine that was evaluated for efficacy in mesothelioma mouse models. Results showed that FP-surv vaccination generated significant immune responses, which led to delayed tumor growth and improved animal survival. We have extended those previous findings in the current study, which involves the pre-clinical development of an optimized version of FP-surv designed for human immunization (HIvax). Survivin-derived peptides for the most common haplotypes in the human population were identified and their immunogenicity confirmed in co-culture experiments using dendritic cells and T cells isolated from healthy donors. Peptides confirmed to induce CD8(+) and CD4(+) T cells activation in humans were then included in 2 transgenes optimized for presentation of processed peptides on MHC-I (HIvax1) and MHC-II (HIvax2). Fowlpox vectors expressing the HIvax transgenes were then generated and their efficacy was evaluated with subsequent co-culture experiments to measure interferon-γ and granzyme B secretion. In these experiments, both antigen specific CD4(+) and CD8(+) T cells were activated by HIvax vaccines with resultant cytotoxic activity against survivin-overexpressing mesothelioma cancer cells. These results provide a rationale for clinical testing of HIvax1 and HIvax2 vaccines in patients with survivin-expressing cancers.

Complete Genome Sequence of Burkholderia cepacia Strain LO6.

Burkholderia cepacia strain LO6 is a betaproteobacterium that was isolated from a cystic fibrosis patient. Here we report the 6.4 Mb draft genome sequence assembled into 2 contigs. This genome sequence will aid the transcriptomic profiling of this bacterium and help us to better understand the mechanisms specific to pulmonary infections.

Demystifying computer science for molecular ecologists.

In this age of data-driven science and high-throughput biology, computational thinking is becoming an increasingly important skill for tackling both new and long-standing biological questions. However, despite its obvious importance and conspicuous integration into many areas of biology, computer science is still viewed as an obscure field that has, thus far, permeated into only a few of the biology curricula across the nation. A national survey has shown that lack of computational literacy in environmental sciences is the norm rather than the exception [Valle & Berdanier (2012) Bulletin of the Ecological Society of America, 93, 373-389]. In this article, we seek to introduce a few important concepts in computer science with the aim of providing a context-specific introduction aimed at research biologists. Our goal was to help biologists understand some of the most important mainstream computational concepts to better appreciate bioinformatics methods and trade-offs that are not obvious to the uninitiated.

Complete Genome Sequence of Pseudoalteromonas sp. Strain OCN003, Isolated from Kane'ohe Bay, O'ahu, Hawaii.

Pseudoalteromonas sp. strain OCN003 is a marine gammaproteobacterium that was isolated from a diseased colony of the common Hawaiian reef coral, Montipora capitata, found on a reef surrounding Moku o Lo'e in Kane'ohe Bay, Hawaii. Here, we report the complete genome of Pseudoalteromonas sp. strain OCN003.

Are we missing half of the viruses in the ocean?

Viruses are abundant in the ocean and a major driving force in plankton ecology and evolution. It has been assumed that most of the viruses in seawater contain DNA and infect bacteria, but RNA-containing viruses in the ocean, which almost exclusively infect eukaryotes, have never been quantified. We compared the total mass of RNA and DNA in the viral fraction harvested from seawater and using data on the mass of nucleic acid per RNA- or DNA-containing virion, estimated the abundances of each. Our data suggest that the abundance of RNA viruses rivaled or exceeded that of DNA viruses in samples of coastal seawater. The dominant RNA viruses in the samples were marine picorna-like viruses, which have small genomes and are at or below the detection limit of common fluorescence-based counting methods. If our results are typical, this means that counts of viruses and the rate measurements that depend on them, such as viral production, are significantly underestimated by current practices. As these RNA viruses infect eukaryotes, our data imply that protists contribute more to marine viral dynamics than one might expect based on their relatively low abundance. This conclusion is a departure from the prevailing view of viruses in the ocean, but is consistent with earlier theoretical predictions.

Chimeric piggyBac transposases for genomic targeting in human cells.

Integrating vectors such as viruses and transposons insert transgenes semi-randomly and can potentially disrupt or deregulate genes. For these techniques to be of therapeutic value, a method for controlling the precise location of insertion is required. The piggyBac (PB) transposase is an efficient gene transfer vector active in a variety of cell types and proven to be amenable to modification. Here we present the design and validation of chimeric PB proteins fused to the Gal4 DNA binding domain with the ability to target transgenes to pre-determined sites. Upstream activating sequence (UAS) Gal4 recognition sites harbored on recipient plasmids were preferentially targeted by the chimeric Gal4-PB transposase in human cells. To analyze the ability of these PB fusion proteins to target chromosomal locations, UAS sites were randomly integrated throughout the genome using the Sleeping Beauty transposon. Both N- and C-terminal Gal4-PB fusion proteins but not native PB were capable of targeting transposition nearby these introduced sites. A genome-wide integration analysis revealed the ability of our fusion constructs to bias 24% of integrations near endogenous Gal4 recognition sequences. This work provides a powerful approach to enhance the properties of the PB system for applications such as genetic engineering and gene therapy.