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

Network-based analysis of heterogeneous patient-matched brain and extracranial melanoma metastasis pairs reveals three homogeneous subgroups

Melanoma, the deadliest form of skin cancer, can metastasize to different organs. Molecular differences between brain and extracranial melanoma metastases are poorly understood. Here, promoter methylation and gene expression of 11 heterogeneous patient-matched pairs of brain and extracranial metastases were analyzed using melanoma-specific gene regulatory networks learned from public transcriptome and methylome data followed by network-based impact propagation of patient-specific alterations. This innovative data analysis strategy allowed to predict potential impacts of patient-specific driver candidate genes on other genes and pathways. The patient-matched metastasis pairs clustered into three robust subgroups with specific downstream targets with known roles in cancer, including melanoma (SG1: RBM38, BCL11B, SG2: GATA3, FES, SG3: SLAMF6, PYCARD). Patient subgroups and ranking of target gene candidates were confirmed in a validation cohort. Summarizing, computational network-based impact analyses of heterogeneous metastasis pairs predicted individual regulatory differences in melanoma brain metastases, cumulating into three consistent subgroups with specific downstream target genes.

Single-cell RNA-sequencing provides new insights into the cell-specific expression patterns and transcriptional regulation of photosynthetic genes in bermudagrass leaf blades

As an important warm-season turfgrass species, bermudagrass (Cynodon dactylon L.) flourishes in warm areas around the world due to the existence of the C4 photosynthetic pathway. However, how C4 photosynthesis operates in bermudagrass leaves is still poorly understood. In this study, we performed single-cell RNA-sequencing on 5296 cells from bermudagrass leaf blades. Eight cell clusters corresponding to mesophyll, bundle sheath, epidermis and vascular bundle cells were successfully identified using known cell marker genes. Expression profiling indicated that genes encoding NADP-dependent malic enzymes (NADP-MEs) were highly expressed in bundle sheath cells, whereas NAD-ME genes were weakly expressed in all cell types, suggesting C4 photosynthesis of bermudagrass leaf blades might be NADP-ME type rather than NAD-ME type. The results also indicated that starch synthesis-related genes showed preferential expression in bundle sheath cells, whereas starch degradation-related genes were highly expressed in mesophyll cells, which agrees with the observed accumulation of starch-filled chloroplasts in bundle sheath cells. Gene co-expression analysis further revealed that different families of transcription factors were co-expressed with multiple C4 photosynthesis-related genes, suggesting a complex transcription regulatory network of C4 photosynthesis might exist in bermudagrass leaf blades. These findings collectively provided new insights into the cell-specific expression patterns and transcriptional regulation of photosynthetic genes in bermudagrass.

Single-cell resolution characterization of myeloid-derived cell states with implication in cancer outcome

Tumor-associated myeloid-derived cells (MDCs) significantly impact cancer prognosis and treatment responses due to their remarkable plasticity and tumorigenic behaviors. Here, we integrate single-cell RNA-sequencing data from different cancer types, identifying 29 MDC subpopulations within the tumor microenvironment. Our analysis reveals abnormally expanded MDC subpopulations across various tumors and distinguishes cell states that have often been grouped together, such as TREM2+ and FOLR2+ subpopulations. Using deconvolution approaches, we identify five subpopulations as independent prognostic markers, including states co-expressing TREM2 and PD-1, and FOLR2 and PDL-2. Additionally, TREM2 alone does not reliably predict cancer prognosis, as other TREM2+ macrophages show varied associations with prognosis depending on local cues. Validation in independent cohorts confirms that FOLR2-expressing macrophages correlate with poor clinical outcomes in ovarian and triple-negative breast cancers. This comprehensive MDC atlas offers valuable insights and a foundation for futher analyses, advancing strategies for treating solid cancers.

Genetic and epigenetic alterations in precursor lesions of endometrial endometrioid carcinoma

The hyperplasia-carcinoma sequence is a stepwise tumourigenic programme towards endometrial cancer in which normal endometrial epithelium becomes neoplastic through non-atypical endometrial hyperplasia (NAEH) and atypical endometrial hyperplasia (AEH), under the influence of unopposed oestrogen. NAEH and AEH are known to exhibit polyclonal and monoclonal cell growth, respectively; yet, aside from focal PTEN protein loss, the genetic and epigenetic alterations that occur during the cellular transition remain largely unknown. We sought to explore the potential molecular mechanisms that promote the NAEH-AEH transition and identify molecular markers that could help to differentiate between these two states. We conducted target-panel sequencing on the coding exons of 596 genes, including 96 endometrial cancer driver genes, and DNA methylome microarrays for 48 NAEH and 44 AEH lesions that were separately collected via macro- or micro-dissection from the endometrial tissues of 30 cases. Sequencing analyses revealed acquisition of the PTEN mutation and the clonal expansion of tumour cells in AEH samples. Further, across the transition, alterations to the DNA methylome were characterised by hypermethylation of promoter/enhancer regions and CpG islands, as well as hypo- and hyper-methylation of DNA-binding regions for transcription factors relevant to endometrial cell differentiation and/or tumourigenesis, including FOXA2, SOX17, and HAND2. The identified DNA methylation signature distinguishing NAEH and AEH lesions was reproducible in a validation cohort with modest discriminative capability. These findings not only support the concept that the transition from NAEH to AEH is an essential step within neoplastic cell transformation of endometrial epithelium but also provide deep insight into the molecular mechanism of the tumourigenic programme. © 2024 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Self-assembling viral histones are evolutionary intermediates between archaeal and eukaryotic nucleosomes

Nucleosomes are DNA-protein complexes composed of histone proteins that form the basis of eukaryotic chromatin. The nucleosome was a key innovation during eukaryotic evolution, but its origin from histone homologues in Archaea remains unclear. Viral histone repeats, consisting of multiple histone paralogues within a single protein, may reflect an intermediate state. Here we examine the diversity of histones encoded by Nucleocytoviricota viruses. We identified 258 histones from 168 viral metagenomes with variable domain configurations including histone singlets, doublets, triplets and quadruplets, the latter comprising the four core histones arranged in series. Viral histone repeats branch phylogenetically between Archaea and eukaryotes and display intermediate functions in Escherichia coli, self-assembling into eukaryotic-like nucleosomes that stack into archaeal-like oligomers capable of impacting genomic activity and condensing DNA. Histone linkage also facilitates nucleosome formation, promoting eukaryotic histone assembly in E. coli. These data support the hypothesis that viral histone repeats originated in stem-eukaryotes and that nucleosome evolution proceeded through histone repeat intermediates.

Evolutionary and biomedical implications of sex differences in the primate brain transcriptome

Humans exhibit sex differences in the prevalence of many neurodevelopmental disorders and neurodegenerative diseases. Here, we generated one of the largest multi-brain-region bulk transcriptional datasets for the rhesus macaque and characterized sex-biased gene expression patterns to investigate the translatability of this species for sex-biased neurological conditions. We identify patterns similar to those in humans, which are associated with overlapping regulatory mechanisms, biological processes, and genes implicated in sex-biased human disorders, including autism. We also show that sex-biased genes exhibit greater genetic variance for expression and more tissue-specific expression patterns, which may facilitate rapid evolution of sex-biased genes. Our findings provide insights into the biological mechanisms underlying sex-biased disease and support the rhesus macaque model for the translational study of these conditions.

Spatial distribution of pesticide use based on crop rotation data in La Plata River basin: a case study from an agricultural region of Uruguay

The intensive global use of pesticides presents an escalating threat to human health, ecosystems, and water quality. To develop national and local environmental management strategies for mitigating pollution caused by pesticides, it is essential to understand the quantities, timing, and location of their application. This study aims to estimate the spatial distribution of pesticide use in an agricultural region of La Plata River basin in Uruguay. Estimates of pesticide use were made by surveying doses applied to each crop. This information was spatialized through identifying agricultural rotations using remote sensing techniques. The study identified the 60 major agricultural rotations in the region and mapped the use and application amount of the nine most significant active ingredients (glyphosate, 2,4-dichlorophenoxyacetic acid, flumioxazin, S-metolachlor, clethodim, flumetsulam, triflumuron, chlorantraniliprole, and fipronil). The results reveal that glyphosate is the most extensively used pesticide (53.5% of the area) and highest amount of use (> 1.44 kg/ha). Moreover, in 19% of the area, at least seven active ingredients are applied in crop rotations. This study marks the initial step in identifying rotations and estimating pesticide applications with high spatial resolution at a regional scale in agricultural regions of La Plata River basin. The results improve the understanding of pesticide spatial distribution based on data obtained from agronomists, technicians, and producers and provide a replicable methodological approach for other geographic and productive contexts. Generating baseline information is key to environmental management and decision making, towards the design of more robust monitoring systems and human exposure assessment.

Chloroplast Genomes Evolution and Phylogenetic Relationships of Caragana species

Caragana sensu lato (s.l.) includes approximately 100 species that are mainly distributed in arid and semi-arid regions. Caragana species are ecologically valuable for their roles in windbreaking and sand fixation. However, the taxonomy and phylogenetic relationships of the genus Caragana are still unclear. In this study, we sequenced and assembled the chloroplast genomes of representative species of Caragana and reconstructed robust phylogenetic relationships at the section level. The Caragana chloroplast genome has lost the inverted repeat region and wascategorized in the inverted repeat loss clade (IRLC). The chloroplast genomes of the eight species ranged from 128,458 bp to 135,401 bp and contained 110 unique genes. All the Caragana chloroplast genomes have a highly conserved structure and gene order. The number of long repeats and simple sequence repeats (SSRs) showed significant variation among the eight species, indicating heterogeneous evolution in Caragana. Selective pressure analysis of the genes revealed that most of the protein-coding genes evolved under purifying selection. The phylogenetic analyses indicated that each section forms a clade, except the section Spinosae, which was divided into two clades. This study elucidated the evolution of the chloroplast genome within the widely distributed genus Caragana. The detailed information obtained from this study can serve as a valuable resource for understanding the molecular dynamics and phylogenetic relationships within Caragana.

Functional architecture of the forebrain cholinergic system in rodents

The basal forebrain cholinergic system (BFCS) participates in functions that are global across the brain, such as sleep-wake cycles, but also participates in capacities that are more behaviorally and anatomically specific, including sensory perception. To better understand the underlying organization principles of the BFCS, more and higher quality anatomical data and analysis is needed. Here, we created a "virtual Basal Forebrain", combining data from numerous rats with cortical retrograde tracer injections into a common 3D reference coordinate space and developed a "spatial density correlation" methodology to analyze patterns in BFCS cortical projection targets, revealing that the BFCS is organized into three principal networks: somatosensory-motor, auditory, and visual. Within each network, clusters of cholinergic cells with increasing complexity innervate cortical targets. These networks represent hierarchically organized building blocks that may enable the BFCS to coordinate spatially selective signaling, including parallel modulation of multiple functionally interconnected yet diverse groups of cortical areas.

Transcriptional immune suppression and up-regulation of double-stranded DNA damage and repair repertoires in ecDNA-containing tumors

Extrachromosomal DNA is a common cause of oncogene amplification in cancer. The non-chromosomal inheritance of ecDNA enables tumors to rapidly evolve, contributing to treatment resistance and poor outcome for patients. The transcriptional context in which ecDNAs arise and progress, including chromosomally-driven transcription, is incompletely understood. We examined gene expression patterns of 870 tumors of varied histological types, to identify transcriptional correlates of ecDNA. Here, we show that ecDNA-containing tumors impact four major biological processes. Specifically, ecDNA-containing tumors up-regulate DNA damage and repair, cell cycle control, and mitotic processes, but down-regulate global immune regulation pathways. Taken together, these results suggest profound alterations in gene regulation in ecDNA-containing tumors, shedding light on molecular processes that give rise to their development and progression.

Clustering plasma concentration-time curves: applications of unsupervised learning in pharmacogenomics

Pharmaceutical researchers are continually searching for techniques to improve both drug development processes and patient outcomes. An area of recent interest is the potential for machine learning (ML) applications within pharmacology. One such application not yet given close study is the unsupervised clustering of plasma concentration-time curves, hereafter, pharmacokinetic (PK) curves. In this paper, we present our findings on how to cluster PK curves by their similarity. Specifically, we find clustering to be effective at identifying similar-shaped PK curves and informative for understanding patterns within each cluster of PK curves. Because PK curves are time series data objects, our approach utilizes the extensive body of research related to the clustering of time series data as a starting point. As such, we examine many dissimilarity measures between time series data objects to find those most suitable for PK curves. We identify Euclidean distance as generally most appropriate for clustering PK curves, and we further show that dynamic time warping, Fréchet, and structure-based measures of dissimilarity like correlation may produce unexpected results. As an illustration, we apply these methods in a case study with 250 PK curves used in a previous pharmacogenomic study. Our case study finds that an unsupervised ML clustering with Euclidean distance, without any subject genetic information, is able to independently validate the same conclusions as the reference pharmacogenomic results. To our knowledge, this is the first such demonstration. Further, the case study demonstrates how the clustering of PK curves may generate insights that could be difficult to perceive solely with population level summary statistics of PK metrics.

Effects of individual traits vs. trait syndromes on assemblages of various herbivore guilds associated with central European Salix

Plants employ diverse anti-herbivore defences that can covary to form syndromes consisting of multiple traits. Such syndromes are hypothesized to impact herbivores more than individual defences. We studied 16 species of lowland willows occurring in central Europe and explored if their chemical and physical traits form detectable syndromes. We tested for phylogenetic trends in the syndromes and explored whether three herbivore guilds (i.e., generalist leaf-chewers, specialist leaf-chewers, and gallers) are affected more by the detected syndromes or individual traits. The recovered syndromes showed low phylogenetic signal and were mainly defined by investment in concentration, richness, or uniqueness of structurally related phenolic metabolites. Resource acquisition traits or inducible volatile organic compounds exhibited a limited correlation with the syndromes. Individual traits composing the syndromes showed various correlations to the assemblages of herbivores from the three studied guilds. In turn, we found some support for the hypothesis that defence syndromes are composed of traits that provide defence against various herbivores. However, individual traits rather than trait syndromes explained more variation for all studied herbivore assemblages. The detected negative correlations between various phenolics suggest that investment trade-offs may occur primarily among plant metabolites with shared metabolic pathways that may compete for their precursors. Moreover, several traits characterizing the recovered syndromes play additional roles in willows other than defence from herbivory. Taken together, our findings suggest that the detected syndromes did not solely evolve as an anti-herbivore defence.

biotextgraph: graphical summarization of functional similarities from textual information

SUMMARY

Functional interpretation of biological entities such as differentially expressed genes is one of the fundamental analyses in bioinformatics. The task can be addressed by using biological pathway databases with enrichment analysis (EA). However, textual description of biological entities in public databases is less explored and integrated in existing tools and it has a potential to reveal new mechanisms. Here, we present a new R package biotextgraph for graphical summarization of omics' textual description data which enables assessment of functional similarities of the lists of biological entities. We illustrate application examples of annotating gene identifiers in addition to EA. The results suggest that the visualization based on words and inspection of biological entities with text can reveal a set of biologically meaningful terms that could not be obtained by using biological pathway databases alone. The results suggest the usefulness of the package in the routine analysis of omics-related data. The package also offers a web-based application for convenient querying.

AVAILABILITY AND IMPLEMENTATION

The package, documentation, and web server are available at: https://github.com/noriakis/biotextgraph.

The Yin and Yang of pathogens and probiotics: interplay between Salmonella enterica sv. Typhimurium and Bifidobacterium infantis during co-infection

Probiotic bacteria have been proposed as an alternative to antibiotics for the control of antimicrobial resistant enteric pathogens. The mechanistic details of this approach remain unclear, in part because pathogen reduction appears to be both strain and ecology dependent. Here we tested the ability of five probiotic strains, including some from common probiotic genera Lactobacillus and Bifidobacterium, to reduce binding of Salmonella enterica sv. Typhimurium to epithelial cells in vitro. Bifidobacterium longum subsp. infantis emerged as a promising strain; however, S. Typhimurium infection outcome in epithelial cells was dependent on inoculation order, with B. infantis unable to rescue host cells from preceding or concurrent infection. We further investigated the complex mechanisms underlying this interaction between B. infantis, S. Typhimurium, and epithelial cells using a multi-omics approach that included gene expression and altered metabolism via metabolomics. Incubation with B. infantis repressed apoptotic pathways and induced anti-inflammatory cascades in epithelial cells. In contrast, co-incubation with B. infantis increased in S. Typhimurium the expression of virulence factors, induced anaerobic metabolism, and repressed components of arginine metabolism as well as altering the metabolic profile. Concurrent application of the probiotic and pathogen notably generated metabolic profiles more similar to that of the probiotic alone than to the pathogen, indicating a central role for metabolism in modulating probiotic-pathogen-host interactions. Together these data imply crosstalk via small molecules between the epithelial cells, pathogen and probiotic that consistently demonstrated unique molecular mechanisms specific probiotic/pathogen the individual associations.

Personalized identification and characterization of genome-wide gene expression differences between patient-matched intracranial and extracranial melanoma metastasis pairs

Melanoma is the most serious type of skin cancer that frequently spreads to other organs of the human body. Especially melanoma metastases to the brain (intracranial metastases) are hard to treat and a major cause of death of melanoma patients. Little is known about molecular alterations and altered mechanisms that distinguish intra- from extracranial melanoma metastases. So far, almost all existing studies compared intracranial metastases from one set of patients to extracranial metastases of an another set of melanoma patients. This neglects the important facts that each melanoma is highly individual and that intra- and extracranial melanoma metastases from the same patient are more similar to each other than to melanoma metastases from other patients in the same organ. To overcome this, we compared the gene expression profiles of 16 intracranial metastases to their corresponding 21 patient-matched extracranial metastases in a personalized way using a three-state Hidden Markov Model (HMM) to identify altered genes for each individual metastasis pair. This enabled three major findings by considering the predicted gene expression alterations across all patients: (i) most frequently altered pathways include cytokine-receptor interaction, calcium signaling, ECM-receptor interaction, cAMP signaling, Jak-STAT and PI3K/Akt signaling, (ii) immune-relevant signaling pathway genes were downregulated in intracranial metastases, and (iii) intracranial metastases were associated with a brain-like phenotype gene expression program. Further, the integration of all differentially expressed genes across the patient-matched melanoma metastasis pairs led to a set of 103 genes that were consistently down- or up-regulated in at least 11 of the 16 of the patients. This set of genes contained many genes involved in the regulation of immune responses, cell growth, cellular signaling and transport processes. An analysis of these genes in the TCGA melanoma cohort showed that the expression behavior of 11 genes was significantly associated with survival. Moreover, a comparison of the 103 genes to three closely related melanoma metastasis studies revealed a core set of eight genes that were consistently down- or upregulated in intra- compared to extracranial metastases in at least two of the three related studies (down: CILP, DPT, FGF7, LAMP3, MEOX2, TMEM119; up: GLDN, PMP2) including FGF7 that was also significantly associated with survival. Our findings contribute to a better characterization of genes and pathways that distinguish intra- from extracranial melanoma metastasis and provide important hints for future experimental studies to identify potential targets for new therapeutic approaches.

Gradients of bacteria in the oceanic water column reveal finely-resolved vertical distributions

Bacterial communities directly influence ecological processes in the ocean, and depth has a major influence due to the changeover in primary energy sources between the sunlit photic zone and dark ocean. Here, we examine the abundance and diversity of bacteria in Monterey Bay depth profiles collected from the surface to just above the sediments (e.g., 2000 m). Bacterial abundance in these Pacific Ocean samples decreased by >1 order of magnitude, from 1.22 ±0.69 ×106 cells ml-1 in the variable photic zone to 1.44 ± 0.25 ×105 and 6.71 ± 1.23 ×104 cells ml-1 in the mesopelagic and bathypelagic, respectively. V1-V2 16S rRNA gene profiling showed diversity increased sharply between the photic and mesopelagic zones. Weighted Gene Correlation Network Analysis clustered co-occurring bacterial amplicon sequence variants (ASVs) into seven subnetwork modules, of which five strongly correlated with depth-related factors. Within surface-associated modules there was a clear distinction between a 'copiotrophic' module, correlating with chlorophyll and dominated by e.g., Flavobacteriales and Rhodobacteraceae, and an 'oligotrophic' module dominated by diverse Oceanospirillales (such as uncultured JL-ETNP-Y6, SAR86) and Pelagibacterales. Phylogenetic reconstructions of Pelagibacterales and SAR324 using full-length 16S rRNA gene data revealed several additional subclades, expanding known microdiversity within these abundant lineages, including new Pelagibacterales subclades Ia.B, Id, and IIc, which comprised 4-10% of amplicons depending on the subclade and depth zone. SAR324 and Oceanospirillales dominated in the mesopelagic, with SAR324 clade II exhibiting its highest relative abundances (17±4%) in the lower mesopelagic (300-750 m). The two newly-identified SAR324 clades showed highest relative abundances in the photic zone (clade III), while clade IV was extremely low in relative abundance, but present across dark ocean depths. Hierarchical clustering placed microbial communities from 900 m samples with those from the bathypelagic, where Marinimicrobia was distinctively relatively abundant. The patterns resolved herein, through high resolution and statistical replication, establish baselines for marine bacterial abundance and taxonomic distributions across the Monterey Bay water column, against which future change can be assessed.

The geographic extent of bird populations affected by renewable-energy development

Bird populations are declining globally. Wind and solar energy can reduce emissions of fossil fuels that drive anthropogenic climate change, yet renewable-energy production represents a potential threat to bird species. Surveys to assess potential effects at renewable-energy facilities are exclusively local, and the geographic extent encompassed by birds killed at these facilities is largely unknown, which creates challenges for minimizing and mitigating the population-level and cumulative effects of these fatalities. We performed geospatial analyses of stable hydrogen isotope data obtained from feathers of 871 individuals of 24 bird species found dead at solar- and wind-energy facilities in California (USA). Most species had individuals with a mix of origins, ranging from 23% to 98% nonlocal. Mean minimum distances to areas of likely origin for nonlocal individuals were as close as 97 to >1250 km, and these minimum distances were larger for species found at solar-energy facilities in deserts than at wind-energy facilities in grasslands (Cohen's d = 6.5). Fatalities were drawn from an estimated 30-100% of species' desingated ranges, and this percentage was significantly smaller for species with large ranges found at wind facilities (Pearson's r = -0.67). Temporal patterns in the geographic origin of fatalities suggested that migratory movements and nonmigratory movements, such as dispersal and nomadism, influence exposure to fatality risk for these birds. Our results illustrate the power of using stable isotope data to assess the geographic extent of renewable-energy fatalities on birds. As the buildout of renewable-energy facilities continues, accurate assessment of the geographic footprint of wildlife fatalities can be used to inform compensatory mitigation for their population-level and cumulative effects.

Consistency of data reporting in fluid responsiveness studies in the critically ill setting: the CODEFIRE consensus from the Cardiovascular Dynamic section of the European Society of Intensive Care Medicine

PURPOSE

To provide consensus recommendations regarding hemodynamic data reporting in studies investigating fluid responsiveness and fluid challenge (FC) use in the intensive care unit (ICU).

METHODS

The Executive Committee of the European Society of Intensive Care Medicine (ESICM) commissioned and supervised the project. A panel of 18 international experts and a methodologist identified main domains and items from a systematic literature, plus 2 ancillary domains. A three-step Delphi process based on an iterative approach was used to obtain the final consensus. In the Delphi 1 and 2, the items were selected with strong (≥ 80% of votes) or week agreement (70-80% of votes), while the Delphi 3 generated recommended (≥ 90% of votes) or suggested (80-90% of votes) items (RI and SI, respectively).

RESULTS

We identified 5 main domains initially including 117 items and the consensus finally resulted in 52 recommendations or suggestions: 18 RIs and 2 SIs statements were obtained for the domain "ICU admission", 11 RIs and 1 SI for the domain "mechanical ventilation", 5 RIs for the domain "reason for giving a FC", 8 RIs for the domain pre- and post-FC "hemodynamic data", and 7 RIs for the domain "pre-FC infused drugs". We had no consensus on the use of echocardiography, strong agreement regarding the volume (4 ml/kg) and the reference variable (cardiac output), while weak on administration rate (within 10 min) of FC in this setting.

CONCLUSION

This consensus found 5 main domains and provided 52 recommendations for data reporting in studies investigating fluid responsiveness in ICU patients.

A lineage-resolved cartography of microRNA promoter activity in C. elegans empowers multidimensional developmental analysis

Elucidating the expression of microRNAs in developing single cells is critical for functional discovery. Here, we construct scCAMERA (single-cell cartography of microRNA expression based on reporter assay), utilizing promoter-driven fluorescent reporters in conjunction with imaging and lineage tracing. The cartography delineates the transcriptional activity of 54 conserved microRNAs in lineage-resolved single cells throughout C. elegans embryogenesis. The combinatorial expression of microRNAs partitions cells into fine clusters reflecting their function and anatomy. Notably, the expression of individual microRNAs exhibits high cell specificity and divergence among family members. Guided by cellular expression patterns, we identify developmental functions of specific microRNAs, including miR-1 in pharynx development and physiology, miR-232 in excretory canal morphogenesis by repressing NHR-25/NR5A, and a functional synergy between miR-232 and miR-234 in canal development, demonstrating the broad utility of scCAMERA. Furthermore, integrative analysis reveals that tissue-specific fate determinants activate microRNAs to repress protein production from leaky transcripts associated with alternative, especially neuronal, fates, thereby enhancing the fidelity of developmental fate differentiation. Collectively, our study offers rich opportunities for multidimensional expression-informed analysis of microRNA biology in metazoans.

Effect of a probiotic and an antibiotic on the mobilome of the porcine microbiota

Introduction: To consider the growing health issues caused by antibiotic resistance from a "one health" perspective, the contribution of meat production needs to be addressed. While antibiotic resistance is naturally present in microbial communities, the treatment of farm animals with antibiotics causes an increase in antibiotic resistance genes (ARG) in the gut microbiome. Pigs are among the most prevalent animals in agriculture; therefore, reducing the prevalence of antibiotic-resistant bacteria in the pig gut microbiome could reduce the spread of antibiotic resistance. Probiotics are often studied as a way to modulate the microbiome and are, therefore, an interesting way to potentially decrease antibiotic resistance. Methods: To assess the efficacy of a probiotic to reduce the prevalence of ARGs in the pig microbiome, six pigs received either treatment with antibiotics (tylvalosin), probiotics (Pediococcus acidilactici MA18/5M; Biopower® PA), or a combination of both. Their faeces and ileal digesta were collected and DNA was extracted for whole genome shotgun sequencing. The reads were compared with taxonomy and ARG databases to identify the taxa and resistance genes in the samples. Results: The results showed that the ARG profiles in the faeces of the antibiotic and combination treatments were similar, and both were different from the profiles of the probiotic treatment (p < 0.05). The effects of the treatments were different in the digesta and faeces. Many macrolide resistance genes were detected in a higher proportion in the microbiome of the pigs treated with antibiotics or the combination of probiotics and antibiotics. Resistance-carrying conjugative plasmids and horizontal transfer genes were also amplified in faeces samples for the antibiotic and combined treatments. There was no effect of treatment on the short chain fatty acid content in the digesta or the faeces. Conclusion: There is no positive effect of adding probiotics to an antibiotic treatment when these treatments are administered simultaneously.

Genomic evolution shapes prostate cancer disease type

The development of cancer is an evolutionary process involving the sequential acquisition of genetic alterations that disrupt normal biological processes, enabling tumor cells to rapidly proliferate and eventually invade and metastasize to other tissues. We investigated the genomic evolution of prostate cancer through the application of three separate classification methods, each designed to investigate a different aspect of tumor evolution. Integrating the results revealed the existence of two distinct types of prostate cancer that arise from divergent evolutionary trajectories, designated as the Canonical and Alternative evolutionary disease types. We therefore propose the evotype model for prostate cancer evolution wherein Alternative-evotype tumors diverge from those of the Canonical-evotype through the stochastic accumulation of genetic alterations associated with disruptions to androgen receptor DNA binding. Our model unifies many previous molecular observations, providing a powerful new framework to investigate prostate cancer disease progression.

Cellular and clinical impact of protein phosphatase enzyme epigenetic silencing in multiple cancer tissues

BACKGROUND

Protein Phosphatase Enzymes (PPE) and protein kinases simultaneously control phosphorylation mechanisms that tightly regulate intracellular signalling pathways and stimulate cellular responses. In human malignancies, PPE and protein kinases are frequently mutated resulting in uncontrolled kinase activity and PPE suppression, leading to cell proliferation, migration and resistance to anti-cancer therapies. Cancer associated DNA hypermethylation at PPE promoters gives rise to transcriptional silencing (epimutations) and is a hallmark of cancer. Despite recent advances in sequencing technologies, data availability and computational capabilities, only a fraction of PPE have been reported as transcriptionally inactive as a consequence of epimutations.

METHODS

In this study, we examined promoter-associated DNA methylation profiles in Protein Phosphatase Enzymes and their Interacting Proteins (PPEIP) in a cohort of 705 cancer patients in five tissues (Large intestine, Oesophagus, Lung, Pancreas and Stomach) in three cell models (primary tumours, cancer cell lines and 3D embedded cancer cell cultures). As a subset of PPEIP are known tumour suppressor genes, we analysed the impact of PPEIP promoter hypermethylation marks on gene expression, cellular networks and in a clinical setting.

RESULTS

Here, we report epimutations in PPEIP are a frequent occurrence in the cancer genome and manifest independent of transcriptional activity. We observed that different tumours have varying susceptibility to epimutations and identify specific cellular signalling networks that are primarily affected by epimutations. Additionally, RNA-seq analysis showed the negative impact of epimutations on most (not all) Protein Tyrosine Phosphatase transcription. Finally, we detected novel clinical biomarkers that inform on patient mortality and anti-cancer treatment sensitivity.

CONCLUSIONS

We propose that DNA hypermethylation marks at PPEIP frequently contribute to the pathogenesis of malignancies and within the precision medicine space, hold promise as biomarkers to inform on clinical features such as patient survival and therapeutic response.

Transposition of HOPPLA in siRNA-deficient plants suggests a limited effect of the environment on retrotransposon mobility in Brachypodium distachyon

Long terminal repeat retrotransposons (LTR-RTs) are powerful mutagens regarded as a major source of genetic novelty and important drivers of evolution. Yet, the uncontrolled and potentially selfish proliferation of LTR-RTs can lead to deleterious mutations and genome instability, with large fitness costs for their host. While population genomics data suggest that an ongoing LTR-RT mobility is common in many species, the understanding of their dual role in evolution is limited. Here, we harness the genetic diversity of 320 sequenced natural accessions of the Mediterranean grass Brachypodium distachyon to characterize how genetic and environmental factors influence plant LTR-RT dynamics in the wild. When combining a coverage-based approach to estimate global LTR-RT copy number variations with mobilome-sequencing of nine accessions exposed to eight different stresses, we find little evidence for a major role of environmental factors in LTR-RT accumulations in B. distachyon natural accessions. Instead, we show that loss of RNA polymerase IV (Pol IV), which mediates RNA-directed DNA methylation in plants, results in high transcriptional and transpositional activities of RLC_BdisC024 (HOPPLA) LTR-RT family elements, and that these effects are not stress-specific. This work supports findings indicating an ongoing mobility in B. distachyon and reveals that host RNA-directed DNA methylation rather than environmental factors controls their mobility in this wild grass model.

Microbiomes of Thalassia testudinum throughout the Atlantic Ocean, Caribbean Sea, and Gulf of Mexico are influenced by site and region while maintaining a core microbiome

Plant microbiomes are known to serve several important functions for their host, and it is therefore important to understand their composition as well as the factors that may influence these microbial communities. The microbiome of Thalassia testudinum has only recently been explored, and studies to-date have primarily focused on characterizing the microbiome of plants in a single region. Here, we present the first characterization of the composition of the microbial communities of T. testudinum across a wide geographical range spanning three distinct regions with varying physicochemical conditions. We collected samples of leaves, roots, sediment, and water from six sites throughout the Atlantic Ocean, Caribbean Sea, and the Gulf of Mexico. We then analyzed these samples using 16S rRNA amplicon sequencing. We found that site and region can influence the microbial communities of T. testudinum, while maintaining a plant-associated core microbiome. A comprehensive comparison of available microbial community data from T. testudinum studies determined a core microbiome composed of 14 ASVs that consisted mostly of the family Rhodobacteraceae. The most abundant genera in the microbial communities included organisms with possible plant-beneficial functions, like plant-growth promoting taxa, disease suppressing taxa, and nitrogen fixers.

The Mediterranean mussel Mytilus galloprovincialis: a novel model for developmental studies in mollusks

A model organism in developmental biology is defined by its experimental amenability and by resources created for the model system by the scientific community. For the most powerful invertebrate models, the combination of both has already yielded a thorough understanding of developmental processes. However, the number of developmental model systems is still limited, and their phylogenetic distribution heavily biased. Members of one of the largest animal lineages, the Spiralia, for example, have long been neglected. In order to remedy this shortcoming, we have produced a detailed developmental transcriptome for the bivalve mollusk Mytilus galloprovincialis, and have expanded the list of experimental protocols available for this species. Our high-quality transcriptome allowed us to identify transcriptomic signatures of developmental progression and to perform a first comparison with another bivalve mollusk: the Pacific oyster Crassostrea gigas. To allow co-labelling studies, we optimized and combined protocols for immunohistochemistry and hybridization chain reaction to create high-resolution co-expression maps of developmental genes. The resources and protocols described here represent an enormous boost for the establishment of Mytilus galloprovincialis as an alternative model system in developmental biology.

Galactinol Regulates JA Biosynthesis to Enhance Tomato Cold Tolerance

Low temperatures can inhibit plant growth and development and reduce fruit yield. This study demonstrated that the expression of AnGolS1 from Ammopiptanthus nanus (A. nanus) encoding a galactinol synthase enhanced tomato cold tolerance. In AnGolS1-overexpressing plants, the jasmonic acid (JA) biosynthesis substrates 13-hydroperoxylinolenicacid and 12,13-epoxylinolenicacid were significantly accumulated, and the expression levels of the ethylene response factor (SlERF4-7) and serine protease inhibitor (SlSPI5) were increased. We speculated that there may be correlations among galactinol, ethylene signaling, the protease inhibitor, protease, and JA levels. The expression levels of SlERF4-7 and SlSPI5 as well as the JA content were significantly increased under exogenous galactinol treatment. Additionally, the expression of SlSPI5 was reduced in SlERF4-7-silenced plants, and SlERF4-7 was confirmed to bind to the dehydration-responsive element (DRE) of the SlSPI5 promoter. These results suggest that SlSPI5 is a target gene of the SlERF4-7 transcription factor. In addition, SlSPI5 interacted with cysteine protease (SlCPase), while SlCPase interacted with lipoxygenase (SlLOX5) and allene oxide synthase (SlAOS2). When SlCPase was silenced, JA levels increased and plant cold tolerance was enhanced. Therefore, galactinol regulates JA biosynthesis to enhance tomato cold tolerance through the SlERF4-7-SlSPI5-SlCPase-SlLOX5/SlAOS2 model. Overall, our study provides new perspectives on the role of galactinol in the JA regulatory network in plant adaptation to low-temperature stress.

Direct reprogramming of non-limb fibroblasts to cells with properties of limb progenitors

The early limb bud consists of mesenchymal limb progenitors derived from the lateral plate mesoderm (LPM). The LPM also gives rise to the mesodermal components of the flank and neck. However, the cells at these other levels cannot produce the variety of cell types found in the limb. Taking advantage of a direct reprogramming approach, we find a set of factors (Prdm16, Zbtb16, and Lin28a) normally expressed in the early limb bud and capable of imparting limb progenitor-like properties to mouse non-limb fibroblasts. The reprogrammed cells show similar gene expression profiles and can differentiate into similar cell types as endogenous limb progenitors. The further addition of Lin41 potentiates the proliferation of the reprogrammed cells. These results suggest that these same four factors may play pivotal roles in the specification of endogenous limb progenitors.

DendroX: multi-level multi-cluster selection in dendrograms

BACKGROUND

Cluster heatmaps are widely used in biology and other fields to uncover clustering patterns in data matrices. Most cluster heatmap packages provide utility functions to divide the dendrograms at a certain level to obtain clusters, but it is often difficult to locate the appropriate cut in the dendrogram to obtain the clusters seen in the heatmap or computed by a statistical method. Multiple cuts are required if the clusters locate at different levels in the dendrogram.

RESULTS

We developed DendroX, a web app that provides interactive visualization of a dendrogram where users can divide the dendrogram at any level and in any number of clusters and pass the labels of the identified clusters for functional analysis. Helper functions are provided to extract linkage matrices from cluster heatmap objects in R or Python to serve as input to the app. A graphic user interface was also developed to help prepare input files for DendroX from data matrices stored in delimited text files. The app is scalable and has been tested on dendrograms with tens of thousands of leaf nodes. As a case study, we clustered the gene expression signatures of 297 bioactive chemical compounds in the LINCS L1000 dataset and visualized them in DendroX. Seventeen biologically meaningful clusters were identified based on the structure of the dendrogram and the expression patterns in the heatmap. We found that one of the clusters consisting of mostly naturally occurring compounds is not previously reported and has its members sharing broad anticancer, anti-inflammatory and antioxidant activities.

CONCLUSIONS

DendroX solves the problem of matching visually and computationally determined clusters in a cluster heatmap and helps users navigate among different parts of a dendrogram. The identification of a cluster of naturally occurring compounds with shared bioactivities implicates a convergence of biological effects through divergent mechanisms.

Methane diffusion affects characteristics of benthic communities in and around microbial mat-covered sediments in the northeastern Japan sea

We investigated relationships between features of benthic macrofaunal communities and geochemical parameters in and around microbial mat-covered sediments associated with a methane seepage on Sakata Knoll in the eastern Japan Sea. A depression on top of the knoll corresponds to a gas-hydrate-bearing area with seepage of methane-rich fluid, and microbial mats cover the seafloor sediments. Sediment cores were collected at three sites for this study: one within a microbial mat, a second a few meters outside of the microbial mat, and a third from a reference site outside the gas-hydrate-bearing areas. Morphological analysis showed that the site inside the microbial mat had higher macrofaunal density and biomass compared with the other sites. 18S rRNA gene analysis showed that annelids were dominant in the surface sediment inside the microbial mat with the possible occurrence of microbial anaerobic oxidation of methane (AOM), whereas in the surface sediments outside the microbial mat and at the reference site the predominant species belonged to phylum Cercozoa. Morphological analysis also showed that the surface sediment inside the microbial mat noticeably favored annelids, with dorvilleid Ophryotrocha sp. and ampharetid Neosabellides sp. identified as major constituents. Statistical analysis showed that sulfidic sediment conditions with concentrations of H2S up to 121 μM resulting from AOM likely resulted in the predominance of annelids with tolerance to sulfide. Both the 18S rRNA genes and macrofaunal characteristics showed that benthic biodiversity among the three sites was greatest outside the microbial mat. The site outside the microbial mat may represent geochemical transition conditions, including a lower rate of upward methane gas-flow compared with the site inside the microbial mat. The high biodiversity there might result from the presence of species specifically suited to the transition zone as well as species also found in photosynthesis-based communities of the background environment.

Fluid chemistry alters faunal trophodynamics but not composition on the deep-sea Capelinhos hydrothermal edifice (Lucky Strike vent field, Mid-Atlantic Ridge)

The recently discovered deep-sea Capelinhos hydrothermal edifice, ~ 1.5 km of the main Lucky Strike (LS) vent field (northern Mid-Atlantic Ridge), contrasts with the other LS edifices in having poorly-altered end-member hydrothermal fluids with low pH and chlorine, and high metal concentrations. Capelinhos unique chemistry and location offer the opportunity to test the effects of local abiotic filters on faunal community structure while avoiding the often-correlated influence of dispersal limitation and depth. In this paper, we characterize for the first time the distribution patterns of the Capelinhos faunal communities, and analyze the benthic invertebrates (> 250 µm) inhabiting diffusive-flow areas and their trophic structures (δ13C, δ15N and δ34S). We hypothesized that faunal communities would differ from those of the nearest LS vent edifices, showing an impoverished species subset due to the potential toxicity of the chemical environment. Conversely, our results show that: (1) community distribution resembles that of other LS edifices, with assemblages visually dominated by shrimps (close to high-temperature focused-fluid areas) and mussels (at low-temperature diffuse flow areas); (2) most species from diffuse flow areas are well-known LS inhabitants, including the bed-forming and chemosymbiotic mussel Bathymodiolus azoricus and (3) communities are as diverse as those of the most diverse LS edifices. On the contrary, stable isotopes suggest different trophodynamics at Capelinhos. The high δ15N and, especially, δ13C and δ34S values suggest an important role of methane oxidation (i.e., methanotrophy), rather than the sulfide oxidation (i.e., thiotrophy) that predominates at most LS edifices. Our results indicate that Capelinhos shows unique environmental conditions, trophic structure and trophodynamics, yet similar fauna, compared to other LS edifices, which suggest a great environmental and trophic plasticity of the vent faunal communities at the LS.

Multipotent progenitors instruct ontogeny of the superior colliculus

The superior colliculus (SC) in the mammalian midbrain is essential for multisensory integration and is composed of a rich diversity of excitatory and inhibitory neurons and glia. However, the developmental principles directing the generation of SC cell-type diversity are not understood. Here, we pursued systematic cell lineage tracing in silico and in vivo, preserving full spatial information, using genetic mosaic analysis with double markers (MADM)-based clonal analysis with single-cell sequencing (MADM-CloneSeq). The analysis of clonally related cell lineages revealed that radial glial progenitors (RGPs) in SC are exceptionally multipotent. Individual resident RGPs have the capacity to produce all excitatory and inhibitory SC neuron types, even at the stage of terminal division. While individual clonal units show no pre-defined cellular composition, the establishment of appropriate relative proportions of distinct neuronal types occurs in a PTEN-dependent manner. Collectively, our findings provide an inaugural framework at the single-RGP/-cell level of the mammalian SC ontogeny.

Drivers of flower visit and resource sharing between the honeybee and native bees in Neotropical coastal sand dunes

The honeybee (Apis mellifera) is one of the most important pollinator species because it can gather resources from a vast variety of plant species, including both natives and introduced, across its geographical distribution. Although A. mellifera interacts with a large diversity of plants and shares resources with other pollinators, there are some plant species with which it interacts more frequently than others. Here, we evaluated the plant traits (i.e., plant length, abundance of bloomed individuals, number of open flowers, and stamen length) that would affect the honeybee visit frequencies to the flowers in a coastal environment in the Gulf of Mexico. Moreover, we evaluated which native bee species (and their body size) overlap floral resource with A. mellifera. We registered 998 plant-bee interactions between 35 plant species and 47 bee species. We observed that plant species with low height and with high abundances of bloomed individuals are positively related to a high frequency of visits by A. mellifera. Moreover, we found that A. mellifera tends to share a higher number of plant species with other bee species with a similar or smaller body size than with bigger species, which makes them a competitor for the resource with honeybees. Our results highlight that the impacts of A. mellifera on plants and native bees could be anticipated based on its individual's characteristics (i.e., plant height and abundance of bloomed individuals) and body size, respectively.

Dual RNA-sequencing of Fusarium head blight resistance in winter wheat

Fusarium head blight (FHB) is a devastating fungal disease responsible for significant yield losses in wheat and other cereal crops across the globe. FHB infection of wheat spikes results in grain contamination with mycotoxins, reducing both grain quality and yield. Breeding strategies have resulted in the production of FHB-resistant cultivars, however, the underlying molecular mechanisms of resistance in the majority of these cultivars are still poorly understood. To improve our understanding of FHB-resistance, we performed a transcriptomic analysis of FHB-resistant AC Emerson, FHB-moderately resistant AC Morley, and FHB-susceptible CDC Falcon in response to Fusarium graminearum. Wheat spikelets located directly below the point of inoculation were collected at 7-days post inoculation (dpi), where dual RNA-sequencing was performed to explore differential expression patterns between wheat cultivars in addition to the challenging pathogen. Differential expression analysis revealed distinct defense responses within FHB-resistant cultivars including the enrichment of physical defense through the lignin biosynthesis pathway, and DON detoxification through the activity of UDP-glycosyltransferases. Nucleotide sequence variants were also identified broadly between these cultivars with several variants being identified within differentially expressed putative defense genes. Further, F. graminearum demonstrated differential expression of mycotoxin biosynthesis pathways during infection, leading to the identification of putative pathogenicity factors.

Genome-wide investigation to assess copy number variants in the Italian local chicken population

BACKGROUND

Copy number variants (CNV) hold significant functional and evolutionary importance. Numerous ongoing CNV studies aim to elucidate the etiology of human diseases and gain insights into the population structure of livestock. High-density chips have enabled the detection of CNV with increased resolution, leading to the identification of even small CNV. This study aimed to identify CNV in local Italian chicken breeds and investigate their distribution across the genome.

RESULTS

Copy number variants were mainly distributed across the first six chromosomes and primarily associated with loss type CNV. The majority of CNV in the investigated breeds were of types 0 and 1, and the minimum length of CNV was significantly larger than that reported in previous studies. Interestingly, a high proportion of the length of chromosome 16 was covered by copy number variation regions (CNVR), with the major histocompatibility complex being the likely cause. Among the genes identified within CNVR, only those present in at least five animals across breeds (n = 95) were discussed to reduce the focus on redundant CNV. Some of these genes have been associated to functional traits in chickens. Notably, several CNVR on different chromosomes harbor genes related to muscle development, tissue-specific biological processes, heat stress resistance, and immune response. Quantitative trait loci (QTL) were also analyzed to investigate potential overlapping with the identified CNVR: 54 out of the 95 gene-containing regions overlapped with 428 QTL associated to body weight and size, carcass characteristics, egg production, egg components, fat deposition, and feed intake.

CONCLUSIONS

The genomic phenomena reported in this study that can cause changes in the distribution of CNV within the genome over time and the comparison of these differences in CNVR of the local chicken breeds could help in preserving these genetic resources.

Genetic and Phenotypic Consequences of Local Transitions between Sexual and Parthenogenetic Reproduction in the Wild

AbstractTransitions from sexual to asexual reproduction have occurred in numerous lineages, but it remains unclear why asexual populations rarely persist. In facultatively parthenogenetic animals, all-female populations can arise when males are absent or become extinct, and such populations could help to understand the genetic and phenotypic changes that occur in the initial stages of transitions to asexuality. We investigated a naturally occurring spatial mosaic of mixed-sex and all-female populations of the facultatively parthenogenetic Australian phasmid Megacrania batesii. Analysis of single-nucleotide polymorphisms indicated multiple independent transitions between reproductive modes. All-female populations had much lower heterozygosity and allelic diversity than mixed-sex populations, but we found few consistent differences in fitness-related traits between population types. All-female populations exhibited more frequent and severe deformities in their (flight-incapable) wings but did not show higher rates of appendage loss. All-female populations also harbored more ectoparasites in swamp (but not beach) habitats. Reproductive mode explained little variation in female body size, fecundity, or egg hatch rate. Our results suggest that transitions to parthenogenetic reproduction can lead to dramatic genetic changes with little immediate effect on performance. All-female M. batesii populations appear to consist of high-fitness genotypes that might be able to thrive for many generations in relatively constant and benign environments but could be vulnerable to environmental challenges, such as increased parasite abundance.

Our fragile future under the cumulative cultural evolution of two technologies

We derive and analyse a model with unusual features characterizing human activities over the long-run. First, human population dynamics draw heavily on consumer-resource modelling in ecology in that humans must consume biological resources to produce new humans. Second, the model also draws heavily from economic growth theory in that humans do not simply consume biological resources; they also produce the resources they consume. Finally, humans use two types of technology. Consumption technology affects the rate at which humans can extract resources. Production technology controls how effectively humans convert labour into new resources. The dynamics of both types of technology are subject to cumulative cultural evolutionary processes that allow both technological progress and regress. The resulting model exhibits a wide range of dynamical regimes. That said, the system is routinely sensitive to initial conditions, with wildly different outcomes given the same parameter values. Moreover, the system exhibits a basic fragility in the sense that human activities often lead to the endogenous extinction of the human species. This can happen gently, or it can follow periods of explosive human activity with super-exponential growth that ends in collapse. This article is part of the theme issue 'Evolution and sustainability: gathering the strands for an Anthropocene synthesis'.

Sex diversity in the 21st century: Concepts, frameworks, and approaches for the future of neuroendocrinology

Sex is ubiquitous and variable throughout the animal kingdom. Historically, scientists have used reductionist methodologies that rely on a priori sex categorizations, in which two discrete sexes are inextricably linked with gamete type. However, this binarized operationalization does not adequately reflect the diversity of sex observed in nature. This is due, in part, to the fact that sex exists across many levels of biological analysis, including genetic, molecular, cellular, morphological, behavioral, and population levels. Furthermore, the biological mechanisms governing sex are embedded in complex networks that dynamically interact with other systems. To produce the most accurate and scientifically rigorous work examining sex in neuroendocrinology and to capture the full range of sex variability and diversity present in animal systems, we must critically assess the frameworks, experimental designs, and analytical methods used in our research. In this perspective piece, we first propose a new conceptual framework to guide the integrative study of sex. Then, we provide practical guidance on research approaches for studying sex-associated variables, including factors to consider in study design, selection of model organisms, experimental methodologies, and statistical analyses. We invite fellow scientists to conscientiously apply these modernized approaches to advance our biological understanding of sex and to encourage academically and socially responsible outcomes of our work. By expanding our conceptual frameworks and methodological approaches to the study of sex, we will gain insight into the unique ways that sex exists across levels of biological organization to produce the vast array of variability and diversity observed in nature.

Methyltransferase-like (METTL) homologues participate in Nicotiana benthamiana antiviral responses

Methyltransferase (MTase) enzymes catalyze the addition of a methyl group to a variety of biological substrates. MTase-like (METTL) proteins are Class I MTases whose enzymatic activities contribute to the epigenetic and epitranscriptomic regulation of multiple cellular processes. N⁶-adenosine methylation (m⁶A) is a common chemical modification of eukaryotic and viral RNA whose abundance is jointly regulated by MTases and METTLs, demethylases, and m⁶A binding proteins. m⁶A affects various cellular processes including RNA degradation, post-transcriptional processing, and antiviral immunity. Here, we used Nicotiana benthamiana and plum pox virus (PPV), an RNA virus of the Potyviridae family, to investigated the roles of MTases in plant-virus interaction. RNA sequencing analysis identified MTase transcripts that are differentially expressed during PPV infection; among these, accumulation of a METTL gene was significantly downregulated. Two N. benthamiana METTL transcripts (NbMETTL1 and NbMETTL2) were cloned and further characterized. Sequence and structural analyses of the two encoded proteins identified a conserved S-adenosyl methionine (SAM) binding domain, showing they are SAM-dependent MTases phylogenetically related to human METTL16 and Arabidopsis thaliana FIONA1. Overexpression of NbMETTL1 and NbMETTL2 caused a decrease of PPV accumulation. In sum, our results indicate that METTL homologues participate in plant antiviral responses.

Effects of Age and Exposure Duration on the Sensitivity of Early Life Stage Fathead Minnow (Pimephales promelas) to Waterborne Propranolol Exposure

Propranolol is a heavily prescribed, nonspecific beta-adrenoceptor (bAR) antagonist frequently found in wastewater effluents, prompting concern over its potential to adversely affect exposed organisms. In the present study, the transcriptional responses of 4, 5, and 6 days postfertilization (dpf) ±1 h fathead minnow, exposed for 6, 24, or 48 h to 0.66 or 3.3 mg/L (nominal) propranolol were characterized using RNA sequencing. The number of differentially expressed genes (DEGs) was used as an estimate of sensitivity. A trend toward increased sensitivity with age was observed; fish >7 dpf at the end of exposure were particularly sensitive to propranolol. The DEGs largely overlapped among treatment groups, suggesting a highly consistent response that was independent of age. Cluster analysis was performed using normalized count data for unexposed and propranolol-exposed fish. Control fish clustered tightly by age, with fish ≥7 dpf clustering away from younger fish, reflecting developmental differences. When clustering was conducted using exposed fish, in cases where propranolol induced a minimal or no transcriptional response, the results mirrored those of the control fish and did not appreciably cluster by treatment. In treatment groups that displayed a more robust transcriptional response, the effects of propranolol were evident; however, fish <7 dpf clustered away from older fish, despite having similar numbers of DEGs. Increased sensitivity at 7 dpf coincided with developmental milestones with the potential to alter propranolol pharmacokinetics or pharmacodynamics, such as the onset of exogenous feeding and gill functionality as well as increased systemic expression of bAR. These results may have broader implications because toxicity testing often utilizes fish <4 dpf, prior to the onset of these potentially important developmental milestones, which may result in an underestimation of risk for some chemicals. Environ Toxicol Chem 2024;00:1-14. Published 2023. This article is a U.S. Government work and is in the public domain in the USA.

Characterization of urban sources of antibiotics and antibiotic-resistance genes in a Dutch sewer catchment

A one year study was conducted in the city of Nijmegen, The Netherlands, to characterize various urban sources of antibiotics and antibiotic resistant genes (ARGs) in wastewater within a single sewer catchment. Prevalence of ermB, tet(W), sul1, sul2, intl1, and 16S rRNA gene was determined at 10 locations within the city. Sampling locations included a nursing home, a student residence, a hospital and an industrial area, among others. Wastewater concentrations of 23 antibiotics were measured using passive sampling. Additionally, excreted loads of 22 antibiotics were estimated based on ambulatory prescription and clinical usage data. Genes sul1 and intl1 were most abundant across most locations. Ciprofloxacin and amoxicillin together contributed over 92 % of the total estimated antibiotic selective pressure at all sampling points. The present study highlights the prominent role that hospitals can have in the prevalence and proliferation of ARGs in urban wastewater. Furthermore, results suggest that even short-term changes in the therapeutic regimen prescribed in hospitals may translate into shifting ARG abundance patterns in hospital wastewater. The methods applied present an opportunity to identify emission hotspots and prioritize intervention options to limit ARG spread from urban wastewater to the environment.

Analysis of geometric morphometrics and molecular phylogeny for Anopheles species in the Republic of Korea

Human malaria, transmitted by Anopheles mosquitoes, is the most predominant mosquito-borne disease that is responsible for hundreds of thousands of deaths worldwide each year. In the Republic of Korea (ROK), there are currently several hundred malaria cases annually, mostly near the demilitarized zone (DMZ). Eight species of Anopheles mosquitoes are currently known to be present in the ROK. Similar to other major malaria vectors in Africa and India, it is very challenging to morphologically differentiate Anopheles mosquitoes in the ROK due to their extremely similar morphology. In this study, wing geometric morphometrics (WGM) were used to differentiate the eight Anopheles species collected at six locations near the DMZ, Seoul and Pyeongtaek from April-October 2021. Phylogenetic analysis was also performed using cytochrome c oxidase subunit 1 (COI), internal transcribed spacer 2 (ITS2), and tyrosine hydroxylase (TH) genes for comparison with WGM analysis and to infer evolutionary relationships. The results of cross-validation (overall accuracy = 74.8%) demonstrated that species identification using WGM alone was not possible with a high accuracy for all eight species. While phylogenetic analyses based on the COI region could not clearly distinguish some species, the analysis based on ITS2 and TH was more useful for resolving the phylogenetic correlation of the eight species. Our results may improve Anopheles species identification strategies for effective identification and control of malaria vectors in the ROK.

Early life social complexity shapes adult neural processing in the communal spiny mouse Acomys cahirinus

RATIONALE

Early life social rearing has profound consequences on offspring behavior and resilience. Yet, most studies examining early life development in rodents use species whose young are born immobile and do not produce complex social behavior until later in development. Furthermore, models of rearing under increased social complexity, rather than deprivation, are needed to provide alternative insight into the development of social neural circuitry.

OBJECTIVES

To understand precocial offspring social development, we manipulated early life social complexity in the communal spiny mouse Acomys cahirinus and assessed long-term consequences on offspring social behavior, exploration, and neural responses to novel social stimuli.

METHODS

Spiny mouse pups were raised in the presence or absence of a non-kin breeding group. Upon adulthood, subjects underwent social interaction tests, an open field test, and a novel object test. Subjects were then exposed to a novel conspecific and novel group and neural responses were quantified via immunohistochemical staining in brain regions associated with social behavior.

RESULTS

Early life social experience did not influence behavior in the test battery, but it did influence social processing. In animals exposed to non-kin during development, adult lateral septal neural responses toward a novel conspecific were weaker and hypothalamic neural responses toward a mixed-sex group were stronger.

CONCLUSIONS

Communal species may exhibit robust behavioral resilience to the early life social environment. But the early life environment can affect how novel social information is processed in the brain during adulthood, with long-term consequences that are likely to shape their behavioral trajectory.

Assessing over decadal biomass burning influence on particulate matter composition in subequatorial Amazon: literature review, remote sensing, chemical speciation and machine learning application

A study on aerosols in the Brazilian subequatorial Amazon region, Tangará da Serra (TS) and Alta Floresta (AF) was conducted and compared to findings in an additional site with background characteristics (Manaus, AM). TS and AF counties suffer from intense biomass burning periods in the dry season, and it accounts for high levels of particles in the atmosphere. Chemical characterization of fine and coarse particulate matter (PM) was performed to quantify water-soluble ions (WSI) and black carbon (BC). The importance of explanatory variables was assessed using three machine learning techniques. Average concentrations of PM in AF and TS were similar (PM2.0, 17±10 µg m-3 (AF) and 16±11 µg m-3 (TS) and PM10-2.0, 13±5 µg m-3 (AF) and 11±7 µg m-3 (TS)), but higher than the background site. BC and SO4 2- were the prevalent components as they represented 27%-68% of particulates chemical composition. The combination of the machine learning techniques provided a further understanding of the pathways for PM concentration variability, and the results highlighted the influence of biomass burning for key sample groups and periods. PM2.0, BC, and most WSI presented higher concentrations in the dry season, providing further support for the influence of biomass burning.

Deep learning with citizen science data enables estimation of species diversity and composition at continental extents

Effective solutions to conserve biodiversity require accurate community- and species-level information at relevant, actionable scales and across entire species' distributions. However, data and methodological constraints have limited our ability to provide such information in robust ways. Herein we employ a Deep-Reasoning Network implementation of the Deep Multivariate Probit Model (DMVP-DRNets), an end-to-end deep neural network framework, to exploit large observational and environmental data sets together and estimate landscape-scale species diversity and composition at continental extents. We present results from a novel year-round analysis of North American avifauna using data from over nine million eBird checklists and 72 environmental covariates. We highlight the utility of our information by identifying critical areas of high species diversity for a single group of conservation concern, the North American wood warblers, while capturing spatiotemporal variation in species' environmental associations and interspecific interactions. In so doing, we demonstrate the type of accurate, high-resolution information on biodiversity that deep learning approaches such as DMVP-DRNets can provide and that is needed to inform ecological research and conservation decision-making at multiple scales.

Multi-omic responses of fish exposed to complex chemical mixtures in the Shenandoah River watershed

To evaluate relationships between different anthropogenic impacts, contaminant occurrence, and fish health, we conducted in situ fish exposures across the Shenandoah River watershed at five sites with different land use. Exposure water was analyzed for over 500 chemical constituents, and organismal, metabolomic, and transcriptomic endpoints were measured in fathead minnows. Adverse reproductive outcomes were observed in fish exposed in the upper watershed at both wastewater treatment plant (WWTP) effluent- and agriculture-impacted sites, including decreased gonadosomatic index and altered secondary sex characteristics. This was accompanied with increased mortality at the site most impacted by agricultural activities. Molecular biomarkers of estrogen exposure were unchanged and consistent with low or non-detectable concentrations of common estrogens, indicating that alternative mechanisms were involved in organismal adverse outcomes. Hepatic metabolomic and transcriptomic profiles were altered in a site-specific manner, consistent with variation in land use and contaminant profiles. Integrated biomarker response data were useful for evaluating mechanistic linkages between contaminants and adverse outcomes, suggesting that reproductive endocrine disruption, altered lipid processes, and immunosuppression may have been involved in these organismal impacts. This study demonstrated linkages between human-impact, contaminant occurrence, and exposure effects in the Shenandoah River watershed and showed increased risk of adverse outcomes in fathead minnows exposed to complex mixtures at sites impacted by municipal wastewater discharges and agricultural practices.

Is apomixis occurring in walnut (Juglans regia L.)? New data from progeny molecular tests and cytological investigations shed light on its reproductive system

Introduction

Persian walnut (Juglans regia) is an economically important nut fruit species cultivated worldwide for its nutritious kernel and timber quality wood. Walnut trees are mostly hetero-dichogamous and, depending on the genotype, some cultivars are protogynous, while others are protandrous. Although selfing is possible when male and female blooms overlap, the dichogamy of the species promotes outcrossing. In addition to sexual reproduction, some reports indicate that elements of apomixis may occur in commercial orchards of walnut varieties and in the last two decades, nut production by apomixis has been reported in walnut. However, there are no reliable studies on the occurrence of apomictic reproduction based on cytoembryological observations and/or molecular marker-progeny tests. This study addresses the combined use of molecular and cytological analyses to gain new insights into the population genetics and reproduction systems of J. regia.

Methods

We systematically analyzed the reproductive origin of individual progeny plants from 8 different cultivated walnut genotypes using microsatellite genotyping and carried out cytohistological investigations of 5 cultivated walnut genotypes arising seed sets from isolated flowers, to shed light on the mode of reproduction.

Results and discussion

These cytometric and genotyping analyses did not support any asexual mode of reproduction or asexual propagation by seed and all individuals studied were identified as zygotic plants produced by crossing. Likewise, the cytological findings did not confirm completely the first component of apomixis, namely apomeiosis. On the other hand, according to histological evidence, adventitious embryony seems to take place at low frequency. Overall, our findings suggest that the occurrence of gametophytic apomixis is unlikely in J. regia, but sporophytic apomixis cannot be completely ruled out.

Transcriptional immune suppression and upregulation of double stranded DNA damage and repair repertoires in ecDNA-containing tumors

Extrachromosomal DNA is a common cause of oncogene amplification in cancer. The non-chromosomal inheritance of ecDNA enables tumors to rapidly evolve, contributing to treatment resistance and poor outcome for patients. The transcriptional context in which ecDNAs arise and progress, including chromosomally-driven transcription, is incompletely understood. We examined gene expression patterns of 870 tumors of varied histological types, to identify transcriptional correlates of ecDNA. Here we show that ecDNA containing tumors impact four major biological processes. Specifically, ecDNA containing tumors upregulate DNA damage and repair, cell cycle control, and mitotic processes, but downregulate global immune regulation pathways. Taken together, these results suggest profound alterations in gene regulation in ecDNA containing tumors, shedding light on molecular processes that give rise to their development and progression.

Resolving network clusters disparity based on dissimilarity measurements with nonmetric analysis of variance

Classic ANOVA (cA) tests the explanatory power of a partitioning on a set of objects. More fit for clusters proximity analysis, nonparametric ANOVA (npA) extends to a case where instead of the object values themselves, their mutual distances are available. However, extending the cA applicability, the metric conditions in npA are limiting. Based on the central limit theorem (CLT), here we introduce nonmetric ANOVA (nmA) that by relaxing the metric properties between objects, allows an ANOVA-like statistical testing of a network clusters disparity. We present a parametric test statistic which under the null hypothesis of no differences between the competing clusters means, follows an exact F-distribution. We apply our method on three diverse biological examples, discuss its parallel performance, and note the specific use of each method tailored by the inherent data properties. The R code is provided at github.com/AmiryousefiLab/nmANOVA.

16S metabarcoding analysis reveals the influence of organic and conventional farming practices on bacterial communities from the rhizospheric of Coffea arabica L

Coffea sp. is cultivated in many tropical countries. Brazil has always adopted intensive agricultural practices, but organic coffee farming is an alternative system based on the non-use of agrochemicals and the rational management of soils. Metabarcoding 16S analysis using next-generation sequencing has been developed to identify and compare the diversity of the Coffea arabica L. rhizospheric bacterial community in two farming areas in São Paulo, Brazil. Dourado uses conventional farming, while Ribeirão Corrente uses organic. We found broad taxonomic composition, with sequences from 24 phyla, 55 classes, 61 orders, 146 families, and 337genus. The three most abundant phyla were Proteobacteria (38.27%), Actinobacteria (15.56%), and Acidobacteria (16.10%). In organic farming, the top 3 were the family Sphingomonadaceae, order Rhizobiales, genus Nocardioides, and Gp6. The genus Gp2 and the phylum Candidatus Saccharibacteria were the most abundant OTUs exclusively present in conventional farming. In the organic farming practice, Proteobacteria, Actinobacteria, and Acidobacteria were also present among the exclusive OTUs; we also found OTUs belonging to Bacteroidetes, Firmicutes, and Verrucomicrobia. Our study indicates a positive effect of organic farming on microbial communities. Fertilization may directly affect soil microbiota, suggesting that a large and active microbial community low in functional diversity might not adapt to new climatic conditions. A diverse community could provide better resilience to environmental changes, improving the productivity of this important crop.

Activation patterns of dopaminergic cell populations reflect different learning scenarios in a cichlid fish, Pseudotropheus zebra

Dopamine is present in all vertebrates and the functional roles of the subsystems are assumed to be similar. Whereas the effect of dopaminergic modulation is well investigated in different target systems, less is known about the factors that are causing the modulation of dopaminergic cells. Using the zebra mbuna, Pseudotropheus zebra, a cichlid fish from Lake Malawi as a model system, we investigated the activation of specific dopaminergic cell populations detected by double-labeling with TH and pS6 antibodies while the animals were solving different learning tasks. Specifically, we compared an intense avoidance learning situation, an instrumental learning task, and a non-learning isolated group and found strong activation of different dopaminergic cell populations. Preoptic-hypothalamic cell populations respond to the stress component in the avoidance task, and the forced movement/locomotion may be responsible for activation in the posterior tubercle. The instrumental learning task had little stress component, but the activation of the raphe superior in this group may be correlated with attention or arousal during the training sessions. At the same time, the weaker activation of the nucleus of the posterior commissure may be related to positive reward acting onto tectal circuits. Finally, we examined the co-activation patterns across all dopaminergic cell populations and recovered robust differences across experimental groups, largely driven by hypothalamic, posterior tubercle, and brain stem regions possibly encoding the valence and salience associated with stressful stimuli. Taken together, our results offer some insights into the different functions of the dopaminergic cell populations in the brain of a non-mammalian vertebrate in correlation with different behavioral conditions, extending our knowledge for a more comprehensive view of the mechanisms of dopaminergic modulation in vertebrates.