Bayesian functional mapping of dynamic quantitative traits

Functional mapping has emerged as a powerful tool for mapping quantitative trait loci (QTL) that control developmental patterns of complex dynamic traits. Original functional mapping has been constructed within the context of simple interval mapping, without consideration of separate multiple linked QTL for a dynamic trait. In this article, we present a statistical framework for mapping QTL that affect dynamic traits by capitalizing on the strengths of functional mapping and composite interval mapping. Within this so-called composite functional-mapping framework, functional mapping models the time-dependent genetic effects of a QTL tested within a marker interval using a biologically meaningful parametric function, whereas composite interval mapping models the time-dependent genetic effects of the markers outside the test interval to control the genome background using a flexible nonparametric approach based on Legendre polynomials. Such a semiparametric framework was formulated by a maximum-likelihood model and implemented with the EM algorithm, allowing for the estimation and the test of the mathematical parameters that define the QTL effects and the regression coefficients of the Legendre polynomials that describe the marker effects. Simulation studies were performed to investigate the statistical behavior of composite functional mapping and compare its advantage in separating multiple linked QTL as compared to functional mapping. We used the new mapping approach to analyze a genetic mapping example in rice, leading to the identification of multiple QTL, some of which are linked on the same chromosome, that control the developmental trajectory of leaf age.

This study provides the first mapping of mycoprotein functionality, composition and structure throughout the Quorn fermentation process. The fermentation broth, RNA-reduced broth (RNA-broth), centrate and their centrifugation deposits and supernatants were characterised. The broth, RNA-broth and their deposits displayed high concentrations of fungal filaments, which contributed to their high gelling properties (with a 5,320 Pa elastic modulus reported for RNA-broth deposits gels). Foams prepared with RNA-broth and centrate supernatants via frothing exhibited high stability (380 min), with high concentrations of a foam-positive cerato-platanin reported in these samples. Emulsions prepared with the broth and broth supernatant showed high emulsifying activity and stability indexes (12.80 m2/g and 15.84 mins for the broth supernatant) and low oil droplet sizes (18.09 µm for the broth). This study identified previously unreported gelling, foaming and/or emulsifying properties for the different Quorn streams, highlighting opportunities to develop novel sustainable alternatives to animal-derived functional ingredients using mycoprotein material.

Background. Clinical trials in cancer frequently include cancer-specific measures of health but not preference-based measures such as the EQ-5D that are suitable for economic evaluation. Mapping functions have been developed to predict EQ-5D values from these measures, but there is considerable uncertainty about the most appropriate model to use, and many existing models are poor at predicting EQ-5D values. This study aims to investigate a range of potential models to develop mapping functions from 2 widely used cancer-specific measures (FACT-G and EORTC-QLQ-C30) and to identify the best model. Methods. Mapping models are fitted to predict EQ-5D-3L values using ordinary least squares (OLS), tobit, 2-part models, splining, and to EQ-5D item-level responses using response mapping from the FACT-G and QLQ-C30. A variety of model specifications are estimated. Model performance and predictive ability are compared. Analysis is based on 530 patients with various cancers for the FACT-G and 771 patients with multiple myeloma, breast cancer, and lung cancer for the QLQ-C30. Results. For FACT-G, OLS models most accurately predict mean EQ-5D values with the best predicting model using FACT-G items with similar results using tobit. Response mapping has low predictive ability. In contrast, for the QLQ-C30, response mapping has the most accurate predictions using QLQ-C30 dimensions. The QLQ-C30 has better predicted EQ-5D values across the range of possible values; however, few respondents in the FACT-G data set have low EQ-5D values, which reduces the accuracy at the severe end. Conclusions. OLS and tobit mapping functions perform well for both instruments. Response mapping gives the best model predictions for QLQ-C30. The generalizability of the FACT-G mapping function is limited to populations in moderate to good health.

The biological and statistical advantages of functional mapping result from joint modeling of the mean-covariance structures for developmental trajectories of a complex trait measured at a series of time points. While an increased number of time points can better describe the dynamic pattern of trait development, significant difficulties in performing functional mapping arise from prohibitive computational times required as well as from modeling the structure of a high-dimensional covariance matrix. In this article, we develop a statistical model for functional mapping of quantitative trait loci (QTL) that govern the developmental process of a quantitative trait on the basis of wavelet dimension reduction. By breaking an original signal down into a spectrum by taking its averages (smooth coefficients) and differences (detail coefficients), we used the discrete Haar wavelet shrinkage technique to transform an inherently high-dimensional biological problem into its tractable low-dimensional representation within the framework of functional mapping constructed by a Gaussian mixture model. Unlike conventional nonparametric modeling of wavelet shrinkage, we incorporate mathematical aspects of developmental trajectories into the smooth coefficients used for QTL mapping, thus preserving the biological relevance of functional mapping in formulating a number of hypothesis tests at the interplay between gene actions/interactions and developmental patterns for complex phenotypes. This wavelet-based parametric functional mapping has been statistically examined and compared with full-dimensional functional mapping through simulation studies. It holds great promise as a powerful statistical tool to unravel the genetic machinery of developmental trajectories with large-scale high-dimensional data.

Background: Plants have evolved the dual capacity for maximizing light assimilation through stem growth (phototropism) and maximizing water and nutrient absorption through root growth (gravitropism). Previous studies have revealed the physiological and molecular mechanisms of these two processes, but the genetic basis for how gravitropism and phototropism interact and coordinate with one another to determine plant growth remains poorly understood. Methods: We designed a seed germination experiment using a full-sib F1 family of Populus euphratica to simultaneously monitor the gravitropic growth of the radicle and the phototropic growth of the plumule throughout seedling ontogeny. We implemented three functional mapping models to identify quantitative trait loci (QTLs) that regulate gravitropic and phototropic growth. Univariate functional mapping dissected each growth trait separately, bivariate functional mapping mapped two growth traits simultaneously, and composite functional mapping mapped the sum of gravitropic and phototropic growth as a main axis. Results: Bivariate model detected 8 QTLs for gravitropism and phototropism (QWRF, GLUR, F-box, PCFS4, UBQ, TAF12, BHLH95, TMN8), composite model detected 7 QTLs for growth of main axis (ATL8, NEFH, PCFS4, UBQ, SOT16, MOR1, PCMP-H), of which, PCFS4 and UBQ were pleiotropically detected with the both model. Many of these QTLs are situated within the genomic regions of candidate genes. Conclusions: The results from our models provide new insight into the mechanisms of genetic control of gravitropism and phototropism in a desert tree, and will stimulate our understanding of the relationships between gravity and light signal transduction pathways and tree adaptation to arid soil.

Various genetic map functions have been proposed to infer the unobservable genetic distance between two loci from the observable recombination fraction between them. Some map functions were found to fit data better than others. When there are more than three markers, multilocus recombination probabilities cannot be uniquely determined by the defining property of map functions, and different methods have been proposed to permit the use of map functions to analyze multilocus data. If for a given map function, there is a probability model for recombination that can give rise to it, then joint recombination probabilities can be deduced from this model. This provides another way to use map functions in multilocus analysis. In this paper we show that stationary renewal processes give rise to most of the map functions in the literature. Furthermore, we show that the interevent distributions of these renewal processes can all be approximated quite well by gamma distributions.

SummaryQuantitative trait loci (QTLs) may affect not only the mean of a trait but also its variability. A special aspect is the variability between multiple measured traits of genotyped animals, such as the within-litter variance of piglet birth weights. The sample variance of repeated measurements is assigned as an observation for every genotyped individual. It is shown that the conditional distribution of the non-normally distributed trait can be approximated by a gamma distribution. To detect QTL effects in the daughter design, a generalized linear model with the identity link function is applied. Suitable test statistics are constructed to test the null hypothesis H0: No QTL with effect on the within-litter variance is segregating versus HA: There is a QTL with effect on the variability of birth weight within litter. Furthermore, estimates of the QTL effect and the QTL position are introduced and discussed. The efficiency of the presented tests is compared with a test based on weighted regression. The error probability of the first type as well as the power of QTL detection are discussed and compared for the different tests.

Adaptation of domesticated species to diverse agroclimatic regions has led to abundant trait diversity. However, the resulting population structure and genetic heterogeneity confounds association mapping of adaptive traits. To address this challenge in sorghum [Sorghum bicolor (L.) Moench]-a widely adapted cereal crop-we developed a nested association mapping (NAM) population using 10 diverse global lines crossed with an elite reference line RTx430. We characterized the population of 2214 recombinant inbred lines at 90,000 SNPs using genotyping-by-sequencing. The population captures ∼70% of known global SNP variation in sorghum, and 57,411 recombination events. Notably, recombination events were four- to fivefold enriched in coding sequences and 5' untranslated regions of genes. To test the power of the NAM population for trait dissection, we conducted joint linkage mapping for two major adaptive traits, flowering time and plant height. We precisely mapped several known genes for these two traits, and identified several additional QTL. Considering all SNPs simultaneously, genetic variation accounted for 65% of flowering time variance and 75% of plant height variance. Further, we directly compared NAM to genome-wide association mapping (using panels of the same size) and found that flowering time and plant height QTL were more consistently identified with the NAM population. Finally, for simulated QTL under strong selection in diversity panels, the power of QTL detection was up to three times greater for NAM vs. association mapping with a diverse panel. These findings validate the NAM resource for trait mapping in sorghum, and demonstrate the value of NAM for dissection of adaptive traits.

To provide the theoretical basis for researching growth, development, and molecular marker-assisted breeding of the economically important Yellow River carp (Cyprinus carpio haematopterus) using dynamic quantitative trait locus (QTL) mapping, we constructed three genetic linkage maps from 207 progeny using a new modified genotyping-by-sequencing method. The three maps contained 16,886, 16,548, and 7482 single nucleotide polymorphism markers, respectively, with an average interval of 0.36 cM, 0.45 cM, and 1.00 cM. We identified 148 QTLs related to four growth traits that were located on 25 chromosomes from three growth stages of Yellow River carp. A total of 32, 36, 43, and 37 QTLs were associated with body length, height, width, and weight, respectively. Among them, 47 QTLs were detected for only one growth trait in one stage, but all of the other QTLs were co-localized. Of the 14 main QTLs, 13 were located on chromosome 12, which suggests the presence of growth-related genes on this chromosome. We then detected 17 candidate genes within 50 K upstream and downstream of the 14 main QTLs. This is the first report of the dynamic QTL mapping of growth traits of Yellow River carp, and the results can be used in future studies of growth, development, and molecular-assisted breeding of this species.

A complex binary trait is a character that has a dichotomous expression but with a polygenic genetic background. Mapping quantitative trait loci (QTL) for such traits is difficult because of the discrete nature and the reduced variation in the phenotypic distribution. Bayesian statistics are proved to be a powerful tool for solving complicated genetic problems, such as multiple QTL with nonadditive effects, and have been successfully applied to QTL mapping for continuous traits. In this study, we show that Bayesian statistics are particularly useful for mapping QTL for complex binary traits. We model the binary trait under the classical threshold model of quantitative genetics. The Bayesian mapping statistics are developed on the basis of the idea of data augmentation. This treatment allows an easy way to generate the value of a hypothetical underlying variable (called the liability) and a threshold, which in turn allow the use of existing Bayesian statistics. The reversible jump Markov chain Monte Carlo algorithm is used to simulate the posterior samples of all unknowns, including the number of QTL, the locations and effects of identified QTL, genotypes of each individual at both the QTL and markers, and eventually the liability of each individual. The Bayesian mapping ends with an estimation of the joint posterior distribution of the number of QTL and the locations and effects of the identified QTL. Utilities of the method are demonstrated using a simulated outbred full-sib family. A computer program written in FORTRAN language is freely available on request.

Most of quantitative trait loci (QTL) underlying soybean seed isoflavone contents were derived from the harvest stage of plant development, which uncover the genetic effects that were expressed in earlier seed developmental stages. The aim of this study was to detect conditional QTL associated with isoflavone accumulation during the entire seed development. A total of 112 recombinant inbred lines developed from the cross between ‘Zhongdou 27’ (higher seed isoflavone content) and ‘Jiunong 20’ (lower seed isoflavone content) were used for the conditional QTL analysis of daidzein (DZ), genistein (GT), glycitein (GC) and total isoflavone (TI) accumulations through composite interval mapping with mixed genetic model. The results indicated that the number and type of QTL and their additive effects for individual and total isoflavone accumulations were different among R3 to R8 developmental stages. Three unconditional QTL and six conditional QTL for DZ, four unconditional QTL and five conditional QTL for GT, six unconditional QTL and five conditional QTL for GC, six unconditional QTL and seven conditional QTL for TI were identified at different developmental stages, respectively. Unconditional and conditional QTL that affect individual and total isoflavone accumulations exhibited multiple expression patterns, implying that some QTL are active for long period and others are transient. Two genomic regions, Satt144‐Satt569 in linkage group F (LG F; chromosome 13, chr 13) for DZ, GC, GT and TI accumulations and Satt540‐Sat_240 in LG M (chr 07) for TI and GC accumulations, were found to significantly affect individual and total isoflavone accumulations in multiple developmental stages, suggesting that the accumulation of soybean seed isoflavones is governed by time‐dependent gene expression.

BackgroundThis study aimed to obtain health utility parameters among Chinese breast cancer patients in different disease states for subsequent health economics model. In addition, we aimed to explore the feasibility of establishing a breast cancer health utility mapping model in China.MethodsMultiple patient-reported health attributes were assessed, including quality of life, which was measured by the Functional Assessment of Cancer Therapy-Breast (FACT-B) instrument; health utility and self-rated health, which were measured by the EuroQol-5 Dimension-5 Level (EQ-5D-5L) questionnaire. Multivariate regression models, including a linear regression model, an ordinal logistic regression model and a Tobit model, were employed to analyze health differences among 446 breast cancer patients. Subgroup analyses were performed to examine differences in multiple dimensions of health derived from the FACT-B and EQ-5D-5L instruments. A mapping function was used to estimate health utility from quality of life. Rank correlation analyses were employed to examine the correlation between estimated and observed health utility values.ResultsA total of 446 breast cancer patients with different disease states were analyzed. The health utility values of breast cancer patients in the P state (without cancer recurrence and metastasis), R state (with cancer recurrence within a year), S state (with primary and recurrent breast cancer for the second year and above), and M state (metastatic cancer) were 0.81 (SD ± 0.23), 0.90 (SD ± 0.12), 0.78 (SD ± 0.31), and 0.74 (SD ± 0.27), respectively. There were positive correlations between all scores, including every domain of the FACT-B instrument (p < 0.001). Results from multivariate analysis suggested that patients in the R and M states had lower scores for overall quality of life (R, β = − 9.45, p < 0.01; M, β = − 6.72, p < 0.05). Patients in the M state had lower health utility values than patients in the P state (β = − 0.11, p < 0.05). Estimated health utility values, which were derived from quality of life by using a mapping function, were significantly correlated with directly measured health utility values (p < 0.001).ConclusionsWe obtained the health utility and health-related quality of life (HRQoL) scores of Chinese breast cancer patients in different disease states. Mapping health utility values from quality of life using four disease states could be feasible in health economic modelling, but the mapping function may need further revision.

BackgroundWith the plummeting cost of the next-generation sequencing technologies, high-density genetic linkage maps could be constructed in a forest hybrid F1 population. However, based on such genetic maps, quantitative trait loci (QTL) mapping cannot be directly conducted with traditional statistical methods or tools because the linkage phase and segregation pattern of molecular markers are not always fixed as in inbred lines.ResultsWe implemented the traditional composite interval mapping (CIM) method to multivariate trait data in forest trees and developed the corresponding software, mvqtlcim. Our method not only incorporated the various segregations and linkage phases of molecular markers, but also applied Takeuchi’s information criterion (TIC) to discriminate the QTL segregation type among several possible alternatives. QTL mapping was performed in a hybrid F1 population of Populus deltoides and P. simonii, and 12 QTLs were detected for tree height over 6 time points. The software package allowed many options for parameters as well as parallel computing for permutation tests. The features of the software were demonstrated with the real data analysis and a large number of Monte Carlo simulations.ConclusionsWe provided a powerful tool for QTL mapping of multiple or longitudinal traits in an outbred F1 population, in which the traditional software for QTL mapping cannot be used. This tool will facilitate studying of QTL mapping and thus will accelerate molecular breeding programs especially in forest trees. The tool package is freely available from https://github.com/tongchf /mvqtlcim.

PTPRD pleiotropy, genetically driven childbirth timing, and corpus callosum microstructure as potential mechanisms underlying ADHD-RLS comorbidity.

While the association between smoking and accelerated facial aging is well documented, the specific pathways underlying this association remain poorly understood. To investigate the shared genetic architecture between smoking and facial aging, we performed genetic analyses based on genome-wide association studies (GWAS) data. These analyses included linkage disequilibrium score regression (LDSC), pleiotropic analysis under composite null hypothesis (PLACO), functional mapping and annotation (FUMA), and multi-marker analysis of genomic annotation (MAGMA). To further explore the shared target genes, we utilized expression quantitative trait loci (eQTLs) and mediation Mendelian randomization (MR) analysis, with subsequent validation conducted through in vitro experiments using NIH/3T3 cells. Additionally, we carried out pan-cancer correlation analyses to assess the broader implications of the identified genes in cancer biology. Through pleiotropy and colocalization analyses, IREB2, along with CHRNA5 and AARS1, were identified as having strong evidence linking smoking and facial aging. Functional enrichment, tissue-specific analyses, and gene co-expression network were conducted to further elucidate the functions of these genes. Following eQTLs and mediation analyses, IREB2 was identified as a potential mediator connecting smoking to facial aging. Cellular experiments demonstrated that exposure to cigarette smoke particles induces cellular senescence and downregulates IREB2 expression. The pan-cancer analysis highlighted IREB2's role in shaping the tumor microenvironment and influencing immune processes. This study identifies IREB2 as a critical factor in the molecular mechanisms by which smoking accelerates facial aging, while also contributing to tumor development and immune evasion. Further functional exploration of IREB2 could uncover new therapeutic avenues to address these interconnected conditions.