Systematic review of genetic risk score in coronary heart disease and other diseases

In order to better take advantage of the abundant results from large-scale genomic association studies, investigators are turning to a genetic risk score (GRS) method in order to combine the information from common modest-effect risk alleles into an efficient risk assessment statistic. The statistical properties of these GRSs are poorly understood. As a first step toward a better understanding of GRSs, a systematic analysis of recent investigations using a GRS was undertaken. GRS studies were searched in the areas of coronary heart disease (CHD), cancer, and other common diseases using bibliographic databases and by hand-searching reference lists and journals. Twenty-one independent case-control studies, cohort studies, and simulation studies (12 in CHD, 9 in other diseases) were identified. The underlying statistical assumptions of the GRS using the experience of the Framingham risk score were investigated. Improvements in the construction of a GRS guided by the concept of composite indicators are discussed. The GRS will be a promising risk assessment tool to improve prediction and diagnosis of common diseases.^

A systematic analysis of INO80 in DNA replication

ATP-dependent chromatin remodeling has been shown to be critical for transcription and DNA repair. However, the involvement of ATP-dependent chromatin remodeling in DNA replication remains poorly defined. Interestingly, we found that the INO80 chromatin-remodeling complex is directly involved in the DNA damage tolerance pathways activated during DNA replication. DNA damage tolerance is important for genomic stability and is controlled by formation of either mono-ubiquitinated or multi-ubiquitinated PCNA, which respectively induce error prone or error-free replication bypass of the lesions. In addition, homologous recombination (HR) mediated by the Rad51 pathway is also involved in the DNA damage tolerance pathways. ^ We found that INO80 is specifically recruited to replication origins during S phase in a genome-wide fashion. In addition, DNA combing analysis shows INO80 is required for the resumption of replication at stalled forks induced by methyl methane-sulfonate (MMS). Mechanistically, we find that INO80 is required for PCNA ubiquitination as well as for Rad51 mediated processing of replication forks after MMS treatment. Furthermore, chromatin immunoprecipitation at specific ARSs indicates INO80 is necessary for Rad18 and Rad51 recruitment to replication forks after MMS treatment. Moreover, 2D gel analysis shows INO80 is necessary to process Rad51 mediated intermediates at impeded replication forks. ^ In conclusion, our findings establish a novel role of a chromatin-remodeling complex in DNA damage tolerance pathways and suggest that chromatin remodeling is fundamentally important to ensure faithful replication of DNA and genome stability in eukaryotes. ^

Population genetics of VNTR loci and their applications in evolutionary studies

Variable number of tandem repeats (VNTR) are genetic loci at which short sequence motifs are found repeated different numbers of times among chromosomes. To explore the potential utility of VNTR loci in evolutionary studies, I have conducted a series of studies to address the following questions: (1) What are the population genetic properties of these loci? (2) What are the mutational mechanisms of repeat number change at these loci? (3) Can DNA profiles be used to measure the relatedness between a pair of individuals? (4) Can DNA fingerprint be used to measure the relatedness between populations in evolutionary studies? (5) Can microsatellite and short tandem repeat (STR) loci which mutate stepwisely be used in evolutionary analyses?^ A large number of VNTR loci typed in many populations were studied by means of statistical methods developed recently. The results of this work indicate that there is no significant departure from Hardy-Weinberg expectation (HWE) at VNTR loci in most of the human populations examined, and the departure from HWE in some VNTR loci are not solely caused by the presence of population sub-structure.^ A statistical procedure is developed to investigate the mutational mechanisms of VNTR loci by studying the allele frequency distributions of these loci. Comparisons of frequency distribution data on several hundreds VNTR loci with the predictions of two mutation models demonstrated that there are differences among VNTR loci grouped by repeat unit sizes.^ By extending the ITO method, I derived the distribution of the number of shared bands between individuals with any kinship relationship. A maximum likelihood estimation procedure is proposed to estimate the relatedness between individuals from the observed number of shared bands between them.^ It was believed that classical measures of genetic distance are not applicable to analysis of DNA fingerprints which reveal many minisatellite loci simultaneously in the genome, because the information regarding underlying alleles and loci is not available. I proposed a new measure of genetic distance based on band sharing between individuals that is applicable to DNA fingerprint data.^ To address the concern that microsatellite and STR loci may not be useful for evolutionary studies because of the convergent nature of their mutation mechanisms, by a theoretical study as well as by computer simulation, I conclude that the possible bias caused by the convergent mutations can be corrected, and a novel measure of genetic distance that makes the correction is suggested. In summary, I conclude that hypervariable VNTR loci are useful in evolutionary studies of closely related populations or species, especially in the study of human evolution and the history of geographic dispersal of Homo sapiens. (Abstract shortened by UMI.) ^

A cladistic analysis of apolipoprotein B gene variation and plasma lipid and apolipoprotein B levels, using familial data

Any functionally important mutation is embedded in an evolutionary matrix of other mutations. Cladistic analysis, based on this, is a method of investigating gene effects using a haplotype phylogeny to define a set of tests which localize causal mutations to branches of the phylogeny. Previous implementations of cladistic analysis have not addressed the issue of analyzing data from related individuals, though in human studies, family data are usually needed to obtain unambiguous haplotypes. In this study, a method of cladistic analysis is described in which haplotype effects are parameterized in a linear model which accounts for familial correlations. The method was used to study the effect of apolipoprotein (Apo) B gene variation on total-, LDL-, and HDL-cholesterol, triglyceride, and Apo B levels in 121 French families. Five polymorphisms defined Apo B haplotypes: the signal peptide Insertion/deletion, Bsp 1286I, XbaI, MspI, and EcoRI. Eleven haplotypes were found, and a haplotype phylogeny was constructed and used to define a set of tests of haplotype effects on lipid and apo B levels.^ This new method of cladistic analysis, the parametric method, found significant effects for single haplotypes for all variables. For HDL-cholesterol, 3 clusters of evolutionarily-related haplotypes affecting levels were found. Haplotype effects accounted for about 10% of the genetic variance of triglyceride and HDL-cholesterol levels. The results of the parametric method were compared to those of a method of cladistic analysis based on permutational testing. The permutational method detected fewer haplotype effects, even when modified to account for correlations within families. Simulation studies exploring these differences found evidence of systematic errors in the permutational method due to the process by which haplotype groups were selected for testing.^ The applicability of cladistic analysis to human data was shown. The parametric method is suggested as an improvement over the permutational method. This study has identified candidate haplotypes for sequence comparisons in order to locate the functional mutations in the Apo B gene which may influence plasma lipid levels. ^

The association of genetic groups 1, 2, and 3 of Mycobacterium tuberculosis with pulmonary and extra-pulmonary tuberculosis

Objective. To investigate the association of the three major genetic groups of Mycobacterium tuberculosis with pulmonary and extra-pulmonary tuberculosis in clustered and non-clustered TB cases in the Houston area. ^ Study design. Secondary analysis of an ambi-directional study. ^ Study population. Three hundred fifty-eight confirmed cases of tuberculosis in the Houston that occurred between October 1995 and May 1997, who had been interviewed by the Houston T13 Initiative staff at Baylor College of Medicine, and whose isolates have had their DNA fingerprint and genetic group determined. ^ Exclusions. Individuals whose mycobacterial genotype was unknown, or whose data variables were unavailable. ^ Source of data. Laboratory results, patient interviews, and medical records at clinics and hospitals of the study population. ^ Results. In clustered cases, the majority of both, pulmonary and extra-pulmonary TB cases were caused by genetic group 1. Independent factors were assessed to determine the interactions that may influence the site of infection or increase the risk for one site or another. HIV negative males were protected against extra-pulmonary TB compared to HIV negative females. Individuals ages 1–14 years were at higher risk of having extra-pulmonary TB. Group 3 organisms were found less frequently in the total population in general, especially in extra-pulmonary disease. This supports the evidence in previous studies that this group is the least virulent and genetically distinct from the other two groups. Group 1 was found more frequently among African Americans than other ethnic groups, a trend for future investigations. ^ Among the non-clustered cases, group 2 organisms were the majority of the organisms found in both sites. They were also the majority of organisms found in African Americans, Caucasians, and Hispanics causing the majority of the infections at both sites. However, group 1 organisms were the overwhelming majority found in Asian/Pacific Islander individuals, which may indicate these organisms are either endemic to that area, or that there is an ethnic biological factor involved. This may also be due to a systematic bias, since isolates from individuals from that geographic region lack adequate copies of the insertion sequence IS6110, which leads to their placement in the non-clustered population. ^ The three genetic groups of Mycobacterium tuberculosis were not found equally distributed between sites of infection in both clustered and non-clustered cases. Furthermore, these groups were not distributed in the same patterns among the clustered and non-clustered cases, but rather in distinct patterns. ^

Familial aggregation, candidate genes and genome scans: Analyzing the role of genetics in stroke

The genetic etiology of stroke likely reflects the influence of multiple loci with small effects, each modulating different pathophysiological processes. This research project utilized three analytical strategies to address the paucity of information related to the identification and characterization of genetic variation associated with stroke in the general population. ^ First, the general contribution of familial factors to stroke susceptibility was evaluated in a population-based sample of unrelated individuals. Increased risk of subclinical cerebral infarction was observed among individuals with a positive parental history of stroke. This association did not appear to be mediated by established stroke risk factors, specifically blood pressure levels or hypertension status. ^ The need to identify specific gene variation associated with stroke in the general population was addressed by evaluating seven candidate gene polymorphisms in a population-based sample of unrelated individuals. Three polymorphisms were significantly associated with increased subclinical cerebral infarction or incident clinical ischemic stroke risk. These relationships include the G-protein β3 subunit 825C/T polymorphism and clinical stroke in Whites, the lipoprotein lipase S/X447 polymorphism and subclinical and clinical stroke in men, and the angiotensin I-converting enzyme Ins/Del polymorphism and subclinical stroke in White men. These associations did not appear to be obfuscated by the stroke risk factors adjusted for in the analysis models specifically blood pressure levels or anti-hypertensive medication use. ^ The final research strategy considered, on a genome-wide scale, the idea that genetic variation may contribute to the occurrence of hypertension or stroke through a common etiologic pathway. Genomic regions were identified for which significant evidence of heterogeneity was observed among hypertensive sibpairs stratified by family history of stroke information. Regions identified on chromosome 15 in African Americans, and chromosome 13 in Whites and African Americans, suggest the presence of genes influencing hypertension and stroke susceptibility. ^ Insight into the role of genetics in stroke is useful for the potential early identification of individuals at increased risk for stroke and improved understanding of the etiology of the disease. The ultimate goal of these endeavors is to guide the development of therapeutic intervention and informed prevention to provide a lasting and positive impact on public health. ^

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

Genetic epidemiology of childhood brain tumors

A rare familial cancer syndrome involving childhood brain tumors (CBT), breast cancer, sarcomas and an array of other tumors has been described (Li and Fraumeni 1969, 1975, 1982, 1987). A survey of CBT identified through the Connnecticut Tumor Registry in 1984 revealed a high frequency of CBT, leukemia and other childhood cancer in siblings of CBT patients (Farwell and Flannery, 1984). Other syndromes such as neurofibromatosis and nevoid basal cell carcinoma syndrome have also been associated with CBT; however, no systematic family studies have been conducted to determine the extent to which cancer aggregates in family members of CBT patients. This family study was designed to determine the frequency of cancer aggregation overall or at specific sites, to determine the frequency of known or potentially hereditary syndromes in families of CBT patients, and to determine a genetic model to characterize familial cancer syndromes and to identify specific kindreds to which such a model(s) might apply. This study includes 244 confirmed CBT patients referred to the University of Texas M. D. Anderson Cancer Center between the years 1944 and 1983, diagnosed under the age of 15 years and resident in the U.S. or Canada. Family histories were obtained on the proband's first (parents, siblings and offspring) and second degree (proband's aunts, uncles and grandparents) relatives following sequential sampling scheme rules. To determine if cancer aggregates in families, we compared the cancer experience in the population to that expected in the general population using Connecticut Tumor Registry calendar year, age, race and sex-specific rates. The standardized incidence ratio (SIR) for cancer overall was 0.91 (41 observed (O) and 44.94 expected (E); 95% Confidence Interval (CI) = 0.65-1.24). We observed a significant excess of colon cancer among the proband's first degree relatives (O/E = 5/1.64; 95% CI = 1.01-7.65), in particular those under age 45 year. Segregation analysis showed evidence for multifactorial inheritance in the small percentage (N = 5) of the families. ^

Molecular and evolutionary genetics of the X -linked visual pigment genes in humans and New World monkeys

Normal humans have one red and at least one green visual pigment genes. These genes are tightly linked as tandem repeats on the X chromosome and each of them has six exons. There is only one X-linked visual pigment gene in New World monkeys (NWMs) but the locus has three polymorphic alleles encoding red, yellow and green visual pigments, respectively. The spectral properties of the squirrel monkey and the marmoset (both NWMs) have been studied and partial sequences of the three alleles are available. To study the evolutionary history of these X-linked opsin genes in humans and NWMs, coding and intron sequences of the three squirrel monkey alleles and the three marmoset alleles were amplified by PCR followed by subcloning and sequencing. Introns 2 and 4 of the human red and green pigment genes were also sequenced. The results obtained are as follows: (1) The sequences of introns 2 and 4 of the human red and green opsin genes are significantly more similar between the two genes than are coding sequences, contrary to the usual situation where coding regions are better conserved in evolution than are introns. The high similarities in the two introns are probably due to recent gene conversion events during evolution of the human lineage. (2) Phylogenetic analysis of both intron and exon sequences indicates that the phylogenetic tree of the available primate opsin genes is the same as the species tree. The two human genes were derived from a gene duplication event after the divergence of the human and NWM lineages. The three alleles in each of the two NWM species diverged after the split of the two NWMs but have persisted in the population for at least 5 million years. (3) Allelic gene conversion might have occurred between the three squirrel monkey alleles. (4) A model of additive effect of hydroxyl-bearing amino acids on spectral tuning is proposed by treating some unknown variables as groups. Under the assumption that some residues have no effect, it is found that at least five amino acid residues, at positions 178 (3 nm), 180 (5 nm), 230 ($-$4 nm), 277 (9 nm) and 285 (13 nm), have linear spectral tuning effects. (5) Adaptive evolution of the opsin genes to different spectral peaks was observed at four residues that are important for spectral tuning. ^

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

GENETIC ANALYSIS OF THE FUNCTION OF THE DROSOPHILA DOUBLESEX-RELATED FACTOR dmrt93B

DMRT (Doublesex and Mab-3 related transcription factor) proteins generally associated with sexual differentiation in many organisms share a common DNA binding domain and are often expressed in reproductive tissues. Aside from doublesex, which is a central factor in the regulation of sex determination, Drosophila possesses three different dmrt genes that are of unknown function. Because the association with sexual differentiation and reproduction is not universal and some DMRT proteins have been found to play other developmental roles we chose to further characterize one of these Drosophila genes. We carried out genetic analysis of dmrt93B, which was previously found to be expressed sex-specifically in the developing somatic gonad and to affect testis morphogenesis in RNAi knockdowns. In order to disrupt this gene, the GAL4 yeast transcriptional activator followed by a polyadenylation signal was inserted after the dmrt93B start codon and introduced into the genome by homologous recombination. Analysis of the knock-in mutation as well as a small deletion removing all dmrt93B sequence demonstrate that loss of function causes partial lethality at the late pupal stage. Surprisingly, these mutations have no significant effect on gonad formation or male fertility. Analysis of GAL4-driven GFP reporter expression indicates that the dmrt93B promoter activity is highly specific to neurons in the suboesophageal and proventricular ganglion in larva and adult of both sexes suggesting a possible role in digestive tract function. Using the Capillary Feeder (CAFÉ) assay to measure daily food intake we find that reduction in this gene’s function leads to an increase in food consumption. These results suggest dmrt93 plays an important role in the formation or maintenance of neurons that affect feeding and support the idea that dmrt genes may not be restricted to roles in sexual differentiation.

THE CONSEQUENCES OF DISRUPTING THE MDM2-P53 BALANCE IN HEMATOPOIESIS

The bone marrow accommodates hematopoietic stem cells and progenitors. These cells provide an indispensible resource for replenishing the blood constituents throughout an organism’s life. A tissue with such a high turn-over rate mandates intact cycling checkpoint and apoptotic pathways to avoid inappropriate cell proliferation and ultimately the development of leukemias. p53, a major tumor suppressor, is a transcription factor that regulates cell cycle, and induces apoptosis and senescence. Mice inheriting a hypomorphic p53 allele in the absence of Mdm2, a p53 inhibitor, have elevated p53 cell cycle activity and die by postnatal day 13 due to hematopoietic failure. Hematopoiesis progresses normally during embryogenesis until it moves to the bone marrow in late development. Increased oxidative stress in the bone marrow compartment postnatally is the impediment for normal hematopoiesis via activation of p53. p53 in turn stimulates the generation of more reactive oxygen species and depletes bone marrow cellularity. Also, p53 exerts various defects on the hematopoietic niche by increasing mesenchymal lineage populations and their differentiation. Hematopoietic defects are rescued with antioxidants or when cells are cultured at low oxygen levels. Deletion of p16 partially rescues bone marrow cellularity and progenitors via a p53-independent pathway. Thus, although p53 is required to inhibit tumorigenesis, Mdm2 is required to control ROS-induced p53 levels for sustainable hematopoiesis and survival during homeostasis.

CELL POLARITY REGULATES ORGAN GROWTH THROUGH THE HIPPO PATHWAY

Defects in apical-basal cell polarity and abnormal expression of cell polarity determinants are linked to human cancer. Loss of polarity is highly correlated with malignancy. In Drosophila, perturbation of apical-basal polarity, including overexpressing the apical determinant Crumbs, can lead to uncontrolled tissue growth. Cells mutant for the basolateral determinant scribble overproliferate and can form neoplastic tumors. Interestingly, scribble mutant clones that arise in wild-type tissues are eliminated and therefore do not manifest their tumorigenic potential. However, the mechanisms by which cell polarity coordinates with growth control pathways in developing organs to achieve appropriate organ size remain obscure. To investigate the function of apical determinants in growth regulation, I investigated the mechanism by which the apical determinant Crumbs affects growth in Drosophila imaginal discs. I found that crumbs gain and loss of function cause overgrowth and induction of Hippo target genes. In addition, Crumbs is required for the proper localization of Expanded, an upstream component of the Hippo pathway. Furthermore, we uncoupled the cell polarity and growth control function of Crb through structure-functional analysis. Taken together, our data identify a role of Crb in growth regulation specifically through modulation of the Hippo pathway. To further explore the role of polarity in growth control, I investigated how cells mutant for basolateral determinants are eliminated by using patches of cells mutant for scribble (scribble mutant clones) as a model system. We found that competitive cell-cell interactions eliminate tumorigenic scribble cells by modulation of the Hippo pathway. The regulation of Hippo signaling is required and sufficient to restrain the tumorous growth of scribble mutant cells. Artificially increasing the relative fitness of scribble mutant cells unleashes their tumorigenic potential. Therefore, we have identified a novel tumor-suppression mechanism that depends on signaling between normal and tumorigenic cells. These data identify evasion of cell competition as a critical step toward malignancy and illustrate a role for wild-type tissue in eliminating abnormal cells and preventing the formation of tumors.

Gene Discovery in Nonsyndromic Cleft Lip with or without Cleft Palate

Nonsyndromic cleft lip with or without cleft palate (NSCLP), a common, complex orofacial birth defect that affects approximately 4,000 newborns each year in the United States, is caused by both genetic and environmental factors. Orofacial clefts affect the mouth and nose, causing severe deformity of the face, which require medical, dental and speech therapies. Despite having substantial genetic liability, less than 25% of the genetic contribute to NSCLP has been identified. The studies described in this thesis were performed to identify genes that contribute to NSCLP and to demonstrate the role of these genes in normal craniofacial development. Using genome scan and candidate gene approaches, novel associations with NSCLP were identified. These include MYH9 (7 SNPs, 0.009≤p<0.05), Wnt3A (4 SNPs, 0.001≤p≤0.005), Wnt11 (2 SNPs, 0.001≤p≤0.01) and CRISPLD2 (4 SNPs, 0.001≤p<0.05). The most interesting findings were for CRISPLD2. This gene is expressed in the fused mouse palate at E17.5. In zebrafish, crispld2 localized to the craniofacial region by one day post fertilization. Morpholino knockdown of crispld2 resulted in a lower survival rates and altered neural crest cell (NCC) clustering. Because NCCs form the tissues that populate the craniofacies, this NCC abnormality resulted in cartilage abnormalities of the jaw including fewer ceratobranchial cartilages forming the lower jaw (three pairs compared to five) and broader craniofacies compared to wild-type zebrafish. These findings suggest that the CRISPLD2 gene plays an important role in normal craniofacial development and perturbation of this gene in humans contributes to orofacial clefting. Overall, these results are important because they contribute to our understanding of normal craniofacial development and orofacial clefting etiology, information that can be used to develop better methods to diagnose, counsel and potentially treat NSCLP patients.