Pvclust: an R package for assessing the uncertainty in hierarchical clustering

Viral community analysis in a marine oxygen minimum zone indicates increased potential for viral manipulation of microbial physiological state

Microbial communities in oxygen minimum zones (OMZs) are known to have significant impacts on global biogeochemical cycles, but viral influence on microbial processes in these regions are much less studied. Here we provide baseline ecological patterns using microscopy and viral metagenomics from the Eastern Tropical North Pacific (ETNP) OMZ region that enhance our understanding of viruses in these climate-critical systems. While extracellular viral abundance decreased below the oxycline, viral diversity and lytic infection frequency remained high within the OMZ, demonstrating that viral influences on microbial communities were still substantial without the detectable presence of oxygen. Viral community composition was strongly related to oxygen concentration, with viral populations in low-oxygen portions of the water column being distinct from their surface layer counterparts. However, this divergence was not accompanied by the expected differences in viral-encoded auxiliary metabolic genes (AMGs) relating to nitrogen and sulfur metabolisms that are known to be performed by microbial communities in these low-oxygen and anoxic regions. Instead, several abundant AMGs were identified in the oxycline and OMZ that may modulate host responses to low-oxygen stress. We hypothesize that this is due to selection for viral-encoded genes that influence host survivability rather than modulating host metabolic reactions within the ETNP OMZ. Together, this study shows that viruses are not only diverse throughout the water column in the ETNP, including the OMZ, but their infection of microorganisms has the potential to alter host physiological state within these biogeochemically important regions of the ocean.

Comparative molecular genomic analyses of a spontaneous rhesus macaque model of mismatch repair-deficient colorectal cancer

Colorectal cancer (CRC) remains the third most common cancer in the US with 15% of cases displaying Microsatellite Instability (MSI) secondary to Lynch Syndrome (LS) or somatic hypermethylation of the MLH1 promoter. A cohort of rhesus macaques from our institution developed spontaneous mismatch repair deficient (MMRd) CRC with a notable fraction harboring a pathogenic germline mutation in MLH1 (c.1029C Collapse

Key Words Collapse

MESH Headings

• Animals
• Colonic Neoplasms/genetics
• Colorectal Neoplasms/genetics
• Colorectal Neoplasms/pathology
• Colorectal Neoplasms, Hereditary Nonpolyposis/genetics
• Colorectal Neoplasms, Hereditary Nonpolyposis/pathology
• DNA Methylation/genetics
• DNA Mismatch Repair/genetics
• Genomics
• Humans
• Macaca mulatta/genetics
• Mice
• Microsatellite Instability
• MutL Protein Homolog 1/genetics

Collapse

Grants

• P30 CA016672 NCI NIH HHS
• R24 OD011173 NIH HHS
• Feinberg Family Foundation
• MDACC Institutional Research Grant Program
• Cattleman for Cancer Research
• National Institutes of Health
• National Cancer Institute

Collapse

Affiliation(s)

Nejla Ozirmak Lermi Department of Clinical Cancer Prevention, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America School of Health Professions, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Stanton B. Gray Michale E. Keeling Center for Comparative Medicine and Research, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Charles M. Bowen Department of Clinical Cancer Prevention, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Laura Reyes-Uribe Department of Clinical Cancer Prevention, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Beth K. Dray Charles River Laboratories, Ashland, Ohio, United States of America
Nan Deng Department of Clinical Cancer Prevention, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
R. Alan Harris Human Genome Sequencing Center and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States of America
Muthuswamy Raveendran Human Genome Sequencing Center and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States of America
Fernando Benavides Department of Epigenetics and Molecular Carcinogenesis, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Carolyn L. Hodo Michale E. Keeling Center for Comparative Medicine and Research, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Melissa W. Taggart Department of Pathology, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Karen Colbert Maresso Department of Clinical Cancer Prevention, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Krishna M. Sinha Department of Clinical Cancer Prevention, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Jeffrey Rogers Human Genome Sequencing Center and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States of America
Eduardo Vilar Department of Clinical Cancer Prevention, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America Clinical Cancer Genetics Program, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America

Collapse

Multiplex component-based allergen macroarray test is useful to predict clinical reactivity to tree nuts in children

BACKGROUND

In tree nut (TN) allergy, singleplex tests showed the diagnostic utility of rAna o 3, rCor a 14/nCor a 9, and nJug r 1/nJug r 4 for cashew/pistachio, hazelnut, and walnut allergies, respectively. However, disadvantages of the tests include high costs and excessive blood sampling in multi-sensitized patients, and a limited number of components. We investigated the utility of a multiplex macroarray (i.e., the ALEX² test) in TN allergy.

METHODS

In 169 children, skin prick test, the component- and extract-specific IgEs of TNs were investigated for clinical reactivity and tolerance.

RESULTS

The predictors (AUC = 0.962-0.749) of clinical reactivity to cashew, pistachio, hazelnut, and walnut were rPis v 1/rAna o 3, rPis v 1/rAna o 3/nPis v 2/nPis v 3, rCor a 14/nCor a 11/nCor a 9, and nJug r 1/nJug r 2/nJug r 6/nJug r 4, respectively. More than 93% of the patients with clinical reactivity to pistachio/cashew, hazelnut and walnut had positivity of (≥0.3 kU_A/L) rPis v 1/rAna o 3, rCor a 14 and nJug r 1/nJug r 2, respectively. The highest accuracies of clinical reactivity to culprit nut were obtained with combination of rPis v 1, sIgE and SPT positivities for cashew/pistachio, rPis v 1 ≥ 1.0 kU_A/L for pistachio, rCor a 14 ≥ 1.0 kU_A/L for hazelnut and combination of nJug r 1 and nJug r 2 positivities for walnut, respectively. Also, higher concentrations of rPis v 1 (≥15.0 kU_A/L), rCor a 14 (≥5.0 kU_A/L) and nJug r 1/nJug r 2 (≥15.0 kU_A/L) had %100 specificity and PPV in predicting clinical reactivity to cashew, hazelnut and walnut, respectively.

CONCLUSIONS

Multiplex macroarray test is useful and reliable in the diagnosis of TN allergy in children, confirms and expands existing knowledge, and can be used as a stand-alone tool in the bottom-up diagnostic approach.

Forest Fire Causes and Motivations in Southern and South-Eastern Europe through the Perception of Experts: Contribution to Enhance the Current Policies

Forest fires causes and motivations are poorly understood in southern and south-eastern Europe. This research aims to identify how experts perceive the different causes of forest fires as defined in the classification proposed by the European Commission in 2013. A panel of experts (N = 271) was gathered from the EU Southern Member States (France, Greece, Italy, Portugal, and Spain) and from Central (Switzerland) and south-eastern Europe (Croatia, Serbia, Bosnia and Herzegovina, Republic of North Macedonia, and Turkey). Experts were asked to answer a questionnaire to score the importance of the 29 fire causes using a five point (1–5) Likert Scale. Agricultural burnings received the highest score, followed by Deliberate fire for profit, and Vegetation management. Most of the events stem from negligence, whereas malicious fire setting is arguably overestimated although there are differences among the countries. This research demonstrates the importance of different techniques to enhance the knowledge of the causes of the complex anthropogenic phenomenon of forest fire occurrence.

The Singular Evolution of Olea Genome Structure

The current view of plant genome evolution proposes that genome size has mainly been determined by polyploidisation and amplification/loss of transposons, with a minor role played by other repeated sequences, such as tandem repeats. In cultivated olive (Olea europaea subsp. europaea var. europaea), available data suggest a singular model of genome evolution, in which a massive expansion of tandem-repeated sequences accompanied changes in nuclear architecture. This peculiar scenario highlights the importance of focusing on Olea genus evolution, to shed light on mechanisms that led to its present genomic structure. Next-generation sequencing technologies, bioinformatics and in situ hybridisation were applied to study the genomic structure of five related Olea taxa, which originated at different times from their last common ancestor. On average, repetitive DNA in the Olea taxa ranged from ~59% to ~73% of the total genome, showing remarkable differences in terms of composition. Among repeats, we identified 11 major families of tandem repeats, with different abundances in the analysed taxa, five of which were novel discoveries. Interestingly, overall tandem repeat abundance was inversely correlated to that of retrotransposons. This trend might imply a competition in the proliferation of these repeat classes. Indeed, O. paniculata, the species closest to the Olea common ancestor, showed very few tandem-repeated sequences, while it was rich in long terminal repeat retrotransposons, suggesting that the amplification of tandem repeats occurred after its divergence from the Olea ancestor. Furthermore, some tandem repeats were physically localised in closely related O. europaea subspecies (i.e., cultivated olive and O. europaea subsp. cuspidata), which showed a significant difference in tandem repeats abundance. For 4 tandem repeats families, a similar number of hybridisation signals were observed in both subspecies, apparently indicating that, after their dissemination throughout the olive genome, these tandem repeats families differentially amplified maintaining the same positions in each genome. Overall, our research identified the temporal dynamics shaping genome structure during Olea speciation, which represented a singular model of genome evolution in higher plants.

Gene Expression over Time during Cell Transformation Due to Non-Genotoxic Carcinogen Treatment of Bhas 42 Cells

The Bhas 42 cell transformation assay (Bhas 42 CTA) is the first Organization for Economic Cooperation and Development (OECD)-certificated method used as a specific tool for the detection of the cell-transformation potential of tumor-promoting compounds, including non-genotoxic carcinogens (NGTxCs), as separate from genotoxic carcinogens. This assay offers the great advantage of enabling the phenotypic detection of oncotransformation. A key benefit of using the Bhas 42 CTA in the study of the cell-transformation mechanisms of tumor-promoting compounds, including non-genotoxic carcinogens, is that the cell-transformation potential of the chemical can be detected directly without treatment with a tumor-initiating compound since Bhas 42 cell line was established by transfecting the v-Ha-ras gene into a mouse fibroblast cloned cell line. Here, we analyzed the gene expression over time, using DNA microarrays, in Bhas 42 cells treated with the tumor-promoting compound 12-O-tetradecanoylphorbol-13-acetate (TPA), and NGTxC, with a total of three repeat experiments. This is the first paper to report on gene expression over time during the process of cell transformation with only a tumor-promoting compound. Pathways that were activated or inactivated during the process of cell transformation in the Bhas 42 cells treated with TPA were related not only directly to RAS but also to various pathways in the hallmarks of cancer.

Distinct Genomic Profiles Are Associated with Treatment Response and Survival in Ovarian Cancer

The majority of patients with ovarian cancer ultimately develop recurrent chemotherapy-resistant disease. Treatment stratification is mainly based on histological subtype and stage, prior response to platinum-based chemotherapy, and time to recurrent disease. Here, we integrated clinical treatment, treatment response, and survival data with whole-genome sequencing profiles of 132 solid tumor biopsies of metastatic epithelial ovarian cancer to explore genome-informed stratification opportunities. Samples from primary and recurrent disease harbored comparable numbers of single nucleotide variants and structural variants. Mutational signatures represented platinum exposure, homologous recombination deficiency, and aging. Unsupervised hierarchical clustering based on genomic input data identified specific ovarian cancer subgroups, characterized by homologous recombination deficiency, genome stability, and duplications. The clusters exhibited distinct response rates and survival probabilities which could thus potentially be used for genome-informed therapy stratification for more personalized ovarian cancer treatment.

Phytoextraction of heavy metals after application of bottom ash and municipal sewage sludge considering the risk of environmental pollution

Waste management to reduce the loss of natural resources has become a basis of sustainable development and a circular economy. When using waste, the heavy metal (HM) concentration must be taken into account since HMs can be potentially released to the environment, posing a toxicity threat. The aim of the study was thus to estimate the availability for plants of Cr, Ni, Cu, Zn, Cd, and Pb introduced into the soil with waste. We hypothesized that the prepared waste mixtures containing coal or biomass ash and municipal sewage sludge would reduce the environmental risk compared to the studied waste used separately. The research was conducted during a 6-year field experiment with grasses and legumes. HM concentration in soil, waste, and plant biomass; tolerance index; and uptake of HMs by plants were measured. The ash-sludge mixtures had a more favourable effect on the soil in terms of pH_KCl, TOC, total nitrogen, and total exchangeable bases than the waste used separately. This provided beneficial conditions for plant growth and development. Consequently, the ash-sludge mixtures increased the plant yield as compared to ash alone, while the mixture containing the biomass ash also enhanced the yield in relation to the sewage sludge. The study showed that the mixtures allowed for a reduction of environmental risk arising from the HM input with waste to the soil. It was proven that HM availability for plants could be beneficially modified by mixing waste. Combining the coal ash with the sewage sludge is particularly recommended, owing to the unfavourable properties of coal ash for plants. The application of the higher dose of the coal ash-sludge mixture showed a better effect than the lower dose, while the influence of both doses of the biomass ash-sludge mixture was similar. Under the ash-sludge treatment, plants took up more HM than under the ash used separately, and the HM concentration in the obtained biomass did not generally exceed that observed under single wastes. This should reduce the accumulation of HMs in the soil during a long-term use of the waste and facilitates the utilisation of the produced biomass.

BinMat: A molecular genetics tool for processing binary data obtained from fragment analysis in R

Processing and visualising trends in the binary data (presence or absence of electropherogram peaks), obtained from fragment analysis methods in molecular biology, can be a time-consuming and often cumbersome process. Scoring and analysing binary data (from methods, such as AFLPs, ISSRs and RFLPs) entail complex workflows that require a high level of computational and bioinformatic skills. The application presented here (BinMat) is a free, open-source and user-friendly R Shiny programme (https://clarkevansteenderen.shinyapps.io/BINMAT/) that automates the analysis pipeline on one platform. It is also available as an R package on the Comprehensive R Archive Network (CRAN) (https://cran.r-project.org/web/packages/BinMat/index.html). BinMat consolidates replicate sample pairs of binary data into consensus reads, produces summary statistics and allows the user to visualise their data as ordination plots and clustering trees without having to use multiple programmes and input files or rely on previous programming experience.

Genomic and Functional Variation of the Chlorophyll d-Producing Cyanobacterium Acaryochloris marina

The Chlorophyll d-producing cyanobacterium Acaryochloris marina is widely distributed in marine environments enriched in far-red light, but our understanding of its genomic and functional diversity is limited. Here, we take an integrative approach to investigate A. marina diversity for 37 strains, which includes twelve newly isolated strains from previously unsampled locations in Europe and the Pacific Northwest of North America. A genome-wide phylogeny revealed both that closely related A. marina have migrated within geographic regions and that distantly related A. marina lineages can co-occur. The distribution of traits mapped onto the phylogeny provided evidence of a dynamic evolutionary history of gene gain and loss during A. marina diversification. Ancestral genes that were differentially retained or lost by strains include plasmid-encoded sodium-transporting ATPase and bidirectional NiFe-hydrogenase genes that may be involved in salt tolerance and redox balance under fermentative conditions, respectively. The acquisition of genes by horizontal transfer has also played an important role in the evolution of new functions, such as nitrogen fixation. Together, our results resolve examples in which genome content and ecotypic variation for nutrient metabolism and environmental tolerance have diversified during the evolutionary history of this unusual photosynthetic bacterium.

Effects of Competitive ELISA-Positive Results of Piroplasmosis on the Performance of Endurance Horses

Simple Summary

Equine piroplasmosis (EP) is a common infectious disease in southern Europe. To better understand the impact and influence of EP on the performance of endurance horses, we collected blood samples from national elite horses during different endurance competitions. The horses were tested against piroplasmosis, and several blood parameters related to performance were also evaluated. It seems that horses without clinical signs of piroplasmosis can participate without performance impairment in competitions of up to 80 km. Although it is recommended that longer distance competitions should be further evaluated, this is the first step for decision-making by organizers and participants in this sport.

Abstract

Endurance is an increasingly popular equestrian sport. However, in southern Europe, there is a high prevalence of horses that are asymptomatic carriers of equine piroplasmosis (EP), a tick-borne disease that could affect their performance. This study aimed to evaluate the impact and influence of EP on the performance of endurance horses. Blood samples were collected from 40 horses in Extremadura, Spain, before and after a race, in different national elite horse endurance competitions. Hematological and biochemical parameters and EP seroprevalence were analysed by competitive enzyme-linked immunosorbent assay. The global seroprevalence of EP was 70%, with 27 horses testing positive for Theileria equi (67.5%) and three (7.5%) for Babesia caballi, with two of these horses (5%) positive for both. Approximately 82.5% of the horses (33 of 40) completed the competition, with no influence on performance or position achieved in those with subclinical parasitosis. There were also no significant differences in hematological or biochemical values between seropositive and seronegative horses. The data suggest that horses without clinical signs of EP can participate without performance impairment in competitions of up to 80 km. Although it is recommended that longer distance competitions should be further evaluated, this is the first step for decision-making by organizers and participants in this sport.

Genome-wide association study reveals white lupin candidate gene involved in anthracnose resistance

GWAS identifies candidate gene controlling resistance to anthracnose disease in white lupin. White lupin (Lupinus albus L.) is a promising grain legume to meet the growing demand for plant-based protein. Its cultivation, however, is severely threatened by anthracnose disease caused by the fungal pathogen Colletotrichum lupini. To dissect the genetic architecture for anthracnose resistance, genotyping by sequencing was performed on white lupin accessions collected from the center of domestication and traditional cultivation regions. GBS resulted in 4611 high-quality single-nucleotide polymorphisms (SNPs) for 181 accessions, which were combined with resistance data observed under controlled conditions to perform a genome-wide association study (GWAS). Obtained disease phenotypes were shown to highly correlate with overall three-year disease assessments under Swiss field conditions (r > 0.8). GWAS results identified two significant SNPs associated with anthracnose resistance on gene Lalb_Chr05_g0216161 encoding a RING zinc-finger E3 ubiquitin ligase which is potentially involved in plant immunity. Population analysis showed a remarkably fast linkage disequilibrium decay, weak population structure and grouping of commercial varieties with landraces, corresponding to the slow domestication history and scarcity of modern breeding efforts in white lupin. Together with 15 highly resistant accessions identified in the resistance assay, our findings show promise for further crop improvement. This study provides the basis for marker-assisted selection, genomic prediction and studies aimed at understanding anthracnose resistance mechanisms in white lupin and contributes to improving breeding programs worldwide.

The pangenome of banana highlights differences between genera and genomes

Banana (Musaceae family) has a complex genetic history and includes a genus Musa with a variety of cultivated clones with edible fruits, Ensete species that are grown for their edible corm, and monospecific Musella whose generic status has been questioned. The most commonly exported banana cultivars belong to Cavendish, a subgroup of Musa triploid cultivars, which is under threat by fungal pathogens, though there are also related species M. balbisiana Colla (B genome), M. textilis Née (T genome), and M. schizocarpa N. W. Simmonds (S genome), along with hybrids of these genomes, which potentially host genes of agronomic interest. Here we present the first cross-genus pangenome of banana, which contains representatives of the Musa and Ensete genera. Clusters based on gene presence-absence variation (PAV) clearly separate Musa and Ensete, while Musa is split further based on species. These results present the first pangenome study across genus boundaries and identifies genes that differentiate between Musaceae species, information that may support breeding programs in these crops.

Large-Scale Meta-Longitudinal Microbiome Data with a Known Batch Factor

Data contamination in meta-approaches where multiple biological samples are combined considerably affects the results of subsequent downstream analyses, such as differential abundance tests comparing multiple groups at a fixed time point. Little has been thoroughly investigated regarding the impact of the lurking variable of various batch sources, such as different days or different laboratories, in more complicated time series experimental designs, for instance, repeatedly measured longitudinal data and metadata. We highlight that the influence of batch factors is significant on subsequent downstream analyses, including longitudinal differential abundance tests, by performing a case study of microbiome time course data with two treatment groups and a simulation study of mimic microbiome longitudinal counts.

Butterfly eyespots evolved via cooption of an ancestral gene-regulatory network that also patterns antennae, legs, and wings

Where do butterfly eyespots come from? One of the long-standing questions in the field of evolution concerns addressing where novel complex traits come from. Here we show that butterfly eyespots, a novel complex trait, likely originated from the redeployment of a preexisting gene-regulatory network regulating antennae, legs, and wings, to novel locations on the wing.

Butterfly eyespots are beautiful novel traits with an unknown developmental origin. Here we show that eyespots likely originated via cooption of parts of an ancestral appendage gene-regulatory network (GRN) to novel locations on the wing. Using comparative transcriptome analysis, we show that eyespots cluster most closely with antennae, relative to multiple other tissues. Furthermore, three genes essential for eyespot development, Distal-less (Dll), spalt (sal), and Antennapedia (Antp), share similar regulatory connections as those observed in the antennal GRN. CRISPR knockout of cis-regulatory elements (CREs) for Dll and sal led to the loss of eyespots, antennae, legs, and also wings, demonstrating that these CREs are highly pleiotropic. We conclude that eyespots likely reused an ancient GRN for their development, a network also previously implicated in the development of antennae, legs, and wings.

Metagenomic Assessment of DNA Viral Diversity in Freshwater Sponges, Baikalospongia bacillifera

Sponges (type Porifera) are multicellular organisms that give shelter to a variety of microorganisms: fungi, algae, archaea, bacteria, and viruses. The studies concerning the composition of viral communities in sponges have appeared rather recently, and the diversity and role of viruses in sponge holobionts remain largely undisclosed. In this study, we assessed the diversity of DNA viruses in the associated community of the Baikal endemic sponge, Baikalospongia bacillifera, using a metagenomic approach, and compared the virome data from samples of sponges and Baikal water (control sample). Significant differences in terms of taxonomy, putative host range of identified scaffolds, and functional annotation of predicted viral proteins were revealed in viromes of sponge B. bacillifera and the Baikal water. This is the evidence in favor of specificity of viral communities in sponges. The diversity shift of viral communities in a diseased specimen, in comparison with a visually healthy sponge, probably reflects the changes in the composition of microbial communities in affected sponges. We identified many viral genes encoding the proteins with metabolic functions; therefore, viruses in Baikal sponges regulate the number and diversity of their associated community, and also take a part in the vital activity of the holobiont, and this is especially significant in the case of damage (or disease) of these organisms in unfavorable conditions. When comparing the Baikal viromes with similar datasets of marine sponge (Ianthella basta), in addition to significant differences in the taxonomic and functional composition of viral communities, we revealed common scaffolds/virotypes in the cross-assembly of reads, which may indicate the presence of some closely related sponge-specific viruses in marine and freshwater sponges.

Bermudagrass Cultivars with Different Tolerance to Nematode Damage Are Characterized by Distinct Fungal but Similar Bacterial and Archaeal Microbiomes

Turfgrass landscapes have expanded rapidly in recent decades and are a major vegetation type in urbanizing ecosystems. While turfgrass areas provide numerous ecosystem services in urban environments, ecological side effects from intensive management are raising concerns regarding their sustainability. One potentially promising approach to ameliorate the ecological impact and decrease the use of agricultural chemicals is to take advantage of naturally evolved turfgrass-associated microbes by harnessing beneficial services provided by microbiomes. Unfortunately, especially compared to agricultural crops, the microbiomes of turfgrasses are not well understood. Here, we analyzed microbial communities inhabiting the leaf and root endospheres as well as soil in two bermudagrass cultivars, 'Latitude 36' and 'TifTuf', which exhibit distinct tolerance to nematode damage, with the goal of identifying potential differences in the microbiomes that might explain their distinct phenotype. We used 16S rRNA gene V4 and ITS2 amplicon sequencing to characterize the microbiomes in combination with microbial cultivation efforts to identify potentially beneficial endophytic fungi and bacteria. Our results show that Latitude 36 and TifTuf showed markedly different fungal microbiomes, each harboring unique taxa from Ascomycota and Glomeromycota, respectively. In contrast, less difference was observed from bacterial and archaeal microbiomes, which were dominated by Bacteroidetes and Thaumarchaeota, respectively. The TifTuf microbiomes exhibited lower microbial diversity compared to Latitude 36. Many sequences could not be classified to a higher taxonomic resolution, indicating a relatively high abundance of hitherto undescribed microorganisms. Our results provide new insights into the structure and composition of turfgrass microbiomes but also raise important questions regarding the functional attributes of key taxa.

Norway spruce deploys tissue-specific responses during acclimation to cold

Climate change in the conifer-dominated boreal forest is expected to lead to warmer but more dynamic winter air temperatures, reducing the depth and duration of snow cover and lowering winter soil temperatures. To gain insight into the mechanisms that have enabled conifers to dominate extreme cold environments, we performed genome-wide RNA-Seq analysis from needles and roots of non-dormant two-year Norway spruce (Picea abies (L.) H. Karst), and contrasted these response to herbaceous model Arabidopsis We show that the main transcriptional response of Norway spruce needles exposed to cold was delayed relative to Arabidopsis, and this delay was associated with slower development of freezing tolerance. Despite this difference in timing, Norway spruce principally utilizes early response transcription factors (TFs) belonging to the same gene families as Arabidopsis, indicating broad evolutionary conservation of cold response networks. In keeping with their different metabolic and developmental states, needles and root of Norway spruce showed contrasting results. Regulatory network analysis identified both conserved TFs with known roles in cold acclimation (e.g. homologs of ICE1, AKS3, and of the NAC and AP2/ERF superfamilies), but also a root-specific bHLH101 homolog, providing functional insights into cold stress response strategies in Norway spruce.

Highly diverse flavobacterial phages isolated from North Sea spring blooms

It is generally recognized that phages are a mortality factor for their bacterial hosts. This could be particularly true in spring phytoplankton blooms, which are known to be closely followed by a highly specialized bacterial community. We hypothesized that phages modulate these dense heterotrophic bacteria successions following phytoplankton blooms. In this study, we focused on Flavobacteriia, because they are main responders during these blooms and have an important role in the degradation of polysaccharides. A cultivation-based approach was used, obtaining 44 lytic flavobacterial phages (flavophages), representing twelve new species from two viral realms. Taxonomic analysis allowed us to delineate ten new phage genera and ten new families, from which nine and four, respectively, had no previously cultivated representatives. Genomic analysis predicted various life styles and genomic replication strategies. A likely eukaryote-associated host habitat was reflected in the gene content of some of the flavophages. Detection in cellular metagenomes and by direct-plating showed that part of these phages were actively replicating in the environment during the 2018 spring bloom. Furthermore, CRISPR/Cas spacers and re-isolation during two consecutive years suggested that, at least part of the new flavophages are stable components of the microbial community in the North Sea. Together, our results indicate that these diverse flavophages have the potential to modulate their respective host populations.

Mapping co-regulatory interactions among ligand-binding sites in ryanodine receptor 1

Ryanodine receptor 1 (RyR1) is an intracellular calcium ion (Ca²⁺ ) release channel required for skeletal muscle contraction. Although cryo-electron microscopy identified binding sites of three coactivators Ca²⁺ , ATP, and caffeine (CFF), the mechanism of co-regulation and synergy of these activators is unknown. Here, we report allosteric connections among the three ligand-binding sites and pore region in (i) Ca²⁺ bound-closed, (ii) ATP/CFF bound-closed, (iii) Ca²⁺ /ATP/CFF bound-closed, and (iv) Ca²⁺ /ATP/CFF bound-open RyR1 states. We identified two dominant networks of interactions that mediate communication between the Ca²⁺ -binding site and pore region in Ca²⁺ bound-closed state, which partially overlapped with the pore communications in ATP/CFF bound-closed RyR1 state. In Ca²⁺ /ATP/CFF bound-closed and -open RyR1 states, co-regulatory interactions were analogous to communications in the Ca²⁺ bound-closed and ATP/CFF bound-closed states. Both ATP- and CFF-binding sites mediate communication between the Ca²⁺ -binding site and the pore region in Ca²⁺ /ATP/CFF bound-open RyR1 structure. We conclude that Ca²⁺ , ATP, and CFF propagate their effects to the pore region through a network of overlapping interactions that mediate allosteric control and molecular synergy in channel regulation.

Modular evolution of secretion systems and virulence plasmids in a bacterial species complex

BACKGROUND

Many named species as defined in current bacterial taxonomy correspond to species complexes. Uncertainties regarding the organization of their genetic diversity challenge research efforts. We utilized the Agrobacterium tumefaciens species complex (a.k.a. Agrobacterium biovar 1), a taxon known for its phytopathogenicity and applications in transformation, as a study system and devised strategies for investigating genome diversity and evolution of species complexes.

RESULTS

We utilized 35 genome assemblies, including 14 newly generated ones, to achieve a phylogenetically balanced sampling of A. tumefaciens. Our genomic analysis suggested that the 10 genomospecies described previously are distinct biological species and supported a quantitative guideline for species delineation. Furthermore, our inference of gene content and core-genome phylogeny allowed for investigations of genes critical in fitness and ecology. For the type VI secretion system (T6SS) involved in interbacterial competition and thought to be conserved, we detected multiple losses and one horizontal gene transfer. For the tumor-inducing plasmids (pTi) and pTi-encoded type IV secretion system (T4SS) that are essential for agrobacterial phytopathogenicity, we uncovered novel diversity and hypothesized their involvement in shaping this species complex. Intriguingly, for both T6SS and T4SS, genes encoding structural components are highly conserved, whereas extensive diversity exists for genes encoding effectors and other proteins.

CONCLUSIONS

We demonstrate that the combination of a phylogeny-guided sampling scheme and an emphasis on high-quality assemblies provides a cost-effective approach for robust analysis in evolutionary genomics. We show that the T6SS VgrG proteins involved in specific effector binding and delivery can be classified into distinct types based on domain organization. The co-occurrence patterns of VgrG-associated domains and the neighboring genes that encode different chaperones/effectors can be used to infer possible interacting partners. Similarly, the associations between plant host preference and the pTi type among these strains can be used to infer phenotype-genotype correspondence. Our strategies for multi-level investigations at scales that range from whole genomes to intragenic domains and phylogenetic depths from between- to within-species are applicable to other bacteria. Furthermore, modularity observed in the molecular evolution of genes and domains is useful for inferring functional constraints and informing experimental works.

Large-Scale Structure-Based Screening of Potential T Cell Cross-Reactivities Involving Peptide-Targets From BCG Vaccine and SARS-CoV-2

Although not being the first viral pandemic to affect humankind, we are now for the first time faced with a pandemic caused by a coronavirus. The Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has been responsible for the COVID-19 pandemic, which caused more than 4.5 million deaths worldwide. Despite unprecedented efforts, with vaccines being developed in a record time, SARS-CoV-2 continues to spread worldwide with new variants arising in different countries. Such persistent spread is in part enabled by public resistance to vaccination in some countries, and limited access to vaccines in other countries. The limited vaccination coverage, the continued risk for resistant variants, and the existence of natural reservoirs for coronaviruses, highlight the importance of developing additional therapeutic strategies against SARS-CoV-2 and other coronaviruses. At the beginning of the pandemic it was suggested that countries with Bacillus Calmette-Guérin (BCG) vaccination programs could be associated with a reduced number and/or severity of COVID-19 cases. Preliminary studies have provided evidence for this relationship and further investigation is being conducted in ongoing clinical trials. The protection against SARS-CoV-2 induced by BCG vaccination may be mediated by cross-reactive T cell lymphocytes, which recognize peptides displayed by class I Human Leukocyte Antigens (HLA-I) on the surface of infected cells. In order to identify potential targets of T cell cross-reactivity, we implemented an in silico strategy combining sequence-based and structure-based methods to screen over 13,5 million possible cross-reactive peptide pairs from BCG and SARS-CoV-2. Our study produced (i) a list of immunogenic BCG-derived peptides that may prime T cell cross-reactivity against SARS-CoV-2, (ii) a large dataset of modeled peptide-HLA structures for the screened targets, and (iii) new computational methods for structure-based screenings that can be used by others in future studies. Our study expands the list of BCG peptides potentially involved in T cell cross-reactivity with SARS-CoV-2-derived peptides, and identifies multiple high-density "neighborhoods" of cross-reactive peptides which could be driving heterologous immunity induced by BCG vaccination, therefore providing insights for future vaccine development efforts.

New Hyperspectral Procedure to Discriminate Intertidal Macroalgae

The recent development and miniaturization of hyperspectral sensors embedded in drones has allowed the acquisition of hyperspectral images with high spectral and spatial resolution. The characteristics of both the embedded sensors and drones (viewing angle, flying altitude, resolution) create opportunities to consider the use of hyperspectral imagery to map and monitor macroalgae communities. In general, the overflight of the areas to be mapped is conconmittently associated accompanied with measurements carried out in the field to acquire the spectra of previously identified objects. An alternative to these simultaneous acquisitions is to use a hyperspectral library made up of pure spectra of the different species in place, that would spare field acquisition of spectra during each flight. However, the use of such a technique requires developed appropriate procedure for testing the level of species classification that can be achieved, as well as the reproducibility of the classification over time. This study presents a novel classification approach based on the use of reflectance spectra of macroalgae acquired in controlled conditions. This overall approach developed is based on both the use of the spectral angle mapper (SAM) algorithm applied on first derivative hyperspectral data. The efficiency of this approach has been tested on a hyperspectral library composed of 16 macroalgae species, and its temporal reproducibility has been tested on a monthly survey of the spectral response of different macro-algae species. In addition, the classification results obtained with this new approach were also compared to the results obtained through the use of the most recent and robust procedure published. The classification obtained shows that the developed approach allows to perfectly discriminate the different phyla, whatever the period. At the species level, the classification approach is less effective when the individuals studied belong to phylogenetically close species (i.e., Fucus spiralis and Fucus serratus).

Genome survey sequencing and characterization of simple sequence repeat (SSR) markers in Platostoma palustre (Blume) A.J.Paton (Chinese mesona)

Platostoma palustre (Blume) A.J.Paton is an annual herbaceous persistent plant of the Labiatae family. However, there is a lack of genomic data for this plant, which severely restricts its genetic improvement. In this study, we performed genome survey sequencing of P. palustre and developed simple sequence repeat (SSR) markers based on the resulting sequence. K-mer analysis revealed that the assembled genome size was approximately 1.21 Gb. A total of 15,498 SSR motifs were identified and characterized in this study; among them, dinucleotide, and hexanucleotide repeats had the highest and lowest, respectively. Among the dinucleotide repeat motifs, AT/TA repeat motifs were the most abundant, and GC/CG repeat motifs were rather rare, accounting for 44.28% and 0.63%, respectively. Genetic similarity coefficient analysis by the UPMGA methods clustered 12 clones, of P. palustre and related species into two subgroups. These results provide helpful information for further research on P. palustre resources and variety improvements.

Unveiling biogeographical patterns of the ichthyofauna in the Tuichi basin, a biodiversity hotspot in the Bolivian Amazon, using environmental DNA

To date, more than 2400 valid fish species have been recorded in the Amazon basin. However, some regions remain poorly documented. This is the case in the Beni basin and in particular in one of its main sub-basins, the Tuichi, an Andean foothills rivers flowing through the Madidi National Park in the Bolivian Amazonia. The knowledge of its ichthyological diversity is, however, essential for the management and protection of aquatic ecosystems, which are threatened by the development of infrastructures (dams, factories and cities), mining and deforestation. Environmental DNA (eDNA) has been relatively little used so far in the Amazon basin. We sampled eDNA from water in 34 sites in lakes and rivers in the Beni basin including 22 sites in the Tuichi sub-basin, during the dry season. To assess the biogeographical patterns of the amazonian ichthyofauna, we implemented a metabarcoding approach using two pairs of specific primers designed and developed in our laboratory to amplify two partially overlapping CO1 fragments, one of 185bp and another of 285bp. We detected 252 fish taxa (207 at species level) among which 57 are newly identified for the Beni watershed. Species compositions are significantly different between lakes and rivers but also between rivers according to their hydrographic rank and altitude. Furthermore, the diversity patterns are related to the different hydro-ecoregions through which the Tuichi flows. The eDNA approach makes it possible to identify and complete the inventory of the ichthyofauna in this still poorly documented Amazon basin. However, taxonomic identification remains constrained by the lack of reference barcodes in public databases and does not allow the assignment of all OTUs. Our results can be taken into account in conservation and management strategies and could serve as a baseline for future studies, including on other Andean tributaries.

Interactome and evolutionary conservation of Dictyostelid small GTPases and their direct regulators

GTP binding proteins known as small GTPases make up one of the largest groups of regulatory proteins and control almost all functions of living cells. Their activity is under, respectively, positive and negative regulation by guanine nucleotide exchange factors (GEFs) and GTPase activating proteins (GAPs), which together with their upstream regulators and the downstream targets of the small GTPases form formidable signalling networks. While genomics has revealed the large size of the GTPase, GEF and GAP repertoires, only a small fraction of their interactions and functions have yet been experimentally explored. Dictyostelid social amoebas have been particularly useful in unravelling the roles of many proteins in the Rac-Rho and Ras-Rap families of GTPases in directional cell migration and regulation of the actin cytoskeleton. Genomes and cell-type specific and developmental transcriptomes are available for Dictyostelium species that span the 0.5 billion years of evolution of the group from their unicellular ancestors. In this work, we identified all GTPases, GEFs and GAPs from genomes representative of the four major taxon groups and investigated their phylogenetic relationships and evolutionary conservation and changes in their functional domain architecture and in their developmental and cell-type specific expression. We performed a hierarchical cluster analysis of the expression profiles of the ~2000 analysed genes to identify putative interacting sets of GTPases, GEFs and GAPs, which highlight sets known to interact experimentally and many novel combinations. This work represents a valuable resource for research into all fields of cellular regulation.

Dysregulated BMP2 in the Placenta May Contribute to Early-Onset Preeclampsia by Regulating Human Trophoblast Expression of Extracellular Matrix and Adhesion Molecules

Many pregnancy disorders, including early-onset preeclampsia (EOPE), are associated with defects in placental trophoblast cell invasion and differentiation during early placental development. Bone morphogenetic protein 2 (BMP2) belongs to the TGF-β superfamily and controls various physiological and developmental processes. However, the expression of BMP2 in the placenta and underlying molecular mechanisms of how BMP2 regulates trophoblast function remain unclear. In this study, we analyzed several publicly available microarray and RNA-seq datasets and revealed differences in expression of TGF-β superfamily members between gestational age-matched non-preeclamptic control and EOPE placentas. Importantly, BMP2 levels were significantly reduced in EOPE placentas compared with controls, and RNAscope in situ hybridization further demonstrated BMP2 expression was disrupted in EOPE placental villi. To explore the molecular mechanisms of BMP2-regulated early trophoblast differentiation, we examined BMP2 expression in first-trimester human placenta and found it to be localized to all subtypes of trophoblasts and the decidua. RNA-seq analysis on control and BMP2-treated primary human trophoblast cells identified 431 differentially expressed genes, including several canonical TGF-β/BMP signaling targets (BAMBI, ID1, INHBA, IGFBP3). Gene ontology annotations revealed that differentially expressed genes were involved in cell adhesion and extracellular matrix organization. Furthermore, we identified adhesion molecule with IgG-like domain 2 (AMIGO2) as a novel target for BMP2 that contributed to BMP2-induced trophoblast invasion and endothelial-like tube formation. Overall, our findings provide insight into the molecular processes controlled by BMP2 during early placental development that may contribute to the pathogenesis of EOPE.

Implications of divergence of methionine adenosyltransferase in archaea

Methionine adenosyltransferase (MAT) catalyzes the biosynthesis of S-adenosyl methionine from l-methionine and ATP. MAT enzymes are ancient, believed to share a common ancestor, and are highly conserved in all three domains of life. However, the sequences of archaeal MATs show considerable divergence compared with their bacterial and eukaryotic counterparts. Furthermore, the structural significance and functional significance of this sequence divergence are not well understood. In the present study, we employed structural analysis and ancestral sequence reconstruction to investigate archaeal MAT divergence. We observed that the dimer interface containing the active site (which is usually well conserved) diverged considerably between the bacterial/eukaryotic MATs and archaeal MAT. A detailed investigation of the available structures supports the sequence analysis outcome: The protein domains and subdomains of bacterial and eukaryotic MAT are more similar than those of archaea. Finally, we resurrected archaeal MAT ancestors. Interestingly, archaeal MAT ancestors show substrate specificity, which is lost during evolution. This observation supports the hypothesis of a common MAT ancestor for the three domains of life. In conclusion, we have demonstrated that archaeal MAT is an ideal system for studying an enzyme family that evolved differently in one domain compared with others while maintaining the same catalytic activity.

Molecular prevalence of resistance determinants, virulence factors and capsular serotypes among colistin resistance carbapenemase producing Klebsiella pneumoniae: a multi-centric retrospective study

The emergence of colistin-carbapenem-resistant Klebsiella pneumoniae (CCR-Kp) in bloodstream infection results in high mortality, and virulence factor contributes further to the difficulty of treatment. A total of 158 carbapenem-resistant K. pneumoniae (CRKP) isolates causing bloodstream infection were collected from three Indian tertiary care hospitals during the 9-month study period, of which 27 isolates exhibited resistance to both colistin and carbapenem antibiotics. In this study, all the strains were characterized for antimicrobial resistance, virulence factors and capsular serotypes that facilitate the development of colistin and carbapenem-resistant K.pneumoniae (CCR-Kp) in bloodstream infection. Fourteen isolates displayed extremely drug resistance (XDR), susceptible only to tigecycline, and the remaining 13 isolates displayed multidrug resistance (MDR). The gene prevalence analysis for CCR-Kp isolates showed the predominance of bla KPC (81.48%) followed by bla NDM (62.96%), bla VIM (37.03%) and bla IMP (18.51%) genes. The distribution of virulence genes was found to be fimH (81.48%), wabG (59.25%), mrkD (55.56%), entB (48.15%), irp1 (33.33%), and rmpA (18.52%). The capsular serotypes K1, K2, K5 and K54 have been identified in 16 isolates. The absence of plasmid-mediated colistin resistance (mcr) genes implies the involvement of other mechanisms. The ERIC and (GTG)5 molecular typing methods detected 18 and 22 distinct clustering patterns among the CCR-Kp isolates, respectively. A strong correlation between ERIC and (GTG)5 genotyping method was established with antimicrobial resistance patterns and virulence determinants at P < 0.05, while no correlation was found with capsular serotyping. Similar virulence and resistance typing among the isolates suggest hospital-acquired infection in a health care setup. These outcomes will advance our awareness of CCR-Kp outbreaks associated with tertiary care hospitals and help forecast their occurrence in the near future.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-03056-4.

Genetic Diversity and Differentiation of Eleven Medicago Species from Campania Region Revealed by Nuclear and Chloroplast Microsatellites Markers

The species belonging to the genus Medicago are considered a very important genetic resource at global level both for planet’s food security and for sustainable rangelands management. The checklist of the Italian flora (2021) includes a total number of 40 Medicago species for Italy, and 27 for Campania region, with a number of doubtful records or related to species no more found in the wild. In this study, 10 Medicago species native to Campania region, and one archaeophyte (M. sativa), identified by means of morphological diagnostic characters, were analyzed in a blind test to assay the efficacy of nine microsatellite markers (five cp-SSRs and four n-SSRs). A total number of 33 individuals from 6 locations were sampled and genotyped. All markers were polymorphic, 40 alleles were obtained with n-SSRs ranging from 8–12 alleles per locus with an average of 10 alleles per marker, PIC values ranged from 0.672 to 0.847, and the most polymorphic SSR was MTIC 564. The cp-SSRs markers were highly polymorphic too; PIC values ranged from 0.644 to 0.891 with an average of 0.776, the most polymorphic cp-SSR was CCMP10. 56 alleles were obtained with cp-SSRs ranging from 7 to 17 alleles per locus with an average of 11. AMOVA analysis with n-SSR markers highlighted a great level of genetic differentiation among the 11 species, with a statistically significant fixation index (FST). UPGMA clustering and Bayesian-based population structure analysis assigned these 11 species to two main clusters, but the distribution of species within clusters was not the same for the two analyses. In conclusion, our results demonstrated that the combination of the used SSRs well distinguished the 11 Medicago species. Moreover, our results demonstrated that the use of a limited number of SSRs might be considered for further genetic studies on other Medicago species.

Simple and robust methylation test for risk stratification of patients with juvenile myelomonocytic leukemia

Juvenile myelomonocytic leukemia (JMML) is a rare myelodysplastic/myeloproliferative neoplasm that develops during infancy and early childhood. The array-based international consensus definition of DNA methylation has recently classified patients with JMML into the following 3 groups: high (HM), intermediate (IM), and low methylation (LM). To develop a simple and robust methylation clinical test, 137 patients with JMML were analyzed using the Digital Restriction Enzyme Analysis of Methylation (DREAM), which is a next-generation sequencing-based methylation analysis. Unsupervised consensus clustering of the discovery cohort (n = 99) using DREAM data identified HM (HM_DREAM; n = 35) and LM subgroups (LM_DREAM; n = 64). Of the 98 cases that could be compared with the international consensus classification, 90 HM (n = 30) and LM (n = 60) cases had 100% concordance with DREAM clustering results. Of the remaining 8 cases comprising the IM group, 4 were classified as belonging to the HM_DREAM group and 4 to the LM_DREAM group. A machine-learning classifier was successfully constructed using a support vector machine (SVM), which divided the validation cohort (n = 38) into HM (HM_SVM, n = 18) and LM (LM_SVM; n = 20) groups. Patients with the HM_SVM profile had a significantly poorer 5-year overall survival rate than those with the LM_SVM profile. In conclusion, we developed a robust methylation test using DREAM for patients with JMML. This simple and straightforward test can be easily incorporated into diagnosis to generate a methylation classification for patients so they can receive risk-adapted treatment in the context of future clinical trials.

Conservation and over-representation of G-quadruplex sequences in regulatory regions of mitochondrial DNA across distinct taxonomic sub-groups

G-quadruplexes have important regulatory roles in the nuclear genome but their distribution and potential roles in mitochondrial DNA (mtDNA) are poorly understood. We analysed 11883 mtDNA sequences from 18 taxonomic sub-groups and identified their frequency and location within mtDNA. Large differences in both the frequency and number of putative quadruplex-forming sequences (PQS) were observed amongst all the organisms and PQS frequency was negatively correlated with an increase in evolutionary age. PQS were over-represented in the 3'UTRs, D-loops, replication origins, and stem loops, indicating regulatory roles for quadruplexes in mtDNA. Variations of the G-quadruplex-forming sequence in the conserved sequence block II (CSBII) region of the human D-loop were conserved amongst other mammals, amphibians, birds, reptiles, and fishes. This D-loop PQS was conserved in the duplicated control regions of some birds and reptiles, indicating its importance to mitochondrial function. The guanine tracts in these PQS also displayed significant length heterogeneity and the length of these guanine tracts were generally longest in bird mtDNA. This information provides further insights into how G4s may contribute to the regulation and function of mtDNA and acts as a database of information for future studies investigating mitochondrial G4s in organisms other than humans.

Saline lakes on the Qinghai-Tibet Plateau harbor unique viral assemblages mediating microbial environmental adaption

The highest plateau on Earth, Qinghai-Tibet Plateau, contains thousands of lakes with broad salinity and diverse and unique microbial communities. However, little is known about their co-occurring viruses. Herein, we identify 4,560 viral Operational Taxonomic Units (vOTUs) from six viromes of three saline lakes on Qinghai-Tibet Plateau, with less than 1% that could be classified. Most of the predicted vOTUs were associated with the dominant bacterial and archaeal phyla. Virus-encoded auxiliary metabolic genes suggest that viruses influence microbial metabolisms of carbon, nitrogen, sulfur, and lipid; the antibiotic resistance mediation; and their salinity adaption. The six viromes clustered together with the ice core viromes and bathypelagic ocean viromes and might represent a new viral habitat. This study has revealed the unique characteristics and potential ecological roles of DNA viromes in the lakes of the highest plateau and established a foundation for the recognition of the viral roles in plateau lake ecosystems.

Development of SSR Databases Available for Both NGS and Capillary Electrophoresis in Apple, Pear and Tea

Developing new varieties in fruit and tea breeding programs is very costly and labor-intensive. Thus, establishing a variety discrimination system is important for protecting breeders’ rights and producers’ profits. Simple sequence repeat (SSR) databases that can be utilized for both next-generation sequencing (SSR-GBS) and polymerase chain reaction–capillary electrophoresis (PCR-CE) would be very useful in variety discrimination. In the present study, SSRs with tri-, tetra- and pentanucleotide repeats were examined in apple, pear and tea. Out of 37 SSRs that showed clear results in PCR-CE, 27 were suitable for SSR-GBS. Among the remaining markers, there was allele dropout for some markers that caused differences between the results of PCR-CE and SSR-GBS. For the selected 27 markers, the alleles detected by SSR-GBS were comparable to those detected by PCR-CE. Furthermore, we developed a computational pipeline for automated genotyping using SSR-GBS by setting a value “α” for each marker, a criterion whether a genotype is homozygous or heterozygous based on allele frequency. The set of 27 markers contains 10, 8 and 9 SSRs for apple, pear and tea, respectively, that are useful for both PCR-CE and SSR-GBS and suitable for automation. The databases help researchers discriminate varieties in various ways depending on sample size, markers and methods.

Metabolic drug survey highlights cancer cell dependencies and vulnerabilities

Interrogation of cellular metabolism with high-throughput screening approaches can unravel contextual biology and identify cancer-specific metabolic vulnerabilities. To systematically study the consequences of distinct metabolic perturbations, we assemble a comprehensive metabolic drug library (CeMM Library of Metabolic Drugs; CLIMET) covering 243 compounds. We, next, characterize it phenotypically in a diverse panel of myeloid leukemia cell lines and primary patient cells. Analysis of the drug response profiles reveals that 77 drugs affect cell viability, with the top effective compounds targeting nucleic acid synthesis, oxidative stress, and the PI3K/mTOR pathway. Clustering of individual drug response profiles stratifies the cell lines into five functional groups, which link to specific molecular and metabolic features. Mechanistic characterization of selective responses to the PI3K inhibitor pictilisib, the fatty acid synthase inhibitor GSK2194069, and the SLC16A1 inhibitor AZD3965, bring forth biomarkers of drug response. Phenotypic screening using CLIMET represents a valuable tool to probe cellular metabolism and identify metabolic dependencies at large.

Helicobacter pylori virulence factors: relationship between genetic variability and phylogeographic origin

Background

Helicobacter pylori is a pathogenic bacteria that colonize the gastrointestinal tract from human stomachs and causes diseases including gastritis, peptic ulcers, gastric lymphoma (MALT), and gastric cancer, with a higher prevalence in developing countries. Its high genetic diversity among strains is caused by a high mutation rate, observing virulence factors (VFs) variations in different geographic lineages. This study aimed to postulate the genetic variability associated with virulence factors present in the Helicobacter pylori strains, to identify the relationship of these genes with their phylogeographic origin.

Methods

The complete genomes of 135 strains available in NCBI, from different population origins, were analyzed using bioinformatics tools, identifying a high rate; as well as reorganization events in 87 virulence factor genes, divided into seven functional groups, to determine changes in position, number of copies, nucleotide identity and size, contrasting them with their geographical lineage and pathogenic phenotype.

Results

Bioinformatics analyses show a high rate of gene annotation errors in VF. Analysis of genetic variability of VFs shown that there is not a direct relationship between the reorganization and geographic lineage. However, regarding the pathogenic phenotype demonstrated in the analysis of many copies, size, and similarity when dividing the strains that possess and not the cag pathogenicity island (cagPAI), having a higher risk of developing gastritis and peptic ulcer was evidenced. Our data has shown that the analysis of the overall genetic variability of all VFs present in each strain of H. pylori is key information in understanding its pathogenic behavior.

A comparison of food sources of nudibranch mollusks at different depths off the Kuril Islands using fatty acid trophic markers

Gastropod molluscs such as nudibranchs are important members of deep-sea benthic ecosystems. However, data on the trophic ecology and feeding specialization of these animals are limited to date. The method of fatty acid trophic markers (FATM) was applied to determine the dietary preferences of nudibranchs off the Kuril Islands. Fatty acid (FA) compositions of Dendronotus sp., Tritonia tetraquetra, and Colga pacifica collected from deep waters were analyzed and compared with those of Aeolidia papillosa and Coryphella verrucosa from the offshore zone. The high level of FATM such as 22:5n-6 and C₂₀ monounsaturated FAs indicated that Dendronotus sp. preys on sea anemones and/or anthoathecates hydroids similarly to that of shallow-water species A. papillosa and C. verrucosa. The high percentage of tetracosapolyenoic acids and the ratio 24:6n-3/24:5n-6 indicated that T. tetraquetra preys on soft corals such as Gersemia and/or Acanella at a depth of 250 m, but soft corals of the family Primnoidae may be the main item in the diet of T. tetraquetra at a depth of 500 m. The high content of Δ 7,13-22:2 and 22:6n-3 shows that C. pacifica can feed on bryozoans. In C. pacifica, 22:5n-6 may be synthesized intrinsically by the mollusks, whereas odd-chain and branched saturated FAs originate from associated bacteria.

Functional constraints on insect immune system components govern their evolutionary trajectories

Roles of constraints in shaping evolutionary outcomes are often considered in the contexts of developmental biology and population genetics, in terms of capacities to generate new variants and how selection limits or promotes consequent phenotypic changes. Comparative genomics also recognizes the role of constraints, in terms of shaping evolution of gene and genome architectures, sequence evolutionary rates, and gene gains or losses, as well as on molecular phenotypes. Characterizing patterns of genomic change where putative functions and interactions of system components are relatively well described offers opportunities to explore whether genes with similar roles exhibit similar evolutionary trajectories. Using insect immunity as our test case system, we hypothesize that characterizing gene evolutionary histories can define distinct dynamics associated with different functional roles. We develop metrics that quantify gene evolutionary histories, employ these to characterize evolutionary features of immune gene repertoires, and explore relationships between gene family evolutionary profiles and their roles in immunity to understand how different constraints may relate to distinct dynamics. We identified three main axes of evolutionary trajectories characterized by gene duplication and synteny, maintenance/stability and sequence conservation, and loss and sequence divergence, highlighting similar and contrasting patterns across these axes amongst subsets of immune genes. Our results suggest that where and how genes participate in immune responses limit the range of possible evolutionary scenarios they exhibit. The test case study system of insect immunity highlights the potential of applying comparative genomics approaches to characterize how functional constraints on different components of biological systems govern their evolutionary trajectories.

A Study on the Clustering of Extra Virgin Olive Oils Extracted from Cultivars Growing in Four Ionian Islands (Greece) by Multivariate Analysis of Their Phenolic Profile, Antioxidant Activity and Genetic Markers

Background: The phenolic fraction of extra virgin olive oil (EVOO) has disease preventive and health-promoting properties which are supported by numerous studies. As such, EVOO is defined as a functional food. The aim of the present study was to characterize the phenolic profile of olive oil from cultivars farmed in the Ionian Islands (Zakynthos, Kefalonia, Lefkada, and Kerkyra) and to investigate the association of phenols to antioxidant activity, which is central to its functionality. Furthermore, the study investigates whether multivariate analyses on the concentration of individual biophenolic compounds and genetic population diversity could classify the olive oil samples based on their geographic origin. Methods: Phenols were determined in 103 samples from different Ionian Island tree populations by ¹H nuclear magnetic resonance (NMR), and sample antioxidant activity was measured by their capacity to reduce the free radical 2,2-diphenyl-1-picrylhydrazyl) (DPPH). Genetic diversity was measured by estimating Nei’s population genetic distance using 15 reproducible bands from random amplified polymorphic DNA (RAPD) genotyping. Results: Principal component analysis (PCA) of the secoiridoid concentrations clustered samples according to cultivar. Clustering based on genetic distances is not concordant with phenolic clustering. A cultivar effect was also demonstrated in the association between the concentration of individual phenols with DPPH reducing activity. Conclusions: Taken together, the study shows that the olive oil phenolic content defines “cultivar-specific phenolic profiles” and that environmental factors other than agronomic conditions contribute more to phenotype variance than genetics.

Observable variations in human sex ratio at birth

The human sex ratio at birth (SRB), defined as the ratio between the number of newborn boys to the total number of newborns, is typically slightly greater than 1/2 (more boys than girls) and tends to vary across different geographical regions and time periods. In this large-scale study, we sought to validate previously-reported associations and test new hypotheses using statistical analysis of two very large datasets incorporating electronic medical records (EMRs). One of the datasets represents over half (∼ 150 million) of the US population for over 8 years (IBM Watson Health MarketScan insurance claims) while another covers the entire Swedish population (∼ 9 million) for over 30 years (the Swedish National Patient Register). After testing more than 100 hypotheses, we showed that neither dataset supported models in which the SRB changed seasonally or in response to variations in ambient temperature. However, increased levels of a diverse array of air and water pollutants, were associated with lower SRBs, including increased levels of industrial and agricultural activity, which served as proxies for water pollution. Moreover, some exogenous factors generally considered to be environmental toxins turned out to induce higher SRBs. Finally, we identified new factors with signals for either higher or lower SRBs. In all cases, the effect sizes were modest but highly statistically significant owing to the large sizes of the two datasets. We suggest that while it was unlikely that the associations have arisen from sex-specific selection mechanisms, they are still useful for the purpose of public health surveillance if they can be corroborated by empirical evidences.

The human sex ratio at birth (SRB), usually slightly greater than 1/2, have been reported to vary in response to a wide array of exogenous factors. In the literature, many such factors have been posited to be associated with higher or lower SRBs, but the studies conducted so far have focused on no more than a few factors at a time and used far smaller datasets, thus prone to generating spurious correlations. We performed a series of statistical tests on 2 large, country-wide health datasets representing the United States and Sweden to investigate associations between putative exogenous factors and the SRB, and were able to validate a set of previously-reported associations while also discovering new signals. We propose to interpret these results simply as public health indicators awaiting further empirical confirmation rather than as implicated in (adaptive) sexual selection mechanisms.

Why Environmental Biomarkers Work: Transcriptome-Proteome Correlations and Modeling of Multistressor Experiments in the Marine Bacterium Trichodesmium

Ocean microbial communities are important contributors to the global biogeochemical reactions that sustain life on Earth. The factors controlling these communities are being increasingly explored using metatranscriptomic and metaproteomic environmental biomarkers. Using published proteomes and transcriptomes from the abundant colony-forming cyanobacterium Trichodesmium (strain IMS101) grown under varying Fe and/or P limitation in low and high CO₂, we observed robust correlations of stress-induced proteins and RNAs (i.e., involved in transport and homeostasis) that yield useful information on the nutrient status under low and/or high CO₂. Conversely, transcriptional and translational correlations of many other central metabolism pathways exhibit broad discordance. A cellular RNA and protein production/degradation model demonstrates how biomolecules with small initial inventories, such as environmentally responsive proteins, achieve large increases in fold-change units as opposed to those with a higher basal expression and inventory such as metabolic systems. Microbial cells, due to their immersion in the environment, tend to show large adaptive responses in both RNA and protein that result in transcript-protein correlations. These observations and model results demonstrate multi-omic coherence for environmental biomarkers and provide the underlying mechanism for those observations, supporting the promise for global application in detecting responses to environmental stimuli in a changing ocean.

Evidence for selection and spatially distinct patterns found in a putative zona pellucida gene in Pacific cod, and implications for management

Genetic differentiation has been observed in marine species even when no obvious barriers to gene flow exist, and understanding such differentiation is essential for effective fisheries management. Highly differentiated outlier loci can provide information on how genetic variation might not only contribute to local adaptation but may also be affected by historical demographic events. A locus which aligned to a predicted zona pellucida sperm-binding protein 3 gene (ZP3) in Atlantic cod (Gadus morhua) was previously identified as the highest outlier based on F ST in a RADseq study of Pacific cod (Gadus macrocephalus) across the West Coast of North America. However, because of the limited length of the RAD sequence and restricted geographic area of sampling, no conclusion on the functional significance of the observed variation was possible. In other marine species, ZP3 is involved in reproductive isolation, local adaptation, and has neofunctionalized as an antifreeze gene, and so it may provide important insights in functional population structure of Pacific cod. Here, we sequenced a 544-bp region of ZP3 in 230 Pacific cod collected from throughout their geographic range. We observed striking patterns of spatial structuring of ZP3 haplotypes, with a sharp break near Kodiak, Alaska, USA where populations within ~200 km of each other are nearly fixed for different haplotypes, contrasting a pattern of isolation by distance at other genetic markers in this region (F ST = 0.003). Phylogenetic analysis of ZP3 haplotypes revealed that the more southern haplotypes appear to be ancestral, with the northern haplotype evolving more recently, potentially in response to a novel selective pressure as Pacific cod recolonized northern latitudes after glaciation. The sharp break in haplotype frequencies suggests strong selective pressures are operating on small spatial scales and illustrates that selection can create high divergence even in marine species with ample opportunities for gene flow.

Surface Roughness of Titanium Orthopedic Implants Alters the Biological Phenotype of Human Mesenchymal Stromal Cells

Metal orthopedic implants are largely biocompatible and generally achieve long-term structural fixation. However, some orthopedic implants may loosen over time even in the absence of infection. In vivo fixation failure is multifactorial, but the fundamental biological defect is cellular dysfunction at the host-implant interface. Strategies to reduce the risk of short- and long-term loosening include surface modifications, implant metal alloy type, and adjuvant substances such as polymethylmethacrylate cement. Surface modifications (e.g., increased surface rugosity) can increase osseointegration and biological ingrowth of orthopedic implants. However, the localized responses of cells to implant surface modifications need to be better characterized. As an in vitro model for investigating cellular responses to metallic orthopedic implants, we cultured mesenchymal stromal/stem cells on clinical-grade titanium disks (Ti6Al4V) that differed in surface roughness as high (porous structured), medium (grit blasted), and low (bead blasted). Topological characterization of clinically relevant titanium (Ti) materials combined with differential mRNA expression analyses (RNA-seq and real-time quantitative polymerase chain reaction) revealed alterations to the biological phenotype of cells cultured on titanium structures that favor early extracellular matrix production and observable responses to oxidative stress and heavy metal stress. These results provide a descriptive model for the interpretation of cellular responses at the interface between native host tissues and three-dimensionally printed modular orthopedic implants, and will guide future studies aimed at increasing the long-term retention of such materials after total joint arthroplasty. Impact statement Using an in vitro model of implant-to-cell interactions by culturing mesenchymal stromal cells (MSCs) on clinically relevant titanium materials of varying topological roughness, we identified mRNA expression patterns consistent with early extracellular matrix (ECM) production and responses to oxidative/heavy metal stress. Implants with high surface roughness may delay the differentiation and ECM formation of MSCs and alter the expression of genes sensitive to reactive oxygen species and protein kinases. In combination with ongoing animal studies, these results will guide future studies aimed at increasing the long-term retention of widely used titanium materials after total joint arthroplasty.

Single-cell atlas of tumor cell evolution in response to therapy in hepatocellular carcinoma and intrahepatic cholangiocarcinoma

BACKGROUND & AIMS

Intratumor molecular heterogeneity is a key feature of tumorigenesis and is linked to treatment failure and patient prognosis. Herein, we aimed to determine what drives tumor cell evolution by performing single-cell transcriptomic analysis.

METHODS

We analyzed 46 hepatocellular carcinoma (HCC) and intrahepatic cholangiocarcinoma (iCCA) biopsies from 37 patients enrolled in interventional studies at the NIH Clinical Center, with 16 biopsies collected before and after treatment from 7 patients. We developed a novel machine learning-based consensus clustering approach to track cellular states of 57,000 malignant and non-malignant cells including tumor cell transcriptome-based functional clonality analysis. We determined tumor cell relationships using RNA velocity and reverse graph embedding. We also studied longitudinal samples from 4 patients to determine tumor cellular state and its evolution. We validated our findings in bulk transcriptomic data from 488 patients with HCC and 277 patients with iCCA.

RESULTS

Using transcriptomic clusters as a surrogate for functional clonality, we observed an increase in tumor cell state heterogeneity which was tightly linked to patient prognosis. Furthermore, increased functional clonality was accompanied by a polarized immune cell landscape which included an increase in pre-exhausted T cells. We found that SPP1 expression was tightly associated with tumor cell evolution and microenvironmental reprogramming. Finally, we developed a user-friendly online interface as a knowledge base for a single-cell atlas of liver cancer.

CONCLUSIONS

Our study offers insight into the collective behavior of tumor cell communities in liver cancer as well as potential drivers of tumor evolution in response to therapy.

LAY SUMMARY

Intratumor molecular heterogeneity is a key feature of tumorigenesis that is linked to treatment failure and patient prognosis. In this study, we present a single-cell atlas of liver tumors from patients treated with immunotherapy and describe intratumoral cell states and their hierarchical relationship. We suggest osteopontin, encoded by the gene SPP1, as a candidate regulator of tumor evolution in response to treatment.

A High-Throughput Short Sequence Typing Scheme for Serratia marcescens Pure Culture and Environmental DNA

Molecular typing methods are used to characterize the relatedness between bacterial isolates involved in infections. These approaches rely mostly on discrete loci or whole-genome sequencing (WGS) analyses of pure cultures. On the other hand, their application to environmental DNA profiling to evaluate epidemiological relatedness among patients and environments has received less attention. We developed a specific, high-throughput short sequence typing (HiSST) method for the opportunistic human pathogen Serratia marcescens. Genes displaying the highest polymorphism were retrieved from the core genome of 60 S. marcescens strains. Bioinformatics analyses showed that use of only three loci (within bssA, gabR, and dhaM) distinguished strains with a high level of efficiency. This HiSST scheme was applied to an epidemiological survey of S. marcescens in a neonatal intensive care unit (NICU). In a first case study, a strain responsible for an outbreak in the NICU was found in a sink drain of this unit, by using HiSST scheme and confirmed by WGS. The HiSST scheme was also applied to environmental DNA extracted from sink-environment samples. Diversity of S. marcescens was modest, with 11, 6, and 4 different sequence types (ST) of gabR, bssA, and dhaM loci among 19 sink drains, respectively. Epidemiological relationships among sinks were inferred on the basis of pairwise comparisons of ST profiles. Further research aimed at relating ST distribution patterns to environmental features encompassing sink location, utilization, and microbial diversity is needed to improve the surveillance and management of opportunistic pathogens. IMPORTANCE Serratia marcescens is an important opportunistic human pathogen, often multidrug resistant and involved in outbreaks of nosocomial infections in neonatal intensive care units. Here, we propose a quick and user-friendly method to select the best typing scheme for nosocomial outbreaks in relating environmental and clinical sources. This method, named high-throughput short sequence typing (HiSST), allows to distinguish strains and to explore the diversity profile of nonculturable S. marcescens. The application of HiSST profile analysis for environmental DNA offers new possibilities to track opportunistic pathogens, identify their origin, and relate their distribution pattern with environmental features encompassing sink location, utilization, and microbial diversity. Adaptation of the method to other opportunistic pathogens is expected to improve knowledge regarding their ecology, which is of significant interest for epidemiological risk assessment and elaborate outbreak mitigation strategies.

TargetTox: A Feature Selection Pipeline for Identifying Predictive Targets Associated with Drug Toxicity

In silico assessment of drug toxicity is becoming a critical step in drug development. Conventional ligand-based models are limited by low accuracy and lack of interpretability. Further, they often fail to explain cellular mechanisms underlying structure-toxicity associations. We addressed these limitations by incorporating target profile as an intermediate connecting structure to toxicity. To accommodate for high-dimensional feature space, we developed a pipeline named TargetTox that can identity a subset of predictive features. We implemented TargetTox to study 569 targets and 815 adverse events. The features identified by TargetTox comprise less than 10% of the original feature space; nevertheless, they accurately predicted binding outcomes for 377 targets and toxicity outcomes for 36 adverse events. We demonstrated that predictive targets tend to be differentially expressed in the tissue of toxicity. We also rediscovered key cellular functions associated with cardiotoxicity from the predictive targets, as well as markers of skin and liver diseases. Furthermore, we found evidence supporting diagnostic and therapeutic applications of some predictive targets in hepatotoxicity and nephrotoxicity. Our findings highlighted the critical role of predictive targets in cellular mechanisms leading to toxicity. In general, our study improved the interpretability of toxicity prediction without sacrificing accuracy. Our novel pipeline may benefit future studies of high-dimensional data sets.