Large registry-based analysis of genetic predisposition to tuberculosis identifies genetic risk factors at HLA.

Tuberculosis is a significant public health concern resulting in the death of over 1 million individuals each year worldwide. While treatment options and vaccines exist, a substantial number of infections still remain untreated or are caused by treatment resistant strains. Therefore, it is important to identify mechanisms that contribute to risk and prognosis of tuberculosis as this may provide tools to understand disease mechanisms and provide novel treatment options for those with severe infection. Our goal was to identify genetic risk factors that contribute to the risk of tuberculosis and to understand biological mechanisms and causality behind the risk of tuberculosis. A total of 1895 individuals in the FinnGen study had International Classification of Diseases-based tuberculosis diagnosis. Genome-wide association study analysis identified genetic variants with statistically significant association with tuberculosis at the human leukocyte antigen (HLA) region (P < 5e-8). Fine mapping of the HLA association provided evidence for one protective haplotype tagged by HLA DQB1*05:01 (P = 1.82E-06, OR = 0.81 [CI 95% 0.74-0.88]), and predisposing alleles tagged by HLA DRB1*13:02 (P = 0.00011, OR = 1.35 [CI 95% 1.16-1.57]). Furthermore, genetic correlation analysis showed association with earlier reported risk factors including smoking (P < 0.05). Mendelian randomization supported smoking as a risk factor for tuberculosis (inverse-variance weighted P < 0.05, OR = 1.83 [CI 95% 1.15-2.93]) with no significant evidence of pleiotropy. Our findings indicate that specific HLA alleles associate with the risk of tuberculosis. In addition, lifestyle risk factors such as smoking contribute to the risk of developing tuberculosis.

Genome-wide association studies highlight novel risk loci for septal defects and left-sided congenital heart defects.

BACKGROUND: Congenital heart defects (CHD) are structural defects of the heart affecting approximately 1% of newborns. They exhibit low penetrance and non-Mendelian patterns of inheritance as varied and complex traits. While genetic factors are known to play an important role in the development of CHD, the specific genetics remain unknown for the majority of patients. To elucidate the underlying genetic risk, we performed a genome wide association study (GWAS) of CHDs in general and specific CHD subgroups using the FinnGen Release 10 (R10) (N > 393,000), followed by functional fine-mapping through eQTL and co-localization analyses using the GTEx database. RESULTS: We discovered three genome-wide significant loci associated with general CHD. Two of them were located in chromosome 17: 17q21.32 (rs2316327, intronic: LRRC37A2, Odds ratio (OR) [95% Confidence Interval (CI)] = 1.17[1.12-1.23], p = 1.5 × 10-9) and 17q25.3 (rs1293973611, nearest: BAHCC1, OR[95%CI] = 4.48[2.80-7.17], p = 7.0 × 10-10), respectively, and in addition to general CHD, the rs1293973611 locus was associated with the septal defect subtype. The third locus was in band 1p21.2 (rs35046143, nearest: PALMD, OR[95%CI] = 1.15[1.09-1.21], p = 7.1 × 10-9), and it was associated with general CHD and left-sided lesions. In the subgroup analysis, two additional loci were associated with septal defects (rs75230966 and rs6824295), and one with left-sided lesions (rs1305393195). In the eQTL analysis the variants rs2316327 (general CHD), and rs75230966 (septal defects) both located in 17q21.32 (with a LD r2 of 0.41) were both predicted to significantly associate with the expression of WNT9B in the atrial appendage tissue category. This effect was further confirmed by co-localization analysis, which also implicated WNT3 expression in the atrial appendage. A meta-analysis of general CHD together with the UK Biobank (combined N = 881,678) provided a different genome-wide significant locus in LRRC37A2; rs16941382 (OR[95%CI] = 1.15[1.11-1.20], p = 1.5 × 10-9) which is in significant LD with rs2316327. CONCLUSIONS: Our results of general CHD and different CHD subcategories identified a complex risk locus on chromosome 17 near BAHCC1 and LRRC37A2, interacting with the genes WNT9B, WNT3 and MYL4, may constitute potential novel CHD risk associated loci, warranting future experimental tests to determine their role.

Host-microbiota-insect interactions drive emergent virulence in a complex tree disease.

Forest declines caused by climate disturbance, insect pests and microbial pathogens threaten the global landscape, and tree diseases are increasingly attributed to the emergent properties of complex ecological interactions between the host, microbiota and insects. To address this hypothesis, we combined reductionist approaches (single and polyspecies bacterial cultures) with emergentist approaches (bacterial inoculations in an oak infection model with the addition of insect larvae) to unravel the gene expression landscape and symptom severity of host-microbiota-insect interactions in the acute oak decline (AOD) pathosystem. AOD is a complex decline disease characterized by predisposing abiotic factors, inner bark lesions driven by a bacterial pathobiome, and larval galleries of the bark-boring beetle Agrilus biguttatus. We identified expression of key pathogenicity genes in Brenneria goodwinii, the dominant member of the AOD pathobiome, tissue-specific gene expression profiles, cooperation with other bacterial pathobiome members in sugar catabolism, and demonstrated amplification of pathogenic gene expression in the presence of Agrilus larvae. This study highlights the emergent properties of complex host-pathobiota-insect interactions that underlie the pathology of diseases that threaten global forest biomes.

Short oligogalacturonides induce pathogen resistance-associated gene expression in Arabidopsis thaliana.

BACKGROUND: Oligogalacturonides (OGs) are important components of damage-associated molecular pattern (DAMP) signaling and influence growth regulation in plants. Recent studies have focused on the impact of long OGs (degree of polymerization (DP) from 10-15), demonstrating the induction of plant defense signaling resulting in enhanced defenses to necrotrophic pathogens. To clarify the role of trimers (trimeric OGs, DP3) in DAMP signaling and their impact on plant growth regulation, we performed a transcriptomic analysis through the RNA sequencing of Arabidopsis thaliana exposed to trimers. RESULTS: The transcriptomic data from trimer-treated Arabidopsis seedlings indicate a clear activation of genes involved in defense signaling, phytohormone signaling and a down-regulation of genes involved in processes related to growth regulation and development. This is further accompanied with improved defenses against necrotrophic pathogens triggered by the trimer treatment, indicating that short OGs have a clear impact on plant responses, similar to those described for long OGs. CONCLUSIONS: Our results demonstrate that trimers are indeed active elicitors of plant defenses. This is clearly indicated by the up-regulation of genes associated with plant defense signaling, accompanied with improved defenses against necrotrophic pathogens. Moreover, trimers simultaneously trigger a clear down-regulation of genes and gene sets associated with growth and development, leading to stunted seedling growth in Arabidopsis.

Revised phylogeny and novel horizontally acquired virulence determinants of the model soft rot phytopathogen Pectobacterium wasabiae SCC3193.

Soft rot disease is economically one of the most devastating bacterial diseases affecting plants worldwide. In this study, we present novel insights into the phylogeny and virulence of the soft rot model Pectobacterium sp. SCC3193, which was isolated from a diseased potato stem in Finland in the early 1980s. Genomic approaches, including proteome and genome comparisons of all sequenced soft rot bacteria, revealed that SCC3193, previously included in the species Pectobacterium carotovorum, can now be more accurately classified as Pectobacterium wasabiae. Together with the recently revised phylogeny of a few P. carotovorum strains and an increasing number of studies on P. wasabiae, our work indicates that P. wasabiae has been unnoticed but present in potato fields worldwide. A combination of genomic approaches and in planta experiments identified features that separate SCC3193 and other P. wasabiae strains from the rest of soft rot bacteria, such as the absence of a type III secretion system that contributes to virulence of other soft rot species. Experimentally established virulence determinants include the putative transcriptional regulator SirB, two partially redundant type VI secretion systems and two horizontally acquired clusters (Vic1 and Vic2), which contain predicted virulence genes. Genome comparison also revealed other interesting traits that may be related to life in planta or other specific environmental conditions. These traits include a predicted benzoic acid/salicylic acid carboxyl methyltransferase of eukaryotic origin. The novelties found in this work indicate that soft rot bacteria have a reservoir of unknown traits that may be utilized in the poorly understood latent stage in planta. The genomic approaches and the comparison of the model strain SCC3193 to other sequenced Pectobacterium strains, including the type strain of P. wasabiae, provides a solid basis for further investigation of the virulence, distribution and phylogeny of soft rot bacteria and, potentially, other bacteria as well.

Genetics of sleep medication purchases suggests causality from sleep problems to psychiatric traits.

STUDY OBJECTIVES: Over 10% of the population in Europe and in the United States use sleep medication to manage sleep problems. Our objective was to elucidate genetic risk factors and clinical correlates that contribute to sleep medication purchase and estimate the comorbid impact of sleep problems. METHODS: We performed epidemiological analysis for psychiatric diagnoses, and genetic association studies of sleep medication purchase in 797 714 individuals from FinnGen Release 7 (N = 311 892) and from the UK Biobank (N = 485 822). Post-association analyses included genetic correlation, co-localization, Mendelian randomization (MR), and polygenic risk estimation. RESULTS: In a GWAS we identified 27 genetic loci significantly associated with sleep medication, located in genes associated with sleep; AUTS2, CACNA1C, MEIS1, KIRREL3, PAX8, GABRA2, psychiatric traits; CACNA1C, HIST1H2BD, NUDT12. TOPAZ1 and TSNARE1. Co-localization and expression analysis emphasized effects on the KPNA2, GABRA2, and CACNA1C expression in the brain. Sleep medications use was epidemiologically related to psychiatric traits in FinnGen (OR [95% (CI)] = 3.86 [3.78 to 3.94], p < 2 × 10-16), and the association was accentuated by genetic correlation and MR; depression (rg = 0.55 (0.027), p = 2.86 × 10-89, p MR = 4.5 × 10-5), schizophrenia (rg = 0.25 (0.026), p = 2.52 × 10-21, p MR = 2 × 10-4), and anxiety (rg = 0.44 (0.047), p = 2.88 × 10-27, p MR = 8.6 × 10-12). CONCLUSIONS: These results demonstrate the genetics behind sleep problems and the association between sleep problems and psychiatric traits. Our results highlight the scientific basis for sleep management in treating the impact of psychiatric diseases.

Nightmares share genetic risk factors with sleep and psychiatric traits.

Nightmares are vivid, extended, and emotionally negative or negative dreams that awaken the dreamer. While sporadic nightmares and bad dreams are common and generally harmless, frequent nightmares often reflect underlying pathologies of emotional regulation. Indeed, insomnia, depression, anxiety, or alcohol use have been associated with nightmares in epidemiological and clinical studies. However, the connection between nightmares and their comorbidities are poorly understood. Our goal was to examine the genetic risk factors for nightmares and estimate correlation or causality between nightmares and comorbidities. We performed a genome-wide association study (GWAS) in 45,255 individuals using a questionnaire-based assessment on the frequency of nightmares during the past month and genome-wide genotyping data. While the GWAS did not reveal individual risk variants, heritability was estimated at 5%. In addition, the genetic correlation analysis showed a robust correlation (rg > 0.4) of nightmares with anxiety (rg = 0.671, p = 7.507e-06), depressive (rg = 0.562, p = 1.282e-07) and posttraumatic stress disorders (rg = 0.4083, p = 0.0152), and personality trait neuroticism (rg = 0.667, p = 4.516e-07). Furthermore, Mendelian randomization suggested causality from insomnia to nightmares (beta = 0.027, p = 0.0002). Our findings suggest that nightmares share genetic background with psychiatric traits and that insomnia may increase an individual's liability to experience frequent nightmares. Given the significant correlations with psychiatric and psychological traits, it is essential to grow awareness of how nightmares affect health and disease and systematically collect information about nightmares, especially from clinical samples and larger cohorts.

Genetic Analysis of Obstructive Sleep Apnea and Its Relationship with Severe COVID-19.

Rationale: Although patients with obstructive sleep apnea (OSA) have a higher risk for coronavirus disease (COVID-19) hospitalization, the causal relationship has remained unexplored. Objectives: To understand the causal relationship between OSA and COVID-19 by leveraging data from vaccination and electronic health records, genetic risk factors from genome-wide association studies, and Mendelian randomization. Methods: We elucidated genetic risk factors for OSA using FinnGen (total N = 377,277), performing genome-wide association. We used the associated variants as instruments for univariate and multivariate Mendelian randomization (MR) analyses and computed absolute risk reduction against COVID-19 hospitalization with or without vaccination. Results: We identified nine novel loci for OSA and replicated our findings in the Million Veteran Program. Furthermore, MR analysis showed that OSA was a causal risk factor for severe COVID-19 (P = 9.41 × 10-4). Probabilistic modeling showed that the strongest genetic risk factor for OSA at the FTO locus reflected a signal of higher body mass index (BMI), whereas BMI-independent association was seen with the earlier reported SLC9A4 locus and a MECOM locus, which is a transcriptional regulator with 210-fold enrichment in the Finnish population. Similarly, multivariate MR analysis showed that the causality for severe COVID-19 was driven by BMI (multivariate MR P = 5.97 × 10-6, ß = 0.47). Finally, vaccination reduced the risk for COVID-19 hospitalization more in the patients with OSA than in the non-OSA controls, with respective absolute risk reductions of 13.3% versus 6.3%. Conclusions: Our analysis identified novel genetic risk factors for OSA and showed that OSA is a causal risk factor for severe COVID-19. The effect is predominantly explained by higher BMI and suggests BMI-dependent effects at the level of individual variants and at the level of comorbid causality.

The public health impact of poor sleep on severe COVID-19, influenza and upper respiratory infections.

BACKGROUND: Poor sleep is associated with an increased risk of infections and all-cause mortality but the causal direction between poor sleep and respiratory infections has remained unclear. We examined if poor sleep contributes as a causal risk factor to respiratory infections. METHODS: We used data on insomnia, influenza and upper respiratory infections (URIs) from primary care and hospital records in the UK Biobank (N ≈ 231,000) and FinnGen (N ≈ 392,000). We computed logistic regression to assess association between poor sleep and infections, disease free survival hazard ratios, and performed Mendelian randomization analyses to assess causality. FINDINGS: Utilizing 23 years of registry data and follow-up, we discovered that insomnia diagnosis associated with increased risk for infections (FinnGen influenza Cox's proportional hazard (CPH) HR = 4.34 [3.90, 4.83], P = 4.16 × 10-159, UK Biobank influenza CPH HR = 1.54 [1.37, 1.73], P = 2.49 × 10-13). Mendelian randomization indicated that insomnia causally predisposed to influenza (inverse-variance weighted (IVW) OR = 1.65, P = 5.86 × 10-7), URI (IVW OR = 1.94, P = 8.14 × 10-31), COVID-19 infection (IVW OR = 1.08, P = 0.037) and risk of hospitalization from COVID-19 (IVW OR = 1.47, P = 4.96 × 10-5). INTERPRETATION: Our findings indicate that chronic poor sleep is a causal risk factor for contracting respiratory infections, and in addition contributes to the severity of respiratory infections. These findings highlight the role of sleep in maintaining sufficient immune response against pathogens. FUNDING: Instrumentarium Science Foundation, Academy of Finland, Signe and Ane Gyllenberg Foundation, National Institutes of Health.

SIRPα controls CD47-dependent platelet clearance in mice and humans.

Over the last decade, more data has revealed that increased surface expression of the "don't eat me" CD47 protein on cancer cells plays a role in immune evasion and tumor progression, with CD47 blockade emerging as a new therapy in immuno-oncology. CD47 is critical in regulating cell homeostasis and clearance, as binding of CD47 to the inhibitory receptor SIRPα can prevent phagocytosis and macrophage-mediated cell clearance. The purpose of this study was to examine the role of the CD47-SIRPα signal in platelet homeostasis and clearance. Therapeutic reagents targeting the CD47-SIRPα axis are very promising for treatment of hematologic malignancies and solid tumors, but lead to transient anemia or thrombocytopenia in a subset of patients. We found that platelet homeostatic clearance is regulated through the CD47-SIRPα axis and that therapeutic blockade to disrupt this interaction in mice and in humans has a significant impact on platelet levels. Furthermore, we identified genetic variations at the SIRPA locus that impact platelet levels in humans such that higher SIRPA gene expression is associated with higher platelet levels. SIRPA expression at either end of the normal range may affect clinical outcomes of treatment with anti-CD47 therapy.

The public health impact of poor sleep on severe COVID-19, influenza and upper respiratory infections.

Background: Poor sleep is associated with an increased risk of infections and all-cause mortality, and acute sleep loss and disruption have been linked with inflammation and poorer immune control. Previous studies, however, have been unable to evidence causality between the chronic effects of poor sleep and respiratory infection risk. In light of the ongoing COVID-19 pandemic and potential future disease outbreaks, understanding the risk factors for these infections is of great importance. Aim: Our goal was to understand if chronic poor sleep could be identified as a causal risk factor for respiratory infections including influenza, upper respiratory infections and COVID-19. Methods: We used population cohorts from the UK Biobank (N ≈ 231,000) and FinnGen (N ≈ 327,000) with ICD-10 based electronic health records and obtained diagnoses of insomnia, influenza and upper respiratory infections (URIs) from primary care and hospital settings. We computed logistic regression to assess association between poor sleep and infections, disease free survival hazard ratios, and used summary statistics from genome-wide association studies of insomnia, influenza, URI and COVID-19 to perform Mendelian randomization analyses and assess causality. Findings: Utilizing 23 years of registry data and follow-up, we saw that insomnia diagnosis associated with increased risk for infections in FinnGen and in UK Biobank (FinnGen influenza HR = 5.32 [4.09, 6.92], P = 1.02×10-35, UK Biobank influenza HR = 1.54 [1.37, 1.73], P = 2.49×10-13). Mendelian randomization indicated that insomnia causally predisposed to influenza (OR = 1.59, P = 6.23×10-4), upper respiratory infections (OR = 1.71, P = 7.60×10-13), COVID-19 infection (OR = 1.08, P = 0.037) and risk of hospitalization from COVID-19 (OR = 1.47, P = 4.96×10-5). Conclusions: Our findings indicate that chronic poor sleep is a causal risk factor for contracting respiratory infections, and in addition contributes to the severity of respiratory infections. These findings highlight the role of sleep in maintaining sufficient immune response against pathogens as suggested by earlier work. As the current COVID-19 pandemic has increased the number of people suffering from poor sleep, safe interventions such as sleep management and treating individuals with insomnia could be promoted to reduce infections and save lives.

Cellular shape reinforces niche to stem cell signaling in the small intestine.

Niche-derived factors regulate tissue stem cells, but apart from the mechanosensory pathways, the effect of niche geometry is not well understood. We used organoids and bioengineered tissue culture platforms to demonstrate that the conical shape of Lgr5+ small intestinal stem cells (ISCs) facilitate their self-renewal and function. Inhibition of non-muscle myosin II (NM II)-driven apical constriction altered ISC shape and reduced niche curvature and stem cell capacity. Niche curvature is decreased in aged mice, suggesting that suboptimal interactions between old ISCs and their niche develop with age. We show that activation of NM IIC or physical restriction to young topology improves in vitro regeneration by old epithelium. We propose that the increase in lateral surface area of ISCs induced by apical constriction promotes interactions between neighboring cells, and the curved topology of the intestinal niche has evolved to maximize signaling between ISCs and neighboring cells.

From Stem Cells to Populations-Using hiPSC, Next-Generation Sequencing, and GWAS to Explore the Genetic and Molecular Mechanisms of Congenital Heart Defects.

Congenital heart defects (CHD) are developmental malformations affecting the heart and the great vessels. Early heart development requires temporally regulated crosstalk between multiple cell types, signaling pathways, and mechanical forces of early blood flow. While both genetic and environmental factors have been recognized to be involved, identifying causal genes in non-syndromic CHD has been difficult. While variants following Mendelian inheritance have been identified by linkage analysis in a few families with multiple affected members, the inheritance pattern in most familial cases is complex, with reduced penetrance and variable expressivity. Furthermore, most non-syndromic CHD are sporadic. Improved sequencing technologies and large biobank collections have enabled genome-wide association studies (GWAS) in non-syndromic CHD. The ability to generate human to create human induced pluripotent stem cells (hiPSC) and further differentiate them to organotypic cells enables further exploration of genotype-phenotype correlations in patient-derived cells. Here we review how these technologies can be used in unraveling the genetics and molecular mechanisms of heart development.

Mendelian randomization highlights insomnia as a risk factor for pain diagnoses.

STUDY OBJECTIVE: Insomnia has been linked to acute and chronic pain conditions; however, it is unclear whether such relationships are causal. Recently, a large number of genetic variants have been discovered for both insomnia and pain through genome-wide association studies (GWASs) providing a unique opportunity to examine the evidence for causal relationships through the use of the Mendelian randomization paradigm. METHODS: To elucidate the causality between insomnia and pain, we performed bidirectional Mendelian randomization analysis in FinnGen, where clinically diagnosed ICD-10 categories of pain had been evaluated. In addition, we used measures of self-reported insomnia symptoms. We used endpoints for pain in the FinnGen Release 5 (R5) (N = 218,379), and a non-overlapping sample for insomnia (UK Biobank (UKBB) and 23andMe, N = 1,331,010 or UKBB alone N = 453,379). We assessed the robustness of results through conventional Mendelian randomization sensitivity analyses. RESULTS: Genetic liability to insomnia symptoms increased the odds of reporting pain (odds ratio (OR) [95% confidence interval (CI)] = 1.47 [1.38-1.58], p = 4.12 × 10-28). Manifested pain had a small effect on increased risk for insomnia (OR [95% CI] = 1.04 [1.01-1.07], p < 0.05). Results were consistent in sensitivity analyses. CONCLUSIONS: Our findings support a bidirectional causal relationship between insomnia and pain. These data support a further clinical investigation into the utility of insomnia treatment as a strategy for pain management and vice versa.

Sleep apnoea is a risk factor for severe COVID-19.

BACKGROUND: Obstructive sleep apnoea (OSA) is associated with higher body mass index (BMI), diabetes, older age and male gender, which are all risk factors for severe COVID-19.We aimed to study if OSA is an independent risk factor for COVID-19 infection or for severe COVID-19. METHODS: OSA diagnosis and COVID-19 infection were extracted from the hospital discharge, causes of death and infectious diseases registries in individuals who participated in the FinnGen study (n=260 405). Severe COVID-19 was defined as COVID-19 requiring hospitalisation. Multivariate logistic regression model was used to examine association. Comorbidities for either COVID-19 or OSA were selected as covariates. We performed a meta-analysis with previous studies. RESULTS: We identified 445 individuals with COVID-19, and 38 (8.5%) of them with OSA of whom 19 out of 91 (20.9%) were hospitalised. OSA associated with COVID-19 hospitalisation independent from age, sex, BMI and comorbidities (p-unadjusted=5.13×10-5, OR-adjusted=2.93 (95% CI 1.02 to 8.39), p-adjusted=0.045). OSA was not associated with the risk of contracting COVID-19 (p=0.25). A meta-analysis of OSA and severe COVID-19 showed association across 15 835 COVID-19 positive controls, and n=1294 patients with OSA with severe COVID-19 (OR=2.37 (95% 1.14 to 4.95), p=0.021). CONCLUSION: Risk for contracting COVID-19 was the same for patients with OSA and those without OSA. In contrast, among COVID-19 positive patients, OSA was associated with higher risk for hospitalisation. Our findings are in line with earlier works and suggest OSA as an independent risk factor for severe COVID-19.

Comparative genomics highlights the importance of drug efflux transporters during evolution of mycoparasitism in Clonostachys subgenus Bionectria (Fungi, Ascomycota, Hypocreales).

Various strains of the mycoparasitic fungal species Clonostachys rosea are used commercially as biological control agents for the control of fungal plant diseases in agricultural crop production. Further improvements of the use and efficacy of C. rosea in biocontrol require a mechanistic understanding of the factors that determines the outcome of the interaction between C. rosea and plant pathogenic fungi. Here, we determined the genome sequences of 11 Clonostachys strains, representing five species in Clonostachys subgenus Bionectria, and performed a comparative genomic analysis with the aim to identify gene families evolving under selection for gene gains or losses. Several gene families predicted to encode proteins involved in biosynthesis of secondary metabolites, including polyketide synthases, nonribosomal peptide syntethases and cytochrome P450s, evolved under selection for gene gains (p ≤ .05) in the Bionectria subgenus lineage. This was accompanied with gene copy number increases (p ≤ .05) in ATP-binding cassette (ABC) transporters and major facilitator superfamily (MFS) transporters predicted to contribute to drug efflux. Most Clonostachys species were also characterized by high numbers of auxiliary activity (AA) family 9 lytic polysaccharide monooxygenases, AA3 glucose-methanol-choline oxidoreductases and additional carbohydrate-active enzyme gene families with putative activity (or binding) towards xylan and rhamnose/pectin substrates. Particular features of the C. rosea genome included expansions (p ≤ .05) of the ABC-B4 multidrug resistance transporters, the ABC-C5 multidrug resistance-related transporters and the 2.A.1.3 drug:H + antiporter-2 MFS drug resistance transporters. The ABC-G1 pleiotropic drug resistance transporter gene abcG6 in C. rosea was induced (p ≤ .009) by exposure to the antifungal Fusarium mycotoxin zearalenone (1121-fold) and various fungicides. Deletion of abcG6 resulted in mutants with reduced (p < .001) growth rates on media containing the fungicides boscalid, fenhexamid and iprodione. Our results emphasize the role of biosynthesis of, and protection against, secondary metabolites in Clonostachys subgenus Bionectria.

LysM Proteins Regulate Fungal Development and Contribute to Hyphal Protection and Biocontrol Traits in Clonostachys rosea.

Lysin motif (LysM) modules are approximately 50 amino acids long and bind to peptidoglycan, chitin and its derivatives. Certain LysM proteins in plant pathogenic and entomopathogenic fungi are shown to scavenge chitin oligosaccharides and thereby dampen host defense reactions. Other LysM proteins can protect the fungal cell wall against hydrolytic enzymes. In this study, we investigated the biological function of LysM proteins in the mycoparasitic fungus Clonostachys rosea. The C. rosea genome contained three genes coding for LysM-containing proteins and gene expression analysis revealed that lysm1 and lysm2 were induced during mycoparasitic interaction with Fusarium graminearum and during colonization of wheat roots. Lysm1 was suppressed in germinating conidia, while lysm2 was induced during growth in chitin or peptidoglycan-containing medium. Deletion of lysm1 and lysm2 resulted in mutants with increased levels of conidiation and conidial germination, but reduced ability to control plant diseases caused by F. graminearum and Botrytis cinerea. The Δlysm2 strain showed a distinct, accelerated mycelial disintegration phenotype accompanied by reduced biomass production and hyphal protection against hydrolytic enzymes including chitinases, suggesting a role of LYSM2 in hyphal protection against chitinases. The Δlysm2 and Δlysm1Δlysm2 strains displayed reduced ability to colonize wheat roots, while only Δlysm1Δlysm2 failed to suppress expression of the wheat defense response genes PR1 and PR4. Based on our data, we propose a role of LYSM1 as a regulator of fungal development and of LYSM2 in cell wall protection against endogenous hydrolytic enzymes, while both are required to suppress plant defense responses. Our findings expand the understanding of the role of LysM proteins in fungal-fungal interactions and biocontrol.

Natural variation of root lesion nematode antagonism in the biocontrol fungus Clonostachys rosea and identification of biocontrol factors through genome-wide association mapping.

Biological control is a promising approach to reduce plant diseases caused by nematodes to ensure high productivity in agricultural production. Large-scale analyses of genetic variation in fungal species used for biocontrol can generate knowledge regarding interaction mechanisms that can improve efficacy of biocontrol applications. In this study, we performed a genome-wide association study (GWAS) for in vitro antagonism against the root lesion nematode Pratylenchus penetrans in 53 previously genome re-sequenced strains of the biocontrol fungus Clonostachys rosea. Nematode mortality in C. rosea potato dextrose broth (PDB) culture filtrates was highly variable and showed continuous variation (p < .001) between strains, indicating a polygenic inheritance. Twenty-one strains produced culture filtrates with higher (p ≤ .05) nematode mortality compared with the PDB control treatment, while ten strains lowered (p ≤ .05) the mortality. The difference in in vitro antagonism against P. penetrans correlated with antagonism against the soybean cyst nematode Heterodera glycines, indicating lack of host specificity in C. rosea. An empirical Bayesian multiple hypothesis testing approach identified 279 single nucleotide polymorphism markers significantly (local false sign rate < 10-10) associated with the trait. Genes present in the genomic regions associated with nematicidal activity included several membrane transporters, a chitinase and genes encoding proteins predicted to biosynthesize secondary metabolites. Gene deletion strains of the predicted nonribosomal peptide synthetase genes nps4 and nps5 were generated and showed increased (p ≤ .001) fungal growth and conidiation rates compared to the wild type. Deletion strains also exhibited reduced (p < .001) nematicidal activity and reduced (p ≤ .05) biocontrol efficacy against nematode root disease and against fusarium foot rot on wheat. In summary, we show that the GWAS approach can be used to identify biocontrol factors in C. rosea, specifically the putative nonribosomal peptide synthetases NPS4 and NPS5.

Extraction of Microbial and Host DNA, RNA, and Proteins from Oak Bark Tissue.

The application of high-throughput nucleic acid and protein sequencing technologies is transforming our understanding of plant microbiomes and their interactions with their hosts in health and disease. However, progress in studying host-microbiome interactions in above-ground compartments of the tree (the phyllosphere) has been hampered due to high concentrations of phenolic compounds, lignin, and other compounds in tree bark that severely limit the success of DNA, RNA, and protein extraction. Here we present modified sample-preparation and kit-based protocols for the extraction of host and microbiome DNA and RNA from oak (Quercus robus and Quercus petraea) bark tissue for subsequent high-throughput sequencing. In addition, reducing the quantity of bark tissue used for an established protein extraction protocol yielded high quality protein for parallel analysis of the oak-microbiota metaproteome. These procedures demonstrate the successful extraction of nucleic acids and proteins from oak tissue using as little as 50 mg of sample input, producing sufficient quantities for nucleic acid sequencing and protein mass spectrometry of tree stem tissues and their associated microbiota.

Integrated multi-omic analysis of host-microbiota interactions in acute oak decline.

BACKGROUND: Britain's native oak species are currently under threat from acute oak decline (AOD), a decline-disease where stem bleeds overlying necrotic lesions in the inner bark and larval galleries of the bark-boring beetle, Agrilus biguttatus, represent the primary symptoms. It is known that complex interactions between the plant host and its microbiome, i.e. the holobiont, significantly influence the health status of the plant. In AOD, necrotic lesions are caused by a microbiome shift to a pathobiome consisting predominantly of Brenneria goodwinii, Gibbsiella quercinecans, Rahnella victoriana and potentially other bacteria. However, the specific mechanistic processes of the microbiota causing tissue necrosis, and the host response, have not been established and represent a barrier to understanding and managing this decline. RESULTS: We profiled the metagenome, metatranscriptome and metaproteome of inner bark tissue from AOD symptomatic and non-symptomatic trees to characterise microbiota-host interactions. Active bacterial virulence factors such as plant cell wall-degrading enzymes, reactive oxygen species defence and flagella in AOD lesions, along with host defence responses including reactive oxygen species, cell wall modification and defence regulators were identified. B. goodwinii dominated the lesion microbiome, with significant expression of virulence factors such as the phytopathogen effector avrE. A smaller proportion of microbiome activity was attributed to G. quercinecans and R. victoriana. In addition, we describe for the first time the potential role of two previously uncharacterised Gram-positive bacteria predicted from metagenomic binning and identified as active in the AOD lesion metatranscriptome and metaproteome, implicating them in lesion formation. CONCLUSIONS: This multi-omic study provides novel functional insights into microbiota-host interactions in AOD, a complex arboreal decline disease where polymicrobial-host interactions result in lesion formation on tree stems. We present the first descriptions of holobiont function in oak health and disease, specifically, the relative lesion activity of B. goodwinii, G. quercinecans, Rahnella victoriana and other bacteria. Thus, the research presented here provides evidence of some of the mechanisms used by members of the lesion microbiome and a template for future multi-omic research into holobiont characterisation, plant polymicrobial diseases and pathogen defence in trees.