Standardization and validation of an induced ovulation model system in buffalo cows: Characterization of gene expression changes in the periovulatory follicle

In bovines characterization of biochemical and molecular determinants of the dominant follicle before and during different time intervals after gonadotrophin surge requires precise identification of the dominant follicle from a follicular wave. The objectives of the present study were to standardize an experimental model in buffalo cows for accurately identifying the dominant follicle of the first wave of follicular growth and characterize changes in follicular fluid hormone concentrations as well as expression patterns of various genes associated with the process of ovulation. From the day of estrus (day 0), animals were subjected to blood sampling and ultrasonography for monitoring circulating progesterone levels and follicular growth. On day 7 of the cycle, animals were administered a PGF2α analogue (Tiaprost Trometamol, 750 μg i.m.) followed by an injection of hCG (2000 IU i.m.) 36 h later. Circulating progesterone levels progressively increased from day 1 of the cycle to 2.26 ± 0.17 ng/ml on day 7 of the cycle, but declined significantly after PGF2α injection. A progressive increase in the size of the dominant follicle was observed by ultrasonography. The follicular fluid estradiol and progesterone concentrations in the dominant follicle were 600 ± 16.7 and 38 ± 7.6 ng/ml, respectively, before hCG injection and the concentration of estradiol decreased to 125.8 ± 25.26 ng/ml, but concentration of progesterone increased to 195 ± 24.6 ng/ml, 24 h post-hCG injection. Inh-α and Cyp19A1 expressions in granulosa cells were maximal in the dominant follicle and declined in response to hCG treatment. Progesterone receptor, oxytocin and cycloxygenase-2 expressions in granulosa cells, regarded as markers of ovulation, were maximal at 24 h post-hCG. The expressions of genes belonging to the super family of proteases were also examined; Cathepsin L expression decreased, while ADAMTS 3 and 5 expressions increased 24 h post-hCG treatment. The results of the current study indicate that sequential treatments of PGF2α and hCG during early estrous cycle in the buffalo cow leads to follicular growth that culminates in ovulation. The model system reported in the present study would be valuable for examining temporo-spatial changes in the periovulatory follicle immediately before and after the onset of gonadotrophin surge.

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.

HDAC inhibitor valproic acid enhances tumor cell kill in adenovirus-HSVtk mediated suicide gene therapy in HNSCC xenograft mouse model

Safety, efficacy and enhanced transgene expression are the primary concerns while using any vector for gene therapy. One of the widely used vectors in clinical. trials is adenovirus which provides a safe way to deliver the therapeutic gene. However, adenovirus has poor transduction efficiency in vivo since most tumor cells express low coxsackie and adenovirus receptors. Similarly transgene expression remains low, possibly because of the chromatization of adenoviral genome upon infection in eukaryotic cells, an effect mediated by histone deacetylases (HDACs). Using a recombinant adenovirus (Ad-HSVtk) carrying the herpes simplex thymidine kinase (HSVtk) and GFP genes we demonstrate that HDAC inhibitor valproic acid can bring about an increase in CAR expression on host cells and thereby enhanced Ad-HSVtk infectivity. It also resulted in an increase in transgene (HSVtk and GFP) expression. This, in turn, resulted in increased cell kill of HNSCC cells, following ganciclovir treatment in vitro as well as in vivo in a xenograft nude mouse model.

Drosophila indirect flight muscle specific Act88F actin mutants as a model system for studying congenital myopathies of the human ACTA1 skeletal muscle actin gene

Most human ACTA1 skeletal actin gene mutations cause dominant, congenital myopathies often with severely reduced muscle function and neonatal mortality. High sequence conservation of actin means many mutated ACTA1 residues are identical to those in the Drosophila Act88F, an indirect flight muscle specific sarcomeric actin. Four known Act88F mutations occur at the same actin residues mutated in ten ACTA1 nemaline mutations, A138D/P, R256H/L, G268C/D/R/S and R372C/S. These Act88F mutants were examined for similar muscle phenotypes. Mutant homozygotes show phenotypes ranging from a lack of myofibrils to almost normal sarcomeres at eclosion. Aberrant Z-disc-like structures and serial Z-disc arrays, ‘zebra bodies’, are observed in homozygotes and heterozygotes of all four Act88F mutants. These electron-dense structures show homologies to human nemaline bodies/rods, but are much smaller than those typically found in the human myopathy. We conclude that the Drosophila indirect flight muscles provide a good model system for studying ACTA1 mutations.

NETGEM: Network Embedded Temporal GEnerative Model for gene expression data

Background: Temporal analysis of gene expression data has been limited to identifying genes whose expression varies with time and/or correlation between genes that have similar temporal profiles. Often, the methods do not consider the underlying network constraints that connect the genes. It is becoming increasingly evident that interactions change substantially with time. Thus far, there is no systematic method to relate the temporal changes in gene expression to the dynamics of interactions between them. Information on interaction dynamics would open up possibilities for discovering new mechanisms of regulation by providing valuable insight into identifying time-sensitive interactions as well as permit studies on the effect of a genetic perturbation. Results: We present NETGEM, a tractable model rooted in Markov dynamics, for analyzing the dynamics of the interactions between proteins based on the dynamics of the expression changes of the genes that encode them. The model treats the interaction strengths as random variables which are modulated by suitable priors. This approach is necessitated by the extremely small sample size of the datasets, relative to the number of interactions. The model is amenable to a linear time algorithm for efficient inference. Using temporal gene expression data, NETGEM was successful in identifying (i) temporal interactions and determining their strength, (ii) functional categories of the actively interacting partners and (iii) dynamics of interactions in perturbed networks. Conclusions: NETGEM represents an optimal trade-off between model complexity and data requirement. It was able to deduce actively interacting genes and functional categories from temporal gene expression data. It permits inference by incorporating the information available in perturbed networks. Given that the inputs to NETGEM are only the network and the temporal variation of the nodes, this algorithm promises to have widespread applications, beyond biological systems. The source code for NETGEM is available from https://github.com/vjethava/NETGEM

Reaction dynamics under confinement: an exact path integral treatment of a two-stage model of stochastic gene expression

Gene expression in living systems is inherently stochastic, and tends to produce varying numbers of proteins over repeated cycles of transcription and translation. In this paper, an expression is derived for the steady-state protein number distribution starting from a two-stage kinetic model of the gene expression process involving p proteins and r mRNAs. The derivation is based on an exact path integral evaluation of the joint distribution, P(p, r, t), of p and r at time t, which can be expressed in terms of the coupled Langevin equations for p and r that represent the two-stage model in continuum form. The steady-state distribution of p alone, P(p), is obtained from P(p, r, t) (a bivariate Gaussian) by integrating out the r degrees of freedom and taking the limit t -> infinity. P(p) is found to be proportional to the product of a Gaussian and a complementary error function. It provides a generally satisfactory fit to simulation data on the same two-stage process when the translational efficiency (a measure of intrinsic noise levels in the system) is relatively low; it is less successful as a model of the data when the translational efficiency (and noise levels) are high.

Genome-Wide Gene Expression Analysis Reveals a Dynamic Interplay between Luteotropic and Luteolytic Factors in the Regulation of Corpus Luteum Function in the Bonnet Monkey (Macaca radiata)

Although LH is essential for survival and function of the corpus luteum (CL) in higher primates, luteolysis occurs during nonfertile cycles without a discernible decrease in circulating LH levels. Using genome-wide expression analysis, several experiments were performed to examine the processes of luteolysis and rescue of luteal function in monkeys. Induced luteolysis with GnRH receptor antagonist (Cetrorelix) resulted in differential regulation of 3949 genes, whereas replacement with exogenous LH (Cetrorelix plus LH) led to regulation of 4434 genes (1563 down-regulation and 2871 up-regulation). A model system for prostaglandin (PG) F-2 alpha-induced luteolysis in the monkey was standardized and demonstrated that PGF(2 alpha) regulated expression of 2290 genes in the CL. Analysis of the LH-regulated luteal transcriptome revealed that 120 genes were regulated in an antagonistic fashion by PGF(2 alpha). Based on the microarray data, 25 genes were selected for validation by real-time RT-PCR analysis, and expression of these genes was also examined in the CL throughout the luteal phase and from monkeys treated with human chorionic gonadotropin (hCG) to mimic early pregnancy. The results indicated changes in expression of genes favorable to PGF(2 alpha) action during the late to very late luteal phase, and expressions of many of these genes were regulated in an opposite manner by exogenous hCG treatment. Collectively, the findings suggest that curtailment of expression of downstream LH-target genes possibly through PGF(2 alpha) action on the CL is among the mechanisms underlying cross talk between the luteotropic and luteolytic signaling pathways that result in the cessation of luteal function, but hCG is likely to abrogate the PGF(2 alpha)-responsive gene expression changes resulting in luteal rescue crucial for the maintenance of early pregnancy. (Endocrinology 150: 1473-1484, 2009)

A model for mammalian male determination based on a passive Y chromosome

A model is suggested for mammalian male determination based on interactions postulated to occur among an autosomal repressor gene, an X-linked male-determining gene termed Tdx, and multiple copies of certain DNA sequences on the Y chromosome that do not code for any protein. The repressor, synthesised in limited amounts, has higher affinity for the Y-linked sequences than for Tdx and its affinity for Tdx is greater than that of RNA polymerase. In XY cells the Y effectively binds all available repressor, permitting transcription of Tdx to occur. In XX cells, since competition from the Y-linked high-affinity sequences is absent, the repressor binds to Tdx and prevents transcription. As a result of this competition between Tdx and the Y-linked high-affinity sites for limiting concentrations of the autosomal repressor, the product of the Tdx gene (TDX) is synthesized in the male but not in the female. It is suggested that in determination of the male sex, the role of the Y chromosome is to serve as a sink for the Tdx repressor. The proposed interactions provide a plausible explanation for the genetic properties of several anomalies of sexual development in mouse, man, and other mammals. The model suggests that the postulated multiple, highaffinity sequences on the Y chromosome of the mouse are included among the DNA sequences referred to as the Sxr-Bkm sequences.

Identification of key DNA elements involved in promoter recognition by Mxr1p, a master regulator of methanol utilization pathway in Pichia pastoris

Mxr1p (methanol expression regulator 1) functions as a key regulator of methanol metabolism in the methylotrophic yeast Pichia pastoris. In this study, a recombinant Mxr1p protein containing the N-terminal zinc finger DNA binding domain was overexpressed and purified from E coli cells and its ability to bind to promoter sequences of AOXI encoding alcohol oxidase was examined. In the AOXI promoter, Mxr1p binds at six different regions. Deletions encompassing these regions result in a significant decrease in AOXI promoter activity in vivo. Based on the analysis of AOXI promoter sequences, a consensus sequence for Mxr1p binding consisting of a core 5' CYCC 3' motif was identified. When the core CYCC sequence is mutated to CYCA, CYCT or CYCM (M = 5-methylcytosine), Mxr1p binding is abolished. Though Mxr1p is the homologue of Saccharomyces cerevisiae Adr1p transcription factor, it does not bind to Adr1p binding site of S. cerevisiae alcohol dehydrogenase promoter (ADH2UAS1). However, two point mutations convert ADH2UAS1 into an Mxr1p binding site. The identification of key DNA elements involved in promoter recognition by Mxr1p is an important step in understanding its function as a master regulator of the methanol utilization pathway in P. pastoris.

Positional Preferences of Polypurine/Polypyrimidine Tracts in Saccharomyces cerevisiae Genome: Implications for cis Regulation of Gene Expression

The complete genome of the baker's yeast S. cerevisiae was analyzed for the presence of polypurine/polypyrimidine (poly[pu/py]) repeats and their occurrences were classified on the basis of their location within and outside open reading frames (ORFs). The analysis reveals that such sequence motifs are present abundantly both in coding as well as noncoding regions. Clear positional preferences are seen when these tracts occur in noncoding regions. These motifs appear to occur predominantly at a unit nucleosomal length both upstream and downstream of ORFs. Moreover, there is a biased distribution of polypurines in the coding strands when these motifs occur within open reading frames. The significance of the biased distribution is discussed with reference to the occurrence of these motifs in other known mRNA sequences and expressed sequence tags. A model for cis regulation of gene expression is proposed based on the ability of these motifs to form an intermolecular triple helix structure when present within the coding region and/or to modulate nucleosome positioning via enhanced histone affinity when present outside coding regions.

Gene expression profiling during Forskolin induced differentiation of BeWo cells by differential display RT-PCR

The differentiation of cytotrophoblasts into syncytiotrophoblasts in the placenta has been employed as a model to investigate stage specific expression as well as regulation of genes during this process. While the cytotrophoblasts are highly invasive and proliferative with relatively less capacity to synthesize pregnancy related proteins, the multinucleated syncytiotrophoblasts are non-proliferative and non-invasive. However, syncytiotrophoblasts are the site of synthesis of a variety of protein, peptide and steroid hormones as well as several growth factors. Both the freshly isolated cytotrophoblasts from human placenta as well as the BeWo cell, a choriocarcinoma cell line model which retain several characteristic of cytotrophoblasts has been employed by us to study regulation of differentiation. In the present study, we have employed the differential display RT-PCR analysis (DD-RT-PCR) to evaluate gene expression changes during Forskolin induced in vitro differentiation of BeWo cells. We have identified several genes which are differentially expressed during differentiation and the differential expression of 10 transcripts was confirmed by Northern blot analysis. Based on the identity of the transcripts an attempt has been made to relate the known function of the gene products, to changes observed during differentiation. Of the several transcripts, one of the transcripts, namely Secretory Leukocyte Protease Inhibitor (SLPI) which is known to have multiple functions was found to increase 15-fold in the syntiotrophoblast.

Proteomic profiling of forskolin-induced differentiated BeWo cells: an in-vitro model of cytotrophoblast differentiation

Placental trophoblastic differentiation is characterized by the fusion of monolayer cytotrophoblasts into syncytiotrophoblasts. During this process of differentiation, several morphological and biochemical changes are known to occur, and this model has been employed to investigate the changes that occur at the gene and protein level during differentiation. Using the sensitive technique of proteomics [two-dimensional gel electrophoresis (2DGE)], changes in protein profile were evaluated in the control and forskolin-induced differentiated cells of trophoblastic choriocarcinoma BeWo cell line. Several proteins were differentially expressed in control and differentiated cells. Four major proteins were up-regulated as assessed by silver staining, and were further characterized as c-h-ras p 21 (phosphorylated), retinoblastoma susceptibility protein I and integrase interactor protein 1. These proteins are known to play an important role in growth arrest of cells, and thus may play a role in initiating the process of differentiation.

Stochastic simulations of the origins and implications of long-tailed distributions in gene expression

Gene expression noise results in protein number distributions ranging from long-tailed to Gaussian. We show how long-tailed distributions arise from a stochastic model of the constituent chemical reactions and suggest that, in conjunction with cooperative switches, they lead to more sensitive selection of a subpopulation of cells with high protein number than is possible with Gaussian distributions. Single-cell-tracking experiments are presented to validate some of the assumptions of the stochastic simulations. We also examine the effect of DNA looping on the shape of protein distributions. We further show that when switches are incorporated in the regulation of a gene via a feedback loop, the distributions can become bimodal. This might explain the bimodal distribution of certain morphogens during early embryogenesis.

Involvement of heme in the transcriptional activation of CYPIIB1/B2 gene by phenobarbitone in rat liver—Studies with succinylacetone

Earlier studies in this laboratory had implicated heme to function as a positive modulator of phenobarbitonemediated activation of CYPIIB1/B2 gene transcription in rat liver. However, recent reports have indicated that succinylacetone, a specific inhibitor of δ-aminolevulinate dehydrase, does not affect this process. The present studies indicate that succinylacetone does inhibit the phenobarbitone-mediated increase in CYPIIB1/B2 mRNAs and their transcription in rat liver at early time points (45 min to 3 h), but the inhibition is not pronounced at later time points (16 h). Succinylacetone is a weaker inhibitor of heme biosynthesis than CoCl2, 3-amino-1,2,4-triazole, or thioacetamide used earlier in this laboratory. Succinylacetone induces δ-aminolevulinate synthase, whereas the other compounds depress the levels of the enzyme. There is a good correlation between the amount of freshly synthesized nuclear heme pool and the activation of CYPIIB1/B2 transcription by phenobarbitone. A model implicating a nuclear heme pool regulating the transcription of δ-aminolevulinate synthase, CYPIIB1/ B2, and heme oxygenase genes is proposed.

Uncoupling Protein 2 and 3 Gene Polymorphisms and Their Association with Type 2 Diabetes in Asian Indians

Background: This study examined the association of -866G/A, Ala55Val, 45bpI/D, and -55C/T polymorphisms at the uncoupling protein (UCP) 3-2 loci with type 2 diabetes in Asian Indians. Methods: A case-control study was performed among 1,406 unrelated subjects (487 with type 2 diabetes and 919 normal glucose-tolerant NGT]), chosen from the Chennai Urban Rural Epidemiology Study, an ongoing population-based study in Southern India. The polymorphisms were genotyped using polymerase chain reaction-restriction fragment length polymorphism and direct sequencing. Haplotype frequencies were estimated using an expectation-maximization algorithm. Linkage disequilibrium was estimated from the estimates of haplotypic frequencies. Results: The genotype (P = 0.00006) and the allele (P = 0.00007) frequencies of Ala55Val of the UCP2 gene showed a significant protective effect against the development of type 2 diabetes. The odds ratios (adjusted for age, sex, and body mass index) for diabetes for individuals carrying Ala/Val was 0.72, and that for individuals carrying Val/Val was 0.37. Homeostasis insulin resistance model assessment and 2-h plasma glucose were significantly lower among Val-allele carriers compared to the Ala/Ala genotype within the NGT group. The genotype (P = 0.02) and the allele (P = 0.002) frequencies of -55C/T of the UCP3 gene showed a significant protective effect against the development of diabetes. The odds ratio for diabetes for individuals carrying CT was 0.79, and that for individuals carrying TT was 0.61. The haplotype analyses further confirmed the association of Ala55Val with diabetes, where the haplotypes carrying the Ala allele were significantly higher in the cases compared to controls. Conclusions: Ala55Val and -55C/T polymorphisms at the UCP3-2 loci are associated with a significantly reduced risk of developing type 2 diabetes in Asian Indians.