An epigenetic toolbox for conservation biologists.

Ongoing climatic shifts and increasing anthropogenic pressures demand an efficient delineation of conservation units and accurate predictions of populations' resilience and adaptive potential. Molecular tools involving DNA sequencing are nowadays routinely used for these purposes. Yet, most of the existing tools focusing on sequence-level information have shortcomings in detecting signals of short-term ecological relevance. Epigenetic modifications carry valuable information to better link individuals, populations, and species to their environment. Here, we discuss a series of epigenetic monitoring tools that can be directly applied to various conservation contexts, complementing already existing molecular monitoring frameworks. Focusing on DNA sequence-based methods (e.g. DNA methylation, for which the applications are readily available), we demonstrate how (a) the identification of epi-biomarkers associated with age or infection can facilitate the determination of an individual's health status in wild populations; (b) whole epigenome analyses can identify signatures of selection linked to environmental conditions and facilitate estimating the adaptive potential of populations; and (c) epi-eDNA (epigenetic environmental DNA), an epigenetic-based conservation tool, presents a non-invasive sampling method to monitor biological information beyond the mere presence of individuals. Overall, our framework refines conservation strategies, ensuring a comprehensive understanding of species' adaptive potential and persistence on ecologically relevant timescales.

Hormesis as an adaptive response to infection.

Hormesis is a phenomenon whereby low-level stress can improve cellular, organ, or organismal fitness in response to a subsequent similar or other stress insult. Whereas hormesis is thought to contribute to the fitness benefits arising from symbiotic host-microbe interactions, the putative benefits of hormesis in host-pathogen interactions have yet to be explored. Hormetic responses have nonetheless been reported in experimental models of infection, a common feature of which is regulation of host mitochondrial function. We propose that these mitohormetic responses could be harnessed therapeutically to limit the severity of infectious diseases.

Genomic signatures of hybridization between Ixodes ricinus and Ixodes persulcatus in natural populations.

Identifying hybridization between common pathogen vectors is essential due to the major public health implications through risks associated with hybrid's enhanced pathogen transmission potential. The hard-ticks Ixodes ricinus and Ixodes persulcatus are the two most common vectors of tick-borne pathogens that affect human and animal health in Europe. Ixodes ricinus is a known native species in Finland with a well-known distribution, whereas I. persulcatus has expanded in range and abundance over the past 60 years, and currently it appears the most common tick species in certain areas in Finland. Here we used double-digest restriction site-associated DNA (ddRAD) sequencing on 186 ticks (morphologically identified as 92 I. ricinus, and 94 I. persulcatus) collected across Finland to investigate whether RAD generated single nucleotide polymorphisms (SNPs) can discriminate tick species and identify potential hybridization events. Two different clustering methods were used to assign specific species based on how they clustered and identified hybrids among them. We were able to discriminate between the two tick species and identified 11 putative hybrids with admixed genomic proportions ranging from approximately 24 to 76 percent. Four of these hybrids were morphologically identified as I. ricinus while the remaining seven were identified as I. persulcatus. Our results thus indicate that RAD SNPs are robust in identifying both species of the ticks as well as putative hybrids. These results further suggest ongoing hybridization between I. ricinus and I. persulcatus in their natural populations in Finland. The unique ability of RAD markers to discriminate between tick species and hybrids adds a useful aspect to tick evolutionary studies. Our findings align with previous studies and suggest a shared evolutionary history between the species, with instances of individuals possessing a considerable proportion of the other species' genome. This study is a significant step in understanding the formation of hybridization zones due to range expansion potentially associated with climate change.

Ferritin heavy chain supports stability and function of the regulatory T cell lineage.

Regulatory T (TREG) cells develop via a program orchestrated by the transcription factor forkhead box protein P3 (FOXP3). Maintenance of the TREG cell lineage relies on sustained FOXP3 transcription via a mechanism involving demethylation of cytosine-phosphate-guanine (CpG)-rich elements at conserved non-coding sequences (CNS) in the FOXP3 locus. This cytosine demethylation is catalyzed by the ten-eleven translocation (TET) family of dioxygenases, and it involves a redox reaction that uses iron (Fe) as an essential cofactor. Here, we establish that human and mouse TREG cells express Fe-regulatory genes, including that encoding ferritin heavy chain (FTH), at relatively high levels compared to conventional T helper cells. We show that FTH expression in TREG cells is essential for immune homeostasis. Mechanistically, FTH supports TET-catalyzed demethylation of CpG-rich sequences CNS1 and 2 in the FOXP3 locus, thereby promoting FOXP3 transcription and TREG cell stability. This process, which is essential for TREG lineage stability and function, limits the severity of autoimmune neuroinflammation and infectious diseases, and favors tumor progression. These findings suggest that the regulation of intracellular iron by FTH is a stable property of TREG cells that supports immune homeostasis and limits the pathological outcomes of immune-mediated inflammation.

Targeting circulating labile heme as a defense strategy against malaria.

Severe presentations of malaria emerge as Plasmodium (P.) spp. parasites invade and lyse red blood cells (RBC), producing extracellular hemoglobin (HB), from which labile heme is released. Here, we tested whether scavenging of extracellular HB and/or labile heme, by haptoglobin (HP) and/or hemopexin (HPX), respectively, counter the pathogenesis of severe presentations of malaria. We found that circulating labile heme is an independent risk factor for cerebral and non-cerebral presentations of severe P. falciparum malaria in children. Labile heme was negatively correlated with circulating HP and HPX, which were, however, not risk factors for severe P. falciparum malaria. Genetic Hp and/or Hpx deletion in mice led to labile heme accumulation in plasma and kidneys, upon Plasmodium infection This was associated with higher incidence of mortality and acute kidney injury (AKI) in ageing but not adult Plasmodium-infected mice, and was corroborated by an inverse correlation between heme and HPX with serological markers of AKI in P. falciparum malaria. In conclusion, HP and HPX act in an age-dependent manner to prevent the pathogenesis of severe presentation of malaria in mice and presumably in humans.

Reframing sepsis immunobiology for translation: towards informative subtyping and targeted immunomodulatory therapies.

Sepsis is a common and deadly condition. Within the current model of sepsis immunobiology, the framing of dysregulated host immune responses into proinflammatory and immunosuppressive responses for the testing of novel treatments has not resulted in successful immunomodulatory therapies. Thus, the recent focus has been to parse observable heterogeneity into subtypes of sepsis to enable personalised immunomodulation. In this Personal View, we highlight that many fundamental immunological concepts such as resistance, disease tolerance, resilience, resolution, and repair are not incorporated into the current sepsis immunobiology model. The focus for addressing heterogeneity in sepsis should be broadened beyond subtyping to encompass the identification of deterministic molecular networks or dominant mechanisms. We explicitly reframe the dysregulated host immune responses in sepsis as altered homoeostasis with pathological disruption of immune-driven resistance, disease tolerance, resilience, and resolution mechanisms. Our proposal highlights opportunities to identify novel treatment targets and could enable successful immunomodulation in the future.

A Review about the Mycoremediation of Soil Impacted by War-like Activities: Challenges and Gaps.

(1) Background: The frequency and intensity of war-like activities (war, military training, and shooting ranges) worldwide cause soil pollution by metals, metalloids, explosives, radionuclides, and herbicides. Despite this environmentally worrying scenario, soil decontamination in former war zones almost always involves incineration. Nevertheless, this practice is expensive, and its efficiency is suitable only for organic pollutants. Therefore, treating soils polluted by wars requires efficient and economically viable alternatives. In this sense, this manuscript reviews the status and knowledge gaps of mycoremediation. (2) Methods: The literature review consisted of searches on ScienceDirect and Web of Science for articles (1980 to 2023) on the mycoremediation of soils containing pollutants derived from war-like activities. (3) Results: This review highlighted that mycoremediation has many successful applications for removing all pollutants of war-like activities. However, the mycoremediation of soils in former war zones and those impacted by military training and shooting ranges is still very incipient, with most applications emphasizing explosives. (4) Conclusion: The mycoremediation of soils from conflict zones is an entirely open field of research, and the main challenge is to optimize experimental conditions on a field scale.

Genomic basis of melanin-associated phenotypes suggests colour-specific environmental adaptations in tawny owls.

Feathers comprise a series of evolutionary innovations but also harbour colour, a key biological trait known to co-vary with life history or complex traits. Those relationships are particularly true in melanin-based pigmentation species due to known pleiotropic effects of the melanocortin pathway - originating from melanin-associated phenotypes. Here, we explore the molecular basis of melanin colouration and expected co-variation at the molecular level in the melanin-based, colour polymorphic system of the tawny owl (Strix aluco). An extensive body of literature has revealed that grey and brown tawny owl colour morphs differ in a series of life history and behavioural traits. Thus, it is plausible to expect co-variation also at molecular level between colour morphs. To investigate this possibility, we assembled the first draft genome of the species against which we mapped ddRADseq reads from 220 grey and 150 brown morphs - representing 10 years of pedigree data from a population in Southern Finland - and explored genome-wide associations with colour phenotype. Our results revealed putative molecular signatures of cold adaptation strongly associated with the grey phenotype, namely, a non-synonymous substitution in MCHR1, plus 2 substitutions in non-coding regions of FTCD and FAM135A whose genotype combinations obtained a predictive power of up to 100% (predicting grey colour). These suggest a molecular basis of cold environment adaptations predicted to be grey-morph specific. Our results potentially reveal part of the molecular machinery of melanin-associated phenotypes and provide novel insights towards understanding the functional genomics of colour polymorphism in melanin-based pigmented species.

Detection of Brain Metastases by Contrast-Enhanced MRI: Comparison of Gadopiclenol and Gadobenate in a Mouse Model.

OBJECTIVES: The aim of this study was to evaluate the capacity of gadopiclenol, a high-relaxivity gadolinium-based contrast agent to detect brain metastases in mice as a function of dose (0.08 mmol/kg or 0.1 mmol/kg) compared with gadobenate at 0.1 mmol/kg. MATERIALS AND METHODS: Brain metastases were induced by ultrasound-guided intracardiac implantation of 1.10 5 MDA-MB-231Br cells in the left ventricle of 18 anesthetized Balb/c Nude nu/nu female mice. At day 28 ± 3 after cell injection, each mouse received 2 crossover intravenous injections at 24-hour intervals, randomly selected from 2 doses of gadopiclenol (0.08 mmol/kg or 0.1 mmol/kg) and gadobenate (0.1 mmol/kg) with n = 6 mice/group (3 groups). Brain magnetic resonance imaging sessions were performed at 4 weeks on a 2.35 T magnet with a 3-dimensional T1-weighted high-resolution gradient echo sequence, before and after each injection. Images were blindly and randomly analyzed to detect enhancing lesions. Contrast-to-noise ratio between the metastases and the surrounding healthy parenchyma was calculated, based on region-of-interest signal measurements. In 2 animals per group, an early time point was added to the protocol (day 22 ± 3) to evaluate the sensitivity of detection as a function of time. After the last imaging session, the presence and location of whole-brain metastases were confirmed by histology in 4 mice. RESULTS: After gadopiclenol, approximately twice as many metastases were detected compared with gadobenate, regardless of the dose. Contrast-to-noise ratios of the detected metastases were 2.3 and 3.3 times higher with gadopiclenol at 0.08 mmol/kg and 0.1 mmol/kg, respectively, compared with gadobenate at 0.1 mmol/kg ( P < 0.0001). Gadopiclenol at the dose of 0.1 mmol/kg resulted in a 1.4-fold higher contrast compared with gadopiclenol at 0.08 mmol/kg ( P < 0.02). In a subset of mice that were imaged 1 week earlier, 2 metastases were detected with gadopiclenol and not with gadobenate. CONCLUSIONS: The high-relaxivity macrocyclic gadolinium-based contrast agent gadopiclenol allowed higher diagnostic performance for detecting brain enhancing metastases in terms of contrast-to-noise ratio and number of detected metastases compared with gadobenate, at both equal (0.1 mmol/kg) dose and 20% lower Gd dose (0.08 mmol/kg). Tumor detection was higher after gadopiclenol at the dose of 0.1 mmol/kg compared with 0.08 mmol/kg.

Biological efficacy of simulated radiolabeled Lipiodol® ultra-fluid and microspheres for various beta emitters: study based on VX2 tumors.

BACKGROUND: Radioembolization is one therapeutic option for the treatment of locally early-stage hepatocellular carcinoma. The aim of this study was to evaluate the distribution of Lipiodol® ultra-fluid and microspheres and to simulate their effectiveness with different beta emitters (90Y, 188Re, 32P, 166Ho, 131I, and 177Lu) on VX2 tumors implanted in the liver of 30 New Zealand rabbits. RESULTS: Twenty-three out of 30 rabbits had exploitable data: 14 in the group that received Lipiodol® ultra-fluid (group L), 6 in the group that received microspheres (group M), and 3 in the control group (group C). The histologic analysis showed that the Lipiodol® ultra-fluid distributes homogeneously in the tumor up to 12 days after injection. The X-ray µCT images showed that Lipiodol® ultra-fluid has a more distal penetration in the tumor than microspheres. The entropy (disorder of the system) in the L group was significantly higher than in the M group (4.06 vs 2.67, p = 0.01). Equivalent uniform biological effective doses (EUBED) for a tumor-absorbed dose of 100 Gy were greater in the L group but without statistical significance except for 177Lu (p = 0.03). The radionuclides ranking by EUBED (from high to low) was 90Y, 188Re, 32P, 166Ho, 131I, and 177Lu. CONCLUSIONS: This study showed a higher ability of Lipiodol® ultra-fluid to penetrate the tumor that translated into a higher EUBED. This study confirms 90Y as a good candidate for radioembolization, although 32P, 166Ho, and 188Re can achieve similar results.

Impact of conservation tillage on wheat performance and its microbiome.

Winter wheat is an important cereal consumed worldwide. However, current management practices involving chemical fertilizers, irrigation, and intensive tillage may have negative impacts on the environment. Conservation agriculture is often presented as a sustainable alternative to maintain wheat production, favoring the beneficial microbiome. Here, we evaluated the impact of different water regimes (rainfed and irrigated), fertilization levels (half and full fertilization), and tillage practices (occasional tillage and no-tillage) on wheat performance, microbial activity, and rhizosphere- and root-associated microbial communities of four winter wheat genotypes (Antequera, Allez-y, Apache, and Cellule) grown in a field experiment. Wheat performance (i.e., yield, plant nitrogen concentrations, and total nitrogen uptake) was mainly affected by irrigation, fertilization, and genotype, whereas microbial activity (i.e., protease and alkaline phosphatase activities) was affected by irrigation. Amplicon sequencing data revealed that habitat (rhizosphere vs. root) was the main factor shaping microbial communities and confirmed that the selection of endophytic microbial communities takes place thanks to specific plant-microbiome interactions. Among the experimental factors applied, the interaction of irrigation and tillage influenced rhizosphere- and root-associated microbiomes. The findings presented in this work make it possible to link agricultural practices to microbial communities, paving the way for better monitoring of these microorganisms in the context of agroecosystem sustainability.

Single-cell RNA sequencing and analysis of rodent blood stage Plasmodium.

Bulk RNA sequencing of Plasmodium spp., the causative parasite of malaria, fails to discriminate developmental-stage-specific gene regulation. Here, we provide a protocol that uses single-cell RNA sequencing of FACS-sorted Plasmodium-chabaudi-chabaudi-AS-infected red blood cells (iRBCs) to characterize developmental-stage-specific modulation of gene expression during malaria blood stage. We describe steps for infecting mice, monitoring disease progression, preparing iRBCs, and single-cell sequencing iRBCs. We then detail procedures for analyzing scRNA-seq data. For complete details on the use and execution of this protocol, please refer to Ramos et al.1.

Seascape genomics reveals limited dispersal and suggests spatially varying selection among European populations of sea lamprey (Petromyzon marinus).

Sea lamprey Petromyzon marinus is an anadromous and semelparous fish without homing behaviors. Despite being a freshwater, free-living organism for a large part of their life cycle, its adulthood is spent as a parasite of marine vertebrates. In their native European range, while it is well-established that sea lampreys comprise a single nearly-panmictic population, few studies have further explored the evolutionary history of natural populations. Here, we performed the first genome-wide characterization of sea lamprey's genetic diversity in their European natural range. The objectives were to investigate the connectivity among river basins and explore evolutionary processes mediating dispersal during the marine phase, with the sequencing of 186 individuals from 8 locations spanning the North Eastern Atlantic coast and the North Sea with double-digest RAD-sequencing, obtaining a total of 30,910 bi-allelic SNPs. Population genetic analyses reinforced the existence of a single metapopulation encompassing freshwater spawning sites within the North Eastern Atlantic and the North Sea, though the prevalence of private alleles at northern latitudes suggested some limits to the species' dispersal. Seascape genomics suggested a scenario where oxygen concentration and river runoffs impose spatially varying selection across their distribution range. Exploring associations with the abundance of potential hosts further suggested that hake and cod could also impose selective pressures, although the nature of such putative biotic interactions was unresolved. Overall, the identification of adaptive seascapes in a panmictic anadromous species could contribute to conservation practices by providing information for restoration activities to mitigate local extinctions on freshwater sites.

Renal control of life-threatening malarial anemia.

Iron recycling prevents the development of anemia under homeostatic conditions. Whether iron recycling was co-opted as a defense strategy to prevent the development of anemia in response to infection is unclear. We find that in severe Plasmodium falciparum malaria, the onset of life-threatening anemia is associated with acute kidney injury (AKI), irrespective of parasite load. Using a well-established experimental rodent model of malaria anemia, we identify a transcriptional response that endows renal proximal tubule epithelial cells (RPTECs) with the capacity to store and recycle iron during P. chabaudi chabaudi (Pcc) infection. This response encompasses the induction of ferroportin 1/SLC40A1, which exports iron from RPTECs and counteracts AKI while supporting compensatory erythropoiesis and preventing the onset of life-threatening malarial anemia. Iron recycling by myeloid cells is dispensable to this protective response, suggesting that RPTECs provide an iron-recycling salvage pathway that prevents the pathogenesis of life-threatening malarial anemia.

Bringing to light nuclear-mitochondrial insertions in the genomes of nocturnal predatory birds.

Mito-nuclear insertions, or NUMTs, relate to genetic material of mitochondrial origin that have been transferred to the nuclear DNA molecule. The increasing amounts of genomic data currently being produced presents an opportunity to investigate this type of patterns in genome evolution of non-model organisms. Identifying NUMTs across a range of closely related taxa allows one to generalize patterns of insertion and maintenance in autosomes, which is ultimately relevant to the understanding of genome biology and evolution. Here we collected existing pairwise genome-mitogenome data of the order Strigiformes, a group that includes all the nocturnal bird predators. We identified NUMTs by applying percent similarity thresholds after blasting mitochondrial genomes against nuclear genome assemblies. We identified NUMTsin all genomes with numbers ranging from 4 in Bubo bubo to 24 in Ciccaba nigrolineata. Statistical analyses revealed NUMT size to negatively correlate with NUMT's sequence similarity to with original mtDNA region. Lastly, characterizing these nuclear insertions of mitochondrial origin in a comparative genomics framework produced variable phylogenetic patterns, suggesting in some cases that insertions might pre-date speciation events within Strigiformes.

Antitumoral effect of local injection of TLR-9 agonist emulsified in Lipiodol with systemic anti-PD-1 in a murine model of colorectal carcinoma.

Introduction: Local treatments of cancer, including transarterial chemoembolization, could enhance responses to systemic immune checkpoint inhibitors such as anti-PD-1 antibodies. Lipiodol, a radiopaque oil, is widely used for transarterial chemoembolization as a tumor-targeting drug carrier and could be used in emulsion with immunomodulators. This study aimed at evaluating the antitumoral effect of intra-tumoral injection of Lipiodol-immunomodulator emulsions combined with systemic anti-PD-1 therapy in a murine model of colorectal carcinoma. Method: Mice (male BALB/c) with anti-PD-1-resistant subcutaneous CT26 tumors were injected with immunomodulators, emulsified or not with Lipiodol (N=10-12/group). Results: The TLR-9 agonist CpG displayed antitumor effects, while Poly I:C and QS21 did not. The Lipiodol-CpG emulsion appeared to be stable and maintained CpG within tumors for a longer time. Repeated intra-tumoral injections, combined with anti-PD-1, induced responses towards the tumor as well as to a distant metastatic-like nodule. This treatment was associated with an increase in proliferative CD8+ T cells and of IFN-γ expression, a decrease in proliferative regulatory T cells but also, surprisingly, an increase in myeloid derived suppressor cells. Conclusions: Local administration of CpG emulsified with Lipiodol led to an effective antitumoral effect when combined to systemic anti-PD-1 therapy. Lipiodol, apart from its radiopaque properties, is an efficient drug-delivery system. The formulated oil-in-water emulsion allows efficient loading and control release of CpG, which induces favorable immune modifications in this murine tumor model.

DNA damage independent inhibition of NF-κB transcription by anthracyclines.

Anthracyclines are among the most used and effective anticancer drugs. Their activity has been attributed to DNA double-strand breaks resulting from topoisomerase II poisoning and to eviction of histones from select sites in the genome. Here, we show that the extensively used anthracyclines Doxorubicin, Daunorubicin, and Epirubicin decrease the transcription of nuclear factor kappa B (NF-κB)-dependent gene targets, but not interferon-responsive genes in primary mouse (Mus musculus) macrophages. Using an NMR-based structural approach, we demonstrate that anthracyclines disturb the complexes formed between the NF-κB subunit RelA and its DNA-binding sites. The anthracycline variants Aclarubicin, Doxorubicinone, and the newly developed Dimethyl-doxorubicin, which share anticancer properties with the other anthracyclines but do not induce DNA damage, also suppressed inflammation, thus uncoupling DNA damage from the effects on inflammation. These findings have implications for anticancer therapy and for the development of novel anti-inflammatory drugs with limited side effects for life-threatening conditions such as sepsis.

Predictive value of common genetic variants in idiopathic pulmonary fibrosis survival.

Idiopathic pulmonary fibrosis (IPF) is a chronic, progressive, fibrosing interstitial pneumonia of unknown etiology. The role of genetic risk factors has been the focus of numerous studies probing for associations of genetic variants with IPF. We aimed to determine whether single-nucleotide polymorphisms (SNPs) of four candidate genes are associated with IPF susceptibility and survival in a Portuguese population. A retrospective case-control study was performed with 64 IPF patients and 74 healthy controls. Ten single-nucleotide variants residing in the MUC5B, TOLLIP, SERPINB1, and PLAU genes were analyzed. Single- and multi-locus analyses were performed to investigate the predictive potential of specific variants in IPF susceptibility and survival. Multifactor dimensionality reduction (MDR) was employed to uncover predictive multi-locus interactions underlying IPF susceptibility. The MUC5B rs35705950 SNP was significantly associated with IPF: T allele carriers were significantly more frequent among IPF patients (75.0% vs 20.3%, P < 1.0 × 10-6). Genotypic and allelic distributions of TOLLIP, PLAU, and SERPINB1 SNPs did not differ significantly between groups. However, the MUC5B-TOLLIP T-C-T-C haplotype, defined by the rs35705950-rs111521887-rs5743894-rs5743854 block, emerged as an independent protective factor in IPF survival (HR = 0.37, 95% CI 0.17-0.78, P = 0.009, after adjustment for FVC). No significant multi-locus interactions correlating with disease susceptibility were detected. MUC5B rs35705950 was linked to an increased risk for IPF, as reported for other populations, but not to disease survival. A haplotype incorporating SNPs of the MUC5B-TOLLIP locus at 11p15.5 seems to predict better survival and could prove useful for prognostic purposes and IPF patient stratification. KEY MESSAGES : The MUC5B rs35705950 minor allele is associated with IPF risk in the Portuguese. No predictive multi-locus interactions of IPF susceptibility were identified by MDR. A haplotype defined by MUC5B and TOLLIP SNPs is a protective factor in IPF survival. The haplotype may be used as a prognostic tool for IPF patient stratification.

Fungal Communities of the Pine Wilt Disease Complex: Studying the Interaction of Ophiostomatales With Bursaphelenchus xylophilus.

Considered one of the most devastating plant-parasitic nematodes worldwide, Bursaphelenchus xylophilus (commonly known as pinewood nematode, PWN) is the causal agent of the pine wilt disease in the Eurasian coniferous forests. This migratory parasitic nematode is carried by an insect vector (Monochamus spp.) into the host tree (Pinus species), where it can feed on parenchymal cells and reproduce massively, resulting in the tree wilting. In declining trees, PWN populations are strongly dependent on fungal communities colonizing the host (predominantly ophiostomatoid fungi known to cause sapwood blue-staining, the blue-stain fungi), which not only influence their development and life cycle but also the number of individuals carried by the insect vector into a new host. Our main aim is to understand if PWN-associated mycobiota plays a key role in the development of PWD, in interaction with the PWN and the insect vector, and to what extent it can be targeted to disrupt the disease cycle. For this purpose, we characterized the fungal communities of Pinus pinaster trees infected and non-infected with PWN in three collection sites in Continental Portugal with different PWD temporal incidences. Our results showed that non-infected P. pinaster mycoflora is more diverse (in terms of abundance and fungal richness) than PWN-infected pine trees in the most recent PWD foci, as opposed to the fungal communities of long-term PWD history sites. Then, due to their ecological importance for PWN survival, representatives of the main ophiostomatoid fungi isolated (Ophiostoma, Leptographium, and Graphilbum) were characterized for their adaptative response to temperature, competition in-between taxa, and as food source for PWN. Under the conditions studied, Leptographium isolates showed promising results for PWN control. They could outcompete the other species, especially O. ips, and significantly reduce the development of PWN populations when compared to Botrytis cinerea (routinely used for PWN lab culturing), suggesting this to be a natural antagonist not only for the other blue-stain species but also for the PWN.