Meixia Ye

A Model Framework for Identifying Genes that Guide the Evolution of Heterochrony

Heterochrony, the phylogenic change in the time of developmental events or rate of development, has been thought to play an important role in producing phenotypic novelty during evolution. Increasing evidence suggests that specific genes are implicated in heterochrony, guiding the process of developmental divergence, but no quantitative models have been instrumented to map such heterochrony genes. Here, we present a computational framework for genetic mapping by which to characterize and locate quantitative trait loci (QTLs) that govern heterochrony described by four parameters, the timing of the inflection point, the timing of maximum acceleration of growth, the timing of maximum deceleration of growth, and the length of linear growth. The framework was developed from functional mapping, a dynamic model derived to map QTLs for the overall process and pattern of development. By integrating an optimality algorithm, the framework allows the so-called heterochrony QTLs (hQTLs) to be tested and quantified. Specific pipelines are given for testing how hQTLs control the onset and offset of developmental events, the rate of development, and duration of a particular developmental stage. Computer simulation was performed to examine the statistical properties of the model and demonstrate its utility to characterize the effect of hQTLs on population diversification due to heterochrony. By analyzing a genetic mapping data in rice, the framework identified an hQTL that controls the timing of maximum growth rate and duration of linear growth stage in plant height growth. The framework provides a tool to study how genetic variation translates into phenotypic innovation, leading a lineage to evolve, through heterochrony.

Integrating Evolutionary Game Theory into Mechanistic Genotype–Phenotype Mapping

Natural selection has shaped the evolution of organisms toward optimizing their structural and functional design. However, how this universal principle can enhance genotype-phenotype mapping of quantitative traits has remained unexplored. Here we show that the integration of this principle and functional mapping through evolutionary game theory gains new insight into the genetic architecture of complex traits. By viewing phenotype formation as an evolutionary system, we formulate mathematical equations to model the ecological mechanisms that drive the interaction and coordination of its constituent components toward population dynamics and stability. Functional mapping provides a procedure for estimating the genetic parameters that specify the dynamic relationship of competition and cooperation and predicting how genes mediate the evolution of this relationship during trait formation.

2HiGWAS: a unifying high-dimensional platform to infer the global genetic architecture of trait development

Whole-genome search of genes is an essential approach to dissecting complex traits, but a marginal one-single-nucleotide polymorphism (SNP)/one-phenotype regression analysis widely used in current genome-wide association studies fails to estimate the net and cumulative effects of SNPs and reveal the developmental pattern of interplay between genes and traits. Here we describe a computational framework, which we refer to as two-side high-dimensional genome-wide association studies (2HiGWAS), to associate an ultrahigh dimension of SNPs with a high dimension of developmental trajectories measured across time and space. The model is implemented with a dual dimension-reduction procedure for both predictors and responses to select a sparse but full set of significant loci from an extremely large pool of SNPs and estimate their net time-varying effects on trait development. The model can not only help geneticists to precisely identify an entire set of genes underlying complex traits but also allow them to elucidate a global picture of how genes control developmental and dynamic processes of trait formation. We investigated the statistical properties of the model via extensive simulation studies. With the increasing availability of GWAS in various organisms, 2HiGWAS will have important implications for genetic studies of developmental compelx traits.

How trees allocate carbon for optimal growth: insight from a game-theoretic model

How trees allocate photosynthetic products to primary height growth and secondary radial growth reflects their capacity to best use environmental resources. Despite substantial efforts to explore tree height-diameter relationship empirically and through theoretical modeling, our understanding of the biological mechanisms that govern this phenomenon is still limited. By thinking of stem woody biomass production as an ecological system of apical and lateral growth components, we implement game theory to model and discern how these two components cooperate symbiotically with each other or compete for resources to determine the size of a tree stem. This resulting allometry game theory is further embedded within a genetic mapping and association paradigm, allowing the genetic loci mediating the carbon allocation of stemwood growth to be characterized and mapped throughout the genome. Allometry game theory was validated by analyzing a mapping data of stem height and diameter growth over perennial seasons in a poplar tree. Several key quantitative trait loci were found to interpret the process and pattern of stemwood growth through regulating the ecological interactions of stem apical and lateral growth. The application of allometry game theory enables the prediction of the situations in which the cooperation, competition or altruism is an optimal decision of a tree to fully use the environmental resources it owns.

Computational identification of genes modulating stem height–diameter allometry

Summary The developmental variation in stem height with respect to stem diameter is related to a broad range of ecological and evolutionary phenomena in trees, but the underlying genetic basis of this variation remains elusive. We implement a dynamic statistical model, functional mapping, to formulate a general procedure for the computational identification of quantitative trait loci (QTLs) that control stem height–diameter allometry during development. Functional mapping integrates the biological principles underlying trait formation and development into the association analysis of DNA genotype and endpoint phenotype, thus providing an incentive for understanding the mechanistic interplay between genes and development. Built on the basic tenet of functional mapping, we explore two core ecological scenarios of how stem height and stem diameter covary in response to environmental stimuli: (i) trees pioneer sunlit space by allocating more growth to stem height than diameter and (ii) trees maintain their competitive advantage through an inverse pattern. The model is equipped to characterize ‘pioneering’ QTLs (pi QTLs) and ‘maintaining’ QTLs (mi QTLs) which modulate these two ecological scenarios, respectively. In a practical application to a mapping population of full‐sib hybrids derived from two Populus species, the model has well proven its versatility by identifying several pi QTLs that promote height growth at a cost of diameter growth and several mi QTLs that benefit radial growth at a cost of height growth. Judicious application of functional mapping may lead to improved strategies for studying the genetic control of the formation mechanisms underlying trade‐offs among quantities of assimilates allocated to different growth parts.

Functional mapping of seasonal transition in perennial plants

Unlike annuals, all perennial plants undergo seasonal transitions during ontogeny. As an adaptive response to seasonal changes in climate, the seasonal pattern of growth is likely to be under genetic control, although its underlying genetic basis remains unknown. Here, we develop a computational model that can map specific quantitative trait loci (QTLs) responsible for seasonal transitions of growth in perennials. The model is founded on functional mapping, a statistical framework to map developmental dynamics, which is reformed to integrate a seasonally adjusted growth function. The new model is equipped with a capacity to characterize the genetic effects of QTLs on seasonal alternation at different ages and then to better elucidate the genetic architecture of development. The model is implemented with a series of testing procedures, including (i) how a QTL controls an overall ontogenetic growth curve, (ii) how the QTL determines seasonal trajectories of growth within years and (iii) how it determines the dynamic nature of age-specific season response. The model was validated through computer simulation. The extension of season adjustment to other types of biological curves is statistically straightforward, facilitating a wider variety of genetic studies into ontogenetic growth and development in perennial plants.

A multi-Poisson dynamic mixture model to cluster developmental patterns of gene expression by RNA-seq

Dynamic changes of gene expression reflect an intrinsic mechanism of how an organism responds to developmental and environmental signals. With the increasing availability of expression data across a time-space scale by RNA-seq, the classification of genes as per their biological function using RNA-seq data has become one of the most significant challenges in contemporary biology. Here we develop a clustering mixture model to discover distinct groups of genes expressed during a period of organ development. By integrating the density function of multivariate Poisson distribution, the model accommodates the discrete property of read counts characteristic of RNA-seq data. The temporal dependence of gene expression is modeled by the first-order autoregressive process. The model is implemented with the Expectation-Maximization algorithm and model selection to determine the optimal number of gene clusters and obtain the estimates of Poisson parameters that describe the pattern of time-dependent expression of genes from each cluster. The model has been demonstrated by analyzing a real data from an experiment aimed to link the pattern of gene expression to catkin development in white poplar. The usefulness of the model has been validated through computer simulation. The model provides a valuable tool for clustering RNA-seq data, facilitating our global view of expression dynamics and understanding of gene regulation mechanisms.

A high-dimensional linkage analysis model for characterizing crossover interference

Linkage analysis has played an important role in understanding genome structure and evolution. However, two-point linkage analysis widely used for genetic map construction can rarely chart a detailed picture of genome organization because it fails to identify the dependence of crossovers distributed along the length of a chromosome, a phenomenon known as crossover interference. Multi-point analysis, proven to be more advantageous in gene ordering and genetic distance estimation for dominant markers than two-point analysis, is equipped with a capacity to discern and quantify crossover interference. Here, we review a statistical model for four-point analysis, which, beyond three-point analysis, can characterize crossover interference that takes place not only between two adjacent chromosomal intervals, but also over multiple successive intervals. This procedure provides an analytical tool to elucidate the detailed landscape of crossover interference over the genome and further infer the evolution of genome structure and organization.

Evo-Devo-EpiR: a genome-wide search platform for epistatic control on the evolution of development

Evo-devo is a theory proposed to study how phenotypes evolve by comparing the developmental processes of different organisms or the same organism experiencing changing environments. It has been recognized that nonallelic interactions at different genes or quantitative trait loci, known as epistasis, may play a pivotal role in the evolution of development, but it has proven difficult to quantify and elucidate this role into a coherent picture. We implement a high-dimensional genome-wide association study model into the evo-devo paradigm and pack it into the R-based Evo-Devo-EpiR, aimed at facilitating the genome-wide landscaping of epistasis for the diversification of phenotypic development. By analyzing a high-throughput assay of DNA markers and their pairs simultaneously, Evo-Devo-EpiR is equipped with a capacity to systematically characterize various epistatic interactions that impact on the pattern and timing of development and its evolution. Enabling a global search for all possible genetic interactions for developmental processes throughout the whole genome, Evo-Devo-EpiR provides a computational tool to illustrate a precise genotype-phenotype map at interface between epistasis, development and evolution.

AlloMap6: an R package for genetic linkage analysis in allohexaploids

Allopolyploids are a group of polyploids with more than two sets of chromosomes derived from different species. Previous linkage analysis of allopolyploids is based on the assumption that different chromosomes pair randomly during meiosis. A more sophisticated model to relax this assumption has been developed for allotetraploids by incorporating the preferential pairing behavior of homologous over homoeologous chromosomes. Here, we show that the basic principle of this model can be extended to perform linkage analysis of higher-ploidy allohexaploids, where multiple preferential pairing factors are used to characterize chromosomal-pairing meiotic features between different constituent species. We implemented the extended model into an R package, called AlloMap6, allowing the recombination fractions and preferential pairing factors to be estimated simultaneously. Allomap6 has two major functionalities, computer simulation and real-data analysis. By analyzing a real data from a full-sib family of allohexaploid persimmon, we tested and validated the usefulness and utility of this package. AlloMap6 lays a foundation for allohexaploid genetic mapping and provides a new horizon to explore the chromosomal kinship of allohexaploids.