IDENTIFYING GENETIC VARIANTS AND CHARACTERIZING THEIR ROLE IN CLUBFOOT

Clubfoot is a common, complex birth defect affecting 4,000 newborns in the United States and 135,000 world-wide each year. The clubfoot deformity is characterized by inward and rigid downward displacement of one or both feet, along with persistent calf muscle hypoplasia. Despite strong evidence for a genetic liability, there is a limited understanding of the genetic and environmental factors contributing to the etiology of clubfoot. The studies described in this dissertation were performed to identify variants and/or genes associated with clubfoot. Genome-wide linkage scan performed on ten multiplex clubfoot families identified seven new chromosomal regions that provide new areas to search for clubfoot genes. Troponin C (TNNC2) the strongest candidate gene, located in 20q12-q13.11, is involved in muscle contraction. Exon sequencing of TNNC2 did not identify any novel coding variants. Interrogation of fifteen muscle contraction genes found strong associations with SNPs located in potential regulatory regions of TPM1 (rs4075583 and rs3805965), TPM2 (rs2025126 and rs2145925) and TNNC2 (rs383112 and rs437122). In previous studies, a strong association was found with rs3801776 located in the basal promoter of HOXA9, a gene also involved in muscle development and patterning. Altogether, this data suggests that SNPs located in potential regulatory regions of genes involved in muscle development and function could alter transcription factor binding leading to changes in gene expression. Functional analysis of 3801776/HOXA9, rs2025126/TPM2 and rs2145925/TPM2 showed altered protein binding, which significantly influenced promoter activity. Although the ancestral allele (G) of rs4075583/TPM1 creates a DNA-protein complex, it did not affect TPM1 promoter activity. However and importantly, in the context of a haplotype, rs4075583/G significantly decreased TPM1 promoter activity. These results suggest dysregulation of multiple skeletal muscle genes, TPM1, TPM2, TNNC2 and HOXA9, working in concert may contribute to clubfoot. However, specific allelic combinations involving these four regulatory SNPs did not confer a significantly higher risk for clubfoot. Other combinations of these variants are being evaluated. Moreover, these variants may interact with yet to be discovered variants in other genes to confer a higher clubfoot risk. Collectively, we show novel evidence for the role of skeletal muscle genes in clubfoot indicating that there are multiple genetic factors contributing to this complex birth defect.

Molecular and genetic analysis of Myxococcus xanthus early *development

To identify more mutations that can affect the early development of Myxococcus xanthus, the synthetic transposon TnT41 was designed and constructed. By virtue of its special features, it can greatly facilitate the processes of mutation screening/selection, mapping, cloning and DNA sequencing. In addition, it allows for the systematic discovery of genes in regulatory hierarchies using their target promoters. In this study, the minimal regulatory region of the early developmentally regulated gene 4521 was used as a reporter in the TnT41 mutagenesis. Both positive (P) mutations and negative (N) mutations were isolated based on their effects on 4521 expression.^ Four of these mutations, i.e. N1, N2, P52 and P54 were analyzed in detail. Mutations N1 and N2 are insertion mutations in a gene designated sasB. The sasB gene is also identified in this study by genetic and molecular analysis of five UV-generated 4521 suppressor mutations. The sasB gene encodes a protein without meaningful homology in the databases. The sasB gene negatively regulates 4521 expression possibly through the SasS-SasR two component system. A wild-type sasB gene is required for normal M. xanthus fruiting body formation and sporulation.^ Cloning and sequencing analysis of the P52 mutation led to the identification of an operon that encodes the M. xanthus high-affinity branched-chain amino acid transporter system. This liv operon consists of five genes designated livK, livH, livM, livC, and livF, respectively. The Liv proteins are highly similar to their counterparts from other bacteria in both amino acid sequences, functional motifs and predicted secondary structures. This system is required for development since liv null mutations cause abnormality in fruiting body formation and a 100-fold decrease in sporulation efficiency.^ Mutation P54 is a TnT41 insertion in the sscM gene of the ssc chemotaxis system, which has been independently identified by Dr. Shi's lab. The sscM gene encodes a MCP (methyl-accepting chemotaxis protein) homologue. The SscM protein is predicted to contain two transmembrane domains, a signaling domain and at least one putative methylation site. Null mutations of this gene abolish the aggregation of starving cells at a very early stage, though the sporulation levels of the mutant can reach 10% that of wild-type cells. ^

Expression and function of $\beta$-tubulin isotypes in Chinese hamster ovary (CHO) cells

The availability of isotype specific antisera for $\beta$-tubulin, coupled with genetic and biochemical analysis, has allowed the determination of $\beta$-tubulin isotype expression and distribution in Chinese hamster ovary (CHO) cells. Using genetic manipulations involving selection for colcemid resistance followed by reversion and reselection for drug resistance, we have succeeded in isolating cell lines that exhibit three major and one minor $\beta$-tubulin spots by two-dimensional gel electrophoresis. In concert with isotype specific antibodies, analysis of these mutants demonstrates that CHO cells express two copies of isotype I, at least one copy of isotype IV, and very small amounts of isotype V. Their stoichiometry is approximately 1:1:0.7:0.2. All three isotypes assemble into both cytoplasmic and spindle microtubules, and are similar in their responses to cold, colcemid, and calcium induced depolymerization. They have comparable turnover rates and are equally sensitive to depression of synthesis upon colchicine treatment. These results suggest that $\beta$-tubulin isotypes are used interchangeably to assemble microtubule structures in CHO cells. However, of 18 colcemid resistant mutants with a demonstrable alteration in $\beta$-tubulin, all were found to have the alteration in isotype I, thus leaving open the possibility that subtle differences in isotype properties may exist. Under various conditions of the cell growth, the relative proportion of each expressed isotype does not significantly seem to change except in the early G1 phase of the cell cycle. At this time the synthesis of isotype V increases more than two fold relative to isotype I and IV, while at the same time, total $\beta$-tubulin synthesis is decreased about 60-70%. ^

Genetic predisposition to prostate cancer

Prostate cancer (PC) is a significant economic and health burden in the U.S. and Europe but its causes are largely unknown. The most significant risk factors (after gender) are age and family history of the disease. A gene with high penetrance but low frequency on chromosome 1q, HPC 1, has been suggested to cause a proportion of the familial aggregation of PC but other more common genes, conferring less risk, are also thought to contribute to disease predisposition. We have pursued a strategy to study both types of genetic risk in PC. To identify high penetrance genes, affected men from thirteen families have been genotyped for genetic linkage analysis at six microsatellite markers spanning 45 cM of 1q24-25. Both LOD score and non-parametric statistics provide no significant support for HPC1 in this genomic region, although 3 of the families did combine to produce a LOD score of 0.9. These families will be included in a genome wide search for other PC predisposition genes as part of a multinational collaboration.^ For study of common genetic factors in PC development, leukocyte DNA samples from an unselected series of 55 patients and 67 controls have been examined for genetic differences in two other candidate genes, the androgen receptor gene, hAR, at Xq11-12, and the vitamin D receptor gene, hVDR, at 12q12-14. hAR was typed for two trinucleotide repeat length polymorphisms, (CAG)$\rm\sb{n}$ and (GGC)$\rm\sb{n},$ encoding polyglutamine and polyglycine tracts, respectively, which have been implicated in PC susceptibility. These data, combined with similarly processed patients and controls from the U.K. show no consistent association of allele length with PC risk. A novel finding, however, has been a significant association between the number of GGC repeats and the length of time between diagnosis and relapse in stage T1-T4 Caucasian patients irrespective of therapy and age of the patient. Of 49 patients who relapsed out of 108 entering the study, those with 16 or fewer GGC repeats had an average relapse-free-period of 101 (+/$-$7.7) months while for those with more than 16 repeats the period averaged 48 (+/$-$2.9) months, a difference of 2.1 fold or 4.4 years.^ The second gene, hVDR, was genotyped at two polymorphisms, a synonymous C/T substitution in exon 9 identified by differential TaqI enzymatic digestion and a variable length polyA tract in the 3$\sp\prime$ UTR. Although these polymorphisms are in strong linkage disequilibrium only the polyA region showed a possible association with PC risk. Men homozygous for alleles with fewer than 18 A's had an increased risk (OR = 3.0, p = 0.0578) compared to controls. This result is opposite to the findings of others and may either indicate off-setting random errors which together balance out to no significant overall effect or reflect more complex genetic and/or environmental associations.^ Overall, this research suggests that single gene familial predisposition may be less prominent in PC than in other cancers and that the characteristics of PC pathology may be useful in identifying the effects of common genetic factors. ^

Functions and intramolecular interactions of ribosomal RNA in decoding of termination codons

Two regions in the 3$\prime$ domain of 16S rRNA (the RNA of the small ribosomal subunit) have been implicated in decoding of termination codons. Using segment-directed PCR random mutagenesis, I isolated 33 translational suppressor mutations in the 3$\prime$ domain of 16S rRNA. Characterization of the mutations by both genetic and biochemical methods indicated that some of the mutations are defective in UGA-specific peptide chain termination and that others may be defective in peptide chain termination at all termination codons. The studies of the mutations at an internal loop in the non-conserved region of helix 44 also indicated that this structure, in a non-conserved region of 16S rRNA, is involved in both peptide chain termination and assembly of 16S rRNA.^ With a suppressible trpA UAG nonsense mutation, a spontaneously arising translational suppressor mutation was isolated in the rrnB operon cloned into a pBR322-derived plasmid. The mutation caused suppression of UAG at two codon positions in trpA but did not suppress UAA or UGA mutations at the same trpA positions. The specificity of the rRNA suppressor mutation suggests that it may cause a defect in UAG-specific peptide chain termination. The mutation is a single nucleotide deletion (G2484$\Delta$) in helix 89 of 23S rRNA (the large RNA of the large ribosomal subunit). The result indicates a functional interaction between two regions of 23S rRNA. Furthermore, it provides suggestive in vivo evidence for the involvement of the peptidyl-transferase center of 23S rRNA in peptide chain termination. The $\Delta$2484 and A1093/$\Delta$2484 (double) mutations were also observed to alter the decoding specificity of the suppressor tRNA lysT(U70), which has a mutation in its acceptor stem. That result suggests that there is an interaction between the stem-loop region of helix 89 of 23S rRNA and the acceptor stem of tRNA during decoding and that the interaction is important for the decoding specificity of tRNA.^ Using gene manipulation procedures, I have constructed a new expression vector to express and purify the cellular protein factors required for a recently developed, realistic in vitro termination assay. The gene for each protein was cloned into the newly constructed vector in such a way that expression yielded a protein with an N-terminal affinity tag, for specific, rapid purification. The amino terminus was engineered so that, after purification, the unwanted N-terminal tag can be completely removed from the protein by thrombin cleavage, yielding a natural amino acid sequence for each protein. I have cloned the genes for EF-G and all three release factors into this new expression vector and the genes for all the other protein factors into a pCAL-n expression vector. These constructs will allow our laboratory group to quickly and inexpensively purify all the protein factors needed for the new in vitro termination assay. (Abstract shortened by UMI.) ^

Prevalence and determinants of type 2 diabetes among Filipino Americans in Houston metropolitan statistical area, Texas

Diabetes in adults (type 2) has emerged as a world health problem. Prevalence and risk factors have been found to vary in different populations. The wide range of prevalence rates worldwide indicates the importance of genetic and environmental factors in the etiology of the disease. The few available studies suggest that Filipinos are among the higher-risk groups for developing diabetes. This cross-sectional study estimated the overall prevalence rate of type 2 diabetes among Filipino Americans, ages 20–74 years and residents of Houston Metropolitan Statistical Area, Texas, to be 16.1%. The observed high prevalence was associated with age, sex, family history of diabetes, obesity, region of birth; and, in women, gestational diabetes and income. The diabetic Filipino Americans had a higher proportion of parental history of diabetes, medical history of hypertension, and history of smoking; were physically less active, but generally non-obese, compared with the United States diabetic population. ^

Persisting abdominal symptoms after travelers' diarrhea: Is there a genetic predisposition?

Functional gastrointestinal disorders (FGIDs) are defined as ailments of the mid or lower gastrointestinal tract which are not attributable to any discernable anatomic or biochemical defects.1 FGIDs include functional bowel disorders, also known as persisting abdominal symptoms (PAS). Irritable bowel syndrome (IBS) is one of the most common illnesses classified under PAS.2,3 This is the first prospective study that looks at the etiology and pathogenesis of post-infectious PAS in the context of environmental exposure and genetic susceptibility in a cohort of US travelers to Mexico. Our objective was to identify infectious, genetic and environmental factors that predispose to post infectious PAS. ^ Methods. This is a secondary data analysis of a prospective study on a cohort of 704 healthy North American tourists to Cuernavaca, Morelos and Guadalajara, Jalisco in Mexico. The subjects at risk for Travelers' diarrhea were assessed for chronic abdominal symptoms on enrollment and six months after the return to the US. ^ Outcomes. PAS was defined as disturbances of mid and lower gastrointestinal system without any known pathological or radiological abnormalities, or infectious, or metabolic causes. It refers to functional bowel disease, category C of functional gastrointestinal diseases as defined by the Rome II criterion. PAS was sub classified into Irritable bowel syndrome (IBS) and functional abdominal disease (FAD). ^ IBS is defined as recurrent abdominal pain or discomfort present at least 25% and associated with improvement with defecation, change in frequency and form of stool. FAD encompasses other abdominal symptoms of chronic nature that do not meet the criteria for IBS. It includes functional diarrhea, functional constipation, functional bloating: and unspecified bowel symptoms. ^ Results. Among the 704 travelers studied, there were 202 cases of PAS. The PAS cases included 175 cases of FAD and 27 cases of IBS. PAS was more frequent among subjects who developed traveler's diarrhea in Mexico compared to travelers who remained healthy during the short term visit to Mexico (52 vs. 38; OR = 1.8; CI, 1.3–2.5, P < 0.001). A statistically significant difference was noted in the mean age of subjects with PAS compared to healthy controls (28 vs. 34 yrs; OR = 0.97, CI, 0.95–0.98; P < 0.001). Travelers who experienced multiple episodes, a later onset of diarrhea in Mexico and passed greater numbers of unformed stools were more likely to be identified in PAS group at six months. Participants who developed TD caused by enterotoxigenic E.coli in Mexico showed a 2.6 times higher risk of developing FAD (P = 0.003). Infection with Providencia ssp. also demonstrated a greater risk to developing PAS. Subjects who sought treatment for diarrhea while in Mexico also displayed a significantly lower frequency of IBS at six months follow up (OR = 0.30; CI, 0.10–0.80; P = 0.02). ^ Forty six SNPs belonging to 14 genes were studied. Seven SNPs were associated with PAS at 6 months. These included four SNPs from the Caspase Recruitment Domain-Containing Protein 15 gene (CARD15), two SNPs from Surfactant Pulmonary-Associated Protein D gene (SFTPD) and one from Decay-Accelerating Factor For Complement gene (CD55). A genetic risk score (GRS) was composed based on the 7 SNPs that showed significant association with PAS. A 20% greater risk for PAS was noted for every unit increase in GRS. The risk increased by 30% for IBS. The mean GRS was high for IBS (2.2) and PAS (1.1) compared to healthy controls (0.51). These data suggests a role for these genetic polymorphisms in defining the susceptibility to PAS. ^ Conclusions. The study allows us to identify individuals at risk for developing post infectious IBS (PI-IBS) and persisting abdominal symptoms after an episode of TD. The observations in this study will be of use in developing measures to prevent and treat post-infectious irritable bowel syndrome among travelers including pre-travel counseling, the use of vaccines, antibiotic prophylaxis or the initiation of early antimicrobial therapy. This study also provides insights into the pathogenesis of post infectious PAS and IBS. (Abstract shortened by UMI.)^

Relationship between dietary patterns and body mass index in Hispanic children ages 4--19 years of age

Background. In over 30 years, the prevalence of overweight for children and adolescents has increased across the United States (Barlow et al., 2007; Ogden, Flegal, Carroll, & Johnson, 2002). Childhood obesity is linked with adverse physiological and psychological issues in youth and affects ethnic/minority populations in disproportionate rates (Barlow et al., 2007; Butte et al., 2006; Butte, Cai, Cole, Wilson, Fisher, Zakeri, Ellis, & Comuzzie, 2007). More importantly, overweight in children and youth tends to track into adulthood (McNaughton, Ball, Mishra, & Crawford, 2008; Ogden et al., 2002). Childhood obesity affects body functions such as the cardiovascular, respiratory, gastrointestinal, and endocrine systems, including emotional health (Barlow et al., 2007, Ogden et al., 2002). Several dietary factors have been associated with the development of obesity in children; however, these factors have not been fully elucidated, especially in ethnic/minority children. In particular, few studies have been done to determine the effects of different meal patterns on the development of obesity in children. Purpose. The purpose of the study is to examine the relationships between daily proportions of energy consumed and energy derived from fat across breakfast, lunch, dinner, and snack, and obesity among Hispanic children and adolescents. Methods. A cross-sectional design was used to evaluate the relationship between dietary patterns and overweight status in Hispanic children and adolescents 4-19 years of age who participated in the Viva La Familia Study. The goal of the Viva La Familia Study was to evaluate genetic and environmental factors affecting childhood obesity and its co-morbidities in the Hispanic population (Butte et al., 2006, 2007). The study enrolled 1030 Hispanic children and adolescents from 319 families and examined factors related to increased body weight by focusing on a multilevel analysis of extensive sociodemographic, genetic, metabolic, and behavioral data. Baseline dietary intakes of the children were collected using 24-hour recalls, and body mass index was calculated from measured height and weight, and classified using the CDC standards. Dietary data were analyzed using a GEE population-averaged panel-data model with a cluster variable family identifier to include possible correlations within related data sets. A linear regression model was used to analyze associations of dietary patterns using possible covariates, and to examine the percentage of daily energy coming from breakfast, lunch, dinner, and snack while adjusting for age, sex, and BMI z-score. Random-effects logistic regression models were used to determine the relationship of the dietary variables with obesity status and to understand if the percent energy intake (%EI) derived from fat from all meals (breakfast, lunch, dinner, and snacks) affected obesity. Results. Older children (age 4-19 years) consumed a higher percent of energy at lunch and dinner and less percent energy from snacks compared to younger children. Age was significantly associated with percentage of total energy intake (%TEI) for lunch, as well as dinner, while no association was found by gender. Percent of energy consumed from dinner significantly differed by obesity status, with obese children consuming more energy at dinner (p = 0.03), but no associations were found between percent energy from fat and obesity across all meals. Conclusions. Information from this study can be used to develop interventions that target dietary intake patterns in obesity prevention programs for Hispanic children and adolescents. In particular, intervention programs for children should target dietary patterns with energy intake that is spread throughout the day and earlier in the day. These results indicate that a longitudinal study should be used to further explore the relationship of dietary patterns and BMI in this and other populations (Dubois et al., 2008; Rodriquez & Moreno, 2006; Thompson et al., 2005; Wilson et al., in review, 2008). ^

Identification and characterization of virulence determinants in the Lyme disease spirochete Borrelia burgdorferi

Borrelia burgdorferi is the etiological agent of Lyme disease, the most common tick-borne disease in the United States. Although the most frequently reported symptom is arthritis, patients can also experience severe cardiac, neurologic, and dermatologic abnormalities. The identification of virulence determinants in infectious B. burgdorferi strains has been limited by their slow growth rate, poor transformability, and general lack of genetic tools. The present study demonstrates the use of transposon mutagenesis for the identification of infectivity-related factors in infectious B. burgdorferi, examines the potential role for chemotaxis in mammalian infection, and describes the development of a novel method for the analysis of recombination events at the Ids antigenic variation locus. A pool of Himar1 mutants was isolated using an infectious B. burgdorferi clone and the transposon vector pMarGent. Clones exhibiting reduced infectivity in mice possessed insertions in virulence determinants putatively involved in host survival and dissemination. These results demonstrated the feasibility of extensive transposon mutagenesis studies for the identification of additional infectivity-related factors. mcp-5 mutants were chosen for further study to determine the role of chemotaxis during infection. Animal studies indicated that mcp-5 mutants exhibited a reduced infectivity potential, and suggested a role for mcp-5 during the early stages of infection. An in vitro phenotype for an mcp-5 mutant was not detected. Genetic complementation of an mcp-5 mutant resulted in restoration of Mcp-5 expression in the complemented clone, as demonstrated by western blotting, but the organisms were not infectious in mice. We believe this result is a consequence of differences in expression between genes located on the linear chromosome and genes present on the circular plasmid used for trans-complementation. Overall, this work implicates mcp-5 as an important determinant of mammalian infectivity. Finally, the development of a computer-assisted method for the analysis of recombination events occurring at the B. burgdorferi vls antigenic variation locus has proven highly valuable for the detailed examination of vls gene conversion. The studies described here provide evidence for the importance of chemotaxis during infection in mice and demonstrate advances in both genetic and computational approaches for the further characterization of the Lyme disease spirochete. ^

Caspase 8: Mediating the effects of a novel proteasome inhibitor, NPI-0052

Targeting the proteasome with the sole FDA approved proteasome inhibitor (PI), bortezomib, has been fruitful in specific cancers. Its success has generated an interest in next-generation PIs that might have a therapeutic advantage in cancers, such as leukemia, where bortezomib monotherapy was less effective. This study focuses on a novel, clinically relevant PI, NPI-0052. Experiments show that NPI-0052 targets chymotrypsin- and caspase-like activities more potently than the trypsin-like activity in leukemia cells. NPI-0052 induced apoptosis, as determined by caspase-3 activation and DNA fragmentation. Using caspase inhibitors and caspase-8 (I9.2) or FADD (I2.1) deficient cells revealed that caspase-8 was essential for NPI-0052-induced apoptosis. NPI-0052 killed cells via a caspase-8-tBid-mitochondrial pathway, relying on caspase-8, whereas bortezomib relies on several caspases. NPI-0052 increased reactive oxygen species (ROS) levels, which contributed towards cytotoxicity since an antioxidant conferred protection. To improve the clinical efficacy of PIs, NPI-0052 was combined with epigenetic anti-cancer agents, histone deacetylase inhibitors (HDACi). NPI-0052 with MS-275 or vorinostat (FDA approved HDACi), synergistically induced apoptosis more effectively than an HDACi/bortezomib regimen in Jurkat cells. Caspase-8 and ROS contributed towards NPI-0052/HDACi cytotoxicity and caspase-8 mediated superoxide production by NPI-0052 or NPI-0052/HDACi. The proximal targets of these agents: proteasome activity and histone acetylation were examined to determine if they contributed towards synergistic effects. HDACi targeted proteasomal β subunits and corresponding catalytic activities responsible for degrading proteins. Immunoblotting showed increases in histone-H3 expression and its acetylation with NPI-0052 or NPI-0052/HDACi in Jurkat and primary cells. Importantly, the hyper-acetylation by NPI-0052 was not detected with bortezomib, suggesting that this effect may be unique to NPI-0052. An antioxidant attenuated histone-H3 expression and acetylation induced by NPI-0052 alone or with HDACi. Furthermore, the hyper-acetylation by NPI-0052 relied on caspase-8. These novel results show that a PI is eliciting classical epigenetic alterations, demonstrated by hyper-acetylation of histone-H3. This alteration was oxidant and caspase-8 dependent. Overall, results reveal that caspase-8 mediates many effects induced by NPI-0052. Data show overlapping activities by NPI-0052 and HDACi which are contributing, along with caspase-8 activation and oxidative stress, to cytotoxic interactions in leukemia cells, reinforcing the potential clinical utility of combining these two compounds. ^

A systematic review and analysis of the association between endometriosis and ovarian cancer

Ovarian cancer is the most lethal of the gynecologic malignancies. The development of endometriosis has been shown to increase one's risk of ovarian cancer. Numerous studies have investigated this association, yet none have synthesized the available data. In a pooled analysis of cohort and case-control studies, the association between endometriosis and ovarian cancer was strengthened. Women who developed endometriosis-associated ovarian cancer were more likely to develop an early stage clear cell or endometrioid ovarian cancer histotypes and were more likely to have a better overall prognosis. The prognostic differences between endometriosis-associated ovarian cancer and ovarian cancer without an associated endometriosis may indicate genetic and environmental differences between groups.^

GENETIC INTERACTIONS OF FACTORS THAT REGULATE ALTERNATIVE RNA SPLICING IN THE MALE GERM LINE OF DROSOPHILA

Alternative RNA splicing is a critical process that contributes variety to protein functions, and further controls cell differentiation and normal development. Although it is known that most eukaryotic genes produce multiple transcripts in which splice site selection is regulated, how RNA binding proteins cooperate to activate and repress specific splice sites is still poorly understood. In addition how the regulation of alternative splicing affects germ cell development is also not well known. In this study, Drosophila Transformer 2 (Tra2) was used as a model to explore both the mechanism of its repressive function on its own pre-mRNA splicing, and the effect of the splicing regulation on spermatogenesis in testis. Half-pint (Hfp), a protein known as splicing activator, was identified in an S2 cell-based RNAi screen as a co-repressor that functions in combination with Tra2 in the splicing repression of the M1 intron. Its repressive splicing function is found to be sequence specific and is dependent on both the weak 3’ splice site and an intronic splicing silencer within the M1 intron. In addition we found that in vivo, two forms of Hfp are expressed in a cell type specific manner. These alternative forms differ at their amino terminus affecting the presence of a region with four RS dipeptides. Using assays in Drosophila S2 cells, we determined that the alternative N terminal domain is necessary in repression. This difference is probably due to differential localization of the two isoforms in the nucleus and cytoplasm. Our in vivo studies show that both Hfp and Tra2 are required for normal spermatogenesis and cooperate in repression of M1 splicing in spermatocytes. But interestingly, Tra2 and Hfp antagonize each other’s function in regulating germline specific alternative splicing of Taf1 (TBP associated factor 1). Genetic and cytological studies showed that mutants of Hfp and Taf1 both cause similar defects in meiosis and spermatogenesis. These results suggest Hfp regulates normal spermatogenesis partially through the regulation of taf1 splicing. These observations indicate that Hfp regulates tra2 and taf1 activity and play an important role in germ cell differentiation of male flies.

BAYESIAN STATISTICAL METHODS IN GENE-ENVIRONMENT AND GENE-GENE INTERACTION STUDIES

Complex diseases such as cancer result from multiple genetic changes and environmental exposures. Due to the rapid development of genotyping and sequencing technologies, we are now able to more accurately assess causal effects of many genetic and environmental factors. Genome-wide association studies have been able to localize many causal genetic variants predisposing to certain diseases. However, these studies only explain a small portion of variations in the heritability of diseases. More advanced statistical models are urgently needed to identify and characterize some additional genetic and environmental factors and their interactions, which will enable us to better understand the causes of complex diseases. In the past decade, thanks to the increasing computational capabilities and novel statistical developments, Bayesian methods have been widely applied in the genetics/genomics researches and demonstrating superiority over some regular approaches in certain research areas. Gene-environment and gene-gene interaction studies are among the areas where Bayesian methods may fully exert its functionalities and advantages. This dissertation focuses on developing new Bayesian statistical methods for data analysis with complex gene-environment and gene-gene interactions, as well as extending some existing methods for gene-environment interactions to other related areas. It includes three sections: (1) Deriving the Bayesian variable selection framework for the hierarchical gene-environment and gene-gene interactions; (2) Developing the Bayesian Natural and Orthogonal Interaction (NOIA) models for gene-environment interactions; and (3) extending the applications of two Bayesian statistical methods which were developed for gene-environment interaction studies, to other related types of studies such as adaptive borrowing historical data. We propose a Bayesian hierarchical mixture model framework that allows us to investigate the genetic and environmental effects, gene by gene interactions (epistasis) and gene by environment interactions in the same model. It is well known that, in many practical situations, there exists a natural hierarchical structure between the main effects and interactions in the linear model. Here we propose a model that incorporates this hierarchical structure into the Bayesian mixture model, such that the irrelevant interaction effects can be removed more efficiently, resulting in more robust, parsimonious and powerful models. We evaluate both of the 'strong hierarchical' and 'weak hierarchical' models, which specify that both or one of the main effects between interacting factors must be present for the interactions to be included in the model. The extensive simulation results show that the proposed strong and weak hierarchical mixture models control the proportion of false positive discoveries and yield a powerful approach to identify the predisposing main effects and interactions in the studies with complex gene-environment and gene-gene interactions. We also compare these two models with the 'independent' model that does not impose this hierarchical constraint and observe their superior performances in most of the considered situations. The proposed models are implemented in the real data analysis of gene and environment interactions in the cases of lung cancer and cutaneous melanoma case-control studies. The Bayesian statistical models enjoy the properties of being allowed to incorporate useful prior information in the modeling process. Moreover, the Bayesian mixture model outperforms the multivariate logistic model in terms of the performances on the parameter estimation and variable selection in most cases. Our proposed models hold the hierarchical constraints, that further improve the Bayesian mixture model by reducing the proportion of false positive findings among the identified interactions and successfully identifying the reported associations. This is practically appealing for the study of investigating the causal factors from a moderate number of candidate genetic and environmental factors along with a relatively large number of interactions. The natural and orthogonal interaction (NOIA) models of genetic effects have previously been developed to provide an analysis framework, by which the estimates of effects for a quantitative trait are statistically orthogonal regardless of the existence of Hardy-Weinberg Equilibrium (HWE) within loci. Ma et al. (2012) recently developed a NOIA model for the gene-environment interaction studies and have shown the advantages of using the model for detecting the true main effects and interactions, compared with the usual functional model. In this project, we propose a novel Bayesian statistical model that combines the Bayesian hierarchical mixture model with the NOIA statistical model and the usual functional model. The proposed Bayesian NOIA model demonstrates more power at detecting the non-null effects with higher marginal posterior probabilities. Also, we review two Bayesian statistical models (Bayesian empirical shrinkage-type estimator and Bayesian model averaging), which were developed for the gene-environment interaction studies. Inspired by these Bayesian models, we develop two novel statistical methods that are able to handle the related problems such as borrowing data from historical studies. The proposed methods are analogous to the methods for the gene-environment interactions on behalf of the success on balancing the statistical efficiency and bias in a unified model. By extensive simulation studies, we compare the operating characteristics of the proposed models with the existing models including the hierarchical meta-analysis model. The results show that the proposed approaches adaptively borrow the historical data in a data-driven way. These novel models may have a broad range of statistical applications in both of genetic/genomic and clinical studies.

Sources of familial aggregation and co-aggregation of fasting blood glucose and nutritional factors among Mexican -Americans in Starr County, Texas

The purpose of this study was to determine the effects of nutrient intake, genetic factors and common household environmental factors on the aggregation of fasting blood glucose among Mexican-Americans in Starr County, Texas. This study was designed to determine: (a) the proportion of variation of fasting blood glucose concentration explained by unmeasured genetic and common household environmental effects; (b) the degree of familial aggregation of measures of nutrient intake; and (c) the extent to which the familial aggregation of fasting blood glucose is explained by nutrient intake and its aggregation. The method of path analysis was employed to determine these various effects.^ Genes play an important role in fasting blood glucose: Genetic variation was found to explain about 40% of the total variation in fasting blood glucose. Common household environmental effects, on the other hand, explained less than 3% of the variation in fasting blood glucose levels among individuals. Common household effects, however, did have significant effects on measures of nutrient intake, though it explained only about 10% of the total variance in nutrient intake. Finally, there was significant familial aggregation of nutrient intake measures, but their aggregation did not contribute significantly to the familial aggregation of fasting blood glucose. These results imply that similarities among relatives for fasting blood glucose are not due to similarities in nutrient intake among relatives. ^