Functional data analysis approaches for genotype-phenotype association studies from next-generation sequencing

Next-generation DNA sequencing platforms can effectively detect the entire spectrum of genomic variation and is emerging to be a major tool for systematic exploration of the universe of variants and interactions in the entire genome. However, the data produced by next-generation sequencing technologies will suffer from three basic problems: sequence errors, assembly errors, and missing data. Current statistical methods for genetic analysis are well suited for detecting the association of common variants, but are less suitable to rare variants. This raises great challenge for sequence-based genetic studies of complex diseases.^ This research dissertation utilized genome continuum model as a general principle, and stochastic calculus and functional data analysis as tools for developing novel and powerful statistical methods for next generation of association studies of both qualitative and quantitative traits in the context of sequencing data, which finally lead to shifting the paradigm of association analysis from the current locus-by-locus analysis to collectively analyzing genome regions.^ In this project, the functional principal component (FPC) methods coupled with high-dimensional data reduction techniques will be used to develop novel and powerful methods for testing the associations of the entire spectrum of genetic variation within a segment of genome or a gene regardless of whether the variants are common or rare.^ The classical quantitative genetics suffer from high type I error rates and low power for rare variants. To overcome these limitations for resequencing data, this project used functional linear models with scalar response to develop statistics for identifying quantitative trait loci (QTLs) for both common and rare variants. To illustrate their applications, the functional linear models were applied to five quantitative traits in Framingham heart studies. ^ This project proposed a novel concept of gene-gene co-association in which a gene or a genomic region is taken as a unit of association analysis and used stochastic calculus to develop a unified framework for testing the association of multiple genes or genomic regions for both common and rare alleles. The proposed methods were applied to gene-gene co-association analysis of psoriasis in two independent GWAS datasets which led to discovery of networks significantly associated with psoriasis.^

IDENTIFICATION AND CHARACTERIZATION OF DISTINCT POPULATIONS OF CLONOGENIC BONE MARROW STROMAL CELLS CAPABLE OF TRANSFERRING THE HEMATOPOIETIC MICROENVIRONMENT IN VIVO AND SUPPORTING LT-HSCs IN VITRO

Bone marrow (BM) stromal cells are ascribed two key functions, 1) stem cells for non-hematopoietic tissues (MSC) and 2) as components of the hematopoietic stem cell niche. Current approaches studying the stromal cell system in the mouse are complicated by the low yield of clonogenic progenitors (CFU-F). Given the perivascular location of MSC in BM, we developed an alternative methodology to isolate MSC from mBM. An intact ‘plug’ of bone marrow is expelled from bones and enzymatically disaggregated to yield a single cell suspension. The recovery of CFU-F (1917.95+199) reproducibly exceeds that obtained using the standard BM flushing technique (14.32+1.9) by at least 2 orders of magnitude (P<0.001; N = 8) with an accompanying 196-fold enrichment of CFU-F frequency. Purified BM stromal and vascular endothelial cell populations are readily obtained by FACS. A detailed immunophenotypic analysis of lineage depleted BM identified PDGFRαβPOS stromal cell subpopulations distinguished by their expression of CD105. Both subpopulations retained their original phenotype of CD105 expression in culture and demonstrate MSC properties of multi-lineage differentiation and the ability to transfer the hematopoietic microenvironment in vivo. To determine the capacity of either subpopulation to support long-term multi-lineage reconstituting HSCs, we fractionated BM stromal cells into either the LinNEGPDGFRαβPOSCD105POS and LINNEGPDGFRαβPOSCD105LOW/- populations and tested their capacity to support LT-HSC by co-culturing each population with either 1 or 10 HSCs for 10 days. Following the 10 day co-culture period, both populations supported transplantable HSCs from 10 cells/well co-cultures demonstrating high levels of donor repopulation with an average of 65+23.6% chimerism from CD105POS co-cultures and 49.3+19.5% chimerism from the CD105NEG co-cultures. However, we observed a significant difference when mice were transplanted with the progeny of a single co-cultured HSC. In these experiments, CD105POS co-cultures (100%) demonstrated long-term multi- lineage reconstitution, while only 4 of 8 mice (50%) from CD105NEG -single HSC co-cultures demonstrated long-term reconstitution, suggesting a more limited expansion of functional stem cells. Taken together, these results demonstrate that the PDGFRαβCD105POS stromal cell subpopulation is distinguished by a unique capacity to support the expansion of long-term reconstituting HSCs in vitro.

Inhibitor of cytochrome c messenger RNA -protein interaction in stimulated rat skeletal muscle

The purpose of this work was to examine the possible mechanisms for the regulation of cytochrome c gene expression in response to increased contractile activity in rat skeletal muscle. The working hypothesis was that increased contractile activity enhances cytochrome c gene expression through a cis-element. A 110% increase in cytochrome c mRNA concentration was observed in tibialis anterior (TA) muscle after 9 days of chronic stimulation. Similar difference (120%) exists between soleus (SO) muscle of higher contractile activity and white vastus lateralis (WV) muscle of lower contractile activity. These results suggest that the endogenous cytochrome c gene expression is regulated by contractile activity. Cytochrome c-reporter genes were injected into skeletal muscles to identify the cis-element that is responsible for the regulation. Although the data was inconclusive, part of it suggested the importance of the 3$\sp\prime$-untranslated region (3$\sp\prime$-UTR) in mediating the response to increased contractile activity.^ RNA gel mobility shift (GMSA) and ultraviolet (UV) cross-linking assays revealed specific RNA-protein interaction in a 50-nucleotide region of the 3$\sp\prime$-UTR in unstimulated TA muscle. Computer analysis predicted a stem-loop structure of 17 nucleotides, which provides a structural basis for RNA-protein interaction. These 17 nucleotides are 100% conserved among rat, mouse and human cytochrome c genes and their 13 pseudogenes, suggesting a functional role for this region. The RNA-protein interaction was significantly less in highly active SO muscle than in inactive WV muscle and was dramatically decreased in stimulated TA muscle due to a protein inhibitor(s) associated with ribosome. It is possible that cytochrome c mRNAs undergoing translation are subject to a compartmentalized regulatory influence.^ The conclusion from these results is that increases in contractile activity induce or activate a protein inhibitor(s) associated with ribosome in rat skeletal muscle. The inhibitor decreases RNA-protein interaction in the 3$\sp\prime$-UTR of cytochrome c mRNA, which may result in increased mRNA stability and/or translation. ^

Cloning and molecular analysis of paraxis: A {\it trans\/}-acting factor required for somite maturation

During vertebrate embryogenesis, cells from the paraxial mesoderm coalesce in a rostral-to-caudal progression to form the somites. Subsequent compartmentalization of the somites yields the sclerotome, myotome and dermatome, which give rise to the axial skeleton, axial musculature, and dermis, respectively. Recently, we cloned a novel basic-Helix-Loop-Helix (bHLH) protein, called scleraxis, which is expressed in the sclerotome, in mesenchymal precursors of bone and cartilage, and in connective tissues. This dissertation focuses on the cloning, expression and functional analysis of a bHLH protein termed paraxis, which is nearly identical to scleraxis within the bHLH region but diverges in both its amino and carboxyl termini. During the process of mouse embryogenesis, paraxis transcripts are first detected at about day 7.5 post coitum within the primitive mesoderm lying posterior to the head and heart primordia. Subsequently, paraxis expression progresses caudally through the paraxial mesoderm, immediately preceding somite formation. Paraxis is expressed at high levels in newly formed somites before the first detectable expression of the myogenic bHLH genes, and as the somite becomes compartmentalized, paraxis becomes downregulated within the myotome.^ To determine the function of paraxis during mammalian embryogenesis, mice were generated with a null mutation in the paraxis locus. Paraxis null mice survived until birth, but exhibited severe foreshortening along the anteroposterior axis due to the absence of vertebrae caudal to the midthoracic region. The phenotype also included axial skeletal defects, particularly shortened bifurcated ribs which were detached from the vertebral column, fused vertebrae and extensive truncation and disorganization caudal to the hindlimbs. Mutant neonates also lacked normal levels of trunk muscle and exhibited defects in the dermis as well as the stratification of the epidermis. Analysis of paraxis -/- mutant embryos has revealed a failure of the somites to both properly epithelialize and compartmentalize, resulting in defects in somite-derived cell lineages. These results suggest that paraxis is an essential component of the genetic pathway regulating somitogenesis. ^