Rega Institute for Medical Research

Bayesian inference of phylogeny using Markov chain Monte Carlo (MCMC) plays a central role in understanding evolutionary history from molecular sequence data. Visualizing and analyzing the MCMC-generated samples from the posterior distribution is a key step in any non-trivial Bayesian inference. We present the software package Tracer (version 1.7) for visualizing and analyzing the MCMC trace files generated through Bayesian phylogenetic inference. Tracer provides kernel density estimation, multivariate visualization, demographic trajectory reconstruction, conditional posterior distribution summary, and more. Tracer is open-source and available at http://beast.community/tracer.

The Bayesian Evolutionary Analysis by Sampling Trees (BEAST) software package has become a primary tool for Bayesian phylogenetic and phylodynamic inference from genetic sequence data. BEAST unifies molecular phylogenetic reconstruction with complex discrete and continuous trait evolution, divergence-time dating, and coalescent demographic models in an efficient statistical inference engine using Markov chain Monte Carlo integration. A convenient, cross-platform, graphical user interface allows the flexible construction of complex evolutionary analyses.

Microbes are dominant drivers of biogeochemical processes, yet drawing a global picture of functional diversity, microbial community structure, and their ecological determinants remains a grand challenge. We analyzed 7.2 terabases of metagenomic data from 243 Tara Oceans samples from 68 locations in epipelagic and mesopelagic waters across the globe to generate an ocean microbial reference gene catalog with >40 million nonredundant, mostly novel sequences from viruses, prokaryotes, and picoeukaryotes. Using 139 prokaryote-enriched samples, containing >35,000 species, we show vertical stratification with epipelagic community composition mostly driven by temperature rather than other environmental factors or geography. We identify ocean microbial core functionality and reveal that >73% of its abundance is shared with the human gut microbiome despite the physicochemical differences between these two ecosystems.

Marine plankton support global biological and geochemical processes. Surveys of their biodiversity have hitherto been geographically restricted and have not accounted for the full range of plankton size. We assessed eukaryotic diversity from 334 size-fractionated photic-zone plankton communities collected across tropical and temperate oceans during the circumglobal Tara Oceans expedition. We analyzed 18S ribosomal DNA sequences across the intermediate plankton-size spectrum from the smallest unicellular eukaryotes (protists, >0.8 micrometers) to small animals of a few millimeters. Eukaryotic ribosomal diversity saturated at ~150,000 operational taxonomic units, about one-third of which could not be assigned to known eukaryotic groups. Diversity emerged at all taxonomic levels, both within the groups comprising the ~11,200 cataloged morphospecies of eukaryotic plankton and among twice as many other deep-branching lineages of unappreciated importance in plankton ecology studies. Most eukaryotic plankton biodiversity belonged to heterotrophic protistan groups, particularly those known to be parasites or symbiotic hosts.

Deep sequencing of the gut microbiomes of 1135 participants from a Dutch population-based cohort shows relations between the microbiome and 126 exogenous and intrinsic host factors, including 31 intrinsic factors, 12 diseases, 19 drug groups, 4 smoking categories, and 60 dietary factors. These factors collectively explain 18.7% of the variation seen in the interindividual distance of microbial composition. We could associate 110 factors to 125 species and observed that fecal chromogranin A (CgA), a protein secreted by enteroendocrine cells, was exclusively associated with 61 microbial species whose abundance collectively accounted for 53% of microbial composition. Low CgA concentrations were seen in individuals with a more diverse microbiome. These results are an important step toward a better understanding of environment-diet-microbe-host interactions.

Short chain fatty acids (SCFAs), including acetate, propionate, and butyrate, are produced at high concentration by bacteria in the gut and subsequently released in the bloodstream. Basal acetate concentrations in the blood (about 100 microm) can further increase to millimolar concentrations following alcohol intake. It was known previously that SCFAs can activate leukocytes, particularly neutrophils. In the present work, we have identified two previously orphan G protein-coupled receptors, GPR41 and GPR43, as receptors for SCFAs. Propionate was the most potent agonist for both GPR41 and GPR43. Acetate was more selective for GPR43, whereas butyrate and isobutyrate were more active on GPR41. The two receptors were coupled to inositol 1,4,5-trisphosphate formation, intracellular Ca2+ release, ERK1/2 activation, and inhibition of cAMP accumulation. They exhibited, however, a differential coupling to G proteins; GPR41 coupled exclusively though the Pertussis toxin-sensitive Gi/o family, whereas GPR43 displayed a dual coupling through Gi/o and Pertussis toxin-insensitive Gq protein families. The broad expression profile of GPR41 in a number of tissues does not allow us to infer clear hypotheses regarding its biological functions. In contrast, the highly selective expression of GPR43 in leukocytes, particularly polymorphonuclear cells, suggests a role in the recruitment of these cell populations toward sites of bacterial infection. The pharmacology of GPR43 matches indeed the effects of SCFAs on neutrophils, in terms of intracellular Ca2+ release and chemotaxis. Such a neutrophil-specific SCFA receptor is potentially involved in the development of a variety of diseases characterized by either excessive or inefficient neutrophil recruitment and activation, such as inflammatory bowel diseases or alcoholism-associated immune depression. GPR43 might therefore constitute a target allowing us to modulate immune responses in these pathological situations.

Cases of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection in Manaus, Brazil, resurged in late 2020 despite previously high levels of infection. Genome sequencing of viruses sampled in Manaus between November 2020 and January 2021 revealed the emergence and circulation of a novel SARS-CoV-2 variant of concern. Lineage P.1 acquired 17 mutations, including a trio in the spike protein (K417T, E484K, and N501Y) associated with increased binding to the human ACE2 (angiotensin-converting enzyme 2) receptor. Molecular clock analysis shows that P.1 emergence occurred around mid-November 2020 and was preceded by a period of faster molecular evolution. Using a two-category dynamical model that integrates genomic and mortality data, we estimate that P.1 may be 1.7- to 2.4-fold more transmissible and that previous (non-P.1) infection provides 54 to 79% of the protection against infection with P.1 that it provides against non-P.1 lineages. Enhanced global genomic surveillance of variants of concern, which may exhibit increased transmissibility and/or immune evasion, is critical to accelerate pandemic responsiveness.

Since the first antiviral drug, idoxuridine, was approved in 1963, 90 antiviral drugs categorized into 13 functional groups have been formally approved for the treatment of the following 9 human infectious diseases: (i) HIV infections (protease inhibitors, integrase inhibitors, entry inhibitors, nucleoside reverse transcriptase inhibitors, nonnucleoside reverse transcriptase inhibitors, and acyclic nucleoside phosphonate analogues), (ii) hepatitis B virus (HBV) infections (lamivudine, interferons, nucleoside analogues, and acyclic nucleoside phosphonate analogues), (iii) hepatitis C virus (HCV) infections (ribavirin, interferons, NS3/4A protease inhibitors, NS5A inhibitors, and NS5B polymerase inhibitors), (iv) herpesvirus infections (5-substituted 2'-deoxyuridine analogues, entry inhibitors, nucleoside analogues, pyrophosphate analogues, and acyclic guanosine analogues), (v) influenza virus infections (ribavirin, matrix 2 protein inhibitors, RNA polymerase inhibitors, and neuraminidase inhibitors), (vi) human cytomegalovirus infections (acyclic guanosine analogues, acyclic nucleoside phosphonate analogues, pyrophosphate analogues, and oligonucleotides), (vii) varicella-zoster virus infections (acyclic guanosine analogues, nucleoside analogues, 5-substituted 2'-deoxyuridine analogues, and antibodies), (viii) respiratory syncytial virus infections (ribavirin and antibodies), and (ix) external anogenital warts caused by human papillomavirus infections (imiquimod, sinecatechins, and podofilox). Here, we present for the first time a comprehensive overview of antiviral drugs approved over the past 50 years, shedding light on the development of effective antiviral treatments against current and emerging infectious diseases worldwide.

The microbiota of the human gut is gaining broad attention owing to its association with a wide range of diseases, ranging from metabolic disorders (e.g. obesity and type 2 diabetes) to autoimmune diseases (such as inflammatory bowel disease and type 1 diabetes), cancer and even neurodevelopmental disorders (e.g. autism). Having been increasingly used in biomedical research, mice have become the model of choice for most studies in this emerging field. Mouse models allow perturbations in gut microbiota to be studied in a controlled experimental setup, and thus help in assessing causality of the complex host-microbiota interactions and in developing mechanistic hypotheses. However, pitfalls should be considered when translating gut microbiome research results from mouse models to humans. In this Special Article, we discuss the intrinsic similarities and differences that exist between the two systems, and compare the human and murine core gut microbiota based on a meta-analysis of currently available datasets. Finally, we discuss the external factors that influence the capability of mouse models to recapitulate the gut microbiota shifts associated with human diseases, and investigate which alternative model systems exist for gut microbiota research.

In this study, we demonstrate that the CXC family of chemokines displays disparate angiogenic activity depending upon the presence or absence of the ELR motif. CXC chemokines containing the ELR motif (ELR-CXC chemokines) were found to be potent angiogenic factors, inducing both in vitro endothelial chemotaxis and in vivo corneal neovascularization. In contrast, the CXC chemokines lacking the ELR motif, platelet factor 4, interferon γ-inducible protein 10, and monokine induced by γ-interferon, not only failed to induce significant in vitro endothelial cell chemotaxis or in vivo corneal neovacularization but were found to be potent angiostatic factors in the presence of either ELR-CXC chemokines or the unrelated angiogenic factor, basic fibroblast growth factor. Additionally, mutant interleukin-8 proteins lacking the ELR motif demonstrated potent angiostatic effects in the presence of either ELR-CXC chemokines or basic fibroblast growth factor. In contrast, a mutant of monokine induced by γ-interferon containing the ELR motif was found to induce in vivo angiogenic activity. These findings suggest a functional role of the ELR motif in determining the angiogenic or angiostatic potential of CXC chemokines, supporting the hypothesis that the net biological balance between angiogenic and angiostatic CXC chemokines may play an important role in regulating overall angiogenesis. In this study, we demonstrate that the CXC family of chemokines displays disparate angiogenic activity depending upon the presence or absence of the ELR motif. CXC chemokines containing the ELR motif (ELR-CXC chemokines) were found to be potent angiogenic factors, inducing both in vitro endothelial chemotaxis and in vivo corneal neovascularization. In contrast, the CXC chemokines lacking the ELR motif, platelet factor 4, interferon γ-inducible protein 10, and monokine induced by γ-interferon, not only failed to induce significant in vitro endothelial cell chemotaxis or in vivo corneal neovacularization but were found to be potent angiostatic factors in the presence of either ELR-CXC chemokines or the unrelated angiogenic factor, basic fibroblast growth factor. Additionally, mutant interleukin-8 proteins lacking the ELR motif demonstrated potent angiostatic effects in the presence of either ELR-CXC chemokines or basic fibroblast growth factor. In contrast, a mutant of monokine induced by γ-interferon containing the ELR motif was found to induce in vivo angiogenic activity. These findings suggest a functional role of the ELR motif in determining the angiogenic or angiostatic potential of CXC chemokines, supporting the hypothesis that the net biological balance between angiogenic and angiostatic CXC chemokines may play an important role in regulating overall angiogenesis.

The potential for using Adeno-associated virus (AAV) as a vector for human gene therapy has stimulated interest in the Dependovirus genus. Serologic data suggest that AAV infections are prevalent in humans, although analyses of viruses and viral sequences from clinical samples are extremely limited. Molecular techniques were used in this study to successfully detect endogenous AAV sequences in 18% of all human tissues screened, with the liver and bone marrow being the most predominant sites. Sequence characterization of rescued AAV DNAs indicated a diverse array of molecular forms which segregate into clades whose members share functional and serologic similarities. One of the most predominant human clades is a hybrid of two previously described AAV serotypes, while another clade was found in humans and several species of nonhuman primates, suggesting a cross-species transmission of this virus. These data provide important information regarding the biology of parvoviruses in humans and their use as gene therapy vectors.

OBJECTIVE: The assessment of potentially confounding factors affecting colon microbiota composition is essential to the identification of robust microbiome based disease markers. Here, we investigate the link between gut microbiota variation and stool consistency using Bristol Stool Scale classification, which reflects faecal water content and activity, and is considered a proxy for intestinal colon transit time. DESIGN: Through 16S rDNA Illumina profiling of faecal samples of 53 healthy women, we evaluated associations between microbiome richness, Bacteroidetes:Firmicutes ratio, enterotypes, and genus abundance with self-reported, Bristol Stool Scale-based stool consistency. Each sample's microbiota growth potential was calculated to test whether transit time acts as a selective force on gut bacterial growth rates. RESULTS: Stool consistency strongly correlates with all known major microbiome markers. It is negatively correlated with species richness, positively associated to the Bacteroidetes:Firmicutes ratio, and linked to Akkermansia and Methanobrevibacter abundance. Enterotypes are distinctly distributed over the BSS-scores. Based on the correlations between microbiota growth potential and stool consistency scores within both enterotypes, we hypothesise that accelerated transit contributes to colon ecosystem differentiation. While shorter transit times can be linked to increased abundance of fast growing species in Ruminococcaceae-Bacteroides samples, hinting to a washout avoidance strategy of faster replication, this trend is absent in Prevotella-enterotyped individuals. Within this enterotype adherence to host tissue therefore appears to be a more likely bacterial strategy to cope with washout. CONCLUSIONS: The strength of the associations between stool consistency and species richness, enterotypes and community composition emphasises the crucial importance of stool consistency assessment in gut metagenome-wide association studies.

Species interaction networks are shaped by abiotic and biotic factors. Here, as part of the Tara Oceans project, we studied the photic zone interactome using environmental factors and organismal abundance profiles and found that environmental factors are incomplete predictors of community structure. We found associations across plankton functional types and phylogenetic groups to be nonrandomly distributed on the network and driven by both local and global patterns. We identified interactions among grazers, primary producers, viruses, and (mainly parasitic) symbionts and validated network-generated hypotheses using microscopy to confirm symbiotic relationships. We have thus provided a resource to support further research on ocean food webs and integrating biological components into ocean models.

Programs that monitor local, national, and regional levels of transmitted HIV-1 drug resistance inform treatment guidelines and provide feedback on the success of HIV-1 treatment and prevention programs. To accurately compare transmitted drug resistance rates across geographic regions and times, the World Health Organization has recommended the adoption of a consensus genotypic definition of transmitted HIV-1 drug resistance. In January 2007, we outlined criteria for developing a list of mutations for drug-resistance surveillance and compiled a list of 80 RT and protease mutations meeting these criteria (surveillance drug resistance mutations; SDRMs). Since January 2007, several new drugs have been approved and several new drug-resistance mutations have been identified. In this paper, we follow the same procedures described previously to develop an updated list of SDRMs that are likely to be useful for ongoing and future studies of transmitted drug resistance. The updated SDRM list has 93 mutations including 34 NRTI-resistance mutations at 15 RT positions, 19 NNRTI-resistance mutations at 10 RT positions, and 40 PI-resistance mutations at 18 protease positions.

The matrix metalloproteinases (MMPs) form an enzyme family of which gelatinase B (MMP-9) represents the largest and most complex member. We focus here on the biochemical properties, regulation, and functions of gelatinase B. The tight regulation of gelatinase B activity is highly complex and is established at five different levels. The transcription of the gelatinase B-gene depends on various cis-elements in its gene promotor and is induced or repressed by a large variety of soluble factors, including cytokines, growth factors, and hormones and by cellular contacts acting through specific signaling pathways. The specific regulation of its secretion occurs in cells storing gelatinase B in granules. After secretion, progelatinase B must be activated through an activation network. The enzyme activity is further regulated by inhibition and by other mechanisms, such as fine-tuning and stabilization by glycosylation. The ability of gelatinase B to degrade components of the extracellular matrix and to regulate the activity of a number of soluble proteins confers an important role in various physiological and pathological processes. These include reproduction, growth, development, inflammation, and vascular and proliferative diseases.

The liver is the largest solid organ in the body, yet it remains incompletely characterized. Here we present a spatial proteogenomic atlas of the healthy and obese human and murine liver combining single-cell CITE-seq, single-nuclei sequencing, spatial transcriptomics, and spatial proteomics. By integrating these multi-omic datasets, we provide validated strategies to reliably discriminate and localize all hepatic cells, including a population of lipid-associated macrophages (LAMs) at the bile ducts. We then align this atlas across seven species, revealing the conserved program of bona fide Kupffer cells and LAMs. We also uncover the respective spatially resolved cellular niches of these macrophages and the microenvironmental circuits driving their unique transcriptomic identities. We demonstrate that LAMs are induced by local lipid exposure, leading to their induction in steatotic regions of the murine and human liver, while Kupffer cell development crucially depends on their cross-talk with hepatic stellate cells via the evolutionarily conserved ALK1-BMP9/10 axis.

Group A rotavirus classification is currently based on the molecular properties of the two outer layer proteins, VP7 and VP4, and the middle layer protein, VP6. As reassortment of all the 11 rotavirus gene segments plays a key role in generating rotavirus diversity in nature, a classification system that is based on all the rotavirus gene segments is desirable for determining which genes influence rotavirus host range restriction, replication, and virulence, as well as for studying rotavirus epidemiology and evolution. Toward establishing such a classification system, gene sequences encoding VP1 to VP3, VP6, and NSP1 to NSP5 were determined for human and animal rotavirus strains belonging to different G and P genotypes in addition to those available in databases, and they were used to define phylogenetic relationships among all rotavirus genes. Based on these phylogenetic analyses, appropriate identity cutoff values were determined for each gene. For the VP4 gene, a nucleotide identity cutoff value of 80% completely correlated with the 27 established P genotypes. For the VP7 gene, a nucleotide identity cutoff value of 80% largely coincided with the established G genotypes but identified four additional distinct genotypes comprised of murine or avian rotavirus strains. Phylogenetic analyses of the VP1 to VP3, VP6, and NSP1 to NSP5 genes showed the existence of 4, 5, 6, 11, 14, 5, 7, 11, and 6 genotypes, respectively, based on nucleotide identity cutoff values of 83%, 84%, 81%, 85%, 79%, 85%, 85%, 85%, and 91%, respectively. In accordance with these data, a revised nomenclature of rotavirus strains is proposed. The novel classification system allows the identification of (i) distinct genotypes, which probably followed separate evolutionary paths; (ii) interspecies transmissions and a plethora of reassortment events; and (iii) certain gene constellations that revealed (a) a common origin between human Wa-like rotavirus strains and porcine rotavirus strains and (b) a common origin between human DS-1-like rotavirus strains and bovine rotaviruses. These close evolutionary links between human and animal rotaviruses emphasize the need for close simultaneous monitoring of rotaviruses in animals and humans.

Disruption of healthy microbial communities has been linked to numerous diseases, yet microbial interactions are little understood. This is due in part to the large number of bacteria, and the much larger number of interactions (easily in the millions), making experimental investigation very difficult at best and necessitating the nascent field of computational exploration through microbial correlation networks. We benchmark the performance of eight correlation techniques on simulated and real data in response to challenges specific to microbiome studies: fractional sampling of ribosomal RNA sequences, uneven sampling depths, rare microbes and a high proportion of zero counts. Also tested is the ability to distinguish signals from noise, and detect a range of ecological and time-series relationships. Finally, we provide specific recommendations for correlation technique usage. Although some methods perform better than others, there is still considerable need for improvement in current techniques.

Interleukin-6 (IL-6), also called 26-kd protein, hybridoma plasmacytoma growth factor, beta 2-interferon, or B cell stimulatory factor 2, is a recently described human cytokine with multiple growth and differentiation activities. Using a very sensitive bioassay based on the growth factor activity of this protein for B cell hybridomas, we found that IL-6 activity was significantly elevated in synovial fluid from patients with rheumatoid arthritis (RA) or other inflammatory arthritides, as compared with that in a group of patients with osteoarthritis. Moreover, IL-6 was detected in about one-third of the serum samples from patients with RA. In the latter group, we found a significant correlation between serum IL-6 activity and serum levels of C-reactive protein, alpha 1-acid glycoprotein, alpha 1-antitrypsin, fibrinogen, and haptoglobin, which indicates that IL-6 is related to disease activity in patients with RA.

Viruses influence ecosystems by modulating microbial population size, diversity, metabolic outputs, and gene flow. Here, we use quantitative double-stranded DNA (dsDNA) viral-fraction metagenomes (viromes) and whole viral community morphological data sets from 43 Tara Oceans expedition samples to assess viral community patterns and structure in the upper ocean. Protein cluster cataloging defined pelagic upper-ocean viral community pan and core gene sets and suggested that this sequence space is well-sampled. Analyses of viral protein clusters, populations, and morphology revealed biogeographic patterns whereby viral communities were passively transported on oceanic currents and locally structured by environmental conditions that affect host community structure. Together, these investigations establish a global ocean dsDNA viromic data set with analyses supporting the seed-bank hypothesis to explain how oceanic viral communities maintain high local diversity.