IDENTIFICATION OF GENES CONTROLLING BIOLOGICAL PROCESSES AND PATHWAYS THROUGH STATISTICAL ANALYSIS AND NETWORK RECONSTRUCTION

The developmental processes and functions of an organism are controlled by the genes and the proteins that are derived from these genes. The identification of key genes and the reconstruction of gene networks can provide a model to help us understand the regulatory mechanisms for the initiation and progression of biological processes or functional abnormalities (e.g. diseases) in living organisms. In this dissertation, I have developed statistical methods to identify the genes and transcription factors (TFs) involved in biological processes, constructed their regulatory networks, and also evaluated some existing association methods to find robust methods for coexpression analyses. Two kinds of data sets were used for this work: genotype data and gene expression microarray data. On the basis of these data sets, this dissertation has two major parts, together forming six chapters. The first part deals with developing association methods for rare variants using genotype data (chapter 4 and 5). The second part deals with developing and/or evaluating statistical methods to identify genes and TFs involved in biological processes, and construction of their regulatory networks using gene expression data (chapter 2, 3, and 6). For the first part, I have developed two methods to find the groupwise association of rare variants with given diseases or traits. The first method is based on kernel machine learning and can be applied to both quantitative as well as qualitative traits. Simulation results showed that the proposed method has improved power over the existing weighted sum method (WS) in most settings. The second method uses multiple phenotypes to select a few top significant genes. It then finds the association of each gene with each phenotype while controlling the population stratification by adjusting the data for ancestry using principal components. This method was applied to GAW 17 data and was able to find several disease risk genes. For the second part, I have worked on three problems. First problem involved evaluation of eight gene association methods. A very comprehensive comparison of these methods with further analysis clearly demonstrates the distinct and common performance of these eight gene association methods. For the second problem, an algorithm named the bottom-up graphical Gaussian model was developed to identify the TFs that regulate pathway genes and reconstruct their hierarchical regulatory networks. This algorithm has produced very significant results and it is the first report to produce such hierarchical networks for these pathways. The third problem dealt with developing another algorithm called the top-down graphical Gaussian model that identifies the network governed by a specific TF. The network produced by the algorithm is proven to be of very high accuracy.

Parentage analysis with genetic markers in natural populations. I. The expected proportion of offspring with unambiguous paternity.

Recent studies indicate that polymorphic genetic markers are potentially helpful in resolving genealogical relationships among individuals in a natural population. Genetic data provide opportunities for paternity exclusion when genotypic incompatibilities are observed among individuals, and the present investigation examines the resolving power of genetic markers in unambiguous positive determination of paternity. Under the assumption that the mother for each offspring in a population is unambiguously known, an analytical expression for the fraction of males excluded from paternity is derived for the case where males and females may be derived from two different gene pools. This theoretical formulation can also be used to predict the fraction of births for each of which all but one male can be excluded from paternity. We show that even when the average probability of exclusion approaches unity, a substantial fraction of births yield equivocal mother-father-offspring determinations. The number of loci needed to increase the frequency of unambiguous determinations to a high level is beyond the scope of current electrophoretic studies in most species. Applications of this theory to electrophoretic data on Chamaelirium luteum (L.) shows that in 2255 offspring derived from 273 males and 70 females, only 57 triplets could be unequivocally determined with eight polymorphic protein loci, even though the average combined exclusionary power of these loci was 73%. The distribution of potentially compatible male parents, based on multilocus genotypes, was reasonably well predicted from the allele frequency data available for these loci. We demonstrate that genetic paternity analysis in natural populations cannot be reliably based on exclusionary principles alone. In order to measure the reproductive contributions of individuals in natural populations, more elaborate likelihood principles must be deployed.

Morphoproteomic analysis reveals an overexpressed and constitutively activated phospholipase D1-mTORC2 pathway in endometrial carcinoma.

The mammalian target of rapamycin (MTOR) assembles into two distinct complexes: mTOR complex 1 (mTORC1) is predominantly cytoplasmic and highly responsive to rapamycin, whereas mTOR complex 2 (mTORC2) is both cytoplasmic and nuclear, and relatively resistant to rapamycin. mTORC1 and mTORC2 phosphorylatively regulate their respective downstream effectors p70S6K/4EBP1, and Akt. The resulting activated mTOR pathways stimulate protein synthesis, cellular proliferation, and cell survival. Moreover, phospholipase D (PLD) and its product, phosphatidic acid (PA) have been implicated as one of the upstream activators of mTOR signaling. In this study, we investigated the activation status as well as the subcellular distribution of mTOR, and its upstream regulators and downstream effectors in endometrial carcinomas (ECa) and non-neoplastic endometrial control tissue. Our data show that the mTORC2 activity is selectively elevated in endometrial cancers as evidenced by a predominant nuclear localization of the activated form of mTOR (p-mTOR at Ser2448) in malignant epithelium, accompanied by overexpression of nuclear p-Akt (Ser473), as well as overexpression of vascular endothelial growth factor (VEGF)-A isoform, the latter a resultant of target gene activation by mTORC2 signaling via hypoxia-inducible factor (HIF)-2alpha. In addition, expression of PLD1, one of the two major isoforms of PLD in human, is increased in tumor epithelium. In summary, we demonstrate that the PLD1/PA-mTORC2 signal pathway is overactivated in endometrial carcinomas. This suggests that the rapamycin-insensitive mTORC2 pathway plays a major role in endometrial tumorigenesis and that therapies designed to target the phospholipase D pathway and components of the mTORC2 pathway should be efficacious against ECa.

Quantitative analysis of arachidonic acid, endocannabinoids, N-acylethanolamines and steroids in biological samples by LCMS/MS: Fit to purpose.

Arachidonic acid (5Z,8Z,11Z,14Z-eicosatetraenoic acid; C20:4) (arachidonate, AA) is a vital polyunsaturated omega-6 fatty acid (PUFA) without its presence the mammalian brain, muscles, and possibly other organs cannot develop or function [1] and [2]. AA fulfils numerous known and possibly yet unknown functions as integral part of mammalian phospholipid membranes and as free AA which also acts as a precursor of a variety of biologically active lipid mediators generally referred to as eicosanoids (e.g., prostaglandins, leukotrienes). A more recent class of eicosanoids is composed of the endogenous cannabinoids (endocannabinoids) 2-arachidonoyl glycerol (2-AG) and arachidonoyl ethanolamide (anandamide, AEA), which act on cannabinoid CB1 and CB2 receptors but also modulate ion channels and nuclear receptors [3] and [4]. In recent years, the role of endocannabinoids as prominent anti-inflammatory and neuromodulatory eicosanoids has been shown by numerous studies [5].

Can widely used cell type markers predict the suitability of immortalized or primary mammary epithelial cell models?

BACKGROUND Mammary cell cultures are convenient tools for in vitro studies of mammary gland biology. However, the heterogeneity of mammary cell types, e.g., glandular milk secretory epithelial or myoepithelial cells, often complicates the interpretation of cell-based data. The present study was undertaken to determine the relevance of bovine primary mammary epithelial cells isolated from American Holstein (bMECUS) or Swiss Holstein-Friesian (bMECCH) cows, and of primary bovine mammary alveolar epithelial cells stably transfected with simian virus-40 (SV-40) large T-antigen (MAC-T) for in vitro analyses. This was evaluated by testing their expression pattern of cytokeratin (CK) 7, 18, 19, vimentin, and α-smooth muscle actin (α-SMA). RESULTS The expression of the listed markers was assessed using real-time quantitative PCR, flow cytometry and immunofluorescence microscopy. Characteristic markers of the mesenchymal (vimentin), myoepithelial (α-SMA) and glandular secretory cells (CKs) showed differential expression among the studied cell cultures, partly depending on the analytical method used. The relative mRNA expression of vimentin, CK7 and CK19, respectively, was lower (P < 0.05) in immortalized than in primary mammary cell cultures. The stain index (based on flow cytometry) of CK7 and CK19 protein was lower (P < 0.05) in MAC-T than in bMECs, while the expression of α-SMA and CK18 showed an inverse pattern. Immunofluorescence microscopy analysis mostly confirmed the mRNA data, while partly disagreed with flow cytometry data (e.g., vimentin level in MAC-T). The differential expression of CK7 and CK19 allowed discriminating between immortal and primary mammary cultures. CONCLUSIONS The expression of the selected widely used cell type markers in primary and immortalized MEC cells did not allow a clear preference between these two cell models for in vitro analyses studying aspects of milk composition. All tested cell models exhibited to a variable degree epithelial and mesenchymal features. Thus, based on their characterization with widely used cell markers, none of these cultures represent an unequivocal alveolar mammary epithelial cell model. For choosing the appropriate in vitro model additional properties such as the expression profile of specific proteins of interest (e.g., transporter proteins) should equally be taken into account.

SOCIAL AND BIOLOGICAL DETERMINANTS OF SHORT-TERM GROWTH VELOCITY AMONG PRESCHOOLERS OF RURAL GUATEMALA: A PATH ANALYSIS APPROACH (WEIGHT-GAIN, IRON DEFICIENCY, GASTROINTESTINAL SYMPTOMS)

This dissertation develops and tests through path analysis a theoretical model to explain how socioeconomic, socioenvironmental, and biologic risk factors simultaneously influence each other to further produce short-term, depressed growth in preschoolers. Three areas of risk factors were identified: child's proximal environment, maturational stage, and biological vulnerability. The theoretical model represented both the conceptual framework and the nature and direction of the hypotheses. Original research completed in 1978-80 and in 1982 provided the background data. It was analyzed first by nested-analysis of variance, followed by path analysis. The study provided evidence of mild iron deficiency and gastrointestinal symptomatology in the etiology of depressed, short-term weight gain. Also, there was evidence suggesting that family resources for material and social survival significantly contribute to the variability of short-term, age-adjusted growth velocity. These results challenge current views of unifocal intervention, whether for prevention or control. For policy formulations, though, the mechanisms underlying any set of interlaced relationships must be decoded. Theoretical formulations here proposed should be reassessed under a more extensive research design. It is suggested that studies should be undertaken where social changes are actually in progress; otherwise, nutritional epidemiology in developing countries operates somewhere between social reality and research concepts, with little grasp of its real potential. The study stresses that there is a connection between substantive theory, empirical observation, and policy issues. ^

Functional analysis of the sea urchin {\it orthodenticle\/} -related gene, {\it SPOTX\/}, in aboral ectoderm-specific gene activation and aboral ectoderm cell fate specification

An important question in developmental biology is how embryonic cell types are derived from a fertilized egg. To address this question, this thesis investigates the mechanisms by which the aboral ectoderm-specific Spec2a gene is spatially and temporally regulated during sea urchin embryogenesis. The Spec2a gene of the sea urchin Strongylocentratus purpuratus has served as a valuable maker to understand the basis of lineage-specific gene activation and the role of transcription factors in cell fate specification. The hypothesis is that transcription factors responsible for cell type-specific gene activation are key components in the initial cell specification step. The Spec2a gene, which encodes a small cytosolic calcium-binding protein, is expressed exclusively in aboral ectoderm cell lineages. The 1516-bp control region of the Spec2a gene contains a 188-bp enhancer element required for temporal activation and aboral ectoderm/mesenchyme cell expression, while an unidentified element upstream of the enhancer represses expression in mesenchyme cells. Using an enhancer activation assay, combined with site-directed mutagenesis, I showed that three TAATCC/T sites within the enhancer are responsible for enhancer activity. Mutagenizing these sites and a fourth one just upstream abolished all activity from the Spec2a control region. A 77-bp DNA fragment from the Spec2a enhancer containing two of the TAATCC/T sites is sufficient for aboral ectoderm/mesenchyme cell expression. A cDNA encoding SpOtx, an orthodenticle-related protein, was cloned from S. purpuratus and shown to bind with high affinity to the TAATCC/T sequences within the Spec2a control region. SpOtx transcripts were found initially in all cells of the cleaving embryo, but they gradually became restricted to oral ectoderm and endoderm cells, suggesting that SpOtx might play a role in the initial temporal activation of the Spec2a gene and most likely has additional functions in the developing embryo. To reveal the broader biological functions of SpOtx, I injected SpOtx mRNA into living sea urchin eggs to determine what effects overexpressing the SpOtx protein might have on embryo development. SpOtx mRNA-injected embryos displayed dramatic alterations in development. Instead of developing into pluteus larvae with 15 different cell types, uniform epithelia balls were formed. These balls consisted of a thin layer of squamous cells with short cilia highly reminiscent of aboral ectoderm. Immunohistochemical staining and RT-PCR demonstrated that the SpOtx-injected embryoids expressed aboral ectoderm markers uniformly, but showed very weak or no expression of markers for non-aboral ectoderm cell types. These data strongly suggested that overexpression of SpOtx redirected the normal fate of non-aboral ectoderm cells to that of aboral ectoderm. These results show that SpOtx is involved in aboral ectoderm differentiation by activating aboral ectoderm-specific genes and that modulating its expression can lead to changes in cell fate. ^

Percentage of cultured cell-populations that stained positively and/or negatively for apoptotic markers cleaved caspase-3 and cleaved PARP, following DNA damage treatments induced by doxorubicin and TNF-alpha co-treatments

The present dataset contains the source data for Figure 2B of Tentner et al. (2012). The data shows the percentage of cultured cell-populations that stained positively and/or negatively for apoptotic markers cleaved caspase-3 and cleaved PARP, following DNA damage treatments induced by various doses of doxorubicin (0, 2 and 10 µmole/L) in the presence (100 ng/mL) or absence (0 ng/mL) of TNF-alpha co-treatment. For the six treatment conditions investigated, cell counts were made by flow cytometry at times 6, 12, 24, and 48 h following treatment; CULTURE DETAILS: U2OS cells were obtained from ATCC were maintained at 21% oxygen and 5% CO2 in Dulbecco's modified Eagle medium supplemented with 10% fetal bovine serum, penicillin, streptomycin, 2mM L-glutamine, and used within 15-20 passages. The first thymidine block was released by washing the plates three times with PBS, and incubating them in fresh thymidine-free media for 12 h. A second thymidine block was then performed by re-addition of thymidine to 2.5 mM followed by incubation for an additional 18 h. Media was aspirated, plates were washed 3 with PBS, and replaced with fresh media in the presence or absence of 10 mM aphidicolin; ANALYSIS DETAILS: See supplementary journal publication; RESULT: The authors of the supplementary journal publication conclude that TNF enhances dose-dependent cell death following doxorubicin-induced DNA damage with minimal affect on dose-dependent cell-cycle arrest.

Analysis of genetic variability in a sample of the durum wheat (Triticum durum Desf.) Spanish collection based on gliadin markers

In this work gliadin proteins were used to analyse the genetic variability in a sample of the durum wheat Spanish collection conserved at the CRF-INIA. In total 38 different alleles were identified at the loci Gli-A1, Gli-A3, Gli-B5, Gli-B1, Gli-A2 and Gli-B2. All the gliadin loci were polymorphic, possessed large genetic diversity and small and large differentiation within and between varieties, respectively. The Gli-A2 and Gli-B2 loci were the most polymorphic, the most fixed within varieties and the most useful to distinguish among varieties. Alternatively, Gli-B1 locus presented the least genetic variability out of the four main loci Gli-A1, Gli-B1, Gli-A2 and Gli-B2. The Gli-B1 alleles coding for the gliadin γ-45, associated with good quality, had an accumulated frequency of 69.7%, showing that the Spanish germplasm could be a good source for breeding quality. The Spanish landraces studied showed new gliadin alleles not catalogued so far. These new alleles might be associated with specific Spanish environment factors. The large number of new alleles identified also indicates that durum wheat Spanish germplasm is rather unique.

Comparative study of two laparoscopic instrument tracker designs for motion analysis and image-guided surgery: a technical evaluation

Laparoscopic instrument tracking systems are a key element in image-guided interventions, which requires high accuracy to be used in a real surgical scenario. In addition, these systems are a suitable option for objective assessment of laparoscopic technical skills based on instrument motion analysis. This study presents a new approach that improves the accuracy of a previously presented system, which applies an optical pose tracking system to laparoscopic practice. A design enhancement of the artificial markers placed on the laparoscopic instrument as well as an improvement of the calibration process are presented as a means to achieve more accurate results. A technical evaluation has been performed in order to compare the accuracy between the previous design and the new approach. Results show a remarkable improvement in the fluctuation error throughout the measurement platform. Moreover, the accumulated distance error and the inclination error have been improved. The tilt range covered by the system is the same for both approaches, from 90º to 7.5º. The relative position error is better for the new approach mainly at close distances to the camera system

Minisatellite marker analysis of Trypanosoma brucei: Reconciliation of clonal, panmictic, and epidemic population genetic structures

The African trypanosome, Trypanosoma brucei, has been shown to undergo genetic exchange in the laboratory, but controversy exists as to the role of genetic exchange in natural populations. Much of the analysis to date has been derived from isoenzyme or randomly amplified polymorphic DNA data with parasite material from a range of hosts and geographical locations. These markers fail to distinguish between the human infective (T. b. rhodesiense) and nonhuman infective (T. b. brucei) “subspecies” so that parasites derived from hosts other than humans potentially contain both subspecies. To overcome some of the inherent problems with the use of such markers and diverse populations, we have analyzed a well-defined population from a discrete geographical location (Busoga, Uganda) using three recently described minisatellite markers. The parasites were primarily isolated from humans and cattle with the latter isolates further characterized by their ability to resist lysis by human serum (equivalent to human infectivity). The minisatellite markers show high levels of polymorphism, and from the data obtained we conclude that T. b. rhodesiense is genetically isolated from T. b. brucei and can be unambiguously identified by its multilocus genotype. Analysis of the genotype frequencies in the separated T. b. brucei and T. b. rhodesiense populations shows the former has an epidemic population structure whereas the latter is clonal. This finding suggests that the strong linkage disequilibrium observed in previous analyses, where human and nonhuman infective trypanosomes were not distinguished, results from the treatment of two genetically isolated populations as a single population.

High-sensitivity array analysis of gene expression for the early detection of disseminated breast tumor cells in peripheral blood

Early detection is an effective means of reducing cancer mortality. Here, we describe a highly sensitive high-throughput screen that can identify panels of markers for the early detection of solid tumor cells disseminated in peripheral blood. The method is a two-step combination of differential display and high-sensitivity cDNA arrays. In a primary screen, differential display identified 170 candidate marker genes differentially expressed between breast tumor cells and normal breast epithelial cells. In a secondary screen, high-sensitivity arrays assessed expression levels of these genes in 48 blood samples, 22 from healthy volunteers and 26 from breast cancer patients. Cluster analysis identified a group of 12 genes that were elevated in the blood of cancer patients. Permutation analysis of individual genes defined five core genes (P ≤ 0.05, permax test). As a group, the 12 genes generally distinguished accurately between healthy volunteers and patients with breast cancer. Mean expression levels of the 12 genes were elevated in 77% (10 of 13) untreated invasive cancer patients, whereas cluster analysis correctly classified volunteers and patients (P = 0.0022, Fisher's exact test). Quantitative real-time PCR confirmed array results and indicated that the sensitivity of the assay (1:2 × 108 transcripts) was sufficient to detect disseminated solid tumor cells in blood. Expression-based blood assays developed with the screening approach described here have the potential to detect and classify solid tumor cells originating from virtually any primary site in the body.

Half tetrad analysis in alfalfa using multiple restriction fragment length polymorphism markers.

A maximum likelihood approach of half tetrad analysis (HTA) based on multiple restriction fragment length polymorphism (RFLP) markers was developed. This procedure estimates the relative frequencies of 2n gametes produced by mechanisms genetically equivalent to first division restitution (FDR) or second division restitution and simultaneously locates the centromere within a linkage group of RFLP marker loci. The method was applied to the diploid alfalfa clone PG-F9 (2n = 2x = 16) previously selected because of its high frequency of 2n egg production. HTA was based on four RFLP loci for which PG-F9 was heterozygous with codominant alleles that were absent in the tetraploid tester. Models including three linked and one unlinked RFLP loci were developed and tested. Results of the HTA showed that PG-F9 produced 6% FDR and 94% second division restitution 2n eggs. Information from a marker locus belonging to one linkage group was used to more precisely locate the centromere on a different linkage group. HTA, together with previous cytological analysis, indicated that in PG-F9, FDR 2n eggs are likely produced by diplospory, a mechanism common among apomictic species. The occurrence of FDR 2n eggs in plant species and their importance for crop evolution and breeding is discussed together with the potential applicability of multilocus HTA in the study of reproductive mutants.