Oestrogen-modulated gene expression in the human endometrium

To identify key regulatory mechanisms in the growth and development of the human endometrium, microarray analysis was performed on uncultured human endometrium collected during menstruation (M) and the late-proliferative (LATE-P)-phase of the menstrual cycle, as well as after 24 h incubation in the presence of oestradiol (17beta-E2). We demonstrate the expression of novel gene transcripts in the human endometrium. i.e. mucin-9, novel oestrogen-responsive gene transcripts, i.e. gelsolin and flotillin-1, and genes known to be expressed in human endometrium but not yet shown to be oestrogen responsive, i.e. connexin-37 and TFF1/pS2. Genes reported to be expressed during the implantation window and implicated in progesterone action, i.e. secretoglobin family 2A, member 2 (mammaglobin) and homeobox-containing proteins, were up-regulated in uncultured LATE-P-phase endometrium compared to M-phase endometrium. Some gene transcripts are regulated directly by 17beta-E2 alone, others are influenced by the in vivo environment as well. These observations emphasise that the regulation of endometrium maturation by oestrogen entails more then just stimulation of cell proliferation.

Choose wisely: Network, ontology and annotation resources for the analysis of Staphylococcus aureus omics data

Staphylococcus aureus (S. aureus) is a prominent human and livestock pathogen investigated widely using omic technologies. Critically, due to availability, low visibility or scattered resources, robust network and statistical contextualisation of the resulting data is generally under-represented. Here, we present novel meta-analyses of freely-accessible molecular network and gene ontology annotation information resources for S. aureus omics data interpretation. Furthermore, through the application of the gene ontology annotation resources we demonstrate their value and ability (or lack-there-of) to summarise and statistically interpret the emergent properties of gene expression and protein abundance changes using publically available data. This analysis provides simple metrics for network selection and demonstrates the availability and impact that gene ontology annotation selection can have on the contextualisation of bacterial omics data.

Genome-wide scan of healthy human connectome discovers SPON1 gene variant influencing dementia severity

Aberrant connectivity is implicated in many neurological and psychiatric disorders, including Alzheimer's disease and schizophrenia. However, other than a few disease-associated candidate genes, we know little about the degree to which genetics play a role in the brain networks; we know even less about specific genes that influence brain connections. Twin and family-based studies can generate estimates of overall genetic influences on a trait, but genome-wide association scans (GWASs) can screen the genome for specific variants influencing the brain or risk for disease. To identify the heritability of various brain connections, we scanned healthy young adult twins with high-field, highangular resolution diffusion MRI. We adapted GWASs to screen the brain's connectivity pattern, allowing us to discover genetic variants that affect the human brain's wiring. The association of connectivity with the SPON1 variant at rs2618516 on chromosome 11 (11p15.2) reached connectome-wide, genome-wide significance after stringent statistical corrections were enforced, and it was replicated in an independent subsample. rs2618516 was shown to affect brain structure in an elderly population with varying degrees of dementia. Older people who carried the connectivity variant had significantly milder clinical dementia scores and lower risk of Alzheimer's disease. As a posthoc analysis, we conducted GWASs on several organizational and topological network measures derived from the matrices to discover variants in and around genes associated with autism (MACROD2), development (NEDD4), and mental retardation (UBE2A) significantly associated with connectivity. Connectome-wide, genome-wide screening offers substantial promise to discover genes affecting brain connectivity and risk for brain diseases.

Gene regulation by translational inhibition is determined by Dicer partnering proteins

MicroRNAs (miRNAs) are small regulatory RNAs produced by Dicer proteins that regulate gene expression in development and adaptive responses to the environment1,​2,​3,​4. In animals, the degree of base pairing between a miRNA and its target messenger RNA seems to determine whether the regulation occurs through cleavage or translation inhibition1. In contrast, the selection of regulatory mechanisms is independent of the degree of mismatch between a plant miRNA and its target transcript5. However, the components and mechanism(s) that determine whether a plant miRNA ultimately regulates its targets by guiding cleavage or translational inhibition are unknown6. Here we show that the form of regulatory action directed by a plant miRNA is determined by DRB2, a DICER-LIKE1 (DCL1) partnering protein. The dependence of DCL1 on DRB1 for miRNA biogenesis is well characterized7,​8,​9, but we show that it is only required for miRNA-guided transcript cleavage. We found that DRB2 determines miRNA-guided translational inhibition and represses DRB1 expression, thereby allowing the active selection of miRNA regulatory action. Furthermore, our results reveal that the core silencing proteins ARGONAUTE1 (AGO1) and SERRATE (SE) are highly regulated by miRNA-guided translational inhibition. DRB2 has been remarkably conserved throughout plant evolution, raising the possibility that translational repression is the ancient form of miRNA-directed gene regulation in plants, and that Dicer partnering proteins, such as human TRBP, might play a similar role in other eukaryotic systems.

Ankylosing spondylitis is associated with the anthrax toxin receptor 2 gene (ANTXR2)

Objectives ANTXR2 variants have been associated with ankylosing spondylitis (AS) in two previous genome-wide association studies (GWAS) (p∼9×10-8). However, a genome-wide significant association (p<5×10-8) was not observed. We conducted a more comprehensive analysis of ANTXR2 in an independent UK sample to confirm and refine this association. Methods A replication study was carried out with 2978 cases and 8365 controls. Then, these were combined with non-overlapping samples from the two previous GWAS in a meta-analysis. Human leukocyte antigen (HLA)-B27 stratification was also performed to test for ANTXR2-HLA-B27 interaction. Results Out of nine single nucleotide polymorphisms (SNP) in the study, five SNPs were nominally associated (p<0.05) with AS in the replication dataset. In the meta-analysis, eight SNPs showed evidence of association, the strongest being with rs12504282 (OR=0.88, p=6.7×10-9). Seven of these SNPs showed evidence for association in the HLA-B27-positive subgroup, but none was associated with HLA-B27-negative AS. However, no statistically significant interaction was detected between HLA-B27 and ANTXR2 variants. Conclusions ANTXR2 variants are clearly associated with AS. The top SNPs from two previous GWAS (rs4333130 and rs4389526) and this study (rs12504282) are in strong linkage disequilibrium (r2≥0.76). All are located near a putative regulatory region. Further studies are required to clarify the role played by these ANTXR2 variants in AS.

Mapping the regulatory sequences controlling 93 breast cancer-associated miRNA genes leads to the identification of two functional promoters of the Hsa-mir-200b cluster, methylation of which is associated with metastasis or hormone receptor status in advanced breast cancer

MicroRNAs (miRNAs) are small non-coding RNAs of 20 nt in length that are capable of modulating gene expression post-transcriptionally. Although miRNAs have been implicated in cancer, including breast cancer, the regulation of miRNA transcription and the role of defects in this process in cancer is not well understood. In this study we have mapped the promoters of 93 breast cancer-associated miRNAs, and then looked for associations between DNA methylation of 15 of these promoters and miRNA expression in breast cancer cells. The miRNA promoters with clearest association between DNA methylation and expression included a previously described and a novel promoter of the Hsa-mir-200b cluster. The novel promoter of the Hsa-mir-200b cluster, denoted P2, is located 2 kb upstream of the 5′ stemloop and maps within a CpG island. P2 has comparable promoter activity to the previously reported promoter (P1), and is able to drive the expression of miR-200b in its endogenous genomic context. DNA methylation of both P1 and P2 was inversely associated with miR-200b expression in eight out of nine breast cancer cell lines, and in vitro methylation of both promoters repressed their activity in reporter assays. In clinical samples, P1 and P2 were differentially methylated with methylation inversely associated with miR-200b expression. P1 was hypermethylated in metastatic lymph nodes compared with matched primary breast tumours whereas P2 hypermethylation was associated with loss of either oestrogen receptor or progesterone receptor. Hypomethylation of P2 was associated with gain of HER2 and androgen receptor expression. These data suggest an association between miR-200b regulation and breast cancer subtype and a potential use of DNA methylation of miRNA promoters as a component of a suite of breast cancer biomarkers.

Role of TATATAA element in the regulation of tRNA(1)(Gly) gene expression in Bombyx mori is position dependent

Transcription of tRNA genes by RNA polymerase III is controlled by the internal conserved sequences within the coding region and the immediate upstream flanking sequences. A highly transcribed copy of glycyl tRNA gene tRNA1(Gly)-1 from Bombyx mori is down regulated by sequences located much farther upstream in the region -150 to -300 nucleotides (nt), with respect to the +1 nt of tRNA. The negative regulatory effect has been narrowed down to a sequence motif 'TATATAA', a perfect consensus recognised by the TATA binding protein, TBP. This sequence element, when brought closer to the transcription start point, on the other hand, exerts a positive effect by promoting transcription of the gene devoid of other cis regulatory elements. The identity of the nuclear protein interacting with this 'TATATAA' element to TBP has been established by antibody and mutagenesis studies. The 'TATATAA' element thus influences the transcription of tRNA genes positively or negatively in a position-dependent manner either by recruitment or sequestration of TBP from the transcription machinery.

Identification of putative regulatory motifs in the upstream regions of co-expressed functional groups of genes in Plasmodium falciparum

Background: Regulation of gene expression in Plasmodium falciparum (Pf) remains poorly understood. While over half the genes are estimated to be regulated at the transcriptional level, few regulatory motifs and transcription regulators have been found. Results: The study seeks to identify putative regulatory motifs in the upstream regions of 13 functional groups of genes expressed in the intraerythrocytic developmental cycle of Pf. Three motif-discovery programs were used for the purpose, and motifs were searched for only on the gene coding strand. Four motifs – the 'G-rich', the 'C-rich', the 'TGTG' and the 'CACA' motifs – were identified, and zero to all four of these occur in the 13 sets of upstream regions. The 'CACA motif' was absent in functional groups expressed during the ring to early trophozoite transition. For functional groups expressed in each transition, the motifs tended to be similar. Upstream motifs in some functional groups showed 'positional conservation' by occurring at similar positions relative to the translational start site (TLS); this increases their significance as regulatory motifs. In the ribonucleotide synthesis, mitochondrial, proteasome and organellar translation machinery genes, G-rich, C-rich, CACA and TGTG motifs, respectively, occur with striking positional conservation. In the organellar translation machinery group, G-rich motifs occur close to the TLS. The same motifs were sometimes identified for multiple functional groups; differences in location and abundance of the motifs appear to ensure different modes of action. Conclusion: The identification of positionally conserved over-represented upstream motifs throws light on putative regulatory elements for transcription in Pf.

Flux balance analysis of biological systems: applications and challenges

Systems level modelling and simulations of biological processes are proving to be invaluable in obtaining a quantitative and dynamic perspective of various aspects of cellular function. In particular, constraint-based analyses of metabolic networks have gained considerable popularity for simulating cellular metabolism, of which flux balance analysis (FBA), is most widely used. Unlike mechanistic simulations that depend on accurate kinetic data, which are scarcely available, FBA is based on the principle of conservation of mass in a network, which utilizes the stoichiometric matrix and a biologically relevant objective function to identify optimal reaction flux distributions. FBA has been used to analyse genome-scale reconstructions of several organisms; it has also been used to analyse the effect of perturbations, such as gene deletions or drug inhibitions in silico. This article reviews the usefulness of FBA as a tool for gaining biological insights, advances in methodology enabling integration of regulatory information and thermodynamic constraints, and finally addresses the challenges that lie ahead. Various use scenarios and biological insights obtained from FBA, and applications in fields such metabolic engineering and drug target identification, are also discussed. Genome-scale constraint-based models have an immense potential for building and testing hypotheses, as well as to guide experimentation.

Evidence for extensive population structure in the white-spotted eagle ray within the Indo-Pacific inferred from mitochondrial gene sequences

The white-spotted eagle ray Aetobatus narinari is a species complex that occurs circumglobally throughout warm-temperate waters. Aetobatus narinari is semi-pelagic and large (up to 300 cm disc width), suggesting high dispersal capabilities and gene flow on a wide spatial scale. Sequence data from two mitochondrial genes, cytochrome b (cytb) and NADH dehydrogenase subunit 4 (ND4), were used to determine the genetic variability within and among 18 sampling locations in the central Indo-Pacific biogeographical region. Populations in the Indo-Pacific were highly genetically structured with c. 70% of the total genetic variation found among three geographical regions (East China Sea, Southeast Asia and Australia). FST was 0.64 for cytb and 0.53 for ND4, with φST values being even larger, that is, 0.78 for cytb and 0.65 for ND4. This high-level genetic partitioning provides strong evidence against extensive gene flow in A. narinari. The degree of genetic population structuring in the Indo-Pacific was similar to that found on a global scale. Global FST was 0.63 for cytb and 0.57 for ND4, and global φST values were 0.94 for cytb and 0.82 for ND4. This suggests that the A. narinari complex may be more speciose than the two or three species proposed to date. Further sampling and genetic analyses are likely to uncover the ‘evolutionarily significant’ and ‘management’ units that are critical to determine the susceptibilities of individual populations to regional fishing pressures and to provide advice on management options. Network analyses showed a close genetic relationship between haplotypes from the central Indo-Pacific and South Africa, providing support for a proposed dispersal pathway from the possible centre of origin of the A. narinari species complex in the Indo-Pacific into the Atlantic Ocean.

Development of regulatory T cells in the human thymus

The role of the immune system is to protect an organism against pathogens while maintaining tolerance against self. T cells are an essential component of the immune system and they develop in the thymus. The AIRE (autoimmune regulator) gene product plays an important role in T cell development, as it promotes expression of peripheral tissue antigens in the thymus. Developing T cells, thymocytes, which recognize self-antigens with high affinity are deleted. However, this deletion process is not perfect and not all autoreactive T cells are destroyed. When the distinction between self and non-self fails, tolerance breaks and the immune system attacks the host s own tissues. This results in autoimmunity. Regulatory T cells contribute to the maintenance of self-tolerance. They can actively suppress the function of autoreactive cells. Several populations of cells with regulatory properties have been described, but the best characterized population is the natural regulatory T cells (Treg cells), which develop in the thymus and express the transcription factor FOXP3. The thymic development of Treg cells in humans is the subject of this thesis. Thymocytes at different developmental stages were analyzed using flow cytometry. The CD4-CD8- double-negative (DN) thymocytes are the earliest T cell precursors in the T cell lineage. My results show that the Treg cell marker FOXP3 is up-regulated already in a subset of these DN thymocytes. FOXP3+ cells were also found among the more mature CD4+CD8+ double-positive (DP) cells and among the CD4+ and CD8+ single-positive (SP) thymocytes. The different developmental stages of the FOXP3+ thymocytes were isolated and their gene expression examined by quantitative PCR. T cell receptor (TCR) repertoire analysis was used to compare these different thymocyte populations. My data show that in humans commitment to the Treg cell lineage is an early event and suggest that the development of Treg cells follows a linear developmental pathway, FOXP3+ DN precursors evolving through the DP stage to become mature CD4+ Treg cells. Most T cells have only one kind of TCR on their cell surface, but a small fraction of cells expresses two different TCRs. My results show that the expression of two different TCRs is enriched among Treg cells. Furthermore, both receptors were capable of transmitting signals when bound by a ligand. By extrapolating flow cytometric data, it was estimated that the majority of peripheral blood Treg cells are indeed dual-specific. The high frequency of dual-specific cells among human Treg cells suggests that dual-specificity has a role in directing these cells to the Treg cell lineage. It is known that both genetic predisposition and environmental factors influence the development of autoimmunity. It is also known that the dysfunction or absence of Treg cells leads to the development of autoimmune manifestations. APECED (autoimmune polyendocrinopathy-candidiasis-ectodermal dystrophy) is a rare monogenic autoimmune disease, caused by mutations in the AIRE gene. In the absence of AIRE gene product, deletion of self-specific T cells is presumably disturbed and autoreactive T cells escape to the periphery. I examined whether Treg cells are also affected in APECED. I found that the frequency of FOXP3+ Treg cells and the level of FOXP3 expression were significantly lower in APECED patients than in controls. Additionally, when studied in cell cultures, the suppressive capacity of the patients' Treg cells was impaired. Additionally, repertoire analysis showed that the TCR repertoire of Treg cells was altered. These results suggest that AIRE contributes to the development of Treg cells in humans and the selection of Treg cells is impaired in APECED patients. In conclusion, my thesis elucidates the developmental pathway of Treg cells in humans. The differentiation of Tregs begins early during thymic development and both the cells dual-specificity and AIRE probably affect the final commitment of Treg cells.

A versatile gene-based test for genome-wide association studies

We have derived a versatile gene-based test for genome-wide association studies (GWAS). Our approach, called VEGAS (versatile gene-based association study), is applicable to all GWAS designs, including family-based GWAS, meta-analyses of GWAS on the basis of summary data, and DNA-pooling-based GWAS, where existing approaches based on permutation are not possible, as well as singleton data, where they are. The test incorporates information from a full set of markers (or a defined subset) within a gene and accounts for linkage disequilibrium between markers by using simulations from the multivariate normal distribution. We show that for an association study using singletons, our approach produces results equivalent to those obtained via permutation in a fraction of the computation time. We demonstrate proof-of-principle by using the gene-based test to replicate several genes known to be associated on the basis of results from a family-based GWAS for height in 11,536 individuals and a DNA-pooling-based GWAS for melanoma in approximately 1300 cases and controls. Our method has the potential to identify novel associated genes; provide a basis for selecting SNPs for replication; and be directly used in network (pathway) approaches that require per-gene association test statistics. We have implemented the approach in both an easy-to-use web interface, which only requires the uploading of markers with their association p-values, and a separate downloadable application.

Association study of the dystrobrevin-binding gene with schizophrenia in Australian and Indian samples

Numerous studies have reported association between variants in the dystrobrevin binding protein 1 (dysbindin) gene (DTNBP1) and schizophrenia. However, the pattern of results is complex and to date, no specific risk marker or haplotype has been consistently identified. The number of single nucleotide polymorphisms (SNPs) tested in these studies has ranged from 5 to 20. We attempted to replicate previous findings by testing 16 SNPs in samples of 41 Australian pedigrees, 194 Australian cases and 180 controls, and 197 Indian pedigrees. No globally significant evidence for association was observed in any sample, despite power calculations indicating sufficient power to replicate several previous findings. Possible explanations for our results include sample differences in background linkage disequilibrium and/or risk allele effect size, the presence of multiple risk alleles upon different haplotypes, or the presence of a single risk allele upon multiple haplotypes. Some previous associations may also represent false positives. Examination of Caucasian HapMap phase II genotype data spanning the DTNBP1 region indicates upwards of 40 SNPs are required to satisfactorily assess all nonredundant variation within DTNBP1 and its potential regulatory regions for association with schizophrenia. More comprehensive studies in multiple samples will be required to determine whether specific DTNBP1 variants function as risk factors for schizophrenia.

trans gene regulation in adaptive evolution: A genetic algorithm model

This is a continuation of earlier studies on the evolution of infinite populations of haploid genotypes within a genetic algorithm framework. We had previously explored the evolutionary consequences of the existence of indeterminate—“plastic”—loci, where a plastic locus had a finite probability in each generation of functioning (being switched “on”) or not functioning (being switched “off”). The relative probabilities of the two outcomes were assigned on a stochastic basis. The present paper examines what happens when the transition probabilities are biased by the presence of regulatory genes. We find that under certain conditions regulatory genes can improve the adaptation of the population and speed up the rate of evolution (on occasion at the cost of lowering the degree of adaptation). Also, the existence of regulatory loci potentiates selection in favour of plasticity. There is a synergistic effect of regulatory genes on plastic alleles: the frequency of such alleles increases when regulatory loci are present. Thus, phenotypic selection alone can be a potentiating factor in a favour of better adaptation.

Identification and functional characterization of a cis-acting positive DNA element regulating CYP 2B1/B2 gene transcription in rat liver

A positive cis-acting DNA element in the near 5'-upstream region of the CYP2B1/B2 genes in rat liver was found to play an important role in the transcription of these genes. An oligonucleotide covering -69 to -98 nt mimicked the gel mobility shift pattern given by the fragment -179 to +29 nt, which was earlier found adequate to confer the regulatory features of this gene. Two major complexes were seen, of which the slower and faster moving complexes became intense under uninduced and Phenobarbitone-induced conditions respectively. Minigene cloned DNA plasmid covering -179 to +181 nt in pUC 19 and Bal 31 mutants derived from this parent were transcribed in whole nuclei and cell free transcription extracts and mutants containing only upto -75 nt of the upstream were poorly transcribed. Transcription extracts from phenobarbitone-injected rat liver nuclei were significantly more active than extracts from uninduced rats in transcribing the minigene constructs. Addition of the oligonucleotide (-69 to -98nt) specifically inhibited the transcription of the minigene construct (-179 to +181 nt) in the cell free transcription system. It is therefore, concluded that the region -69 to -98 nt acts as a positive cis-acting element in the transcription of the CYP2B1/B2 genes and in mediating the inductive effects of phenobarbitone.