Search

For much of terrestrial biodiversity, the evolutionary pathways of adaptation from marine ancestors are poorly understood and have usually been viewed as a binary trait. True crabs, the decapod crustacean infraorder Brachyura, comprise over 7600 species representing a striking diversity of morphology and ecology, including repeated adaptation to non-marine habitats. Here, we reconstruct the evolutionary history of Brachyura using new and published sequences of 10 genes for 344 tips spanning 88 of 109 brachyuran families. Using 36 newly vetted fossil calibrations, we infer that brachyurans most likely diverged in the Triassic, with family-level splits in the late Cretaceous and early Paleogene. By contrast, the root age is underestimated with automated sampling of 328 fossil occurrences explicitly incorporated into the tree prior, suggesting such models are a poor fit under heterogeneous fossil preservation. We apply recently defined trait-by-environment associations to classify a gradient of transitions from marine to terrestrial lifestyles. We estimate that crabs left the marine environment at least 7 and up to 17 times convergently, and returned to the sea from non-marine environments at least twice. Although the most highly terrestrial- and many freshwater-adapted crabs are concentrated in Thoracotremata, Bayesian threshold models of ancestral state reconstruction fail to identify shifts to higher terrestrial grades due to the degree of underlying change required. Lineages throughout our tree inhabit intertidal and marginal marine environments, corroborating the inference that the early stages of terrestrial adaptation have a lower threshold to evolve. Our framework and extensive new fossil and natural history datasets will enable future comparisons of non-marine adaptation at the morphological and molecular level. Crabs provide an important window into the early processes of adaptation to novel environments, and different degrees of evolutionary constraint that might help predict these pathways. [Brachyura; convergent evolution; crustaceans; divergence times; fossil calibration; molecular phylogeny; terrestrialization; threshold model.]., (© The Author(s) 2023. Published by Oxford University Press on behalf of the Society of Systematic Biologists.)

Escherichia coli is the leading cause of urinary tract infections (UTIs) in children and adults. The gastrointestinal tract is the primary reservoir of uropathogenic E. coli , which can be acquired from a variety of environmental exposures, including retail meat. In the current study, we used a novel statistical-genomic approach to estimate the proportion of pediatric UTIs caused by foodborne zoonotic E. coli strains. E. coli urine isolates were collected from DC residents aged 2 months to 17 years from the Children's National Medical Center Laboratory, 2013-2014. During the same period, E. coli isolates were collected from retail poultry products purchased from 15 sites throughout DC. A total of 52 urine and 56 poultry isolates underwent whole-genome sequencing, core genome phylogenetic analysis, and host-origin prediction by a Bayesian latent class model that incorporated data on the presence of mobile genetic elements (MGEs) among E. coli isolates from multiple vertebrate hosts. A total of 56 multilocus sequence types were identified among the isolates. Five sequence types-ST10, ST38, ST69, ST117, and ST131-were observed among both urine and poultry isolates. Using the Bayesian latent class model, we estimated that 19% (10/52) of the clinical E. coli isolates in our population were foodborne zoonotic strains. These data suggest that a substantial portion of pediatric UTIs in the Washington DC region may be caused by E. coli strains originating in food animals and likely transmitted via contaminated poultry meat.IMPORTANCE Escherichia coli UTIs are a heavy public health burden and can have long-term negative health consequences for pediatric patients. E. coli has an extremely broad host range, including humans, chickens, turkeys, pigs, and cattle. E. coli derived from food animals is a frequent contaminant of retail meat products, but little is known about the risk these strains pose to pediatric populations. Quantifying the proportion of pediatric UTIs caused by food-animal-derived E. coli , characterizing the highest-risk strains, and identifying their primary reservoir species could inform novel intervention strategies to reduce UTI burden in this vulnerable population. Our results suggest that retail poultry meat may be an important vehicle for pediatric exposure to zoonotic E. coli strains capable of causing UTIs. Vaccinating poultry against the highest-risk strains could potentially reduce poultry colonization, poultry meat contamination, and downstream pediatric infections., Competing Interests: The authors declare no conflict of interest.

Given the multitude of challenges Earth is facing, sustainability science is of key importance to our continued existence. Evolution is the fundamental biological process underlying the origin of all biodiversity. This phylogenetic diversity fosters the resilience of ecosystems to environmental change, and provides numerous resources to society, and options for the future. Genetic diversity within species is also key to the ability of populations to evolve and adapt to environmental change. Yet, the value of evolutionary processes and the consequences of their impairment have not generally been considered in sustainability research. We argue that biological evolution is important for sustainability and that the concepts, theory, data, and methodological approaches used in evolutionary biology can, in crucial ways, contribute to achieving the UN Sustainable Development Goals (SDGs). We discuss how evolutionary principles are relevant to understanding, maintaining, and improving Nature Contributions to People (NCP) and how they contribute to the SDGs. We highlight specific applications of evolution, evolutionary theory, and evolutionary biology's diverse toolbox, grouped into four major routes through which evolution and evolutionary insights can impact sustainability. We argue that information on both within-species evolutionary potential and among-species phylogenetic diversity is necessary to predict population, community, and ecosystem responses to global change and to make informed decisions on sustainable production, health, and well-being. We provide examples of how evolutionary insights and the tools developed by evolutionary biology can not only inspire and enhance progress on the trajectory to sustainability, but also highlight some obstacles that hitherto seem to have impeded an efficient uptake of evolutionary insights in sustainability research and actions to sustain SDGs. We call for enhanced collaboration between sustainability science and evolutionary biology to understand how integrating these disciplines can help achieve the sustainable future envisioned by the UN SDGs., Competing Interests: The authors declare that they do not have any conflict of interest., (© 2024 The Author(s). Evolutionary Applications published by John Wiley & Sons Ltd.)

The white spot syndrome virus (WSSV) is a causative agent of white spot disease (WSD) in crustaceans, especially in cultivated black tiger shrimp ( Penaeus monodon ), leading to significant economic losses in the aquaculture sector. The present study describes four whole genome sequences of WSSV obtained from coastal regions of Bangladesh., Competing Interests: The authors declare no conflict of interest.

Introduction: During the COVID-19 Delta variant surge, the CLAIRE cross-sectional study sampled saliva from 120 hospitalized patients, 116 of whom had a positive COVID-19 PCR test. Patients received antibiotics upon admission due to possible secondary bacterial infections, with patients at risk of sepsis receiving broad-spectrum antibiotics (BSA)., Methods: The saliva samples were analyzed with shotgun DNA metagenomics and respiratory RNA virome sequencing. Medical records for the period of hospitalization were obtained for all patients. Once hospitalization outcomes were known, patients were classified based on their COVID-19 disease severity and the antibiotics they received., Results: Our study reveals that BSA regimens differentially impacted the human salivary microbiome and disease progression. 12 patients died and all of them received BSA. Significant associations were found between the composition of the COVID-19 saliva microbiome and BSA use, between SARS-CoV-2 genome coverage and severity of disease. We also found significant associations between the non-bacterial microbiome and severity of disease, with Candida albicans detected most frequently in critical patients. For patients who did not receive BSA before saliva sampling, our study suggests Staphylococcus aureus as a potential risk factor for sepsis., Discussion: Our results indicate that the course of the infection may be explained by both monitoring antibiotic treatment and profiling a patient's salivary microbiome, establishing a compelling link between microbiome and the specific antibiotic type and timing of treatment. This approach can aid with emergency room triage and inpatient management but also requires a better understanding of and access to narrow-spectrum agents that target pathogenic bacteria., Competing Interests: PB and KF were employed by Lifetime Omics. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The authors declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision., (Copyright © 2024 Buendia, Fernandez, Raley, Rahnavard, Crandall and Castro.)

Background: Predicting phenotypes from genetic variation is foundational for fields as diverse as bioengineering and global change biology, highlighting the importance of efficient methods to predict gene functions. Linking genetic changes to phenotypic changes has been a goal of decades of experimental work, especially for some model gene families, including light-sensitive opsin proteins. Opsins can be expressed in vitro to measure light absorption parameters, including λmax-the wavelength of maximum absorbance-which strongly affects organismal phenotypes like color vision. Despite extensive research on opsins, the data remain dispersed, uncompiled, and often challenging to access, thereby precluding systematic and comprehensive analyses of the intricate relationships between genotype and phenotype., Results: Here, we report a newly compiled database of all heterologously expressed opsin genes with λmax phenotypes that we call the Visual Physiology Opsin Database (VPOD). VPOD_1.0 contains 864 unique opsin genotypes and corresponding λmax phenotypes collected across all animals from 73 separate publications. We use VPOD data and deepBreaks to show regression-based machine learning (ML) models often reliably predict λmax, account for nonadditive effects of mutations on function, and identify functionally critical amino acid sites., Conclusion: The ability to reliably predict functions from gene sequences alone using ML will allow robust exploration of molecular-evolutionary patterns governing phenotype, will inform functional and evolutionary connections to an organism's ecological niche, and may be used more broadly for de novo protein design. Together, our database, phenotype predictions, and model comparisons lay the groundwork for future research applicable to families of genes with quantifiable and comparable phenotypes., (© The Author(s) 2024. Published by Oxford University Press GigaScience.)

Sucralose and acesulfame-potassium consumption alters gut microbiota in rodents, with unclear effects in humans. We examined effects of three-times daily sucralose- and acesulfame-potassium-containing diet soda consumption for 1 ( n  = 17) or 8 ( n  = 8) weeks on gut microbiota composition in young adults. After 8 weeks of diet soda consumption, the relative abundance of Proteobacteria, specifically Enterobacteriaceae , increased; and, increased abundance of two Proteobacteria taxa was also observed after 1 week of diet soda consumption compared with sparkling water. In addition, three taxa in the Bacteroides genus increased following 1 week of diet soda consumption compared with sparkling water. The clinical relevance of these findings and effects of sucralose and acesulfame-potassium consumption on human gut microbiota warrant further investigation in larger studies. Clinical trial registration: NCT02877186 and NCT03125356., Competing Interests: None of the authors have any conflicts of interest to disclose.

The bacterial communities of the human skin impact its physiology and homeostasis, hence elucidating the composition and structure of the healthy skin bacteriome is paramount to understand how bacterial imbalance (i.e., dysbiosis) may lead to disease. To obtain an integrated view of the spatial diversity of the skin bacteriome, we surveyed from 2019 to 2023 five skin regions (belly button, behind ears, between toes, calves and forearms) with different physiological characteristics (dry, moist and sebaceous) in 129 healthy adults (579 samples - after data cleaning). Estimating bacterial diversity through 16S rRNA metataxonomics, we identified significant ( p  < 0.0001) differences in the bacterial relative abundance of the four most abundant phyla and 11 genera, alpha- and beta-diversity indices and predicted functional profiles (36 to 400 metabolic pathways) across skin regions and microenvironments. No significant differences, however, were observed across genders, ages, and ethnicities. As previously suggested, dry skin regions (forearms and calves) were more even, richer, and functionally distinct than sebaceous (behind ears) and moist (belly button and between toes) regions. Within skin regions, bacterial alpha- and beta-diversity also varied significantly for some of the years compared, suggesting that skin bacterial stability may be region and subject dependent. Our results, hence, confirm that the skin bacteriome varies systematically across skin regions and microenvironments and provides new insights into the internal and external factors driving bacterial diversity., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision., (Copyright © 2023 Pérez-Losada and Crandall.)

Macrophage-lineage cells are indispensable to immunity and physiology of all vertebrates. Amongst these, amphibians represent a key stage in vertebrate evolution and are facing decimating population declines and extinctions, in large part due to emerging infectious agents. While recent studies indicate that macrophages and related innate immune cells are critically involved during these infections, much remains unknown regarding the ontogeny and functional differentiation of these cell types in amphibians. Accordingly, in this review we coalesce what has been established to date about amphibian blood cell development (hematopoiesis), the development of key amphibian innate immune cells (myelopoiesis) and the differentiation of amphibian macrophage subsets (monopoiesis). We explore the current understanding of designated sites of larval and adult hematopoiesis across distinct amphibian species and consider what mechanisms may lend to these species-specific adaptations. We discern the identified molecular mechanisms governing the functional differentiation of disparate amphibian (chiefly Xenopus laevis) macrophage subsets and describe what is known about the roles of these subsets during amphibian infections with intracellular pathogens. Macrophage lineage cells are at the heart of so many vertebrate physiological processes. Thus, garnering greater understanding of the mechanisms responsible for the ontogeny and functionality of these cells in amphibians will lend to a more comprehensive view of vertebrate evolution., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

Most experiments studying bacterial microbiomes rely on the PCR amplification of all or part of the gene for the 16S rRNA subunit, which serves as a biomarker for identifying and quantifying the various taxa present in a microbiome sample. Several computational methods exist for analyzing 16S amplicon sequencing. However, the most-used bioinformatics tools cannot produce high quality genus-level or species-level taxonomic calls and may underestimate the potential accuracy of these calls. We used 16S sequencing data from mock bacterial communities to evaluate the sensitivity and specificity of several bioinformatics pipelines and genomic reference libraries used for microbiome analyses, concentrating on measuring the accuracy of species-level taxonomic assignments of 16S amplicon reads. We evaluated the tools DADA2, QIIME 2, Mothur, PathoScope 2, and Kraken 2 in conjunction with reference libraries from Greengenes, SILVA, Kraken 2, and RefSeq. Profiling tools were compared using publicly available mock community data from several sources, comprising 136 samples with varied species richness and evenness, several different amplified regions within the 16S rRNA gene, and both DNA spike-ins and cDNA from collections of plated cells. PathoScope 2 and Kraken 2, both tools designed for whole-genome metagenomics, outperformed DADA2, QIIME 2 using the DADA2 plugin, and Mothur, which are theoretically specialized for 16S analyses. Evaluations of reference libraries identified the SILVA and RefSeq/Kraken 2 Standard libraries as superior in accuracy compared to Greengenes. These findings support PathoScope and Kraken 2 as fully capable, competitive options for genus- and species-level 16S amplicon sequencing data analysis, whole genome sequencing, and metagenomics data tools., (© 2023. Springer Nature Limited.)

The clade Pancrustacea, comprising crustaceans and hexapods, is the most diverse group of animals on earth, containing over 80% of animal species and half of animal biomass. It has been the subject of several recent phylogenomic analyses, yet relationships within Pancrustacea show a notable lack of stability. Here, the phylogeny is estimated with expanded taxon sampling, particularly of malacostracans. We show small changes in taxon sampling have large impacts on phylogenetic estimation. By analyzing identical orthologs between two slightly different taxon sets, we show that the differences in the resulting topologies are due primarily to the effects of taxon sampling on the phylogenetic reconstruction method. We compare trees resulting from our phylogenomic analyses with those from the literature to explore the large tree space of pancrustacean phylogenetic hypotheses and find that statistical topology tests reject the previously published trees in favor of the maximum likelihood trees produced here. Our results reject several clades including Caridoida, Eucarida, Multicrustacea, Vericrustacea, and Syncarida. Notably, we find Copepoda nested within Allotriocarida with high support and recover a novel relationship between decapods, euphausiids, and syncarids that we refer to as the Syneucarida. With denser taxon sampling, we find Stomatopoda sister to this latter clade, which we collectively name Stomatocarida, dividing Malacostraca into three clades: Leptostraca, Peracarida, and Stomatocarida. A new Bayesian divergence time estimation is conducted using 13 vetted fossils. We review our results in the context of other pancrustacean phylogenetic hypotheses and highlight 15 key taxa to sample in future studies., (Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution 2023.)

Transposable elements (TE) are typically silenced by DNA methylation and repressive histone modifications in differentiated healthy human tissues. However, TE expression increases in a wide range of cancers and is correlated with global hypomethylation of cancer genomes. We assessed expression and DNA methylation of TEs in fibroblast cells that were serially transduced with hTERT, SV40, and HRASR24C to immortalize and then transform them, modeling the different steps of the tumorigenesis process. RNA sequencing and whole-genome bisulfite sequencing were performed at each stage of transformation. TE expression significantly increased as cells progressed through transformation, with the largest increase in expression after the final stage of transformation, consistent with data from human tumors. The upregulated TEs were dominated by endogenous retroviruses [long terminal repeats (LTR)]. Most differentially methylated regions (DMR) in all stages were hypomethylated, with the greatest hypomethylation in the final stage of transformation. A majority of the DMRs overlapped TEs from the RepeatMasker database, indicating that TEs are preferentially demethylated. Many hypomethylated TEs displayed a concordant increase in expression. Demethylation began during immortalization and continued into transformation, while upregulation of TE transcription occurred in transformation. Numerous LTR elements upregulated in the model were also identified in The Cancer Genome Atlas datasets of breast, colon, and prostate cancer. Overall, these findings indicate that TEs, specifically endogenous retroviruses, are demethylated and transcribed during transformation., Significance: Analysis of epigenetic and transcriptional changes in a transformation model reveals that transposable element expression and methylation are dysregulated during oncogenic transformation., (©2023 American Association for Cancer Research.)

Introduction: Human endogenous retroviruses (HERVs) are transcriptionally-active remnants of ancient retroviral infections that may play a role in Alzheimer's disease., Methods: We combined two, publicly available RNA-Seq datasets with a third, novel dataset for a total cohort of 103 patients with Alzheimer's disease and 45 healthy controls. We use telescope to perform HERV quantification for these samples and simultaneously perform gene expression analysis., Results: We identify differentially expressed genes and differentially expressed HERVs in Alzheimer's disease patients. Differentially expressed HERVs are scattered throughout the genome; many of them are members of the HERV-K superfamily. A number of HERVs are correlated with the expression of dysregulated genes in Alzheimer's and are physically proximal to genes which drive disease pathways., Discussion: Dysregulated expression of ancient retroviral insertions in the human genome are present in Alzheimer's disease and show localization patterns that may explain how these elements drive pathogenic gene expression., Competing Interests: CC has received research support from: GlaxoSmithKline and Eisai. The funders of the study had no role in the collection, analysis, or interpretation of data; in the writing of the report; or in the decision to submit the manuscript for publication. CC is a member of the advisory board of Vivid Genomics and Circular Genomics and owns stocks. TD has received research support from AMPEL BioSolutions. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Dawson, Rentia, Sanford, Cruchaga, Kauwe and Crandall.)

Raccoon roundworm (Baylisascaris procyonis) negatively affects woodrat (Neotoma spp.) populations but is not known to occur in the endemic range of endangered Key Largo woodrats (Neotoma floridana smalli). Rectal swabs from 23 raccoons (Procyon lotor) in Key Largo were screened for raccoon roundworm by PCR. All tests were negative, suggesting continued absence., (© Wildlife Disease Association 2023.)

Regulation of intron retention (IR), a form of alternative splicing, is a newly recognized checkpoint in gene expression. Since there are numerous abnormalities in gene expression in the prototypic autoimmune disease systemic lupus erythematosus (SLE), we sought to determine whether IR was intact in patients with this disease. We, therefore, studied global gene expression and IR patterns of lymphocytes in SLE patients. We analyzed RNA-seq data from peripheral blood T cell samples from 14 patients suffering from systemic lupus erythematosus (SLE) and 4 healthy controls and a second, independent data set of RNA-seq data from B cells from16 SLE patients and 4 healthy controls. We identified intron retention levels from 26,372 well annotated genes as well as differential gene expression and tested for differences between cases and controls using unbiased hierarchical clustering and principal component analysis. We followed with gene-disease enrichment analysis and gene-ontology enrichment analysis. Finally, we then tested for significant differences in intron retention between cases and controls both globally and with respect to specific genes. Overall decreased IR was found in T cells from one cohort and B cells from another cohort of patients with SLE and was associated with increased expression of numerous genes, including those encoding spliceosome components. Different introns within the same gene displayed both up- and down-regulated retention profiles indicating a complex regulatory mechanism. These results indicate that decreased IR in immune cells is characteristic of patients with active SLE and may contribute to the abnormal expression of specific genes in this autoimmune disease., (© 2023. The Author(s).)

Persons with cystic fibrosis (PwCF) suffer from pulmonary exacerbations (PEx) related in part to lung infection. While higher microbial diversity is associated with higher lung function, the data on the impact of short-term antibiotics on changes in microbial diversity is conflicting. Further, Prevotella secretes beta-lactamases, which may influence recovery of lung function. We hypothesize that sub-therapeutic and broad spectrum antibiotic exposure leads to decreasing microbial diversity. Our secondary aim was to evaluate the concerted association of beta-lactam pharmacokinetics (PK), antibiotic spectrum, microbial diversity, and antibiotic resistance on lung function recovery using a pathway analysis. This was a retrospective observational study of persons with CF treated with IV antibiotics for PEx between 2016 and 2020 at Children's National Hospital; respiratory samples and clinical information were collected at hospital admission for PEx (E), end of antibiotic treatment (T), and follow-up (F). Metagenomic sequencing was performed; PathoScope 2.0 and AmrPlusPlus were used for taxonomic assignment of sequences to bacteria and antibiotic resistance genes (ARGs). M/W Pharm was used for PK modeling. Comparison of categorical and continuous variables and pathway analysis were performed in STATA. Twenty-two PwCF experienced 43 PEx. The study cohort had a mean age of 14.6 years. Only 12/43 beta-lactam courses had therapeutic PK, and 18/43 were broad spectrum. A larger decrease in richness between E and T was seen in the therapeutic PK group (sufficient - 20.1 vs. insufficient - 1.59, p = 0.025) and those receiving broad spectrum antibiotics (broad - 14.5 vs. narrow - 2.8, p = 0.030). We did not detect differences in the increase in percent predicted forced expiratory volume in one second (ppFEV1) at end of treatment compared to PEx based on beta-lactam PK (sufficient 13.6% vs. insufficient 15.1%) or antibiotic spectrum (broad 11.5% vs. narrow 16.6%). While both therapeutic beta-lactam PK and broad-spectrum antibiotics decreased richness between PEx and the end of treatment, we did not detect longstanding changes in alpha diversity or an association with superior recovery of lung function compared with subtherapeutic PK and narrow spectrum antimicrobials., (© 2023. The Author(s).)

Introduction: Pulmonary exacerbations (PEx) in persons with cystic fibrosis (CF) are primarily related to acute or chronic inflammation associated with bacterial lung infections, which may be caused by several bacteria that activate similar bacterial genes and produce similar by-products. The goal of our study was to perform a stratified functional analysis of bacterial genes at three distinct time points in the treatment of a PEx in order to determine the role that specific airway microbiome community members may play within each clinical state (i.e., PEx, end of antibiotic treatment, and follow-up). Our secondary goal was to compare the change between clinical states with the metabolic activity of specific airway microbiome community members., Methods: This was a prospective observational study of persons with CF treated with intravenous antibiotics for PEx between 2016 and 2020 at Children's National Hospital. Demographic and clinical information as well as respiratory samples were collected at hospital admission for PEx, end of antibiotic treatment, and follow-up. Metagenomic sequencing was performed; MetaPhlAn3 and HUMANn3 were used to assign sequences to bacterial species and bacterial metabolic genes, respectively., Results: Twenty-two persons with CF, with a mean age of 14.5 (range 7-23) years, experienced 45 PEx during the study period. Two-hundred twenty-one bacterial species were identified in the respiratory samples from the study cohort. Ten bacterial species had differential gene abundance across changes in the clinical state including Staphylococcus aureus , Streptococcus salivarius , and Veillonella atypica (all padj < 0.01 and log2FoldChange > |2|). These corresponded to a differential abundance of bacterial genes, with S. aureus accounting for 81% of the genes more abundant in PEx and S. salivarius accounting for 83% of the genes more abundant in follow-up, all compared to the end of treatment. Lastly, 8,653 metabolic pathways were identified across samples, with again S. aureus and S. salivarius contributing to the differential abundance of pathways (106 in PEx vs. 66 in follow-up, respectively). V. atypica was associated with a single metabolic pathway (UDP- N -acetyl-D-glucosamine biosynthesis) increased in follow-up compared to PEx., Discussion: Taken together, these data suggest that the metabolic potential of bacterial species can provide more insight into changes across clinical states than the relative abundance of the bacteria alone., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Shumyatsky, Burrell, Chaney, Sami, Koumbourlis, Freishtat, Crandall, Zemanick and Hahn.)

Background: Molecular epidemiological approaches provide opportunities to characterize HIV transmission dynamics. We analyzed HIV sequences and virus load (VL) results obtained during routine clinical care, and individual’s zip-code location to determine utility of this approach. Methods: HIV-1 pol sequences aligned using ClustalW were subtyped using REGA. A maximum likelihood (ML) tree was generated using IQTree. Transmission clusters with ≤3% genetic distance (GD) and ≥90% bootstrap support were identified using ClusterPicker. We conducted Bayesian analysis using BEAST to confirm transmission clusters. The proportion of nucleotides with ambiguity ≤0.5% was considered indicative of early infection. Descriptive statistics were applied to characterize clusters and group comparisons were performed using chi-square or t-test. Results: Among 2775 adults with data from 2014−2015, 2589 (93%) had subtype B HIV-1, mean age was 44 years (SD 12.7), 66.4% were male, and 25% had nucleotide ambiguity ≤0.5. There were 456 individuals in 193 clusters: 149 dyads, 32 triads, and 12 groups with ≥ four individuals per cluster. More commonly in clusters were males than females, 349 (76.5%) vs. 107 (23.5%), p < 0.0001; younger individuals, 35.3 years (SD 12.1) vs. 44.7 (SD 12.3), p < 0.0001; and those with early HIV-1 infection by nucleotide ambiguity, 202/456 (44.3%) vs. 442/2133 (20.7%), p < 0.0001. Members of 43/193 (22.3%) of clusters included individuals in different jurisdictions. Clusters ≥ four individuals were similarly found using BEAST. HIV-1 viral load (VL) ≥3.0 log10 c/mL was most common among individuals in clusters ≥ four, 18/21, (85.7%) compared to 137/208 (65.8%) in clusters sized 2−3, and 927/1169 (79.3%) who were not in a cluster (p < 0.0001). Discussion: HIV sequence data obtained for HIV clinical management provide insights into regional transmission dynamics. Our findings demonstrate the additional utility of HIV-1 VL data in combination with phylogenetic inferences as an enhanced contact tracing tool to direct HIV treatment and prevention services. Trans-jurisdictional approaches are needed to optimize efforts to end the HIV epidemic.

Fecal microbiota transplantation (FMT) is a promising therapeutic modality for the treatment and prevention of metabolic disease. We previously conducted a double-blind, randomized, placebo-controlled pilot trial of FMT in obese metabolically healthy patients in which we found that FMT enhanced gut bacterial bile acid metabolism and delayed the development of impaired glucose tolerance relative to the placebo control group. Therefore, we conducted a secondary analysis of fecal samples collected from these patients to assess the potential gut microbial species contributing to the effect of FMT to improve metabolic health and increase gut bacterial bile acid metabolism. Fecal samples collected at baseline and after 4 weeks of FMT or placebo treatment underwent shotgun metagenomic analysis. Ultra-high-performance liquid chromatography-mass spectrometry was used to profile fecal bile acids. FMT-enriched bacteria that have been implicated in gut bile acid metabolism included Desulfovibrio fairfieldensis and Clostridium hylemonae . To identify candidate bacteria involved in gut microbial bile acid metabolism, we assessed correlations between bacterial species abundance and bile acid profile, with a focus on bile acid products of gut bacterial metabolism. Bacteroides ovatus and Phocaeicola dorei were positively correlated with unconjugated bile acids. Bifidobacterium adolescentis , Collinsella aerofaciens , and Faecalibacterium prausnitzii were positively correlated with secondary bile acids. Together, these data identify several candidate bacteria that may contribute to the metabolic benefits of FMT and gut bacterial bile acid metabolism that requires further functional validation.

An important outcome for veterinary education is ensuring that graduates can provide an appropriate level of care to patients and clients by demonstrating core competencies such as communication skills. In addition, accreditation requirements dictate the need to assess learning outcomes and may drive the motivation to incorporate relevant and appropriate methods of entry assessments for incoming students. Predicting the success of Doctor of Veterinary Medicine (DVM) students based on entry assessment performance has been scantly investigated and can be challenging. Specifically, no research presently exists on predicting DVM students' first-year performance in relation to communication skills at the time of program entry. Objectives of this exploratory study were to investigate (a) the relationship between communication skills outcomes from multiple mini-interview (MMI) data and first-year academic performance related to communication and (b) the relationship between communication skills outcomes from MMI data and self-reported first-year communication reflections. A retrospective single-class study was conducted. Data were analyzed using descriptive statistics, correlation statistics, regression models, and paired t -tests to identify relationships among variables. Paired t -tests showed that students felt more prepared to meet second-year expectations over first-year expectations. Spearman's correlation revealed an association between MMI communication scores and one pre-year 1 survey question related to professionalism. Noo relationships were observed between MMI communication scores and marks from a self-reflection assignment in a communications course, or grades from a clinical medicine course that included clinical communication. The merit for further exploration of the relationship between communication competencies and student performance is discussed.

Background: Cystic fibrosis (CF) is characterized by recurrent pulmonary exacerbations (PEx) and lung function decline. PEx are frequently treated with antibiotics. However, little is known about the effects of antibiotics on the airway microbiome of persons with CF over time. The purpose of this study was to evaluate changes in the microbiome and lung function in persons with CF over 1 year following an initial study pulmonary exacerbation (iPEx)., Methods: Twenty children aged ≤18 years with CF were enrolled in the study, which occurred prior to the routine administration of highly effective modulator therapy. Respiratory samples and spirometry were obtained at a minimum of quarterly visits and up to 1 year after an iPEx. Metagenomic sequencing was performed, and bacterial taxa were assigned using MetaPhlAn 2.0. Paired t test, analysis of variance, and generalized least squares regression were used to compare outcome variables., Results: The mean age of study participants at the time of the iPEx was 10.6 years. There were 3 ± 1.6 PEx treated with antibiotics per person during the study period. Bacterial richness was similar at 1 year compared to iPEx (40.3 vs 39.3, P  = .852), whereas the mean Shannon diversity index was significantly higher at 1 year (2.84 vs 1.62, P  < .001). The number of PEx treated with antibiotics was not associated with changes in microbial diversity but was associated with changes in lung function., Conclusions: In our 1-year prospective study, we found that microbial diversity increased despite decreases in lung function associated with repeated PEx events requiring antibiotic therapy., Competing Interests: Potential conflicts of interest. All authors: No reported conflicts., (© The Author(s) 2022. Published by Oxford University Press on behalf of Infectious Diseases Society of America.)

Rhizosphere microbial communities exert critical roles in plant health, nutrient cycling, and soil fertility. Despite the essential functions conferred by microbes, the source and acquisition of the rhizosphere are not entirely clear. Therefore, we investigated microbial community diversity and potential source using the only two native Antarctic plants, Deschampsia antarctica (Da) and Colobanthus quitensis (Cq), as models. We interrogated rhizosphere and bulk soil microbiomes at six locations in the Byers Peninsula, Livingston Island, Antarctica, both individual plant species and their association (Da.Cq). Our results show that host plant species influenced the richness and diversity of bacterial communities in the rhizosphere. Here, the Da rhizosphere showed the lowest richness and diversity of bacteria compared to Cq and Da.Cq rhizospheres. In contrast, for rhizosphere fungal communities, plant species only influenced diversity, whereas the rhizosphere of Da exhibited higher fungal diversity than the Cq rhizosphere. Also, we found that environmental geographic pressures (i.e., sampling site, latitude, and altitude) and, to a lesser extent, biotic factors (i.e., plant species) determined the species turnover between microbial communities. Moreover, our analysis shows that the sources of the bacterial communities in the rhizosphere were local soils that contributed to homogenizing the community composition of the different plant species growing in the same sampling site. In contrast, the sources of rhizosphere fungi were local (for Da and Da.Cq) and distant soils (for Cq). Here, the host plant species have a specific effect in acquiring fungal communities to the rhizosphere. However, the contribution of unknown sources to the fungal rhizosphere (especially in Da and Da.Cq) indicates the existence of relevant stochastic processes in acquiring these microbes. Our study shows that rhizosphere microbial communities differ in their composition and diversity. These differences are explained mainly by the microbial composition of the soils that harbor them, acting together with plant species-specific effects. Both plant species acquire bacteria from local soils to form part of their rhizosphere. Seemingly, the acquisition process is more complex for fungi. We identified a significant contribution from unknown fungal sources due to stochastic processes and known sources from soils across the Byers Peninsula., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Guajardo-Leiva, Alarcón, Gutzwiller, Gallardo-Cerda, Acuña-Rodríguez, Molina-Montenegro, Crandall, Pérez-Losada and Castro-Nallar.)

Host evolutionary history is a key factor shaping the earthworm cast microbiome, although its effect can be shadowed by the earthworm's diet. To untangle dietary from taxon effects, we raised nine earthworm species on a uniform diet of cow manure and compared cast microbiome across species while controlling for diet. Our results showed that, under controlled laboratory conditions, earthworm microbiomes are species-specific, more diverse than that of the controlled diet, and mainly comprised of native bacteria (i.e. not acquired from the diet). Furthermore, diet has a medium to large convergence effect on microbiome composition since earthworms shared 16%-74% of their bacterial amplicon sequence variants (ASV). The interspecies core microbiome included 10 ASVs, while their intraspecies core microbiomes were larger and varied in ASV richness (24%-48%) and sequence abundance across earthworm species. This specificity in core microbiomes and variable degree of similarity in bacterial composition suggest that phylosymbiosis could determine earthworm microbiome assembly. However, lack of congruence between the earthworm phylogeny and the microbiome dendrogram suggests that a consistent diet fed over several generations may have weakened potential phylosymbiotic effects. Thus, cast microbiome assembly in earthworms seem to be the result of an interplay among host phylogeny and diet., (© The Author(s) 2022. Published by Oxford University Press on behalf of FEMS.)

Elimination of human immunodeficiency virus (HIV) reservoirs is a critical endpoint to eradicate HIV. One therapeutic intervention against latent HIV is "shock and kill." This strategy is based on the transcriptional activation of latent HIV with a latency-reversing agent (LRA) with the consequent killing of the reactivated cell by either the cytopathic effect of HIV or the immune system. We have previously found that the small molecule 3-hydroxy-1,2,3-benzotriazin-4(3H)-one (HODHBt) acts as an LRA by increasing signal transducer and activator of transcription (STAT) factor activation mediated by interleukin-15 (IL-15) in cells isolated from aviremic participants. The IL-15 superagonist N-803 is currently under clinical investigation to eliminate latent reservoirs. IL-15 and N-803 share similar mechanisms of action by promoting the activation of STATs and have shown some promise in preclinical models directed toward HIV eradication. In this work, we evaluated the ability of HODHBt to enhance IL-15 signaling in natural killer (NK) cells and the biological consequences associated with increased STAT activation in NK cell effector and memory-like functions. We showed that HODHBt increased IL-15-mediated STAT phosphorylation in NK cells, resulting in increases in the secretion of CXCL-10 and interferon gamma (IFN-γ) and the expression of cytotoxic proteins, including granzyme B, granzyme A, perforin, granulysin, FASL, and TRAIL. This increased cytotoxic profile results in increased cytotoxicity against HIV-infected cells and different tumor cell lines. HODHBt also improved the generation of cytokine-induced memory-like NK cells. Overall, our data demonstrate that enhancing the magnitude of IL-15 signaling with HODHBt favors NK cell cytotoxicity and memory-like generation, and thus, targeting this pathway could be further explored for HIV cure interventions. IMPORTANCE Several clinical trials targeting the HIV latent reservoir with LRAs have been completed. In spite of a lack of clinical benefit, they have been crucial to elucidate hurdles that "shock and kill" strategies have to overcome to promote an effective reduction of the latent reservoir to lead to a cure. These hurdles include low reactivation potential mediated by LRAs, the negative influence of some LRAs on the activity of natural killer and effector CD8 T cells, an increased resistance to apoptosis of latently infected cells, and an exhausted immune system due to chronic inflammation. To that end, finding therapeutic strategies that can overcome some of these challenges could improve the outcome of shock and kill strategies aimed at HIV eradication. Here, we show that the LRA HODHBt also improves IL-15-mediated NK cell effector and memory-like functions. As such, pharmacological enhancement of IL-15-mediated STAT activation can open new therapeutic avenues toward an HIV cure.

Proteins are direct products of the genome and metabolites are functional products of interactions between the host and other factors such as environment, disease state, clinical information, etc. Omics data, including proteins and metabolites, are useful in characterizing biological processes underlying COVID-19 along with patient data and clinical information, yet few methods are available to effectively analyze such diverse and unstructured data. Using an integrated approach that combines proteomics and metabolomics data, we investigated the changes in metabolites and proteins in relation to patient characteristics (e.g., age, gender, and health outcome) and clinical information (e.g., metabolic panel and complete blood count test results). We found significant enrichment of biological indicators of lung, liver, and gastrointestinal dysfunction associated with disease severity using publicly available metabolite and protein profiles. Our analyses specifically identified enriched proteins that play a critical role in responses to injury or infection within these anatomical sites, but may contribute to excessive systemic inflammation within the context of COVID-19. Furthermore, we have used this information in conjunction with machine learning algorithms to predict the health status of patients presenting symptoms of COVID-19. This work provides a roadmap for understanding the biochemical pathways and molecular mechanisms that drive disease severity, progression, and treatment of COVID-19., (© 2022. The Author(s).)

Background: Mastitis pathogenesis involves a wide range of opportunistic and apparently resident microorganims including bacteria, viruses and archaea. In dairy animals, microbes reside in the host, interact with environment and evade the host immune system, providing a potential for host-tropism to favor mastitis pathogenesis. To understand the host-tropism phenomena of bovine-tropic mastitis microbiomes, we developed a cow-to-mouse mastitis model., Methods: A cow-to-mouse mastitis model was established by fecal microbiota transplantation (FMT) and milk microbiota transplantation (MMT) to pregnant mice to assess microbiome dysbiosis and genomic functional perturbations through shotgun whole metagenome sequencing (WMS) along with histopathological changes in mice mammary gland and colon tissues., Results: The cow-to-mouse FMT and MMT from clinical mastitis (CM) cows induced mastitis syndromes in mice as evidenced by histopathological changes in mammary gland and colon tissues. The WMS of 24 samples including six milk (CM = 3, healthy; H = 3), six fecal (CM = 4, H = 2) samples from cows, and six fecal (CM = 4, H = 2) and six mammary tissue (CM = 3, H = 3) samples from mice generating 517.14 million reads (average: 21.55 million reads/sample) mapped to 2191 bacterial, 94 viral and 54 archaeal genomes. The Kruskal-Wallis test revealed significant differences (p = 0.009) in diversity, composition, and relative abundances in microbiomes between CM- and H-metagenomes. These differences in microbiome composition were mostly represented by Pseudomonas aeruginosa, Lactobacillus crispatus, Klebsiella oxytoca, Enterococcus faecalis, Pantoea dispersa in CM-cows (feces and milk), and Muribaculum spp., Duncaniella spp., Muribaculum intestinale, Bifidobacterium animalis, Escherichia coli, Staphylococcus aureus, Massilia oculi, Ralstonia pickettii in CM-mice (feces and mammary tissues). Different species of Clostridia, Bacteroida, Actinobacteria, Flavobacteriia and Betaproteobacteria had a strong co-occurrence and positive correlation as the indicator species of murine mastitis. However, both CM cows and mice shared few mastitis-associated microbial taxa (1.14%) and functional pathways regardless of conservation of mastitis syndromes, indicating the higher discrepancy in mastitis-associated microbiomes among lactating mammals., Conclusions: We successfully induced mastitis by FMT and MMT that resulted in microbiome dysbiosis and genomic functional perturbations in mice. This study induced mastitis in a mouse model through FMT and MMT, which might be useful for further studies- focused on pathogen(s) involved in mastitis, their cross-talk among themselves and the host., (© 2022. The Author(s).)

CRF47&#95;BF is a circulating recombinant form (CRF) of the human immunodeficiency virus type 1 (HIV-1), the etiological agent of AIDS. CRF47&#95;BF represents one of 19 CRFx&#95;BFs and has a geographic focus in Spain, where it was first identified in 2010. Since its discovery, CRF47&#95;BF has expanded considerably in Spain, predominantly through heterosexual contact (∼56% of the infections). Little is known, however, about the origin and diversity of this CRF or its epidemiological correlates, as very few samples have been available so far. This study conducts a phylogenetic analysis with representatives of all CRFx&#95;BF sequence types along with HIV-1 M Group subtypes to validate that the CRF47&#95;BF sequences share a unique evolutionary history. The CRFx&#95;BF sequences cluster into a single, not well supported, clade that includes their dominant parent subtypes (B and F). This clade also includes subtype D and excludes sub-subtype F2. However, the CRF47&#95;BF sequences all share a most recent common ancestor. Further analysis of this clade couples CRF47&#95;BF protease-reverse transcriptase sequences and epidemiological data from an additional 87 samples collected throughout Spain, as well as additional CRF47&#95;BF database sequences from Brazil and Spain to investigate the origin and phylodynamics of CRF47&#95;BF. The Spanish region with the highest proportion of CRF47&#95;BF samples in the data set was the Basque Country (43.7%) with Navarre next highest at 19.5%. We include in our analysis epidemiological data on host sex, mode of transmission, time of collection, and geographic region. The phylodynamic analysis indicates that CRF47&#95;BF originated in Brazil around 1999-2000 and spread to Spain from Brazil in 2002-2003. The virus spread rapidly throughout Spain with an increase in population size from 2011 to 2015 and leveling off more recently. Three strongly supported clusters associated with Spanish regions (Basque Country, Navarre, and Aragon), together comprising 60.8% of the Spanish samples, were identified, one of which was also associated with transmission among men who have sex with men. The expansion in Spain of CRF47&#95;BF, together with that of other CRFs and subtype variants of South American origin, previously reported, reflects the increasing relationship between the South American and European HIV-1 epidemics., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Hill, Pérez-Losada, Delgado, Benito, Montero, Gil, Sánchez, Cañada-García, García-Bodas, Crandall, Thomson and the Spanish Group for the Study of New HIV Diagnoses.)

Earthworms heavily modify the soil microbiome as it passes throughout their guts. However, there are no detailed studies describing changes in the composition, structure and diversity of soil microbiomes during gut transit and once they are released back to the soil as casts. To address this knowledge gap, we used 16S rRNA next-generation sequencing to characterize the microbiomes of soil, gut and casts from the earthworm Aporrectodea caliginosa . We also studied whether these three microbiomes are clearly distinct in composition or can be merged into metacommunities. A large proportion of bacteria was unique to each microbiome-soil (82%), gut (89%) and casts (75%), which indicates that the soil microbiome is greatly modified during gut transit. The three microbiomes also differed in alpha diversity, which peaked during gut transit and decreased in casts. Furthermore, gut transit also modified the structure of the soil microbiome, which clustered away from those of the earthworm gut and cast samples. However, this clustering pattern was not supported by metacommunity analysis, which indicated that soil and gut samples make up one metacommunity and cast samples another. These results have important implications for understanding the dynamics of soil microbial communities and nutrient cycles.

The study of microbial communities or microbiotas in animals and environments is important because of their impact in a broad range of industrial applications, diseases and ecological roles. High throughput sequencing (HTS) is the best strategy to characterize microbial composition and function. Microbial profiles can be obtained either by shotgun sequencing of genomes, or through amplicon sequencing of target genes (e.g., 16S rRNA for bacteria and ITS for fungi). Here, we compared both HTS approaches at assessing taxonomic and functional diversity of bacterial and fungal communities during vermicomposting of white grape marc. We applied specific HTS workflows to the same 12 microcosms, with and without earthworms, sampled at two distinct phases of the vermicomposting process occurring at 21 and 63 days. Metataxonomic profiles were inferred in DADA2, with bacterial metabolic pathways predicted via PICRUSt2. Metagenomic taxonomic profiles were inferred in PathoScope, while bacterial functional profiles were inferred in Humann2. Microbial profiles inferred by metagenomics and metataxonomics showed similarities and differences in composition, structure, and metabolic function at different taxonomic levels. Microbial composition and abundance estimated by both HTS approaches agreed reasonably well at the phylum level, but larger discrepancies were observed at lower taxonomic ranks. Shotgun HTS identified ~1.8 times more bacterial genera than 16S rRNA HTS, while ITS HTS identified two times more fungal genera than shotgun HTS. This is mainly a consequence of the difference in resolution and reference richness between amplicon and genome sequencing approaches and databases, respectively. Our study also revealed great differences and even opposite trends in alpha- and beta-diversity between amplicon and shotgun HTS. Interestingly, amplicon PICRUSt2-imputed functional repertoires overlapped ~50% with shotgun Humann2 profiles. Finally, both approaches indicated that although bacteria and fungi are the main drivers of biochemical decomposition, earthworms also play a key role in plant vermicomposting. In summary, our study highlights the strengths and weaknesses of metagenomics and metataxonomics and provides new insights on the vermicomposting of white grape marc. Since both approaches may target different biological aspects of the communities, combining them will provide a better understanding of the microbiotas under study., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Pérez-Losada, Narayanan, Kolbe, Ramos-Tapia, Castro-Nallar, Crandall and Domínguez.)

Introgression is an important biological process affecting at least 10% of the extant species in the animal kingdom. Introgression significantly impacts inference of phylogenetic species relationships where a strictly binary tree model cannot adequately explain reticulate net-like species relationships. Here, we use phylogenomic approaches to understand patterns of introgression along the evolutionary history of a unique, nonmodel insect system: dragonflies and damselflies (Odonata). We demonstrate that introgression is a pervasive evolutionary force across various taxonomic levels within Odonata. In particular, we show that the morphologically "intermediate" species of Anisozygoptera (one of the three primary suborders within Odonata besides Zygoptera and Anisoptera), which retain phenotypic characteristics of the other two suborders, experienced high levels of introgression likely coming from zygopteran genomes. Additionally, we find evidence for multiple cases of deep inter-superfamilial ancestral introgression. [Gene flow; Odonata; phylogenomics; reticulate evolution.]., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society of Systematic Biologists.)

Objective: In response to COVID-19, the informatics community united to aggregate as much clinical data as possible to characterize this new disease and reduce its impact through collaborative analytics. The National COVID Cohort Collaborative (N3C) is now the largest publicly available HIPAA limited dataset in US history with over 6.4 million patients and is a testament to a partnership of over 100 organizations., Materials and Methods: We developed a pipeline for ingesting, harmonizing, and centralizing data from 56 contributing data partners using 4 federated Common Data Models. N3C data quality (DQ) review involves both automated and manual procedures. In the process, several DQ heuristics were discovered in our centralized context, both within the pipeline and during downstream project-based analysis. Feedback to the sites led to many local and centralized DQ improvements., Results: Beyond well-recognized DQ findings, we discovered 15 heuristics relating to source Common Data Model conformance, demographics, COVID tests, conditions, encounters, measurements, observations, coding completeness, and fitness for use. Of 56 sites, 37 sites (66%) demonstrated issues through these heuristics. These 37 sites demonstrated improvement after receiving feedback., Discussion: We encountered site-to-site differences in DQ which would have been challenging to discover using federated checks alone. We have demonstrated that centralized DQ benchmarking reveals unique opportunities for DQ improvement that will support improved research analytics locally and in aggregate., Conclusion: By combining rapid, continual assessment of DQ with a large volume of multisite data, it is possible to support more nuanced scientific questions with the scale and rigor that they require., (© The Author(s) 2021. Published by Oxford University Press on behalf of the American Medical Informatics Association.)

Background: The barnacles are a group of >2,000 species that have fascinated biologists, including Darwin, for centuries. Their lifestyles are extremely diverse, from free-swimming larvae to sessile adults, and even root-like endoparasites. Barnacles also cause hundreds of millions of dollars of losses annually due to biofouling. However, genomic resources for crustaceans, and barnacles in particular, are lacking., Results: Using 62× Pacific Biosciences coverage, 189× Illumina whole-genome sequencing coverage, 203× HiC coverage, and 69× CHi-C coverage, we produced a chromosome-level genome assembly of the gooseneck barnacle Pollicipes pollicipes. The P. pollicipes genome is 770 Mb long and its assembly is one of the most contiguous and complete crustacean genomes available, with a scaffold N50 of 47 Mb and 90.5% of the BUSCO Arthropoda gene set. Using the genome annotation produced here along with transcriptomes of 13 other barnacle species, we completed phylogenomic analyses on a nearly 2 million amino acid alignment. Contrary to previous studies, our phylogenies suggest that the Pollicipedomorpha is monophyletic and sister to the Balanomorpha, which alters our understanding of barnacle larval evolution and suggests homoplasy in a number of naupliar characters. We also compared transcriptomes of P. pollicipes nauplius larvae and adults and found that nearly one-half of the genes in the genome are differentially expressed, highlighting the vastly different transcriptomes of larvae and adult gooseneck barnacles. Annotation of the genes with KEGG and GO terms reveals that these stages exhibit many differences including cuticle binding, chitin binding, microtubule motor activity, and membrane adhesion., Conclusion: This study provides high-quality genomic resources for a key group of crustaceans. This is especially valuable given the roles P. pollicipes plays in European fisheries, as a sentinel species for coastal ecosystems, and as a model for studying barnacle adhesion as well as its key position in the barnacle tree of life. A combination of genomic, phylogenetic, and transcriptomic analyses here provides valuable insights into the evolution and development of barnacles., (© The Author(s) 2022. Published by Oxford University Press GigaScience.)

Importance: Prior observational studies suggest that aspirin use may be associated with reduced mortality in high-risk hospitalized patients with COVID-19, but aspirin's efficacy in patients with moderate COVID-19 is not well studied., Objective: To assess whether early aspirin use is associated with lower odds of in-hospital mortality in patients with moderate COVID-19., Design, Setting, and Participants: Observational cohort study of 112 269 hospitalized patients with moderate COVID-19, enrolled from January 1, 2020, through September 10, 2021, at 64 health systems in the United States participating in the National Institute of Health's National COVID Cohort Collaborative (N3C)., Exposure: Aspirin use within the first day of hospitalization., Main Outcome and Measures: The primary outcome was 28-day in-hospital mortality, and secondary outcomes were pulmonary embolism and deep vein thrombosis. Odds of in-hospital mortality were calculated using marginal structural Cox and logistic regression models. Inverse probability of treatment weighting was used to reduce bias from confounding and balance characteristics between groups., Results: Among the 2 446 650 COVID-19-positive patients who were screened, 189 287 were hospitalized and 112 269 met study inclusion. For the full cohort, Median age was 63 years (IQR, 47-74 years); 16.1% of patients were African American, 3.8% were Asian, 52.7% were White, 5.0% were of other races and ethnicities, 22.4% were of unknown race and ethnicity. In-hospital mortality occurred in 10.9% of patients. After inverse probability treatment weighting, 28-day in-hospital mortality was significantly lower in those who received aspirin (10.2% vs 11.8%; odds ratio [OR], 0.85; 95% CI, 0.79-0.92; P < .001). The rate of pulmonary embolism, but not deep vein thrombosis, was also significantly lower in patients who received aspirin (1.0% vs 1.4%; OR, 0.71; 95% CI, 0.56-0.90; P = .004). Patients who received early aspirin did not have higher rates of gastrointestinal hemorrhage (0.8% aspirin vs 0.7% no aspirin; OR, 1.04; 95% CI, 0.82-1.33; P = .72), cerebral hemorrhage (0.6% aspirin vs 0.4% no aspirin; OR, 1.32; 95% CI, 0.92-1.88; P = .13), or blood transfusion (2.7% aspirin vs 2.3% no aspirin; OR, 1.14; 95% CI, 0.99-1.32; P = .06). The composite of hemorrhagic complications did not occur more often in those receiving aspirin (3.7% aspirin vs 3.2% no aspirin; OR, 1.13; 95% CI, 1.00-1.28; P = .054). Subgroups who appeared to benefit the most included patients older than 60 years (61-80 years: OR, 0.79; 95% CI, 0.72-0.87; P < .001; >80 years: OR, 0.79; 95% CI, 0.69-0.91; P < .001) and patients with comorbidities (1 comorbidity: 6.4% vs 9.2%; OR, 0.68; 95% CI, 0.55-0.83; P < .001; 2 comorbidities: 10.5% vs 12.8%; OR, 0.80; 95% CI, 0.69-0.93; P = .003; 3 comorbidities: 13.8% vs 17.0%, OR, 0.78; 95% CI, 0.68-0.89; P < .001; >3 comorbidities: 17.0% vs 21.6%; OR, 0.74; 95% CI, 0.66-0.84; P < .001)., Conclusions and Relevance: In this cohort study of US adults hospitalized with moderate COVID-19, early aspirin use was associated with lower odds of 28-day in-hospital mortality. A randomized clinical trial that includes diverse patients with moderate COVID-19 is warranted to adequately evaluate aspirin's efficacy in patients with high-risk conditions.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has evolved into eight fundamental clades with four of these clades (G, GH, GR, and GV) globally prevalent in 2020. To explain plausible epistatic effects of the signature co-occurring mutations of these circulating clades on viral replication and transmission fitness, we proposed a hypothetical model using in silico approach. Molecular docking and dynamics analyses showed the higher infectiousness of a spike mutant through more favorable binding of G 614 with the elastase-2. RdRp mutation p.P323L significantly increased genome-wide mutations (p < 0.0001), allowing for more flexible RdRp (mutated)-NSP8 interaction that may accelerate replication. Superior RNA stability and structural variation at NSP3:C241T might impact protein, RNA interactions, or both. Another silent 5'-UTR:C241T mutation might affect translational efficiency and viral packaging. These four G-clade-featured co-occurring mutations might increase viral replication. Sentinel GH-clade ORF3a:p.Q57H variants constricted the ion-channel through intertransmembrane-domain interaction of cysteine(C81)-histidine(H57). The GR-clade N:p.RG203-204KR would stabilize RNA interaction by a more flexible and hypo-phosphorylated SR-rich region. GV-clade viruses seemingly gained the evolutionary advantage of the confounding factors; nevertheless, N:p.A220V might modulate RNA binding with no phenotypic effect. Our hypothetical model needs further retrospective and prospective studies to understand detailed molecular events and their relationship to the fitness of SARS-CoV-2., (© 2021 Wiley Periodicals LLC.)

Life on Earth has evolved from initial simplicity to the astounding complexity we experience today. Bacteria and archaea have largely excelled in metabolic diversification, but eukaryotes additionally display abundant morphological innovation. How have these innovations come about and what constraints are there on the origins of novelty and the continuing maintenance of biodiversity on Earth? The history of life and the code for the working parts of cells and systems are written in the genome. The Earth BioGenome Project has proposed that the genomes of all extant, named eukaryotes-about 2 million species-should be sequenced to high quality to produce a digital library of life on Earth, beginning with strategic phylogenetic, ecological, and high-impact priorities. Here we discuss why we should sequence all eukaryotic species, not just a representative few scattered across the many branches of the tree of life. We suggest that many questions of evolutionary and ecological significance will only be addressable when whole-genome data representing divergences at all of the branchings in the tree of life or all species in natural ecosystems are available. We envisage that a genomic tree of life will foster understanding of the ongoing processes of speciation, adaptation, and organismal dependencies within entire ecosystems. These explorations will resolve long-standing problems in phylogenetics, evolution, ecology, conservation, agriculture, bioindustry, and medicine., Competing Interests: The authors declare no competing interest., (Copyright © 2022 the Author(s). Published by PNAS.)

A global international initiative, such as the Earth BioGenome Project (EBP), requires both agreement and coordination on standards to ensure that the collective effort generates rapid progress toward its goals. To this end, the EBP initiated five technical standards committees comprising volunteer members from the global genomics scientific community: Sample Collection and Processing, Sequencing and Assembly, Annotation, Analysis, and IT and Informatics. The current versions of the resulting standards documents are available on the EBP website, with the recognition that opportunities, technologies, and challenges may improve or change in the future, requiring flexibility for the EBP to meet its goals. Here, we describe some highlights from the proposed standards, and areas where additional challenges will need to be met., Competing Interests: Competing interest statement: P.F. is a member of the scientific advisory boards of Fabric Genomics, Inc. and Eagle Genomics, Ltd., (Copyright © 2022 the Author(s). Published by PNAS.)

Competing Interests: The authors declare no competing interest.

The microbiota of the nasopharyngeal tract (NT) play a role in host immunity against respiratory infectious diseases. However, scant information is available on interactions of SARS-CoV-2 with the nasopharyngeal microbiome. This study characterizes the effects of SARS-CoV-2 infection on human nasopharyngeal microbiomes and their relevant metabolic functions. Twenty-two (n = 22) nasopharyngeal swab samples (including COVID-19 patients = 8, recovered humans = 7, and healthy people = 7) were collected, and underwent to RNAseq-based metagenomic investigation. Our RNAseq data mapped to 2281 bacterial species (including 1477, 919 and 676 in healthy, COVID-19 and recovered metagenomes, respectively) indicating a distinct microbiome dysbiosis. The COVID-19 and recovered samples included 67% and 77% opportunistic bacterial species, respectively compared to healthy controls. Notably, 79% commensal bacterial species found in healthy controls were not detected in COVID-19 and recovered people. Similar dysbiosis was also found in viral and archaeal fraction of the nasopharyngeal microbiomes. We also detected several altered metabolic pathways and functional genes in the progression and pathophysiology of COVID-19. The nasopharyngeal microbiome dysbiosis and their genomic features determined by our RNAseq analyses shed light on early interactions of SARS-CoV-2 with the nasopharyngeal resident microbiota that might be helpful for developing microbiome-based diagnostics and therapeutics for this novel pandemic disease., (© 2021. The Author(s).)

SARS-CoV-2 (CoV) is the etiological agent of the COVID-19 pandemic and evolves to evade both host immune systems and intervention strategies. We divided the CoV genome into 29 constituent regions and applied novel analytical approaches to identify associations between CoV genomic features and epidemiological metadata. Our results show that nonstructural protein 3 (nsp3) and Spike protein (S) have the highest variation and greatest correlation with the viral whole-genome variation. S protein variation is correlated with nsp3, nsp6, and 3'-to-5' exonuclease variation. Country of origin and time since the start of the pandemic were the most influential metadata associated with genomic variation, while host sex and age were the least influential. We define a novel statistic-coherence-and show its utility in identifying geographic regions (populations) with unusually high (many new variants) or low (isolated) viral phylogenetic diversity. Interestingly, at both global and regional scales, we identify geographic locations with high coherence neighboring regions of low coherence; this emphasizes the utility of this metric to inform public health measures for disease spread. Our results provide a direction to prioritize genes associated with outcome predictors (e.g., health, therapeutic, and vaccine outcomes) and to improve DNA tests for predicting disease status., (© 2021. The Author(s).)

In response to the COVID-19 outbreak, scientists and medical researchers are capturing a wide range of host responses, symptoms and lingering postrecovery problems within the human population. These variable clinical manifestations suggest differences in influential factors, such as innate and adaptive host immunity, existing or underlying health conditions, comorbidities, genetics and other factors-compounding the complexity of COVID-19 pathobiology and potential biomarkers associated with the disease, as they become available. The heterogeneous data pose challenges for efficient extrapolation of information into clinical applications. We have curated 145 COVID-19 biomarkers by developing a novel cross-cutting disease biomarker data model that allows integration and evaluation of biomarkers in patients with comorbidities. Most biomarkers are related to the immune (SAA, TNF-∝ and IP-10) or coagulation (D-dimer, antithrombin and VWF) cascades, suggesting complex vascular pathobiology of the disease. Furthermore, we observe commonality with established cancer biomarkers (ACE2, IL-6, IL-4 and IL-2) as well as biomarkers for metabolic syndrome and diabetes (CRP, NLR and LDL). We explore these trends as we put forth a COVID-19 biomarker resource (https://data.oncomx.org/covid19) that will help researchers and diagnosticians alike., (© The Author(s) 2021. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)