Genome-wide Association Study Analysis Research Articles

A novel locus in CSMD1 gene is associated with increased susceptibility to severe malaria in Malian children

BackgroundPlasmodium falciparum malaria is still a leading cause of child mortality in sub-Saharan Africa. The clinical manifestations of malaria range from asymptomatic infection to severe disease. The variation in clinical presentation is partly attributed to host genetic factors with estimated narrow-sense heritability of 23%. Here, we investigate the associations between candidate gene polymorphisms and the likelihood of severe malaria (SM) in a cohort of Malian children.MethodsBased on our previous genome-wide association studies (GWAS) analysis, candidate genes were selected for in-depth analysis using several criteria including gene-level GWAS scores, functional overlap with malaria pathogenesis, and evidence of association with protection or susceptibility to other infectious or inflammatory diseases. Single Nucleotide Polymorphisms (SNPs) residing within these genes were selected mainly based on p-values from previous severe malaria susceptibility GWAS studies and minor allele frequency (MAF) in West African populations.ResultsOf 182 candidate genes reported in our previous study, 11 genes and 22 SNPs residing in these genes were selected. The selected SNPs were genotyped using KASP technology in 477 DNA samples (87 SM and 390 controls). Logistic regression analysis revealed that a common intron variant, rs13340578 in CUB and Sushi Multi Domain (CSMD1) gene, is associated with increased odds of SM in recessive mode of inheritance (MAF = 0.42, OR = 1.8, 95% CI = [1.78, 1.84], p = 0.029). The SNP is in linkage disequilibrium (LD) with multiple variants with regulatory features.ConclusionTaken together, the current study showed that an intron variant rs13340578, residing in CSMD1 gene, is associated with increased susceptibility to malaria. This finding suggests that modified regulation of complement may contribute to malaria disease severity. Further studies are needed to identify the causal variants and the underlying molecular mechanisms.

Genome-wide association study for seedling heat tolerance under two temperature conditions in bread wheat (Triticum aestivum L.)

BackgroundAs the greenhouse effect intensifies, global temperatures are steadily increasing, posing a challenge to bread wheat (Triticum aestivum L.) production. It is imperative to comprehend the mechanism of high temperature tolerance in wheat and implement breeding programs to identify and develop heat-tolerant wheat germplasm and cultivars.ResultsTo identify quantitative trait loci (QTL) related to heat stress tolerance (HST) at seedling stage in wheat, a panel of 253 wheat accessions which were re-sequenced used to conduct genome-wide association studies (GWAS) using the factored spectrally transformed linear mixed models (FaST-LMM). For most accessions, the growth of seedlings was found to be inhibited under heat stress. Analysis of the phenotypic data revealed that under heat stress conditions, the main root length, total root length, and shoot length of seedlings decreased by 47.46%, 49.29%, and 15.19%, respectively, compared to those in normal conditions. However, 17 varieties were identified as heat stress tolerant germplasm. Through GWAS analysis, a total of 115 QTLs were detected under both heat stress and normal conditions. Furthermore, 15 stable QTL-clusters associated with heat response were identified. By combining gene expression, haplotype analysis, and gene annotation information within the physical intervals of the 15 QTL-clusters, two novel candidate genes, TraesCS4B03G0152700/TaWRKY74-B and TraesCS4B03G0501400/TaSnRK3.15-B, were responsive to temperature and identified as potential regulators of HST in wheat at the seedling stage.ConclusionsThis study conducted a detailed genetic analysis and successfully identified two genes potentially associated with HST in wheat at the seedling stage, laying a foundation to further dissect the regulatory mechanism underlying HST in wheat under high temperature conditions. Our finding could serve as genomic landmarks for wheat breeding aimed at improving adaptation to heat stress in the face of climate change.

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.

GWAS analysis reveals the genetic basis of blast resistance associated with heading date in rice

Rice blast is a destructive fungal disease affecting rice plants at various growth stages, significantly threatening global yield stability. Development of resistant rice cultivars stands as a practical means of disease control. Generally, association mapping with a diversity panel powerfully identifies new alleles controlling trait of interest. On the other hand, utilization of a breeding panel has its advantage that can be directly applied in a breeding program. In this study, we conducted a genome-wide association study (GWAS) for blast resistance using 296 commercial rice cultivars with low population structure but large phenotypic diversity. We attempt to answer the genetic basis behind rice blast resistance among early maturing cultivars by subdividing the population based on its Heading date 1 (Hd1) functionality. Subpopulation-specific GWAS using the mixed linear model (MLM) based on blast nursery screening conducted in three years revealed a total of 26 significant signals, including three nucleotide-binding site leucine-rich repeat (NBS-LRR) genes (Os06g0286500, Os06g0286700, and Os06g0287500) located at Piz locus on chromosome 6, and one at the Pi-ta locus (Os12g0281300) on chromosome 12. Haplotype analysis revealed blast resistance associated with Piz locus was exclusively specific to Type 14 hd1 among japonica rice. Our findings provide valuable insights for breeding blast resistant rice and highlight the applicability of our elite cultivar panel to detect superior alleles associated with important agronomic traits.

Assessment of genetic and clinical factors in T2D susceptibility among patients with hypertension.

Hypertension (HTN) and Type 2 Diabetes (T2D) often coexist, therefore understanding the relationship between both diseases is imperative to guide targeted prevention/therapy. This study aims to explore the relationship between HTN and T2D using genome-wide association study (GWAS) analysis and biochemical data to understand the implication of both clinical and genetic factors in these pathologies. A total of 2,876 patients were enrolled. Using GWAS and biochemical data, patients with both T2D and HTN were compared to patients with only HTN. Specificity was confirmed by testing the detected genetic variants for associations with HTN development in T2D patients, or with HTN in healthy subjects. Regression models were applied to examine the association of T2D in patients with HTN with cardiovascular risk factors. Replication was performed using UK Biobank dataset with 31,170 subjects. Data showed that females with HTN are at higher risk of developing T2D due to dyslipidemia, while males faced higher risk due to high BMI (body mass index) and family history of T2D. GWAS identified Single Nucleotide Polymorphisms (SNPs) linked to T2D in patients with HTN. Notably, rs7865889, rs7756992, and rs10896290 were positively associated with T2D, whereas rs12737517 yielded negative association. Three SNPs were replicated in the UK Biobank (rs10896290, rs7865889, and rs7756992). Incorporating clinical and genetic screening into risk assessment is important for the detection and prevention of T2D in patients with HTN. The detected SNPs (rs7865889, rs12737517, and rs10896290), especially the protective SNP (rs12737517), provide an opportunity for better diagnosis, prevention, and therapy of patients with T2D and HTN.

Genome-wide association study and transcriptome analysis reveal natural variation of key genes regulation flowering time in rapeseed

Genome-wide association study and transcriptome analysis reveal natural variation of key genes regulation flowering time in rapeseed

Genome-wide association analysis explores the genetic loci of amino acid content in duck’s breast muscle

BackgroundAmino acids are the basic components of protein and an important index to evaluate meat quality. With the rapid development of genomics, candidate regions and genes affecting amino acid content in livestock and poultry have been gradually revealed. Hence, genome-wide association study (GWAS) can be used to screen candidate loci associated with amino acid content in duck meat.ResultIn the current study, the content of 16 amino acids was detected in 358 duck breast muscles. The proportion of Glu to the total amino acid content was relatively high, and the proportion was 0.14. However, the proportion of Met content was relatively low, at just 0.03. By comparative analysis, significant differences were found between males and females in 3 amino acids, including Ser, Met, and Phe. In addition, 12 SNPs were significantly correlated with Pro content by GWAS analysis, and these SNPs were annotated by 7 protein-coding genes; 8 significant SNPs were associated with Tyr content, and these SNPs were annotated by 6 protein-coding genes. At the same time, linkage disequilibrium (LD) analysis was performed on these regions with significant signals. The results showed that three SNPs in the 55–56 Mbp region of chromosome 3 were highly correlated with the leader SNP (chr3:55526954) that affected Pro content (r2 > 0.6). Similarly, LD analysis showed that there were three SNPs in the 21.2–21.6 Mbp region of chromosome 13, which were highly correlated with leader SNP (chr13:21421661) (r2 > 0.6). Moreover, Through functional enrichment analysis of all candidate genes. The results of GO enrichment analysis showed that several significant GO items were associated with amino acid transport function, including amino acid transmembrane transport and glutamine transport. The results further indicate that these candidate genes are closely associated with amino acid transport. Among them, key candidate genes include SLC38A1. For KEGG enrichment analysis, CACNA2D3 and CACNA1D genes were covered by significant pathways.ConclusionIn this study, GWAS analysis found a total of 28 significant SNPs affecting amino acid content. Through gene annotation, a total of 20 candidate genes were screened. In addition, Through LD analysis and enrichment analysis, we considered that SERAC1, CACNA2D3 and SLC38A1 genes are important candidate genes affecting amino acid content in duck breast muscle.

GWAS and WGCNA Analysis Uncover Candidate Genes Associated with Oil Content in Soybean.

Soybean vegetable oil is an important source of the human diet. However, the analysis of the genetic mechanism leading to changes in soybean oil content is still incomplete. In this study, a total of 227 soybean materials were applied and analyzed by a genome-wide association study (GWAS). There are 44 quantitative trait nucleotides (QTNs) that were identified as associated with oil content. A total of six, four, and 34 significant QTN loci were identified in Xiangyang, Hulan, and Acheng, respectively. Of those, 26 QTNs overlapped with or were near the known oil content quantitative trait locus (QTL), and 18 new QTNs related to oil content were identified. A total of 594 genes were located near the peak single nucleotide polymorphism (SNP) from three tested environments. These candidate genes exhibited significant enrichment in tropane, piperidine, and pyridine alkaloid biosynthesiss (ko00960), ABC transporters (ko02010), photosynthesis-antenna proteins (ko00196), and betalain biosynthesis (ko00965). Combined with the GWAS and weighted gene co-expression network analysis (WGCNA), four candidate genes (Glyma.18G300100, Glyma.11G221100, Glyma.13G343300, and Glyma.02G166100) that may regulate oil content were identified. In addition, Glyma.18G300100 was divided into two main haplotypes in the studied accessions. The oil content of haplotype 1 is significantly lower than that of haplotype 2. Our research findings provide a theoretical basis for improving the regulatory mechanism of soybean oil content.

Novel insight into the genetic signatures of altitude adaptation related body composition in Tibetans.

The Tibetan population residing in high-altitude (HA) regions has adapted to extreme hypoxic environments. However, there is limited understanding of the genetic basis of body compositions in Tibetan population adapted to HA. We performed a genome-wide association study (GWAS) to identify genetic variants associated with HA and HA-related body composition traits. A total of 755,731 single nucleotide polymorphisms (SNPs) were genotyped using the precision medicine diversity array from 996 Tibetan college students. T-tests and Pearson correlation analysis were used to estimate the association between body compositions and altitude. The mixed linear regression identified the SNPs significantly associated with HA and HA-related body compositions. LASSO regression was used to screen for important SNPs in HA and body compositions. Significant differences were observed in lean body mass (LBW), muscle mass (MM), total body water (TBW), standard weight (SBW), basal metabolic rate (BMR), total protein (TP), and total inorganic salt (Is) in different altitudes stratification. We identified three SNPs in EPAS1 (rs1562453, rs7589621 and rs7583392) that were significantly associated with HA (p < 5 × 10-7). GWAS analysis of 7 HA-related body composition traits, we identified 14 SNPs for LBM, 11 SNPs for TBW, 15 SNPs for MM, 16 SNPs for SBW, 9 SNPs for BMR, 12 SNPs for TP, and 26 SNPs for Is (p < 5.0 × 10-5). These findings provide insight into the genetic basis of body composition in Tibetan college students adapted to HA, and lay the foundation for further investigation into the molecular mechanisms underlying HA adaptation.

Genomic and transcriptome analyses reveal potential contributors to erucic acid biosynthesis in seeds of rapeseed (Brassica napus).

Through comprehensive genomic and transcriptomic analyses, we identified a set of 23 genes that act up- or downstream of erucic acid content (EAC) production in rapeseed seeds. We selected example genes to showcase the distribution of single nucleotide polymorphisms, haplotypes associated with EAC phenotypes, and the creation of molecular markers differentiating low EAC and high EAC genotypes. Erucic acid content (EAC) is a crucial trait in rapeseed, with low LEAC oil recognized for its health benefits and high EA oil holding industrial value. Despite its significance, the genomic consequences of intensive LEAC-cultivar selection and the genetic basis underlying EA regulation remain largely unexplored. To address this knowledge gap, we conducted selective signal analyses, genome-wide association studies (GWAS), and transcriptome analyses. Our investigation unveiled the genetic footprints resulting from LEAC selection in germplasm populations, drawing attention to specific loci that contribute to enriching diversity. By integrating GWAS and transcriptome analyses, we identified a set of 23 genes that play a significant role in determining EAC in seeds or are downstream consequences of EA-level alterations. These genes have emerged as promising candidates for elucidating the potential mechanisms governing EAC in rapeseed. To exemplify the findings, we selected specific genes to demonstrate the distribution of single nucleotide polymorphisms and haplotypes associated with different EAC phenotypes. Additionally, we showcased to develop molecular markers distinguishing between LEAC and high EAC genotypes.

Assessing the causality of interferon‐γ and its receptor 1/2 with systemic lupus erythematosus risk using genetic data

AbstractBackgroundThe interferon‐γ (IFN‐γ) signaling pathway is activated in systemic lupus erythematosus (SLE). This study aimed to assess the causal association between IFN‐γ, IFN‐γ receptor 1 (IFN‐γR1), and IFN‐γR2 and SLE using a bidirectional Mendelian randomization (MR) design.MethodsGenetic instruments for exposure to IFN‐γ, IFN‐γR1, and IFN‐γR2 were derived from a large genome‐wide association study (GWAS) that included a sample size of 3301 participants. Instrumental variables for SLE were selected from another independent GWAS analysis comprising 5201 cases and 6099 controls with European ancestry. Bidirectional two‐sample MR was performed using inverse variance weighting, MR‐Egger regression, and weighted median methods. A series of sensitivity analyses were conducted to assess the robustness of the results.ResultsThe inverse variance weighting showed that IFN‐γ had a positive causal association with the risk of SLE (odd ratio [OR] = 1.24, 95% confidence interval [CI]: 1.03–1.47, p = 0.018). IFN‐γR2 levels were not associated with SLE risk after adjustment for multiple comparisons (OR = 0.85, 95% CI: 0.73–0.99, p = 0.034). No genetic association was also detected between IFN‐γR1 and SLE (OR = 0.97, 95% CI: 0.79–1.19, p = 0.768). Evidence from bidirectional MR did not support reverse causality. The weighted median regression also showed directionally similar estimates.ConclusionHigher levels of IFN‐γ are significantly associated with an increased risk of SLE, providing insights into the pathogenesis of SLE.

Genome-Wide Epistatic Network Analyses of Semantic Fluency in Older Adults.

Semantic fluency impairment has been attributed to a wide range of neurocognitive and psychiatric conditions, especially in the older population. Moderate heritability estimates on semantic fluency were obtained from both twin and family-based studies suggesting genetic contributions to the observed variation across individuals. Currently, effort in identifying the genetic variants underlying the heritability estimates for this complex trait remains scarce. Using the semantic fluency scale and genome-wide SNP genotype data from the Long Life Family Study (LLFS), we performed a genome-wide association study (GWAS) and epistasis network analysis on semantic fluency in 2289 individuals aged over 60 years from the American LLFS cohorts and replicated the findings in 1129 individuals aged over 50 years from the Danish LLFS cohort. In the GWAS, two SNPs with genome-wide significance (rs3749683, p = 2.52 × 10-8; rs880179, p = 4.83 × 10-8) mapped to the CMYAS gene on chromosome 5 were detected. The epistasis network analysis identified five modules as significant (4.16 × 10-5 < p < 7.35 × 10-3), of which two were replicated (p < 3.10 × 10-3). These two modules revealed significant enrichment of tissue-specific gene expression in brain tissues and high enrichment of GWAS catalog traits, e.g., obesity-related traits, blood pressure, chronotype, sleep duration, and brain structure, that have been reported to associate with verbal performance in epidemiological studies. Our results suggest high tissue specificity of genetic regulation of gene expression in brain tissues with epistatic SNP networks functioning jointly in modifying individual verbal ability and cognitive performance.

Large-scale transcriptomic and genomic analyses reveal a novel functional gene SERPINB6 for chicken carcass traits

BackgroundCarcass traits are crucial indicators of meat production efficiency. However, the molecular regulatory mechanisms associated with these traits remain unclear.ResultsIn this study, we conducted comprehensive transcriptomic and genomic analyses on 399 Tiannong partridge chickens to identify key genes and variants associated with carcass traits and to elucidate the underlying regulatory mechanisms. Based on association analyses with the elastic net (EN) model, we identified 12 candidate genes (AMY1A, AP3B2, CEBPG, EEF2, EIF4EBP1, FGFR1, FOXD3, GOLM1, LOC107052698, PABPC1, SERPINB6 and TBC1D16) for 4 carcass-related traits, namely live weight, dressed weight, eviscerated weight, and breast muscle weight. SERPINB6 was identified as the only overlapping gene by 3 analyses, EN model analysis, weighted gene co-expression network analysis and differential expression analysis. Cell-level experiments confirmed that SERPINB6 promotes the proliferation of chicken DF1 cells and primary myoblasts. Further expression genome-wide association study and association analysis indicated that rs317934171 is the critical site that enhances SERPINB6 expression. Furthermore, a dual-luciferase reporter assay proved that gga-miR-1615 targets the 3′UTR of SERPINB6.ConclusionsCollectively, our findings reveal that SERPINB6 serves as a novel gene for chicken carcass traits by promoting fibroblast and myoblast proliferation. Additionally, the downstream variant rs317934171 regulates SERPINB6 expression. These results identify a new target gene and molecular marker for the molecular mechanisms of chicken carcass traits.

Genetic Diversity and Genome-Wide Association Analysis of the Hulled/Naked Trait in a Barley Collection from Shanghai Agricultural Gene Bank.

Barley is one of the most important cereal crops in the world, and its value as a food is constantly being revealed, so the research into and the use of barley germplasm are very important for global food security. Although a large number of barley germplasm samples have been collected globally, their specific genetic compositions are not well understood, and in many cases their origins are even disputed. In this study, 183 barley germplasm samples from the Shanghai Agricultural Gene Bank were genotyped using genotyping-by-sequencing (GBS) technology, SNPs were identified and their genetic parameters were estimated, principal component analysis (PCA) was preformed, and the phylogenetic tree and population structure of the samples were also analyzed. In addition, a genome-wide association study (GWAS) was carried out for the hulled/naked grain trait, and a KASP marker was developed using an associated SNP. The results showed that a total of 181,906 SNPs were identified, and these barley germplasm samples could be roughly divided into three categories according to the phylogenetic analysis, which was generally consistent with the classification of the traits of row type and hulled/naked grain. Population structure analysis showed that the whole barley population could be divided into four sub-populations (SPs), the main difference from previous classifications being that the two-rowed and the hulled genotypes were sub-divided into two SPs. The GWAS analysis of the hulled/naked trait showed that many associated loci were unrelated to the Nud/nud locus, indicating that there might be new loci controlling the trait. A KASP marker was developed for one exon-type SNP on chromosome 7. Genotyping based on the KASP assay was consistent with that based on SNPs, indicating that the gene of this locus might be associated with the hulled/naked trait. The above work not only lays a good foundation for the future utilization of this barley germplasm population but it provides new loci and candidate genes for the hulled/naked trait.

Metalloproteins as risk factors for osteoarthritis: improving and understanding causal estimates using Mendelian randomization.

Osteoarthritis (OA) is one of the most prevalent musculoskeletal disorders and a primary cause of pain and disability among the elderly population. Research on the relationship between metalloproteins (MPs) and OA is limited, and causality remains unclear. Our objective is to utilize Mendelian randomization (MR) to explore the possible causal relationship between MPs and OA. The data on MPs were derived from a Genome-Wide Association Study (GWAS) analysis involving 3301 samples. The GWAS data for OA were obtained from an analysis involving 462,933 European individuals. In this study, a variety of two-sample Mendelian randomization methods (two-sample MR) to evaluate the causal effect of MPs on OA, including inverse variance weighted method (IVW), MR-Egger method, weighted median method (WM), simple mode, weight mode, and Wald ratio. The primary MR analysis using the IVW method reveals a significant negative correlation between Metallothionein-1F (MT-1F), zinc finger protein 134 (ZNF134), calcium/calmodulin-dependent protein kinase type 1D (CAMK1D), and EF-hand calcium-binding domain-containing protein 14 (EFCAB14) with the occurrence of osteoarthritis (OA) (p value < 0.05). However, no causal relationship was observed in the opposite direction between these MPs and OA. Notably, even in combined models accounting for confounding factors, the negative association between these four MPs and OA remained significant. Sensitivity analysis demonstrated no evidence of horizontal pleiotropy or heterogeneity, and leave-one-out analysis confirmed the robustness of the results. In this study, we have established a conspicuous association between four distinct MPs and OA. This discovery augments our understanding of potential avenues for the diagnosis and treatment of this condition. Key Points • The MR method was employed to assess the relationship between MPs and OA. • A total of four types of MPs have demonstrated inhibitory effects on the occurrence of OA.

PLA2G12A protects against diet-induced obesity and insulin resistance by enhancing energy expenditure and clearance of circulating triglycerides.

Secreted phospholipase A2s are involved in the development of obesity, type 2 diabetes mellitus (T2DM) and cardiovascular disease, which have become serious and growing health concerns worldwide. Integration of genome-wide association study and gene co-expression networks analysis showed that the secreted phospholipase A2 group XIIA (PLA2G12A) may participate in hepatic lipids metabolism. Nevertheless, the role of PLA2G12A in lipid metabolism and its potential mechanism remain elusive. Here, we used AAV9 vector carrying human PLA2G12A gene to exogenously express hPLA2G12A in the liver of mice. We demonstrated that the overexpression of hPLA2G12A resulted in a significant decrease in serum lipid levels in wild-type mice fed with chow diet or high-fat diet (HFD). Moreover, hPLA2G12A treatment protected against diet-induced obesity and insulin resistance in mice fed a HFD. Notably, we found that hPLA2G12A treatment confers protection against obesity and hyperlipidemia independent of its enzymatic activity, but rather by increasing physical activity and energy expenditure. Furthermore, we demonstrated that hPLA2G12A treatment induced upregulation of ApoC2 and Cd36 and downregulation of Angptl8, which contributed to the increase in clearance of circulating triglycerides and hepatic uptake of fatty acids without affecting hepatic de novo lipogenesis, very low-density lipoprotein secretion, or intestinal lipid absorption. Our study highlights the potential of PLA2G12A gene therapy as a promising approach for treating obesity, insulin resistance and T2DM.

Discovering new genes for alfalfa (Medicago sativa) growth and biomass resilience in combined salinity and Phoma medicaginis infection through GWAS.

Salinity and Phoma medicaginis infection represent significant challenges for alfalfa cultivation in South Africa, Europe, Australia, and, particularly, Tunisia. These constraints have a severe impact on both yield and quality. The primary aim of this study was to establish the genetic basis of traits associated with biomass and growth of 129 Medicago sativa genotypes through genome-wide association studies (GWAS) under combined salt and P. medicaginis infection stresses. The results of the analysis of variance (ANOVA) indicated that the variation in these traits could be primarily attributed to genotype effects. Among the test genotypes, the length of the main stem, the number of ramifications, the number of chlorotic leaves, and the aerial fresh weight exhibited the most significant variation. The broad-sense heritability (H²) was relatively high for most of the assessed traits, primarily due to genetic factors. Cluster analysis, applied to morpho-physiological traits under the combined stresses, revealed three major groups of accessions. Subsequently, a GWAS analysis was conducted to validate significant associations between 54,866 SNP-filtered single-nucleotide polymorphisms (SNPs) and seven traits. The study identified 27 SNPs that were significantly associated with the following traits: number of healthy leaves (two SNPs), number of chlorotic leaves (five SNPs), number of infected necrotic leaves (three SNPs), aerial fresh weight (six SNPs), aerial dry weight (nine SNPs), number of ramifications (one SNP), and length of the main stem (one SNP). Some of these markers are related to the ionic transporters, cell membrane rigidity (related to salinity tolerance), and the NBS_LRR gene family (associated with disease resistance). These findings underscore the potential for selecting alfalfa genotypes with tolerance to the combined constraints of salinity and P. medicaginis infection.

Abstract PS02-05: Trajectories and predictors of peripheral neuropathy after neoadjuvant chemotherapy in a prospective cohort of 11,014 patients with early breast cancer

Abstract Background Peripheral neuropathy (PN) is a debilitating adverse event in patients with early breast cancer (BC) receiving (neo)adjuvant chemotherapy (CT). We harnessed the CANTO cohort study to detail clinical trajectories and explore clinical and genetic predictors of PN. Methods CANTO (CANcer TOxicities - NCT01993498) prospectively enrolled invasive stage I-III BC patients (pts) of 26 French comprehensive cancer centers. Pts were assessed at diagnosis, 3-6 (M6), 12 (M12), 36 (M36), and 60 (M60) months after treatment, defined as completion of surgery, CT or radiotherapy, whichever comes last. At each time point, PN events including paresthesia, sensory and motor neuropathy were collected according to NCI-CTC v4.0 criteria. A genome-wide association study (GWAS) was conducted to identify genetic predictors of PN, using Illumina GSA BeadChips. Minimac4/1000G was used to impute additional single nucleotide polymorphisms (SNPs). After stringent quality control measures, 1,894,475 SNPs with a minor allele frequency (MAF) &amp;gt;0.05 were evaluated. Longitudinal trajectories of PN events were descriptively examined. Statistical associations between each SNP and PN events measured at different time points, and between SNPs and trajectories, were performed with logistic regression assuming a log-additive genetic model. All analyses were adjusted for key clinical parameters and the first ten axes of principal component analysis of the genetic data to control for population stratification. Results Of 12,012 included pts (data lock Aug. 2022), 11,014 (91.7%) were analyzed. Age was &amp;lt; 50 and &amp;gt;65 in 3407 and 2759 pts (31% and 25%, respectively). Overweight/obesity and diabetes were recorded in 5328 (48%) and 458 pts (4%). A neurologic history was observed in 1360 pts (13%). Stage 0-I and II-III were observed in 5356 (49%) and 5512 pts (51%), and 4011 pts had an axillary dissection (37%). CT was administered in 5790 pts (53%), including a taxane (tax) in 5542 pts (96%). Shortly, all grade and grade 2+ PN was observed in 29% and 12%; 27% and 10%; 20% and 10%; 13% and 10% of pts at M6, M12, M36 and M60, respectively. We derived 5 trajectories of overall PN, deemed “never”, “always”, “appears”, “transient” and “disappears”, in respectively 54%, 10%, 9%, 7% and 19% of the pts. Similar trajectories were built for the 3 categories of PN events. We built 9 predictive nomograms including key clinical parameters and time of analysis. E.g., grade 2+ PN at M6 was independently predicted by a medical history of carpal tunnel syndrome, past neurologic history and exposure to tax. Grade 2+ PN at M36 was predicted by past neurologic history and exposure to tax, and largely by previous PN at M6 or M12. The GWAS analysis included 7633 pts (84%). Four independent SNPs adding predictive value to the clinical nomograms were identified at a suggestive level of association (p &amp;lt; 1e-6), including one SNP in the NCAM1 gene (involved in neurogenesis) for M6 PN (OR=1.40, 95%CI 1.23-1.60, p=4.8e-7, see table) and one SNP in the CLDN11 gene (regulates oligodendrocytes) for M12 sensory PN (OR=1.32; 95%CI 1.17-1.49, p=7.5e-6). For rare ( &amp;lt; 5%) toxicities, 10 independent suggestive SNPs were identified, including one SNP in the NELL1 gene (neural cell growth regulation and differentiation) for the “always” trajectory (OR=1.61, 95%CI 1.31-1.99, p=7.6e-6) and one SNP in the KCNIP1 gene (neuronal sensor) for M36 motor PN (OR=1.5, 95%CI 1.7-3.5, p=3.1e-07). Conclusions Risk of early peripheral neuropathy is associated with previous neurological history and taxane exposure. Patients with PN at M6 or M12 are highly exposed to long-term PN. Some key SNPs may add independent predictive value for specific PN endpoints. Exploratory results of interest were produced for rarer toxicities and typical clinical trajectories. Detailed results will be presented at the meeting. All PN, M6 Multivariate analyses of the risk of peripheral neuropathy, all grades, at M6, without and with GWAS input. Citation Format: Youenn Drouet, Emilie Thomas, Florence Lerebours, Barbara Pistilli, Olivier Trédan, Christelle Jouannaud, Marion Fournier, Philippe Rouanet, Laurence Vanlemmens, CHARLES COUTANT, Asma Dhaini Merimeche, Baptiste Sauterey, Christelle Levy, Mario Campone, Carole Tarpin, Marie-Ange Mouret-Reynier, Olivier Rigal, Thierry Petit, Sophie Guillermet, Antoine Arnaud, Mahmoud Ibrahim, Sylvie Giacchetti, Florence Dalenc, Johanna Wassermann, Olivier Arsène, Ariane Darut-Jouve, Sibille Everhard, Ines Vaz Luis, Anne-Laure Martin, Fabrice André, Jean-François Deleuze, Alain Viari, Matthieu Carton, Paul Cottu. Trajectories and predictors of peripheral neuropathy after neoadjuvant chemotherapy in a prospective cohort of 11,014 patients with early breast cancer [abstract]. In: Proceedings of the 2023 San Antonio Breast Cancer Symposium; 2023 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2024;84(9 Suppl):Abstract nr PS02-05.

Identification of novel SNPs and candidate genes significantly affecting growth in grass carp (Ctenopharyngodon idella) through GWAS analysis

Grass carp (Ctenopharyngodon idella) is one of the most important cultivated fish species in the global freshwater aquaculture with fast growth, low cost of production and delicious meat. A breeding program to improve economically important growth traits in grass carp would benefit from effective molecular markers. Although a genetic linkage map of grass carp has been constructed, and quantitative trait loci (QTL) mapping and a genome-wide association study (GWAS) have been conducted, the effective number of SNP markers and genes currently available for breeding is still very limited. In this study, we conducted a GWAS based on re-sequencing of individuals to identify genomic regions and candidate genes potentially associated with various growth traits of grass carp, including body weight (BW), body length (BL) and body height (BH). A total of 4,854,366 high-quality SNPs were identified and located on 24 chromosomes. Among them, 31 SNPs were identified as being associated with growth traits, including 8 SNPs on chromosomes 1, 2, 10, 13, 14, 17 and 22 associated with BW, 21 SNPs on chromosomes 2, 4, 5, 8, 10, 13, 14, 15, 17, 23 and 24 correlated with BL, and 18 SNPs on chromosomes 2, 9, 10, 13, 14, 17, 23 and 24 related to BH. In addition, 13 of these SNPs were simultaneously associated with two or three growth traits. Within a 48 kb range upstream and downstream of these significant SNPs, 120 candidate genes were identified to be associated with the QTLs. Most of these genes were related to growth and development, fat metabolism, bone mineralization and differentiation, vascular development, or feeding regulation. Notably, the genes pwp1, cers5, chrebp, lsr, ctnnbip1, smoc2, ano6 and cckrb were specifically associated with growth, fat, bone formation and feeding. Our findings would provide a foundation for molecular marker-assisted selection, genomic selection, and functional analysis of growth traits in grass carp.

Genetic control of root morphology in response to nitrogen across rapeseed diversity.

Rapeseed (Brassica napus L.) is an oil-containing crop of great economic value but with considerable nitrogen requirement. Breeding root systems that efficiently absorb nitrogen from the soil could be a driver to ensure genetic gains for more sustainable rapeseed production. The aim of this study is to identify genomic regions that regulate root morphology in response to nitrate availability. The natural variability offered by 300 inbred lines was screened at two experimental locations. Seedlings grew hydroponically with low or elevated nitrate levels. Fifteen traits related to biomass production and root morphology were measured. On average across the panel, a low nitrate level increased the root-to-shoot biomass ratio and the lateral root length. A large phenotypic variation was observed, along with important heritability values and genotypic effects, but low genotype-by-nitrogen interactions. Genome-wide association study and bulk segregant analysis were used to identify loci regulating phenotypic traits. The first approach nominated 319 SNPs that were combined into 80 QTLs. Three QTLs identified on the A07 and C07 chromosomes were stable across nitrate levels and/or experimental locations. The second approach involved genotyping two groups of individuals from an experimental F2 population created by crossing two accessions with contrasting lateral root lengths. These individuals were found in the tails of the phenotypic distribution. Co-localized QTLs found in both mapping approaches covered a chromosomal region on the A06 chromosome. The QTL regions contained some genes putatively involved in root organogenesis and represent selection targets for redesigning the root morphology of rapeseed.