Functional mapping of ontogeny in flowering plants.

All organisms face the problem of how to perform a sequence of developmental changes and transitions during ontogeny. We revise functional mapping, a statistical model originally derived to map genes that determine developmental dynamics, to take into account the entire process of ontogenetic growth from embryo to adult and from the vegetative to reproductive phase. The revised model provides a framework that reconciles the genetic architecture of development at different stages and elucidates a comprehensive picture of the genetic control mechanisms of growth that change gradually from a simple to a more complex level. We use an annual flowering plant, as an example, to demonstrate our model by which to map genes and their interactions involved in embryo and postembryonic growth. The model provides a useful tool to study the genetic control of ontogenetic growth in flowering plants and any other organisms through proper modifications based on their biological characteristics.

A computational algorithm for functional clustering of proteome dynamics during development.

Phenotypic traits, such as seed development, are a consequence of complex biochemical interactions among genes, proteins and metabolites, but the underlying mechanisms that operate in a coordinated and sequential manner remain elusive. Here, we address this issue by developing a computational algorithm to monitor proteome changes during the course of trait development. The algorithm is built within the mixture-model framework in which each mixture component is modeled by a specific group of proteins that display a similar temporal pattern of expression in trait development. A nonparametric approach based on Legendre orthogonal polynomials was used to fit dynamic changes of protein expression, increasing the power and flexibility of protein clustering. By analyzing a dataset of proteomic dynamics during early embryogenesis of the Chinese fir, the algorithm has successfully identified several distinct types of proteins that coordinate with each other to determine seed development in this forest tree commercially and environmentally important to China. The algorithm will find its immediate applications for the characterization of mechanistic underpinnings for any other biological processes in which protein abundance plays a key role.

Genome-Wide Identification and Analysis of the WRKY Gene Family in Asparagus officinalis.

In recent years, the related research of the WRKY gene family has been gradually promoted, which is mainly reflected in the aspects of environmental stress and hormone response. However, to make the study of the WRKY gene family more complete, we also need to focus on the whole-genome analysis and identification of the family. In previous studies, the whole WRKY gene family of Arabidopsis, legumes and other plants has been thoroughly studied. However, since the publication of Asparagus officinalis genome-wide data, there has never been an analysis of the whole WRKY gene family. To understand more broadly the function of the WRKY gene family, the whole genome and salt stress transcriptome data of asparagus were used for comprehensive analysis in this study, including WRKY gene family identification, phylogenetic tree construction, analysis of conserved mods and gene domains, extraction of cis-acting elements, intron/exon analysis, species collinearity analysis, and WRKY expression analysis under salt stress. The results showed that a total of 70 genes were selected and randomly distributed on 10 chromosomes and one undefined chromosome. According to the functional classification of Arabidopsis thaliana, the WRKY family of asparagus was divided into 11 subgroups (C1-C9, U1, U2). It is worth considering that the distribution rules of gene-conserved motifs, gene domains and introns/exons in the same subfamily are similar, which suggests that genes in the same subfamily may regulate similar physiological processes. In this study, 11 cis-acting elements of WRKY family were selected, among which auxin, gibberellin, abscisic acid, salicylic acid and other hormone-regulated induction elements were involved. In addition, environmental stress (such as drought stress and low-temperature response) also accounted for a large proportion. Interestingly, we analyzed a total of two tandem duplicate genes and 13 segmental duplication genes, suggesting that this is related to the amplification of the WRKY gene family. Transcriptome data analysis showed that WRKY family genes could regulate plant growth and development by up-regulating and down-regulating gene expression under salt stress. Volcanic maps showed that 3 and 15 AoWRKY genes were significantly up-regulated or down-regulated in NI&NI+S and AMF&AMF+S, respectively. These results provide a new way to analyze the evolution and function of the WRKY gene family, and can provide a reference for the production and research of asparagus.

Genome-Wide Identification and Analysis of the GRAS Transcription Factor Gene Family in Theobroma cacao.

GRAS genes exist widely and play vital roles in various physiological processes in plants. In this study, to identify Theobroma cacao (T. cacao) GRAS genes involved in environmental stress and phytohormones, we conducted a genome-wide analysis of the GRAS gene family in T. cacao. A total of 46 GRAS genes of T. cacao were identified. Chromosomal distribution analysis showed that all the TcGRAS genes were evenly distributed on ten chromosomes. Phylogenetic relationships revealed that GRAS proteins could be divided into twelve subfamilies (HAM: 6, LISCL: 10, LAS: 1, SCL4/7: 1, SCR: 4, DLT: 1, SCL3: 3, DELLA: 4, SHR: 5, PAT1: 6, UN1: 1, UN2: 4). Of the T. cacao GRAS genes, all contained the GRAS domain or GRAS superfamily domain. Subcellular localization analysis predicted that TcGRAS proteins were located in the nucleus, chloroplast, and endomembrane system. Gene duplication analysis showed that there were two pairs of tandem repeats and six pairs of fragment duplications, which may account for the rapid expansion in T. cacao. In addition, we also predicted the physicochemical properties and cis-acting elements. The analysis of GO annotation predicted that the TcGRAS genes were involved in many biological processes. This study highlights the evolution, diversity, and characterization of the GRAS genes in T. cacao and provides the first comprehensive analysis of this gene family in the cacao genome.

Genome-Wide Identification and Analysis of the R2R3-MYB Gene Family in Theobroma cacao.

The MYB gene family is involved in the regulation of plant growth, development and stress responses. In this paper, to identify Theobroma cacao R2R3-MYB (TcMYB) genes involved in environmental stress and phytohormones, we conducted a genome-wide analysis of the R2R3-MYB gene family in Theobroma cacao (cacao). A total of 116 TcMYB genes were identified, and they were divided into 23 subgroups according to the phylogenetic analysis. Meanwhile, the conserved motifs, gene structures and cis-acting elements of promoters were analyzed. Moreover, these TcMYB genes were distributed on 10 chromosomes. We conducted a synteny analysis to understand the evolution of the cacao R2R3-MYB gene family. A total of 37 gene pairs of TcMYB genes were identified through tandem or segmental duplication events. Additionally, we also predicted the subcellular localization and physicochemical properties. All the studies showed that TcMYB genes have multiple functions, including responding to environmental stresses. The results provide an understanding of R2R3-MYB in Theobroma cacao and lay the foundation for a further functional analysis of TcMYB genes in the growth of cacao.

Genome-Wide Identification and Analysis of the MADS-Box Gene Family in Theobroma cacao.

The MADS-box family gene is a class of transcription factors that have been extensively studied and involved in several plant growth and development processes, especially in floral organ specificity, flowering time and initiation and fruit development. In this study, we identified 69 candidate MADS-box genes and clustered these genes into five subgroups (Mα: 11; Mß: 2; Mγ: 14; Mδ: 9; MIKC: 32) based on their phylogenetical relationships with Arabidopsis. Most TcMADS genes within the same subgroup showed a similar gene structure and highly conserved motifs. Chromosomal distribution analysis revealed that all the TcMADS genes were evenly distributed in 10 chromosomes. Additionally, the cis-acting elements of promoter, physicochemical properties and subcellular localization were also analyzed. This study provides a comprehensive analysis of MADS-box genes in Theobroma cacao and lays the foundation for further functional research.

Genome-wide association studies of callus differentiation for the desert tree, Populus euphratica.

Callus differentiation is a key developmental process in plant regeneration from cells. A better understanding of the genetic architecture of callus differentiation timing can help improve tissue transformation and the efficiency of artificial propagation. In this study, we investigated genotypic variation in callus differentiation capacity among 297 diverse P. euphratica trees sampled from a natural population. We employed a genome-wide association study (GWAS) of binary and growth-based parameters to identify loci and characterize the genetic architecture and genetic network underlying regulation of callus differentiation in P. euphratica. The results of this GWAS experiment suggested potential associations controlling whether the callus could differentiate and the process of callus differentiation. We identified multiple significant quantitative trait loci (QTLs), including the genes LOG1 and LOG7 and a locus containing WOX1. We reconstructed a genetic network that visualizes how each QTL interacts uniquely with other variants, and several core QTLs were detected that are involved in the degree of callus differentiation, providing potential targets for selection. This study represents one of the first to identify genetic variants affecting callus differentiation in a forest tree. Our results suggest that callus differentiation may be a typical qualitative-quantitative trait controlled by a major gene as well as polygenes across the genome of P. euphratica. This GWAS will help to design more complex and specific molecular tools for systematically manipulating organ regeneration.

Genome-Wide Analysis of the NAC Domain Transcription Factor Gene Family in Theobroma cacao.

As a plant-specific transcription factor, the NAC (NAM, ATAF1/2 and CUC2) domain protein plays an important role in plant growth and development, as well as stress resistance. Based on the genomic data of the cacao tree, this study identified 102 cacao NAC genes and named them according to their location within the genome. The phylogeny of the protein sequence of the cacao tree NAC family was analyzed using various bioinformatic methods, and then divided into 12 subfamilies. Then, the amino-acid composition, physicochemical properties, genomic location, gene structure, conserved domains, and promoter cis-acting elements were analyzed. This study provides information on the evolution of the TcNAC gene and its possible functions, laying the foundation for further research on the NAC family.

Lysine-specific histone demethylase 1 inhibition promotes reprogramming by facilitating the expression of exogenous transcriptional factors and metabolic switch.

Lysine-specific histone demethylase 1 (LSD1) regulates histone methylation and influences the epigenetic state of cells during the generation of induced pluripotent stem cells (iPSCs). Here we reported that LSD1 inhibition via shRNA or specific inhibitor, tranylcypromine, promoted reprogramming at early stage via two mechanisms. At early stage of reprogramming, LSD1 inhibition increased the retrovirus-mediated exogenous expression of Oct4, Klf4, and Sox2 by blocking related H3K4 demethylation. Since LSD1 inhibition still promoted reprogramming even when iPSCs were induced with small-molecule compounds in a virus-free system, additional mechanisms should be involved. When RNA-seq was used for analysis, it was found that LSD1 inhibition reversed some gene expression changes induced by OKS, which subsequently promoted reprogramming. For example, by partially rescuing the decreased expression of Hif1α, LSD1 inhibition reversed the up-regulation of genes in oxidative phosphorylation pathway and the down-regulation of genes in glycolysis pathway. Such effects facilitated the metabolic switch from oxidative phosphorylation to glycolysis and subsequently promoted iPSCs induction. In addition, LSD1 inhibition also promoted the conversion from pre-iPSCs to iPSCs by facilitating the similar metabolic switch. Therefore, LSD1 inhibition promotes reprogramming by facilitating the expression of exogenous transcriptional factors and metabolic switch.

Systems mapping for hematopoietic progenitor cell heterogeneity.

Cells with the same genotype growing under the same conditions can show different phenotypes, which is known as "population heterogeneity". The heterogeneity of hematopoietic progenitor cells has an effect on their differentiation potential and lineage choices. However, the genetic mechanisms governing population heterogeneity remain unclear. Here, we present a statistical model for mapping the quantitative trait locus (QTL) that affects hematopoietic cell heterogeneity. This strategy, termed systems mapping, integrates a system of differential equations into the framework for systems mapping, allowing hypotheses regarding the interplay between genetic actions and cell heterogeneity to be tested. A simulation approach based on cell heterogeneity dynamics has been designed to test the statistical properties of the model. This model not only considers the traditional QTLs, but also indicates the methylated QTLs that can illustrate non-genetic individual differences. It has significant implications for probing the molecular, genetic and epigenetic mechanisms of hematopoietic progenitor cell heterogeneity.

Sequential introduction of reprogramming factors reveals a time-sensitive requirement for individual factors and a sequential EMT-MET mechanism for optimal reprogramming.

Present practices for reprogramming somatic cells to induced pluripotent stem cells involve simultaneous introduction of reprogramming factors. Here we report that a sequential introduction protocol (Oct4-Klf4 first, then c-Myc and finally Sox2) outperforms the simultaneous one. Surprisingly, the sequential protocol activates an early epithelial-to-mesenchymal transition (EMT) as indicated by the upregulation of Slug and N-cadherin followed by a delayed mesenchymal-to-epithelial transition (MET). An early EMT induced by 1.5-day TGF-ß treatment enhances reprogramming with the simultaneous protocol, whereas 12-day treatment blocks reprogramming. Consistent results were obtained when the TGF-ß antagonist Repsox was applied in the sequential protocol. These results reveal a time-sensitive role of individual factors for optimal reprogramming and a sequential EMT-MET mechanism at the start of reprogramming. Our studies provide a rationale for further optimizing reprogramming, and introduce the concept of a sequential EMT-MET mechanism for cell fate decision that should be investigated further in other systems, both in vitro and in vivo.

How to compute which genes control drug resistance dynamics.

Increasing evidence shows that genes have a pivotal role in affecting the dynamic pattern of viral loads in the body of a host. By reviewing the biochemical interactions between a virus and host cells as a dynamic system, we outline a computational approach for mapping the genetic control of virus dynamics. The approach integrates differential equations (DEs) to quantify the dynamic origin and behavior of a viral infection system. It enables geneticists to generate various testable hypotheses about the genetic control mechanisms for virus dynamics and infection. The experiment designed according to this approach will also enable researchers to gain insight into the role of genes in limiting virus abundance and the dynamics of viral drug resistance, facilitating the development of personalized medicines to eliminate viral infections.

A statistical design for testing transgenerational genomic imprinting in natural human populations.

Genomic imprinting is a phenomenon in which the same allele is expressed differently, depending on its parental origin. Such a phenomenon, also called the parent-of-origin effect, has been recognized to play a pivotal role in embryological development and pathogenesis in many species. Here we propose a statistical design for detecting imprinted loci that control quantitative traits based on a random set of three-generation families from a natural population in humans. This design provides a pathway for characterizing the effects of imprinted genes on a complex trait or disease at different generations and testing transgenerational changes of imprinted effects. The design is integrated with population and cytogenetic principles of gene segregation and transmission from a previous generation to next. The implementation of the EM algorithm within the design framework leads to the estimation of genetic parameters that define imprinted effects. A simulation study is used to investigate the statistical properties of the model and validate its utilization. This new design, coupled with increasingly used genome-wide association studies, should have an immediate implication for studying the genetic architecture of complex traits in humans.

The C-terminal pentapeptide of Nanog tryptophan repeat domain interacts with Nac1 and regulates stem cell proliferation but not pluripotency.

Overexpression of Nanog in mouse embryonic stem (ES) cells has been shown to abrogate the requirement of leukemia inhibitory factor for self-renewal in culture. Little is known about the molecular mechanism of Nanog function. Here we describe the role of the tryptophan repeat (WR) domain, one of the two transactivators at its C terminus, in regulating stem cell proliferation as well as pluripotency. We first created a supertransactivator, W2W3x10, by duplicating repeats W2W3 10 times and discovered that it can functionally substitute for wild type WR at sustaining pluripotency, albeit with a significantly slower cell cycle, phenocopying Nanog(9W) with the C-terminal pentapeptide (WNAAP) of WR deleted. ES cells carrying both W2W3x10 and Nanog(9W) have a longer G1 phase, a shorter S phase in cell cycle distribution and progression analysis, and a lower level of pAkt(Ser473) compared with wild type Nanog, suggesting that both mutants impact the cell cycle machinery via the phosphatidylinositol 3-kinase/Akt pathway. Both mutants remain competent in dimerizing with Nanog but cannot form a complex with Nac1 efficiently, suggesting that WNAAP may be involved in Nac1 binding. By tagging Gal4DBD with WNAAP, we demonstrated that this pentapeptide is sufficient to confer Nac1 binding. Furthermore, we can rescue W2W3x10 by placing WNAAP at the corresponding locations. Finally, we found that Nanog and Nac1 synergistically up-regulate ERas expression and promote the proliferation of ES cells. These results suggest that Nanog interacts with Nac1 through WNAAP to regulate the cell cycle of ES cells via the ERas/phosphatidylinositol 3-kinase/Akt pathway, but not pluripotency, thus decoupling cell cycle control from pluripotency.

Sorting nexin 33 induces mammalian cell micronucleated phenotype and actin polymerization by interacting with Wiskott-Aldrich syndrome protein.

Sorting nexin 33 (SNX33) is a novel member of the sorting nexin superfamily with three predicted structural domains, SH3-PX-BAR. Very little is known about the cellular function of SNX33. In an effort to analyze its structure/function relationship, we attempted but failed to generate stable cell lines for short hairpin RNA or overexpression SNX33. Transient knockdown of SNX33 induces both HeLa and MCF7 cells to grow multiple long processes, delay the G(1)/M transition, and become more apoptotic, implying that SNX33 may control cell cycle process through influence the cytoskeleton. In vitro cell lineage analysis revealed that cells transfected with SNX33 failed to divide and became micronucleated, suggesting a specific defect in cytokinesis. Further analysis revealed that SNX33 induced the accumulation of actin at the perinuclear space, which might have disabled the cytokinetic machinery. However, SNX33 appears to mediate actin polymerization indirectly, as they do not interact with each other. SNX33 interacts with itself and SNX9. Interestingly, it also interacts with VCA domain of Wiskott-Aldrich syndrome protein (WASp), a protein known to be involved in actin polymerization. Indeed, cells overexpressing WASp failed to divide and form stable colonies as SNX33, consistent with the notion that SNX33 may interfere with cytokinesis. On the other hand, knockdown of WASp alleviates the phenotype induced by SNX33. Taken together, our results suggest that SNX33 plays a role in maintaining cell shape and cell cycle progression through its interaction with WASp.

Regulation of the pluripotency marker Rex-1 by Nanog and Sox2.

Rex-1 (Zfp-42) is a known marker for undifferentiated embryonic stem cells and teratocarcinoma cells. However, the mechanism by which Rex-1 is regulated in pluripotent cells remains unresolved. Here we report that Nanog, an Nk-2 homeodomain protein known for its role in maintaining stem cell pluripotency, is a transcription activator for the Rex-1 promoter. Knockdown of Nanog in embryonic stem cells resulted in a reduction of Rex-1 expression, whereas forced expression of Nanog in P19 stimulated Rex-1 expression. Employing a Rex-1 reporter, we demonstrate that Nanog transactivates Rex-1 directly. Serial deletion studies mapped the Nanog-responsive element between -187 and -286 of the Rex-1 promoter. Although Oct-3/4 and Sox2 can both transactivate Rex-1 promoter, only Sox2 cooperates with Nanog in up-regulating Rex-1. Furthermore, we demonstrate that the C terminus of Nanog is responsible for transactivating the Rex-1 promoter, a function that can be substituted for by a viral transactivator Vp16 efficiently in NIH3T3 cells but less so in P19 cells. Taking these findings together, we conclude that Rex-1 is a direct target of Nanog, which is augmented by Sox2 and Oct-3/4.