  The widely used genetic pleiotropic analysis of multiple phenotypes are often
designed for examining the relationship between common variants and a few
phenotypes. They are not suited for both high dimensional phenotypes and high
dimensional genotype (next-generation sequencing) data. To overcome these
limitations, we develop sparse structural equation models (SEMs) as a general
framework for a new paradigm of genetic analysis of multiple phenotypes. To
incorporate both common and rare variants into the analysis, we extend the
traditional multivariate SEMs to sparse functional SEMs. To deal with high
dimensional phenotype and genotype data, we employ functional data analysis and
the alternative direction methods of multiplier (ADMM) techniques to reduce
data dimension and improve computational efficiency. Using large scale
simulations we showed that the proposed methods have higher power to detect
true causal genetic pleiotropic structure than other existing methods.
Simulations also demonstrate that the gene-based pleiotropic analysis has
higher power than the single variant-based pleiotropic analysis. The proposed
method is applied to exome sequence data from the NHLBI Exome Sequencing
Project (ESP) with 11 phenotypes, which identifies a network with 137 genes
connected to 11 phenotypes and 341 edges. Among them, 114 genes showed
pleiotropic genetic effects and 45 genes were reported to be associated with
phenotypes in the analysis or other cardiovascular disease (CVD) related
phenotypes in the literature.
