Functional mapping of human growth trajectories

Funmap2: an R package for QTL mapping using longitudinal phenotypes

Quantitative trait locus (QTL) mapping has been used as a powerful tool for inferring the complexity of the genetic architecture that underlies phenotypic traits. This approach has shown its unique power to map the developmental genetic architecture of complex traits by implementing longitudinal data analysis. Here, we introduce the R package Funmap2 based on the functional mapping framework, which integrates prior biological knowledge into the statistical model. Specifically, the functional mapping framework is engineered to include longitudinal curves that describe the genetic effects and the covariance matrix of the trait of interest. Funmap2 chooses the type of longitudinal curve and covariance matrix automatically using information criteria. Funmap2 is available for download at https://github.com/wzhy2000/Funmap2.

Plastic expression of heterochrony quantitative trait loci (hQTLs) for leaf growth in the common bean (Phaseolus vulgaris)

Heterochrony, that is, evolutionary changes in the relative timing of developmental events and processes, has emerged as a key concept that links evolution and development. Genes associated with heterochrony encode molecular components of developmental timing mechanisms. However, our understanding of how heterochrony genes alter the expression of heterochrony in response to environmental changes remains very limited. We applied functional mapping to find quantitative trait loci (QTLs) responsible for growth trajectories of leaf area and leaf mass in the common bean (Phaseolus vulgaris) grown in two contrasting environments. We identified three major QTLs pleiotropically expressed under the two environments. Further characterization of the temporal pattern of these QTLs indicates that they are heterochrony QTLs (hQTLs) in terms of their role in influencing four heterochronic parameters: the timing of the inflection point, the timing of maximum acceleration and deceleration, and the duration of linear growth. The pattern of gene action by the hQTLs on each parameter was unique, being environmentally dependent and varying between two allometrically related leaf growth traits. These results provide new insights into the complexity of genetic mechanisms that control trait formation in plants and provide novel findings that will be of use in studying the evolutionary trends.

Go with the flow-biology and genetics of the lactation cycle

Lactation is a dynamic process, which evolved to meet dietary demands of growing offspring. At the same time, the mother's metabolism changes to meet the high requirements of nutrient supply to the offspring. Through strong artificial selection, the strain of milk production on dairy cows is often associated with impaired health and fertility. This led to the incorporation of functional traits into breeding aims to counteract this negative association. Potentially, distributing the total quantity of milk per lactation cycle more equally over time could reduce the peak of physiological strain and improve health and fertility. During lactation many factors affect the production of milk: food intake; digestion, absorption, and transportation of nutrients; blood glucose levels; activity of cells in the mammary gland, liver, and adipose tissue; synthesis of proteins and fat in the secretory cells; and the metabolic and regulatory pathways that provide fatty acids, amino acids, and carbohydrates. Whilst the endocrine regulation and physiology of the dynamic process of milk production seems to be understood, the genetics that underlie these dynamics are still to be uncovered. Modeling of longitudinal traits and estimating the change in additive genetic variation over time has shown that the genetic contribution to the expression of a trait depends on the considered time-point. Such time-dependent studies could contribute to the discovery of missing heritability. Only very few studies have estimated exact gene and marker effects at different time-points during lactation. The most prominent gene affecting milk yield and milk fat, DGAT1, exhibits its main effects after peak production, whilst the casein genes have larger effects in early lactation. Understanding the physiological dynamics and elucidating the time-dependent genetic effects behind dynamically expressed traits will contribute to selection decisions to further improve productive and healthy breeding populations.

Analyzing growth trajectories

Growth trajectories play a central role in life course epidemiology, often providing fundamental indicators of prenatal or childhood development, as well as an array of potential determinants of adult health outcomes. Statistical methods for the analysis of growth trajectories have been widely studied, but many challenging problems remain. Repeated measurements of length, weight and head circumference, for example, may be available on most subjects in a study, but usually only sparse temporal sampling of such variables is feasible. It can thus be challenging to gain a detailed understanding of growth patterns, and smoothing techniques are inevitably needed. Moreover, the problem is exacerbated by the presence of large fluctuations in growth velocity during early infancy, and high variability between subjects. Existing approaches, however, can be inflexible because of a reliance on parametric models, require computationally intensive methods that are unsuitable for exploratory analyses, or are only capable of examining each variable separately. This article proposes some new nonparametric approaches to analyzing sparse data on growth trajectories, with flexibility and ease of implementation being key features. The methods are illustrated using data on participants in the Collaborative Perinatal Project.

Linking transcriptomic and genomic variation to growth in brook charr hybrids (Salvelinus fontinalis, Mitchill)

Hybridization can lead to phenotypic differences arising from changes in gene expression patterns or new allele combinations. Variation in gene expression is thought to be controlled by differences in transcription regulation of parental alleles, either through cis- or trans-regulatory elements. A previous study among brook charr hybrids from different populations (Rupert, Laval, and domestic) showing distinct length at age during early life stages also revealed different patterns in transcription regulation inheritance of transcript abundance. In the present study, transcript abundance using RNA-sequencing and quantitative real-time PCR, single-nucleotide polymorphism (SNP) genotypes and allelic imbalance were assessed in order to understand the molecular mechanisms underlying the observed transcriptomic and differences in length at age among domestic × Rupert hybrids and Laval × domestic hybrids. We found 198 differentially expressed genes between the two hybrid crosses, and allelic imbalance could be analyzed for 69 of them. Among these 69 genes, 36 genes exhibited cis-acting regulatory effects in both of the two crosses, thus confirming the prevalent role of cis-acting regulatory elements in the regulation of differentially expressed genes among intraspecific hybrids. In addition, we detected a significant association between SNP genotypes of three genes and length at age. Our study is thus one of the few that have highlighted some of the molecular mechanisms potentially involved in the differential phenotypic expression in intraspecific hybrids for nonmodel species.

The influence of parental effects on transcriptomic landscape during early development in brook charr (Salvelinus fontinalis, Mitchill)

Parental effects represent an important source of variation in offspring phenotypes. Depending on the specific mechanisms involved, parental effects may be caused to different degrees by either the maternal or the paternal parent, and these effects may in turn act at different stages of development. To detect parental effects acting on gene transcription regulation and length phenotype during ontogeny, the transcriptomic profiles of two reciprocal hybrids from Laval × Rupert and Laval × Domestic populations of brook charr were compared at hatching, yolk sac resorption and 15 weeks after exogenous feeding. Using a salmonid cDNA microarray, our results show that parental effects modulated gene expression among reciprocal hybrids only at the yolk sac resorption stage. In addition, Laval × Domestic and Laval × Rupert reciprocal hybrids differed in the magnitude of theses parental effects, with 199 and 630 differentially expressed transcripts, respectively. This corresponds to a maximum of 18.5% of the analyzed transcripts. These transcripts are functionally related to cell cycle, nucleic acid metabolism and intracellular protein traffic, which is consistent with observed differences associated with embryonic development and growth differences in other fish species. Our results thus illustrate how parental effects on patterns of gene transcription seem dependent on the genetic architecture of the parents. In addition, in absence of transcriptional differences, non-transcript deposits in the yolk sac could contribute to the observed length differences among the reciprocal hybrids before yolk sac resorption.

Genetic effects and correlations between production and fertility traits and their dependency on the lactation-stage in Holstein Friesians

Background

This study focused on the dynamics of genome-wide effects on five milk production and eight fertility traits as well as genetic correlations between the traits. For 2,405 Holstein Friesian bulls, estimated breeding values (EBVs) were used. The production traits were additionally assessed in 10-day intervals over the first 60 lactation days, as this stage is physiologically the most crucial time in milk production.

Results

SNPs significantly affecting the EBVs of the production traits could be separated into three groups according to the development of the size of allele effects over time: 1) increasing effects for all traits; 2) decreasing effects for all traits; and 3) increasing effects for all traits except fat yield. Most of the significant markers were found within 22 haplotypes spanning on average 135,338 bp. The DGAT1 region showed high density of significant markers, and thus, haplotype blocks. Further functional candidate genes are proposed for haplotype blocks of significant SNPs (KLHL8, SICLEC12, AGPAT6 and NID1). Negative genetic correlations were found between yield and fertility traits, whilst content traits showed positive correlations with some fertility traits. Genetic correlations became stronger with progressing lactation. When correlations were estimated within genotype classes, correlations were on average 0.1 units weaker between production and fertility traits when the yield increasing allele was present in the genotype.

Conclusions

This study provides insight into the expression of genetic effects during early lactation and suggests possible biological explanations for the presented time-dependent effects. Even though only three markers were found with effects on fertility, the direction of genetic correlations within genotype classes between production and fertility traits suggests that alleles increasing the milk production do not affect fertility in a more negative way compared to the decreasing allele.

A maximum likelihood approach to functional mapping of longitudinal binary traits

Despite their importance in biology and biomedicine, genetic mapping of binary traits that change over time has not been well explored. In this article, we develop a statistical model for mapping quantitative trait loci (QTLs) that govern longitudinal responses of binary traits. The model is constructed within the maximum likelihood framework by which the association between binary responses is modeled in terms of conditional log odds-ratios. With this parameterization, the maximum likelihood estimates (MLEs) of marginal mean parameters are robust to the misspecification of time dependence. We implement an iterative procedures to obtain the MLEs of QTL genotype-specific parameters that define longitudinal binary responses. The usefulness of the model was validated by analyzing a real example in rice. Simulation studies were performed to investigate the statistical properties of the model, showing that the model has power to identify and map specific QTLs responsible for the temporal pattern of binary traits.

Gene-environment interactions on growth trajectories

It has been suggested that children with larger brains tend to perform better on IQ tests or cognitive function tests. Prenatal head growth and head growth in infancy are two crucial periods for subsequent intelligence. Studies have shown that environmental exposure to air pollutants during pregnancy is associated with fetal growth reduction, developmental delay, and reduced IQ. Meanwhile, genetic polymorphisms may modify the effect of environment on head growth. However, studies on gene-environment or gene-gene interactions on growth trajectories have been quite limited partly due to the difficulty to quantitatively measure interactions on growth trajectories. Moreover, it is known that assessing the significance of gene-environment or gene-gene interactions on cross-sectional outcomes empirically using the permutation procedures may bring substantial errors in the tests. We proposed a score that quantitatively measures interactions on growth trajectories and developed an algorithm with a parametric bootstrap procedure to empirically assess the significance of the interactions on growth trajectories under the likelihood framework. We also derived a Wald statistic to test for interactions on growth trajectories and compared it to the proposed parametric bootstrap procedure. Through extensive simulation studies, we demonstrated the feasibility and power of the proposed testing procedures. We applied our method to a real dataset with head circumference measures from birth to age 7 on a cohort currently being conducted by the Columbia Center for Children's Environmental Health (CCCEH) in Krakow, Poland, and identified several significant gene-environment interactions on head circumference growth trajectories.

How to cluster gene expression dynamics in response to environmental signals

Organisms usually cope with change in the environment by altering the dynamic trajectory of gene expression to adjust the complement of active proteins. The identification of particular sets of genes whose expression is adaptive in response to environmental changes helps to understand the mechanistic base of gene-environment interactions essential for organismic development. We describe a computational framework for clustering the dynamics of gene expression in distinct environments through Gaussian mixture fitting to the expression data measured at a set of discrete time points. We outline a number of quantitative testable hypotheses about the patterns of dynamic gene expression in changing environments and gene-environment interactions causing developmental differentiation. The future directions of gene clustering in terms of incorporations of the latest biological discoveries and statistical innovations are discussed. We provide a set of computational tools that are applicable to modeling and analysis of dynamic gene expression data measured in multiple environments.

Curve-based multivariate distance matrix regression analysis: application to genetic association analyses involving repeated measures

Most, if not all, human phenotypes exhibit a temporal, dosage-dependent, or age effect. Despite this fact, it is rare that data are collected over time or in sequence in relevant studies of the determinants of these phenotypes. The costs and organizational sophistication necessary to collect repeated measurements or longitudinal data for a given phenotype are clearly impediments to this, but greater efforts in this area are needed if insights into human phenotypic expression are to be obtained. Appropriate data analysis methods for genetic association studies involving repeated or longitudinal measures are also needed. We consider the use of longitudinal profiles obtained from fitted functions on repeated data collections from a set of individuals whose similarities are contrasted between sets of individuals with different genotypes to test hypotheses about genetic influences on time-dependent phenotype expression. The proposed approach can accommodate uncertainty of the fitted functions, as well as weighting factors across the time points, and is easily extended to a wide variety of complex analysis settings. We showcase the proposed approach with data from a clinical study investigating human blood vessel response to tyramine. We also compare the proposed approach with standard analytic procedures and investigate its robustness and power via simulation studies. The proposed approach is found to be quite flexible and performs either as well or better than traditional statistical methods.