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.^

DEVELOPMENT OF NOVEL METHODS TO MINIMIZE THE IMPACT OF SEQUENCING ERRORS IN THE NEXT-GENERATION SEQUENCING DATA ANALYSIS

Next-generation sequencing (NGS) technology has become a prominent tool in biological and biomedical research. However, NGS data analysis, such as de novo assembly, mapping and variants detection is far from maturity, and the high sequencing error-rate is one of the major problems. . To minimize the impact of sequencing errors, we developed a highly robust and efficient method, MTM, to correct the errors in NGS reads. We demonstrated the effectiveness of MTM on both single-cell data with highly non-uniform coverage and normal data with uniformly high coverage, reflecting that MTM’s performance does not rely on the coverage of the sequencing reads. MTM was also compared with Hammer and Quake, the best methods for correcting non-uniform and uniform data respectively. For non-uniform data, MTM outperformed both Hammer and Quake. For uniform data, MTM showed better performance than Quake and comparable results to Hammer. By making better error correction with MTM, the quality of downstream analysis, such as mapping and SNP detection, was improved. SNP calling is a major application of NGS technologies. However, the existence of sequencing errors complicates this process, especially for the low coverage (

The use of second-generation antipsychotics and the changes in physical growth in children and adolescents with perinatally acquired HIV.

Second-generation antipsychotics (SGAs) are increasingly prescribed to treat psychiatric symptoms in pediatric patients infected with HIV. We examined the relationship between prescribed SGAs and physical growth in a cohort of youth with perinatally acquired HIV-1 infection. Pediatric AIDS Clinical Trials Group (PACTG), Protocol 219C (P219C), a multicenter, longitudinal observational study of children and adolescents perinatally exposed to HIV, was conducted from September 2000 until May 2007. The analysis included P219C participants who were perinatally HIV-infected, 3-18 years old, prescribed first SGA for at least 1 month, and had available baseline data prior to starting first SGA. Each participant prescribed an SGA was matched (based on gender, age, Tanner stage, baseline body mass index [BMI] z score) with 1-3 controls without antipsychotic prescriptions. The main outcomes were short-term (approximately 6 months) and long-term (approximately 2 years) changes in BMI z scores from baseline. There were 236 participants in the short-term and 198 in the long-term analysis. In linear regression models, youth with SGA prescriptions had increased BMI z scores relative to youth without antipsychotic prescriptions, for all SGAs (short-term increase = 0.192, p = 0.003; long-term increase = 0.350, p < 0.001), and for risperidone alone (short-term = 0.239, p = 0.002; long-term = 0.360, p = 0.001). Participants receiving both protease inhibitors (PIs) and SGAs showed especially large increases. These findings suggest that growth should be carefully monitored in youth with perinatally acquired HIV who are prescribed SGAs. Future research should investigate the interaction between PIs and SGAs in children and adolescents with perinatally acquired HIV infection.

A molecular analysis of the conditional defect in MuSVts 110 viral RNA splicing

Cells infected with MuSVts110 express a viral RNA which contains an inherent conditional defect in RNA splicing. It has been shown previously that splicing of the MuSVts110 primary transcript is essential to morphological transformation of 6m2 cells in vitro. A growth temperature of 33$\sp\circ$C is permissive for viral RNA splicing,and, consequently, 6m2 cells appear morphologically transformed at this temperature. However, 6m2 cells appear phenotypically normal when incubated at 39$\sp\circ$C, the non-permissive temperature for viral RNA splicing.^ After a shift from 39$\sp\circ$C to 33$\sp\circ$C, the coordinate splicing of previously synthesized and newly transcribed MuSVts110 RNA was achieved. By S1 nuclease analysis of total RNA isolated at various times, 5$\sp\prime$ splice site cleavage of the MuSVts110 transcript appeared to occur 60 minutes after the shift to 33$\sp\circ$C, and 30 minutes prior to detectable exon ligation. In addition, consistent with the permissive temperatures and the kinetic timeframe of viral RNA splicing after a shift to 33$\sp\circ$C, four temperature sensitive blockades to primer extension were identified 26-75 bases upstream of the 3$\sp\prime$ splice site. These blockades likely reflect four branchpoint sequences utilized in the formation of MuSVts110 lariat splicing-intermediates.^ The 54-5A4 cell line is a spontaneous revertant of 6m2 cells and appears transformed at all growth temperatures. Primer extension sequence analysis has shown that a five base deletion occurred at the 3$\sp\prime$ splice site in MuSVts110 RNA allowing the expression of a viral transforming protein in 54-5A4 in the absence of RNA splicing, whereas in the parental 6m2 cell line, a splicing event is necessary to generate a similar transforming protein. As a consequence of this deletion, splicing cannot occur and the formation of the four MuSVts110 branched-intermediates were not observed at any temperature in 54-5A4 cells. However, 5$\sp\prime$ splice site cleavage was still detected at 33$\sp\circ$C.^ Finally, we have investigated the role of the 1488 bp deletion which occurred in the generation of MuSVts110 in the activation of temperature sensitive viral RNA splicing. This deletion appears solely responsible for splice site activation. Whether intron size is the crucial factor in MuSVts110 RNA splicing or whether inhibitory sequences were removed by the deletion is currently unknown. (Abstract shortened with permission of author.) ^

Genetic analysis of factors regulating mesenchymal cell fates

Formation of cartilage and bone involves sequential processes in which undifferentiated mesenchyme aggregates into primordial condensations which subsequently grow and differentiate, resulting in morphogenesis of the adult skeleton. While much has been learned about the structural molecules which comprise cartilage and bone, little is known about the nuclear factors which regulate chondrogenesis and osteogenesis. MHox is a homeobox-containing gene which is expressed in the mesenchyme of facial, limb, and vertebral skeletal precursors during mouse embryogenesis. MHox expression has been shown to require epithelial-derived signals, suggesting that MHox may regulate the epithelial-mesenchymal interactions required for skeletal organogenesis. To determine the functions of MHox, we generated a loss-of-function mutation in the MHox gene. Mice homozygous for a mutant MHox allele exhibit defects of skeletogenesis, involving the loss or malformation of craniofacial, limb and vertebral skeletal structures. The affected skeletal elements are derived from the cranial neural crest, as well as somitic and lateral mesoderm. Analysis of the mutant phenotype during ontogeny demonstrated a defect in the formation or growth of chondrogenic and osteogenic precursors. These findings provide evidence that MHox regulates the formation of preskeletal condensations from undifferentiated mesenchyme. In addition, generation of mice doubly mutant for the MHox and S8 homeobox genes reveal that these two genes interact to control formation of the limb and craniofacial skeleton. Mice carrying mutant alleles for S8 and MHox exhibit an exaggeration of the craniofacial and limb phenotypes observed in the MHox mutant mouse. Thus, MHox and S8 are components of a combinatorial genetic code controlling generation of the skeleton of the skull and limbs. ^

ANALYSIS OF THE CORYNEBACTERIUM PARVUM-INDUCED DECLINE OF MURINE NATURAL KILLER CELL-MEDIATED CYTOTOXICITY

Treatment of mice with the immunomodulating agent, Corynebacterium parvum (C. parvum), was shown to result in a severe and long-lasting depression of splenic natural killer (NK) cell-mediated cytotoxicity 5-21 days post-inoculation. Because NK cells have been implicated in immunosurveillance against malignancy (due to their spontaneous occurrence and rapid reactivity to a variety of histological types of tumors), as well as in resistance to established tumors, this decreased activity was of particular concern, since this effect is contrary to that which would be considered therapeutically desirable in cancer treatment (i.e. a potentiation of antitumor effector functions, including NK cell activity, would be expected to lead to a more effective destruction of malignant cells). Therefore, an analysis of the mechanism of this decline of splenic NK cell activity in C.parvum treated mice was undertaken.^ From in vitro co-culturing experiments, it was found that low NK-responsive C. parvum splenocytes were capable of reducing the normally high-reactivity of cells from untreated syngeneic mice to YAC-1 lymphoma, suggesting the presence of NK-directed suppressor cells in C. parvum treated animals. This was further supported by the demonstration of normal levels of cytotoxicity in C. parvum splenocyte preparations following Ficoll-Hypaque separation, which coincided with removal of the NK-suppressive capabilities of these cells. The T cell nature of these regulatory cells was indicated by (1) the failure of C. parvum to cause a reduction of NK cell activity, or the generation of NK-directed suppressor cells in T cell-deficient athymic mice, (2) the removal of C. parvum-induced suppression by T cell-depleting fractionation procedures or treatments, and (3) demonstration of suppression of NK cell activity by T cell-enriched C. parvum splenocytes. These studies suggest, therefore, that the eventual reduction of suppression by T cell elimination and/or inhibition, may result in a promotion of the antitumor effectiveness of C. parvum due to the contribution of "freed" NK effector cell activity.^ However, the temporary suppression of NK cell activity induced by C. parvum (reactivity of treated mice returns to normal levels within 28 days after C. parvum injection), may in fact be favorable in some situations, e.g. in bone marrow transplantation cases, since NK cells have been suggested to play a role also in the process of bone marrow graft rejection.^ Therefore, the discriminate use of agents such as C. parvum may allow for the controlled regulation of NK cell activity suggested to be necessary for the optimalization of therapeutic regimens. ^

Generation of an hTERT-restricted, CMV-augmented oncolytic adenovirus capable of tumor-specific transgene expression: Development, characterization, and optimization

Current shortcomings in cancer therapy require the generation of new, broadly applicable, potent, targeted treatments. Here, an adenovirus is engineered to replicate specifically in cells with active human telomerase promotion using a modified hTERT promoter, fused to a CMV promoter element. The virus was also modified to contain a visible reporter transgene, GFP. The virus, Ad/hTC-GFP-E1 was characterized in vitro and demonstrated tumor specific activity both by dose and over time course experiments in a variety of cell lines. In vivo, Ad/hTC-GFP-E1 was affected at suppressing tumor growth and providing a survival benefit without causing any measurable toxicity. To increase the host range of the vector, the fiber region was modified to contain an RGD-motif. The vector, AdRGD/hTC-GFP-E1, was recharacterized in vitro, revealing heightened levels of infectivity and toxicity however maintaining a therapeutic window between cancer and normal cell toxicity. AdRGD/hTC-GFP-E1 was administered in vivo by limb perfusion and was observed to be tumor specific both in expression and replication. To further enhance the efficacy of viral vectors in lung delivery, asthma medications were investigated for their abilities to enhance transgene delivery and expression. A combination of bronchodilators, mast cell inhibitors, and mucolytic agents was devised which demonstrated fold increases in expression in immunocompetent mouse lungs as single agents and more homogenous, intense levels of expression when done in combination of all agents. To characterize the methods in which some cancers are resistant or may become resistant to oncolytic treatments, several small molecule inhibitors of metabolic pathways were applied in combination with oncolytic infection in vitro. SP600125 and PD 98059, respective JNK and ERK inhibitors, successfully suppressed oncolytic toxicity, however did not affect infectivity or transgene expression of Ad/hTC-GFP-E1. JNK and ERK inhibition did significantly suppress viral replication, however, as analyzed by lysate transfer and titration assays. In contrast, SB 203580, an inhibitor for p38, did not demonstrate any protective effects with infected cells. Flow cytometric analysis indicated a possible correlation with G1 arrest and suppressed viral production, however more compounds must be investigated to clarify this observation. ^

INFANT MORTALITY IN THE MEXICAN AMERICAN COMMUNITY: AN ANALYSIS OF DIFFERENTIALS IN INFANT, NEONATAL, AND POSTNEONATAL MORTALITY BY ETHNICITY AND PARENTAL NATIVITY, 1974-75 SINGLE LIVE BIRTH COHORT, HARRIS COUNTY, TEXAS

The paradoxically low infant mortality rates for Mexican Americans in Texas have been attributed to inaccuracies in vital registration and idiosyncracies in Mexican migration in rural areas along the U.S.-Mexico border. This study examined infant (IMR), neonatal (NMR), and postneonatal (PNMR) mortality rates of Mexican Americans in an urban, non-border setting, using linked birth and death records of the 1974-75 single live birth cohort (N = 68,584) in Harris County, Texas, which includes the city of Houston and is reported to have nearly complete birth and death registration. The use of parental nativity with the traditional Spanish surname criterion made it possible to distinguish infants of Mexican-born immigrants from those of Blacks, Anglos, other Hispanics, and later-generation, more Anglicized Mexican Americans. Mortality rates were analyzed by ethnicity, parental nativity, and cause of death, with respect to birth weight, birth order, maternal age, legitimacy status, and time of first prenatal care.^ While overall IMRs showed Spanish surname rates slightly higher than Anglo rates, infants of Mexican-born immigrants had much lower NMRs than did Anglos, even for moderately low birth weight infants. However, among infants under 1500 grams, presumably unable to be discharged home in the neonatal period, Mexican Americans had the highest NMR. The inconsistency suggested unreported deaths for Mexican American low birth weight infants after hospital discharge. The PNMR of infants of Mexican immigrants was also lower than for Anglos, and the usual mortality differentials were reversed: high-risk categories of high birth order, high maternal age, and late/no prenatal care had the lowest PNMRs. Since these groups' characteristics are congruent with those of low-income migrants, the data suggested the possibility of migration losses. Cause of death analysis suggested that prematurity and birth injuries are greater problems than heretofore recognized among Mexican Americans, and that home births and "shoebox burials" may be unrecorded even in an urban setting.^ Caution is advised in the interpretation of infant mortality rates for a Spanish surname population of Mexican origin, even in an urban, non-border area with reportedly excellent birth and death registration. ^

Network-based analysis of genome-wide scans of natural selection

The genomic era brought by recent advances in the next-generation sequencing technology makes the genome-wide scans of natural selection a reality. Currently, almost all the statistical tests and analytical methods for identifying genes under selection was performed on the individual gene basis. Although these methods have the power of identifying gene subject to strong selection, they have limited power in discovering genes targeted by moderate or weak selection forces, which are crucial for understanding the molecular mechanisms of complex phenotypes and diseases. Recent availability and rapid completeness of many gene network and protein-protein interaction databases accompanying the genomic era open the avenues of exploring the possibility of enhancing the power of discovering genes under natural selection. The aim of the thesis is to explore and develop normal mixture model based methods for leveraging gene network information to enhance the power of natural selection target gene discovery. The results show that the developed statistical method, which combines the posterior log odds of the standard normal mixture model and the Guilt-By-Association score of the gene network in a naïve Bayes framework, has the power to discover moderate/weak selection gene which bridges the genes under strong selection and it helps our understanding the biology under complex diseases and related natural selection phenotypes.^

Genetic analysis of Pitx2 in left -right asymmetry and cardiovascular development

Pitx2, a paired-related homeobox gene that is mutated in human Rieger Syndrome, plays a key role in transferring the early asymmetric signals to individual organs. Pitx2 encodes three isoforms, Pitx2a, Pitx2b and Pitx2c. I found that Pitx2c was the Pitx2 isoform for regulating left-right asymmetry in heart, lung and the predominant isoform in guts. Previous studies suggested that the generation of left-right asymmetry within individual organs is an all or none, random event. Phenotypic analysis of various Pitx2 allelic combinations, that encode graded levels of Pitx2c, reveals an organ-intrinsic mechanism for regulating left-right asymmetric morphogenesis based on differential response to Pitx2c levels. The heart needs low Pitx2c levels, while the lungs and duodenum require higher doses of Pitx2c. In addition, the duodenal rotation is under strict control of Pitx2c activity. Left-right asymmetry development for aortic arch arteries involves complex vascular remodeling. Left-sided expression of Pitx2c in these developing vessels implied its potential function in this process. In order to determine if Pitx2c also can regulate the left-right asymmetry of the aortic arch arteries, a Pitx2c-specific loss of function mutation is generated. Although in wild type mice, the direction of the aortic arch is always oriented toward the left side, the directions of the aortic arches in the mutants were randomized, showing that Pitx2c also determined the left-right asymmetry of these vessels. I have further showed that the cardiac neural crest wasn't involved in this vascular remodeling process. In addition, all mutant embryos had Double Outlet Right Ventricle (DORV), a common congenital heart disease. This study provided insight into the mechanism of Pitx2c-mediated late stages of left-right asymmetry development and identified the roles of Pitx2c in regulation of aortic arch remodeling and heart development. ^