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.

Investigating the role of common cis-regulatory variants in modifying penetrance of putatively damaging, inherited variants in severe neurodevelopmental disorders.

Recent work has revealed an important role for rare, incompletely penetrant inherited coding variants in neurodevelopmental disorders (NDDs). Additionally, we have previously shown that common variants contribute to risk for rare NDDs. Here, we investigate whether common variants exert their effects by modifying gene expression, using multi-cis-expression quantitative trait loci (cis-eQTL) prediction models. We first performed a transcriptome-wide association study for NDDs using 6987 probands from the Deciphering Developmental Disorders (DDD) study and 9720 controls, and found one gene, RAB2A, that passed multiple testing correction (p = 6.7 × 10-7). We then investigated whether cis-eQTLs modify the penetrance of putatively damaging, rare coding variants inherited by NDD probands from their unaffected parents in a set of 1700 trios. We found no evidence that unaffected parents transmitting putatively damaging coding variants had higher genetically-predicted expression of the variant-harboring gene than their child. In probands carrying putatively damaging variants in constrained genes, the genetically-predicted expression of these genes in blood was lower than in controls (p = 2.7 × 10-3). However, results for proband-control comparisons were inconsistent across different sets of genes, variant filters and tissues. We find limited evidence that common cis-eQTLs modify penetrance of rare coding variants in a large cohort of NDD probands.

A wild rice CSSL population facilitated identification of salt tolerance genes and rice germplasm innovation.

Salt stress is one of the major factors that limits rice production. Therefore, identification of salt-tolerant alleles from wild rice is important for rice breeding. In this study, we constructed a set of chromosome segment substitution lines (CSSLs) using wild rice as the donor parent and cultivated rice Nipponbare (Nip) as the recurrent parent. Salt tolerance germinability (STG) was evaluated, and its association with genotypes was determined using this CSSL population. We identified 17 QTLs related to STG. By integrating the transcriptome and genome data, four candidate genes were identified, including the previously reported AGO2 and WRKY53. Compared with Nip, wild rice AGO2 has a structure variation in its promoter region and the expression levels were upregulated under salt treatments; wild rice WRKY53 also has natural variation in its promoter region, and the expression levels were downregulated under salt treatments. Wild rice AGO2 and WRKY53 alleles have combined effects for improving salt tolerance at the germination stage. One CSSL line, CSSL118 that harbors these two alleles was selected. Compared with the background parent Nip, CSSL118 showed comprehensive salt tolerance and higher yield, with improved transcript levels of reactive oxygen species scavenging genes. Our results provided promising genes and germplasm resources for future rice salt tolerance breeding.

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.

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.

Genetic insights into superior grain number traits: a QTL analysis of wheat-Agropyron cristatum derivative pubing3228.

BACKGROUND: Agropyron cristatum (L.) is a valuable genetic resource for expanding the genetic diversity of common wheat. Pubing3228, a novel wheat-A. cristatum hybrid germplasm, exhibits several desirable agricultural traits, including high grain number per spike (GNS). Understanding the genetic architecture of GNS in Pubing3228 is crucial for enhancing wheat yield. This study aims to analyze the specific genetic regions and alleles associated with high GNS in Pubing3228. METHODS: The study employed a recombination inbred line (RIL) population derived from a cross between Pubing3228 and Jing4839 to investigate the genetic regions and alleles linked to high GNS. Quantitative Trait Loci (QTL) analysis and candidate gene investigation were utilized to explore these traits. RESULTS: A total of 40 QTLs associated with GNS were identified across 16 chromosomes, accounting for 4.25-17.17% of the total phenotypic variation. Five QTLs (QGns.wa-1D, QGns.wa-5 A, QGns.wa-7Da.1, QGns.wa-7Da.2 and QGns.wa-7Da.3) accounter for over 10% of the phenotypic variation in at least two environments. Furthermore, 94.67% of the GNS QTL with positive effects originated from Pubing3228. Candidate gene analysis of stable QTLs identified 11 candidate genes for GNS, including a senescence-associated protein gene (TraesCS7D01G148000) linked to the most significant SNP (AX-108,748,734) on chromosome 7D, potentially involved in reallocating nutrients from senescing tissues to developing seeds. CONCLUSION: This study provides new insights into the genetic mechanisms underlying high GNS in Pubing3228, offering valuable resources for marker-assisted selection in wheat breeding to enhance yield.

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.

Genetic determinants of plasma protein levels in the Estonian population.

The proteome holds great potential as an intermediate layer between the genome and phenome. Previous protein quantitative trait locus studies have focused mainly on describing the effects of common genetic variations on the proteome. Here, we assessed the impact of the common and rare genetic variations as well as the copy number variants (CNVs) on 326 plasma proteins measured in up to 500 individuals. We identified 184 cis and 94 trans signals for 157 protein traits, which were further fine-mapped to credible sets for 101 cis and 87 trans signals for 151 proteins. Rare genetic variation contributed to the levels of 7 proteins, with 5 cis and 14 trans associations. CNVs were associated with the levels of 11 proteins (7 cis and 5 trans), examples including a 3q12.1 deletion acting as a hub for multiple trans associations; and a CNV overlapping NAIP, a sensor component of the NAIP-NLRC4 inflammasome which is affecting pro-inflammatory cytokine interleukin 18 levels. In summary, this work presents a comprehensive resource of genetic variation affecting the plasma protein levels and provides the interpretation of identified effects.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Scattered throughout the human genome are variations in the genetic code that make individuals more or less likely to develop certain traits. To identify these variants, scientists carry out Genome-wide association studies (GWAS) which compare the DNA variants of large groups of people with and without the trait of interest. This method has been able to find the underlying genes for many human diseases, but it has limitations. For instance, some variations are linked together due to where they are positioned within DNA, which can result in GWAS falsely reporting associations between genetic variants and traits. This phenomenon, known as linkage equilibrium, can be avoided by analyzing functional genomics which looks at the multiple ways a gene's activity can be influenced by a variation. For instance, how the gene is copied and decoded in to proteins and RNA molecules, and the rate at which these products are generated. Researchers can now use an artificial intelligence technique called deep learning to generate functional genomic data from a particular DNA sequence. Here, Song et al. used one of these deep learning models to calculate the functional genomics of haplotypes, groups of genetic variants inherited from one parent. The approach was applied to DNA samples from over 350 thousand individuals included in the UK BioBank. An activity score, defined as the haplotype function score (or HFS for short), was calculated for at least two haplotypes per individual, and then compared to various complex traits like height or bone density. Song et al. found that the HFS framework was better at finding links between genes and specific traits than existing methods. It also provided more information on the biology that may be underpinning these outcomes. Although more work is needed to reduce the computer processing times required to calculate the HFS, Song et al. believe that their new method has the potential to improve the way researchers identify links between genes and human 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 ).

Trans-eQTL mapping in gene sets identifies network effects of genetic variants.

Nearly all trait-associated variants identified in genome-wide association studies (GWASs) are noncoding. The cis regulatory effects of these variants have been extensively characterized, but how they affect gene regulation in trans has been the subject of fewer studies because of the difficulty in detecting trans-expression quantitative loci (eQTLs). We developed trans-PCO for detecting trans effects of genetic variants on gene networks. Our simulations demonstrate that trans-PCO substantially outperforms existing trans-eQTL mapping methods. We applied trans-PCO to two gene expression datasets from whole blood, DGN (N = 913) and eQTLGen (N = 31,684), and identified 14,985 high-quality trans-eSNP-module pairs associated with 197 co-expression gene modules and biological processes. We performed colocalization analyses between GWAS loci of 46 complex traits and the trans-eQTLs. We demonstrated that the identified trans effects can help us understand how trait-associated variants affect gene regulatory networks and biological pathways.