In silico experimental modeling of cancer treatment

In silico experimental modeling of cancer involves combining findings from biological literature with computer-based models of biological systems in order to conduct investigations of hypotheses entirely in the computer laboratory. In this paper, we discuss the use of in silico modeling as a precursor to traditional clinical and laboratory research, allowing researchers to refine their experimental programs with an aim to reducing costs and increasing research efficiency. We explain the methodology of in silico experimental trials before providing an example of in silico modeling from the biomathematical literature with a view to promoting more widespread use and understanding of this research strategy.

Candidate gene association studies : a comprehensive guide to useful in silico tools

The candidate gene approach has been a pioneer in the field of genetic epidemiology, identifying risk alleles and their association with clinical traits. With the advent of rapidly changing technology, there has been an explosion of in silico tools available to researchers, giving them fast, efficient resources and reliable strategies important to find casual gene variants for candidate or genome wide association studies (GWAS). In this review, following a description of candidate gene prioritisation, we summarise the approaches to single nucleotide polymorphism (SNP) prioritisation and discuss the tools available to assess functional relevance of the risk variant with consideration to its genomic location. The strategy and the tools discussed are applicable to any study investigating genetic risk factors associated with a particular disease. Some of the tools are also applicable for the functional validation of variants relevant to the era of GWAS and next generation sequencing (NGS).

In silico analyses reveal common cellular pathways affected by loss of heterozygosity (LOH) events in the lymphomagenesis of Non-Hodgkin’s lymphoma (NHL)

Background The analysis of cellular networks and pathways involved in oncogenesis has increased our knowledge about the pathogenic mechanisms that underlie tumour biology and has unmasked new molecular targets that may lead to the design of better anti-cancer therapies. Recently, using a high resolution loss of heterozygosity (LOH) analysis, we identified a number of potential tumour suppressor genes (TSGs) within common LOH regions across cases suffering from two of the most common forms of Non-Hodgkin’s lymphoma (NHL), Follicular Lymphoma (FL) and Diffuse Large B-cell Lymphoma (DLBCL). From these studies LOH of the protein tyrosine phosphatase receptor type J (PTPRJ) gene was identified as a common event in the lymphomagenesis of these B-cell lymphomas. The present study aimed to determine the cellular pathways affected by the inactivation of these TSGs including PTPRJ in FL and DLBCL tumourigenesis. Results Pathway analytical approaches identified that candidate TSGs located within common LOH regions participate within cellular pathways, which may play a crucial role in FL and DLBCL lymphomagenesis (i.e., metabolic pathways). These analyses also identified genes within the interactome of PTPRJ (i.e. PTPN11 and B2M) that when inactivated in NHL may play an important role in tumourigenesis. We also detected genes that are differentially expressed in cases with and without LOH of PTPRJ, such as NFATC3 (nuclear factor of activated T-cells, cytoplasmic, calcineurin-dependent 3). Moreover, upregulation of the VEGF, MAPK and ERBB signalling pathways was also observed in NHL cases with LOH of PTPRJ, indicating that LOH-driving events causing inactivation of PTPRJ, apart from possibly inducing a constitutive activation of these pathways by reduction or abrogation of its dephosphorylation activity, may also induce upregulation of these pathways when inactivated. This finding implicates these pathways in the lymphomagenesis and progression of FL and DLBCL. Conclusions The evidence obtained in this research supports findings suggesting that FL and DLBCL share common pathogenic mechanisms. Also, it indicates that PTPRJ can play a crucial role in the pathogenesis of these B-cell tumours and suggests that activation of PTPRJ might be an interesting novel chemotherapeutic target for the treatment of these B-cell tumours.

In silico biology : making the most of parallel computing

Biological systems are typically complex and adaptive, involving large numbers of entities, or organisms, and many-layered interactions between these. System behaviour evolves over time, and typically benefits from previous experience by retaining memory of previous events. Given the dynamic nature of these phenomena, it is non-trivial to provide a comprehensive description of complex adaptive systems and, in particular, to define the importance and contribution of low-level unsupervised interactions to the overall evolution process. In this chapter, the authors focus on the application of the agent-based paradigm in the context of the immune response to HIV. Explicit implementation of lymph nodes and the associated lymph network, including lymphatic chain structure, is a key objective, and requires parallelisation of the model. Steps taken towards an optimal communication strategy are detailed.

In vitro and in vivo studies on protease-activated receptor-2

Protease-activated receptor-2 (PAR2) is a G protein coupled receptor (GPCR) that is activated by proteolytic cleavage of its amino terminal domain by trypsin-like serine proteases. Cleavage of this receptor exposes a neoepitope, termed the tethered ligand (TL), which binds intramolecularly within the receptor to stimulate signal transduction via coupled G proteins. PAR2-mediated signal transduction is also experimentally stimulated by hexapeptides (agonist peptides; APs) that are homologous to the TL sequence. Due to the irreversible nature of PAR2 proteolysis, downstream signal transduction is tightly regulated. Following activation, PAR2 is rapidly uncoupled from downstream signalling by the post-translational modifications phosphorylation and ubiquination which facilitate interactions with â- arrestin. This scaffolding protein couples PAR2 to the internalisation machinery initiating its desensitisation and trafficking through the early and late endosomes followed by receptor degradation. PAR2 is widely expressed in mammalian tissues with key roles for this receptor in cardiovascular, respiratory, nervous and musculoskeletal systems. This receptor has also been linked to pathological states with aberrant expression and signalling noted in several cancers. In prostate cancer, PAR2 signalling induces migration and proliferation of tumour derived cell lines, while elevated receptor expression has been noted in malignant tissues. Importantly, a role for this receptor has also been suggested in prostate cancer bone metastasis as coexpression of PAR2 and a proteolytic activator has been demonstrated by immunohistochemical analysis. Based on these data, the primary focus of this project has been on two aspects of PAR2 biology. The first is characterisation of cellular mechanisms that regulate PAR2 signalling and trafficking. The second aspect is the role of this receptor in prostate cancer bone metastasis. In addition, to permit these studies, it was first necessary to evaluate the specificity of the commercially available anti-PAR2 antibodies SAM11, C17, N19 and H99. The evaluation of the four commercially available antibodies was assessed using four techniques: immunoprecipitation; Western blot analysis; immunofluorescence; and flow cytometry. These approaches demonstrated that three of the antibodies efficiently detect ectopically expressed PAR2 by each of these techniques. A significant finding from this study was that N19 was the only antibody able to specifically detect N-glycosylated endogenous PAR2 by Western blot analysis. This analysis was performed on lysates from prostate cancer derived cell lines and tissue derived from wildtype and PAR2 knockout mice. Importantly, further evaluation demonstrated that this antibody also efficiently detects endogenous PAR2 at the cell surface by flow cytometry. The anti-PAR2 antibody N19 was used to explore the in vitro role of palmitoylation, the post-translational addition of palmitate, in PAR2 signalling, trafficking, cell surface expression and desensitization. Significantly, use of the palmitoylation inhibitor 2-bromopalmitate indicated that palmitate addition is important in trafficking of PAR2 endogenously expressed by prostate cancer cell lines. This was supported by palmitate labelling experiments using two approaches which showed that PAR2 stably expressed by CHO cells is palmitoylated and that palmitoylation occurs on cysteine 361. Another key finding from this study is that palmitoylation is required for optimal PAR2 signalling as Ca2+ flux assays indicated that in response to trypsin agonism, palmitoylation deficient PAR2 is ~9 fold less potent than wildtype receptor with a reduction of about 33% in the maximum signal induced via the mutant receptor. Confocal microscopy, flow cytometry and cell surface biotinylation analyses demonstrated that palmitoylation is required for efficient cell surface expression of PAR2. Importantly, this study also identified that palmitoylation of this receptor within the Golgi apparatus is required for efficient agonist-induced rab11amediated trafficking of PAR2 to the cell surface. Interestingly, palmitoylation is also required for receptor desensitization, as agonist-induced â-arrestin recruitment and receptor degradation were markedly reduced in CHO-PAR2-C361A cells compared with CHO-PAR2 cells. Collectively, these data provide new insights on the life cycle of PAR2 and demonstrate that palmitoylation is critical for efficient signalling, trafficking, cell surface localization and degradation of this receptor. This project also evaluated PAR2 residues involved in ligand docking. Although the extracellular loop (ECL)2 of PAR2 is known to be required for agonist-induced signal transduction, the binding pocket for receptor agonists remains to be determined. In silico homology modelling, based on a crystal structure for the prototypical GPCR rhodopsin, and ligand docking were performed to identify PAR2 transmembrane (TM) amino acids potentially involved in agonist binding. These methods identified 12 candidate residues that were mutated to examine the binding site of the PAR2 TL, revealed by trypsin cleavage, as well as of the soluble ligands 2f-LIGRLO-NH2 and GB110, which are both structurally based on the AP SLIGRLNH2. Ligand binding was evaluated from the impact of the mutated residues on PAR2-mediated calcium mobilisation. An important finding from these experiments was that mutation of residues Y156 and Y326 significantly reduced 2f-LIGRLO-NH2 and GB110 agonist activity. L307 was also important for GB110 activity. Intriguingly, mutation of PAR2 residues did not alter trypsin-induced signalling to the same extent as for the soluble agonists. The reason for this difference remains to be further examined by in silico and in vitro experimentation and, potentially, crystal structure studies. However, these findings identified the importance of TM domains in PAR2 ligand docking and will enhance the design of both PAR2 agonists and potentially agents to inhibit signalling (antagonists). The potential importance of PAR2 in prostate cancer bone metastasis was examined using a mouse model. In patients, prostate cancer bone metastases cause bone growth by disrupting bone homeostasis. In an attempt to mimic prostate cancer growth in bone, PAR2 responsive 22Rv1 prostate cancer cells, which form mixed osteoblastic and osteolytic lesions, were injected into the proximal aspect of mouse tibiae. A role for PAR2 was assessed by treating these mice with the recently developed PAR2 antagonist GB88. As controls, animals bearing intra-tibial tumours were also treated with vehicle (olive oil) or the prostate cancer chemotherapeutic docetaxel. The effect of these treatments on bone was examined radiographically and by micro-CT. Consistent with previous studies, 22Rv1 tumours caused osteoblastic periosteal spicule formation and concurrent osteolytic bone loss. Significantly, blockade of PAR2 signalling reduced the osteoblastic and osteolytic phenotype of 22Rv1 tumours in bone. No bone defects were detected in mice treated with docetaxel. These qualitative data will be followed in the future by quantitative micro-CT analysis as well as histology and histomorphometry analysis of already collected tissues. Nonetheless, these preliminary experiments highlight a potential role for PAR2 in prostate cancer growth in bone. In summary, in vitro studies have defined mechanisms regulating PAR2 activation, downstream signalling and trafficking and in vivo studies point to a potential role for this receptor in prostate cancer bone metastasis. The outcomes of this project are that a greater understanding of the biology of PAR2 may lead to the development of strategies to modulate the function of this receptor in disease.

Genetic polymorphisms in the human tissue kallikrein (KLK) locus and their implication in various malignant and non-malignant diseases

The Kallikrein (KLK) gene locus encodes a family of serine proteases and is the largest contiguous cluster of protease-encoding genes attributed an evolutionary age of 330 million years. The KLK locus has been implicated as a high susceptibility risk loci in numerous cancer studies through the last decade. The KLK3 gene already has established clinical relevance as a biomarker in prostate cancer prognosis through its encoded protein, prostate-specific antigen. Data mined through genome-wide association studies (GWAS) and next-generation sequencing point to many important candidate single nucleotide polymorphisms (SNPs) in KLK3 and other KLK genes. SNPs in the KLK locus have been found to be associated with several diseases including cancer, hypertension, cardiovascular disease and atopic dermatitis. Moreover, introducing a model incorporating SNPs to improve the efficiency of prostate-specific antigen in detecting malignant states of prostate cancer has been recently suggested. Establishing the functional relevance of these newly-discovered SNPs, and their interactions with each other, through in silico investigations followed by experimental validation, can accelerate the discovery of diagnostic and prognostic biomarkers. In this review, we discuss the various genetic association studies on the KLK loci identified either through candidate gene association studies or at the GWAS and post-GWAS front to aid researchers in streamlining their search for the most significant, relevant and therapeutically promising candidate KLK gene and/or SNP for future investigations.

The ancient gene C12orf29 : an exploration of its role in the chordate body plan

The sheep (Ovis aries) is commonly used as a large animal model in skeletal research. Although the sheep genome has been sequenced there are still only a limited number of annotated mRNA sequences in public databases. A complementary DNA (cDNA) library was constructed to provide a generic resource for further exploration of genes that are actively expressed in bone cells in sheep. It was anticipated that the cDNA library would provide molecular tools for further research into the process of fracture repair and bone homeostasis, and add to the existing body of knowledge. One of the hallmarks of cDNA libraries has been the identification of novel genes and in this library the full open reading frame of the gene C12orf29 was cloned and characterised. This gene codes for a protein of unknown function with a molecular weight of 37 kDa. A literature search showed that no previous studies had been conducted into the biological role of C12orf29, except for some bioinformatics studies that suggested a possible link with cancer. Phylogenetic analyses revealed that C12orf29 had an ancient pedigree with a homologous gene found in some bacterial taxa. This implied that the gene was present in the last common eukaryotic ancestor, thought to have existed more than 2 billion years ago. This notion was further supported by the fact that the gene is found in taxa belonging to the two major eukaryotic branches, bikonts and unikonts. In the bikont supergroup a C12orf29-like gene was found in the single celled protist Naegleria gruberi, whereas in the unikont supergroup, encompassing the metazoa, the gene is universal to all chordate and, therefore, vertebrate species. It appears to have been lost to the majority of cnidaria and protostomes taxa; however, C12orf29-like genes have been found in the cnidarian freshwater hydra and the protostome Pacific oyster. The experimental data indicate that C12orf29 has a structural role in skeletal development and tissue homeostasis, whereas in silico analysis of the human C12orf29 promoter region suggests that its expression is potentially under the control of the NOTCH, WNT and TGF- developmental pathways, as well SOX9 and BAPX1; pathways that are all heavily involved in skeletogenesis. Taken together, this investigation provides strong evidence that C12orf29 has a very important role in the chordate body plan, in early skeletal development, cartilage homeostasis, and also a possible link with spina bifida in humans.

Targeting amino acid transport in metastatic castration-resistant prostate cancer : effects on cell cycle, cell growth, and tumor development

Background L-type amino acid transporters (LATs) uptake neutral amino acids including L-leucine into cells, stimulating mammalian target of rapamycin complex 1 signaling and protein synthesis. LAT1 and LAT3 are overexpressed at different stages of prostate cancer, and they are responsible for increasing nutrients and stimulating cell growth. Methods We examined LAT3 protein expression in human prostate cancer tissue microarrays. LAT function was inhibited using a leucine analog (BCH) in androgen-dependent and -independent environments, with gene expression analyzed by microarray. A PC-3 xenograft mouse model was used to study the effects of inhibiting LAT1 and LAT3 expression. Results were analyzed with the Mann-Whitney U or Fisher exact tests. All statistical tests were two-sided. Results LAT3 protein was expressed at all stages of prostate cancer, with a statistically significant decrease in expression after 4–7 months of neoadjuvant hormone therapy (4–7 month mean = 1.571; 95% confidence interval = 1.155 to 1.987 vs 0 month = 2.098; 95% confidence interval = 1.962 to 2.235; P = .0187). Inhibition of LAT function led to activating transcription factor 4–mediated upregulation of amino acid transporters including ASCT1, ASCT2, and 4F2hc, all of which were also regulated via the androgen receptor. LAT inhibition suppressed M-phase cell cycle genes regulated by E2F family transcription factors including critical castration-resistant prostate cancer regulatory genes UBE2C, CDC20, and CDK1. In silico analysis of BCH-downregulated genes showed that 90.9% are statistically significantly upregulated in metastatic castration-resistant prostate cancer. Finally, LAT1 or LAT3 knockdown in xenografts inhibited tumor growth, cell cycle progression, and spontaneous metastasis in vivo. Conclusion Inhibition of LAT transporters may provide a novel therapeutic target in metastatic castration-resistant prostate cancer, via suppression of mammalian target of rapamycin complex 1 activity and M-phase cell cycle genes.

How much information can be obtained from tracking the position of the leading edge in a scratch assay?

Moving cell fronts are an essential feature of wound healing, development and disease. The rate at which a cell front moves is driven, in part, by the cell motility, quantified in terms of the cell diffusivity $D$, and the cell proliferation rate �$\lambda$. Scratch assays are a commonly-reported procedure used to investigate the motion of cell fronts where an initial cell monolayer is scratched and the motion of the front is monitored over a short period of time, often less than 24 hours. The simplest way of quantifying a scratch assay is to monitor the progression of the leading edge. Leading edge data is very convenient since, unlike other methods, it is nondestructive and does not require labeling, tracking or counting individual cells amongst the population. In this work we study short time leading edge data in a scratch assay using a discrete mathematical model and automated image analysis with the aim of investigating whether such data allows us to reliably identify $D$ and $\lambda$�. Using a naıve calibration approach where we simply scan the relevant region of the ($D$;$\lambda$�) parameter space, we show that there are many choices of $D$ and $\lambda$� for which our model produces indistinguishable short time leading edge data. Therefore, without due care, it is impossible to estimate $D$ and $\lambda$� from this kind of data. To address this, we present a modified approach accounting for the fact that cell motility occurs over a much shorter time scale than proliferation. Using this information we divide the duration of the experiment into two periods, and we estimate $D$ using data from the first period, while we estimate �$\lambda$ using data from the second period. We confirm the accuracy of our approach using in silico data and a new set of in vitro data, which shows that our method recovers estimates of $D$ and $\lamdba$� that are consistent with previously-reported values except that that our approach is fast, inexpensive, nondestructive and avoids the need for cell labeling and cell counting.

Carbon nanotube-based super nanotubes: Tunable thermal conductivity in three dimensions

Advances in nanomaterials/nanostructures offer the possibility of fabricating multifunctional materials for use in engineering applications. Carbon nanotube (CNT)-based nanostructures are a representative building block for these multifunctional materials. Based on a series of in silico studies, we investigated the possibility of tuning the thermal conductivity of a three-dimensional CNT-based nanostructure: a single-walled CNT-based super-nanotube. The thermal conductivity of the super-nanotubes was shown to vary with different connecting carbon rings and super-nanotubes with longer constituent single-walled CNTs and larger diameters had a smaller thermal conductivity. The inverse of the thermal conductivity of the super-nanotubes showed a good linear relationship with the inverse of the length. The thermal conductivity was approximately proportional to the inverse of the temperature, but was insensitive to the axial strain as a result of the Poisson ratio. These results provide a fundamental understanding of the thermal conductivity of the super-nanotubes and will guide their future design/fabrication and engineering applications.

The tumour-promoting receptor tyrosine kinase, EphB4, regulates expression of Integrin-β8 in prostate cancer cells

Background The EphB4 receptor tyrosine kinase is overexpressed in many cancers including prostate cancer. The molecular mechanisms by which this ephrin receptor influences cancer progression are complex as there are tumor-promoting ligand-independent mechanisms in place as well as ligand-dependent tumor suppressive pathways. Methods We employed transient knockdown of EPHB4 in prostate cancer cells, coupled with gene microarray analysis, to identify genes that were regulated by EPHB4 and may represent linked tumor-promoting factors. We validated target genes using qRT-PCR and employed functional assays to determine their role in prostate cancer migration and invasion. Results We discovered that over 500 genes were deregulated upon EPHB4 siRNA knockdown, with integrin β8 (ITGB8) being the top hit (29-fold down-regulated compared to negative non-silencing siRNA). Gene ontology analysis found that the process of cell adhesion was highly deregulated and two other integrin genes, ITGA3 and ITGA10, were also differentially expressed. In parallel, we also discovered that over-expression of EPHB4 led to a concomitant increase in ITGB8 expression. In silico analysis of a prostate cancer progression microarray publically available in the Oncomine database showed that both EPHB4 and ITGB8 are highly expressed in prostatic intraepithelial neoplasia, the precursor to prostate cancer. Knockdown of ITGB8 in PC-3 and 22Rv1 prostate cancer cells in vitro resulted in significant reduction of cell migration and invasion. Conclusions These results reveal that EphB4 regulates integrin β8 expression and that integrin β8 plays a hitherto unrecognized role in the motility of prostate cancer cells and thus targeting integrin β8 may be a new treatment strategy for prostate cancer.

Genetic determinants of heel bone properties: Genome-wide association meta-analysis and replication in the GEFOS/GENOMOS consortium

Quantitative ultrasound of the heel captures heel bone properties that independently predict fracture risk and, with bone mineral density (BMD) assessed by X-ray (DXA), may be convenient alternatives for evaluating osteoporosis and fracture risk. We performed a meta-analysis of genome-wide association (GWA) studies to assess the genetic determinants of heel broadband ultrasound attenuation (BUA; n 5 14 260), velocity of sound (VOS; n 5 15 514) and BMD (n 5 4566) in 13 discovery cohorts. Independent replication involved seven cohorts with GWA data (in silico n 5 11 452) and new genotyping in 15 cohorts (de novo n 5 24 902). In combined random effects, meta-analysis of the discovery and replication cohorts, nine single nucleotide polymorphisms (SNPs) had genome-wide significant (P < 5 3 108) associations with heel bone properties. Alongside SNPs within or near previously identified osteoporosis susceptibility genes including ESR1 (6q25.1: rs4869739, rs3020331, rs2982552), SPTBN1 (2p16.2: rs11898505), RSPO3 (6q22.33: rs7741021), WNT16 (7q31.31: rs2908007), DKK1 (10q21.1: rs7902708) and GPATCH1 (19q13.11: rs10416265), we identified a new locus on chromosome 11q14.2 (rs597319 close to TMEM135, a gene recently linked to osteoblastogenesis and longevity) significantly associated with both BUA and VOS (P < 8.23 3 1014). In meta-analyses involving 25 cohorts with up to 14 985 fracture cases, six of 10 SNPs associated with heel bone properties at P < 5 3 106 also had the expected direction of association with any fracture (P < 0.05), including threeSNPswithP < 0.005: 6q22.33 (rs7741021), 7q31.31 (rs2908007) and 10q21.1 (rs7902708). In conclusion, thisGWAstudy reveals the effect of several genescommon to central DXA-derivedBMDand heel ultrasound/DXAmeasures and points to anewgenetic locus with potential implications for better understanding of osteoporosis pathophysiology.

Identification of IDUA and WNT16 phosphorylation-related non-synonymous polymorphisms for bone mineral density in meta-analyses of genome-wide association studies

Protein phosphorylation regulates a wide variety of cellular processes. Thus, we hypothesize that single-nucleotide polymorphisms (SNPs) that may modulate protein phosphorylation could affect osteoporosis risk. Based on a previous conventional genome-wide association (GWA) study, we conducted a three-stage meta-analysis targeting phosphorylation-related SNPs (phosSNPs) for femoral neck (FN)-bone mineral density (BMD), total hip (HIP)-BMD, and lumbar spine (LS)-BMD phenotypes. In stage 1, 9593 phosSNPs were meta-analyzed in 11,140 individuals of various ancestries. Genome-wide significance (GWS) and suggestive significance were defined by α = 5.21 × 10–6 (0.05/9593) and 1.00 × 10–4, respectively. In stage 2, nine stage 1–discovered phosSNPs (based on α = 1.00 × 10–4) were in silico meta-analyzed in Dutch, Korean, and Australian cohorts. In stage 3, four phosSNPs that replicated in stage 2 (based on α = 5.56 × 10–3, 0.05/9) were de novo genotyped in two independent cohorts. IDUA rs3755955 and rs6831280, and WNT16 rs2707466 were associated with BMD phenotypes in each respective stage, and in three stages combined, achieving GWS for both FN-BMD (p = 8.36 × 10–10, p = 5.26 × 10–10, and p = 3.01 × 10–10, respectively) and HIP-BMD (p = 3.26 × 10–6, p = 1.97 × 10–6, and p = 1.63 × 10–12, respectively). Although in vitro studies demonstrated no differences in expressions of wild-type and mutant forms of IDUA and WNT16B proteins, in silico analyses predicts that WNT16 rs2707466 directly abolishes a phosphorylation site, which could cause a deleterious effect on WNT16 protein, and that IDUA phosSNPs rs3755955 and rs6831280 could exert indirect effects on nearby phosphorylation sites. Further studies will be required to determine the detailed and specific molecular effects of these BMD-associated non-synonymous variants. © 2015 American Society for Bone and Mineral Research.

Identification of a novel FGFRL1 MicroRNA target site polymorphism for bone mineral density in meta-analyses of genome-wide association studies

MicroRNAs (miRNAs) are critical post-transcriptional regulators. Based on a previous genome-wide association (GWA) scan, we conducted a polymorphism in microRNAs' Target Sites (poly-miRTS)-centric multistage meta-analysis for lumbar spine (LS)-, total hip (HIP)-, and femoral neck (FN)-bone mineral density (BMD). In stage I, 41,102 poly-miRTSs were meta-analyzed in 7 cohorts with a genome-wide significance (GWS) α=0.05/41,102=1.22×10-6. By applying α=5×10-5 (suggestive significance), 11 poly-miRTSs were selected, with FGFRL1 rs4647940 and PRR5 rs3213550 as top signals for FN-BMD (P-value=7.67×10-6 and 1.58×10-5) in gender-combined sample. In stage II in silico replication (two cohorts), FGFRL1 rs4647940 was the only signal marginally replicated for FN-BMD (P-value=5.08×10-3) at α=0.10/11=9.09×10-3. PRR5 rs3213550 was also selected based on biological significance. In stage III de novo genotyping replication (two cohorts), FGFRL1 rs4647940 was the only signal significantly replicated for FN-BMD (P-value=7.55×10-6) at α=0.05/2=0.025 in gender-combined sample. Aggregating three stages, FGFRL1 rs4647940 was the single stage I-discovered and stages II- and III-replicated signal attaining GWS for FN-BMD (P-value=8.87×10-12). Dual-luciferase reporter assays demonstrated that FGFRL1 3' untranslated region harboring rs4647940 appears to be hsa-miR-140-5p's target site. In a zebrafish microinjection experiment, dre-miR-140-5p is shown to exert a dramatic impact on craniofacial skeleton formation. Taken together, we provided functional evidence for a novel FGFRL1 poly-miRTS rs4647940 in a previously known 4p16.3 locus, and experimental and clinical genetics studies have shown both FGFRL1 and hsa-miR-140-5p are important for bone formation. © The Author 2015. Published by Oxford University Press. All rights reserved.

Chlamydial infertility in women: Diagnosis, epidemiology and immune response

This project primarily focused on the development of a novel, non-invasive peptide-based ELISA to diagnose C. trachomatis-related infertility in women. In this study, an in silico approach was applied to identify and design novel peptide antigens that are specific to women with C. trachomatis-related infertility. The serological assay was developed such that it could potentially replace invasive techniques for early infertility investigation. The thesis further investigated the sero-prevalence of chlamydial infertility in women from a low-resource, high prevalence setting. In addition, the thesis also identified immune markers that were associated with the C. trachomatis-related infertility in women.