GWAS Analysis Research Articles

Identification of significant SNPs and candidate loci for blast disease resistance via GWAS and population structure analysis in ARC panel of Oryza sativa

Identification of significant SNPs and candidate loci for blast disease resistance via GWAS and population structure analysis in ARC panel of Oryza sativa

Natural variation in the chickpea metabolome under drought stress.

Chickpea is the world's fourth largest grown legume crop, which significantly contributes to food security by providing calories and dietary protein globally. However, the increased frequency of drought stress has significantly reduced chickpea production in recent years. Here, we have performed a field experiment with 36 diverse chickpea genotypes to evaluate grain yield, photosynthetic activities and molecular traits related to drought stress. For metabolomics analysis, leaf tissue was collected at three time points representing different pod-filling stages. We identified L-threonic acid, fructose and sugar alcohols involved in chickpea adaptive drought response within the mid-pod-filling stage. A stress susceptibility index for each genotype was calculated to identify tolerance capacity under drought, distributing the 36 genotypes into four categories from best to worst performance. To understand how biochemical mechanisms control different traits for genetic improvement, we performed a differential Jacobian analysis, which unveiled the interplay between various metabolic pathways across three time points, including higher flux towards inositol interconversions, glycolysis for high-performing genotypes, fumarate to malate conversion, and carbon and nitrogen metabolism perturbations. Metabolic GWAS (mGWAS) analysis uncovered gene candidates involved in glycolysis and MEP pathway corroborating with the differential biochemical Jacobian results. Accordingly, this proposed data analysis strategy bridges the gap from pure statistical association to causal biochemical relations by exploiting natural variation. Our study offers new perspectives on the genetic and metabolic understanding of drought tolerance-associated diversity in the chickpea metabolome and led to the identification of metabolic control points that can be also tested in other legume crops.

Genetic architecture of meat traits in Large White sows

Currently, genome-wide association analysis is a modern and reliable method for analyzing genomic information about animals, as well as determining the “genotype — phenotype” relationship. This study aims to use the GWAS method to identify significant SNPs located within or linked to genes for meat traits in Large White pigs — backfat thickness over the 6–7th and 10–12th vertebrae, and loin muscle depth. The conducted GWA analysis revealed 60 genes, of which 17 are associated with biological functionality, annotated using the DAVID program. Three genes were found to have codification in the Pig QTL database. The genes were divided into 10 groups based on gene ontology (GO). Of all the genes, the AUTS2 gene, located on chromosome 3 and predicting the number of corpora lutea in sows, is of greatest interest. The results of this scientific work will contribute to the development of a genetic evaluation system and improvement of meat qualities in pigs.

Genome-Wide Association Studies (GWAS) and Transcriptome Analysis Reveal Male Heterogametic Sex-Determining Regions and Candidate Genes in Northern Snakeheads (Channa argus).

The Northern snakehead (Channa argus) is a significant economic aquaculture species in China. Exhibiting sexual dimorphism in the growth rate between females and males, mono-sex breeding holds substantial value for aquaculture. This study employed GWAS and transcriptome analysis were applied to identify sex determination genomic regions and develop sex-specific markers. A total of 270 single-nucleotide polymorphisms (SNPs) and 31 insertion-deletions (InDels) were identified as being sexually dimorphic through GWAS and fixation index (Fst) scanning. Based on GWAS results, two sex-specific InDel markers were developed, effectively distinguishing genetic sex for XX females, XY males, and YY super-males via (polymerase chain reaction) PCR amplification. A major genomic segment of approximately 115 kb on chromosome 3 (Chr 03) was identified as the sex-determination region. A comparative transcriptome analysis of gonads for three sexes identified 158 overlapping differentially expressed genes (DEGs). Additionally, three sex-related candidate genes were identified near the sex determination region, including id2, sox11, and rnf144a. Further studies are required to elucidate the functions of these genes. Overall, two sex-specific InDel markers support a male heterogametic XX/XY sex-determination system in Northern snakeheads and three candidate genes offer new insights into sex determination and the evolution of sex chromosomes in teleost fish.

The levels of amino acid metabolites in serum induce the pathogenesis of atopic dermatitis by mediating the inflammatory protein S100A12

Atopic dermatitis (AD) is a chronic inflammatory skin disease affecting tens of millions of people globally. The causal relationship between metabolites and AD pathology has not yet been formally indicated, and the mediating mechanism by which metabolites affect AD has not yet been explored. This study aimed to determine the genetic relationship between metabolites and AD and to determine the pathways through which amino acid metabolites affect AD. Meta-analysis integrates the results of multiple GWAS analyses using METAL software. Using bidirectional two-sample Mendelian randomization (MR), we analyzed the causal relationships between metabolites and AD. The principal MR test of causal effects was conducted using inverse-variance weighted regression, and we used reverse MR analysis to exclude reverse causality. We also performed the MR-PRESSO test to detect and correct for possible pleiotropic effects, and used the Cochran Q test to assess heterogeneity. Two-step MR was utilized to analyze the mediating factors between amino acid metabolites and the onset of AD. The correlation between mediating factors (inflammatory protein S100A12) and immune cell infiltration was analyzed using the edgeR and GSVA software packages. Using single-cell sequencing data from skin tissues of patients with AD, we studied the regulatory role of the S100A12 gene in immune cells. Multiple drug databases and macromolecular docking were used to search for S100A12-targeting drugs. Bidirectional two-sample MR analyses indicated that twenty-two metabolites and one inflammatory protein (S100A12) were significantly associated with AD pathogenesis. S100A12 is a mediator of amino acid metabolites (N6-methyllysine; N2-acetyl,N6,N6-dimethyllysine and N6,N6-dimethyllysine) that are genetically associated with AD. S100A12 was positively correlated with the infiltration of multiple immune cell types in lesional AD skin. The amino acid metabolites N6-methyllysine; N2-acetyl,N6,N6-dimethyllysine and N6,N6-dimethyllysine influence AD pathogenesis by mediating S100A12 expression.

Interactions of the CSF3R polymorphism and early stress on future orientation: evidence for the differential model of stress-related growth.

This study aims to explore the concept of future orientation, which encompasses individuals' thoughts about the future, goal-setting, planning, response to challenges and behavioural adjustments in evolving situations. Often viewed as a psychological resource, future orientation is believed to be developed from psychological resilience. The study investigates the curvilinear relationship between childhood maltreatment and future orientation while examining the moderating effects of genotype. A total of 14,675 Chinese adults self-reported their experiences of childhood maltreatment and their future orientation. The influence of genetic polymorphism was evaluated through genome-wide interaction studies (GWIS; genome-wide association study [GWAS] using gene × environment interaction) and a candidate genes approach. Both GWAS and candidate genes analyses consistently indicated that rs4498771 and its linked single-nucleotide polymorphisms, located in the intergenic area surrounding CSF3R, significantly interacted with early trauma to influence future orientation. Nonlinear regression analyses identified a quadratic or cubic association between future orientation and childhood maltreatment across some genotypes. Specifically, as levels of childhood maltreatment increased, future orientation declined for all genotypes. However, upon reaching a certain threshold, future orientation exhibited a rebound in individuals with specific genotypes. The findings suggest that individuals with certain genotypes exhibit greater resilience to childhood maltreatment. Based on these results, we propose a new threshold model of stress-related growth.

A-167 Hematologic setpoints are a deep physiologic phenotype which enable personalized blood count reference intervals

Abstract Background The complete blood count (CBC) is an important screening tool for healthy adults, commonly ordered at periodic physical exams. However, test results are typically interpreted using one-size-fits-all reference intervals (RIs), which don’t account for the low index-of-individuality of most CBC markers. This undermines the goal of precision medicine - to tailor patient care based on individual characteristics. Methods Here we utilize a database of >50,000 healthy outpatients with yearly CBCs, to investigate the value of patient-specific reference intervals (PRIs) for 9 CBC markers. We use Gaussian mixture modelling to calculate each patient’s setpoints, and to construct associated PRIs. Mechanistic drivers of setpoints were investigated through GWAS and via analysis of excess clinical specimens. Clinical outcomes were analyzed by comparing sensitivity of RIs and PRIs for detection of future cytopenia/cythemia and mortality. Results We found healthy adults had patient-specific setpoints that were stable over decades, and distinctly identifiable from 98% of other patients. Setpoints reflected a deep physiologic phenotype, enabling improved detection of genetic determinants of hematologic regulation (Fig1A), including discovery of multiple novel loci via GWAS analysis. Deviation from PRIs was associated with greater risk of cytopenia/cythemia and mortality (Fig1B) than deviation from the RI for all setpoints except MCHC. Setpoints were significantly correlated with long-term all-cause mortality, even when limited to those within the RI (Fig1C). Conclusions We show that hematologic setpoints reflect a deep, patient-specific phenotype. This phenotype can enhance mechanistic studies, and allow for increased diagnostic sensitivity, through calculation of a personalized reference interval. Setpoints are sufficiently stable and patient-specific to help realize the promise of precision medicine for healthy adults.Fig1A.Significance of GWAS hits using setpoints versus single CBCs. Fig1B. 10yr mortality hazard ratio for tests outside the RI and PRI. Fig1C. Association of WBC setpoint with mortality (limited to WBC setpoints inside the RI).

Integrated GWAS and transcriptome analysis reveals key genes associated with muscle fibre and fat traits in Gushi chicken

ABSTRACT 1. In the following experiment meat quality traits of a Gushi-Anka F2 resource population were measured, and their heritability estimated. Intramuscular fat (IMF) had medium heritability (0.35) but leg muscle fibre density (LMD), leg muscle fibre diameter (LMF), breast muscle fibre density (BMD), fresh fat content (FFA), and absolute dry fat content (AFC) had low heritability (0–0.2). The IMF presented the most important genetic additive effect among the poultry meat quality-related traits studied. 2. The phenotypic data of meat quality traits in the Gushi-Anka F2 resource population were combined with genotyping by sequencing (GBS) data to obtain genotype data. Six meat quality traits in 734 birds were analysed by GWAS. Based on these variants, 83 significant (–log10(p) > 4.42) single nucleotide polymorphisms and four quantitative trait loci (QTL) regions corresponding to 175 genes were identified. Further linkage disequilibrium (LD) analysis was conducted on chromosome 13 (Chr13) and chromosome 27 (Chr27) QTL regions. 3. Based on the transcriptome data and GWAS results, 12 shared genes - ITGB3, DNAJC27, ETV4, C7orf50, FKBP1B, G3BP1, IGF2BP1, KCNH6, LOC416263, SCARA5, SMIM5 and TBL1XR1 were identified as candidate genes influencing muscle fibre and fat traits.

Genetics of suicide ideation. A role for inflammation and neuroplasticity?

Suicide is a leading cause of death worldwide. Suicide ideation (SI) is a known risk factor for suicide behaviour (SB). The current psychobiology and genetic predisposition to SI and SB are poorly defined. Despite convincing relevance of a genetic background for SI, there is no current implementable knowledge about the genetic makeup that identifies subjects at risk for it. One of the possible reasons for the absence of a clear-cut evidence is the polygenetic nature of SI along with the very large sample sizes that are needed to observe significant genetic association result. The CATIE sample was instrumental to the analysis. SI was retrieved as measured by the Calgary test. Clinical possible covariates were identified by a nested regression model. A principal component analysis helped in defining the possible genetic stratification factors. A GWAS analysis, polygenic risk score associated with a random forest analysis and a molecular pathway analysis were undertaken to identify the genetic contribution to SI. As a result, 741 Schizophrenic individuals from the CATIE were available for the genetic analysis, including 166,325 SNPs after quality control and pruning. No GWAS significant result was found. The random forest analysis conducted by combining the polygenic risk score and several clinical variables resulted in a possibly overfitting model (OOB error rate < 1%). The molecular pathway analysis revealed several molecular pathways possibly involved in SI, of which those involved in microglia functioning were of particular interest. A medium-small sample of SKZ individuals was analyzed to shed a light on the genetic of SI. As an expected result from the underpowered sample, no GWAS positive result was retrieved, but the molecular pathway analysis indicated a possible role of microglia and neurodevelopment in SI.

Identification of Candidate Genes for Cold Tolerance at Seedling Stage by GWAS in Rice (Oryza sativa L.).

Due to global climate change, cold temperatures have significantly impacted rice production, resulting in reduced yield and quality. In this study, we investigated two traits related to the cold tolerance (CT) of 1992 diverse rice accessions at the seedling stage. Geng accessions exhibited higher levels of CT compared to xian accessions, with the GJ-tmp subgroup displaying the strongest CT. However, extreme CT accessions were also identified within the xian subspecies. Through GWAS analysis based on the survival rate (SR) and leaf score of cold tolerance (SCT), a total of 29 QTLs associated with CT at the seedling stage were identified, among which four QTLs (qSR3.1a, qSR4.1a, qSR11.1x, and qSR12.1a) were found to be important. Furthermore, five candidate genes (LOC_Os03g44760, LOC_Os04g06900, LOC_Os04g07260, LOC_Os11g40610, and LOC_Os12g10710) along with their favorable haplotypes were identified through gene function annotation and haplotype analysis. Pyramiding multiple favorable haplotypes resulted in a significant improvement in CT performance. Subsequently, three selected accessions (CX534, B236, and IRIS_313-8565), carrying different superior alleles for CT, were selected and recommended for molecular breeding for CT using marker-assisted selection (MAS). The findings from this study provide valuable resources for enhancing rice's ability for CT while laying a foundation for the future cloning of novel genes involved in conferring CT.

Abstract A007: Beyond continental ancestry categories: trans-population genomic similarity continuum reveals overlooked somatic-germline interactions across cancer types

Abstract A growing body of research highlighted a strong presence of continuums that describe human genomic diversity. Historical scholarship has also extensively investigated the use of discrete classification systems in human genomics further questioning their validity. To investigate the limiting effect of pre-defined stratification, we performed a wide comparative driver GWAS analysis at the trans-population (i.e. without stratification) vs continent-based population levels across 16 cancer types. Results demonstrate that in the MSK-IMPACT cohort (n = 56,826 samples), the analysis of the trans-population genomic similarity continuum captures significant somatic-germline interactions (p &amp;lt; 5x10-8) in cancer type-specific and pancancer cohorts, with a set of interactions overlooked by the discrete continent-based paradigm. We demonstrate that overlooked germline SNPs have similar allele frequencies and directional relationships with somatic alterations across pre-defined discrete categories, and that earlier stratification prevents their detection due to lowered statistical power. PheWAS analysis of detected germline variants revealed associations with cancer and non-cancer phenotypes reported in previous studies as well as novel hits. To expand the scope of the trans-population lens, we carried out a comparative analysis on 17 whole-exome TCGA cohorts to identify potential cancer drivers based on mutational frequencies. Results similarly revealed known and potential cancer drivers at the trans-population level that would be overlooked when pre-defined continent-based stratification is applied. The lens we adopt further highlights the necessity for new quantitative perspectives to comprehensively detect genomic patterns and to leverage entire genomic datasets without pre-defined stratification schemes. Citation Format: Hussein Mohsen, Tomin Perea-Chamblee, Jian Carrot-Zhang. Beyond continental ancestry categories: trans-population genomic similarity continuum reveals overlooked somatic-germline interactions across cancer types [abstract]. In: Proceedings of the 17th AACR Conference on the Science of Cancer Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; 2024 Sep 21-24; Los Angeles, CA. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2024;33(9 Suppl):Abstract nr A007.

Genome-wide association study of lignin trait in elite spring wheat against spot blotch disease

Spot blotch disease in wheat is one of major yield-limiting factor in warm humid conditions of South Asia. Developing the disease-resistant variety through molecular methods is an economically and environmentally eco-friendly approach. Diseased wheat leaves associated with lignin content increased by two or three times at 72 h compared with control. The genomic region associated with lignin content was investigated in a set of 289 diverse spring wheat genotypes. The genome - wide association studies panel employed a 90K Illumina SNP array and phenotyped in four environments. The GWAS analysis showed a total of 86 marker-trait associations on 1A, 1B, 2A, 2B, 2D, 3A, 3D, 4A, 5A, 5B, 6A, 6D, 7A, 7B, 7D were linked lignin trait. Insilco analysis revealed that significant SNPs were located on import putative candidate genes thirteen in the Transmembrane domain of the protein-kinase family and were involved in the host-pathogen interaction. The identified significant MTAs for further validation and useful for the wheat breeding programs to develop the spot blotch disease resistance.

Elucidating Genetic Mechanisms of Summer Stress Tolerance in Chinese Cabbage through GWAS and Phenotypic Analysis

Elucidating Genetic Mechanisms of Summer Stress Tolerance in Chinese Cabbage through GWAS and Phenotypic Analysis

GWAS Analysis to Identify Candidate Genes Related to Phosphorus Deficiency Tolerance by GWAS in Rice

GWAS Analysis to Identify Candidate Genes Related to Phosphorus Deficiency Tolerance by GWAS in Rice

Genome-wide association study for metabolic syndrome reveals APOA5 single nucleotide polymorphisms with multilayered effects in Koreans

Background and purposeGenome-wide association studies (GWAS) of metabolic syndrome (MetS) have predominantly focused on non-Asian populations, with limited representation from East Asian cohorts. Moreover, previous GWAS analyses have primarily emphasized the significance of top single nucleotide polymorphisms (SNPs), poorly explaining other SNP signals in linkage disequilibrium. This study aimed to reveal the interaction between rs651821 and rs2266788, the principal variants of apolipoprotein A5 (APOA5), within the most significant loci identified through GWAS on MetS.MethodsGWAS on MetS and its components was conducted using the data from the Korean Genome and Epidemiology Study (KoGES) city cohort comprising 58,600 individuals with available biochemical, demographic, lifestyle factors, and the most significant APOA5 locus was analyzed further in depth.ResultsAccording to GWAS of MetS and its diagnostic components, a significant association between the APOA5 SNPs rs651821/rs2266788 and MetS/triglycerides/high-density lipoprotein phenotypes was revealed. However, a conditional analysis employing rs651821 unveiled a reversal in the odds ratio for rs2266788. Therefore, rs651821 and rs2266788 emerged as independent and opposing signals in the extended GWAS analysis, i.e., the multilayered effects. Further gene-environment interaction analyses regarding lifestyle factors such as smoking, alcohol consumption, and physical activity underscored these multilayered effects.ConclusionThis study unveils the intricate interplay between rs651821 and rs2266788 derived from MetS GWAS. Removing the influence of lead SNP reveals an independent protective signal associated with rs2266788, suggesting a multilayered effect between these SNPs. These findings underline the need for novel perspectives in future MetS GWAS.

Developing genomic tools to assist turnip rape [Brassica rapa (L.) subsp.oleifera (DC.) Metzg.] breeding

IntroductionTurnip rape is recognized as an oilseed crop contributing to environmentally sustainable agriculture via integration into crop rotation systems. Despite its various advantages, the crop’s cultivation has declined globally due to a relatively low productivity, giving way to other crops. The use of genomic tools could enhance the breeding process and accelerate genetic gains. Therefore, the present research investigated 170 turnip rape accessions representing its global gene pool to identify SNP markers associated nine phenological and agro-morphological traits and estimate the genomic breeding values (GEBVs) of the germplasm through GWAS and genomic prediction analyses, respectively.MethodsField trials were conducted at two sites in northern and southern Sweden to obtain the phenotypic data while genotyping was conducted via the genotyping-by-sequencing (GBS) method. The traits studied include days to flowering (DTF) and maturity (DTM), plant height (PH), seed yield (YLD), thousand seed weight (TSW), silique length (SL), number of siliques (NS), number of seeds per silique (SS), and pod shattering resistance (PSHR).Results and conclusionAnalysis of variance revealed substantial variation among accessions, with significant genotype-by-environment interaction for most traits. A total of 25, 17, 16, 14, 7, 5, 3, and 3 MTAs were identified for TSW, DTF, PH, PSHR, SL, YLD, SS and DTM, respectively. An 80%–20% training-test set genomic prediction analysis was conducted using the ridge regression – BLUP (RR-BLUP) model. The accuracy of genomic prediction for most traits was high, indicating that these tools may assist turnip rape breeders in accelerating genetic gains. The study highlights the potential of genomic tools to significantly advance breeding programs for turnip rape by identifying pivotal SNP markers and effectively estimating genomic breeding values. Future breeding perspectives should focus on leveraging these genomic insights to enhance agronomic traits and productivity, thereby reinstating turnip rape as a competitive and sustainable crop in Sweden and broader global agriculture.

Multitrait Genetic Analysis Identifies Novel Pleiotropic Loci for Depression and Schizophrenia in East Asians.

While genetic correlations, pleiotropic loci, and shared genetic mechanisms of psychiatric disorders have been extensively studied in European populations, the investigation of these factors in East Asian populations has been relatively limited. To identify novel pleiotropic risk loci for depression and schizophrenia (SCZ) in East Asians. We utilized the most comprehensive dataset available for East Asians and quantified the genetic overlap between depression, SCZ, and their related traits via a multitrait genome-wide association study. Global and local genetic correlations were estimated by LDSC and ρ-HESS. Pleiotropic loci were identified by the multitrait analysis of GWAS (MTAG). Besides the significant correlation between depression and SCZ, our analysis revealed genetic correlations between depression and obesity-related traits, such as weight, BMI, T2D, and HDL. In SCZ, significant correlations were detected with HDL, heart diseases and use of various medications. Conventional meta-analysis of depression and SCZ identified a novel locus at 1q25.2 in East Asians. Further multitrait analysis of depression, SCZ and related traits identified ten novel pleiotropic loci for depression, and four for SCZ. Our findings demonstrate shared genetic underpinnings between depression and SCZ in East Asians, as well as their associated traits, providing novel candidate genes for the identification and prioritization of therapeutic targets specific to this population.

A rice amino acid transporter OsMP1 acts as a quantitative trait locus against blast fungus and leaf-blight bacterium.

Amino acid homeostasis is interconnected with the immune network of plants. During plant-pathogen interaction, amino acid transporters (AATs) have been shown to be involved in plant immune responses. However, the molecular mechanism by which how AATs function in this process remains elusive. In this study, we identify OsMP1 that acts as a quantitative trait locus against blast fungus from a joint analysis of GWAS and QTL mapping in rice. Heterogeneous expression of OsMP1 in yeast supports its function in transporting a wide range of amino acids, including Thr, Ser, Phe, His and Glu. OsMP1 could also mediate 15N-Glu efflux and influx in Xenopus oocyte cells. The expression of OsMP1 is dramatically induced by Magnaporthe oryzae in the resistant landrace Heikezijing, while remaining unresponsive in the susceptible landrace Suyunuo. Overexpressing OsMP1 in Suyunuo enhances disease resistance to blast fungus and leaf-blight bacterium without yield penalty. Furthermore, the overexpression of OsMP1 leads to increased accumulation of Thr, Ser, Phe and His in the leaves. And the heightened levels of these amino acids contribute to reduced disease susceptibility, which is associated with upregulated jasmonic acid pathway. Thus, our results elucidate the pivotal role of OsMP1 in disease resistance and provide a potential target for breeding more resistant rice cultivars without compromising yield.

Mining for QTL controlling maize low-phosphorus response genes combined with deep resequencing of RIL parental genomes and in silico GWAS analysis.

Extensive and comprehensive phenotypic data from a maize RIL population under both low- and normal-Pi treatments were used to conduct QTL mapping. Additionally, we integrated parental resequencing data from the RIL population, GWAS results, and transcriptome data to identify candidate genes associated with low-Pi stress in maize. Phosphorus (Pi) is one of the essential nutrients that greatly affect the maize yield. However, the genes underlying the QTL controlling maize low-Pi response remain largely unknown. In this study, a total of 38 traits at both seedling and maturity stages were evaluated under low- and normal-Pi conditions using a RIL population constructed from X178 (tolerant) and 9782 (sensitive), and most traits varied significantly between low- and normal-Pi treatments. Twenty-nine QTLs specific to low-Pi conditions were identified after excluding those with common intervals under both low- and normal-Pi conditions. Furthermore, 45 additional QTLs were identified based on the index value ((Trait_under_LowPi-Trait_under_NormalPi)/Trait_under_NormalPi) of each trait. These 74 QTLs collectively were classified as Pi-dependent QTLs. Additionally, 39 Pi-dependent QTLs were clustered in nine HotspotQTLs. The Pi-dependent QTL interval contained 19,613 unique genes, 6,999 of which exhibited sequence differences with non-synonymous mutation sites between X178 and 9782. Combined with in silico GWAS results, 277 consistent candidate genes were identified, with 124 genes located within the HotspotQTL intervals. The transcriptome analysis revealed that 21 genes, including the Pi transporter ZmPT7 and the strigolactones pathway-related gene ZmPDR1, exhibited consistent low-Pi stress response patterns across various maize inbred lines or tissues. It is noteworthy that ZmPDR1 in maize roots can be sharply up-regulated by low-Pi stress, suggesting its potential importance as a candidate gene for responding to low-Pi stress through the strigolactones pathway.

Genomic characterization and virulence gene profiling of Erysipelothrix rhusiopathiae isolated from widespread muskox mortalities in the Canadian Arctic Archipelago

BackgroundMuskoxen are important ecosystem components and provide food, economic opportunities, and cultural well-being for Indigenous communities in the Canadian Arctic. Between 2010 and 2021, Erysipelothrix rhusiopathiae was isolated from carcasses of muskoxen, caribou, a seal, and an Arctic fox during multiple large scale mortality events in the Canadian Arctic Archipelago. A single strain (‘Arctic clone’) of E. rhusiopathiae was associated with the mortalities on Banks, Victoria and Prince Patrick Islands, Northwest Territories and Nunavut, Canada (2010–2017). The objectives of this study were to (i) characterize the genomes of E. rhusiopathiae isolates obtained from more recent muskox mortalities in the Canadian Arctic in 2019 and 2021; (ii) identify and compare common virulence traits associated with the core genome and mobile genetic elements (i.e. pathogenicity islands and prophages) among Arctic clone versus other E. rhusiopathiae genomes; and iii) use pan-genome wide association studies (GWAS) to determine unique genetic contents of the Arctic clone that may encode virulence traits and that could be used for diagnostic purposes.ResultsPhylogenetic analyses revealed that the newly sequenced E. rhusiopathiae isolates from Ellesmere Island, Nunavut (2021) also belong to the Arctic clone. Of 17 virulence genes analysed among 28 Arctic clone isolates, four genes – adhesin, rhusiopathiae surface protein-A (rspA), choline binding protein-B (cbpB) and CDP-glycerol glycerophosphotransferase (tagF) – had amino acid sequence variants unique to this clone when compared to 31 other E. rhusiopathiae genomes. These genes encode proteins that facilitate E. rhusiopathiae to attach to the host endothelial cells and form biofilms. GWAS analyses using Scoary found several unique genes to be overrepresented in the Arctic clone.ConclusionsThe Arctic clone of E. rhusiopathiae was associated with multiple muskox mortalities spanning over a decade and multiple Arctic islands with distances over 1000 km, highlighting the extent of its spatiotemporal spread. This clone possesses unique gene content, as well as amino acid variants in multiple virulence genes that are distinct from the other closely related E. rhusiopathiae isolates. This study establishes an essential foundation on which to investigate whether these differences are correlated with the apparent virulence of this specific clone through in vitro and in vivo studies.