Genes controlling hydrolysate toxin tolerance identified by QTL analysis of the natural Saccharomyces cerevisiae BCC39850

Lignocellulosic material can be converted to valorized products such as fuels. Pretreatment is an essential step in conversion, which is needed to increase the digestibility of the raw material for microbial fermentation. However, pretreatment generates by-products (hydrolysate toxins) that are detrimental to microbial growth. In this study, natural Saccharomyces strains isolated from habitats in Thailand were screened for their tolerance to synthetic hydrolysate toxins (synHTs). The Saccharomyces cerevisiae natural strain BCC39850 (toxin-tolerant) was crossed with the laboratory strain CEN.PK2-1C (toxin-sensitive), and quantitative trait locus (QTL) analysis was performed on the segregants using phenotypic scores of growth (OD600) and glucose consumption. VMS1, DET1, KCS1, MRH1, YOS9, SYO1, and YDR042C were identified from QTLs as candidate genes associated with the tolerance trait. CEN.PK2-1C knockouts of the VMS1, YOS9, KCS1, and MRH1 genes exhibited significantly greater hydrolysate toxin sensitivity to growth, whereas CEN.PK2-1C knock-ins with replacement of VMS1 and MRH1 genes from the BCC39850 alleles showed significant increased ethanol production titers compared with the CEN.PK2-1C parental strain in the presence of synHTs. The discovery of VMS1, YOS9, MRH1, and KCS1 genes associated with hydrolysate toxin tolerance in S. cerevisiae indicates the roles of the endoplasmic-reticulum-associated protein degradation pathway, plasma membrane protein association, and the phosphatidylinositol signaling system in this trait. KEY POINTS: • QTL analysis was conducted using a hydrolysate toxin-tolerant S. cerevisiae natural strain • Deletion of VMS1, YOS9, MRH1, and KCS1 genes associated with hydrolysate toxin-sensitivity • Replacement of VMS1 and MRH1 with natural strain alleles increased ethanol production titers in the presence of hydrolysate toxins.

Critical evaluation of the reliability of DNA methylation probes on the Illumina MethylationEPIC v1.0 BeadChip microarrays

DNA methylation (DNAm) plays a crucial role in a number of complex diseases. However, the reliability of DNAm levels measured using Illumina arrays varies across different probes. Previous research primarily assessed probe reliability by comparing duplicate samples between the 450k-450k or 450k-EPIC platforms, with limited investigations on Illumina EPIC v1.0 arrays. We conducted a comprehensive assessment of the EPIC v1.0 array probe reliability using 69 blood DNA samples, each measured twice, generated by the Alzheimer's Disease Neuroimaging Initiative study. We observed higher reliability in probes with average methylation beta values of 0.2 to 0.8, and lower reliability in type I probes or those within the promoter and CpG island regions. Importantly, we found that probe reliability has significant implications in the analyses of Epigenome-wide Association Studies (EWAS). Higher reliability is associated with more consistent effect sizes in different studies, the identification of differentially methylated regions (DMRs) and methylation quantitative trait locus (mQTLs), and significant correlations with downstream gene expression. Moreover, blood DNAm measurements obtained from probes with higher reliability are more likely to show concordance with brain DNAm measurements. Our findings, which provide crucial reliability information for probes on the EPIC v1.0 array, will serve as a valuable resource for future DNAm studies.

Genetic analysis of milk citrate predicted by milk mid-infrared spectra of Holstein cows in early lactation

Milk citrate is regarded as an early biomarker of negative energy balance in dairy cows during early lactation and serves as a suitable candidate phenotype for genomic selection due to its wide availability across a large number of cows through milk mid-infrared spectra prediction. However, its genetic background is not well known. Therefore, the objectives of this study were to (1) analyze the genetic parameters of milk citrate; (2) identify genomic regions associated with milk citrate; and (3) analyze the functional annotation of candidate genes and quantitative trait loci (QTL) related to milk citrate in Walloon Holstein cows. In total, 134,517 test-day milk-citrate phenotypes (mmol/L) collected within the first 50 d in milk on 52,198 Holstein cows were used. These milk-citrate phenotypes, predicted by milk mid-infrared spectra, were divided into 3 traits according to the first (citrate1), second (citrate2), and third to fifth parity (citrate3+). Genomic information for 566,170 SNPs was available for 4,479 animals. A multiple-trait repeatability model was used to estimate genetic parameters. A single-step GWAS was used to identify candidate genes for citrate and post-GWAS analysis was done to investigate the relationship and function of the identified candidate genes. The heritabilities estimated for citrate1, citrate2, and citrate3+ were 0.40, 0.37, and 0.35, respectively. The genetic correlations among the 3 traits ranged from 0.98 to 0.99. The genomic correlations among the 3 traits were also close to 1.00 across the genomic regions (1 Mb) in the whole genome, which means that citrate can be considered as a single trait in the first 5 parities. In total, 603 significant SNPs located on 3 genomic regions (chromosome 7, 68.569-68.575 Mb; chromosome 14, 0.15-1.90 Mb; and chromosome 20, 54.00-64.28 Mb), were identified to be associated with milk citrate. We identified 89 candidate genes including GPT, ANKH, PPP1R16A, and 32 QTL reported in the literature related to the identified significant SNPs. These identified QTL were mainly reported associated with milk fatty acids and metabolic diseases in dairy cows. This study suggests that milk citrate in Holstein cows is highly heritable and has the potential to be used as an early proxy for the negative energy balance of Holstein cows in a breeding objective.

Exploring genomic inbreeding and selection signatures in a commercial Brangus herd through functional annotation

Composite breeds, including Brangus, are widely utilized in subtropical and tropical regions to harness the advantages of both Bos t. taurus and Bos t. indicus breeds. The formation and subsequent selection of composite breeds may result in discernible signatures of selection and shifts in genomic population structure. The objectives of this study were to 1) assess genomic inbreeding, 2) identify signatures of selection, 3) assign functional roles to these signatures in a commercial Brangus herd, and 4) contrast signatures of selection between selected and non-selected cattle from the same year. A total of 4035 commercial Brangus cattle were genotyped using the GGP-F250K array. Runs of Homozygosity (ROH) were used to identify signatures of selection and calculate genomic inbreeding. Quantitative trait loci (QTL) enrichment analysis and literature search identified phenotypic traits linked to ROH islands. Genomic inbreeding averaged 5%, primarily stemming from ancestors five or more generations back. A total of nine ROH islands were identified, QTL enrichment analysis revealed traits related to growth, milk composition, carcass, reproductive, and meat quality traits. Notably, the ROH island on BTA14 encompasses the pleiomorphic adenoma (PLAG1) gene, which has been linked to growth, carcass, and reproductive traits. Moreover, ROH islands associated with milk yield and composition were more pronounced in selected replacement heifers of the population, underscoring the importance of milk traits in cow-calf production. In summary, our research sheds light on the changing genetic landscape of the Brangus breed due to selection pressures and reveals key genomic regions impacting production traits.

Haplotype function score improves biological interpretation and cross-ancestry polygenic prediction of human complex traits

We propose a new framework for human genetic association studies: at each locus, a deep learning model (in this study, Sei) is used to calculate the functional genomic activity score for two haplotypes per individual. This score, defined as the Haplotype Function Score (HFS), replaces the original genotype in association studies. Applying the HFS framework to 14 complex traits in the UK Biobank, we identified 3619 independent HFS-trait associations with a significance of p < 5 × 10-8. Fine-mapping revealed 2699 causal associations, corresponding to a median increase of 63 causal findings per trait compared with single-nucleotide polymorphism (SNP)-based analysis. HFS-based enrichment analysis uncovered 727 pathway-trait associations and 153 tissue-trait associations with strong biological interpretability, including 'circadian pathway-chronotype' and 'arachidonic acid-intelligence'. Lastly, we applied least absolute shrinkage and selection operator (LASSO) regression to integrate HFS prediction score with SNP-based polygenic risk scores, which showed an improvement of 16.1-39.8% in cross-ancestry polygenic prediction. We concluded that HFS is a promising strategy for understanding the genetic basis of human complex traits.

GbyE: an integrated tool for genome widely association study and genome selection based on genetic by environmental interaction

BACKGROUND

The growth and development of organism were dependent on the effect of genetic, environment, and their interaction. In recent decades, lots of candidate additive genetic markers and genes had been detected by using genome-widely association study (GWAS). However, restricted to computing power and practical tool, the interactive effect of markers and genes were not revealed clearly. And utilization of these interactive markers is difficult in the breeding and prediction, such as genome selection (GS).

RESULTS

Through the Power-FDR curve, the GbyE algorithm can detect more significant genetic loci at different levels of genetic correlation and heritability, especially at low heritability levels. The additive effect of GbyE exhibits high significance on certain chromosomes, while the interactive effect detects more significant sites on other chromosomes, which were not detected in the first two parts. In prediction accuracy testing, in most cases of heritability and genetic correlation, the majority of prediction accuracy of GbyE is significantly higher than that of the mean method, regardless of whether the rrBLUP model or BGLR model is used for statistics. The GbyE algorithm improves the prediction accuracy of the three Bayesian models BRR, BayesA, and BayesLASSO using information from genetic by environmental interaction (G × E) and increases the prediction accuracy by 9.4%, 9.1%, and 11%, respectively, relative to the Mean value method. The GbyE algorithm is significantly superior to the mean method in the absence of a single environment, regardless of the combination of heritability and genetic correlation, especially in the case of high genetic correlation and heritability.

CONCLUSIONS

Therefore, this study constructed a new genotype design model program (GbyE) for GWAS and GS using Kronecker product. which was able to clearly estimate the additive and interactive effects separately. The results showed that GbyE can provide higher statistical power for the GWAS and more prediction accuracy of the GS models. In addition, GbyE gives varying degrees of improvement of prediction accuracy in three Bayesian models (BRR, BayesA, and BayesCpi). Whatever the phenotype were missed in the single environment or multiple environments, the GbyE also makes better prediction for inference population set. This study helps us understand the interactive relationship between genomic and environment in the complex traits. The GbyE source code is available at the GitHub website ( https://github.com/liu-xinrui/GbyE ).

QTL mapping for the flag leaf-related traits using RILs derived from Trititrigia germplasm line SN304 and wheat cultivar Yannong15 in multiple environments

BACKGROUND

Developing and enriching genetic resources plays important role in the crop improvement. The flag leaf affects plant architecture and contributes to the grain yield of wheat (Triticum aestivum L.). The genetic improvement of flag leaf traits faces problems such as a limited genetic basis. Among the various genetic resources of wheat, Thinopyrum intermedium has been utilized as a valuable resource in genetic improvement due to its disease resistance, large spikes, large leaves, and multiple flowers. In this study, a recombinant inbred line (RIL) population was derived from common wheat Yannong15 and wheat-Th. intermedium introgression line SN304 was used to identify the quantitative trait loci (QTL) for flag leaf-related traits.

RESULTS

QTL mapping was performed for flag leaf length (FLL), flag leaf width (FLW) and flag leaf area (FLA). A total of 77 QTLs were detected, and among these, 51 QTLs with positive alleles were contributed by SN304. Fourteen major QTLs for flag leaf traits were detected on chromosomes 2B, 3B, 4B, and 2D. Additionally, 28 QTLs and 8 QTLs for flag leaf-related traits were detected in low-phosphorus and drought environments, respectively. Based on major QTLs of positive alleles from SN304, we identified a pair of double-ended anchor primers mapped on chromosome 2B and amplified a specific band of Th. intermedium in SN304. Moreover, there was a major colocated QTL on chromosome 2B, called QFll/Flw/Fla-2B, which was delimited to a physical interval of approximately 2.9 Mb and contained 20 candidate genes. Through gene sequence and expression analysis, four candidate genes associated with flag leaf formation and growth in the QTL interval were identified.

CONCLUSION

These results promote the fine mapping of QFll/Flw/Fla-2B, which have pleiotropic effects, and will facilitate the identification of candidate genes for flag leaf-related traits. Additionally, this work provides a theoretical basis for the application of Th. intermedium in wheat breeding.

Characterization of a wheat stable QTL for spike length and its genetic effects on yield-related traits

Spike length (SL) is one of the most important agronomic traits affecting yield potential and stability in wheat. In this study, a major stable quantitative trait locus (QTL) for SL, i.e., qSl-2B, was detected in multiple environments in a recombinant inbred line (RIL) mapping population, KJ-RILs, derived from a cross between Kenong 9204 (KN9204) and Jing 411 (J411). The qSl-2B QTL was mapped to the 60.06-73.06 Mb region on chromosome 2B and could be identified in multiple mapping populations. An InDel molecular marker in the target region was developed based on a sequence analysis of the two parents. To further clarify the breeding use potential of qSl-2B, we analyzed its genetic effects and breeding selection effect using both the KJ-RIL population and a natural mapping population, which consisted of 316 breeding varieties/advanced lines. The results showed that the qSl-2B alleles from KN9204 showed inconsistent genetic effects on SL in the two mapping populations. Moreover, in the KJ-RILs population, the additive effects analysis of qSl-2B showed that additive effect was higher when both qSl-2D and qSl-5A harbor negative alleles under LN and HN. In China, a moderate selection utilization rate for qSl-2B was found in the Huanghuai winter wheat area and the selective utilization rate for qSl-2B continues to increase. The above findings provided a foundation for the genetic improvement of wheat SL in the future via molecular breeding strategies.

Identification of quantitative trait loci associated with leaf rust resistance in rye by precision mapping

BACKGROUND

Leaf rust (LR) is among the most destructive fungal diseases of rye (Secale cereale L.). Despite intensive research using various analytical and methodological approaches, such as quantitative trait locus (QTL) mapping, candidate gene expression analysis, and transcriptome sequencing, the genetic basis of the rye immune response to LR remains unclear.

RESULTS

A genome-wide association study was employed to detect QTLs controlling the immune response to LR of rye. A mapping population, G38A, was constructed by crossing two inbred lines: 723 (susceptible to LR) and JKI-NIL-Pr3 (a donor of the LR resistance gene Pr3). For genotyping, SNP-DArT and silico-DArT markers were used. Resistance phenotyping was conducted by visual assessment of the infection severity in detached leaf segments inoculated with two isolates of Puccinia recondita f. sp. secalis, namely, 60/17/2.1 (isolate S) in the main experiment and 86/n/2.1_5x (isolate N) in the validation experiment, at 10 and 17 days post-infection (dpi), respectively. In total, 42,773 SNP-DArT and 105,866 silico-DArT markers were included in the main analysis including isolate S, of which 129 and 140 SNP-DArTs and 767 and 776 silico-DArTs were significantly associated (p ≤ 0.001; - log10(p) ≥ 3.0) with the immune response to LR at 10 and 17 dpi, respectively. Most significant markers were mapped to chromosome 1R. The number of common markers from both systems and at both time points occupying common chromosomal positions was 37, of which 21 were positioned in genes, comprising 18 markers located in exons and three in introns. This gene pool included genes encoding proteins with a known function in response to LR (e.g., a NBS-LRR disease resistance protein-like protein and carboxyl-terminal peptidase).

CONCLUSION

This study has expanded and supplemented existing knowledge of the genetic basis of rye resistance to LR by (1) detecting two QTLs associated with the LR immune response of rye, of which one located on the long arm of chromosome 1R is newly detected, (2) assigning hundreds of markers significantly associated with the immune response to LR to genes in the 'Lo7' genome, and (3) predicting the potential translational effects of polymorphisms of SNP-DArT markers located within protein-coding genes.

Genetic mapping reveals the complex genetic architecture controlling slow canopy wilting in soybean

In soybean [Glycine max (L.) Merr.], drought stress is the leading cause of yield loss from abiotic stress in rain-fed US growing areas. Only 10% of the US soybean production is irrigated; therefore, plants must possess physiological mechanisms to tolerate drought stress. Slow canopy wilting is a physiological trait that is observed in a few exotic plant introductions (PIs) and may lead to yield improvement under drought stress. Canopy wilting of 130 recombinant inbred lines (RILs) derived from Hutcheson × PI 471938 grown under drought stress was visually evaluated and genotyped with the SoySNP6K BeadChip. Over four years, field evaluations of canopy wilting were conducted under rainfed conditions at three locations across the US (Georgia, Kansas, and North Carolina). Due to the variation in weather among locations and years, the phenotypic data were collected from seven environments. Substantial variation in canopy wilting was observed among the genotypes in the RIL population across environments. Three QTLs were identified for canopy wilting from the RIL population using composite interval mapping on chromosomes (Chrs) 2, 8, and 9 based on combined environmental analyses. These QTLs inherited the favorable alleles from PI 471938 and accounted for 11, 10, and 14% of phenotypic variation, respectively. A list of 106 candidate genes were narrowed down for these three QTLs based on the published information. The QTLs identified through this research can be used as targets for further investigation to understand the mechanisms of slow canopy wilting. These QTLs could be deployed to improve drought tolerance through a targeted selection of the genomic regions from PI 471938.

TOPMed imputed genomics enhances genomic atlas of the human proteome in brain, cerebrospinal fluid, and plasma

Comprehensive expression quantitative trait loci studies have been instrumental for understanding tissue-specific gene regulation and pinpointing functional genes for disease-associated loci in a tissue-specific manner. Compared to gene expressions, proteins more directly affect various biological processes, often dysregulated in disease, and are important drug targets. We previously performed and identified tissue-specific protein quantitative trait loci in brain, cerebrospinal fluid, and plasma. We now enhance this work by analyzing more proteins (1,300 versus 1,079) and an almost twofold increase in high quality imputed genetic variants (8.4 million versus 4.4 million) by using TOPMed reference panel. We identified 38 genomic regions associated with 43 proteins in brain, 150 regions associated with 247 proteins in cerebrospinal fluid, and 95 regions associated with 145 proteins in plasma. Compared to our previous study, this study newly identified 12 loci in brain, 30 loci in cerebrospinal fluid, and 22 loci in plasma. Our improved genomic atlas uncovers the genetic control of protein regulation across multiple tissues. These resources are accessible through the Online Neurodegenerative Trait Integrative Multi-Omics Explorer for use by the scientific community.

Fine mapping and breeding application of two brown planthopper resistance genes derived from landrace rice

The Brown planthopper (Nilaparvata lugens Stål; BPH) is known to cause significant damage to rice crops in Asia, and the use of host-resistant varieties is an effective and environmentally friendly approach for controlling BPH. However, genes limited resistance genes that are used in insect-resistant rice breeding programs, and landrace rice varieties are materials resources that carry rich and versatile genes for BPH resistance. Two landrace indica rice accessions, CL45 and CL48, are highly resistant to BPH and show obvious antibiosis against BPH. A novel resistance locus linked to markers 12M16.983 and 12M19.042 was identified, mapped to chromosome 12 in CL45, and designated Bph46. It was finely mapped to an interval of 480 kb and Gene 3 may be the resistance gene. Another resistance locus linked to markers RM26567 and 11MA104 was identified and mapped to chromosome 11 in CL48 and designated qBph11.3 according to the nominating rule. It was finely mapped to an interval of 145 kb, and LOC_Os11g29090 and LOC_Os11g29110 may be the resistance genes. Moreover, two markers, 12M16.983 and 11MA104, were developed for CL45 and CL48, respectively, using marker-assisted selection (MAS) and were confirmed by backcrossing individuals and phenotypic detection. Interestingly, we found that the black glume color is closely linked to the BPH resistance gene in CL48 and can effectively assist in the identification of positive individuals for breeding. Finally, several near-isogenic lines with a 9311 or KW genetic background, as well as pyramid lines with two resistance parents, were developed using MAS and exhibited significantly high resistance against BPHs.

DeLIVR: a deep learning approach to IV regression for testing nonlinear causal effects in transcriptome-wide association studies

Transcriptome-wide association studies (TWAS) have been increasingly applied to identify (putative) causal genes for complex traits and diseases. TWAS can be regarded as a two-sample two-stage least squares method for instrumental variable (IV) regression for causal inference. The standard TWAS (called TWAS-L) only considers a linear relationship between a gene's expression and a trait in stage 2, which may lose statistical power when not true. Recently, an extension of TWAS (called TWAS-LQ) considers both the linear and quadratic effects of a gene on a trait, which however is not flexible enough due to its parametric nature and may be low powered for nonquadratic nonlinear effects. On the other hand, a deep learning (DL) approach, called DeepIV, has been proposed to nonparametrically model a nonlinear effect in IV regression. However, it is both slow and unstable due to the ill-posed inverse problem of solving an integral equation with Monte Carlo approximations. Furthermore, in the original DeepIV approach, statistical inference, that is, hypothesis testing, was not studied. Here, we propose a novel DL approach, called DeLIVR, to overcome the major drawbacks of DeepIV, by estimating a related but different target function and including a hypothesis testing framework. We show through simulations that DeLIVR was both faster and more stable than DeepIV. We applied both parametric and DL approaches to the GTEx and UK Biobank data, showcasing that DeLIVR detected additional 8 and 7 genes nonlinearly associated with high-density lipoprotein (HDL) cholesterol and low-density lipoprotein (LDL) cholesterol, respectively, all of which would be missed by TWAS-L, TWAS-LQ, and DeepIV; these genes include BUD13 associated with HDL, SLC44A2 and GMIP with LDL, all supported by previous studies.

The genomic basis of nitrogen utilization efficiency and trait plasticity to improve nutrient stress tolerance in cultivated sunflower

Maintaining crop productivity is challenging as population growth, climate change, and increasing fertilizer costs necessitate expanding crop production to poorer lands whilst reducing inputs. Enhancing crops' nutrient use efficiency is thus an important goal, but requires a better understanding of related traits and their genetic basis. We investigated variation in low nutrient stress tolerance in a diverse panel of cultivated sunflower genotypes grown under high and low nutrient conditions, assessing relative growth rate (RGR) as performance. We assessed variation in traits related to nitrogen utilization efficiency (NUtE), mass allocation, and leaf elemental content. Across genotypes, nutrient limitation generally reduced RGR. Moreover, there was a negative correlation between vigor (RGR in control) and decline in RGR in response to stress. Given this trade-off, we focused on nutrient stress tolerance independent of vigor. This tolerance metric correlated with the change in NUtE, plasticity for a suite of morphological traits, and leaf element content. Genome-wide associations revealed regions associated with variation and plasticity in multiple traits, including two regions with seemingly additive effects on NUtE change. Our results demonstrate potential avenues for improving sunflower nutrient stress tolerance independent of vigor, and highlight specific traits and genomic regions that could play a role in enhancing tolerance.

QTL mapping and candidate gene analysis reveal two major loci regulating green leaf color in non-heading Chinese cabbage

KEY MESSAGE

Two major-effect QTL GlcA07.1 and GlcA09.1 for green leaf color were fine mapped into 170.25 kb and 191.41 kb intervals on chromosomes A07 and A09, respectively, and were validated by transcriptome analysis. Non-heading Chinese cabbage (NHCC) is a leafy vegetable with a wide range of green colors. Understanding the genetic mechanism behind broad spectrum of green may facilitate the breeding of high-quality NHCC. Here, we used F2 and F7:8 recombination inbred line (RIL) population from a cross between Wutacai (dark-green) and Erqing (lime-green) to undertake the genetic analysis and quantitative trait locus (QTL) mapping in NHCC. The genetic investigation of the F2 population revealed that the variation of green leaf color was controlled by two recessive genes. Six pigments associated with green leaf color, including total chlorophyll, chlorophyll a, chlorophyll b, total carotenoids, lutein, and carotene were quantified and applied for QTL mapping in the RIL population. A total of 7 QTL were detected across the whole genome. Among them, two major-effect QTL were mapped on chromosomes A07 (GlcA07.1) and A09 (GlcA09.1) corresponding to two QTL identified in the F2 population. The QTL GlcA07.1 and GlcA09.1 were further fine mapped into 170.25 kb and 191.41 kb genomic regions, respectively. By comparing gene expression level and gene annotation, BraC07g023810 and BraC07g023970 were proposed as the best candidates for GlcA07.1, while BraC09g052220 and BraC09g052270 were suggested for GlcA09.1. Two InDel molecular markers (GlcA07.1-BcGUN4 and GlcA09.1-BcSG1) associated with BraC07gA023810 and BraC09g052220 were developed and could effectively identify leaf color in natural NHCC accessions, suggesting their potential for marker-assisted leaf color selection in NHCC breeding.

Control of false discoveries in grouped hypothesis testing for eQTL data

BACKGROUND

Expression quantitative trait locus (eQTL) analysis aims to detect the genetic variants that influence the expression of one or more genes. Gene-level eQTL testing forms a natural grouped-hypothesis testing strategy with clear biological importance. Methods to control family-wise error rate or false discovery rate for group testing have been proposed earlier, but may not be powerful or easily apply to eQTL data, for which certain structured alternatives may be defensible and may enable the researcher to avoid overly conservative approaches.

RESULTS

In an empirical Bayesian setting, we propose a new method to control the false discovery rate (FDR) for grouped hypotheses. Here, each gene forms a group, with SNPs annotated to the gene corresponding to individual hypotheses. The heterogeneity of effect sizes in different groups is considered by the introduction of a random effects component. Our method, entitled Random Effects model and testing procedure for Group-level FDR control (REG-FDR), assumes a model for alternative hypotheses for the eQTL data and controls the FDR by adaptive thresholding. As a convenient alternate approach, we also propose Z-REG-FDR, an approximate version of REG-FDR, that uses only Z-statistics of association between genotype and expression for each gene-SNP pair. The performance of Z-REG-FDR is evaluated using both simulated and real data. Simulations demonstrate that Z-REG-FDR performs similarly to REG-FDR, but with much improved computational speed.

CONCLUSION

Our results demonstrate that the Z-REG-FDR method performs favorably compared to other methods in terms of statistical power and control of FDR. It can be of great practical use for grouped hypothesis testing for eQTL analysis or similar problems in statistical genomics due to its fast computation and ability to be fit using only summary data.

Genomic wide association study and selective sweep analysis identify genes associated with improved yield under drought in Turkish winter wheat germplasm

A panel comprising of 84 Turkish winter wheat landraces (LR) and 73 modern varieties (MV) was analyzed with genome wide association study (GWAS) to identify genes/genomic regions associated with increased yield under favorable and drought conditions. In addition, selective sweep analysis was conducted to detect signatures of selection in the winter wheat genome driving the differentiation between LR and MV, to gather an understanding of genomic regions linked to adaptation and yield improvement. The panel was genotyped with 25 K wheat SNP array and phenotyped for agronomic traits for two growing seasons (2018 and 2019) in Konya, Turkey. Year 2018 was treated as drought environment due to very low precipitation prior to heading whereas year 2019 was considered as a favorable season. GWAS conducted with SNPs and haplotype blocks using mixed linear model identified 18 genomic regions in the vicinities of known genes i.e., TaERF3-3A, TaERF3-3B, DEP1-5A, FRIZZY PANICLE-2D, TaSnRK23-1A, TaAGL6-A, TaARF12-2A, TaARF12-2B, WAPO1, TaSPL16-7D, TaTGW6-A1, KAT-2B, TaOGT1, TaSPL21-6B, TaSBEIb, trs1/WFZP-A, TaCwi-A1-2A and TaPIN1-7A associated with grain yield (GY) and yield related traits. Haplotype-based GWAS identified five haplotype blocks (H1A-42, H2A-71, H4A-48, H7B-123 and H7B-124), with the favorable haplotypes showing a yield increase of > 700 kg/ha in the drought season. SNP-based GWAS, detected only one larger effect genomic region on chromosome 7B, in common with haplotype-based GWAS. On an average, the percentage variation (PV) explained by haplotypes was 8.0% higher than PV explained by SNPs for all the investigated traits. Selective sweep analysis detected 39 signatures of selection between LR and MV of which 15 were within proximity of known functional genes controlling flowering (PRR-A1, PPR-D1, TaHd1-6B), GY and GY components (TaSus2-2B, TaGS2-B1, AG1-1A/WAG1-1A, DUO-A1, DUO-B1, AG2-3A/WAG2-3A, TaLAX1, TaSnRK210-4A, FBP, TaLAX1, TaPIL1 and AP3-1-7A/WPA3-7A) and 10 regions underlying various transcription factors and regulatory genes. The study outcomes contribute to utilization of LR in breeding winter wheat.

Through the ages: a commentary on 'Diversification of quantitative morphological traits in wheat'

This article comments on:

Yixiang Shan and Colin P. Osborne, Diversification of quantitative morphological traits in wheat, Annals of Botany, Volume 133, Issue 3, 1 March 2024, Pages 413–426 https://doi.org/10.1093/aob/mcad202

Genetic linkage mapping and quantitative trait locus (QTL) analysis of sweet basil (Ocimum basilicum L.) to identify genomic regions associated with cold tolerance and major volatiles

Chilling sensitivity is one of the greatest challenges affecting the marketability and profitability of sweet basil (Ocimum basilicum L.) in the US and worldwide. Currently, there are no sweet basils commercially available with significant chilling tolerance and traditional aroma profiles. This study was conducted to identify quantitative trait loci (QTLs) responsible for chilling tolerance and aroma compounds in a biparental mapping population, including the Rutgers advanced breeding line that served as a chilling tolerant parent, 'CB15', the chilling sensitive parent, 'Rutgers Obsession DMR' and 200 F2 individuals. Chilling tolerance was assessed by percent necrosis using machine learning and aroma profiling was evaluated using gas chromatography (GC) mass spectrometry (MS). Single nucleotide polymorphism (SNP) markers were generated from genomic sequences derived from double digestion restriction-site associated DNA sequencing (ddRADseq) and converted to genotype data using a reference genome alignment. A genetic linkage map was constructed and five statistically significant QTLs were identified in response to chilling temperatures with possible interactions between QTLs. The QTL on LG24 (qCH24) demonstrated the largest effect for chilling response and was significant in all three replicates. No QTLs were identified for linalool, as the population did not segregate sufficiently to detect this trait. Two significant QTLs were identified for estragole (also known as methyl chavicol) with only qEST1 on LG1 being significant in the multiple-QTL model (MQM). QEUC26 was identified as a significant QTL for eucalyptol (also known as 1,8-cineole) on LG26. These QTLs may represent key mechanisms for chilling tolerance and aroma in basil, providing critical knowledge for future investigation of these phenotypic traits and molecular breeding.

The story of a decade: Genomics, functional genomics, and molecular breeding in Brassica napus

Rapeseed (Brassica napus L.) is one of the major global sources of edible vegetable oil and is also used as a feed and pioneer crop and for sightseeing and industrial purposes. Improvements in genome sequencing and molecular marker technology have fueled a boom in functional genomic studies of major agronomic characters such as yield, quality, flowering time, and stress resistance. Moreover, introgression and pyramiding of key functional genes have greatly accelerated the genetic improvement of important traits. Here we summarize recent progress in rapeseed genomics and genetics, and we discuss effective molecular breeding strategies by exploring these findings in rapeseed. These insights will extend our understanding of the molecular mechanisms and regulatory networks underlying agronomic traits and facilitate the breeding process, ultimately contributing to more sustainable agriculture throughout the world.

Genetic dissection of ten photosynthesis-related traits based on InDel- and SNP-GWAS in soybean

KEY MESSAGE

A total of 416 InDels and 112 SNPs were significantly associated with soybean photosynthesis-related traits. GmIWS1 and GmCDC48 might be related to chlorophyll fluorescence and gas-exchange parameters, respectively. Photosynthesis is one of the main factors determining crop yield. A better understanding of the genetic architecture for photosynthesis is of great significance for soybean yield improvement. Our previous studies identified 5,410,112 single nucleotide polymorphisms (SNPs) from the resequencing data of 219 natural soybean accessions. Here, we identified 634,106 insertions and deletions (InDels) from these 219 accessions and used these InDel variations to perform principal component and linkage disequilibrium analysis of this population. The genome-wide association study (GWAS) were conducted on six chlorophyll fluorescence parameters (chlorophyll content, light energy absorbed per reaction center, quantum yield for electron transport, probability that a trapped exciton moves an electron into the electron transport chain beyond primary quinone acceptor, maximum quantum yield of photosystem II primary photochemistry in the dark-adapted state, performance index on absorption basis) and four gas-exchange parameters (intercellular carbon dioxide concentration, stomatal conductance, net photosynthesis rate, transpiration rate) and revealed 416 significant InDels and 112 significant SNPs. Based on GWAS results, GmIWS1 (encoding a transcription elongation factor) and GmCDC48 (encoding a cell division cycle protein) with the highest expression in the mapping region were determined as the candidate genes responsible for chlorophyll fluorescence and gas-exchange parameters, respectively. Further identification of favorable haplotypes with higher photosynthesis, seed weight and seed yield were carried out for GmIWS1 and GmCDC48. Overall, this study revealed the natural variations and candidate genes underlying the photosynthesis-related traits based on abundant phenotypic and genetic data, providing valuable insights into the genetic mechanisms controlling photosynthesis and yield in soybean.

GWAS supported by computer vision identifies large numbers of candidate regulators of in planta regeneration in Populus trichocarpa

Plant regeneration is an important dimension of plant propagation and a key step in the production of transgenic plants. However, regeneration capacity varies widely among genotypes and species, the molecular basis of which is largely unknown. Association mapping methods such as genome-wide association studies (GWAS) have long demonstrated abilities to help uncover the genetic basis of trait variation in plants; however, the performance of these methods depends on the accuracy and scale of phenotyping. To enable a large-scale GWAS of in planta callus and shoot regeneration in the model tree Populus, we developed a phenomics workflow involving semantic segmentation to quantify regenerating plant tissues over time. We found that the resulting statistics were of highly non-normal distributions, and thus employed transformations or permutations to avoid violating assumptions of linear models used in GWAS. We report over 200 statistically supported quantitative trait loci (QTLs), with genes encompassing or near to top QTLs including regulators of cell adhesion, stress signaling, and hormone signaling pathways, as well as other diverse functions. Our results encourage models of hormonal signaling during plant regeneration to consider keystone roles of stress-related signaling (e.g. involving jasmonates and salicylic acid), in addition to the auxin and cytokinin pathways commonly considered. The putative regulatory genes and biological processes we identified provide new insights into the biological complexity of plant regeneration, and may serve as new reagents for improving regeneration and transformation of recalcitrant genotypes and species.

A k-mer-based bulked segregant analysis approach to map seed traits in unphased heterozygous potato genomes

High-throughput sequencing-based methods for bulked segregant analysis (BSA) allow for the rapid identification of genetic markers associated with traits of interest. BSA studies have successfully identified qualitative (binary) and quantitative trait loci (QTLs) using QTL mapping. However, most require population structures that fit the models available and a reference genome. Instead, high-throughput short-read sequencing can be combined with BSA of k-mers (BSA-k-mer) to map traits that appear refractory to standard approaches. This method can be applied to any organism and is particularly useful for species with genomes diverged from the closest sequenced genome. It is also instrumental when dealing with highly heterozygous and potentially polyploid genomes without phased haplotype assemblies and for which a single haplotype can control a trait. Finally, it is flexible in terms of population structure. Here, we apply the BSA-k-mer method for the rapid identification of candidate regions related to seed spot and seed size in diploid potato. Using a mixture of F1 and F2 individuals from a cross between 2 highly heterozygous parents, candidate sequences were identified for each trait using the BSA-k-mer approach. Using parental reads, we were able to determine the parental origin of the loci. Finally, we mapped the identified k-mers to a closely related potato genome to validate the method and determine the genomic loci underlying these sequences. The location identified for the seed spot matches with previously identified loci associated with pigmentation in potato. The loci associated with seed size are novel. Both loci are relevant in future breeding toward true seeds in potato.

Detecting genetic effects on phenotype variability to capture gene-by-environment interactions: a systematic method comparison

Genetically associated phenotypic variability has been widely observed across organisms and traits, including in humans. Both gene-gene and gene-environment interactions can lead to an increase in genetically associated phenotypic variability. Therefore, detecting the underlying genetic variants, or variance Quantitative Trait Loci (vQTLs), can provide novel insights into complex traits. Established approaches to detect vQTLs apply different methodologies from variance-only approaches to mean-variance joint tests, but a comprehensive comparison of these methods is lacking. Here, we review available methods to detect vQTLs in humans, carry out a simulation study to assess their performance under different biological scenarios of gene-environment interactions, and apply the optimal approaches for vQTL identification to gene expression data. Overall, with a minor allele frequency (MAF) of less than 0.2, the squared residual value linear model (SVLM) and the deviation regression model (DRM) are optimal when the data follow normal and non-normal distributions, respectively. In addition, the Brown-Forsythe (BF) test is one of the optimal methods when the MAF is 0.2 or larger, irrespective of phenotype distribution. Additionally, a larger sample size and more balanced sample distribution in different exposure categories increase the power of BF, SVLM, and DRM. Our results highlight vQTL detection methods that perform optimally under realistic simulation settings and show that their relative performance depends on the phenotype distribution, allele frequency, sample size, and the type of exposure in the interaction model underlying the vQTL.

Characterization of adaptation mechanisms in sorghum using a multireference back-cross nested association mapping design and envirotyping

Identifying the genetic factors impacting the adaptation of crops to environmental conditions is of key interest for conservation and selection purposes. It can be achieved using population genomics, and evolutionary or quantitative genetics. Here we present a sorghum multireference back-cross nested association mapping population composed of 3,901 lines produced by crossing 24 diverse parents to 3 elite parents from West and Central Africa-back-cross nested association mapping. The population was phenotyped in environments characterized by differences in photoperiod, rainfall pattern, temperature levels, and soil fertility. To integrate the multiparental and multi-environmental dimension of our data we proposed a new approach for quantitative trait loci (QTL) detection and parental effect estimation. We extended our model to estimate QTL effect sensitivity to environmental covariates, which facilitated the integration of envirotyping data. Our models allowed spatial projections of the QTL effects in agro-ecologies of interest. We utilized this strategy to analyze the genetic architecture of flowering time and plant height, which represents key adaptation mechanisms in environments like West Africa. Our results allowed a better characterization of well-known genomic regions influencing flowering time concerning their response to photoperiod with Ma6 and Ma1 being photoperiod-sensitive and the region of possible candidate gene Elf3 being photoperiod-insensitive. We also accessed a better understanding of plant height genetic determinism with the combined effects of phenology-dependent (Ma6) and independent (qHT7.1 and Dw3) genomic regions. Therefore, we argue that the West and Central Africa-back-cross nested association mapping and the presented analytical approach constitute unique resources to better understand adaptation in sorghum with direct application to develop climate-smart varieties.

Linkage mapping of root shape traits in two carrot populations

This study investigated the genetic basis of carrot root shape traits using composite interval mapping in two biparental populations (n = 119 and n = 128). The roots of carrot F2:3 progenies were grown over 2 years and analyzed using a digital imaging pipeline to extract root phenotypes that compose market class. Broad-sense heritability on an entry-mean basis ranged from 0.46 to 0.80 for root traits. Reproducible quantitative trait loci (QTL) were identified on chromosomes 2 and 6 on both populations. Colocalization of QTLs for phenotypically correlated root traits was also observed and coincided with previously identified QTLs in published association and linkage mapping studies. Individual QTLs explained between 14 and 27% of total phenotypic variance across traits, while four QTLs for length-to-width ratio collectively accounted for up to 73% of variation. Predicted genes associated with the OFP-TRM (OVATE Family Proteins-TONNEAU1 Recruiting Motif) and IQD (IQ67 domain) pathway were identified within QTL support intervals. This observation raises the possibility of extending the current regulon model of fruit shape to include carrot storage roots. Nevertheless, the precise molecular mechanisms through which this pathway operates in roots characterized by secondary growth originating from cambium layers remain unknown.

Genetic improvement of Egyptian cotton (Gossypium barbadense L.) for high yield and fiber quality properties under semi arid conditions

Between 2016 and 2018, the Agriculture Research Center's Sakha Agriculture Research Station conducted two rounds of pedigree selection on a segregating population of cotton (Gossypium barbadense L.) using the F2, F3, and F4 generations resulting from crossing Giza 94 and Suvin. In 2016, the top 5% of plants from the F2 population were selected based on specific criteria. The superior families from the F3 generation were then selected to produce the F4 families in 2017, which were grown in the 2018 summer season in single plant progeny rows and bulk experiments with a randomized complete block design of three replications. Over time, most traits showed increased mean values in the population, with the F2 generation having higher Genotypic Coefficient of Variance (GCV) and Phenotypic Coefficient of Variance (PCV) values compared to the succeeding generations for the studied traits. The magnitude of GCV and PCV in the F3 and F4 generations was similar, indicating that genotype had played a greater role than the environment. Moreover, the mean values of heritability in the broad sense increased from generation to generation. Selection criteria I2, I4, and I5 were effective in improving most of the yield and its component traits, while selection criterion I1 was efficient in improving earliness traits. Most of the yield and its component traits showed a positive and significant correlation with each other, highlighting their importance in cotton yield. This suggests that selecting to improveone or more of these traits would improve the others. Families number 9, 13, 19, 20, and 21 were the best genotypes for relevant yield characters, surpassing the better parent, check variety, and giving the best values for most characters. Therefore, the breeder could continue to use these families in further generations as breeding genotypes to develop varieties with high yields and its components.

Mapping QTL controlling count traits with excess zeros and ones using a zero-and-one-inflated generalized Poisson regression model

The research on the quantitative trait locus (QTL) mapping of count data has aroused the wide attention of researchers. There are frequent problems in applied research that limit the application of the conventional Poisson model in the analysis of count phenotypes, which include the overdispersion and excess zeros and ones. In this article, a novel model, that is, the zero-and-one-inflated generalized Poisson (ZOIGP) model, is proposed to deal with these problems. Based on the proposed model, a score test is performed for the inflation parameter, in which the ZOIGP model with a constant proportion of excess zeros and ones is compared with a standard generalized Poisson model. To illustrate the practicability of the ZOIGP model, we extend it to the QTL interval mapping application that underpins count phenotype with excess zeros and excess ones. The genetic effects are estimated utilizing the expectation-maximization algorithm embedded with the Newton-Raphson algorithm, and the genome-wide scan and likelihood ratio test is performed to map and test the potential QTLs. The statistical properties exhibited by the proposed method are investigated through simulation. Finally, a real data analysis example is used to illustrate the utility of the proposed method for QTL mapping.

A novel OsHB5-OsAPL-OsMADS27/OsWRKY102 regulatory module regulates grain size in rice

Grain size, a crucial trait that determines rice yield and quality, is typically regulated by multiple genes. Although numerous genes controlling grain size have been identified, the precise and dynamic regulatory network governing grain size is still not fully understood. In this study, we unveiled a novel regulatory module composed of OsHB5, OsAPL and OsMADS27/OsWRKY102, which plays a crucial role in modulating grain size in rice. As a positive regulator of grain size, OsAPL has been found to interact with OsHB5 both in vitro and in vivo. Through chromatin immunoprecipitation-sequencing, we successfully mapped two potential targets of OsAPL, namely OsMADS27, a positive regulator in grain size and OsWRKY102, a negative regulator in lignification that is also associated with grain size control. Further evidence from EMSA and chromatin immunoprecipitation-quantitative PCR experiments has shown that OsAPL acts as an upstream transcription factor that directly binds to the promoters of OsMADS27 and OsWRKY102. Moreover, EMSA and dual-luciferase reporter assays have indicated that the interaction between OsAPL and OsHB5 enhances the repressive effect of OsAPL on OsMADS27 and OsWRKY102. Collectively, our findings discovered a novel regulatory module, OsHB5-OsAPL-OsMADS27/OsWRKY102, which plays a significant role in controlling grain size in rice. These discoveries provide potential targets for breeding high-yield and high-quality rice varieties.

Population genetic analysis based on the polymorphisms mediated by transposons in the genomes of pig

Transposable elements (TEs) mobility is capable of generating a large number of structural variants (SVs), which can have considerable potential as molecular markers for genetic analysis and molecular breeding in livestock. Our results showed that the pig genome contains mainly TE-SVs generated by short interspersed nuclear elements (51,873/76.49%), followed by long interspersed nuclear elements (11,131/16.41%), and more than 84% of the common TE-SVs (Minor allele frequency, MAF > 0.10) were validated to be polymorphic. Subsequently, we utilized the identified TE-SVs to gain insights into the population structure, resulting in clear differentiation among the three pig groups and facilitating the identification of relationships within Chinese local pig breeds. In addition, we investigated the frequencies of TEs in the gene coding regions of different pig groups and annotated the respective TE types, related genes, and functional pathways. Through genome-wide comparisons of Large White pigs and Chinese local pigs utilizing the Beijing Black pigs, we identified TE-mediated SVs associated with quantitative trait loci and observed that they were mainly involved in carcass traits and meat quality traits. Lastly, we present the first documented evidence of TE transduction in the pig genome.

Evaluating three strategies of genome-wide association analysis for integrating data from multiple populations

In livestock, genome-wide association studies (GWAS) are usually conducted in a single population (single-GWAS) with limited sample size and detection power. To enhance the detection power of GWAS, meta-analysis of GWAS (meta-GWAS) and mega-analysis of GWAS (mega-GWAS) have been proposed to integrate data from multiple populations at the level of summary statistics or individual data, respectively. However, there is a lack of comparison for these different strategies, which makes it difficult to guide the best practice of GWAS integrating data from multiple study populations. To maximize the comparison of different association analysis strategies across multiple populations, we conducted single-GWAS, meta-GWAS, and mega-GWAS for the backfat thickness of 100 kg (BFT_100) and days to 100 kg (DAYS_100) within each of the three commercial pig breeds (Duroc, Yorkshire, and Landrace). Based on controlling the genome inflation factor to one, we calculated corrected p-values (pC ). In Yorkshire, with the largest sample size, mega-GWAS, meta-GWAS and single-GWAS detected 149, 38 and 20 significant SNPs (pC  < 1E-5) associated with BFT_100, as well as 26, four, and one QTL, respectively. Among them, pC of SNPs from mega-GWAS was the lowest, followed by meta-GWAS and single-GWAS. The correlation of pC among the three GWAS strategies ranged from 0.60 to 0.75 and the correlation of SNP effect values between meta-GWAS and mega-GWAS was 0.74, all showing good agreement. Collectively, even though there are differences in the integration of individual data or summary statistics, integrating data from multiple populations is an effective means of genetic argument for complex traits, especially mega-GWAS versus single-GWAS.

Genome-wide association analysis revealed new QTL and candidate genes affecting the teat number in Dutch Large White pigs

Teat number (TNUM) is an important reproductive trait of sows, which affects the weaning survival rate of piglets. In this study, 1166 Dutch Large White pigs with TNUM phenotype were used as the research object. These pigs were genotyped by 50K SNP chip and the chip data were further imputed to the resequencing level. The estimated heritabilities of left teat number (LTN), right teat number (RTN) and total teat number (TTN) were 0.21, 0.19 and 0.3, respectively. Based on chip data, significant SNPs for RTN on SSC2, SSC5, SSC9 and SSC13 were identified using genome-wide association analysis (GWAS). Significant SNPs for TTN were identified on SSC2, SSC5 and SSC7. Based on imputed data, the GWAS identified a significant SNP (rs329158522) for LTN on SSC17, two significant SNPs (rs342855242 and rs80813115) for RTN on SSC2 and SSC9, and two significant SNPs (rs327003548 and rs326943811) for TTN on SSC5 and SSC6. Among them, four novel QTL were discovered. The Bayesian fine-mapping method was used to fine map the QTL identified in the GWAS of the imputed data, and the confidence intervals of QTL affecting LTN (SSC17: 45.22-46.20 Mb), RTN (SSC9: 122.18-122.80 Mb) and TTN (SSC5: 14.01-15.91 Mb, SSC6: 120.06-121.25 Mb) were detected. A total of 52 candidate genes were obtained. Furthermore, we identified five candidate genes, WNT10B, AQP5, FMNL3, NUAK1 and CKAP4, for the first time, which involved in breast development and other related functions by gene annotation. Overall, this study provides new molecular markers for the breeding of teat number in pigs.

Genome-wide association study of growth curve parameters reveals novel genomic regions and candidate genes associated with metatarsal bone traits in chickens

The growth and development of chicken bones have an enormous impact on the health and production performance of chickens. However, the development pattern and genetic regulation of the chicken skeleton are poorly understood. This study aimed to evaluate metatarsal bone growth and development patterns in chickens via non-linear models, and to identify the genetic determinants of metatarsal bone traits using a genome-wide association study (GWAS) based on growth curve parameters. Data on metatarsal length (MeL) and metatarsal circumference (MeC) were obtained from 471 F2 chickens (generated by crossing broiler sires, derived from a line selected for high abdominal fat, with Baier layer dams) at 4, 6, 8, 10, and 12 weeks of age. Four non-linear models (Gompertz, Logistic, von Bertalanffy, and Brody) were used to fit the MeL and MeC growth curves. Subsequently, the estimated growth curve parameters of the mature MeL or MeC (A), time-scale parameter (b), and maturity rate (K) from the non-linear models were utilized as substitutes for the original bone data in GWAS. The Logistic and Brody models displayed the best goodness-of-fit for MeL and MeC, respectively. Single-trait and multi-trait GWASs based on the growth curve parameters of the Logistic and Brody models revealed 4 618 significant single nucleotide polymorphisms (SNPs), annotated to 332 genes, associated with metatarsal bone traits. The majority of these significant SNPs were located on Gallus gallus chromosome (GGA) 1 (167.433-176.318 Mb), GGA2 (96.791-103.543 Mb), GGA4 (65.003-83.104 Mb) and GGA6 (64.685-95.285 Mb). Notably, we identified 12 novel GWAS loci associated with chicken metatarsal bone traits, encompassing 35 candidate genes. In summary, the combination of single-trait and multi-trait GWASs based on growth curve parameters uncovered numerous genomic regions and candidate genes associated with chicken bone traits. The findings benefit an in-depth understanding of the genetic architecture underlying metatarsal growth and development in chickens.

Mapping and validation of Fusarium wilt race 2 resistance QTL from Citrullus amarus line USVL246-FR2

KEY MESSAGE

Fon race 2 resistant QTLs were identified on chromosomes 8 and 9. Families homozygous for resistance alleles at a haplotype of three KASP markers had 42% lower disease severity than those with susceptible alleles in an independent, interspecific validation population confirming their utility for introgression of Fusarium wilt resistance. Fusarium oxysporum f. sp. niveum (Fon) race 2 causes Fusarium wilt in watermelon and threatens watermelon production worldwide. Chemical management options are not effective, and no resistant edible watermelon cultivars have been released. Implementation of marker-assisted selection to develop resistant cultivars requires identifying sources of resistance and the underlying quantitative trait loci (QTL), developing molecular markers associated with the QTL, and validating marker-phenotype associations with an independent population. An intraspecific Citrullus amarus recombinant inbred line population from a cross of resistant USVL246-FR2 and susceptible USVL114 was used for mapping Fon race 2 resistance QTL. KASP markers were developed (N = 51) for the major QTL on chromosome 9 and minor QTL on chromosomes 1, 6, and 8. An interspecific F2:3 population was developed from resistance donor USVL246-FR2 (C. amarus) and a susceptible cultivar 'Sugar Baby' (Citrullus lanatus) to validate the utility of the markers for introgression of resistance from the wild crop relative into cultivated watermelon. Only 16 KASP markers segregated in the interspecific C. amarus/lanatus validation population. Four markers showed significant differences in the separation of genotypes based on family-mean disease severity, but together explained only 16% of the phenotypic variance. Genotypes that inherited homozygous resistant parental alleles at three KASP markers had 42% lower family-mean disease severity than homozygous susceptible genotypes. Thus, haplotype analysis was more effective at predicting the mean disease severity of families than single markers. The haplotype identified in this study will be valuable for developing Fon race 2 resistant watermelon cultivars.

Omics-based construction of regulatory variants can be applied to help decipher pig liver-related traits

Genetic variants can influence complex traits by altering gene expression through changes to regulatory elements. However, the genetic variants that affect the activity of regulatory elements in pigs are largely unknown, and the extent to which these variants influence gene expression and contribute to the understanding of complex phenotypes remains unclear. Here, we annotate 90,991 high-quality regulatory elements using acetylation of histone H3 on lysine 27 (H3K27ac) ChIP-seq of 292 pig livers. Combined with genome resequencing and RNA-seq data, we identify 28,425 H3K27ac quantitative trait loci (acQTLs) and 12,250 expression quantitative trait loci (eQTLs). Through the allelic imbalance analysis, we validate two causative acQTL variants in independent datasets. We observe substantial sharing of genetic controls between gene expression and H3K27ac, particularly within promoters. We infer that 46% of H3K27ac exhibit a concomitant rather than causative relationship with gene expression. By integrating GWAS, eQTLs, acQTLs, and transcription factor binding prediction, we further demonstrate their application, through metabolites dulcitol, phosphatidylcholine (PC) (16:0/16:0) and published phenotypes, in identifying likely causal variants and genes, and discovering sub-threshold GWAS loci. We provide insight into the relationship between regulatory elements and gene expression, and the genetic foundation for dissecting the molecular mechanism of phenotypes.

TaCRTISO dosage modulates plant height and spike number per plant in wheat

An allelic variation of TaCRTISO is valuable in adjusting spike number per plant and plant height in wheat breeding.

Identification of QTLs for symbiotic nitrogen fixation and related traits in a soybean recombinant inbred line population

KEY MESSAGE

The genetic architecture of symbiotic N fixation and related traits was investigated in the field. QTLs were identified for percent N derived from the atmosphere, shoot [N] and C to N ratio. Soybean [Glycine max (L.) Merr.] is cultivated worldwide and is the most abundant source of plant-based protein. Symbiotic N2 fixation (SNF) in legumes such as soybean is of great importance; however, yields may still be limited by N in both high yielding and stressful environments. To better understand the genetic architecture of SNF and facilitate the development of high yielding cultivars and sustainable soybean production in stressful environments, a recombinant inbred line population consisting of 190 lines, developed from a cross between PI 442012A and PI 404199, was evaluated for N derived from the atmosphere (Ndfa), N concentration ([N]), and C to N ratio (C/N) in three environments. Significant genotype, environment and genotype × environment effects were observed for all three traits. A linkage map was constructed containing 3309 single nucleotide polymorphism (SNP) markers. QTL analysis was performed for additive effects of QTLs, QTL × environment interactions, and QTL × QTL interactions. Ten unique additive QTLs were identified across all traits and environments. Of these, two QTLs were detected for Ndfa and eight for C/N. Of the eight QTLs for C/N, four were also detected for [N]. Using QTL × environment analysis, six QTLs were detected, of which five were also identified in the additive QTL analysis. The QTL × QTL analysis identified four unique epistatic interactions. The results of this study may be used for genomic selection and introgression of favorable alleles for increased SNF, [N], and C/N via marker-assisted selection.

Mapping of a Major-Effect Quantitative Trait Locus for Seed Dormancy in Wheat

The excavation and utilization of dormancy loci in breeding are effective endeavors for enhancing the resistance to pre-harvest sprouting (PHS) of wheat varieties. CH1539 is a wheat breeding line with high-level seed dormancy. To clarify the dormant loci carried by CH1539 and obtain linked molecular markers, in this study, a recombinant inbred line (RIL) population derived from the cross of weak dormant SY95-71 and strong dormant CH1539 was genotyped using the Wheat17K single-nucleotide polymorphism (SNP) array, and a high-density genetic map covering 21 chromosomes and consisting of 2437 SNP markers was constructed. Then, the germination percentage (GP) and germination index (GI) of the seeds from each RIL were estimated. Two QTLs for GP on chromosomes 5A and 6B, and four QTLs for GI on chromosomes 5A, 6B, 6D and 7A were identified. Among them, the QTL on chromosomes 6B controlling both GP and GI, temporarily named QGp/Gi.sxau-6B, is a major QTL for seed dormancy with the maximum phenotypic variance explained of 17.66~34.11%. One PCR-based diagnostic marker Ger6B-3 for QGp/Gi.sxau-6B was developed, and the genetic effect of QGp/Gi.sxau-6B on the RIL population and a set of wheat germplasm comprising 97 accessions was successfully confirmed. QGp/Gi.sxau-6B located in the 28.7~30.9 Mbp physical position is different from all the known dormancy loci on chromosomes 6B, and within the interval, there are 30 high-confidence annotated genes. Our results revealed a novel QTL QGp/Gi.sxau-6B whose CH1539 allele had a strong and broad effect on seed dormancy, which will be useful in further PHS-resistant wheat breeding.

Quantitative trait loci detection for three tiller-related traits and the effects on wheat (Triticum aestivum L.) yields

KEY MESSAGE

A total of 38 putative additive QTLs and 55 pairwise putative epistatic QTLs for tiller-related traits were reported, and the candidate genes underlying qMtn-KJ-5D, a novel major and stable QTL for maximum tiller number, were characterized. Tiller-related traits play an important role in determining the yield potential of wheat. Therefore, it is important to elucidate the genetic basis for tiller number when attempting to use genetic improvement as a tool for enhancing wheat yields. In this study, a quantitative trait locus (QTL) analysis of three tiller-related traits was performed on the recombinant inbred lines (RILs) of a mapping population, referred to as KJ-RILs, that was derived from a cross between the Kenong 9204 (KN9204) and Jing 411 (J411) lines. A total of 38 putative additive QTLs and 55 pairwise putative epistatic QTLs for spike number per plant (SNPP), maximum tiller number (MTN), and ear-bearing tiller rate (EBTR) were detected in eight different environments. Among these QTLs with additive effects, three major and stable QTLs were first documented herein. Almost all but two pairwise epistatic QTLs showed minor interaction effects accounting for no more than 3.0% of the phenotypic variance. The genetic effects of two colocated major and stable QTLs, i.e., qSnpp-KJ-5D.1 and qMtn-KJ-5D, for yield-related traits were characterized. The breeding selection effect of the beneficial allele for the two QTLs was characterized, and its genetic effects on yield-related traits were evaluated. The candidate genes underlying qMtn-KJ-5D were predicted based on multi-omics data, and TraesKN5D01HG00080 was identified as a likely candidate gene. Overall, our results will help elucidate the genetic architecture of tiller-related traits and can be used to develop novel wheat varieties with high yields.

Elucidation of the genetic determination of body weight and size in Chinese local chicken breeds by large-scale genomic analyses

BACKGROUND

Body weight and size are important economic traits in chickens. While many growth-related quantitative trait loci (QTLs) and candidate genes have been identified, further research is needed to confirm and characterize these findings. In this study, we investigate genetic and genomic markers associated with chicken body weight and size. This study provides new insights into potential markers for genomic selection and breeding strategies to improve meat production in chickens.

METHODS

We performed whole-genome resequencing of and Wenshang Barred (WB) chickens (n = 596) and three additional breeds with varying body sizes (Recessive White (RW), WB, and Luxi Mini (LM) chickens; (n = 50)). We then used selective sweeps of mutations coupled with genome-wide association study (GWAS) to identify genomic markers associated with body weight and size.

RESULTS

We identified over 9.4 million high-quality single nucleotide polymorphisms (SNPs) among three chicken breeds/lines. Among these breeds, 287 protein-coding genes exhibited positive selection in the RW and WB populations, while 241 protein-coding genes showed positive selection in the LM and WB populations. Genomic heritability estimates were calculated for 26 body weight and size traits, including body weight, chest breadth, chest depth, thoracic horn, body oblique length, keel length, pelvic width, shank length, and shank circumference in the WB breed. The estimates ranged from 0.04 to 0.67. Our analysis also identified a total of 2,522 genome-wide significant SNPs, with 2,474 SNPs clustered around two genomic regions. The first region, located on chromosome 4 (7.41-7.64 Mb), was linked to body weight after ten weeks and body size traits. LCORL, LDB2, and PPARGC1A were identified as candidate genes in this region. The other region, located on chromosome 1 (170.46-171.53 Mb), was associated with body weight from four to eighteen weeks and body size traits. This region contained CAB39L and WDFY2 as candidate genes. Notably, LCORL, LDB2, and PPARGC1A showed highly selective signatures among the three breeds of chicken with varying body sizes.

CONCLUSION

Overall this study provides a comprehensive map of genomic variants associated with body weight and size in chickens. We propose two genomic regions, one on chromosome 1 and the other on chromosome 4, that could helpful for developing genome selection breeding strategies to enhance meat yield in chickens.

Pangenome-genotyped structural variation improves molecular phenotype mapping in cattle

Expression and splicing quantitative trait loci (e/sQTL) are large contributors to phenotypic variability. Achieving sufficient statistical power for e/sQTL mapping requires large cohorts with both genotypes and molecular phenotypes, and so, the genomic variation is often called from short-read alignments, which are unable to comprehensively resolve structural variation. Here we build a pangenome from 16 HiFi haplotype-resolved cattle assemblies to identify small and structural variation and genotype them with PanGenie in 307 short-read samples. We find high (>90%) concordance of PanGenie-genotyped and DeepVariant-called small variation and confidently genotype close to 21 million small and 43,000 structural variants in the larger population. We validate 85% of these structural variants (with MAF > 0.1) directly with a subset of 25 short-read samples that also have medium coverage HiFi reads. We then conduct e/sQTL mapping with this comprehensive variant set in a subset of 117 cattle that have testis transcriptome data, and find 92 structural variants as causal candidates for eQTL and 73 for sQTL. We find that roughly half of the top associated structural variants affecting expression or splicing are transposable elements, such as SV-eQTL for STN1 and MYH7 and SV-sQTL for CEP89 and ASAH2 Extensive linkage disequilibrium between small and structural variation results in only 28 additional eQTL and 17 sQTL discovered when including SVs, although many top associated SVs are compelling candidates.

Unveiling synergistic QTLs associated with slow wilting in soybean (Glycine max [L.] Merr.)

KEY MESSAGE

A stable QTL qSW_Gm10 works with a novel locus, qSW_Gm01, in a synergistic manner for controlling slow-wilting traits at the early vegetative stage under drought stress in soybean. Drought is one of the major environmental factors which limits soybean yield. Slow wilting is a promising trait that can enhance drought resilience in soybean without additional production costs. Recently, a Korean soybean cultivar SS2-2 was reported to exhibit slow wilting at the early vegetative stages. To find genetic loci responsible for slow wilting, in this study, quantitative trait loci (QTL) analysis was conducted using a recombinant inbred line (RIL) population derived from crossing between Taekwangkong (fast-wilting) and SS2-2 (slow-wilting). Wilting score and leaf moisture content were evaluated at the early vegetative stages for three years. Using the ICIM-MET module, a novel QTL on Chr01, qSW_Gm01 was identified, together with a previously known QTL, qSW_Gm10. These two QTLs were found to work synergistically for slow wilting of the RILs under the water-restricted condition. Furthermore, the SNP markers from the SoySNP50K dataset, located within these QTLs, were associated with the wilting phenotype in 30 diverse soybean accessions. Two genes encoding protein kinase 1b and multidrug resistance-associated protein 4 were proposed as candidate genes for qSW_Gm01 and qSW_Gm10, respectively, based on a comprehensive examination of sequence variation and gene expression differences in the parental lines under drought conditions. These genes may play a role in slow wilting by optimally regulating stomatal aperture. Our findings provide promising genetic resources for improving drought resilience in soybean and give valuable insights into the genetic mechanisms governing slow wilting.

Identification and validation of major and stable quantitative trait locus for falling number in common wheat (Triticum aestivum L.)

KEY MESSAGE

A major and stable QTL, QFn.sau-1B.2, which can explain 13.6% of the PVE in FN and has a positive effect on resistance in SGR, was mapped and validated. The falling number (FN) is considered one of the most important quality traits of wheat grain and is the most important quality evaluation index for wheat trade worldwide. The quantitative trait loci (QTLs) for FN were mapped in three years of experiments. 23, 30, and 58 QTLs were identified using the ICIM-BIP, ICIM-MET, and ICIM-EPI methods, respectively. Among them, seven QTLs were considered stable. QFn.sau-1B.2, which was mapped to the 1BL chromosome, can explain 13.6% of the phenotypic variation on average and is considered a major and stable QTL for FN. This QTL was mapped in a 1 cM interval and is flanked by the markers AX-110409346 and AX-108743901. Epistatic analysis indicated that QFN.sau-1B.2 has a strong influence on FN through both additive and epistatic effects. The Kompetitive Allele-Specific PCR marker KASP-AX-108743901, which is closely linked to QFn.sau-1B.2, was designed. The genetic effect of QFn.sau-1B.2 on FN was successfully confirmed in Chuannong18 × T1208 and CN17 × CN11 populations. Moreover, the results of the additive effects of favorable alleles for FN showed that the QTLs for FN had significant effects not only on FN but also on the resistance to spike germination. Within the interval of QFn.sau-1B.2, 147 high-confidence genes were found. According to the gene annotation and the transcriptome data, four genes might be associated with FN. QFn.sau-1B.2 may provide a new resource for the high-quality breeding of wheat in the future.

Identification of additional body weight QTLs in the Berlin Fat Mouse BFMI861 lines using time series data

The Berlin Fat Mouse Inbred line (BFMI) is a model for obesity and metabolic syndrome. The sublines BFMI861-S1 and BFMI861-S2 differ in weight despite high genetic similarity and a shared obesity-related locus. This study focused on identifying additional body weight quantitative trait loci (QTLs) by analyzing weekly weight measurements in a male population of the advanced intercross line BFMI861-S1 x BFMI861-S2. QTL analysis, utilizing 200 selectively genotyped mice (GigaMUGA) and 197 males genotyped for top SNPs, revealed a genome-wide significant QTL on Chr 15 (68.46 to 81.40 Mb) for body weight between weeks 9 to 20. Notably, this QTL disappeared (weeks 21 to 23) and reappeared (weeks 24 and 25) coinciding with a diet change. Additionally, a significant body weight QTL on Chr 16 (3.89 to 22.79 Mb) was identified from weeks 6 to 25. Candidate genes, including Gpt, Cbx6, Apol6, Apol8, Sun2 (Chr 15) and Trap1, Rrn3, Mapk1 (Chr 16), were prioritized. This study unveiled two additional body weight QTLs, one of which is novel and responsive to diet changes. These findings illuminate genomic regions influencing weight in BFMI and emphasize the utility of time series data in uncovering novel genetic factors.

Genomic dissection of the correlation between milk yield and various health traits using functional and evolutionary information about imputed sequence variants of 34,497 German Holstein cows

BACKGROUND

Over the last decades, it was subject of many studies to investigate the genomic connection of milk production and health traits in dairy cattle. Thereby, incorporating functional information in genomic analyses has been shown to improve the understanding of biological and molecular mechanisms shaping complex traits and the accuracies of genomic prediction, especially in small populations and across-breed settings. Still, little is known about the contribution of different functional and evolutionary genome partitioning subsets to milk production and dairy health. Thus, we performed a uni- and a bivariate analysis of milk yield (MY) and eight health traits using a set of ~34,497 German Holstein cows with 50K chip genotypes and ~17 million imputed sequence variants divided into 27 subsets depending on their functional and evolutionary annotation. In the bivariate analysis, eight trait-combinations were observed that contrasted MY with each health trait. Two genomic relationship matrices (GRM) were included, one consisting of the 50K chip variants and one consisting of each set of subset variants, to obtain subset heritabilities and genetic correlations. In addition, 50K chip heritabilities and genetic correlations were estimated applying merely the 50K GRM.

RESULTS

In general, 50K chip heritabilities were larger than the subset heritabilities. The largest heritabilities were found for MY, which was 0.4358 for the 50K and 0.2757 for the subset heritabilities. Whereas all 50K genetic correlations were negative, subset genetic correlations were both, positive and negative (ranging from -0.9324 between MY and mastitis to 0.6662 between MY and digital dermatitis). The subsets containing variants which were annotated as noncoding related, splice sites, untranslated regions, metabolic quantitative trait loci, and young variants ranked highest in terms of their contribution to the traits` genetic variance. We were able to show that linkage disequilibrium between subset variants and adjacent variants did not cause these subsets` high effect.

CONCLUSION

Our results confirm the connection of milk production and health traits in dairy cattle via the animals` metabolic state. In addition, they highlight the potential of including functional information in genomic analyses, which helps to dissect the extent and direction of the observed traits` connection in more detail.

Meta-analysis of identified genomic regions and candidate genes underlying salinity tolerance in rice (Oryza sativa L.)

Rice output has grown globally, yet abiotic factors are still a key cause for worry. Salinity stress seems to have the more impact on crop production out of all abiotic stresses. Currently one of the most significant challenges in paddy breeding for salinity tolerance with the help of QTLs, is to determine the QTLs having the best chance of improving salinity tolerance with the least amount of background noise from the tolerant parent. Minimizing the size of the QTL confidence interval (CI) is essential in order to primarily include the genes responsible for salinity stress tolerance. By considering that, a genome-wide meta-QTL analysis on 768 QTLs from 35 rice populations published from 2001 to 2022 was conducted to identify consensus regions and the candidate genes underlying those regions responsible for the salinity tolerance, as it reduces the confidence interval (CI) to many folds from the initial QTL studies. In the present investigation, a total of 65 MQTLs were extracted with an average CI reduced from 17.35 to 1.66 cM including the smallest of 0.01 cM. Identification of the MQTLs for individual traits and then classifying the target traits into correlated morphological, physiological and biochemical aspects, resulted in more efficient interpretation of the salinity tolerance, identifying the candidate genes and to understand the salinity tolerance mechanism as a whole. The results of this study have a huge potential to improve the rice genotypes for salinity tolerance with the help of MAS and MABC.

Profiling the selected hotspots for ear traits in two maize-teosinte populations

KEY MESSAGE

Eight selected hotspots related to ear traits were identified from two maize-teosinte populations. Throughout the history of maize cultivation, ear-related traits have been selected. However, little is known about the specific genes involved in shaping these traits from their origins in the wild progenitor, teosinte, to the characteristics observed in modern maize. In this study, five ear traits (kernel row number [KRN], ear length [EL], kernel number per row [KNR], cob diameter [CD], and ear diameter [ED]) were investigated, and eight quantitative trait loci (QTL) hotspots were identified in two maize-teosinte populations. Notably, our findings revealed a significant enrichment of genes showing a selection signature and expressed in the ear in qbdCD1.1, qbdCD5.1, qbpCD2.1, qbdED1.1, qbpEL1.1, qbpEL5.1, qbdKNR1.1, and qbdKNR10.1, suggesting that these eight QTL are selected hotspots involved in shaping the maize ear. By combining the results of the QTL analysis with data from previous genome-wide association study (GWAS) involving two natural panels, we identified eight candidate selected genes related to KRN, KNR, CD, and ED. Among these, considering their expression pattern and sequence variation, Zm00001d025111, encoding a WD40/YVTN protein, was proposed as a positive regulator of KNR. This study presents a framework for understanding the genomic distribution of selected loci crucial in determining ear-related traits.

A new winter wheat genetic resource harbors untapped diversity from synthetic hexaploid wheat

KEY MESSAGE

The NIAB_WW_SHW_NAM population, a large nested association mapping panel, is a useful resource for mapping QTL from synthetic hexaploid wheat that can improve modern elite wheat cultivars. The allelic richness harbored in progenitors of hexaploid bread wheat (Triticum aestivum L.) is a useful resource for addressing the genetic diversity bottleneck in modern cultivars. Synthetic hexaploid wheat (SHW) is created through resynthesis of the hybridisation events between the tetraploid (Triticum turgidum subsp. durum Desf.) and diploid (Aegilops tauschii Coss.) bread wheat progenitors. We developed a large and diverse winter wheat nested association mapping (NAM) population (termed the NIAB_WW_SHW_NAM) consisting of 3241 genotypes derived from 54 nested back-cross 1 (BC1) populations, each formed via back-crossing a different primary SHW into the UK winter wheat cultivar 'Robigus'. The primary SHW lines were created using 15 T. durum donors and 47 Ae. tauschii accessions that spanned the lineages and geographical range of the species. Primary SHW parents were typically earlier flowering, taller and showed better resistance to yellow rust infection (Yr) than 'Robigus'. The NIAB_WW_SHW_NAM population was genotyped using a single nucleotide polymorphism (SNP) array and 27 quantitative trait loci (QTLs) were detected for flowering time, plant height and Yr resistance. Across multiple field trials, a QTL for Yr resistance was found on chromosome 4D that corresponded to the Yr28 resistance gene previously reported in other SHW lines. These results demonstrate the value of the NIAB_WW_SHW_NAM population for genetic mapping and provide the first evidence of Yr28 working in current UK environments and genetic backgrounds. These examples, coupled with the evidence of commercial wheat breeders selecting promising genotypes, highlight the potential value of the NIAB_WW_SHW_NAM to variety improvement.