Search

Background: DNA methylation is a critical regulatory mechanism of gene expression, influencing various human diseases and traits. While traditional expression quantitative trait loci (eQTL) studies have helped elucidate the genetic regulation of gene expression, there is a growing need to explore environmental influences on gene expression. Existing methods such as PrediXcan and FUSION focus on genotype-based associations but overlook the impact of environmental factors. To address this gap, we present MOSES (methylation-based gene association), a novel approach that utilizes DNA methylation to identify environmentally regulated genes associated with traits or diseases without relying on measured gene expression., Results: MOSES involves training, imputation, and association testing. It employs elastic-net penalized regression models to estimate the influence of CpGs and SNPs (if available) on gene expression. We developed and compared four MOSES versions incorporating different methylation and genetic data: (1) cis-DNA methylation within 1 Mb of promoter regions, (2) both cis-SNPs and cis-CpGs, 3) both cis- and a part of trans- CpGs (±5Mb away) from promoter regions), and 4) long-range DNA methylation (±10 Mb away) from promoter regions. Our analysis using nasal epithelium and white blood cell data from the Epigenetic Variation and Childhood Asthma in Puerto Ricans (EVA-PR) study demonstrated that MOSES, particularly the version incorporating long-range CpGs (MOSES-DNAm 10 M), significantly outperformed existing methods like PrediXcan, MethylXcan, and Biomethyl in predicting gene expression. MOSES-DNAm 10 M identified more differentially expressed genes (DEGs) associated with atopic asthma, particularly those involved in immune pathways, highlighting its superior performance in uncovering environmentally regulated genes. Further application of MOSES to lung tissue data from idiopathic pulmonary fibrosis (IPF) patients confirmed its robustness and versatility across different diseases and tissues., Conclusion: MOSES represents an innovative advancement in gene association studies, leveraging DNA methylation to capture the influence of environmental factors on gene expression. By incorporating long-range CpGs, MOSES-DNAm 10 M provides superior predictive accuracy and gene association capabilities compared to traditional genotype-based methods. This novel approach offers valuable insights into the complex interplay between genetics and the environment, enhancing our understanding of disease mechanisms and potentially guiding therapeutic strategies. The user-friendly MOSES R package is publicly available to advance studies in various diseases, including immune-related conditions like asthma., Competing Interests: Declarations Ethical approval and consent to participate The Epigenetic Variation and Childhood Asthma in Puerto Ricans (EVA-PR) study was approved by the institutional review boards of the University of Puerto Rico (San Juan, PR) and the University of Pittsburgh (Pittsburgh, PA). Written parental consent and participant assent were obtained for individuals younger than 18, while consent was obtained from participants 18 years and older. Consent for publication Not applicable. Competing interests The authors declare no competing interests., (© 2024. The Author(s).)

Main Conclusion: Association analysis identified 77 marker-trait associations (MTAs) for PUE traits, of which 10 were high-confidence MTAs. Candidate-gene mining and in-silico expression analysis identified 13 putative candidate genes for PUE traits. Bread wheat (Triticum aestivum L.) is a major cereal crop affected by phosphorus (P) deficiency, which affects root characteristics, plant biomass, and other attributes related to P-use efficiency (PUE). Understanding the genetic mechanisms of PUE traits helps in developing bread wheat cultivars that perform well in low-P environments. With this objective, we evaluated a bread wheat panel comprising 304 accessions for 14 PUE traits with high-throughput phenotyping under low-P and optimum-P treatments and observed a significant genetic variation among germplasm lines for studied traits. Genome-wide association study (GWAS) using 14,025 high-quality single-nucleotide polymorphisms identified 77 marker-trait associations (MTAs), of which 10 were chosen as high-confidence MTAs as they had > 10% phenotypic variation with logarithm of odds (LOD) scores of more than five. Candidate-gene (CG) mining from high-confidence MTAs identified 180 unique gene models, of which 78 were differentially expressed (DEGs) with at least twofold change in expression under low-P over optimum-P. Of the 78-DEGs, 13 were thought to be putative CGs as they exhibited functional relevance to PUE traits. These CGs mainly encode for important proteins and their products involved in regulating root system architecture, P uptake, transport, and utilization. Promoter analysis from 1500 bp upstream of gene start site for 13 putative CGs revealed the presence of light responsive, salicylic-acid responsive, gibberellic-acid (GA)-responsive, auxin-responsive, and cold responsive cis-regulatory elements. High-confidence MTAs and putative CGs identified in this study can be employed in breeding programs to improve PUE traits in bread wheat., Competing Interests: Declarations Conflict of interest On behalf of all authors, the corresponding author states that there is no conflict of interest. Ethical approval and consent to participate Not applicable. Consent to participate Not applicable. Consent for publication Not applicable., (© 2024. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Background: This study aimed to identify DNA methylation biomarkers associated with myopia using summary-data-based Mendelian randomization (SMR)., Methods: A systematic search of the PubMed, Web of Science, Cochrane Library, and Embase databases was conducted up to March 27, 2024. SMR analyses were performed to integrate genome-wide association study (GWAS) with methylation quantitative trait loci (mQTL) and expression quantitative trait loci (eQTL) studies. The heterogeneity in the dependent instrument (HEIDI) test was utilized to distinguish pleiotropic associations from linkage disequilibrium., Results: The systematic review identified 26 DNA methylation biomarkers in five studies, with no overlap observed among those identified by different studies. After integrating GWAS with multi-omics data of mQTL and eQTL, six genes were significantly associated with myopia: PRMT6 (cg00944433 and cg15468180), SH3YL1 (cg03299269, cg11361895, and cg13354988), ZKSCAN4 (cg01192291), GATS (cg17830204), NPAT (cg04826772), and UBE2I (cg03545757 and cg08025960)., Conclusions: We identified six methylation biomarkers associated with the risk of myopia that may be helpful to elucidate the etiology mechanisms of myopia. Further experimental validation studies are required to corroborate these findings., Competing Interests: Declarations Ethics approval and consent to participate Given that this study represents a reevaluation of previously collected and published data, no additional ethical approval or informed consent from the subjects was needed. Competing interests The authors declare no competing interests., (© 2024. The Author(s).)

Balancing the tradeoff between quantity and quality of phenotypic data is critical in omics studies. Measurements below the limit of quantification (BLQ) are often tagged in quality control fields, but these flags are currently underutilized in human genetics studies. Extreme phenotype sampling is advantageous for mapping rare variant effects. We hypothesize that genetic drivers, along with environmental and technical factors, contribute to the presence of BLQ flags. Here, we introduce "hypometric genetics" (hMG) analysis and uncover a genetic basis for BLQ flags, indicating an additional source of genetic signal for genetic discovery, especially from phenotypic extremes. Applying our hMG approach to n = 227,469 UK Biobank individuals with metabolomic profiles, we reveal more than 5% heritability for BLQ flags and report biologically relevant associations, for example, at APOC3, APOA5, and PDE3B loci. For common variants, polygenic scores trained only for BLQ flags predict the corresponding quantitative traits with 91% accuracy, validating the genetic basis. For rare coding variant associations, we find an asymmetric 65.4% higher enrichment of metabolite-lowering associations for BLQ flags, highlighting the impact of putative loss-of-function variants with large effects on phenotypic extremes. Joint analysis of binarized BLQ flags and the corresponding quantitative metabolite measurements improves power in Bayesian rare variant aggregation tests, resulting in an average of 181% more prioritized genes. Our approach is broadly applicable to omics profiling. Overall, our results underscore the benefit of integrating quality control flags and quantitative measurements and highlight the advantage of joint analysis of population-based samples and phenotypic extremes in human genetics studies., Competing Interests: Declaration of interests Massachusetts Institute of Technology filed patent applications based on the findings and regarding the inclusive polygenic score approach used in the study. Y.T. and M.K. are designated as inventors of the applications., (Copyright © 2024 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Background: Russian wildrye (RWR, Psathyrostachys juncea ) is an outcrossing perennial grass that plays a crucial role in foragaing and rangeland restoration due to its tiller producing capabilities, nevertheless, a genetic map has yet to be constructed due to a shortage of efficient and reliable molecular markers. This also limits the identification, localization, and cloning of economically important traits related to tiller density during breeding., Methods: Therefore, this study aimed to create a F 1 mapping population with 147 individual lines and their two parents, which were selected based on varying tiller densities. We then used this mapping population to conduct specific-locus amplified fragment sequencing (SLAF-seq) to generate SLAF markers and discover single nucleotide polymorphisms (SNPs)., Results: Initially, we generated a total of 1,438.38 million pair-end reads with an average sequencing depth of 84.92 in the maternal line, 79.34 in the parental line, and 27.05 in each F 1 individual line, respectively. Following the filtering of low-depth SLAF tags, a total of 558,344 high-quality SLAFs were identified. A total of 1,519,903 SNP markers were obtained, and 62,424 polymorphic SNPs were discovered. From these, 4,644 polymorphic SNPs were selected and used for the construction of a genetic map encompassing seven linkage groups. The genetic map spanned 1,416.60 cM with an average distance of 0.31 cM between adjacent markers. Comparative analysis between the seven linkage groups of RWR SLAF tag and the whole-genome sequences in barley ( Hordeum vulgare L.) revealed homology values ranging from 17.5% to 34.6%, and the collinearity between the RWR linkage groups and the barley homology groups ranged from 0.6787 to 0.9234, with an average value of 0.8158. Additionally, 143 significant quantitative trait locus (QTLs) with Logarithm of Odds (LOD) value greater than 2.5 for five tiller related traits were detected using three consecutive years of phenotypic trait data from the F1 population, further verifying the map's reliability., Competing Interests: The authors declare that they have no competing interests., (© 2024 Ma et al.)

Amyotrophic Lateral Sclerosis (ALS) is a devastating, immensely complex neurodegenerative disease by lack of effective treatments. We developed a network medicine methodology via integrating human brain multi-omics data to prioritize drug targets and repurposable treatments for ALS. We leveraged non-coding ALS loci effects from genome-wide associated studies (GWAS) on human brain expression quantitative trait loci (QTL) (eQTL), protein QTL (pQTL), splicing QTL (sQTL), methylation QTL (meQTL), and histone acetylation QTL (haQTL). Using a network-based deep learning framework, we identified 105 putative ALS-associated genes enriched in known ALS pathobiological pathways. Applying network proximity analysis of predicted ALS-associated genes and drug-target networks under the human protein-protein interactome (PPI) model, we identified potential repurposable drugs (i.e., Diazoxide and Gefitinib) for ALS. Subsequent validation established preclinical evidence for top-prioritized drugs. In summary, we presented a network-based multi-omics framework to identify drug targets and repurposable treatments for ALS and other neurodegenerative disease if broadly applied., (© 2024. The Author(s).)

Flesh firmness is a critical breeding trait that determines consumer selection, shelf life, and transportation. The genetic basis controlling firmness in apple (Malus × domestica Borkh.) remains to be fully elucidated. We aimed to decipher genetic variance for firmness at harvest and develop potential molecular markers for marker-assisted breeding. Maturity firmness for 439 F1 hybrids from a cross of "Cripps Pink" and "Fuji" was determined in 2016 and 2017. The phenotype segregated extensively, with a Gaussian distribution. In a combined bulked segregant analysis (BSA) and RNA-sequencing analysis, 84 differentially expressed genes were screened from the 10 quantitative trait loci regions. Interestingly, next-generation re-sequencing analysis revealed a Harbinger-like transposon element insertion upstream of the candidate gene PECTATE LYASE5 (MdPL5); the genotype was associated with flesh firmness at harvest. The presence of this transposon repressed MdPL5 expression and was closely linked to the extra-hard phenotype. MdPL5 was demonstrated to promote softening in apples and tomatoes. Subsequently, using the MdPL5 promoter as bait, MdNAC1-L was identified as a transcription activator that positively regulates ripening and softening in the developing fruit. We also demonstrated that MdNAC1-L could induce the up-regulation of MdPL5, MdPG1, and the ethylene-related genes MdACS1 and MdACO1. Our findings provide insight into TE-related genetic variation and the PL-mediated regulatory network for the firmness of apple fruit., Competing Interests: Conflict of interest statement. None declared., (© The Author(s) 2024. Published by Oxford University Press on behalf of American Society of Plant Biologists. All rights reserved. For commercial re-use, please contact reprints@oup.com for reprints and translation rights for reprints. All other permissions can be obtained through our RightsLink service via the Permissions link on the article page on our site—for further information please contact journals.permissions@oup.com.)

The worldwide cultivated Cucurbita pepo L. is one of the most diverse species in the plant kingdom. In this study, chilling tolerance over a wide range of cultivars was characterized to discover the allelic variants to improving the postharvest quality of the immature fruit during cold storage. For this purpose, fruits from 126 accessions of worldwide origin have been evaluated for weight loss and chilling injury after 3, 7 and 14 days of cold storage, classifying them into tolerant, partially tolerant, and sensitive accessions. To verify this classification, antioxidant capacity and lipid peroxidation (MDA) of contrasting accessions (tolerant vs. sensitive) were assessed. The antioxidant capacity significantly decreased during cold storage in the sensitive accessions, while it was maintained in tolerant accessions. Additionally, the sensitive accessions presented a higher accumulation of MDA during this period. Finally, a GWAS analysis using GBS data available in CuGenDBv2, combined with weight loss percentage data, led to the identification of a candidate QTL located on chromosome 17 that regulates postharvest cold tolerance in zucchini. The region contains four SNPs whose alternative alleles were significantly associated with lower weight loss percentage and chilling injury indices during cold storage. Two SNPs are located in the 3' UTR region of the gene CpERS1, a gene involved in ethylene perception. The other two SNPs generate missense mutations in the coding region of a Pectin methyl esterase inhibitor gene (CpPMI). The role of this QTL and these variants in chilling tolerance is discussed., (© 2024 The Author(s). Physiologia Plantarum published by John Wiley & Sons Ltd on behalf of Scandinavian Plant Physiology Society.)

Functional genes within genomic loci associated with systemic lupus erythematosus (SLE), as identified by genome-wide association studies, exhibit cell-specific characteristics. This study delves into the impact of genetic variants within SLE loci on gene expression in different types of immune cells, unraveling the complex interplay between genetics and immunopathogenesis. Through the integration of genetic association and single-cell transcriptomic sequencing data, we identified potential cell-specific susceptibility genes for SLE across diverse immune cell subsets. The single-cell eQTL analysis revealed 30,409 associations involving 3583 SLE-associated SNPs. These SNPs exhibited associations with expression levels of 147 genes across 14 distinct cell types. The single-cell summary data-based Mendelian randomization (SMR) analysis identified 119 significant associations between the expression levels of 44 genes and SLE. Notably, myeloid cells exhibited associations solely within the MHC region, while T, B, and natural killer cells showed associations with both MHC and non-MHC genes in relation to SLE. Analysis of single-cell transcriptomic data from 33 children SLE cases and 11 match controls (227,303 cells), as well as 7 adult SLE cases and 5 match controls (78,414 cells) highlights differential expression of key genes. Notably, genetic variants within HLA-DRB1, HLA-DRB5, HLA-DQA1, HLA-DQB1, IRF7, IRF5, BLK and HLA-DPA1 play a pivotal role in mediating immune dysregulation in specific immune cell types. Our study contributes to a comprehensive understanding of the intricate relationships between genetics, gene expression and SLE susceptibility. The findings shed light on the cell-specific impacts of genetic variants within SLE-associated genomic loci., Competing Interests: Declaration of Competing Interest The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Lotus (Nelumbo nucifera G.) is a perennial aquatic horticultural plant with diverse architectures. Distinct plant architecture (PA) has certain attractive and practical qualities, but its genetic morphogenesis in lotus remains elusive. In this study, we employ genome-wide association analysis (GWAS) for the seven traits of petiole length (PLL), leaf length (LL), leaf width (LW), peduncle length (PLF), flower diameter (FD), petal length (PeL), and petal width (PeW) in 301 lotus accessions. A total of 90 loci are identified to associate with these traits across 4 years of trials. Meanwhile, we perform RNA sequencing (RNA-seq) to analyze the differential expression of the gene (DEG) transcripts between large and small PA (LPA and SPA) of lotus stems (peduncles and petioles). As a result, eight key candidate genes are identified that are all primarily involved in plant cell wall remodeling significantly associated with PA traits by integrating the results of DEGs and GWAS. To verify this result, we compare the cell wall compositions and structures of LPA versus SPA in representative lotus germplasms. Intriguingly, compared with the SPA lotus, the LPA varieties have higher content of cellulose and hemicellulose, but less filling substrates of pectin and lignin. Additionally, we verified longer cellulose chains and higher cellulose crystallinity with less interference in LPA varieties. Taken together, our study illustrates how plant cell wall remodeling affects PA in lotus, shedding light on the genetic architecture of this significant ornamental trait and offering a priceless genetic resource for future genomic-enabled breeding., (© 2024 Society for Experimental Biology and John Wiley & Sons Ltd.)

Plant trichomes are an excellent model for studying cell differentiation and development, providing crucial defenses against biotic and abiotic stresses. There is a well-established inverse relationship between trichome density and aphid prevalence, indicating that higher trichome density leads to reduced aphid infestations. Here we present the cloning and characterization of a dominant quantitative trait locus, HIC (hirsute cotton), which significantly enhances cotton trichome density. This enhancement leads to markedly improved resistance against cotton aphids. The HIC encodes an HD-ZIP IV transcriptional activator, crucial for trichome initiation. Overexpression of HIC leads to a substantial increase in trichome density, while knockdown of HIC results in a marked decrease in density, confirming its role in trichome regulation. We identified a variant in the HIC promoter (-810 bp A to C) that increases transcription of HIC and trichome density in hirsute cotton compared with Gossypium hirsutum cultivars with fewer or no trichomes. Interestingly, although the -810 variant in the HIC promoter is the same in G. barbadense and hirsute cotton, the presence of a copia-like retrotransposon insertion in the coding region of HIC in G. barbadense causes premature transcription termination. Further analysis revealed that HIC positively regulates trichome density by directly targeting the EXPANSIN A2 gene, which is involved in cell wall development. Taken together, our results underscore the pivotal function of HIC as a primary regulator during the initial phases of trichome formation, and its prospective utility in enhancing aphid resistance in superior cotton cultivars via selective breeding., (© 2024 Society for Experimental Biology and John Wiley & Sons Ltd.)

Tillers are specialized lateral shoots arising from axillary buds at basal nodes, and are also an important agronomic trait that determines the aboveground biomass and grain yield of various gramineous crops. So far, few genes have been reported to control tiller formation and most have been in the annual crop rice (Oryza sativa). Orchardgrass (Dactylis glomerata) is an important perennial forage crop with great economic and ecological value, but its genes regulating tillering have remained largely unknown. In the present study, we used a natural population of 264 global orchardgrass germplasms to determine genes associated with quantitative variation in tiller number through genome-wide association study analysis. A total of 19 putative loci and 55 genes associated with tiller number were thus identified. Additionally, 26 putative differentially expressed genes with tiller number, including DgCYC-C1, were identified by RNA-seq and genome-wide association study analysis. DgCYC-C1 which is involved in cell division, was overexpressed, revealing that DgCYC-C1 positively regulates tiller number. These results provide some new candidate genes or loci for the improvement of tiller number in crops, which might advance new sustainable strategies to meet global crop production challenges., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024. Published by Elsevier Masson SAS.)

There is an urgent need to improve wheat for upcoming challenges, including biotic and abiotic stresses. Sustainable wheat improvement requires the introduction of new genes and alleles in high-yielding wheat cultivars. Using new approaches, tools, and technologies to identify and introduce new genes in wheat cultivars is critical. High-quality genomes, transcriptomes, and pangenomes provide essential resources and tools to examine wheat closely to identify and manipulate new and targeted genes and alleles. Wheat genomics has improved excellently in the past 5 years, generating multiple genomes, pangenomes, and transcriptomes. Leveraging these resources allows us to accelerate our crop improvement pipelines. This review summarizes the progress made in wheat genomics and trait discovery in the past 5 years., Competing Interests: Declaration of interests The authors declare no conflicts of interest., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Transcriptome-wide association studies (TWAS) have provided valuable insight in identifying genes that may impact cigarette smoking. Most of previous studies, however, mainly focused on European ancestry. Limited TWAS studies have been conducted across multiple ancestries to explore genes that may impact smoking behaviors. In this study, we used cis-eQTL data of cerebral cortex from multiple ancestries in MetaBrain, including European, East Asian, and African samples, as reference panels to perform multi-ancestry TWAS analyses on ancestry-matched GWASs of four smoking behaviors including smoking initiation, smoking cessation, age of smoking initiation, and number of cigarettes per day in GWAS & Sequencing Consortium of Alcohol and Nicotine use (GSCAN). Multiple-ancestry fine-mapping approach was conducted to identify credible gene sets associated with these four traits. Enrichment and module network analyses were further performed to explore the potential roles of these identified gene sets. A total of 719 unique genes were identified to be associated with at least one of the four smoking traits across ancestries. Among those, 249 genes were further prioritized as putative causal genes in multiple ancestry-based fine-mapping approach. Several well-known smoking-related genes, including PSMA4, IREB2, and CHRNA3, showed high confidence across ancestries. Some novel genes, e.g., TSPAN3 and ANK2, were also identified in the credible sets. The enrichment analysis identified a series of critical pathways related to smoking such as synaptic transmission and glutamate receptor activity. Leveraging the power of the latest multi-ancestry GWAS and eQTL data sources, this study revealed hundreds of genes and relevant biological processes related to smoking behaviors. These findings provide new insights for future functional studies on smoking behaviors., (© 2024. The Author(s), under exclusive licence to Springer Nature Limited.)

Psychiatric disorders are highly heritable yet polygenic, potentially involving hundreds of risk genes. Genome-wide association studies have identified hundreds of genomic susceptibility loci with susceptibility to psychiatric disorders; however, the contribution of these loci to the underlying psychopathology and etiology remains elusive. Here we generated deep human brain proteomics data by quantifying 11,608 proteins across 268 subjects using 11-plex tandem mass tag coupled with two-dimensional liquid chromatography-tandem mass spectrometry. Our analysis revealed 788 cis-acting protein quantitative trait loci associated with the expression of 883 proteins at a genome-wide false discovery rate <5%. In contrast to expression at the transcript level and complex diseases that are found to be mainly influenced by noncoding variants, we found protein expression level tends to be regulated by non-synonymous variants. We also provided evidence of 76 shared regulatory signals between gene expression and protein abundance. Mediation analysis revealed that for most (88%) of the colocalized genes, the expression levels of their corresponding proteins are regulated by cis-pQTLs via gene transcription. Using summary data-based Mendelian randomization analysis, we identified 4 proteins and 19 genes that are causally associated with schizophrenia. We further integrated multiple omics data with network analysis to prioritize candidate genes for schizophrenia risk loci. Collectively, our findings underscore the potential of proteome-wide linkage analysis in gaining mechanistic insights into the pathogenesis of psychiatric disorders., (© 2024. The Author(s).)

Cold stress detrimentally influences fruit development, leading to a substantial yield reduction in many fruit-bearing vegetables. Cucumber, a vegetable of subtropical origin, is especially sensitive to cold. Cold-inducible parthenocarpy (CIP) promises fruit yield under cold conditions. Previously, we identified a CIP line EC5 in cucumber, which showed strong parthenocarpy and sustained fruit growth under cold conditions (16°C day/10°C night). However, the candidate gene and genetic mechanism underlying CIP in cucumber remain unknown. In this study, both BSA-seq and conventional QTL mapping strategies were employed on F 2 populations to delve into the genetic control of CIP. A single QTL, CIP5.1, was consistently mapped across two winter seasons in 2021 and 2022. Fine mapping delimited the CIP locus into a 38.3 kb region on chromosome 5, harboring 8 candidate genes. Among these candidates, CsAGL11 (CsaV3&#95;5G040370) was identified, exhibiting multiple deletions/insertions in the promoter and 5'UTR region. The CsAGL11 gene encodes a MADS-box transcription factor protein, which is homologous to the genes previously recognized as negative regulators in ovule and fruit development of Arabidopsis and tomato. Correspondingly, cold treatment resulted in decreased expression of CsAGL11 during the early developmental stage of the fruit in EC5. A promoter activity assay confirmed promoter polymorphisms leading to weak transcriptional activation of CsAGL11 under cold conditions. This study deepens our understanding of the genetic characteristics of CIP and elucidates the potential role of the CsAGL11 gene in developing cucumber cultivars with enhanced fruiting under cold conditions., Competing Interests: Declaration of Competing Interest All the authors declare no conflict of interest., (Copyright © 2024. Published by Elsevier B.V.)

Tocochromanols, collectively known as Vitamin E, serve as natural lipid-soluble antioxidants that are exclusively obtained through dietary intake in humans. Synthesized by all plants, tocochromanols play an important role in protecting polyunsaturated fatty acids in plant seeds from lipid peroxidation. While the genes involved in tocochromanol biosynthesis have been fully elucidated in Arabidopsis thaliana, Oryza sativa and Zea mays, the genetic basis of tocochromanol accumulation in sweet corn remains poorly understood. This gap is a consequence of limited natural genetic diversity and harvest at immature growth stages. In this study, we conducted comprehensive genome-wide association studies (GWAS) on a sweet corn panel of 295 individuals with a high-density molecular marker set. In total, thirteen quantitative trait loci (QTLs) for individual and derived tocochromanol traits were identified. Our analysis identified novel roles for three genes, ZmCS2, Zmshki1 and ZmB4FMV1, in the regulation of α-tocopherol accumulation in sweet corn kernels. We genetically validated the role of Zmshki1 through the generation of a knock-out line using CRISPR-Cas9 technology. Further gene-based GWAS revealed the function of the canonical tyrosine metabolic enzymes ZmCS2 and Zmhppd1 in the regulation of total tocochromanol content. This comprehensive assessment of the genetic basis for variation in vitamin E content establishes a solid foundation for enhancing vitamin E content not only in sweet corn, but also in other cereal crops., Competing Interests: Declaration of Competing Interest The authors declare that they have no conflict of interest., (Copyright © 2024 Elsevier B.V. All rights reserved.)

Blood composition is indicative of health-related traits such as immunity and metabolism. The use of molecular genetics to investigate alterations in these attributes in laying ducks is a novel approach. Our objective was to employ genome - wide association studies (GWAS) and haplotype - sharing analysis to identify genomic regions and potential genes associated with 11 blood components in Shaoxing ducks. Our findings revealed 35 SNPs and 1 SNP associated with low-density lipoprotein cholesterol (LDL) and globulin (GLB), respectively. We identified 36 putative candidate genes for the LDL trait in close proximity to major QTLs and key loci. Based on their biochemical and physiological properties, TRA2A , NPY , ARHGEF26 , DHX36 , and AADAC are the strongest putative candidate genes. Through linkage disequilibrium analysis and haplotype sharing analysis, we identified three haplotypes and one haplotype, respectively, that were significantly linked with LDL and GLB. These haplotypes could be selected as potential candidate haplotypes for molecular breeding of Shaoxing ducks. Additionally, we utilized a bootstrap test to verify the reliability of GWAS with small experimental samples. The test can be accessed at https://github.com/xuwenwu24/Bootstrap-test.

The elucidation of genetic architecture and molecular regulatory networks underlying complex traits remains a significant challenge in life science, largely due to the substantial background effects that arise from epistasis and gene-environment interactions. The chromosome segment substitution line (CSSL) is an ideal material for genetic and molecular dissection of complex traits due to its near-isogenic properties; yet a comprehensive analysis, from the basic identification of substitution segments to advanced regulatory network, is still insufficient. Here, we developed two cotton CSSL populations on the Gossypium hirsutum background, representing wide adaptation and high lint yield, with introgression from G. barbadense, representing superior fibre quality. We sequenced 99 CSSLs that demonstrated significant differences from G. hirsutum in fibre, and characterized 836 dynamic fibre transcriptomes in three crucial developmental stages. We developed a workflow for precise resolution of chromosomal substitution segments; the genome sequencing revealed substitutions collectively representing 87.25% of the G. barbadense genome. Together, the genomic and transcriptomic survey identified 18 novel fibre-quality-related quantitative trait loci with high genetic contributions and the comprehensive landscape of fibre development regulation. Furthermore, analysis determined unique cis-expression patterns in CSSLs to be the driving force for fibre quality alteration; building upon this, the co-expression regulatory network revealed biological relationships among the noted pathways and accurately described the molecular interactions of GhHOX3, GhRDL1 and GhEXPA1 during fibre elongation, along with reliable predictions for their interactions with GhTBA8A5. Our study will enhance more strategic employment of CSSL in crop molecular biology and breeding programmes., (© 2024 The Author(s). Plant Biotechnology Journal published by Society for Experimental Biology and The Association of Applied Biologists and John Wiley & Sons Ltd.)

To understand whether domestication had an impact on susceptibility and responsiveness to arbuscular mycorrhizal fungi (AMF) in tomato (Solanum lycopersicum), we investigated two tomato cultivars ("M82" and "Moneymaker") and a panel of wild relatives including S. neorickii, S. habrochaites and S. pennellii encompassing the whole Lycopersicon clade. Most genotypes revealed good AM colonisation levels when inoculated with the AMF Funneliformis mosseae. By contrast, both S. pennellii accessions analysed showed a very low colonisation, but with normal arbuscule morphology, and a negative response in terms of root and shoot biomass. This behaviour was independent of fungal identity and environmental conditions. Genomic and transcriptomic analyses revealed in S. pennellii the lack of genes identified within QTLs for AM colonisation, a limited transcriptional reprogramming upon mycorrhization and a differential regulation of strigolactones and AM-related genes compared to tomato. Donor plants experiments indicated that the AMF could represent a cost for S. pennellii: F. mosseae could extensively colonise the root only when it was part of a mycorrhizal network, but a higher mycorrhization led to a higher inhibition of plant growth. These results suggest that genetics and functional traits of S. pennellii are responsible for the limited extent of AMF colonisation., (© 2024 The Author(s). Plant, Cell & Environment published by John Wiley & Sons Ltd.)

Background: Rice ECQ (eating and cooking quality) is an important determinant of rice consumption and market expansion. Therefore, improvement of ECQ is one of the primary goals in rice breeding. However, ECQ-related quantitative trait loci (QTL) have not yet been fully revealed. The present study aimed to identify a major effect QTL for rice taste, an important component of ECQ via genotyping-by-sequencing, to reveal the associated molecular mechanisms, and to predict key candidate genes., Results: A population of F 9 recombinant inbred lines resulting from a cross between R668 (national standard of high-quality third class) and R838 (common edible rice) was used to construct a high-density genetic map (2,295.062 cM). The map comprises 639,504 markers distributed on 12 linkage elements with an average genetic distance of 0.004 cM. We detected a major taste-related QTL, qECQ8, which explained 41.4% of phenotypic variance and had LOD values of 4.42-7.73. Using a five-generation NIL population from the backcross of "Ganxiangzhan No. 1" carrying qECQ8 with the recurrent parent R838 (without qECQ8), we narrowed qECQ8 to a 187.5 kb interval between markers M33 and M37 on Chr8. Comparative transcriptomic analysis revealed that photosynthesis, glyoxylate and dicarboxylate metabolism, carbon fixation in photosynthetic organisms, and alpha-linolenic acid metabolism were induced in developing seeds of lines containing qECQ8. Furthermore, we identified two candidate genes in the qECQ8 region, including LOC&#95;Os08g30550 (zinc knuckle family protein), a major candidate for genetic-assisted breeding of high-quality rice., Conclusion: Our findings provide important genetic resources for targeted improvement of rice taste quality and may facilitate the genetic breeding of rice ECQ., (© 2024. The Author(s).)

Lifespan is influenced by complex interactions between genetic and environmental factors. Studying those factors in model organisms of a single genetic background limits their translational value for humans. Here, we mapped lifespan determinants in 85 C. elegans recombinant inbred advanced intercross lines (RIAILs). We assessed molecular profiles-transcriptome, proteome, and lipidome-and life-history traits, including lifespan, development, growth dynamics, and reproduction. RIAILs exhibited large variations in lifespan, which correlated positively with developmental time. We validated three longevity modulators, including rict-1, gfm-1, and mltn-1, among the top candidates obtained from multiomics data integration and quantitative trait locus (QTL) mapping. We translated their relevance to humans using UK Biobank data and showed that variants in GFM1 are associated with an elevated risk of age-related heart failure. We organized our dataset as a resource that allows interactive explorations for new longevity targets., Competing Interests: Declaration of interests J.J.C. is a consultant for Thermo Scientific, Seer, and 908 Devices. E.K., L.M., and M.C. are employees of and J.A. is a shareholder of Nagi Bioscience S.A., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

Background: Frozen shoulder (FS) is characterized by the thickening and fibrosis of the joint capsule, leading to joint contracture and a reduction in joint volume. The precise etiology responsible for these pathological changes remains elusive. Therefore, the primary aim of this study was to explore the potential involvement of pathogenic genes in FS and analyze their underlying roles in the disease progression. Methods: Differential expression analysis and weighted gene co-expression network analysis (WGCNA) were employed to investigate co-expressed genes potentially associated with FS. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses were conducted to elucidate the potential roles of these co-expressed genes. Subsequently, Mendelian randomization (MR) analysis was performed using expression quantitative trait loci datasets for the co-expressed genes, combined with summary statistics from the genome-wide association study of FS, aiming to identify key genes causally associated with FS. The identified key genes were further validated through reverse transcription-quantitative PCR (RT-qPCR). Additionally, a nomogram model and receiver operating characteristic (ROC) curves were established to assess the diagnostic value of the hub genes. Furthermore, the infiltration of immune cells was evaluated using the CIBERSORT algorithm and the relationship between key genes and immune-infiltrating cells was analyzed. Results: 295 overlapping co-expressed genes were identified by intersecting the differentially expressed genes with the hub genes obtained from associated modules identified through WGCNA. Utilizing MR analysis, four key genes, namely ADAMTS1, NR4A2, PARD6G and SMKR1, were found to exhibit positive causal relationships with FS, which were subsequently validated through RT-qPCR analysis. Moreover, the diagnostic value of these four key genes was demonstrated through the development of a nomogram model and the construction of ROC curves. Notably, a causal relationship between ADAMTS1 and immune cell infiltration in FS was observed. Conclusion: Our study suggested genetic predisposition to higher expression levels of ADAMTS1, NR4A2, PARD6G and SMKR1, was associated with an increased risk of FS. Further investigations elucidating the functional roles of these genes will enhance our understanding of the pathogenesis of FS and may facilitate the development of targeted treatment strategies., Competing Interests: Competing Interests: The authors have declared that no competing interest exists., (© The author(s).)

Key Message: Cellulose synthase-like OsCSLD4 plays a pivotal role in regulating diverse agronomic traits, enhancing resistance against bacterial leaf blight, and modulating metabolite indices based on the multi-omics analysis in rice. To delve deeper into this complex network between agronomic traits and metabolites in rice, we have compiled a dataset encompassing genome, phenome, and metabolome, including 524 diverse accessions, 11 agronomic traits, and 841 metabolites, enabling us to pinpoint eight hotspots through GWAS. We later discovered four distinct metabolite categories, encompassing 15 metabolites that are concurrently present on the QTL qC12.1, associated with leaf angle of flag and spikelet length, and finally focused the cellulose synthase-like OsCSLD4, which was pinpointed through a rigorous process encompassing sequence variation, haplotype, ATAC, and differential expression across diverse tissues. Compared to the wild type, csld4 exhibited significant reductions in the plant height, flag leaf length, leaf width, spikelet length, 1000-grain weight, grain width, grain thickness, fertility, yield per plant, and bacterial blight resistance. However, there were significant increase in tiller numbers, degree of leaf rolling, flowering period, growth period, grain length, and empty kernel rate. Furthermore, the content of four polyphenol metabolites, excluding metabolite N-feruloyltyramine (mr1268), notably rose, whereas the levels of the other three polyphenol metabolites, smiglaside C (mr1498), 4-coumaric acid (mr1622), and smiglaside A (mr1925) decreased significantly in mutant csld4. The content of amino acid L-tyramine (mr1446) exhibited a notable increase, whereas the alkaloid trigonelline (mr1188) displayed a substantial decrease among the mutants. This study offered a comprehensive multi-omics perspective to analyze the genetic mechanism of OsCSLD4, and breeders can potentially enhance rice's yield, bacterial leaf blight resistance, and metabolite content, leading to more sustainable and profitable rice production., (© 2024. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Improving the powdery mildew resistance of bitter gourd is highly important for achieving high yield and high quality. To better understand the genetic basis of powdery mildew resistance in bitter gourd, this study analyzed 300 lines of recombinant inbred lines (RILs) formed by hybridizing the powdery mildew-resistant material MC18 and the powdery mildew-susceptible material MC402. A high-density genetic map of 1222.04 cM was constructed via incorporating 1,996,505 SNPs generated by resequencing data from 180 lines, and quantitative trait locus (QTL) positioning was performed using phenotypic data at different inoculation stages. A total of seven QTLs related to powdery mildew resistance were identified on four chromosomes, among which qPm-3-1 was detected multiple times and at multiple stages after inoculation. By selecting 18 KASP markers that were evenly distributed throughout the region, 250 lines and parents were genotyped, and the interval was narrowed to 207.22 kb, which explained 13.91% of the phenotypic variation. Through RNA-seq analysis of the parents, 11,868 differentially expressed genes (DEGs) were screened. By combining genetic analysis, gene coexpression, and sequence comparison analysis of extreme materials, two candidate genes controlling powdery mildew resistance in bitter gourd were identified ( evm.TU.chr3.2934 (C3H) and evm.TU.chr3.2946 (F-box-LRR)). These results represent a step forward in understanding the genetic regulatory network of powdery mildew resistance in bitter gourd and lay a molecular foundation for the genetic improvement in powdery mildew resistance.

A genome-wide association study (GWAS) identifies trait-associated loci, but identifying the causal genes can be a bottleneck, due in part to slow decay of linkage disequilibrium (LD). A transcriptome-wide association study (TWAS) addresses this issue by identifying gene expression-phenotype associations or integrating gene expression quantitative trait loci with GWAS results. Here, we used self-pollinated soybean (Glycine max [L.] Merr.) as a model to evaluate the application of TWAS to the genetic dissection of traits in plant species with slow LD decay. We generated RNA sequencing data for a soybean diversity panel and identified the genetic expression regulation of 29 286 soybean genes. Different TWAS solutions were less affected by LD and were robust to the source of expression, identifing known genes related to traits from different tissues and developmental stages. The novel pod-color gene L2 was identified via TWAS and functionally validated by genome editing. By introducing a new exon proportion feature, we significantly improved the detection of expression variations that resulted from structural variations and alternative splicing. As a result, the genes identified through our TWAS approach exhibited a diverse range of causal variations, including SNPs, insertions or deletions, gene fusion, copy number variations, and alternative splicing. Using this approach, we identified genes associated with flowering time, including both previously known genes and novel genes that had not previously been linked to this trait, providing insights complementary to those from GWAS. In summary, this study supports the application of TWAS for candidate gene identification in species with low rates of LD decay., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Flowering time is an important agronomic trait for canola breeders, as it provides growers with options for minimizing exposure to heat stress during flowering and to more effectively utilize soil moisture. Plants have evolved various systems to control seasonal rhythms in reproductive phenology including an internal circadian clock that responds to environmental signals. In this study, we used canola cultivar 'Westar' as a recurrent parent and canola cultivar 'Surpass 400' as the donor parent to generate a chromosome segment substitution line (CSSL) and to map a flowering time locus on chromosome A10 using molecular marker-assisted selection. This CSSL contains an introgressed 4.6 mega-bases (Mb) segment (between 13 and 17.6 Mb) of Surpass 400, which substantially delayed flowering compared with Westar. To map flowering time gene(s) within this locus, eight introgression lines (ILs) were developed carrying a series of different lengths of introgressed chromosome A10 segments using five co-dominant polymorphic markers located at 13.5, 14.0, 14.5, 15.0, 15.5, and 16.0 Mb. Eight ILs were crossed with Westar reciprocally and flowering time of resultant 16 F 1 hybrids and parents were evaluated in a greenhouse (2021 and 2022). Four ILs (IL005, IL017, IL035, and IL013) showed delayed flowering compared to Westar (P < 0.0001), and their reciprocal crosses displayed a phenotype intermediate in flowering time of both homozygote parents. These results indicated that flowering time is partial or incomplete dominance, and the flowering time locus mapped within a 1 Mb region between two co-dominant polymorphic markers at 14.5-15.5 Mb on chromosome A10. The flowering time locus was delineated to be between 14.60 and 15.5 Mb based on genotypic data at the crossover site, and candidate genes within this region are associated with flowering time in canola and/or Arabidopsis. The co-dominant markers identified on chromosome A10 should be useful for marker assisted selection in breeding programs but will need to be validated to other breeding populations or germplasm accessions of canola., (© 2024. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Background: Vegetable soybean is rich in nutrients and has a unique flavor. It is highly preferred by people because of its pharmacological activities, including those that regulate the intestines and lower blood pressure. The pod color of vegetable soybeans is an important quality that indicates their freshness and has a significant impact on their commercialization., Results: In this study, pod color was evaluated in 301 vegetable soybean accessions collected from various regions. Genome-wide association analysis was carried out using the Mixed linear model (MLM), a total of 18 quantitative trait loci including 117 SNPs were detected. Two significant QTLs located on chromosomes 6 (qGPCL4 /qGPCa1/qGPCb2) and 18 (qGPCL10/qGPCb3) were consistently detected across different variables. Based on gene functional annotation, 30 candidate genes were identified in these two candidate intervals. Combined with transcriptome analysis, Glyma.18g241700 has been identified as a candidate gene for regulating pod color in vegetable soybeans. Glyma.18g241700 encodes a chlorophyll photosystem I subunit XI. which localizes to the chloroplast named GmPsaL, qRT-PCR analysis showed that GmPsaL was specifically highly expressed in developing pods. Furthermore, overexpression of GmPsaL in transgenetic Arabidopsis plants produced dark green pods., Conclusions: These findings may be useful for clarifying the genetic basis of the pod color of vegetable soybeans. The identified candidate genes may be useful for the genetic improvement of the appearance quality of vegetable soybeans., (© 2024. The Author(s).)

Background: Mitochondrial dysfunction (MD) is increasingly recognized as a key pathophysiological contributor in Alzheimer disease (AD). As differential MD genes expression may serve as either a causative factor or a consequence in AD, and expression of these genes could be influenced by epigenetic modifications or interact with inflammatory cytokines, hence, the precise role of MD in AD remains uncertain., Methods: Meta-analysis of brain transcriptome datasets was conducted to pinpoint differentially expressed genes (DEGs) associated with MD in AD. We utilized three-step SMR to analyze the AD genome-wide association study summaries with expression quantitative trait loci (eQTLs) and DNA methylation QTLs from the blood and brain tissues, respectively. Through SMR and colocalization analysis, we further explored the interactions between brain eQTLs and inflammatory cytokines., Results: Five datasets were meta-analyzed to prioritize 825 DEGs in AD from 1339 MD-related genes. Among these, seven genes from blood samples such as NDUFS8 and SPG7 and thirty-two genes from brain tissue including CLU and MAPT were identified as candidate AD-causal MD genes and regulated by methylation level. Furthermore, we revealed 13 MD gene expression-inflammatory pathway pairs involving LDLR, ACE and PTPMT1 along with interleukin-17C, interleukin-18 and hepatocyte growth factor., Conclusions: This study highlighted that the AD-causal MD genes could be regulated by epigenetic changes and interact with inflammatory cytokines, providing evidence for AD prevention and intervention., (© 2024. The Author(s).)

Alternative splicing (AS) plays crucial roles in regulating various biological processes in plants. However, the genetic mechanisms underlying AS and its role in controlling important agronomic traits in rice (Oryza sativa) remain poorly understood. In this study, we explored AS in rice leaves and panicles using the rice minicore collection. Our analysis revealed a high level of transcript isoform diversity, with approximately one-fifth of the potential isoforms acting as major transcripts in both tissues. Regarding the genetic mechanism of AS, we found that the splicing of 833 genes in the leaf and 1,230 genes in the panicle was affected by cis-genetic variation. Twenty-one percent of these AS events could only be explained by large structural variations. Approximately 77.5% of genes with significant splicing quantitative trait loci (sGenes) exhibited tissue-specific regulation, and AS can cause 26.9% (leaf) and 23.6% (panicle) of sGenes to have altered, lost, or gained functional domains. Additionally, through splicing-phenotype association analysis, we identified phosphate-starvation-induced RING-type E3 ligase (OsPIE1; LOC&#95;Os01g72480), whose splicing ratio was significantly associated with plant height. In summary, this study provides an understanding of AS in rice and its contribution to the regulation of important agronomic traits., Competing Interests: Conflict of interest statement. None declared., (© The Author(s) 2024. Published by Oxford University Press on behalf of American Society of Plant Biologists. All rights reserved. For commercial re-use, please contact reprints@oup.com for reprints and translation rights for reprints. All other permissions can be obtained through our RightsLink service via the Permissions link on the article page on our site—for further information please contact journals.permissions@oup.com.)

Objective: N6-methyladenosine (m 6 A) is the most prevalent modification of RNA in eukaryotes which is associated with many cellular processes and diseases. Here, our objective is to explore whether genetic variants in m 6 A modification genes are associated with the risk of non-syndrome orofacial clefts (NSOCs)., Methods: The transmission disequilibrium test (TDT) was performed to calculate the association between single nucleotide polymorphisms (SNPs) in m 6 A modification genes and NSOCs risk in 944 case-parent trios. The function of SNP was predicted by HaploReg, RegulomeDB and histone enrichment data. The expression quantitative trait locus (eQTL) analysis was examined using Genotype-Tissue Expression (GTEx) and eQTLGen. The role of gene in the development of NSOCs was assessed with correlation and enrichment analysis based on gene expression data in mice craniofacial tissue and zebrafish embryo., Results: We identified that rs8078195 (A > C) in METTL16 was suggestively associated with the increased risk of NSOCs (OR = 1.32, p = 1.80E - 03). The region surrounding rs8078195 was subjected to deoxyribonuclease hypersensitivity and enriched with multiple histone modifications. In addition, it had a significant eQTL effect with METTL16 in skin tissue and human peripheral blood, which played an important role in NSOCs development. Bioinformatic analysis indicated that METTL16 contributed to the development of NSOCs probably by regulating cell cycle process., Conclusions: Rs8078195 in METTL16 was associated with the occurrence of NSOCs., (© 2024 Wiley Periodicals LLC.)

Identifying genomic regions under selection is essential for understanding the genetic mechanisms driving species evolution and adaptation. Traditional methods often fall short in detecting complex, spatially varying selection signals. Recent advances in deep learning, however, present promising new approaches for uncovering subtle selection signals that traditional methods might miss. In this study, we utilized the deep learning framework DeepGenomeScan to detect spatially varying selection signatures across 15 bovine populations worldwide. Our analysis uncovered novel insights into selective sweep hotspots within the bovine genome, revealing key genes associated with physiological and adaptive traits that were previously undetected. We identified significant quantitative trait loci linked to milk protein and fat percentages. By comparing the selection signatures identified in this study with those reported in the Bovine Genome Variation Database, we discovered 38 novel genes under selection that were not identified through traditional methods. These genes are primarily associated with milk and meat yield and quality. Our findings enhance our understanding of spatially varying selection's impact on bovine genomic diversity, laying a foundation for future research in genetic improvement and conservation. This is the first deep learning-based study of selection signatures in cattle, offering new insights for evolutionary and livestock genomics research.

Genome-wide association study (GWAS) with single nucleotide polymorphisms (SNPs) has been widely used to explore genetic controls of phenotypic traits. Alternatively, GWAS can use counts of substrings of length k from longer sequencing reads, k-mers, as genotyping data. Using maize cob and kernel color traits, we demonstrated that k-mer GWAS can effectively identify associated k-mers. Co-expression analysis of kernel color k-mers and genes directly found k-mers from known causal genes. Analyzing complex traits of kernel oil and leaf angle resulted in k-mers from both known and candidate genes. A gene encoding a MADS transcription factor was functionally validated by showing that ectopic expression of the gene led to less upright leaves. Evolution analysis revealed most k-mers positively correlated with kernel oil were strongly selected against in maize populations, while most k-mers for upright leaf angle were positively selected. In addition, genomic prediction of kernel oil, leaf angle, and flowering time using k-mer data resulted in a similarly high prediction accuracy to the standard SNP-based method. Collectively, we showed k-mer GWAS is a powerful approach for identifying trait-associated genetic elements. Further, our results demonstrated the bridging role of k-mers for data integration and functional gene discovery., (© 2024 Society for Experimental Biology and John Wiley & Sons Ltd.)

Our study utilized genome-wide association studies (GWAS) to link nucleotide variants to traits in Populus trichocarpa, a species with rapid linkage disequilibrium decay. The aim was to overcome the challenge of interpreting statistical associations at individual loci without sufficient biological context, which often leads to reliance solely on gene annotations from unrelated model organisms. We employed an integrative approach that included GWAS targeting multiple traits using three individual techniques for lignocellulose phenotyping, expression quantitative trait loci (eQTL) analysis to construct transcriptional regulatory networks around each candidate locus and co-expression analysis to provide biological context for these networks, using lignocellulose biosynthesis in Populus trichocarpa as a case study. The research identified three candidate genes potentially involved in lignocellulose formation, including one previously recognized gene (Potri.005G116800/VND1, a critical regulator of secondary cell wall formation) and two genes (Potri.012G130000/AtSAP9 and Potri.004G202900/BIC1) with newly identified putative roles in lignocellulose biosynthesis. Our integrative approach offers a framework for providing biological context to loci associated with trait variation, facilitating the discovery of new genes and regulatory networks., (© 2024 The Author(s). New Phytologist © 2024 New Phytologist Foundation.)

Epistasis refers to nonallelic interaction between genes that cause bias in estimates of genetic parameters for a phenotype with interactions of two or more genes affecting the same trait. Partitioning of epistatic effects allows true estimation of the genetic parameters affecting phenotypes. Multigenic variation plays a central role in the evolution of complex characteristics, among which pleiotropy, where a single gene affects several phenotypic characters, has a large influence. While pleiotropic interactions provide functional specificity, they increase the challenge of gene discovery and functional analysis. Overcoming pleiotropy-based phenotypic trade-offs offers potential for assisting breeding for complex traits. Modelling higher order nonallelic epistatic interaction, pleiotropy and non-pleiotropy-induced variation, and genotype × environment interaction in genomic selection may provide new paths to increase the productivity and stress tolerance for next generation of crop cultivars. Advances in statistical models, software and algorithm developments, and genomic research have facilitated dissecting the nature and extent of pleiotropy and epistasis. We overview emerging approaches to exploit positive (and avoid negative) epistatic and pleiotropic interactions in a plant breeding context, including developing avenues of artificial intelligence, novel exploitation of large-scale genomics and phenomics data, and involvement of genes with minor effects to analyse epistatic interactions and pleiotropic quantitative trait loci, including missing heritability., (© 2024 The Author(s). Plant Biotechnology Journal published by Society for Experimental Biology and The Association of Applied Biologists and John Wiley & Sons Ltd.)

Developing early maturing lentil has the potential to minimize yield losses, mainly during terminal drought. Whole-genome resequencing (WGRS) based QTL-seq identified the loci governing earliness in lentil. The genetic analysis for maturity duration provided a good fit to 3:1 segregation (F 2 ), indicating earliness as a recessive trait. WGRS of Globe Mutant (late parent), late-flowering, and early-flowering bulks (from RILs) has generated 1124.57, 1052.24 million raw and clean reads, respectively. The QTL-Seq identified three QTLs (LcqDTF3.1, LcqDTF3.2, and LcqDTF3.3) on chromosome 3 having 246244 SNPs and 15577 insertions/deletions (InDels) and 13 flowering pathway genes. Of these, 11 exhibited sequence variations between bulks and validation (qPCR) revealed a significant difference in the expression of nine candidate genes (LcGA20oxG, LcFRI, LcLFY, LcSPL13a, Lcu.2RBY.3g060720, Lcu.2RBY.3g062540, Lcu.2RBY.3g062760, LcELF3a, and LcEMF1). Interestingly, the LcELF3a gene showed significantly higher expression in late-flowering genotype and exhibited substantial involvement in promoting lateness. Subsequently, an InDel marker (I-SP-383.9; LcELF3a gene) developed from LcqDTF3.2 QTL region showed 82.35% PVE (phenotypic variation explained) for earliness. The cloning, sequencing, and comparative analysis of the LcELF3a gene from both parents revealed 23 SNPs and InDels. Interestingly, a 52 bp deletion was recorded in the LcELF3a gene of L4775, predicted to cause premature termination of protein synthesis after 4 missense amino acids beyond the 351st amino acid due to the frameshift during translation. The identified InDel marker holds significant potential for breeding early maturing lentil varieties., (© 2024 The Author(s). Plant Biotechnology Journal published by Society for Experimental Biology and The Association of Applied Biologists and John Wiley & Sons Ltd.)

Castor (Ricinus communis L.) is an importance crop cultivated for its oil and economic value. Seed size is a crucial factor that determines crop yield. Gaining insight into the molecular regulatory processes of seed development is essential for the genetic enhancement and molecular breeding of castor. Here, we successfully fine-mapped a major QTL related to seed size, qSS3, to a 180 kb interval on chromosome 03 using F2 populations (DL01×WH11). A 17.6-kb structural variation (SV) was detected through genomic comparison between DL01 and WH11. Analysis of haplotypes showed that the existence of the complete 17.6 kb structural variant may lead to the small seed characteristic in castor. In addition, we found that qSS3 contains the microRNA396b (miR396b) sequence, which is situated within the 17.6 kb SV. The results of our experiment offer additional evidence that miR396-Growth Regulating Factor 4 (GRF4) controls seed size by impacting the growth and multiplication of seed coat and endosperm cells. Furthermore, we found that RcGRF4 activates the expression of YUCCA6 (YUC6), facilitating the production of IAA in seeds and thereby impacting the growth of castor seeds. Our research has discovered a crucial functional module that controls seed size, offering a fresh understanding of the mechanism underlying seed size regulation in castor., Competing Interests: Conflict of interest statement. The authors declare that there is no conflict of interest., (© The Author(s) 2024. Published by Oxford University Press on behalf of American Society of Plant Biologists. All rights reserved. For commercial re-use, please contact reprints@oup.com for reprints and translation rights for reprints. All other permissions can be obtained through our RightsLink service via the Permissions link on the article page on our site—for further information please contact journals.permissions@oup.com.)

Wheat (Triticum aestivum L.) is one of the most important crops worldwide and a major source of human cadmium (Cd) intake. Limiting grain Cd concentration (Gr&#95;Cd&#95;Conc) in wheat is necessary to ensure food safety. However, the genetic factors associated with Cd uptake, translocation and distribution and Gr&#95;Cd&#95;Conc in wheat are poorly understood. Here, we mapped quantitative trait loci (QTLs) for Gr&#95;Cd&#95;Conc and its related transport pathway using a recombinant inbred line (RIL) population derived from 2 Polish wheat varieties (RIL&#95;DT; dwarf Polish wheat [DPW] and tall Polish wheat [TPW]). We identified 29 novel major QTLs for grain and tissue Cd concentration; 14 novel major QTLs for Cd uptake, translocation, and distribution; and 27 major QTLs for agronomic traits. We also analyzed the pleiotropy of these QTLs. Six novel QTLs (QGr&#95;Cd&#95;Conc-1A, QGr&#95;Cd&#95;Conc-3A, QGr&#95;Cd&#95;Conc-4B, QGr&#95;Cd&#95;Conc-5B, QGr&#95;Cd&#95;Conc-6A, and QGr&#95;Cd&#95;Conc-7A) for Gr&#95;Cd&#95;Conc explained 8.16% to 17.02% of the phenotypic variation. QGr&#95;Cd&#95;Conc-3A, QGr&#95;Cd&#95;Conc-6A, and QGr&#95;Cd&#95;Conc-7A pleiotropically regulated Cd transport; 3 other QTLs were organ-specific for Gr&#95;Cd&#95;Conc. We fine-mapped the locus of QGr&#95;Cd&#95;Conc-4B and identified the candidate gene as Cation/Ca exchanger 2 (TpCCX2-4B), which was differentially expressed in DPW and TPW. It encodes an endoplasmic reticulum membrane/plasma membrane-localized Cd efflux transporter in yeast. Overexpression of TpCCX2-4B reduced Gr&#95;Cd&#95;Conc in rice. The average Gr&#95;Cd&#95;Conc was significantly lower in TpCCX2-4BDPW genotypes than in TpCCX2-4BTPW genotypes of the RIL&#95;DT population and 2 other natural populations, based on a Kompetitive allele-specific PCR marker derived from the different promoter sequences between TpCCX2-4BDPW and TpCCX2-4BTPW. Our study reveals the genetic mechanism of Cd accumulation in wheat and provides valuable resources for genetic improvement of low-Cd-accumulating wheat cultivars., Competing Interests: Conflict of interest statement. None declared., (© The Author(s) 2024. Published by Oxford University Press on behalf of American Society of Plant Biologists. All rights reserved. For commercial re-use, please contact reprints@oup.com for reprints and translation rights for reprints. All other permissions can be obtained through our RightsLink service via the Permissions link on the article page on our site—for further information please contact journals.permissions@oup.com.)

Rapeseed (Brassica napus) silique is the major carbohydrate source for seed development, and the final silique length has attracted great attention from breeders. However, no studies had focused on the dynamic character of silique elongation length (SEL). Here, the dynamic SEL investigation in a natural population including 588 lines over two years indicate that dynamic SEL during 0-20 days after flowering was the most essential stage associated with seed number per silique (SPS) and thousand seed weight (TSW). Then, nine loci were identified to be associated with SEL based on GWAS analysis, among which five SNPs (over 50%) distributed on the A02 chromosome within 6.08 to 6.48 Mb. Subsequently, we screened 5078 differentially expressed genes between two extreme materials. An unknown protein, BnaA02.SE, was identified combining with GWAS and RNA-Seq analysis. Subcellular localization and expression profiles analysis demonstrated that BnaA02.SE is a chloroplast- and nucleus-localized protein mainly expressed in pericarps and leaves. Furthermore, transgenic verification and dynamic cytological observation reveal that overexpressed BnaA02.SE can promote silique elongation by regulating JA and IAA contents, affecting cell proliferation and expansion, respectively, and finally enhance seed yield by influencing SPS and TSW. Haplotype analysis reveal that the homologs of BnaA02.SE may also be involved in silique elongation regulation. Our findings provided comprehensive insights into a newly SEL trait, and cloned the first gene (BnaA02.SE) controlling silique elongation in B. napus. The identified BnaA02.SE and its homologs can offer a valuable target for improving B. napus yield., (© 2024 The Author(s). Plant Biotechnology Journal published by Society for Experimental Biology and The Association of Applied Biologists and John Wiley & Sons Ltd.)

Chronotype is a proxy sleep measure that has been associated with neuropsychiatric disorders. By investigating how chronotype influences risk for neuropsychiatric disorders and vice versa, we may identify modifiable risk factors for each phenotype. Here we used Mendelian randomization (MR), to explore causal effects by (1) studying the causal relationships between neuropsychiatric disorders and chronotype and (2) characterizing the genetic components of these phenotypes. Firstly, we investigated if a causal role exists between five neuropsychiatric disorders and chronotype using the largest genome-wide association studies (GWAS) available. Secondly, we integrated data from expression quantitative trait loci (eQTLs) to investigate the role of gene expression alterations on these phenotypes. Evening chronotype was causal for increased risk of schizophrenia and autism spectrum disorder and schizophrenia was causal for a tendency toward evening chronotype. We identified 12 eQTLs where gene expression changes in brain or blood were causal for one of the phenotypes, including two eQTLs for SNX19 in hippocampus and hypothalamus that were causal for schizophrenia. These findings provide important evidence for the complex, bidirectional relationship that exists between a sleep-based phenotype and neuropsychiatric disorders, and use gene expression data to identify causal roles for genes at associated loci., (© 2024 The Authors. American Journal of Medical Genetics Part B: Neuropsychiatric Genetics published by Wiley Periodicals LLC.)

Rice grain number is a crucial agronomic trait impacting yield. In this study, we characterized a quantitative trait locus (QTL), GRAIN NUMBER 1.1 (GN1.1), which encodes a Flowering Locus T-like1 (FT-L1) protein and acts as a negative regulator of grain number in rice. The elite allele GN1.1 B , derived from the Oryza indica variety, BF3-104, exhibits a 14.6% increase in grain yield compared with the O. japonica variety, Nipponbare, based on plot yield tests. We demonstrated that GN1.1 interacted with and enhanced the stability of ADP-ribosylation factor (Arf)-GTPase-activating protein (Gap), OsZAC. Loss of function of OsZAC results in increased grain number. Based on our data, we propose that GN1.1 B facilitates the elevation of auxin content in young rice panicles by affecting polar auxin transport (PAT) through interaction with OsZAC. Our study unveils the pivotal role of the GN1.1 locus in rice panicle development and presents a novel, promising allele for enhancing rice grain yield through genetic improvement., (© 2024 Institute of Botany, Chinese Academy of Sciences.)

Background: Understanding the genetic control of pod shatter resistance and its association with pod length is crucial for breeding improved pod shatter resistance and reducing pre-harvest yield losses due to extensive shattering in cultivars of Brassica species. In this study, we evaluated a doubled haploid (DH) mapping population derived from an F 1 cross between two Brassica carinata parental lines Y-BcDH64 and W-BcDH76 (YWDH), originating from Ethiopia and determined genetic bases of variation in pod length and pod shatter resistance, measured as rupture energy. The YWDH population, its parental lines and 11 controls were grown across three years for genetic analysis., Results: By using three quantitative trait loci (QTL) analytic approaches, we identified nine genomic regions on B02, B03, B04, B06, B07 and C01 chromosomes for rupture energy that were repeatedly detected across three growing environments. One of the QTL on chromosome B07, flanked with DArTseq markers 100,046,735 and 100,022,658, accounted for up to 27.6% of genetic variance in rupture energy. We observed no relationship between pod length and rupture energy, suggesting that pod length does not contribute to variation in pod shatter resistance. Comparative mapping identified six candidate genes; SHP1 on B6, FUL and MAN on chromosomes B07, IND and NST2 on B08, and MAN7 on C07 that mapped within 0.2 Mb from the QTL for rupture energy., Conclusion: The results suggest that favourable alleles of stable QTL on B06, B07, B08 and C01 for pod shatter resistance can be incorporated into the shatter-prone B. carinata and its related species to improve final seed yield at harvest., (© 2024. Crown.)

Background: Castor (Ricinus communis L., 2n = 2x = 20) is an important industrial crop, due to its oil is very important to the global special chemical industry. Seed size and seed weight are fundamentally important in determining castor yield, while little is known about it. In this study, QTL analysis and candidate gene mining of castor seed size and seed weight were conducted with composite interval mapping (CIM), inclusive composite interval mapping (ICIM) and marker enrichment strategy in 4 populations, i.e., populations F 2 , BC 1 , S 1 -1 and S 1 -2, derived from 2 accessions with significant phenotypic differences., Results: In the QTL primary mapping, 2 novel QTL clusters were detected in marker intervals RCM520-RCM76 and RCM915-RCM950. In order to verify their accuracy and to narrow their intervals, QTL remapping was carried out in populations F 2 and BC 1 . Among them, 44 and 30 QTLs underlying seed size and seed weight were detected in F 2 population using methods CIM and ICIM-ADD respectively, including 4-9 and 3-5 ones conferring each trait were identified with a phenotypic variation explained ranged from 37.92 to 115.81% and 32.86-45.98% respectively. The remapping results in BC 1 population were consistent with those in F 2 population. Importantly, 3 QTL clusters (i.e. QTL-cluster1, QTL-cluster2 and QTL-cluster3) were found in marker intervals RCM74-RCM76 (37.1 kb), RCM930-RCM950 (259.8 kb) and RCM918-RCM920 (172.9 kb) respectively; in addition, all of them were detected again, the former one was found in the S 1 -2 population, and the latter two were found simultaneously in the populations S 1 -1 and S 1 -2. Finally, 6 candidate genes (i.e. LOC8266555, LOC8281168, LOC8281151, LOC8259066, LOC8258591 and LOC8270077) were screened in the above QTL clusters, they were differentially expressed in multiple seed tissues of both parents, signifying the potential role in regulating seed size and seed weight., Conclusion: The above results not only provide new insights into the genetic structure of seed size and seed weight in castor, but also lay the foundation for the functional identification of these candidate genes., (© 2024. The Author(s).)

COVID-19 casualties vary among different ancestral groups due to a variety of factors. Here, we present a protocol for analyzing publicly available genome-wide association studies (GWASs) to search for ancestry-specific genetic factors related to severe COVID-19. We describe steps for downloading and comparing two COVID-19 GWASs, calculating expression quantitative trait loci, and single-cell gene expression analysis. We demonstrate this approach using GWASs from Host Genetics Initiative; however, it is applicable to other databases such as the UK Biobank. For complete details on the use and execution of this protocol, please refer to Cheng et al. 1 ., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

Background: Seed vigor recognized as a quantitative trait is of particular importance for agricultural production. However, limited knowledge is available for understanding genetic basis of wheat seed vigor., Methods: The aim of this study was to identify quantitative trait loci (QTL) responsible for 10 seed vigor-related traits representing multiple aspects of seed-vigor dynamics during artificial aging with 6 different treatment times (0, 24, 36, 48, 60, and 72 h) under controlled conditions (48 °C, 95% humidity, and dark). The mapping populations were two wheat introgression lines (IL-1 and IL-2) derived from recipient parent (Lumai 14) and donor parent (Shaanhan 8675 or Jing 411)., Results: A total of 26 additive QTLs and 72 pairs of epistatic QTLs were detected for wheat seed-vigor traits. Importantly, chromosomes 1B and 7B contained several co-located QTLs, and chromosome 2A had a QTL-rich region near the marker Xwmc667, indicating that these QTLs may affect wheat seed vigor with pleiotropic effects. Furthermore, several possible consistent QTLs (hot-spot regions) were examined by comparison analysis of QTLs detected in this study and reported previously. Finally, a set of candidate genes for wheat seed vigor were predicted to be involved in transcription regulation, carbohydrate and lipid metabolism., Conclusion: The present findings lay new insights into the mechanism underlying wheat seed vigor, providing valuable information for wheat genetic improvement especially marker-assisted breeding to increase seed vigor and consequently achieve high grain yield despite of further investigation required., Competing Interests: The authors declare there are no competing interests., (©2024 Yang et al.)

Multilocus genome-wide association study has become the state-of-the-art tool for dissecting the genetic architecture of complex and multiomic traits. However, most existing multilocus methods require relatively long computational time when analyzing large datasets. To address this issue, in this study, we proposed a fast mrMLM method, namely, best linear unbiased prediction multilocus random-SNP-effect mixed linear model (BLUPmrMLM). First, genome-wide single-marker scanning in mrMLM was replaced by vectorized Wald tests based on the best linear unbiased prediction (BLUP) values of marker effects and their variances in BLUPmrMLM. Then, adaptive best subset selection (ABESS) was used to identify potentially associated markers on each chromosome to reduce computational time when estimating marker effects via empirical Bayes. Finally, shared memory and parallel computing schemes were used to reduce the computational time. In simulation studies, BLUPmrMLM outperformed GEMMA, EMMAX, mrMLM, and FarmCPU as well as the control method (BLUPmrMLM with ABESS removed), in terms of computational time, power, accuracy for estimating quantitative trait nucleotide positions and effects, false positive rate, false discovery rate, false negative rate, and F1 score. In the reanalysis of two large rice datasets, BLUPmrMLM significantly reduced the computational time and identified more previously reported genes, compared with the aforementioned methods. This study provides an excellent multilocus model method for the analysis of large-scale and multiomic datasets. The software mrMLM v5.1 is available at BioCode (https://ngdc.cncb.ac.cn/biocode/tool/BT007388) or GitHub (https://github.com/YuanmingZhang65/mrMLM)., (© The Author(s) 2024. Published by Oxford University Press and Science Press on behalf of the Beijing Institute of Genomics, Chinese Academy of Sciences / China National Center for Bioinformation and Genetics Society of China.)

Main Conclusion: The Ustilaginoidea virens -rice pathosystem has been used as a model for flower-infecting fungal pathogens. The molecular biology of the interactions between U. virens and rice, with an emphasis on the attempt to get a deeper comprehension of the false smut fungus's genomes, proteome, host range, and pathogen biology, has been investigated. Meta-QTL analysis was performed to identify potential QTL hotspots for use in marker-assisted breeding. The Rice False Smut (RFS) caused by the fungus Ustilaginoidea virens currently threatens rice cultivators across the globe. RFS infects rice panicles, causing a significant reduction in grain yield. U. virens can also parasitize other hosts though they play only a minor role in its life cycle. Furthermore, because it produces mycotoxins in edible rice grains, it puts both humans and animals at risk of health problems. Although fungicides are used to control the disease, some fungicides have enabled the pathogen to develop resistance, making its management challenging. Several QTLs have been reported but stable gene(s) that confer RFS resistance have not been discovered yet. This review offers a comprehensive overview of the pathogen, its virulence mechanisms, the genome and proteome of U. virens, and its molecular interactions with rice. In addition, information has been compiled on reported resistance QTLs, facilitating the development of a consensus genetic map using meta-QTL analysis for identifying potential QTL hotspots. Finally, this review highlights current developments and trends in U. virens-rice pathosystem research while identifying opportunities for future investigations., (© 2024. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Soluble sugar and organic acids are key determinants of fruit organoleptic quality and directly affect the commodity value and economic returns of fruit crops. We performed whole-genome sequencing of the apple varieties Gala and Xiahongrou, along with their F1 hybrids, to construct a high-density bin map. Our quantitative genetic analysis pinpointed 53 quantitative trait loci (QTLs) related to 11 sugar and acid traits. We identified a candidate gene, MdNADP-ME, responsible for malate degradation, in a stable QTL on linkage group 15. Sequence analysis revealed an A/C SNP in the promoter region (MEp-799) that influences binding of the MdMYB2 transcription factor, thereby affecting MdNADP-ME expression. In our study of various apple genotypes, this SNP has been demonstrated to be linked to malate and fructose levels. We also developed a dCAPS marker associated with fruit fructose content. These results substantiate the role of MdNADP-ME in maintaining the equilibrium between sugar and acid contents in apple fruits., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

Grain and flag leaf size are two important agronomic traits that influence grain yield in rice (Oryza sativa). Many quantitative trait loci (QTLs) and genes that regulate these traits individually have been identified, however, few QTLs and genes that simultaneously control these two traits have been identified. In this study, we conducted a genome-wide association analysis in rice and detected a major locus, WIDTH OF LEAF AND GRAIN (WLG), that was associated with both grain and flag leaf width. WLG encodes a RING-domain E3 ubiquitin ligase. WLGhap.B, which possesses five single nucleotide polymophysim (SNP) variations compared to WLGhap.A, encodes a protein with enhanced ubiquitination activity that confers increased rice leaf width and grain size, whereas mutation of WLG leads to narrower leaves and smaller grains. Both WLGhap.A and WLGhap.B interact with LARGE2, a HETC-type E3 ligase, however, WLGhap.B exhibits stronger interaction with LARGE2, thus higher ubiquitination activity toward LARGE2 compared with WLGhap.A. Lysine1021 is crucial for the ubiquitination of LARGE2 by WLG. Loss-of-function of LARGE2 in wlg-1 phenocopies large2-c in grain and leaf width, suggesting that WLG acts upstream of LARGE2. These findings reveal the genetic and molecular mechanism by which the WLG-LARGE2 module mediates grain and leaf size in rice and suggest the potential of WLGhap.B in improving rice yield., Competing Interests: Conflict of interest statement. None declared., (© The Author(s) 2024. Published by Oxford University Press on behalf of American Society of Plant Biologists. All rights reserved. For commercial re-use, please contact reprints@oup.com for reprints and translation rights for reprints. All other permissions can be obtained through our RightsLink service via the Permissions link on the article page on our site—for further information please contact journals.permissions@oup.com.)

Plant aquaporins are involved in numerous physiological processes, such as cellular homeostasis, tissue hydraulics, transpiration, and nutrient supply, and are key players of the response to environmental cues. While varying expression patterns of aquaporin genes have been described across organs, developmental stages, and stress conditions, the underlying regulation mechanisms remain elusive. Hence, this work aimed to shed light on the expression variability of 4 plasma membrane intrinsic protein (PIP) genes in maize (Zea mays) leaves, and its genetic causes, through expression quantitative trait locus (eQTL) mapping across a 252-hybrid diversity panel. Significant genetic variability in PIP transcript abundance was observed to different extents depending on the isoforms. The genome-wide association study mapped numerous eQTLs, both local and distant, thus emphasizing the existing natural diversity of PIP gene expression across the studied panel and the potential to reveal regulatory actors and mechanisms. One eQTL associated with PIP2;5 expression variation was characterized. Genomic sequence comparison and in vivo reporter assay attributed, at least partly, the local eQTL to a transposon-containing polymorphism in the PIP2;5 promoter. This work paves the way to the molecular understanding of PIP gene regulation and its possible integration into larger networks regulating physiological and stress adaptation processes., Competing Interests: Conflict of interest statement. None declared., (© The Author(s) 2024. Published by Oxford University Press on behalf of American Society of Plant Biologists.)