Genome-wide Association Study Approach Research Articles

A multi-trait GWAS identifies novel genes influencing albuminuria.

Albuminuria is common and associated with increased risks of end-stage kidney disease and cardiovascular diseases, yet its underlying mechanism remains obscure. Previous genome-wide association studies (GWAS) for albuminuria did not consider gene pleiotropy and primarily focused on European ancestry populations. This study adopted a multi-trait analysis of GWAS (MTAG) approach to jointly analyze two vital kidney traits, estimated glomerular filtration rate (eGFR) and urine albumin-to-creatinine ratio (UACR) to identify and prioritize the genes associated with UACR. Data from the Taiwan Biobank from 2012 to 2023 were analyzed. GWAS of UACR and eGFR were performed separately and the summary statistics from these GWAS were jointly analyzed using MTAG. The polygenic risk scores (PRS) of UACR were constructed for validation. The UACR-associated loci were further fine-mapped and prioritized based on their deleteriousness, eQTL associations, and relatedness to Mendelian kidney diseases. MTAG analysis of the UACR revealed 15 genetic loci, including 12 novel loci. The PRS for UACR was significantly associated with urinary albumin level (P<0.001) and microalbuminuria (P=0.001 ∼ 0.045). A list of priority genes was generated. Twelve genes with high priority included the albumin endocytic receptor gene LRP2 and ciliary genes IFT172. The findings of this multi-trait GWAS suggest that primary cilia play a role in sensing mechanical stimuli, leading to albumin endocytosis. The priority list of genes warrants further translational investigation to reduce albuminuria.

Genome sequence analyses identify novel risk loci for multiple system atrophy.

Multiple system atrophy (MSA) is an adult-onset, sporadic synucleinopathy characterized by parkinsonism, cerebellar ataxia, and dysautonomia. The genetic architecture of MSA is poorly understood, and treatments are limited to supportive measures. Here, we performed a comprehensive analysis of whole genome sequence data from 888 European-ancestry MSA cases and 7,128 controls to systematically investigate the genetic underpinnings of this understudied neurodegenerative disease. We identified four significantly associated risk loci using a genome-wide association study approach. Transcriptome-wide association analyses prioritized USP38-DT, KCTD7, and lnc-KCTD7-2 as novel susceptibility genes for MSA within these loci, and single-nucleus RNA sequence analysis found that the associated variants acted as cis-expression quantitative trait loci for multiple genes across neuronal and glial cell types. In conclusion, this study highlights the role of genetic determinants in the pathogenesis of MSA, and the publicly available data from this study represent a valuable resource for investigating synucleinopathies.

Inter-individual variation in ovarian reserve after gonadotoxic treatment in female childhood cancer survivors – a genome-wide association study: results from PanCareLIFE

ObjectiveWe aimed to discover new variants associated with low ovarian reserve after gonadotoxic treatment among adult female childhood cancer survivors using a genome-wide association study approach. DesignGenome-wide association study. SubjectsA discovery cohort of adult female childhood cancer survivors, from the pan‐European PanCareLIFE cohort (n=743; median age: 25.8 years), excluding those who received bilateral ovarian irradiation, bilateral oophorectomy, central nerve system or total body irradiation, or stem cell transplantation. Replication was attempted in the USA‐based St. Jude Lifetime Cohort (n=391; median age: 31.3 years). ExposureFemale childhood cancer survivors are at risk of therapy‐related gonadal impairment. Alkylating agents are well‐established risk factors, and the inter‐individual variability in gonadotoxicity may be explained by genetic polymorphisms. Data were collected in real-life conditions and cyclophosphamide equivalent dose was used to quantify alkylation agent exposure. InterventionNo intervention was performed. Main Outcome MeasureAnti‐Müllerian hormone (AMH) levels served as a proxy for ovarian function and findings were combined in a meta‐analysis. ResultsThree genome-wide significant (<5.0x10-8) and 16 genome-wide suggestive (<5.0x10-6) loci were associated with log-transformed AMH levels, adjusted for cyclophosphamide equivalent dose of alkylating agents, age at diagnosis, and age at study in the PanCareLIFE cohort. Based on effect allele frequency (EAF) (>0.01 if not genome-wide significant), p-value (<5.0×10−6), and biological relevance, 15 SNPs were selected for replication. None of the SNPs were statistically significantly associated with AMH levels. A meta-analysis indicated that rs78861946 was associated at borderline genome-wide statistical significance (Reference/effect allele: C/T; EAF: 0.04, Beta (SE): −0.484 (0.091), p‐value= 9.39×10−8). ConclusionThis study found no genetic variants associated with a lower ovarian reserve after gonadotoxic treatment, as the findings of this GWAS were not statistically significant replicated in the replication cohort. Suggestive evidence for potential importance of one variant is briefly discussed, but the lack of statistical significance calls for larger cohort sizes. As the population of childhood cancer survivors is increasing, large-scale and systematic research is needed to identify genetic variants that could aid predictive risk models of gonadotoxicity and as well as fertility preservation options for childhood cancer survivors.

Genome-wide association study and haplotype analysis reveal novel candidate genes for resistance to powdery mildew in soybean.

Powdery mildew disease (PMD) is caused by the obligate biotrophic fungus Microsphaera diffusa Cooke & Peck (M. diffusa) and results in significant yield losses in soybean (Glycine max (L.) Merr.) crops. By identifying disease-resistant genes and breeding soybean accessions with enhanced resistance, we can effectively mitigate the detrimental impact of PMD on soybeans. We analyzed PMD resistance in a diversity panel of 315 soybean accessions in two locations over 3 years, and candidate genes associated with PMD resistance were identified through genome-wide association studies (GWAS), haplotype analysis, qRT-PCR, and EMS mutant analysis. Based on the GWAS approach, we identified a region on chromosome 16 (Chr16) in which 21 genes form a gene cluster that is highly correlated with PMD resistance. In order to validate and refine these findings, we conducted haplotype analysis of 21 candidate genes and indicated there are single nucleotide polymorphisms (SNPs) and insertion-deletions (InDels) variations of Glyma.16G214000, Glyma.16G214200, Glyma.16G215100 and Glyma.16G215300 within the coding and promoter regions that exhibit a strong association with resistance against PMD. Subsequent structural analysis of candidate genes within this cluster revealed that in 315 accessions, the majority of accessions exhibited resistance to PMD when Glyma.16G214300, Glyma.16G214800 and Glyma.16G215000 were complete; however, they demonstrated susceptibility to PMD when these genes were incomplete. Quantitative real-time PCR assays (qRT-PCR) of possible candidate genes showed that 14 candidate genes (Glyma.16G213700, Glyma.16G213800, Glyma.16G213900, Glyma.16G214000, Glyma.16G214200, Glyma.16G214300, Glyma.16G214500, Glyma.16G214585, Glyma.16G214669, Glyma.16G214700, Glyma.16G214800, Glyma.16G215000, Glyma.16G215100 and Glyma.16G215300) were involved in PMD resistance. Finally, we evaluated the PMD resistance of mutant lines from the Williams 82 EMS mutations library, which revealed that mutants of Glyma.16G214000, Glyma.16G214200, Glyma.16G214300, Glyma.16G214800, Glyma.16G215000, Glyma.16G215100 and Glyma.16G215300, exhibited sensitivity to PMD. Combined with the analysis results of GWAS, haplotypes, qRT-PCR and mutants, the genes Glyma.16G214000, Glyma.16G214200, Glyma.16G214300, Glyma.16G214800, Glyma.16G215000, Glyma.16G215100 and Glyma.16G215300 were identified as highly correlated with PMD resistance. The candidate genes identified above are all NLR family genes, and these discoveries deepen our understanding of the molecular basis of PMD resistance in soybeans and will be useful for guiding breeding strategies.

Genome-Wide Genetic Analysis of Dropout in a Controlled Exercise Intervention in Sedentary Adults With Overweight or Obesity and Cardiometabolic Disease.

Despite the benefits of exercise, many individuals are unable or unwilling to adopt an exercise intervention. The purpose of this analysis was to identify putative genetic variants associated with dropout from exercise training interventions among individuals in the STRRIDE trials. We used a genome-wide association study approach to identify genetic variants in 603 participants initiating a supervised exercise intervention. Exercise intervention dropout occurred when a subject withdrew from further participation in the study or was otherwise lost to follow-up. Exercise intervention dropout was associated with a cluster of single-nucleotide polymorphisms with the top candidate being rs722069 (T/C, risk allele = C) (unadjusted p = 2.2 × 10-7, odds ratio = 2.23) contained within a linkage disequilibrium block on chromosome 16. In Genotype-Tissue Expression, rs722069 is an expression quantitative trait locus of the EARS2, COG7, and DCTN5 genes in skeletal muscle tissue. In subsets of the STRRIDE genetic cohort with available muscle gene expression (n = 37) and metabolic data (n = 82), at baseline the C allele was associated with lesser muscle expression of EARS2 (p < .002) and COG7 (p = .074) as well as lesser muscle concentrations of C2- and C3-acylcarnitines (p = .026). Our observations imply that exercise intervention dropout is genetically moderated through alterations in gene expression and metabolic pathways in skeletal muscle. Individual genetic traits may allow the development of a biomarker-based approach for identifying individuals who may benefit from more intensive counseling and other interventions to optimize exercise intervention adoption. STRRIDE I = NCT00200993; STRRIDE AT/RT = NCT00275145; STRRIDE-PD = NCT00962962.

Intergenomic signatures of coevolution between Tasmanian devils and an infectious cancer

Coevolution is common and frequently governs host-pathogen interaction outcomes. Phenotypes underlying these interactions often manifest as the combined products of the genomes of interacting species, yet traditional quantitative trait mapping approaches ignore these intergenomic interactions. Devil facial tumor disease (DFTD), an infectious cancer afflicting Tasmanian devils (Sarcophilus harrisii), has decimated devil populations due to universal host susceptibility and a fatality rate approaching 100%. Here, we used a recently developed joint genome-wide association study (i.e., co-GWAS) approach, 15 y of mark-recapture data, and 960 genomes to identify intergenomic signatures of coevolution between devils and DFTD. Using a traditional GWA approach, we found that both devil and DFTD genomes explained a substantial proportion of variance in how quickly susceptible devils became infected, although genomic architectures differed across devils and DFTD; the devil genome had fewer loci of large effect whereas the DFTD genome had a more polygenic architecture. Using a co-GWA approach, devil-DFTD intergenomic interactions explained ~3× more variation in how quickly susceptible devils became infected than either genome alone, and the top genotype-by-genotype interactions were significantly enriched for cancer genes and signatures of selection. A devil regulatory mutation was associated with differential expression of a candidate cancer gene and showed putative allele matching effects with two DFTD coding sequence variants. Our results highlight the need to account for intergenomic interactions when investigating host-pathogen (co)evolution and emphasize the importance of such interactions when considering devil management strategies.

Bone health index in the assessment of bone health: The Generation R Study

Bone Health Index (BHI) has been proposed as a useful instrument for assessing bone health in children. However, its relationship with fracture risk remains unknown. We aimed to investigate whether BHI is associated with bone mineral density (BMD) and prevalent fracture odds in children from the Generation R Study. We also implemented genome-wide association study (GWAS) and polygenic score (PGS) approaches to improve our understanding of BHI and its potential. In total, 4150 children (49.4 % boys; aged 9.8 years) with genotyped data and bone assessments were included in this study. BMD was measured across the total body (less head following ISCD guidelines) using a GE-Lunar iDXA densitometer; and BHI was determined from the hand DXA scans using BoneXpert®. Fractures were self-reported collected with home questionnaires. The association of BHI with BMD and fractures was evaluated using linear models corrected for age, sex, ethnicity, height, and weight. We observed a positive correlation between BHI and BMD (ρ = 0.32, p-value<0.0001). Further, every SD decrease in BHI was associated with an 11 % increased risk of prevalent fractures (OR:1.11, 95 % CI 1.00–1.24, p-value = 0.05). Our BHI GWAS identified variants (lead SNP rs1404264-A, p-value = 2.61 × 10−14) mapping to the ING3/CPED1/WNT16 locus. Children in the extreme tails of the BMD PGS presented a difference in BHI values of −0.10 standard deviations (95% CI -0.14 to −0.07; p-value<0.0001). On top of the demonstrated epidemiological association of BHI with both BMD and fracture risk, our results reveal a partially shared biological background between BHI and BMD. These findings highlight the potential value of using BHI to screen children at risk of fracture.

Association of FGF4L1 Retrogene Insertion with Prolapsed Gland of the Nictitans (Cherry Eye) in Dogs.

Cherry eye is the common name for prolapse of the nictitans gland, a tear-producing gland situated under the third eyelid of dogs. Cherry eye is characterized by a red fleshy protuberance in the corner of the eye, resembling a cherry. This protrusion is a displacement of the normal gland of the third eyelid, thought to be caused by a defect in the connective tissue that secures the gland in place. Options for treatment may include anti-inflammatory medications in mild cases, but surgical replacement of the gland is usually indicated. Cherry eye is most often seen in dogs under the age of two years, with certain breeds having a higher incidence, suggesting a potential genetic association. Integration of panel genetic testing into routine clinical practice allows for the generation of large numbers of genotyped individuals paired with clinical records and enables the investigation of common disorders using a genome-wide association study (GWAS) approach at scale. In this investigation, several thousand cases and controls for cherry eye in both purebred dogs and mixed breeds are used for a large-scale GWAS, revealing a single peak of genome-wide significance on canine chromosome 18, directly at the location of the previously identified FGF4 insertion known to cause chondrodysplasia in several breeds.

Glutaredoxin regulation of primary root growth is associated with early drought stress tolerance in pearl millet

Seedling root traits impact plant establishment under challenging environments. Pearl millet is one of the most heat and drought tolerant cereal crops that provides a vital food source across the sub-Saharan Sahel region. Pearl millet’s early root system features a single fast-growing primary root which we hypothesize is an adaptation to the Sahelian climate. Using crop modeling, we demonstrate that early drought stress is an important constraint in agrosystems in the Sahel where pearl millet was domesticated. Furthermore, we show that increased pearl millet primary root growth is correlated with increased early water stress tolerance in field conditions. Genetics including genome-wide association study and quantitative trait loci (QTL) approaches identify genomic regions controlling this key root trait. Combining gene expression data, re-sequencing and re-annotation of one of these genomic regions identified a glutaredoxin-encoding gene PgGRXC9 as the candidate stress resilience root growth regulator. Functional characterization of its closest Arabidopsis homolog AtROXY19 revealed a novel role for this glutaredoxin (GRX) gene clade in regulating cell elongation. In summary, our study suggests a conserved function for GRX genes in conferring root cell elongation and enhancing resilience of pearl millet to its Sahelian environment.

Identification of a Branch Number Locus in Soybean Using BSA-Seq and GWAS Approaches.

The determination of the soybean branch number plays a pivotal role in plant morphogenesis and yield components. This polygenic trait is subject to environmental influences, and despite its significance, the genetic mechanisms governing the soybean branching number remain incompletely understood. To unravel these mechanisms, we conducted a comprehensive investigation employing a genome-wide association study (GWAS) and bulked sample analysis (BSA). The GWAS revealed 18 SNPs associated with the soybean branch number, among which qGBN3 on chromosome 2 emerged as a consistently detected locus across two years, utilizing different models. In parallel, a BSA was executed using an F2 population derived from contrasting cultivars, Wandou35 (low branching number) and Ruidou1 (high branching number). The BSA results pinpointed a significant quantitative trait locus (QTL), designated as qBBN1, located on chromosome 2 by four distinct methods. Importantly, both the GWAS and BSA methods concurred in co-locating qGBN3 and qBBN1. In the co-located region, 15 candidate genes were identified. Through gene annotation and RT-qPCR analysis, we predicted that Glyma.02G125200 and Glyma.02G125600 are candidate genes regulating the soybean branch number. These findings significantly enhance our comprehension of the genetic intricacies regulating the branch number in soybeans, offering promising candidate genes and materials for subsequent investigations aimed at augmenting the soybean yield. This research represents a crucial step toward unlocking the full potential of soybean cultivation through targeted genetic interventions.

Performing highly parallelized and reproducible GWAS analysis on biobank-scale data.

Genome-wide association studies (GWAS) are transforming genetic research and enable the detection of novel genotype-phenotype relationships. In the last two decades, over 60 000 genetic associations across thousands of traits have been discovered using a GWAS approach. Due to increasing sample sizes, researchers are increasingly faced with computational challenges. A reproducible, modular and extensible pipeline with a focus on parallelization is essential to simplify data analysis and to allow researchers to devote their time to other essential tasks. Here we present nf-gwas, a Nextflow pipeline to run biobank-scale GWAS analysis. The pipeline automatically performs numerous pre- and post-processing steps, integrates regression modeling from the REGENIE package and supports single-variant, gene-based and interaction testing. It includes an extensive reporting functionality that allows to inspect thousands of phenotypes and navigate interactive Manhattan plots directly in the web browser. The pipeline is tested using the unit-style testing framework nf-test, a crucial requirement in clinical and pharmaceutical settings. Furthermore, we validated the pipeline against published GWAS datasets and benchmarked the pipeline on high-performance computing and cloud infrastructures to provide cost estimations to end users. nf-gwas is a highly parallelized, scalable and well-tested Nextflow pipeline to perform GWAS analysis in a reproducible manner.

A Genome-Wide Association Study of Metabolic Syndrome in the Taiwanese Population.

The purpose of this study was to investigate genetic factors associated with metabolic syndrome (MetS) by conducting a large-scale genome-wide association study (GWAS) in Taiwan, addressing the limited data on Asian populations compared to Western populations. Using data from the Taiwan Biobank, comprehensive clinical and genetic information from 107,230 Taiwanese individuals was analyzed. Genotyping data from the TWB1.0 and TWB2.0 chips, including over 650,000 single nucleotide polymorphisms (SNPs), were utilized. Genotype imputation using the 1000 Genomes Project was performed, resulting in more than 9 million SNPs. MetS was defined based on a modified version of the Adult Treatment Panel III criteria. Among all participants (mean age: 50 years), 23% met the MetS definition. GWAS analysis identified 549 SNPs significantly associated with MetS, collectively mapping to 10 genomic risk loci. Notable risk loci included rs1004558, rs3812316, rs326, rs4486200, rs2954038, rs10830963, rs662799, rs62033400, rs183130, and rs34342646. Gene-set analysis revealed 22 associated genes: CETP, LPL, APOA5, SIK3, ZPR1, APOC1, BUD13, MLXIPL, TOMM40, GCK, YKT6, RPS6KB1, FTO, VMP1, TUBD1, BCL7B, C19orf80 (ANGPTL8), SIDT2, SENP7, PAFAH1B2, DOCK6, and FOXA2. This study identified genomic risk loci for MetS in a large Taiwanese population through a comprehensive GWAS approach. These associations provide novel insights into the genetic basis of MetS and hold promise for the potential discovery of clinical biomarkers.

Single and longitudinal genome-wide association studies for dairy traits available in goats with three recorded lactations.

Milk yield and composition phenotypes are systematically recorded across several lactations in goats, but the majority of genome-wide association studies (GWAS) performed so far have rather ignored the longitudinal nature of such data. Here, we have used two different GWAS approaches to analyse data from three lactations recorded in Murciano-Granadina goats. In Analysis 1, independent GWAS have been carried out for each trait and lactation, while a single longitudinal GWAS, jointly considering all data, has been performed in Analysis 2. In both analyses, genome-wide significant QTL for lactose percentage on chromosome 2 (129.77-131.01 Mb) and for milk protein percentage on the chromosome 6 (74.8-94.6 Mb) casein gene cluster region were detected. In Analysis 1, several QTL were not replicated in all three lactations, possibly due to the existence of lactation-specific genetic determinants. In Analysis 2, we identified several genome-wide significant QTL related to milk yield and protein content that were not uncovered in Analysis 1. The increased number of QTL identified in Analysis 2 suggests that the longitudinal GWAS is particularly well suited for the genetic analysis of dairy traits. Moreover, our data confirm that variability within or close to the casein complex is the main genetic determinant of milk protein percentage in Murciano-Granadina goats.

Not just females and males: Unravelling the complex sex determinism of the hemp palm, Trachycarpus fortunei.

The ornamental Asian palm Trachycarpus fortunei (Arecaceae: Coryphoideae) is widely planted in temperate regions. In Europe, it has spread outside of gardens, particularly on the southern side of the Alps. Sexual expression in the species is complex, varying from dioecy to polygamy. This study investigated (1) sexual floral development and (2) genetic markers implicated in sex determinism. The morphology and anatomy of floral organs at different developmental stages were studied using SEM observations and anatomical section. Sex determinism was explored using a genome-wide association study approach, searching for correlations between 31,000 single-nucleotide polymorphisms and sex affiliation of 122 palms from 21 wild populations. We observed that sexual differentiation appears late in floral development of T. fortunei. Morpho-anatomical characters of flowers conducive to panmixia were observed, such as well-differentiated septal nectaries that are thought to promote cross-pollination. At the molecular level, homozygous and heterozygous allelic systems with closely linked regions were found for sex determinism in individuals with female and "dominant-male" phenotypes, respectively. Through our wide sampling in the southern Alps, the closely linked genetic regions in males suggest that at least fifteen percent of wild palms are the direct offspring of "males" that can also produce fertile pistillate flowers. Trachycarpus fortunei is a further example of unstable sexual expression found in the family Arecaceae and represents an evolutionary path towards an XY genetic system. Our structural and genetic results may explain the high species dispersal ability in the southern Alps.

Genome-Wide Association Studies (GWAS) Approaches for the Detection of Genetic Variants Associated with Antibiotic Resistance: A Systematic Review.

Antibiotic resistance is a significant threat to public health worldwide. Genome-wide association studies (GWAS) have emerged as a powerful tool to identify genetic variants associated with this antibiotic resistance. By analyzing large datasets of bacterial genomes, GWAS can provide valuable insights into the resistance mechanisms and facilitate the discovery of new drug targets. The present study aimed to undertake a systematic review of different GWAS approaches used for detecting genetic variants associated with antibiotic resistance. We comprehensively searched the PubMed and Scopus databases to identify relevant studies published from 2013 to February 2023. A total of 40 studies met our inclusion criteria. These studies explored a wide range of bacterial species, antibiotics, and study designs. Notably, most of the studies were centered around human pathogens such as Mycobacterium tuberculosis, Escherichia coli, Neisseria gonorrhoeae, and Staphylococcus aureus. The review seeks to explore the several GWAS approaches utilized to investigate the genetic mechanisms associated with antibiotic resistance. Furthermore, it examines the contributions of GWAS approaches in identifying resistance-associated genetic variants through binary and continuous phenotypes. Overall, GWAS holds great potential to enhance our understanding of bacterial resistance and improve strategies to combat infectious diseases.

KGWASflow: a modular, flexible, and reproducible Snakemake workflow for k-mers-based GWAS.

Genome-wide association studies (GWAS) have been widely used to identify genetic variation associated with complex traits. Despite its success and popularity, the traditional GWAS approach comes with a variety of limitations. For this reason, newer methods for GWAS have been developed, including the use of pan-genomes instead of a reference genome and the utilization of markers beyond single-nucleotide polymorphisms, such as structural variations and k-mers. The k-mers-based GWAS approach has especially gained attention from researchers in recent years. However, these new methodologies can be complicated and challenging to implement. Here, we present kGWASflow, a modular, user-friendly, and scalable workflow to perform GWAS using k-mers. We adopted an existing kmersGWAS method into an easier and more accessible workflow using management tools like Snakemake and Conda and eliminated the challenges caused by missing dependencies and version conflicts. kGWASflow increases the reproducibility of the kmersGWAS method by automating each step with Snakemake and using containerization tools like Docker. The workflow encompasses supplemental components such as quality control, read-trimming procedures, and generating summary statistics. kGWASflow also offers post-GWAS analysis options to identify the genomic location and context of trait-associated k-mers. kGWASflow can be applied to any organism and requires minimal programming skills. kGWASflow is freely available on GitHub (https://github.com/akcorut/kGWASflow) and Bioconda (https://anaconda.org/bioconda/kgwasflow).

Combined reference-free and multi-reference based GWAS uncover cryptic variation underlying rapid adaptation in a fungal plant pathogen.

Microbial pathogens often harbor substantial functional diversity driven by structural genetic variation. Rapid adaptation from such standing variation threatens global food security and human health. Genome-wide association studies (GWAS) provide a powerful approach to identify genetic variants underlying recent pathogen adaptation. However, the reliance on single reference genomes and single nucleotide polymorphisms (SNPs) obscures the true extent of adaptive genetic variation. Here, we show quantitatively how a combination of multiple reference genomes and reference-free approaches captures substantially more relevant genetic variation compared to single reference mapping. We performed reference-genome based association mapping across 19 reference-quality genomes covering the diversity of the species. We contrasted the results with a reference-free (i.e., k-mer) approach using raw whole-genome sequencing data in a panel of 145 strains collected across the global distribution range of the fungal wheat pathogen Zymoseptoria tritici. We mapped the genetic architecture of 49 life history traits including virulence, reproduction and growth in multiple stressful environments. The inclusion of additional reference genome SNP datasets provides a nearly linear increase in additional loci mapped through GWAS. Variants detected through the k-mer approach explained a higher proportion of phenotypic variation than a reference genome-based approach and revealed functionally confirmed loci that classic GWAS approaches failed to map. The power of GWAS in microbial pathogens can be significantly enhanced by comprehensively capturing structural genetic variation. Our approach is generalizable to a large number of species and will uncover novel mechanisms driving rapid adaptation of pathogens.

Dissecting genomic regions and underlying sheath blight resistance traits in rice (Oryza sativa L.) using a genome-wide association study.

The productivity of rice is greatly affected by the infection of the plant pathogenic fungus Rhizoctonia solani, which causes a significant grain yield reduction globally. There exist a limited number of rice accessions that are available to develop sheath blight resistance (ShB). Our objective was to identify a good source of the ShB resistance, understand the heritability, and trait interactions, and identify the genomic regions for ShB resistance traits by genome-wide association studies (GWAS). In the present study, a set of 330 traditional landraces and improved rice varieties were evaluated for ShB resistance and created a core panel of 192 accessions used in the GWAS. This panel provides a more considerable amount of genetic variance and found a significant phenotypic variation among the panel of rice accessions for all the agro-morphological and disease-resistance traits over the seasons. The infection rate of ShB and disease reaction were calculated as percent disease index (PDI) and area under the disease progress curve (AUDPC). The correlation analysis showed a significant positive association between PDIs and AUPDC and a negative association between PDI and plant height, flag leaf length, and grain yield. The panel was genotyped with 133 SSR microsatellite markers, resulting in a genome coverage of 314.83 Mb, and the average distance between markers is 2.53 Mb. By employing GLM and MLM (Q + K) models, 30 marker-trait associations (MTAs) were identified with targeted traits over the seasons. Among these QTLs, eight were found to be novel and located on 2, 4, 8, 10, and 12 chromosomes, which explained the phenotypic variation ranging from 5% to 15%. With the GWAS approach, six candidate genes were identified. Os05t0566400, Os08t0155900, and Os09t0567300 were found to be associated with defense mechanisms against ShB. These findings provided insights into the novel donors of IC283139, IC 277248, Sivappuchithirai Kar, and Bowalia. The promising genomic regions on 10 of 12 chromosomes associated with ShB would be useful in developing rice varieties with durable disease resistance.

Genome resources for the elite bread wheat cultivar Aikang 58 and mining of elite homeologous haplotypes for accelerating wheat improvement

Despite recent progress in crop genomics studies, the genomic changes brought about by modern breeding selection are still poorly understood, thus hampering genomics-assisted breeding, especially in polyploid crops with compound genomes such as common wheat (Triticum aestivum). In this work, we constructed genome resources for the modern elite common wheat variety Aikang 58 (AK58). Comparative genomics between AK58 and the landrace cultivar Chinese Spring shed light on genomic changes that occurred through recent varietal improvement. We also explored subgenome diploidization and divergence in common wheat and developed a homoeologous locus–based genome-wide association study (HGWAS) approach, which was more effective than single homoeolog–based GWAS in unraveling agronomic trait–associated loci. A total of 123 major HGWAS loci were detected using a genetic population derived from AK58 and Chinese Spring. Elite homoeologous haplotypes, formed by combinations of subgenomic homoeologs of the associated loci, were found in both parents and progeny, and many could substantially improve wheat yield and related traits. We built a website (available at https://triticeae.henau.edu.cn/aikang58/) where users can download genome assembly sequence and annotation data for AK58, perform blast analysis, and run JBrowse. Our work enriches genome resources for wheat, provides new insights into genomic changes during modern wheat improvement, and suggests that efficient mining of elite homoeologous haplotypes can make a substantial contribution to genomics-assisted breeding in common wheat and other polyploid crops.

Genome-Wide Association Study and Genomic Prediction of Fusarium Wilt Resistance in Common Bean Core Collection.

The common bean (Phaseolus vulgaris L.) is a globally cultivated leguminous crop. Fusarium wilt (FW), caused by Fusarium oxysporum f. sp. phaseoli (Fop), is a significant disease leading to substantial yield loss in common beans. Disease-resistant cultivars are recommended to counteract this. The objective of this investigation was to identify single nucleotide polymorphism (SNP) markers associated with FW resistance and to pinpoint potential resistant common bean accessions within a core collection, utilizing a panel of 157 accessions through the Genome-wide association study (GWAS) approach with TASSEL 5 and GAPIT 3. Phenotypes for Fop race 1 and race 4 were matched with genotypic data from 4740 SNPs of BARCBean6K_3 Infinium Bea Chips. After ranking the 157-accession panel and revealing 21 Fusarium wilt-resistant accessions, the GWAS pinpointed 16 SNPs on chromosomes Pv04, Pv05, Pv07, Pv8, and Pv09 linked to Fop race 1 resistance, 23 SNPs on chromosomes Pv03, Pv04, Pv05, Pv07, Pv09, Pv10, and Pv11 associated with Fop race 4 resistance, and 7 SNPs on chromosomes Pv04 and Pv09 correlated with both Fop race 1 and race 4 resistances. Furthermore, within a 30 kb flanking region of these associated SNPs, a total of 17 candidate genes were identified. Some of these genes were annotated as classical disease resistance protein/enzymes, including NB-ARC domain proteins, Leucine-rich repeat protein kinase family proteins, zinc finger family proteins, P-loopcontaining nucleoside triphosphate hydrolase superfamily, etc. Genomic prediction (GP) accuracy for Fop race resistances ranged from 0.26 to 0.55. This study advanced common bean genetic enhancement through marker-assisted selection (MAS) and genomic selection (GS) strategies, paving the way for improved Fop resistance.