Qualimap: evaluating next-generation sequencing alignment data

lncRNA LINC00941 modulates MTA2/NuRD occupancy to suppress premature human epidermal differentiation

Numerous long non-coding RNAs (lncRNAs) were shown to have a functional impact on cellular processes such as human epidermal homeostasis. However, the mechanism of action for many lncRNAs remains unclear to date. Here, we report that lncRNA LINC00941 regulates keratinocyte differentiation on an epigenetic level through association with the NuRD complex, one of the major chromatin remodelers in cells. We find that LINC00941 interacts with NuRD-associated MTA2 and CHD4 in human primary keratinocytes. LINC00941 perturbation changes MTA2/NuRD occupancy at bivalent chromatin domains in close proximity to transcriptional regulator genes, including the EGR3 gene coding for a transcription factor regulating epidermal differentiation. Notably, LINC00941 depletion resulted in reduced NuRD occupancy at the EGR3 gene locus, increased EGR3 expression in human primary keratinocytes, and increased abundance of EGR3-regulated epidermal differentiation genes in cells and human organotypic epidermal tissues. Our results therefore indicate a role of LINC00941/NuRD in repressing EGR3 expression in non-differentiated keratinocytes, consequentially preventing premature differentiation of human epidermal tissues.

Exploring the gene expression network involved in the heat stress response of a thermotolerant tomato genotype

BACKGROUND

The increase in temperatures due to the current climate change dramatically affects crop cultivation, resulting in yield losses and altered fruit quality. Tomato is one of the most extensively grown and consumed horticultural products, and although it can withstand a wide range of climatic conditions, heat stress can affect plant growth and development specially on the reproductive stage, severely influencing the final yield. In the present work, the heat stress response mechanisms of one thermotolerant genotype (E42) were investigated by exploring its regulatory gene network. This was achieved through a promoter analysis based on the identification of the heat stress elements (HSEs) mapping in the promoters, combined with a gene co-expression network analysis aimed at identifying interactions among heat-related genes.

RESULTS

Results highlighted 82 genes presenting HSEs in the promoter and belonging to one of the 52 gene networks obtained by the GCN analysis; 61 of these also interact with heat shock factors (Hsfs). Finally, a list of 13 candidate genes including two Hsfs, nine heat shock proteins (Hsps) and two GDSL esterase/lipase (GELPs) were retrieved by focusing on those E42 genes exhibiting HSEs in the promoters, interacting with Hsfs and showing variants, compared to Heinz reference genome, with HIGH and/or MODERATE impact on the translated protein. Among these, the Gene Ontology annotation analysis evidenced that only LeHsp100 (Solyc02g088610) belongs to a network specifically involved in the response to heat stress.

CONCLUSIONS

As a whole, the combination of bioinformatic analyses carried out on genomic and trascriptomic data available for tomato, together with polymorphisms detected in HS-related genes of the thermotolerant E42 allowed to determine a subset of candidate genes involved in the HS response in tomato. This study provides a novel approach in the investigation of abiotic stress response mechanisms and further studies will be conducted to validate the role of the highlighted genes.

Prediction models of persistent taxane-induced peripheral neuropathy among breast cancer survivors using whole-exome sequencing

Persistent taxane-induced peripheral neuropathy (TIPN) is highly prevalent among early-stage breast cancer survivors (ESBCS) and has detrimental effect on quality of life. We leveraged logistic regression models to develop and validate polygenic prediction models to estimate the risk of persistent PN symptoms in a training cohort and validation cohort taking clinical risk factors into account. Based on 337 whole-exome sequenced ESBCS two of five prediction models for individual PN symptoms obtained AUC results above 60% when validated. Using the model for numbness in feet (35 SNVs) in the test cohort, 73% survivors were correctly predicted. For tingling in feet (55 SNVs) 70% were correctly predicted. Both models included SNVs from the ADAMTS20, APT6V0A2, CCDC88C, CYP2C8, EPHA5, NR1H3, PSKH2/APTV0D2, and SCN10A genes. For cramps in feet, difficulty climbing stairs and difficulty opening a jar the validation was unsuccessful. Polygenic prediction models including clinical risk factors can estimate the risk of persistent taxane-induced numbness in feet and tingling in feet in ESBCS.

Overcoming microbial inhibition of S. sonnei through the exploitation of genomically-predicted antibiotic resistance profiles for the development of food enrichment media

Linking outbreaks of Shigella spp. to specific foods is challenging due to poor selectivity of current enrichment media. We have previously shown that enrichment media, tailored to the genomically-predicted antimicrobial resistance (AMR) of Shiga toxigenic E. coli strains, enhances their isolation from foods. This study investigates the application of this approach for Shigella isolation. The AMR gene profiles of 21,908 published S. sonnei genomes indicated a high prevalence of genes conferring resistance to streptomycin (aadA, aph(3")-Ib, aph(6)-Id , 92.8%), sulfonamides (sul1, sul2, 74.8%) and/or trimethoprim (dfrA, 96.2%). Genomic analysis and antibiotic susceptibility testing conducted with a panel of 17 outbreak-associated S. sonnei strains confirmed the correlation of AMR gene detection with resistance phenotypes. Supplementation of Shigella Broth (SB) with up to 400 µg/mL of trimethoprim or sulfadiazine did not supress growth of sensitive strains, whereas 100 µg/mL of streptomycin increased the selectivity of this broth. All three antibiotics increased the selectivity of modified Tryptone Soya Broth (mTSB). Based on these results, supplemented media formulations were developed and assessed by measuring the relative growth of S. sonnei in cultures co-inoculated with a strain of bacteriocin-producing E. coli that is inhibitory to Shigella growth. S. sonnei was not recovered from co-cultures grown in SB or mTSB without antibiotics. In contrast, media supplemented with streptomycin at 50 and 100 µg/mL, trimethoprim at 25 and 50 µg/mL, and sulfadiazine at 100 µg/mL increased the relative proportion of S. sonnei in post-enrichment cultures. The enhanced recovery of resistant S. sonnei strains achieved in this study indicates that, in cases where genomic data is available for clinical S. sonnei isolates, customization of selective enrichment media based on AMR gene detection could be a valuable tool for supporting the investigation of foodborne shigellosis outbreaks.

Colonial-driven extinction of the blue antelope despite genomic adaptation to low population size

Low genomic diversity is generally indicative of small population size and is considered detrimental by decreasing long-term adaptability.1,2,3,4,5,6 Moreover, small population size may promote gene flow with congeners and outbreeding depression.7,8,9,10,11,12,13 Here, we examine the connection between habitat availability, effective population size (Ne), and extinction by generating a 40× nuclear genome from the extinct blue antelope (Hippotragus leucophaeus). Historically endemic to the relatively small Cape Floristic Region in southernmost Africa,14,15 populations were thought to have expanded and contracted across glacial-interglacial cycles, tracking suitable habitat.16,17,18 However, we found long-term low Ne, unaffected by glacial cycles, suggesting persistence with low genomic diversity for many millennia prior to extinction in ∼AD 1800. A lack of inbreeding, alongside high levels of genetic purging, suggests adaptation to this long-term low Ne and that human impacts during the colonial era (e.g., hunting and landscape transformation), rather than longer-term ecological processes, were central to its extinction. Phylogenomic analyses uncovered gene flow between roan (H. equinus) and blue antelope, as well as between roan and sable antelope (H. niger), approximately at the time of divergence of blue and sable antelope (∼1.9 Ma). Finally, we identified the LYST and ASIP genes as candidates for the eponymous bluish pelt color of the blue antelope. Our results revise numerous aspects of our understanding of the interplay between genomic diversity and evolutionary history and provide the resources for uncovering the genetic basis of this extinct species' unique traits.

Parasympathetic neurons derived from human pluripotent stem cells model human diseases and development

Autonomic parasympathetic neurons (parasymNs) control unconscious body responses, including "rest-and-digest." ParasymN innervation is important for organ development, and parasymN dysfunction is a hallmark of autonomic neuropathy. However, parasymN function and dysfunction in humans are vastly understudied due to the lack of a model system. Human pluripotent stem cell (hPSC)-derived neurons can fill this void as a versatile platform. Here, we developed a differentiation paradigm detailing the derivation of functional human parasymNs from Schwann cell progenitors. We employ these neurons (1) to assess human autonomic nervous system (ANS) development, (2) to model neuropathy in the genetic disorder familial dysautonomia (FD), (3) to show parasymN dysfunction during SARS-CoV-2 infection, (4) to model the autoimmune disease Sjögren's syndrome (SS), and (5) to show that parasymNs innervate white adipocytes (WATs) during development and promote WAT maturation. Our model system could become instrumental for future disease modeling and drug discovery studies, as well as for human developmental studies.

Characterization of MLST-99 Salmonella Typhimurium and the monophasic variant I:4,[5],12:i:- isolated from Canadian Atlantic coast shellfish

Salmonella enterica subsp. enterica Typhimurium and its monophasic variant I 1;4,[5],12:i:- (MVST) are responsible for thousands of reported cases of salmonellosis each year in Canada, and countries worldwide. We investigated S. Typhimurium and MVST isolates recovered from raw shellfish harvested in Atlantic Canada by the Canadian Food Inspection Agency (CFIA) over the past decade, to assess the potential impact of these isolates on human illness and to explore possible routes of shellfish contamination. Whole-genome sequence analysis was performed on 210 isolates of S. Typhimurium and MVST recovered from various food sources, including shellfish. The objective was to identify genetic markers linked to ST-99, a sequence type specifically associated with shellfish, which could explain their high prevalence in shellfish. We also investigated the genetic similarity amongst CFIA ST-99 isolates recovered in different years and geographical locations. Finally, the study aimed to enhance the molecular serotyping of ST-99 isolates, as they are serologically classified as MVST but are frequently misidentified as S. Typhimurium through sequence analysis. To ensure recovery of ST-99 from shellfish was not due to favourable growth kinetics, we measured the growth rates of these isolates relative to other Salmonella and determined that ST-99 did not have a faster growth rate and/or shorter lag phase than other Salmonella evaluated. The CFIA ST-99 isolates from shellfish were highly clonal, with up to 81 high-quality single nucleotide variants amongst isolates. ST-99 isolates both within the CFIA collection and those isolated globally carried numerous unique deletions, insertions and mutations in genes, including some considered important for virulence, such as gene deletions in the type VI secretion system. Interestingly, several of these genetic characteristics appear to be unique to North America. Most notably was a large genomic region showing a high prevalence in genomes from Canadian isolates compared to those from the USA. Although the functions of the majority of the proteins encoded within this region remain unknown, the genes umuC and umuD, known to be protective against UV light damage, were present. While this study did not specifically examine the effects of mutations and insertions, results indicate that these isolates may be adapted to survive in specific environments, such as ocean water, where wild birds and/or animals serve as the natural hosts. Our hypothesis is reinforced by a global phylogenetic analysis, which indicates that isolates obtained from North American shellfish and wild birds are infrequently connected to isolates from human sources. These findings suggest a distinct ecological niche for ST-99, potentially indicating their specialization and adaptation to non-human hosts and environments, such as oceanic habitats.

Effects of Chronic Inflammatory Activation of Murine and Human Arterial Endothelial Cells at Normal Lipoprotein and Cholesterol Levels In Vivo and In Vitro

The activation of endothelial cells is crucial for immune defense mechanisms but also plays a role in the development of atherosclerosis. We have previously shown that inflammatory stimulation of endothelial cells on top of elevated lipoprotein/cholesterol levels accelerates atherogenesis. The aim of the current study was to investigate how chronic endothelial inflammation changes the aortic transcriptome of mice at normal lipoprotein levels and to compare this to the inflammatory response of isolated endothelial cells in vitro. We applied a mouse model expressing constitutive active IκB kinase 2 (caIKK2)-the key activator of the inflammatory NF-κB pathway-specifically in arterial endothelial cells and analyzed transcriptomic changes in whole aortas, followed by pathway and network analyses. We found an upregulation of cell death and mitochondrial beta-oxidation pathways with a predicted increase in endothelial apoptosis and necrosis and a simultaneous reduction in protein synthesis genes. The highest upregulated gene was ACE2, the SARS-CoV-2 receptor, which is also an important regulator of blood pressure. Analysis of isolated human arterial and venous endothelial cells supported these findings and also revealed a reduction in DNA replication, as well as repair mechanisms, in line with the notion that chronic inflammation contributes to endothelial dysfunction.

Microsatellite instability in mismatch repair proficient colorectal cancer: clinical features and underlying molecular mechanisms

BACKGROUND

Both defects in mismatch repair (dMMR) and high microsatellite instability (MSI-H) have been recognised as crucial biomarkers that guide treatment strategies and disease management in colorectal cancer (CRC). As MMR and MSI tests are being widely conducted, an increasing number of MSI-H tumours have been identified in CRCs with mismatch repair proficiency (pMMR). The objective of this study was to assess the clinical features of patients with pMMR/MSI-H CRC and elucidate the underlying molecular mechanism in these cases.

METHODS

From January 2015 to December 2018, 1684 cases of pMMR and 401 dMMR CRCs were enrolled. Of those patients, 93 pMMR/MSI-H were identified. The clinical phenotypes and prognosis were analysed. Frozen and paraffin-embedded tissue were available in 35 patients with pMMR/MSI-H, for which comprehensive genomic and transcriptomic analyses were performed.

FINDINGS

In comparison to pMMR/MSS CRCs, pMMR/MSI-H CRCs exhibited significantly less tumour progression and better long-term prognosis. The pMMR/MSI-H cohorts displayed a higher presence of CD8+ T cells and NK cells when compared to the pMMR/MSS group. Mutational signature analysis revealed that nearly all samples exhibited deficiencies in MMR genes, and we also identified deleterious mutations in MSH3-K383fs.

INTERPRETATION

This study revealed pMMR/MSI-H CRC as a distinct subgroup within CRC, which manifests diverse clinicopathological features and long-term prognostic outcomes. Distinct features in the tumour immune-microenvironment were observed in pMMR/MSI-H CRCs. Pathogenic deleterious mutations in MSH3-K383fs were frequently detected, suggesting another potential biomarker for identifying MSI-H.

FUNDING

This work was supported by the Science and Technology Commission of Shanghai Municipality (20DZ1100101).

Complete genome sequence of Pseudomonas canadensis strain Pcan-CK-23 isolated from Zophobas morio larvae

We present the complete genome sequence of Pseudomonas canadensis. The strain (Pcan-CK-23) was isolated from Zophobas morio (superworm) larvae. The genome consisted of a 6,424,469 bp chromosome with a GC content of 60.3% and 5,973 genes. Pcan-CK-23 can be used as a reference genome for further studies with P. canadensis.

Droplet based whole genome amplification for sequencing minute amounts of purified Mycobacterium tuberculosis DNA

Implementation of whole genome sequencing (WGS) for patient care is hindered by limited Mycobacterium tuberculosis (Mtb) in clinical specimens and slow Mtb growth. We evaluated droplet multiple displacement amplification (dMDA) for amplification of minute amounts of Mtb DNA to enable WGS as an alternative to other Mtb enrichment methods. Purified genomic Mtb-DNA (0.1, 0.5, 1, and 5 pg) was encapsulated and amplified using the Samplix Xdrop-instrument and sequenced alongside a control sample using standard Illumina protocols followed by MAGMA-analysis. The control and 5 pg input dMDA samples underwent nanopore sequencing followed by Nanoseq and TB-profiler analysis. dMDA generated 105-2400 ng DNA from the 0.1-5 pg input DNA, respectively. Followed by Illumina WGS, dMDA raised mean sequencing depth from 7 × for 0.1 pg input DNA to ≥ 60 × for 5 pg input and the control sample. Bioinformatic analysis revealed a high number of false positive and false negative variants when amplifying ≤ 0.5 pg input DNA. Nanopore sequencing of the 5 pg dMDA sample presented excellent coverage depth, breadth, and accurate strain characterization, albeit elevated false positive and false negative variants compared to Illumina-sequenced dMDA sample with identical Mtb DNA input. dMDA coupled with Illumina WGS for samples with ≥ 5 pg purified Mtb DNA, equating to approximately 1000 copies of the Mtb genome, offers precision for drug resistance, phylogeny, and transmission insights.

The transcription factor STE12 influences growth on several carbon sources and production of dehydroacetic acid (DHAA) in Trichoderma reesei

The filamentous ascomycete Trichoderma reesei, known for its prolific cellulolytic enzyme production, recently also gained attention for its secondary metabolite synthesis. Both processes are intricately influenced by environmental factors like carbon source availability and light exposure. Here, we explore the role of the transcription factor STE12 in regulating metabolic pathways in T. reesei in terms of gene regulation, carbon source utilization and biosynthesis of secondary metabolites. We show that STE12 is involved in regulating cellulase gene expression and growth on carbon sources associated with iron homeostasis. STE12 impacts gene regulation in a light dependent manner on cellulose with modulation of several CAZyme encoding genes as well as genes involved in secondary metabolism. STE12 selectively influences the biosynthesis of the sorbicillinoid trichodimerol, while not affecting the biosynthesis of bisorbibutenolide, which was recently shown to be regulated by the MAPkinase pathway upstream of STE12 in the signaling cascade. We further report on the biosynthesis of dehydroacetic acid (DHAA) in T. reesei, a compound known for its antimicrobial properties, which is subject to regulation by STE12. We conclude, that STE12 exerts functions beyond development and hence contributes to balance the energy distribution between substrate consumption, reproduction and defense.

Multicellular ecotypes shape progression of lung adenocarcinoma from ground-glass opacity toward advanced stages

Lung adenocarcinoma is a type of cancer that exhibits a wide range of clinical radiological manifestations, from ground-glass opacity (GGO) to pure solid nodules, which vary greatly in terms of their biological characteristics. Our current understanding of this heterogeneity is limited. To address this gap, we analyze 58 lung adenocarcinoma patients via machine learning, single-cell RNA sequencing (scRNA-seq), and whole-exome sequencing, and we identify six lung multicellular ecotypes (LMEs) correlating with distinct radiological patterns and cancer cell states. Notably, GGO-associated neoantigens in early-stage cancers are recognized by CD8+ T cells, indicating an immune-active environment, while solid nodules feature an immune-suppressive LME with exhausted CD8+ T cells, driven by specific stromal cells such as CTHCR1+ fibroblasts. This study also highlights EGFR(L858R) neoantigens in GGO samples, suggesting potential CD8+ T cell activation. Our findings offer valuable insights into lung adenocarcinoma heterogeneity, suggesting avenues for targeted therapies in early-stage disease.

Comparison of qPCR and metagenomic sequencing methods for quantifying antibiotic resistance genes in wastewater

Surveillance methods of circulating antibiotic resistance genes (ARGs) are of utmost importance in order to tackle what has been described as one of the greatest threats to humanity in the 21st century. In order to be effective, these methods have to be accurate, quickly deployable, and scalable. In this study, we compare metagenomic shotgun sequencing (TruSeq DNA sequencing) of wastewater samples with a state-of-the-art PCR-based method (Resistomap HT-qPCR) on four wastewater samples that were taken from hospital, industrial, urban and rural areas. ARGs that confer resistance to 11 antibiotic classes have been identified in these wastewater samples using both methods, with the most abundant observed classes of ARGs conferring resistance to aminoglycoside, multidrug-resistance (MDR), macrolide-lincosamide-streptogramin B (MLSB), tetracycline and beta-lactams. In comparing the methods, we observed a strong correlation of relative abundance of ARGs obtained by the two tested methods for the majority of antibiotic classes. Finally, we investigated the source of discrepancies in the results obtained by the two methods. This analysis revealed that false negatives were more likely to occur in qPCR due to mutated primer target sites, whereas ARGs with incomplete or low coverage were not detected by the sequencing method due to the parameters set in the bioinformatics pipeline. Indeed, despite the good correlation between the methods, each has its advantages and disadvantages which are also discussed here. By using both methods together, a more robust ARG surveillance program can be established. Overall, the work described here can aid wastewater treatment plants that plan on implementing an ARG surveillance program.

Effects of nitrous oxide and ketamine on electrophysiological and molecular responses in the prefrontal cortex of mice: A comparative study

Nitrous oxide (N2O; laughing gas) has recently reported to produce rapid antidepressant effects, but little is known about the underlying mechanisms. We performed transcriptomics, in situ hybridization, and electrophysiological studies to examine the potential shared signatures induced by 1 h inhalation of 50% N2O and a single subanesthetic dose of ketamine (10 mg/kg, i.p.) in the medial prefrontal cortex (mPFC) in adult mice. Both treatments similarly affected the transcription of several negative regulators of mitogen-activated protein kinases (MAPKs), namely, dual specificity phosphatases (DUSPs). The effects were primarily located in the pyramidal cells. Notably, the overall effects of N2O on mRNA expression were much more prominent and widespread compared to ketamine. Ketamine caused an elevation of the spiking frequency of putative pyramidal neurons and increased gamma activity (30-100 Hz) of cortical local field potentials. However, N2O produced no such effects. Spiking amplitudes and spike-to-local field potential phase locking of putative pyramidal neurons and interneurons in this brain area showed no uniform changes across treatments. Our findings suggest that N2O and subanesthetic-dose ketamine target MAPK pathway in the mPFC but produce varying acute electrophysiological responses.

Multimodal integration of neuroimaging and genetic data for the diagnosis of mood disorders based on computer vision models

Mood disorders, particularly major depressive disorder (MDD) and bipolar disorder (BD), are often underdiagnosed, leading to substantial morbidity. Harnessing the potential of emerging methodologies, we propose a novel multimodal fusion approach that integrates patient-oriented brain structural magnetic resonance imaging (sMRI) scans with DNA whole-exome sequencing (WES) data. Multimodal data fusion aims to improve the detection of mood disorders by employing established deep-learning architectures for computer vision and machine-learning strategies. We analyzed brain imaging genetic data of 321 East Asian individuals, including 147 patients with MDD, 78 patients with BD, and 96 healthy controls. We developed and evaluated six fusion models by leveraging common computer vision models in image classification: Vision Transformer (ViT), Inception-V3, and ResNet50, in conjunction with advanced machine-learning techniques (XGBoost and LightGBM) known for high-dimensional data analysis. Model validation was performed using a 10-fold cross-validation. Our ViT ⊕ XGBoost fusion model with MRI scans, genomic Single Nucleotide polymorphism (SNP) data, and unweighted polygenic risk score (PRS) outperformed baseline models, achieving an incremental area under the curve (AUC) of 0.2162 (32.03% increase) and 0.0675 (+8.19%) and incremental accuracy of 0.1455 (+25.14%) and 0.0849 (+13.28%) compared to SNP-only and image-only baseline models, respectively. Our findings highlight the opportunity to refine mood disorder diagnostics by demonstrating the transformative potential of integrating diverse, yet complementary, data modalities and methodologies.

Treatment with the apoptosis inhibitor Asunercept reduces clone sizes in patients with lower risk Myelodysplastic Neoplasms

In low-risk Myelodysplastic Neoplasms (MDS), increased activity of apoptosis-promoting factors such as tumor necrosis factor (TNFα) and pro-apoptotic Fas ligand (CD95L) have been described as possible pathomechanisms leading to impaired erythropoiesis. Asunercept (APG101) is a novel therapeutic fusion protein blocking CD95, which has previously shown partial efficacy in reducing transfusion requirement in a clinical phase I trial for low-risk MDS patients (NCT01736436; 2012-11-26). In the current study we aimed to evaluate the effect of Asunercept therapy on the clonal bone marrow composition to identify potential biomarkers to predict response. Bone marrow samples of n = 12 low-risk MDS patients from the above referenced clinical trial were analyzed by serial deep whole exome sequencing in a total of n = 58 time points. We could distinguish a mean of 3.5 molecularly defined subclones per patient (range 2-6). We observed a molecular response defined as reductions of dominant clone sizes by a variant allele frequency (VAF) decrease of at least 10% (mean 20%, range: 10.5-39.2%) in dependency of Asunercept treatment in 9 of 12 (75%) patients. Most of this decline in clonal populations was observed after completion of 12 weeks treatment. Particularly early and pronounced reductions of clone sizes were found in subclones driven by mutations in genes involved in regulation of methylation (n = 1 DNMT3A, n = 1 IDH2, n = 1 TET2). Our results suggest that APG101 could be efficacious in reducing clone sizes of mutated hematopoietic cells in the bone marrow of Myelodysplastic Neoplasms, which warrants further investigation.

Whole-exome sequencing of Nigerian benign prostatic hyperplasia reveals increased alterations in apoptotic pathways

BACKGROUND

Through whole-exome sequencing of 60 formalin-fixed paraffin-embedded Nigerian (NGRn) benign prostatic hyperplasia (BPH) samples, we identified germline and somatic alterations in apoptotic pathways impacting BPH development and progression. Prostate enlargement is a common occurrence in male aging; however, this enlargement can lead to lower urinary tract symptoms that negatively impact quality of life. This impact is disproportionately present in men of African ancestry. BPH pathophysiology is poorly understood and studies examining non-European populations are lacking.

METHODS

In this study, NGRn BPH, normal prostate, and prostate cancer (PCa) tumor samples were sequenced and compared to characterize genetic alterations in NGRn BPH.

RESULTS

Two hundred and two nonbenign, ClinVar-annotated germline variants were present in NGRn BPH samples. Six genes [BRCA1 (92%), HSD3B1 (85%), TP53 (37%), PMS2 (23%), BARD1 (20%), and BRCA2 (17%)] were altered in at least 10% of samples; however, compared to NGRn normal and tumor, the frequency of alterations in BPH samples showed no significant differences at the gene or variant level. BRCA2_rs11571831 and TP53_rs1042522 germline alterations had a statistically significant co-occurrence interaction in BPH samples. In at least two BPH samples, 173 genes harbored somatic variants known to be clinically actionable. Three genes (COL18A1, KIF16B, and LRP1) showed a statistically significant (p < 0.05) higher frequency in BPH. NGRn BPH also had five gene pairs (PKD1/KIAA0100, PKHD1/PKD1, DNAH9/LRP1B, NWD1/DCHS2, and TCERG1/LMTK2) with statistically significant co-occurring interactions. Two hundred and seventy-nine genes contained novel somatic variants in NGRn BPH. Three genes (CABP1, FKBP1C, and RP11-595B24.2) had a statistically significant (p < 0.05) higher alteration frequency in NGRn BPH and three were significantly higher in NGRn tumor (CACNA1A, DMKN, and CACNA2D2). Pairwise Fisher's exact tests showed 14 gene pairs with statistically significant (p < 0.05) interactions and four interactions approaching significance (p < 0.10). Mutational patterns in NGRn BPH were similar to COSMIC (Catalog of Somatic Mutations in Cancer) signatures associated with aging and dysfunctional DNA damage repair.

CONCLUSIONS

NGRn BPH contained significant germline alteration interactions (BRCA2_rs11571831 and TP53_rs1042522) and increased somatic alteration frequencies (LMTK2, LRP1, COL18A1, CABP1, and FKBP1C) that impact apoptosis. Normal prostate development is maintained by balancing apoptotic and proliferative activity. Dysfunction in either mechanism can lead to abnormal prostate growth. This work is the first to examine genomic sequencing in NGRn BPH and provides data that fill known gaps in the understanding BPH and how it impacts men of African ancestry.

Rescue of Mycobacterium bovis DNA Obtained from Cultured Samples during Official Surveillance of Animal TB: Key Steps for Robust Whole Genome Sequence Data Generation

Epidemiological surveillance of animal tuberculosis (TB) based on whole genome sequencing (WGS) of Mycobacterium bovis has recently gained track due to its high resolution to identify infection sources, characterize the pathogen population structure, and facilitate contact tracing. However, the workflow from bacterial isolation to sequence data analysis has several technical challenges that may severely impact the power to understand the epidemiological scenario and inform outbreak response. While trying to use archived DNA from cultured samples obtained during routine official surveillance of animal TB in Portugal, we struggled against three major challenges: the low amount of M. bovis DNA obtained from routinely processed animal samples; the lack of purity of M. bovis DNA, i.e., high levels of contamination with DNA from other organisms; and the co-occurrence of more than one M. bovis strain per sample (within-host mixed infection). The loss of isolated genomes generates missed links in transmission chain reconstruction, hampering the biological and epidemiological interpretation of data as a whole. Upon identification of these challenges, we implemented an integrated solution framework based on whole genome amplification and a dedicated computational pipeline to minimize their effects and recover as many genomes as possible. With the approaches described herein, we were able to recover 62 out of 100 samples that would have otherwise been lost. Based on these results, we discuss adjustments that should be made in official and research laboratories to facilitate the sequential implementation of bacteriological culture, PCR, downstream genomics, and computational-based methods. All of this in a time frame supporting data-driven intervention.

Mitochondrial and Nuclear DNA Variants in Amyotrophic Lateral Sclerosis: Enrichment in the Mitochondrial Control Region and Sirtuin Pathway Genes in Spinal Cord Tissue

Amyotrophic Lateral Sclerosis (ALS) is a progressive disease with prevalent mitochondrial dysfunctions affecting both upper and lower motor neurons in the motor cortex, brainstem, and spinal cord. Despite mitochondria having their own genome (mtDNA), in humans, most mitochondrial genes are encoded by the nuclear genome (nDNA). Our study aimed to simultaneously screen for nDNA and mtDNA genomes to assess for specific variant enrichment in ALS compared to control tissues. Here, we analysed whole exome (WES) and whole genome (WGS) sequencing data from spinal cord tissues, respectively, of 6 and 12 human donors. A total of 31,257 and 301,241 variants in nuclear-encoded mitochondrial genes were identified from WES and WGS, respectively, while mtDNA reads accounted for 73 and 332 variants. Despite technical differences, both datasets consistently revealed a specific enrichment of variants in the mitochondrial Control Region (CR) and in several of these genes directly associated with mitochondrial dynamics or with Sirtuin pathway genes within ALS tissues. Overall, our data support the hypothesis of a variant burden in specific genes, highlighting potential actionable targets for therapeutic interventions in ALS.

Isolation of Burkholderia pseudomallei from a goat in New Caledonia: implications for animal and human health monitoring and serological tool comparison

BACKGROUND

Melioidosis is a serious bacterial infection caused by Burkholderia pseudomallei, a gram-negative bacterium commonly found in soil and water. It can affect both humans and animals, and is endemic in regions such as Southeast Asia and Northern Australia. In recent years, there have been reports of an emergence of human melioidosis in other areas, including New Caledonia.

RESULTS

During standard laboratory analysis in New Caledonia in 2021, a strain of B. pseudomallei was isolated from a goat. The strain was characterized using both MLST and WGS techniques and was found to cluster with previously described local human strains from the area. In parallel, several serological tests (CFT, ELISA, Luminex (Hcp1, GroEL, BPSS1840), arrays assay and a latex agglutination test) were performed on animals from the farm where the goat originated, and/or from three other neighboring farms. Using two commercial ELISA kits, seropositive animals were found only on the farm where the infected goat originated and tests based on recombinant proteins confirmed the usefulness of the Hcp1 protein for the diagnosis of melioidosis in animals.

CONCLUSIONS

Despite the regular reports of human cases, this is the first confirmed case of melioidosis in an animal in New Caledonia. These results confirm the presence of the bacterium in the region and highlight the importance of vigilance for both animal and human health. It is critical that all health partners, including breeders, veterinarians, and biologists, work together to monitor and prevent the spread of the disease.

Shared and distinct interactions of type 1 and type 2 Epstein-Barr Nuclear Antigen 2 with the human genome

BACKGROUND

There are two major genetic types of Epstein-Barr Virus (EBV): type 1 (EBV-1) and type 2 (EBV-2). EBV functions by manipulating gene expression in host B cells, using virus-encoded gene regulatory proteins including Epstein-Barr Nuclear Antigen 2 (EBNA2). While type 1 EBNA2 is known to interact with human transcription factors (hTFs) such as RBPJ, EBF1, and SPI1 (PU.1), type 2 EBNA2 shares only ~ 50% amino acid identity with type 1 and thus may have distinct binding partners, human genome binding locations, and functions.

RESULTS

In this study, we examined genome-wide EBNA2 binding in EBV-1 and EBV-2 transformed human B cells to identify shared and unique EBNA2 interactions with the human genome, revealing thousands of type-specific EBNA2 ChIP-seq peaks. Computational predictions based on hTF motifs and subsequent ChIP-seq experiments revealed that both type 1 and 2 EBNA2 co-occupy the genome with SPI1 and AP-1 (BATF and JUNB) hTFs. However, type 1 EBNA2 showed preferential co-occupancy with EBF1, and type 2 EBNA2 preferred RBPJ. These differences in hTF co-occupancy revealed possible mechanisms underlying type-specific gene expression of known EBNA2 human target genes: MYC (shared), CXCR7 (type 1 specific), and CD21 (type 2 specific). Both type 1 and 2 EBNA2 binding events were enriched at systemic lupus erythematosus (SLE) and multiple sclerosis (MS) risk loci, while primary biliary cholangitis (PBC) risk loci were specifically enriched for type 2 peaks.

CONCLUSIONS

This study reveals extensive type-specific EBNA2 interactions with the human genome, possible differences in EBNA2 interaction partners, and a possible new role for type 2 EBNA2 in autoimmune disorders. Our results highlight the importance of considering EBV type in the control of human gene expression and disease-related investigations.

Draft genomes of meropenem-resistant Pseudomonads from the cloacae of turtles in Lake Opinicon, Canada

We report the draft genome sequences of seven meropenem-resistant bacterial isolates from the cloacae of painted and musk turtles at Lake Opinicon, Canada. This study improves our understanding of the spread of meropenem-resistant bacteria in the wild.

Acetyl-CoA production by Mediator-bound 2-ketoacid dehydrogenases boosts de novo histone acetylation and is regulated by nitric oxide

Histone-modifying enzymes depend on the availability of cofactors, with acetyl-coenzyme A (CoA) being required for histone acetyltransferase (HAT) activity. The discovery that mitochondrial acyl-CoA-producing enzymes translocate to the nucleus suggests that high concentrations of locally synthesized metabolites may impact acylation of histones and other nuclear substrates, thereby controlling gene expression. Here, we show that 2-ketoacid dehydrogenases are stably associated with the Mediator complex, thus providing a local supply of acetyl-CoA and increasing the generation of hyper-acetylated histone tails. Nitric oxide (NO), which is produced in large amounts in lipopolysaccharide-stimulated macrophages, inhibited the activity of Mediator-associated 2-ketoacid dehydrogenases. Elevation of NO levels and the disruption of Mediator complex integrity both affected de novo histone acetylation within a shared set of genomic regions. Our findings indicate that the local supply of acetyl-CoA generated by 2-ketoacid dehydrogenases bound to Mediator is required to maximize acetylation of histone tails at sites of elevated HAT activity.

Whole-genome sequencing and comparative genomics reveal candidate genes associated with quality traits in Dioscorea alata

BACKGROUND

Quality traits are essential determinants of consumer preferences. Dioscorea alata (Greater Yam), is a starchy tuber crop in tropical regions. However, a comprehensive understanding of the genetic basis underlying yam tuber quality remains elusive. To address this knowledge gap, we employed population genomics and candidate gene association approaches to unravel the genetic factors influencing the quality attributes of boiled yam.

METHODS AND RESULTS

Comparative genomics analysis of 45 plant species revealed numerous novel genes absent in the existing D. alata gene annotation. This approach, adding 48% more genes, significantly enhanced the functional annotation of three crucial metabolic pathways associated with boiled yam quality traits: pentose and glucuronate interconversions, starch and sucrose metabolism, and flavonoid biosynthesis. In addition, the whole-genome sequencing of 127 genotypes identified 27 genes under selection and 22 genes linked to texture, starch content, and color through a candidate gene association analysis. Notably, five genes involved in starch content and cell wall composition, including 1,3-beta Glucan synthase, β-amylase, and Pectin methyl esterase, were common to both approaches and their expression levels were assessed by transcriptomic data.

CONCLUSIONS

The analysis of the whole-genome of 127 genotypes of D. alata and the study of three specific pathways allowed the identification of important genes for tuber quality. Our findings provide insights into the genetic basis of yam quality traits and will help the enhancement of yam tuber quality through breeding programs.

Evaluation of aphid resistance on different rose cultivars and transcriptome analysis in response to aphid infestation

BACKGROUND

The rose is one of the most important ornamental flowers in the world for its aesthetic beauty but can be attacked by many pests such as aphids. Aphid infestation causes tremendous damage on plant tissues leading to harmed petals and leaves. Rose cultivars express different levels of resistance to aphid infestation yet the information remains unclear. Not only that, studies about the transcriptional analysis on defending mechanisms against aphids in rose are limited so far.

RESULTS

In this study, the aphid resistance of 20 rose cultivars was evaluated, and they could be sorted into six levels based on the number ratio of aphids. And then, a transcriptome analysis was conducted after aphid infestation in one high resistance (R, Harmonie) and one highly susceptibility (S, Carefree Wonder) rose cultivar. In open environment the majority of rose cultivars had the highest aphid number at May 6th or May 15th in 2020 and the resistance to infestation could be classified into six levels. Differential expression analysis revealed that there were 1,626 upregulated and 767 downregulated genes in the R cultivar and 481 upregulated and 63 downregulated genes in the S cultivar after aphid infestation. Pathway enrichment analysis of the differentially expressed genes revealed that upregulated genes in R and S cultivars were both enriched in defense response, biosynthesis of secondary metabolites (phenylpropanoid, alkaloid, and flavonoid), carbohydrate metabolism (galactose, starch, and sucrose metabolism) and lipid processing (alpha-linolenic acid and linolenic acid metabolism) pathways. In the jasmonic acid metabolic pathway, linoleate 13S-lipoxygenase was specifically upregulated in the R cultivar, while genes encoding other crucial enzymes, allene oxide synthase, allene oxide cyclase, and 12-oxophytodienoate reductase were upregulated in both cultivars. Transcription factor analysis and transcription factor binding search showed that WRKY transcription factors play a pivotal role during aphid infestation in the R cultivar.

CONCLUSIONS

Our study indicated the potential roles of jasmonic acid metabolism and WRKY transcription factors during aphid resistance in rose, providing clues for future research.

Enhancer mutations modulate the severity of chemotherapy-induced myelosuppression

Non-small cell lung cancer is often diagnosed at advanced stages, and many patients are still treated with classical chemotherapy. The unselective nature of chemotherapy often results in severe myelosuppression. Previous studies showed that protein-coding mutations could not fully explain the predisposition to myelosuppression. Here, we investigate the possible role of enhancer mutations in myelosuppression susceptibility. We produced transcriptome and promoter-interaction maps (using HiCap) of three blood stem-like cell lines treated with carboplatin or gemcitabine. Taking advantage of publicly available enhancer datasets, we validated HiCap results in silico and in living cells using epigenetic CRISPR technology. We also developed a network approach for interactome analysis and detection of differentially interacting genes. Differential interaction analysis provided additional information on relevant genes and pathways for myelosuppression compared with differential gene expression analysis at the bulk level. Moreover, we showed that enhancers of differentially interacting genes are highly enriched for variants associated with differing levels of myelosuppression. Altogether, our work represents a prominent example of integrative transcriptome and gene regulatory datasets analysis for the functional annotation of noncoding mutations.

Karyotypic stasis and swarming influenced the evolution of viral tolerance in a species-rich bat radiation

The emergence of COVID-19 and severe acute respiratory syndrome (SARS) has prioritized understanding bats' viral tolerance. Myotis bats are exceptionally species rich and have evolved viral tolerance. They also exhibit swarming, a cryptic behavior where large, multi-species assemblages gather for mating, which has been hypothesized to promote interspecific hybridization. To resolve the coevolution of genome architecture and their unusual antiviral tolerance, we undertook a phylogenomic analysis of 60 Old World Myotis genomes. We demonstrate an extensive history of introgressive hybridization that has replaced the species phylogeny across 17%-93% of the genome except for pericentromeric regions of macrochromosomes. Introgression tracts were enriched on microchromosome regions containing key antiviral pathway genes overexpressed during viral challenge experiments. Together, these results suggest that the unusual Myotis karyotype may have evolved to selectively position immune-related genes in high recombining genomic regions prone to introgression of divergent alleles, including a diversity of interleukin loci responsible for the release of pro-inflammatory cytokines.

Ehf and Fezf2 regulate late medullary thymic epithelial cell and thymic tuft cell development

Thymic epithelial cells are indispensable for T cell maturation and selection and the induction of central immune tolerance. The self-peptide repertoire expressed by medullary thymic epithelial cells is in part regulated by the transcriptional regulator Aire (Autoimmune regulator) and the transcription factor Fezf2. Due to the high complexity of mTEC maturation stages (i.e., post-Aire, Krt10+ mTECs, and Dclk1+ Tuft mTECs) and the heterogeneity in their gene expression profiles (i.e., mosaic expression patterns), it has been challenging to identify the additional factors complementing the transcriptional regulation. We aimed to identify the transcriptional regulators involved in the regulation of mTEC development and self-peptide expression in an unbiased and genome-wide manner. We used ATAC footprinting analysis as an indirect approach to identify transcription factors involved in the gene expression regulation in mTECs, which we validated by ChIP sequencing. This study identifies Fezf2 as a regulator of the recently described thymic Tuft cells (i.e., Tuft mTECs). Furthermore, we identify that transcriptional regulators of the ELF, ESE, ERF, and PEA3 subfamily of the ETS transcription factor family and members of the Krüppel-like family of transcription factors play a role in the transcriptional regulation of genes involved in late mTEC development and promiscuous gene expression.

Insights from Melipona bicolor hybrid genome assembly: a stingless bee genome with chromosome-level scaffold

BACKGROUND

The highly eusocial stingless bees are crucial pollinators of native and agricultural ecosystems. Nevertheless, genomic studies within this bee tribe remain scarce. We present the genome assembly of the stingless bee Melipona bicolor. This bee is a remarkable exception to the typical single-queen colony structure, since in this species, multiple queens may coexist and share reproductive duties, resulting in genetically diverse colonies with weak kinship connections. As the only known genuinely polygynous bee, M. bicolor's genome provides a valuable resource for investigating sociality beyond kin selection.

RESULTS

The genome was assembled employing a hybrid approach combining short and long reads, resulting in 241 contigs spanning 259 Mb (N50 of 6.2 Mb and 97.5% complete BUSCOs). Comparative analyses shed light on some evolutionary aspects of stingless bee genomics, including multiple chromosomal rearrangements in Melipona. Additionally, we explored the evolution of venom genes in M. bicolor and other stingless bees, revealing that, apart from two genes, the conserved repertoire of venom components remains under purifying selection in this clade.

CONCLUSION

This study advances our understanding of stingless bee genomics, contributing to the conservation efforts of these vital pollinators and offering insights into the evolutionary mechanisms driving their unique adaptations.

A pairwise cytokine code explains the organism-wide response to sepsis

Sepsis is a systemic response to infection with life-threatening consequences. Our understanding of the molecular and cellular impact of sepsis across organs remains rudimentary. Here, we characterize the pathogenesis of sepsis by measuring dynamic changes in gene expression across organs. To pinpoint molecules controlling organ states in sepsis, we compare the effects of sepsis on organ gene expression to those of 6 singles and 15 pairs of recombinant cytokines. Strikingly, we find that the pairwise effects of tumor necrosis factor plus interleukin (IL)-18, interferon-gamma or IL-1β suffice to mirror the impact of sepsis across tissues. Mechanistically, we map the cellular effects of sepsis and cytokines by computing changes in the abundance of 195 cell types across 9 organs, which we validate by whole-mouse spatial profiling. Our work decodes the cytokine cacophony in sepsis into a pairwise cytokine message capturing the gene, cell and tissue responses of the host to the disease.

Three Rounds of Read Correction Significantly Improve Eukaryotic Protein Detection in ONT Reads

BACKGROUND

Eukaryotes' whole-genome sequencing is crucial for species identification, gene detection, and protein annotation. Oxford Nanopore Technology (ONT) is an affordable and rapid platform for sequencing eukaryotes; however, the relatively higher error rates require computational and bioinformatic efforts to produce more accurate genome assemblies. Here, we evaluated the effect of read correction tools on eukaryote genome completeness, gene detection and protein annotation.

METHODS

Reads generated by ONT of four eukaryotes, C. albicans, C. gattii, S. cerevisiae, and P. falciparum, were assembled using minimap2 and underwent three rounds of read correction using flye, medaka and racon. The generates consensus FASTA files were compared for total length (bp), genome completeness, gene detection, and protein-annotation by QUAST, BUSCO, BRAKER1 and InterProScan, respectively.

RESULTS

Genome completeness was dependent on the assembly method rather than on the read correction tool; however, medaka performed better than flye and racon. Racon significantly performed better than flye and medaka in gene detection, while both racon and medaka significantly performed better than flye in protein-annotation.

CONCLUSION

We show that three rounds of read correction significantly affect gene detection and protein annotation, which are dependent on assembly quality in preference to assembly completeness.

Altered gene expression linked to germline dysfunction following exposure of Caenorhabditis elegans to DEET

N,N-diethyl-meta-toluamide (DEET) is a commonly used synthetic insect repellent. Although the neurological effects of DEET have been widely investigated, its effects on the germline are less understood. Here, we show that exposure of the nematode Caenorhabditis elegans, which is highly predictive of mammalian reprotoxicity, resulting in internal DEET levels within the range detected in human biological samples, causes activation of p53/CEP-1-dependent germ cell apoptosis, altered meiotic recombination, chromosome abnormalities, and missegregation. RNA-sequencing analysis links DEET-induced alterations in the expression of genes related to redox processes and chromatin structure to reduced mitochondrial function, impaired DNA double-strand break repair progression, and defects during early embryogenesis. We propose that Caenorhabditis elegans exposure to DEET interferes with gene expression, leading to increased oxidative stress and altered chromatin structure, resulting in germline effects that pose a risk to reproductive health.

UV RESISTANCE LOCUS 8-Mediated UV-B Response Is Required Alongside CRYPTOCHROME 1 for Plant Survival in Sunlight under Field Conditions

As sessile, photoautotrophic organisms, plants are subjected to fluctuating sunlight that includes potentially detrimental ultraviolet-B (UV-B) radiation. Experiments under controlled conditions have shown that the UV-B photoreceptor UV RESISTANCE LOCUS 8 (UVR8) controls acclimation and tolerance to UV-B in Arabidopsis thaliana; however, its long-term impact on plant fitness under naturally fluctuating environments remain poorly understood. Here, we quantified the survival and reproduction of different Arabidopsis mutant genotypes under diverse field and laboratory conditions. We found that uvr8 mutants produced more fruits than wild type when grown in growth chambers under artificial low-UV-B conditions but not under natural field conditions, indicating a fitness cost in the absence of UV-B stress. Importantly, independent double mutants of UVR8 and the blue light photoreceptor gene CRYPTOCHROME 1 (CRY1) in two genetic backgrounds showed a drastic reduction in fitness in the field. Experiments with UV-B attenuation in the field and with supplemental UV-B in growth chambers demonstrated that UV-B caused the cry1 uvr8 conditional lethal phenotype. Using RNA-seq data of field-grown single and double mutants, we explicitly identified genes showing significant statistical interaction of UVR8 and CRY1 mutations in the presence of UV-B in the field. They were enriched in Gene Ontology categories related to oxidative stress, photoprotection and DNA damage repair in addition to UV-B response. Our study demonstrates the functional importance of the UVR8-mediated response across life stages in natura, which is partially redundant with that of cry1. Moreover, these data provide an integral picture of gene expression associated with plant responses under field conditions.

Differential Gene Expression in Human Fibroblasts Simultaneously Exposed to Ionizing Radiation and Simulated Microgravity

During future space missions, astronauts will be exposed to cosmic radiation and microgravity (μG), which are known to be health risk factors. To examine the differentially expressed genes (DEG) and their prevalent biological processes and pathways as a response to these two risk factors simultaneously, 1BR-hTERT human fibroblast cells were cultured under 1 gravity (1G) or simulated μG for 48 h in total and collected at 0 (sham irradiated), 3 or 24 h after 1 Gy of X-ray or Carbon-ion (C-ion) irradiation. A three-dimensional clinostat was used for the simulation of μG and the simultaneous radiation exposure of the samples. The RNA-seq method was used to produce lists of differentially expressed genes between different environmental conditions. Over-representation analyses were performed and the enriched biological pathways and targeting transcription factors were identified. Comparing sham-irradiated cells under simulated μG and 1G conditions, terms related to response to oxygen levels and muscle contraction were identified. After irradiation with X-rays or C-ions under 1G, identified DEGs were found to be involved in DNA damage repair, signal transduction by p53 class mediator, cell cycle arrest and apoptosis pathways. The same enriched pathways emerged when cells were irradiated under simulated μG condition. Nevertheless, the combined effect attenuated the transcriptional response to irradiation which may pose a subtle risk in space flights.

Human Sterols Are Overproduced, Stored and Excreted in Yeasts

Sterols exert a profound influence on numerous cellular processes, playing a crucial role in both health and disease. However, comprehending the effects of sterol dysfunction on cellular physiology is challenging. Consequently, numerous processes affected by impaired sterol biosynthesis still elude our complete understanding. In this study, we made use of yeast strains that produce cholesterol instead of ergosterol and investigated the cellular response mechanisms on the transcriptome as well as the lipid level. The exchange of ergosterol for cholesterol caused the downregulation of phosphatidylethanolamine and phosphatidylserine and upregulation of phosphatidylinositol and phosphatidylcholine biosynthesis. Additionally, a shift towards polyunsaturated fatty acids was observed. While the sphingolipid levels dropped, the total amounts of sterols and triacylglycerol increased, which resulted in 1.7-fold enlarged lipid droplets in cholesterol-producing yeast cells. In addition to internal storage, cholesterol and its precursors were excreted into the culture supernatant, most likely by the action of ABC transporters Snq2, Pdr12 and Pdr15. Overall, our results demonstrate that, similarly to mammalian cells, the production of non-native sterols and sterol precursors causes lipotoxicity in K. phaffii, mainly due to upregulated sterol biosynthesis, and they highlight the different survival and stress response mechanisms on multiple, integrative levels.

High-Quality Genome Assemblies Reveal Evolutionary Dynamics of Repetitive DNA and Structural Rearrangements in the Drosophila virilis Subgroup

High-quality genome assemblies across a range of nontraditional model organisms can accelerate the discovery of novel aspects of genome evolution. The Drosophila virilis group has several attributes that distinguish it from more highly studied species in the Drosophila genus, such as an unusual abundance of repetitive elements and extensive karyotype evolution, in addition to being an attractive model for speciation genetics. Here, we used long-read sequencing to assemble five genomes of three virilis group species and characterized sequence and structural divergence and repetitive DNA evolution. We find that our contiguous genome assemblies allow characterization of chromosomal arrangements with ease and can facilitate analysis of inversion breakpoints. We also leverage a small panel of resequenced strains to explore the genomic pattern of divergence and polymorphism in this species and show that known demographic histories largely predicts the extent of genome-wide segregating polymorphism. We further find that a neo-X chromosome in Drosophila americana displays X-like levels of nucleotide diversity. We also found that unusual repetitive elements were responsible for much of the divergence in genome composition among species. Helitron-derived tandem repeats tripled in abundance on the Y chromosome in D. americana compared to Drosophila novamexicana, accounting for most of the difference in repeat content between these sister species. Repeats with characteristics of both transposable elements and satellite DNAs expanded by 3-fold, mostly in euchromatin, in both D. americana and D. novamexicana compared to D. virilis. Our results represent a major advance in our understanding of genome biology in this emerging model clade.

New indicators for monitoring genetic diversity applied to alpine brown trout populations using whole genome sequence data

International policy recently adopted commitments to maintain genetic diversity in wild populations to secure their adaptive potential, including metrics to monitor temporal trends in genetic diversity - so-called indicators. A national programme for assessing trends in genetic diversity was recently initiated in Sweden. Relating to this effort, we systematically assess contemporary genome-wide temporal trends (40 years) in wild populations using the newly adopted indicators and whole genome sequencing (WGS). We use pooled and individual WGS data from brown trout (Salmo trutta) in eight alpine lakes in protected areas. Observed temporal trends in diversity metrics (nucleotide diversity, Watterson's ϴ and heterozygosity) lie within proposed acceptable threshold values for six of the lakes, but with consistently low values in lakes above the tree line and declines observed in these northern-most lakes. Local effective population size is low in all lakes, highlighting the importance of continued protection of interconnected systems to allow genetic connectivity for long-term viability of these populations. Inbreeding (FROH ) spans 10%-30% and is mostly represented by ancient (<1 Mb) runs of homozygosity, with observations of little change in mutational load. We also investigate adaptive dynamics over evolutionarily short time frames (a few generations); identifying putative parallel selection across all lakes within a gene pertaining to skin pigmentation as well as candidates of selection unique to specific lakes and lake systems involved in reproduction and immunity. We demonstrate the utility of WGS for systematic monitoring of natural populations, a priority concern if genetic diversity is to be protected.

Metabolic memory of Δ9-tetrahydrocannabinol exposure in pluripotent stem cells and primordial germ cells-like cells

Cannabis, the most consumed illicit psychoactive drug in the world, is increasingly used by pregnant women. However, while cannabinoid receptors are expressed in the early embryo, the impact of phytocannabinoids exposure on early embryonic processes is lacking. Here, we leverage a stepwise in vitro differentiation system that captures the early embryonic developmental cascade to investigate the impact of exposure to the most abundant phytocannabinoid, Δ9-tetrahydrocannabinol (Δ9-THC). We demonstrate that Δ9-THC increases the proliferation of naive mouse embryonic stem cells (ESCs) but not of their primed counterpart. Surprisingly, this increased proliferation, dependent on the CB1 receptor binding, is only associated with moderate transcriptomic changes. Instead, Δ9-THC capitalizes on ESCs' metabolic bivalence by increasing their glycolytic rates and anabolic capabilities. A memory of this metabolic rewiring is retained throughout differentiation to Primordial Germ Cell-Like Cells in the absence of direct exposure and is associated with an alteration of their transcriptional profile. These results represent the first in-depth molecular characterization of the impact of Δ9-THC exposure on early stages of germline development.

Direct long read visualization reveals metabolic interplay between two antimalarial drug targets

Increases in the copy number of large genomic regions, termed genome amplification, are an important adaptive strategy for malaria parasites. Numerous amplifications across the Plasmodium falciparum genome contribute directly to drug resistance or impact the fitness of this protozoan parasite. During the characterization of parasite lines with amplifications of the dihydroorotate dehydrogenase (DHODH) gene, we detected increased copies of an additional genomic region that encompassed 3 genes (~5 kb) including GTP cyclohydrolase I (GCH1 amplicon). While this gene is reported to increase the fitness of antifolate resistant parasites, GCH1 amplicons had not previously been implicated in any other antimalarial resistance context. Here, we further explored the association between GCH1 and DHODH copy number. Using long read sequencing and single read visualization, we directly observed a higher number of tandem GCH1 amplicons in parasites with increased DHODH copies (up to 9 amplicons) compared to parental parasites (3 amplicons). While all GCH1 amplicons shared a consistent structure, expansions arose in 2-unit steps (from 3 to 5 to 7, etc copies). Adaptive evolution of DHODH and GCH1 loci was further bolstered when we evaluated prior selection experiments; DHODH amplification was only successful in parasite lines with pre-existing GCH1 amplicons. These observations, combined with the direct connection between metabolic pathways that contain these enzymes, lead us to propose that the GCH1 locus is beneficial for the fitness of parasites exposed to DHODH inhibitors. This finding highlights the importance of studying variation within individual parasite genomes as well as biochemical connections of drug targets as novel antimalarials move towards clinical approval.

TWAS revealed significant causal loci for milk production and its composition in Murrah buffaloes

Milk yield is the most complex trait in dairy animals, and mapping all causal variants even with smallest effect sizes has been difficult with the genome-wide association study (GWAS) sample sizes available in geographical regions with small livestock holdings such as Indian sub-continent. However, Transcriptome-wide association studies (TWAS) could serve as an alternate for fine mapping of expression quantitative trait loci (eQTLs). This is a maiden attempt to identify milk production and its composition related genes using TWAS in Murrah buffaloes (Bubalus bubalis). TWAS was conducted on a test (N = 136) set of Murrah buffaloes genotyped through ddRAD sequencing. Their gene expression level was predicted using reference (N = 8) animals having both genotype and mammary epithelial cell (MEC) transcriptome information. Gene expression prediction was performed using Elastic-Net and Dirichlet Process Regression (DPR) model with fivefold cross-validation and without any cross-validation. DPR model without cross-validation predicted 80.92% of the total genes in the test group of Murrah buffaloes which was highest compared to other methods. TWAS in test individuals based on predicted gene expression, identified a significant association of one unique gene for Fat%, and two for SNF% at Bonferroni corrected threshold. The false discovery rates (FDR) corrected P-values of the top ten SNPs identified through GWAS were comparatively higher than TWAS. Gene ontology of TWAS-identified genes was performed to understand the function of these genes, it was revealed that milk production and composition genes were mainly involved in Relaxin, AMPK, and JAK-STAT signaling pathway, along with CCRI, and several key metabolic processes. The present study indicates that TWAS offers a lower false discovery rate and higher significant hits than GWAS for milk production and its composition traits. Hence, it is concluded that TWAS can be effectively used to identify genes and cis-SNPs in a population, which can be used for fabricating a low-density genomic chip for predicting milk production in Murrah buffaloes.

Evidence for early domestic yak, taurine cattle, and their hybrids on the Tibetan Plateau

Domestic yak, cattle, and their hybrids are fundamental to herder survival at high altitudes on the Tibetan Plateau. However, little is known about their history. Bos remains are uncommon in this region, and ancient domestic yak have not been securely identified. To identify Bos taxa and investigate their initial management, we conducted zooarchaeological analyses of 193 Bos specimens and sequenced five nuclear genomes from recently excavated assemblages at Bangga. Morphological data indicated that more cattle than yak were present. Ancient mitochondrial DNA and nuclear genome sequences identified taurine cattle and provided evidence for domestic yak and yak-cattle hybridization ~2500 years ago. Reliance on diverse Bos species and their hybrid has increased cattle adaptation and herder resilience to plateau conditions. Ancient cattle and yak at Bangga were closely related to contemporary livestock, indicating early herder legacies and the continuity of cattle and yak husbandry on the Tibetan Plateau.

OvaPrint-A Cell-free DNA Methylation Liquid Biopsy for the Risk Assessment of High-grade Serous Ovarian Cancer

PURPOSE

High-grade serous ovarian carcinoma (HGSOC) is the most lethal epithelial ovarian cancer (EOC) and is often diagnosed at late stage. In women with a known pelvic mass, surgery followed by pathologic assessment is the most reliable way to diagnose EOC and there are still no effective screening tools in asymptomatic women. In the current study, we developed a cell-free DNA (cfDNA) methylation liquid biopsy for the risk assessment of early-stage HGSOC.

EXPERIMENTAL DESIGN

We performed reduced representation bisulfite sequencing to identify differentially methylated regions (DMR) between HGSOC and normal ovarian and fallopian tube tissue. Next, we performed hybridization probe capture for 1,677 DMRs and constructed a classifier (OvaPrint) on an independent set of cfDNA samples to discriminate HGSOC from benign masses. We also analyzed a series of non-HGSOC EOC, including low-grade and borderline samples to assess the generalizability of OvaPrint. A total of 372 samples (tissue n = 59, plasma n = 313) were analyzed in this study.

RESULTS

OvaPrint achieved a positive predictive value of 95% and a negative predictive value of 88% for discriminating HGSOC from benign masses, surpassing other commercial tests. OvaPrint was less sensitive for non-HGSOC EOC, albeit it may have potential utility for identifying low-grade and borderline tumors with higher malignant potential.

CONCLUSIONS

OvaPrint is a highly sensitive and specific test that can be used for the risk assessment of HGSOC in symptomatic women. Prospective studies are warranted to validate OvaPrint for HGSOC and further develop it for non-HGSOC EOC histotypes in both symptomatic and asymptomatic women with adnexal masses.

PAX3-FOXO1 uses its activation domain to recruit CBP/P300 and shape RNA Pol2 cluster distribution

Activation of oncogenic gene expression from long-range enhancers is initiated by the assembly of DNA-binding transcription factors (TF), leading to recruitment of co-activators such as CBP/p300 to modify the local genomic context and facilitate RNA-Polymerase 2 (Pol2) binding. Yet, most TF-to-coactivator recruitment relationships remain unmapped. Here, studying the oncogenic fusion TF PAX3-FOXO1 (P3F) from alveolar rhabdomyosarcoma (aRMS), we show that a single cysteine in the activation domain (AD) of P3F is important for a small alpha helical coil that recruits CBP/p300 to chromatin. P3F driven transcription requires both this single cysteine and CBP/p300. Mutants of the cysteine reduce aRMS cell proliferation and induce cellular differentiation. Furthermore, we discover a profound dependence on CBP/p300 for clustering of Pol2 loops that connect P3F to its target genes. In the absence of CBP/p300, Pol2 long range enhancer loops collapse, Pol2 accumulates in CpG islands and fails to exit the gene body. These results reveal a potential novel axis for therapeutic interference with P3F in aRMS and clarify the molecular relationship of P3F and CBP/p300 in sustaining active Pol2 clusters essential for oncogenic transcription.

Transcriptomic Changes in Response to Form of Selenium on the Interferon-Tau Signaling Mechanism in the Caruncular Tissue of Beef Heifers at Maternal Recognition of Pregnancy

We have reported that selenium (Se) provided to grazing beef cattle in an inorganic (ISe) form versus a 1:1 mixture (MIX) of inorganic and organic (OSe) forms affects cholesterol biosynthesis in the corpus luteum (CL), the abundance of interferon tau (IFNτ) and progesterone (P4)-induced mRNAs in the caruncular (CAR) tissue of the endometrium, and conceptus length at maternal recognition of pregnancy (MRP). In this study, beef heifers were supplemented with a vitamin-mineral mix containing 35 ppm Se as ISe or MIX to achieve a Se-adequate status. Inseminated heifers were killed at MRP (d 17, n = 6 per treatment) for tissue collection. In CAR samples from MIX versus ISe heifers, qPCR revealed that mRNA encoding the thyroid regulating DIO2 and DIO3 was decreased (p < 0.05) and a complete transcriptomic analysis revealed effects on the interferon JAK-STAT1/2 pathway, including decreased expression of mRNAs encoding the classical interferon stimulated genes IFIT1, IFIT2, IFIT3, IRF1, IRF9, ISG15, OAS2, and RSAD2 (p < 0.05). Treatment also affected the abundance of mRNAs contributing to the immunotolerant environment (p < 0.05). In combination, these findings suggest more advanced preparation of the CAR and developing conceptus for implantation and to evade immune rejection by the maternal system in MIX- vs. ISe-treated heifers.

DNA methylation extends lifespan in the bumblebee Bombus terrestris

Epigenetic alterations are a primary hallmark of ageing. In mammals, age-related epigenetic changes alter gene expression profiles, disrupt cellular homeostasis and physiological functions and, therefore, promote ageing. It remains unclear whether ageing is also driven by epigenetic mechanisms in invertebrates. Here, we used a pharmacological hypomethylating agent (RG108) to evaluate the effects of DNA methylation (DNAme) on lifespan in an insect-the bumblebee Bombus terrestris. RG108 extended mean lifespan by 43% and induced the differential methylation of genes involved in hallmarks of ageing, including DNA damage repair and chromatin organization. Furthermore, the longevity gene sirt1 was overexpressed following the treatment. Functional experiments demonstrated that SIRT1 protein activity was positively associated with lifespan. Overall, our study indicates that epigenetic mechanisms are conserved regulators of lifespan in both vertebrates and invertebrates and provides new insights into how DNAme is involved in the ageing process in insects.

Identification of intergenerational epigenetic inheritance by whole genome DNA methylation analysis in trios

Genome-wide association studies have identified thousands of loci associated with common diseases and traits. However, a large fraction of heritability remains unexplained. Epigenetic modifications, such as the observed in DNA methylation have been proposed as a mechanism of intergenerational inheritance. To investigate the potential contribution of DNA methylation to the missing heritability, we analysed the methylomes of four healthy trios (two parents and one offspring) using whole genome bisulphite sequencing. Of the 1.5 million CpGs (19%) with over 20% variability between parents in at least one family and compatible with a Mendelian inheritance pattern, only 3488 CpGs (0.2%) lacked correlation with any SNP in the genome, marking them as potential sites for intergenerational epigenetic inheritance. These markers were distributed genome-wide, with some preference to be located in promoters. They displayed a bimodal distribution, being either fully methylated or unmethylated, and were often found at the boundaries of genomic regions with high/low GC content. This analysis provides a starting point for future investigations into the missing heritability of simple and complex traits.

Adaptive laboratory evolution for acetic acid-tolerance matches sourdough challenges with yeast phenotypes

Acetic acid tolerance of Saccharomyces cerevisiae is an important trait in sourdough fermentation processes, where the accumulation of acid by the growth of lactic acid bacteria reduces the yeast metabolic activity. In this work, we have carried out adaptive laboratory evolution (ALE) experiments in two sourdough isolates of S. cerevisiae exposed to acetic acid, or alternatively to acetic acid and myriocin, an inhibitor of sphingolipid biosynthesis that sped-up the evolutionary adaptation. Evolution approaches resulted in acetic tolerance, and surprisingly, increased lactic susceptibility. Four evolved clones, one from each parental strain and evolutionary scheme, were selected on the basis of their potential for CO2 production in sourdough conditions. Among them, two showed phenotypic instability characterized by strong lactic sensitivity after several rounds of growth under unstressed conditions, while two others, displayed increased constitutive acetic tolerance with no loss of growth in lactic medium. Genome sequencing and ploidy level analysis of all strains revealed aneuploidies, which could account for phenotypic heterogeneity. In addition, copy number variations (CNVs), affecting specially to genes involved in ion transport or flocculation, and single nucleotide polymorphisms (SNPs) were identified. Mutations in several genes, ARG82, KEX1, CTK1, SPT20, IRA2, ASG1 or GIS4, were confirmed as involved in acetic and/or lactic tolerance, and new determinants of these phenotypes, MSN5 and PSP2, identified.

Nanopore sequencing reveals methylation changes associated with obesity in circulating cell-free DNA from Göttingen Minipigs

Profiling of circulating cell-free DNA (cfDNA) by tissue-specific base modifications, such as 5-methylcytosines (5mC), may enable the monitoring of ongoing pathophysiological processes. Nanopore sequencing allows genome-wide 5mC detection in cfDNA without bisulphite conversion. The aims of this study were: i) to find differentially methylated regions (DMRs) of cfDNA associated with obesity in Göttingen minipigs using Nanopore sequencing, ii) to validate a subset of the DMRs using methylation-specific PCR (MSP-PCR), and iii) to compare the cfDNA DMRs with those from whole blood genomic DNA (gDNA). Serum cfDNA and gDNA were obtained from 10 lean and 7 obese Göttingen Minipigs both with experimentally induced myocardial infarction and sequenced using Oxford Nanopore MinION. A total of 1,236 cfDNA DMRs (FDR<0.01) were associated with obesity. In silico analysis showed enrichment of the adipocytokine signalling, glucagon signalling, and cellular glucose homoeostasis pathways. A strong cfDNA DMR was discovered in PPARGC1B, a gene linked to obesity and type 2 diabetes. The DMR was validated using MSP-PCR and correlated significantly with body weight (P < 0.05). No DMRs intersected between cfDNA and gDNA, suggesting that cfDNA originates from body-wide shedding of DNA. In conclusion, nanopore sequencing detected differential methylation in minute quantities (0.1-1 ng/µl) of cfDNA. Future work should focus on translation into human and comparing 5mC from somatic tissues to pinpoint the exact location of pathology.

The Next-Generation Oral Selective Estrogen Receptor Degrader Camizestrant (AZD9833) Suppresses ER+ Breast Cancer Growth and Overcomes Endocrine and CDK4/6 Inhibitor Resistance

Oral selective estrogen receptor degraders (SERD) could become the backbone of endocrine therapy (ET) for estrogen receptor-positive (ER+) breast cancer, as they achieve greater inhibition of ER-driven cancers than current ETs and overcome key resistance mechanisms. In this study, we evaluated the preclinical pharmacology and efficacy of the next-generation oral SERD camizestrant (AZD9833) and assessed ER-co-targeting strategies by combining camizestrant with CDK4/6 inhibitors (CDK4/6i) and PI3K/AKT/mTOR-targeted therapy in models of progression on CDK4/6i and/or ET. Camizestrant demonstrated robust and selective ER degradation, modulated ER-regulated gene expression, and induced complete ER antagonism and significant antiproliferation activity in ESR1 wild-type (ESR1wt) and mutant (ESR1m) breast cancer cell lines and patient-derived xenograft (PDX) models. Camizestrant also delivered strong antitumor activity in fulvestrant-resistant ESR1wt and ESR1m PDX models. Evaluation of camizestrant in combination with CDK4/6i (palbociclib or abemaciclib) in CDK4/6-naive and -resistant models, as well as in combination with PI3Kαi (alpelisib), mTORi (everolimus), or AKTi (capivasertib), indicated that camizestrant was active with CDK4/6i or PI3K/AKT/mTORi and that antitumor activity was further increased by the triple combination. The response was observed independently of PI3K pathway mutation status. Overall, camizestrant shows strong and broad antitumor activity in ER+ breast cancer as a monotherapy and when combined with CDK4/6i and PI3K/AKT/mTORi.

SIGNIFICANCE

Camizestrant, a next-generation oral SERD, shows promise in preclinical models of ER+ breast cancer alone and in combination with CDK4/6 and PI3K/AKT/mTOR inhibitors to address endocrine resistance, a current barrier to treatment.