![155980-Thumbnail Image.png](https://d1rbsgppyrdqq4.cloudfront.net/s3fs-public/styles/width_400/public/2021-08/155980-Thumbnail%20Image.png?versionId=P3TbjB2qKDdWomyodp6JuAmi9boL4ykh&X-Amz-Content-Sha256=UNSIGNED-PAYLOAD&X-Amz-Algorithm=AWS4-HMAC-SHA256&X-Amz-Credential=AKIASBVQ3ZQ42ZLA5CUJ/20240615/us-west-2/s3/aws4_request&X-Amz-Date=20240615T002512Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=19cbb97cf10296c9dc24cf952b21b406d5aeea3edfc1ab9e795009e817a5468b&itok=ql6JWeXj)
In this dissertation, I evaluate the ecological drivers and fitness consequences of non-kin queen cooperation, by comparing the reproduction of mature single-queen versus polygynous harvester ant (Pogonomyrmex californicus) colonies in the field. I captured and quantified the total number and biomass of reproductives across multiple mating seasons, comparing between populations that vary in the proportion of single queen versus polygynous colonies, to assess the fitness outcomes of queen cooperation. Colonies in a mainly polygynous site had lower reproductive investment than those in sites with predominantly single-queen colonies. The site dominated by polygyny had higher colony density and displayed evidence of resource limitation, pressures that may drive the evolution of queen cooperation.
I also used microsatellite markers to examine how polygynous queens share worker and reproductive production with nest-mate queens. The majority of queens fairly contribute to worker production and equally share reproductive output. However, there is a low frequency of queens that under-produce workers and over-produce reproductive offspring. This suggests that cheating by reproducing queens is possible, but uncommon. Competitive pressure from neighboring colonies could reduce the success of colonies that contain cheaters and maintain a low frequency of this phenotype in the population.
![155994-Thumbnail Image.png](https://d1rbsgppyrdqq4.cloudfront.net/s3fs-public/styles/width_400/public/2021-08/155994-Thumbnail%20Image.png?versionId=sTLdy1U3u6bsNuT8vAaAtKnJSfBtDAXJ&X-Amz-Content-Sha256=UNSIGNED-PAYLOAD&X-Amz-Algorithm=AWS4-HMAC-SHA256&X-Amz-Credential=AKIASBVQ3ZQ42ZLA5CUJ/20240616/us-west-2/s3/aws4_request&X-Amz-Date=20240616T023934Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=f106e78fc94e744e362ab9ab2b7e59380c1eec5b5572007fd0153ea31afdbcbc&itok=3SbKzrRr)
This dissertation proposes two PageRank-based analytical methods, Pathways of Topological Rank Analysis (PoTRA) and miR2Pathway, discussed in Chapter 1 and Chapter 2, respectively. PoTRA focuses on detecting pathways with an altered number of hub genes in corresponding pathways between two phenotypes. The basis for PoTRA is that the loss of connectivity is a common topological trait of cancer networks, as well as the prior knowledge that a normal biological network is a scale-free network whose degree distribution follows a power law where a small number of nodes are hubs and a large number of nodes are non-hubs. However, from normal to cancer, the process of the network losing connectivity might be the process of disrupting the scale-free structure of the network, namely, the number of hub genes might be altered in cancer compared to that in normal samples. Hence, it is hypothesized that if the number of hub genes is different in a pathway between normal and cancer, this pathway might be involved in cancer. MiR2Pathway focuses on quantifying the differential effects of miRNAs on the activity of a biological pathway when miRNA-mRNA connections are altered from normal to disease and rank disease risk of rewired miRNA-mediated biological pathways. This dissertation explores how rewired gene-gene interactions and rewired miRNA-mRNA interactions lead to aberrant activity of biological pathways, and rank pathways for their disease risk. The two methods proposed here can be used to complement existing genomics analysis methods to facilitate the study of biological mechanisms behind disease at the systems-level.
![157966-Thumbnail Image.png](https://d1rbsgppyrdqq4.cloudfront.net/s3fs-public/styles/width_400/public/2021-09/157966-Thumbnail%20Image.png?versionId=YQ70Qb938z6qwGzju5mwwe9lBMkJkS5I&X-Amz-Content-Sha256=UNSIGNED-PAYLOAD&X-Amz-Algorithm=AWS4-HMAC-SHA256&X-Amz-Credential=AKIASBVQ3ZQ42ZLA5CUJ/20240617/us-west-2/s3/aws4_request&X-Amz-Date=20240617T050226Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=897fbcee377bc0c26b6139bb9ee1c0990f4b377107760f9e4dbf896986174035&itok=_eZXnAFq)
designing personalized treatments and improving clinical outcomes of cancers. Such
investigations require accurate delineation of the subclonal composition of a tumor, which
to date can only be reliably inferred from deep-sequencing data (>300x depth). The
resulting algorithm from the work presented here, incorporates an adaptive error model
into statistical decomposition of mixed populations, which corrects the mean-variance
dependency of sequencing data at the subclonal level and enables accurate subclonal
discovery in tumors sequenced at standard depths (30-50x). Tested on extensive computer
simulations and real-world data, this new method, named model-based adaptive grouping
of subclones (MAGOS), consistently outperforms existing methods on minimum
sequencing depth, decomposition accuracy and computation efficiency. MAGOS supports
subclone analysis using single nucleotide variants and copy number variants from one or
more samples of an individual tumor. GUST algorithm, on the other hand is a novel method
in detecting the cancer type specific driver genes. Combination of MAGOS and GUST
results can provide insights into cancer progression. Applications of MAGOS and GUST
to whole-exome sequencing data of 33 different cancer types’ samples discovered a
significant association between subclonal diversity and their drivers and patient overall
survival.
![158771-Thumbnail Image.png](https://d1rbsgppyrdqq4.cloudfront.net/s3fs-public/styles/width_400/public/2021-09/158771-Thumbnail%20Image.png?versionId=64.R10SUWEL_AxkWDD4lICRxoR_gU2eo&X-Amz-Content-Sha256=UNSIGNED-PAYLOAD&X-Amz-Algorithm=AWS4-HMAC-SHA256&X-Amz-Credential=AKIASBVQ3ZQ42ZLA5CUJ/20240617/us-west-2/s3/aws4_request&X-Amz-Date=20240617T183038Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=0cd6f09631e6c77d366906657b7d96ef85145be3c91b3ebb7df4d44640abb240&itok=pzfZs3rU)
![161916-Thumbnail Image.png](https://d1rbsgppyrdqq4.cloudfront.net/s3fs-public/styles/width_400/public/2021-11/161916-Thumbnail%20Image.png?versionId=nHWAScIfIC3rTQR0bK77xR9glYO7cZ6T&X-Amz-Content-Sha256=UNSIGNED-PAYLOAD&X-Amz-Algorithm=AWS4-HMAC-SHA256&X-Amz-Credential=AKIASBVQ3ZQ42ZLA5CUJ/20240617/us-west-2/s3/aws4_request&X-Amz-Date=20240617T042505Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=8b98e17d4dbfe0b39f1d8aa1953a6fa18758ad9c5da77627b8ba41a3b7b77cae&itok=_1j9bU6J)
![129370-Thumbnail Image.png](https://d1rbsgppyrdqq4.cloudfront.net/s3fs-public/styles/width_400/public/2021-04/129370-Thumbnail%20Image.png?versionId=YV4F288JP9C07bqu7WYS2oPetI4EKvUm&X-Amz-Content-Sha256=UNSIGNED-PAYLOAD&X-Amz-Algorithm=AWS4-HMAC-SHA256&X-Amz-Credential=AKIASBVQ3ZQ42ZLA5CUJ/20240615/us-west-2/s3/aws4_request&X-Amz-Date=20240615T081655Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=d6cadd3f6e6c6ec2a323b8669a803db9bd46311e8db2a226d2496722d9ff2601&itok=1H-XKn8y)
Adaptation requires genetic variation, but founder populations are generally genetically depleted. Here we sequence two populations of an inbred ant that diverge in phenotype to determine how variability is generated. Cardiocondyla obscurior has the smallest of the sequenced ant genomes and its structure suggests a fundamental role of transposable elements (TEs) in adaptive evolution. Accumulations of TEs (TE islands) comprising 7.18% of the genome evolve faster than other regions with regard to single-nucleotide variants, gene/exon duplications and deletions and gene homology. A non-random distribution of gene families, larvae/adult specific gene expression and signs of differential methylation in TE islands indicate intragenomic differences in regulation, evolutionary rates and coalescent effective population size. Our study reveals a tripartite interplay between TEs, life history and adaptation in an invasive species.
![128995-Thumbnail Image.png](https://d1rbsgppyrdqq4.cloudfront.net/s3fs-public/styles/width_400/public/2021-04/128995-Thumbnail%20Image.png?versionId=_y7EfM094PXEGa1kmvpHtpWfZBnjQK97&X-Amz-Content-Sha256=UNSIGNED-PAYLOAD&X-Amz-Algorithm=AWS4-HMAC-SHA256&X-Amz-Credential=AKIASBVQ3ZQ42ZLA5CUJ/20240615/us-west-2/s3/aws4_request&X-Amz-Date=20240615T054722Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=31c547944c16445bdcfb9d73970c4a1edcf2f7c2d43431e18df0602a620d9a0e&itok=KddQAOWN)
Background: Obesity is a metabolic disease caused by environmental and genetic factors. However, the epigenetic mechanisms of obesity are incompletely understood. The aim of our study was to investigate the role of skeletal muscle DNA methylation in combination with transcriptomic changes in obesity.
Results: Muscle biopsies were obtained basally from lean (n = 12; BMI = 23.4 ± 0.7 kg/m[superscript 2]) and obese (n = 10; BMI = 32.9 ± 0.7 kg/m[superscript 2]) participants in combination with euglycemic-hyperinsulinemic clamps to assess insulin sensitivity. We performed reduced representation bisulfite sequencing (RRBS) next-generation methylation and microarray analyses on DNA and RNA isolated from vastus lateralis muscle biopsies. There were 13,130 differentially methylated cytosines (DMC; uncorrected P < 0.05) that were altered in the promoter and untranslated (5' and 3'UTR) regions in the obese versus lean analysis. Microarray analysis revealed 99 probes that were significantly (corrected P < 0.05) altered. Of these, 12 genes (encompassing 22 methylation sites) demonstrated a negative relationship between gene expression and DNA methylation. Specifically, sorbin and SH3 domain containing 3 (SORBS3) which codes for the adapter protein vinexin was significantly decreased in gene expression (fold change −1.9) and had nine DMCs that were significantly increased in methylation in obesity (methylation differences ranged from 5.0 to 24.4 %). Moreover, differentially methylated region (DMR) analysis identified a region in the 5'UTR (Chr.8:22,423,530–22,423,569) of SORBS3 that was increased in methylation by 11.2 % in the obese group. The negative relationship observed between DNA methylation and gene expression for SORBS3 was validated by a site-specific sequencing approach, pyrosequencing, and qRT-PCR. Additionally, we performed transcription factor binding analysis and identified a number of transcription factors whose binding to the differentially methylated sites or region may contribute to obesity.
Conclusions: These results demonstrate that obesity alters the epigenome through DNA methylation and highlights novel transcriptomic changes in SORBS3 in skeletal muscle.
![128932-Thumbnail Image.png](https://d1rbsgppyrdqq4.cloudfront.net/s3fs-public/styles/width_400/public/2021-04/128932-Thumbnail%20Image.png?versionId=exwtcOP1SXrEX0aA2Fb0qbd9fhzB.J9M&X-Amz-Content-Sha256=UNSIGNED-PAYLOAD&X-Amz-Algorithm=AWS4-HMAC-SHA256&X-Amz-Credential=AKIASBVQ3ZQ42ZLA5CUJ/20240617/us-west-2/s3/aws4_request&X-Amz-Date=20240617T193416Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=836ac100aadef9836c30ff5d50c201920d9ad12d3f8e5c5154f6c51751247cad&itok=2OFNADC5)
We have previously hypothesized a biological pathway of activity-dependent synaptic plasticity proteins that addresses the dual genetic and environmental contributions to schizophrenia. Accordingly, variations in the immediate early gene EGR3, and its target ARC, should influence schizophrenia susceptibility. We used a pooled Next-Generation Sequencing approach to identify variants across these genes in U.S. populations of European (EU) and African (AA) descent. Three EGR3 and one ARC SNP were selected and genotyped for validation, and three SNPs were tested for association in a replication cohort. In the EU group of 386 schizophrenia cases and 150 controls EGR3 SNP rs1877670 and ARC SNP rs35900184 showed significant associations (p = 0.0078 and p = 0.0275, respectively). In the AA group of 185 cases and 50 controls, only the ARC SNP revealed significant association (p = 0.0448). The ARC SNP did not show association in the Han Chinese (CH) population. However, combining the EU, AA, and CH groups revealed a highly significant association of ARC SNP rs35900184 (p = 2.353 x 10-7; OR [95% CI] = 1.54 [1.310–1.820]). These findings support previously reported associations between EGR3 and schizophrenia. Moreover, this is the first report associating an ARC SNP with schizophrenia and supports recent large-scale GWAS findings implicating the ARC complex in schizophrenia risk. These results support the need for further investigation of the proposed pathway of environmentally responsive, synaptic plasticity-related, schizophrenia genes.
![128958-Thumbnail Image.png](https://d1rbsgppyrdqq4.cloudfront.net/s3fs-public/styles/width_400/public/2021-04/128958-Thumbnail%20Image.png?versionId=tvshD1VebMDZkXhDMmh4L7CJo3ypcNv5&X-Amz-Content-Sha256=UNSIGNED-PAYLOAD&X-Amz-Algorithm=AWS4-HMAC-SHA256&X-Amz-Credential=AKIASBVQ3ZQ42ZLA5CUJ/20240615/us-west-2/s3/aws4_request&X-Amz-Date=20240615T220909Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=c5deb9610a0009fde48dbaf759a579b196e5260ad6a2c127a29b5bcbb9fefd57&itok=qLLcjpvx)
Background: Immunosignaturing is a new peptide microarray based technology for profiling of humoral immune responses. Despite new challenges, immunosignaturing gives us the opportunity to explore new and fundamentally different research questions. In addition to classifying samples based on disease status, the complex patterns and latent factors underlying immunosignatures, which we attempt to model, may have a diverse range of applications.
Methods: We investigate the utility of a number of statistical methods to determine model performance and address challenges inherent in analyzing immunosignatures. Some of these methods include exploratory and confirmatory factor analyses, classical significance testing, structural equation and mixture modeling.
Results: We demonstrate an ability to classify samples based on disease status and show that immunosignaturing is a very promising technology for screening and presymptomatic screening of disease. In addition, we are able to model complex patterns and latent factors underlying immunosignatures. These latent factors may serve as biomarkers for disease and may play a key role in a bioinformatic method for antibody discovery.
Conclusion: Based on this research, we lay out an analytic framework illustrating how immunosignatures may be useful as a general method for screening and presymptomatic screening of disease as well as antibody discovery.
![129066-Thumbnail Image.png](https://d1rbsgppyrdqq4.cloudfront.net/s3fs-public/styles/width_400/public/2021-04/129066-Thumbnail%20Image.png?versionId=cQihiHzCwAvA7C2wHXkQT7beDCah_B_t&X-Amz-Content-Sha256=UNSIGNED-PAYLOAD&X-Amz-Algorithm=AWS4-HMAC-SHA256&X-Amz-Credential=AKIASBVQ3ZQ42ZLA5CUJ/20240617/us-west-2/s3/aws4_request&X-Amz-Date=20240617T064652Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=784d7fa31b166f719b2c23f06b83a91b9f5073d0135521037a01307eb457276d&itok=XA7zXwDv)
Background: Glioblastoma is the most aggressive primary central nervous tumor and carries a very poor prognosis. Invasion precludes effective treatment and virtually assures tumor recurrence. In the current study, we applied analytical and bioinformatics approaches to identify a set of microRNAs (miRs) from several different human glioblastoma cell lines that exhibit significant differential expression between migratory (edge) and migration-restricted (core) cell populations. The hypothesis of the study is that differential expression of miRs provides an epigenetic mechanism to drive cell migration and invasion.
Results: Our research data comprise gene expression values for a set of 805 human miRs collected from matched pairs of migratory and migration-restricted cell populations from seven different glioblastoma cell lines. We identified 62 down-regulated and 2 up-regulated miRs that exhibit significant differential expression in the migratory (edge) cell population compared to matched migration-restricted (core) cells. We then conducted target prediction and pathway enrichment analysis with these miRs to investigate potential associated gene and pathway targets. Several miRs in the list appear to directly target apoptosis related genes. The analysis identifies a set of genes that are predicted by 3 different algorithms, further emphasizing the potential validity of these miRs to promote glioblastoma.
Conclusions: The results of this study identify a set of miRs with potential for decreased expression in invasive glioblastoma cells. The verification of these miRs and their associated targeted proteins provides new insights for further investigation into therapeutic interventions. The methodological approaches employed here could be applied to the study of other diseases to provide biomedical researchers and clinicians with increased opportunities for therapeutic interventions.