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.

The Tudor-domain protein TDRD7, mutated in congenital cataract, controls the heat shock protein HSPB1 (HSP27) and lens fiber cell morphology.

Mutations of the RNA granule component TDRD7 (OMIM: 611258) cause pediatric cataract. We applied an integrated approach to uncover the molecular pathology of cataract in Tdrd7-/- mice. Early postnatal Tdrd7-/- animals precipitously develop cataract suggesting a global-level breakdown/misregulation of key cellular processes. High-throughput RNA sequencing integrated with iSyTE-bioinformatics analysis identified the molecular chaperone and cytoskeletal modulator, HSPB1, among high-priority downregulated candidates in Tdrd7-/- lens. A protein fluorescence two-dimensional difference in-gel electrophoresis (2D-DIGE)-coupled mass spectrometry screen also identified HSPB1 downregulation, offering independent support for its importance to Tdrd7-/- cataractogenesis. Lens fiber cells normally undergo nuclear degradation for transparency, posing a challenge: how is their cell morphology, also critical for transparency, controlled post-nuclear degradation? HSPB1 functions in cytoskeletal maintenance, and its reduction in Tdrd7-/- lens precedes cataract, suggesting cytoskeletal defects may contribute to Tdrd7-/- cataract. In agreement, scanning electron microscopy (SEM) revealed abnormal fiber cell morphology in Tdrd7-/- lenses. Further, abnormal phalloidin and wheat germ agglutinin (WGA) staining of Tdrd7-/- fiber cells, particularly those exhibiting nuclear degradation, reveals distinct regulatory mechanisms control F-actin cytoskeletal and/or membrane maintenance in post-organelle degradation maturation stage fiber cells. Indeed, RNA immunoprecipitation identified Hspb1 mRNA in wild-type lens lysate TDRD7-pulldowns, and single-molecule RNA imaging showed co-localization of TDRD7 protein with cytoplasmic Hspb1 mRNA in differentiating fiber cells, suggesting that TDRD7-ribonucleoprotein complexes may be involved in optimal buildup of key factors. Finally, Hspb1 knockdown in Xenopus causes eye/lens defects. Together, these data uncover TDRD7's novel upstream role in elevation of stress-responsive chaperones for cytoskeletal maintenance in post-nuclear degradation lens fiber cells, perturbation of which causes early-onset cataracts.

Unraveling Fe(II)-Oxidizing Mechanisms in a Facultative Fe(II) Oxidizer, Sideroxydans lithotrophicus Strain ES-1, via Culturing, Transcriptomics, and Reverse Transcription-Quantitative PCR.

Sideroxydans lithotrophicus ES-1 grows autotrophically either by Fe(II) oxidation or by thiosulfate oxidation, in contrast to most other isolates of neutrophilic Fe(II)-oxidizing bacteria (FeOB). This provides a unique opportunity to explore the physiology of a facultative FeOB and constrain the genes specific to Fe(II) oxidation. We compared the growth of S. lithotrophicus ES-1 on Fe(II), thiosulfate, and both substrates together. While initial growth rates were similar, thiosulfate-grown cultures had higher yield with or without Fe(II) present, which may give ES-1 an advantage over obligate FeOB. To investigate the Fe(II) and S oxidation pathways, we conducted transcriptomics experiments, validated with reverse transcription-quantitative PCR (RT-qPCR). We explored the long-term gene expression response at different growth phases (over days to a week) and expression changes during a short-term switch from thiosulfate to Fe(II) (90 min). The dsr and sox sulfur oxidation genes were upregulated in thiosulfate cultures. The Fe(II) oxidase gene cyc2 was among the top expressed genes during both Fe(II) and thiosulfate oxidation, and addition of Fe(II) to thiosulfate-grown cells caused an increase in cyc2 expression. These results support the role of Cyc2 as the Fe(II) oxidase and suggest that ES-1 maintains readiness to oxidize Fe(II), even in the absence of Fe(II). We used gene expression profiles to further constrain the ES-1 Fe(II) oxidation pathway. Notably, among the most highly upregulated genes during Fe(II) oxidation were genes for alternative complex III, reverse electron transport, and carbon fixation. This implies a direct connection between Fe(II) oxidation and carbon fixation, suggesting that CO2 is an important electron sink for Fe(II) oxidation. IMPORTANCE Neutrophilic FeOB are increasingly observed in various environments, but knowledge of their ecophysiology and Fe(II) oxidation mechanisms is still relatively limited. Sideroxydans isolates are widely observed in aquifers, wetlands, and sediments, and genome analysis suggests metabolic flexibility contributes to their success. The type strain ES-1 is unusual among neutrophilic FeOB isolates, as it can grow on either Fe(II) or a non-Fe(II) substrate, thiosulfate. Almost all our knowledge of neutrophilic Fe(II) oxidation pathways comes from genome analyses, with some work on metatranscriptomes. This study used culture-based experiments to test the genes specific to Fe(II) oxidation in a facultative FeOB and refine our model of the Fe(II) oxidation pathway. We gained insight into how facultative FeOB like ES-1 connect Fe, S, and C biogeochemical cycling in the environment and suggest a multigene indicator would improve understanding of Fe(II) oxidation activity in environments with facultative FeOB.

Deletion of Fibroblast growth factor 9 globally and in skeletal muscle results in enlarged tuberosities at sites of deltoid tendon attachments.

BACKGROUND: The growth of most bony tuberosities, like the deltoid tuberosity (DT), rely on the transmission of muscle forces at the tendon-bone attachment during skeletal growth. Tuberosities distribute muscle forces and provide mechanical leverage at attachment sites for joint stability and mobility. The genetic factors that regulate tuberosity growth remain largely unknown. In mouse embryos with global deletion of fibroblast growth factor 9 (Fgf9), the DT size is notably enlarged. In this study, we explored the tissue-specific regulation of DT size using both global and targeted deletion of Fgf9. RESULTS: We showed that cell hypertrophy and mineralization dynamics of the DT, as well as transcriptional signatures from skeletal muscle but not bone, were influenced by the global loss of Fgf9. Loss of Fgf9 during embryonic growth led to increased chondrocyte hypertrophy and reduced cell proliferation at the DT attachment site. This endured hypertrophy and limited proliferation may explain the abnormal mineralization patterns and locally dysregulated expression of markers of endochondral development in Fgf9null attachments. We then showed that targeted deletion of Fgf9 in skeletal muscle leads to postnatal enlargement of the DT. CONCLUSION: Taken together, we discovered that Fgf9 may play an influential role in muscle-bone cross-talk during embryonic and postnatal development.

Bacteriophages against Vibrio coralliilyticus and Vibrio tubiashii: Isolation, Characterization, and Remediation of Larval Oyster Mortalities.

Vibrio coralliilyticus and Vibrio tubiashii are pathogens responsible for high larval oyster mortality rates in shellfish hatcheries. Bacteriophage therapy was evaluated to determine its potential to remediate these mortalities. Sixteen phages against V. coralliilyticus and V. tubiashii were isolated and characterized from Hawaiian seawater. Fourteen isolates were members of the Myoviridae family, and two were members of the Siphoviridae In proof-of-principle trials, a cocktail of five phages reduced mortalities of larval Eastern oysters (Crassostrea virginica) and Pacific oysters (Crassostrea gigas) by up to 91% 6 days after challenge with lethal doses of V. coralliilyticus Larval survival depended on the oyster species, the quantities of phages and vibrios applied, and the species and strain of Vibrio A later-generation cocktail, designated VCP300, was formulated with three lytic phages subsequently named Vibrio phages vB_VcorM-GR7B, vB_VcorM-GR11A, and vB_VcorM-GR28A (abbreviated 7B, 11A, and 28A, respectively). Together, these three phages displayed host specificity toward eight V. coralliilyticus strains and a V. tubiashii strain. Larval C. gigas mortalities from V. coralliilyticus strains RE98 and OCN008 were significantly reduced by >90% (P < 0.0001) over 6 days with phage treatment compared to those of untreated controls. Genomic sequencing of phages 7B, 11A, and 28A revealed 207,758-, 194,800-, and 154,046-bp linear DNA genomes, respectively, with the latter showing 92% similarity to V. coralliilyticus phage YC, a strain from the Great Barrier Reef, Australia. Phage 7B and 11A genomes showed little similarity to phages in the NCBI database. This study demonstrates the promising potential for phage therapy to reduce larval oyster mortalities in oyster hatcheries.IMPORTANCE Shellfish hatcheries encounter episodic outbreaks of larval oyster mortalities, jeopardizing the economic stability of hatcheries and the commercial shellfish industry. Shellfish pathogens like Vibrio coralliilyticus and Vibrio tubiashii have been recognized as major contributors of larval oyster mortalities in U.S. East and West Coast hatcheries for many years. This study isolated, identified, and characterized bacteriophages against these Vibrio species and demonstrated their ability to reduce mortalities from V. coralliilyticus in larval Pacific oysters and from both V. coralliilyticus and V. tubiashii in larval Eastern oysters. Phage therapy offers a promising approach for stimulating hatchery production to ensure the well-being of hatcheries and the commercial oyster trade.

Multiple Megaplasmids Confer Extremely High Levels of Metal Tolerance in Alteromonas Strains.

Alteromonas is a widely distributed genus of marine Gammaproteobacteria, with representatives shown to be key players in diverse processes, including biogeochemical cycling and biofouling of marine substrata. While Alteromonas spp. are early colonizers of copper-based antifouling paints on marine vessels, their mechanism of tolerance is poorly understood. PacBio whole-genome sequencing of Alteromonas macleodii strains CUKW and KCC02, isolated from Cu/Ni alloy test coupons submerged in oligotrophic coastal waters, indicated the presence of multiple megaplasmids (ca. 200 kb) in both. A pulsed-field gel electrophoresis method was developed and used to confirm the presence of multiple megaplasmids in these two strains; it was then used to screen additional Alteromonas strains for which little to no sequencing data exist. Plasmids were not detected in any of the other strains. Bioinformatic analysis of the CUKW and KCC02 plasmids identified numerous genes associated with metal resistance. Copper resistance orthologs from both the Escherichia coli Cue and Cus and Pseudomonas syringae Cop systems were present, at times as multiple copies. Metal growth assays in the presence of copper, cobalt, manganese, and zinc performed with 10 Alteromonas strains demonstrated the ability of CUKW and KCC02 to grow at metal concentrations inhibitory to all the other strains tested. This study reports multiple megaplasmids in Alteromonas strains. Bioinformatic analysis of the CUKW and KCC02 plasmids indicate that they harbor elements of the Tra system conjugation apparatus, although their type of mobility remains to be experimentally verified.IMPORTANCE Copper is commonly used as an antifouling agent on ship hulls. Alteromonas spp. are early colonizers of copper-based antifouling paint, but their mechanism of tolerance is poorly understood. Sequencing of A. macleodii strains isolated from copper test materials for marine ships indicated the presence of multiple megaplasmids. Plasmids serve as key vectors in horizontal gene transfer and confer traits such as metal resistance, detoxification, ecological interaction, and antibiotic resistance. Bioinformatic analysis identified many metal resistance genes and genes associated with mobility. Understanding the molecular mechanisms and capacity for gene transfer within marine biofilms provides a platform for the development of novel antifouling solutions targeting genes involved in copper tolerance and biofilm formation.

High-fat-diet-mediated dysbiosis promotes intestinal carcinogenesis independently of obesity.

Several features common to a Western lifestyle, including obesity and low levels of physical activity, are known risk factors for gastrointestinal cancers. There is substantial evidence suggesting that diet markedly affects the composition of the intestinal microbiota. Moreover, there is now unequivocal evidence linking dysbiosis to cancer development. However, the mechanisms by which high-fat diet (HFD)-mediated changes in the microbial community affect the severity of tumorigenesis in the gut remain to be determined. Here we demonstrate that an HFD promotes tumour progression in the small intestine of genetically susceptible, K-ras(G12Dint), mice independently of obesity. HFD consumption, in conjunction with K-ras mutation, mediated a shift in the composition of the gut microbiota, and this shift was associated with a decrease in Paneth-cell-mediated antimicrobial host defence that compromised dendritic cell recruitment and MHC class II molecule presentation in the gut-associated lymphoid tissues. When butyrate was administered to HFD-fed K-ras(G12Dint) mice, dendritic cell recruitment in the gut-associated lymphoid tissues was normalized, and tumour progression was attenuated. Importantly, deficiency in MYD88, a signalling adaptor for pattern recognition receptors and Toll-like receptors, blocked tumour progression. The transfer of faecal samples from HFD-fed mice with intestinal tumours to healthy adult K-ras(G12Dint) mice was sufficient to transmit disease in the absence of an HFD. Furthermore, treatment with antibiotics completely blocked HFD-induced tumour progression, suggesting that distinct shifts in the microbiota have a pivotal role in aggravating disease. Collectively, these data underscore the importance of the reciprocal interaction between host and environmental factors in selecting a microbiota that favours carcinogenesis, and they suggest that tumorigenesis is transmissible among genetically predisposed individuals.

Light Modulates the Physiology of Nonphototrophic Actinobacteria.

Light is a source of energy and an environmental cue that is available in excess in most surface environments. In prokaryotic systems, conversion of light to energy by photoautotrophs and photoheterotrophs is well understood, but the conversion of light to information and the cellular response to that information have been characterized in only a few species. Our goal was to explore the response of freshwater Actinobacteria, which are ubiquitous in illuminated aquatic environments, to light. We found that Actinobacteria without functional photosystems grow faster in the light, likely because sugar transport and metabolism are upregulated in the light. Based on the action spectrum of the growth effect and comparisons of the genomes of three Actinobacteria with this growth rate phenotype, we propose that the photosensor in these strains is a putative CryB-type cryptochrome. The ability to sense light and upregulate carbohydrate transport during the day could allow these cells to coordinate their time of maximum organic carbon uptake with the time of maximum organic carbon release by primary producers.IMPORTANCE Sunlight provides information about both place and time. In sunlit aquatic environments, primary producers release organic carbon and nitrogen along with other growth factors during the day. The ability of Actinobacteria to coordinate organic carbon uptake and utilization with production of photosynthate enables them to grow more efficiently in the daytime, and it potentially gives them a competitive advantage over heterotrophs that constitutively produce carbohydrate transporters, which is energetically costly, or produce transporters only after detection of the substrate(s), which delays their response. Understanding how light cues the transport of organic carbon and its conversion to biomass is key to understanding biochemical mechanisms within the carbon cycle, the fluxes through it, and the variety of mechanisms by which light enhances growth.

A reference genome of the Chinese hamster based on a hybrid assembly strategy.

Accurate and complete genome sequences are essential in biotechnology to facilitate genome-based cell engineering efforts. The current genome assemblies for Cricetulus griseus, the Chinese hamster, are fragmented and replete with gap sequences and misassemblies, consistent with most short-read-based assemblies. Here, we completely resequenced C. griseus using single molecule real time sequencing and merged this with Illumina-based assemblies. This generated a more contiguous and complete genome assembly than either technology alone, reducing the number of scaffolds by >28-fold, with 90% of the sequence in the 122 longest scaffolds. Most genes are now found in single scaffolds, including up- and downstream regulatory elements, enabling improved study of noncoding regions. With >95% of the gap sequence filled, important Chinese hamster ovary cell mutations have been detected in draft assembly gaps. This new assembly will be an invaluable resource for continued basic and pharmaceutical research.

Computational clustering for viral reference proteomes.

MOTIVATION: The enormous number of redundant sequenced genomes has hindered efforts to analyze and functionally annotate proteins. As the taxonomy of viruses is not uniformly defined, viral proteomes pose special challenges in this regard. Grouping viruses based on the similarity of their proteins at proteome scale can normalize against potential taxonomic nomenclature anomalies. RESULTS: We present Viral Reference Proteomes (Viral RPs), which are computed from complete virus proteomes within UniProtKB. Viral RPs based on 95, 75, 55, 35 and 15% co-membership in proteome similarity based clusters are provided. Comparison of our computational Viral RPs with UniProt's curator-selected Reference Proteomes indicates that the two sets are consistent and complementary. Furthermore, each Viral RP represents a cluster of virus proteomes that was consistent with virus or host taxonomy. We provide BLASTP search and FTP download of Viral RP protein sequences, and a browser to facilitate the visualization of Viral RPs. AVAILABILITY AND IMPLEMENTATION: http://proteininformationresource.org/rps/viruses/ CONTACT: chenc@udel.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Ribonucleotide reductases reveal novel viral diversity and predict biological and ecological features of unknown marine viruses.

Virioplankton play a crucial role in aquatic ecosystems as top-down regulators of bacterial populations and agents of horizontal gene transfer and nutrient cycling. However, the biology and ecology of virioplankton populations in the environment remain poorly understood. Ribonucleotide reductases (RNRs) are ancient enzymes that reduce ribonucleotides to deoxyribonucleotides and thus prime DNA synthesis. Composed of three classes according to O2 reactivity, RNRs can be predictive of the physiological conditions surrounding DNA synthesis. RNRs are universal among cellular life, common within viral genomes and virioplankton shotgun metagenomes (viromes), and estimated to occur within >90% of the dsDNA virioplankton sampled in this study. RNRs occur across diverse viral groups, including all three morphological families of tailed phages, making these genes attractive for studies of viral diversity. Differing patterns in virioplankton diversity were clear from RNRs sampled across a broad oceanic transect. The most abundant RNRs belonged to novel lineages of podoviruses infecting α-proteobacteria, a bacterial class critical to oceanic carbon cycling. RNR class was predictive of phage morphology among cyanophages and RNR distribution frequencies among cyanophages were largely consistent with the predictions of the "kill the winner-cost of resistance" model. RNRs were also identified for the first time to our knowledge within ssDNA viromes. These data indicate that RNR polymorphism provides a means of connecting the biological and ecological features of virioplankton populations.

Polygenic inheritance of cryptorchidism susceptibility in the LE/orl rat.

STUDY HYPOTHESIS: Susceptibility to inherited cryptorchidism in the LE/orl rat may be associated with genetic loci that influence developmental patterning of the gubernaculum by the fetal testis. STUDY FINDING: Cryptorchidism in the LE/orl rat is associated with a unique combination of homozygous minor alleles at multiple loci, and the encoded proteins are co-localized with androgen receptor (AR) and Leydig cells in fetal gubernaculum and testis, respectively. WHAT IS KNOWN ALREADY: Prior studies have shown aberrant perinatal gubernacular migration, muscle patterning defects and reduced fetal testicular testosterone in the LE/orl strain. In addition, altered expression of androgen-responsive, cytoskeletal and muscle-related transcripts in the LE/orl fetal gubernaculum suggest a role for defective AR signaling in cryptorchidism susceptibility. STUDY DESIGN, SAMPLES/MATERIALS, METHODS: The long-term LE/orl colony and short-term colonies of outbred Crl:LE and Crl:SD, and inbred WKY/Ncrl rats were maintained for studies. Animals were intercrossed (LE/orl X WKY/Ncrl), and obligate heterozygotes were reciprocally backcrossed to LE/orl rats to generate 54 F2 males used for genotyping and/or linkage analysis. At least five fetuses per gestational time point from two or more litters were used for quantitative real-time RT-PCR (qRT-PCR) and freshly harvested embryonic (E) day 17 gubernaculum was used to generate conditionally immortalized cell lines. We completed genotyping and gene expression analyses using genome-wide microsatellite markers and single nucleotide polymorphism (SNP) arrays, PCR amplification, direct sequencing, restriction enzyme digest with fragment analysis, whole genome sequencing (WGS), and qRT-PCR. Linkage analysis was performed in Haploview with multiple testing correction, and qRT-PCR data were analyzed using ANOVA after log transformation. Imaging was performed using custom and commercial antibodies directed at candidate proteins in gubernaculum and testis tissues, and gubernaculum cell lines. MAIN RESULTS AND THE ROLE OF CHANCE: LE/orl rats showed reduced fertility and fecundity, and higher risk of perinatal death as compared with Crl:LE rats, but there were no differences in breeding outcomes between normal and unilaterally cryptorchid males. Linkage analysis identified multiple peaks, and with selective breeding of outbred Crl:LE and Crl:SD strains for alleles within two of the most significant (P < 0.003) peaks on chromosomes 6 and 16, we were able to generate a non-LE/orl cryptorchid rat. Associated loci contain potentially functional minor alleles (0.25-0.36 in tested rat strains) including an exonic deletion in Syne2, a large intronic insertion in Ncoa4 (an AR coactivator) and potentially deleterious variants in Solh/Capn15, Ankrd28, and Hsd17b2. Existing WGS data indicate that homozygosity for these combined alleles does not occur in any other sequenced rat strain. We observed a modifying effect of the Syne2(del) allele on expression of other candidate genes, particularly Ncoa4, and for muscle and hormone-responsive transcripts. The selected candidate genes/proteins are highly expressed, androgen-responsive and/or co-localized with developing muscle and AR in fetal gubernaculum, and co-localized with Leydig cells in fetal testis. LIMITATIONS, REASONS FOR CAUTION: The present study identified multiple cryptorchidism-associated linkage peaks in the LE/orl rat, containing potentially causal alleles. These are strong candidate susceptibility loci, but further studies are needed to demonstrate functional relevance to the phenotype. WIDER IMPLICATIONS OF THE FINDINGS: Association data from both human and rat models of spontaneous, nonsyndromic cryptorchidism support a polygenic etiology of the disease. Both the present study and a human genome-wide association study suggest that common variants with weak effects contribute to susceptibility, and may exist in genes encoding proteins that participate in AR signaling in the developing gubernaculum. These findings have potential implications for the gene-environment interaction in the etiology of cryptorchidism. LARGE SCALE DATA: Sequences were deposited in the Rat Genome Database (RGD, http://rgd.mcw.edu/). STUDY FUNDING AND COMPETING INTERESTS: This work was supported by: R01HD060769 from the Eunice Kennedy Shriver National Institute for Child Health and Human Development (NICHD), 2P20GM103446 and P20GM103464 from the National Institute of General Medical Sciences (NIGMS), and Nemours Biomedical Research. The authors have no competing interests to declare.

Direct assessment of viral diversity in soils by random PCR amplification of polymorphic DNA.

Viruses are the most abundant and diverse biological entities within soils, yet their ecological impact is largely unknown. Defining how soil viral communities change with perturbation or across environments will contribute to understanding the larger ecological significance of soil viruses. A new approach to examining the composition of soil viral communities based on random PCR amplification of polymorphic DNA (RAPD-PCR) was developed. A key methodological improvement was the use of viral metagenomic sequence data for the design of RAPD-PCR primers. This metagenomically informed approach to primer design enabled the optimization of RAPD-PCR sensitivity for examining changes in soil viral communities. Initial application of RAPD-PCR viral fingerprinting to soil viral communities demonstrated that the composition of autochthonous soil viral assemblages noticeably changed over a distance of meters along a transect of Antarctic soils and across soils subjected to different land uses. For Antarctic soils, viral assemblages segregated upslope from the edge of dry valley lakes. In the case of temperate soils at the Kellogg Biological Station, viral communities clustered according to land use treatment. In both environments, soil viral communities changed along with environmental factors known to shape the composition of bacterial host communities. Overall, this work demonstrates that RAPD-PCR fingerprinting is an inexpensive, high-throughput means for addressing first-order questions of viral community dynamics within environmental samples and thus fills a methodological gap between narrow single-gene approaches and comprehensive shotgun metagenomic sequencing for the analysis of viral community diversity.

Gallionellaceae pangenomic analysis reveals insight into phylogeny, metabolic flexibility, and iron oxidation mechanisms.

The iron-oxidizing Gallionellaceae drive a wide variety of biogeochemical cycles through their metabolisms and biominerals. To better understand the environmental impacts of Gallionellaceae, we need to improve our knowledge of their diversity and metabolisms, especially any novel iron oxidation mechanisms. Here, we used a pangenomic analysis of 103 genomes to resolve Gallionellaceae phylogeny and explore the range of genomic potential. Using a concatenated ribosomal protein tree and key gene patterns, we determined Gallionellaceae has four genera, divided into two groups-iron-oxidizing bacteria (FeOB) Gallionella, Sideroxydans, and Ferriphaselus with known iron oxidases (Cyc2, MtoA) and nitrite-oxidizing bacteria (NOB) Candidatus Nitrotoga with nitrite oxidase (Nxr). The FeOB and NOB have similar electron transport chains, including genes for reverse electron transport and carbon fixation. Auxiliary energy metabolisms including S oxidation, denitrification, and organotrophy were scattered throughout the Gallionellaceae FeOB. Within FeOB, we found genes that may represent adaptations for iron oxidation, including a variety of extracellular electron uptake (EEU) mechanisms. FeOB genomes encoded more predicted c-type cytochromes overall, notably more multiheme c-type cytochromes (MHCs) with >10 CXXCH motifs. These include homologs of several predicted outer membrane porin-MHC complexes, including MtoAB and Uet. MHCs are known to efficiently conduct electrons across longer distances and function across a wide range of redox potentials that overlap with mineral redox potentials, which can help expand the range of usable iron substrates. Overall, the results of pangenome analyses suggest that the Gallionellaceae genera Gallionella, Sideroxydans, and Ferriphaselus are primarily iron oxidizers, capable of oxidizing dissolved Fe2+ as well as a range of solid iron or other mineral substrates.

Metagenomic Analysis of a Concrete Bridge Reveals a Microbial Community Dominated by Halophilic Bacteria and Archaea.

Concrete hosts a small but diverse microbiome that changes over time. Shotgun metagenomic sequencing would enable assessment of both the diversity and function of the microbial community in concrete, but a number of unique challenges make this difficult for concrete samples. The high concentration of divalent cations in concrete interferes with nucleic acid extraction, and the extremely low biomass in concrete means that DNA from laboratory contamination may be a large fraction of the sequence data. Here, we develop an improved method for DNA extraction from concrete, with higher yield and lower laboratory contamination. To show that this method provides DNA of sufficient quality and quantity to do shotgun metagenomic sequencing, DNA was extracted from a sample of concrete obtained from a road bridge and sequenced with an Illumina MiSeq system. This microbial community was dominated by halophilic Bacteria and Archaea, with enriched functional pathways related to osmotic stress responses. Although this was a pilot-scale effort, we demonstrate that metagenomic sequencing can be used to characterize microbial communities in concrete and that older concrete structures may host different microbes than recently poured concrete. IMPORTANCE Prior work on the microbial communities of concrete focused on the surfaces of concrete structures such as sewage pipes or bridge pilings, where thick biofilms were easy to observe and sample. Because the biomass inside concrete is so low, more recent analyses of the microbial communities inside concrete used amplicon sequencing methods to describe those communities. However, to understand the activity and physiology of microbes in concrete, or to develop living infrastructure, we must develop more direct methods of community analysis. The method developed here for DNA extraction and metagenomic sequencing can be used for analysis of microbial communities inside concrete and can likely be adapted for other cementitious materials.

Gallionellaceae pangenomic analysis reveals insight into phylogeny, metabolic flexibility, and iron oxidation mechanisms.

IMPORTANCE: Neutrophilic iron-oxidizing bacteria (FeOB) produce copious iron (oxyhydr)oxides that can profoundly influence biogeochemical cycles, notably the fate of carbon and many metals. To fully understand environmental microbial iron oxidation, we need a thorough accounting of iron oxidation mechanisms. In this study, we show the Gallionellaceae FeOB genomes encode both characterized iron oxidases as well as uncharacterized multiheme cytochromes (MHCs). MHCs are predicted to transfer electrons from extracellular substrates and likely confer metabolic capabilities that help Gallionellaceae occupy a range of different iron- and mineral-rich niches. Gallionellaceae appear to specialize in iron oxidation, so it would be advantageous for them to have multiple mechanisms to oxidize various forms of iron, given the many iron minerals on Earth, as well as the physiological and kinetic challenges faced by FeOB. The multiple iron/mineral oxidation mechanisms may help drive the widespread ecological success of Gallionellaceae.

miR-31-mediated local translation at the mitotic spindle is important for early development.

miR-31 is a highly conserved microRNA that plays critical roles in cell proliferation, migration, and differentiation. We discovered miR-31 and some of its validated targets are enriched on the mitotic spindle of the dividing sea urchin embryo and mammalian cells. Using the sea urchin embryo, we found that miR-31 inhibition led to developmental delay correlated with increased cytoskeleton and chromosomal defects. We identified miR-31 to directly suppress several actin remodeling transcripts, ß-actin, Gelsolin, Rab35 and Fascin, which were localized to the mitotic spindle. miR-31 inhibition leads to increased newly translated Fascin at the spindles. Forced ectopic localization of Fascin transcripts to the cell membrane and translation led to significant developmental and chromosomal segregation defects, leading to our hypothesis that miR-31 regulates local translation at the mitotic spindle to ensure proper cell division. Furthermore, miR-31-mediated post-transcriptional regulation at the mitotic spindle may be an evolutionarily conserved regulatory paradigm of mitosis.

Functional Annotation Routines Used by ABRF Bioinformatics Core Facilities - Observations, Comparisons, and Considerations.

The functional annotation of gene lists is a common analysis routine required for most genomics experiments, and bioinformatics core facilities must support these analyses. In contrast to methods such as the quantitation of RNA-Seq reads or differential expression analysis, our research group noted a lack of consensus in our preferred approaches to functional annotation. To investigate this observation, we selected 4 experiments that represent a range of experimental designs encountered by our cores and analyzed those data with 6 tools used by members of the Association of Biomolecular Resource Facilities (ABRF) Genomic Bioinformatics Research Group (GBIRG). To facilitate comparisons between tools, we focused on a single biological result for each experiment. These results were represented by a gene set, and we analyzed these gene sets with each tool considered in our study to map the result to the annotation categories presented by each tool. In most cases, each tool produces data that would facilitate identification of the selected biological result for each experiment. For the exceptions, Fisher's exact test parameters could be adjusted to detect the result. Because Fisher's exact test is used by many functional annotation tools, we investigated input parameters and demonstrate that, while background set size is unlikely to have a significant impact on the results, the numbers of differentially expressed genes in an annotation category and the total number of differentially expressed genes under consideration are both critical parameters that may need to be modified during analyses. In addition, we note that differences in the annotation categories tested by each tool, as well as the composition of those categories, can have a significant impact on results.

Spontaneously Produced Lysogenic Phages Are an Important Component of the Soybean Bradyrhizobium Mobilome.

The ability of Bradyrhizobium spp. to nodulate and fix atmospheric nitrogen in soybean root nodules is critical to meeting humanity's nutritional needs. The intricacies of soybean bradyrhizobia-plant interactions have been studied extensively; however, bradyrhizobial ecology as influenced by phages has received somewhat less attention, even though these interactions may significantly impact soybean yield. In batch culture, four soybean bradyrhizobia strains, Bradyrhizobium japonicum S06B (S06B-Bj), B. japonicum S10J (S10J-Bj), Bradyrhizobium diazoefficiens USDA 122 (USDA 122-Bd), and Bradyrhizobium elkanii USDA 76T (USDA 76-Be), spontaneously (without apparent exogenous chemical or physical induction) produced tailed phages throughout the growth cycle; for three strains, phage concentrations exceeded cell numbers by ~3-fold after 48 h of incubation. Phage terminase large-subunit protein phylogeny revealed possible differences in phage packaging and replication mechanisms. Bioinformatic analyses predicted multiple prophage regions within each soybean bradyrhizobia genome, preventing accurate identification of spontaneously produced prophage (SPP) genomes. A DNA sequencing and mapping approach accurately delineated the boundaries of four SPP genomes within three of the soybean bradyrhizobia chromosomes and suggested that the SPPs were capable of transduction. In addition to the phages, S06B-Bj and USDA 76-Be contained three to four times more insertion sequences (IS) and large, conjugable, broad host range plasmids, both of which are known drivers of horizontal gene transfer (HGT) in soybean bradyrhizobia. These factors indicate that SPP along with IS and plasmids participate in HGT, drive bradyrhizobia evolution, and play an outsized role in bradyrhizobia ecology. IMPORTANCE Previous studies have shown that IS and plasmids mediate HGT of symbiotic nodulation (nod) genes in soybean bradyrhizobia; however, these events require close cell-to-cell contact, which could be limited in soil environments. Bacteriophage-assisted gene transduction through spontaneously produced prophages provides a stable means of HGT not limited by the constraints of proximal cell-to-cell contact. These phage-mediated HGT events may shape soybean bradyrhizobia population ecology, with concomitant impacts on soybean agriculture.

Ubiquitous, B12-dependent virioplankton utilizing ribonucleotide-triphosphate reductase demonstrate interseasonal dynamics and associate with a diverse range of bacterial hosts in the pelagic ocean.

Through infection and lysis of their coexisting bacterial hosts, viruses impact the biogeochemical cycles sustaining globally significant pelagic oceanic ecosystems. Currently, little is known of the ecological interactions between lytic viruses and their bacterial hosts underlying these biogeochemical impacts at ecosystem scales. This study focused on populations of lytic viruses carrying the B12-dependent Class II monomeric ribonucleotide reductase (RNR) gene, ribonucleotide-triphosphate reductase (Class II RTPR), documenting seasonal changes in pelagic virioplankton and bacterioplankton using amplicon sequences of Class II RTPR and the 16S rRNA gene, respectively. Amplicon sequence libraries were analyzed using compositional data analysis tools that account for the compositional nature of these data. Both virio- and bacterioplankton communities responded to environmental changes typically seen across seasonal cycles as well as shorter term upwelling-downwelling events. Defining Class II RTPR-carrying viral populations according to major phylogenetic clades proved a more robust means of exploring virioplankton ecology than operational taxonomic units defined by percent sequence homology. Virioplankton Class II RTPR populations showed positive associations with a broad phylogenetic diversity of bacterioplankton including dominant taxa within pelagic oceanic ecosystems such as Prochlorococcus and SAR11. Temporal changes in Class II RTPR virioplankton, occurring as both free viruses and within infected cells, indicated possible viral-host pairs undergoing sustained infection and lysis cycles throughout the seasonal study. Phylogenetic relationships inferred from Class II RTPR sequences mirrored ecological patterns in virio- and bacterioplankton populations demonstrating possible genome to phenome associations for an essential viral replication gene.