Receptor-specific in vivo desensitization by the G protein-coupled receptor kinase-5 in transgenic mice.

Transgenic mice were generated with cardiac-specific overexpression of the G protein-coupled receptor kinase-5 (GRK5), a serine/threonine kinase most abundantly expressed in the heart compared with other tissues. Animals overexpressing GRK5 showed marked beta-adrenergic receptor desensitization in both the anesthetized and conscious state compared with nontransgenic control mice, while the contractile response to angiotensin II receptor stimulation was unchanged. In contrast, the angiotensin II-induced rise in contractility was significantly attenuated in transgenic mice overexpressing the beta-adrenergic receptor kinase-1, another member of the GRK family. These data suggest that myocardial overexpression of GRK5 results in selective uncoupling of G protein-coupled receptors and demonstrate that receptor specificity of the GRKs may be important in determining the physiological phenotype.

Myocardial expression of a constitutively active alpha 1B-adrenergic receptor in transgenic mice induces cardiac hypertrophy.

Transgenic mice were generated by using the alpha-myosin heavy chain promoter coupled to the coding sequence of a constitutively active mutant alpha 1B-adrenergic receptor (AR). These transgenic animals demonstrated cardiac-specific expression of this alpha 1-AR with resultant activation of phospholipase C as shown by increased myocardial diacylglycerol content. A phenotype consistent with cardiac hypertrophy developed in adult transgenic mice with increased heart/body weight ratios, myocyte cross-sectional areas, and ventricular atrial natriuretic factor mRNA levels relative to nontransgenic controls. These transgenic animals may provide insight into the biochemical triggers that induce hypertrophy in cardiac disease and serve as a convenient experimental model for studies of this condition.

The receptor kinase family: primary structure of rhodopsin kinase reveals similarities to the beta-adrenergic receptor kinase.

Light-dependent deactivation of rhodopsin as well as homologous desensitization of beta-adrenergic receptors involves receptor phosphorylation that is mediated by the highly specific protein kinases rhodopsin kinase (RK) and beta-adrenergic receptor kinase (beta ARK), respectively. We report here the cloning of a complementary DNA for RK. The deduced amino acid sequence shows a high degree of homology to beta ARK. In a phylogenetic tree constructed by comparing the catalytic domains of several protein kinases, RK and beta ARK are located on a branch close to, but separate from the cyclic nucleotide-dependent protein kinase and protein kinase C subfamilies. From the common structural features we conclude that both RK and beta ARK are members of a newly delineated gene family of guanine nucleotide-binding protein (G protein)-coupled receptor kinases that may function in diverse pathways to regulate the function of such receptors.

cAMP stimulates transcription of the beta 2-adrenergic receptor gene in response to short-term agonist exposure.

In addition to conveying cellular responses to an effector molecule, receptors are often themselves regulated by their effectors. We have demonstrated that epinephrine modulates both the rate of transcription of the beta 2-adrenergic receptor (beta 2AR) gene and the steady-state level of beta 2AR mRNA in DDT1MF-2 cells. Short-term (30 min) exposure to epinephrine (100 nM) stimulates the rate of beta 2AR gene transcription, resulting in a 3- to 4-fold increase in steady-state beta 2AR mRNA levels. These effects are mimicked by 1 mM N6,O2'-dibutyryladenosine 3',5'-cyclic monophosphate (Bt2cAMP) or foskolin but not by phorbol esters. The half-life of the beta 2AR mRNA after addition of actinomycin D (46.7 +/- 10.2 min; mean +/- SEM; n = 5) remained unchanged after 30 min of epinephrine treatment (46.8 +/- 10.6 min; mean +/- SEM; n = 4), indicating that a change in transcription rate is the predominant factor responsible for the increase of beta 2AR mRNA. Whereas brief exposure to epinephrine or Bt2cAMP does not significantly affect the total number of cellular beta 2ARs (assessed by ligand binding), continued exposure results in a gradual decline in beta 2AR number to approximately 20% (epinephrine) or approximately 45% (Bt2cAMP) of the levels in control cells by 24 hr. Similar decreases in agonist-stimulated adenylyl cyclase activity are observed. This loss of receptors with prolonged agonist exposure is accompanied by a 50% reduction in beta 2AR mRNA. Transfection of the beta 2AR promoter region cloned onto a reporter gene (bacterial chloramphenicol acetyltransferase) allowed demonstration of a 2- to 4-fold induction of transcription by agents that elevate cAMP levels, such as forskolin or phosphodiesterase inhibitors. These results establish the presence of elements within the proximal promoter region of the beta 2AR gene responsible for the transcriptional enhancing activity of cAMP and demonstrate that beta 2AR gene expression is regulated by a type of feedback mechanism involving the second messenger cAMP.

WNT3 is a biomarker capable of predicting the definitive endoderm differentiation potential of hESCs.

Generation of functional cells from human pluripotent stem cells (PSCs) through in vitro differentiation is a promising approach for drug screening and cell therapy. However, the observed large and unavoidable variation in the differentiation potential of different human embryonic stem cell (hESC)/induced PSC (iPSC) lines makes the selection of an appropriate cell line for the differentiation of a particular cell lineage difficult. Here, we report identification of WNT3 as a biomarker capable of predicting definitive endoderm (DE) differentiation potential of hESCs. We show that the mRNA level of WNT3 in hESCs correlates with their DE differentiation efficiency. In addition, manipulations of hESCs through WNT3 knockdown or overexpression can respectively inhibit or promote DE differentiation in a WNT3 level-dependent manner. Finally, analysis of several hESC lines based on their WNT3 expression levels allowed accurate prediction of their DE differentiation potential. Collectively, our study supports the notion that WNT3 can serve as a biomarker for predicting DE differentiation potential of hESCs.

Human-chimpanzee differences in a FZD8 enhancer alter cell-cycle dynamics in the developing neocortex.

The human neocortex differs from that of other great apes in several notable regards, including altered cell cycle, prolonged corticogenesis, and increased size [1-5]. Although these evolutionary changes most likely contributed to the origin of distinctively human cognitive faculties, their genetic basis remains almost entirely unknown. Highly conserved non-coding regions showing rapid sequence changes along the human lineage are candidate loci for the development and evolution of uniquely human traits. Several studies have identified human-accelerated enhancers [6-14], but none have linked an expression difference to a specific organismal trait. Here we report the discovery of a human-accelerated regulatory enhancer (HARE5) of FZD8, a receptor of the Wnt pathway implicated in brain development and size [15, 16]. Using transgenic mice, we demonstrate dramatic differences in human and chimpanzee HARE5 activity, with human HARE5 driving early and robust expression at the onset of corticogenesis. Similar to HARE5 activity, FZD8 is expressed in neural progenitors of the developing neocortex [17-19]. Chromosome conformation capture assays reveal that HARE5 physically and specifically contacts the core Fzd8 promoter in the mouse embryonic neocortex. To assess the phenotypic consequences of HARE5 activity, we generated transgenic mice in which Fzd8 expression is under control of orthologous enhancers (Pt-HARE5::Fzd8 and Hs-HARE5::Fzd8). In comparison to Pt-HARE5::Fzd8, Hs-HARE5::Fzd8 mice showed marked acceleration of neural progenitor cell cycle and increased brain size. Changes in HARE5 function unique to humans thus alter the cell-cycle dynamics of a critical population of stem cells during corticogenesis and may underlie some distinctive anatomical features of the human brain.

Nutritional control of mRNA isoform expression during developmental arrest and recovery in C. elegans.

Nutrient availability profoundly influences gene expression. Many animal genes encode multiple transcript isoforms, yet the effect of nutrient availability on transcript isoform expression has not been studied in genome-wide fashion. When Caenorhabditis elegans larvae hatch without food, they arrest development in the first larval stage (L1 arrest). Starved larvae can survive L1 arrest for weeks, but growth and post-embryonic development are rapidly initiated in response to feeding. We used RNA-seq to characterize the transcriptome during L1 arrest and over time after feeding. Twenty-seven percent of detectable protein-coding genes were differentially expressed during recovery from L1 arrest, with the majority of changes initiating within the first hour, demonstrating widespread, acute effects of nutrient availability on gene expression. We used two independent approaches to track expression of individual exons and mRNA isoforms, and we connected changes in expression to functional consequences by mining a variety of databases. These two approaches identified an overlapping set of genes with alternative isoform expression, and they converged on common functional patterns. Genes affecting mRNA splicing and translation are regulated by alternative isoform expression, revealing post-transcriptional consequences of nutrient availability on gene regulation. We also found that phosphorylation sites are often alternatively expressed, revealing a common mode by which alternative isoform expression modifies protein function and signal transduction. Our results detail rich changes in C. elegans gene expression as larvae initiate growth and post-embryonic development, and they provide an excellent resource for ongoing investigation of transcriptional regulation and developmental physiology.

Sensitive and precise quantification of insulin-like mRNA expression in Caenorhabditis elegans.

Insulin-like signaling regulates developmental arrest, stress resistance and lifespan in the nematode Caenorhabditis elegans. However, the genome encodes 40 insulin-like peptides, and the regulation and function of individual peptides is largely uncharacterized. We used the nCounter platform to measure mRNA expression of all 40 insulin-like peptides as well as the insulin-like receptor daf-2, its transcriptional effector daf-16, and the daf-16 target gene sod-3. We validated the platform using 53 RNA samples previously characterized by high density oligonucleotide microarray analysis. For this set of genes and the standard nCounter protocol, sensitivity and precision were comparable between the two platforms. We optimized conditions of the nCounter assay by varying the mass of total RNA used for hybridization, thereby increasing sensitivity up to 50-fold and reducing the median coefficient of variation as much as 4-fold. We used deletion mutants to demonstrate specificity of the assay, and we used optimized conditions to assay insulin-like gene expression throughout the C. elegans life cycle. We detected expression for nearly all insulin-like genes and find that they are expressed in a variety of distinct patterns suggesting complexity of regulation and specificity of function. We identified insulin-like genes that are specifically expressed during developmental arrest, larval development, adulthood and embryogenesis. These results demonstrate that the nCounter platform provides a powerful approach to analyzing insulin-like gene expression dynamics, and they suggest hypotheses about the function of individual insulin-like genes.

Polymorphisms in the kinesin-like factor 1 B gene and risk of epithelial ovarian cancer in Eastern Chinese women.

The kinesin-like factor 1 B (KIF1B) gene plays an important role in the process of apoptosis and the transformation and progression of malignant cells. Genetic variations in KIF1B may contribute to risk of epithelial ovarian cancer (EOC). In this study of 1,324 EOC patients and 1,386 cancer-free female controls, we investigated associations between two potentially functional single nucleotide polymorphisms in KIF1B and EOC risk by the conditional logistic regression analysis. General linear regression model was used to evaluate the correlation between the number of variant alleles and KIF1B mRNA expression levels. We found that the rs17401966 variant AG/GG genotypes were significantly associated with a decreased risk of EOC (adjusted odds ratio (OR) = 0.81, 95 % confidence interval (CI) = 0.68-0.97), compared with the AA genotype, but no associations were observed for rs1002076. Women who carried both rs17401966 AG/GG and rs1002076 AG/AA genotypes of KIF1B had a 0.82-fold decreased risk (adjusted 95 % CI = 0.69-0.97), compared with others. Additionally, there was no evidence of possible interactions between about-mentioned co-variants. Further genotype-phenotype correlation analysis indicated that the number of rs17401966 variant G allele was significantly associated with KIF1B mRNA expression levels (P for GLM = 0.003 and 0.001 in all and Chinese subjects, respectively), with GG carriers having the lowest level of KIF1B mRNA expression. Taken together, the rs17401966 polymorphism likely regulates KIF1B mRNA expression and thus may be associated with EOC risk in Eastern Chinese women. Larger, independent studies are warranted to validate our findings.

Pax3 expression enhances PDGF-B-induced brainstem gliomagenesis and characterizes a subset of brainstem glioma.

High-grade Brainstem Glioma (BSG), also known as Diffuse Intrinsic Pontine Glioma (DIPG), is an incurable pediatric brain cancer. Increasing evidence supports the existence of regional differences in gliomagenesis such that BSG is considered a distinct disease from glioma of the cerebral cortex (CG). In an effort to elucidate unique characteristics of BSG, we conducted expression analysis of mouse PDGF-B-driven BSG and CG initiated in Nestin progenitor cells and identified a short list of expression changes specific to the brainstem gliomagenesis process, including abnormal upregulation of paired box 3 (Pax3). In the neonatal mouse brain, Pax3 expression marks a subset of brainstem progenitor cells, while it is absent from the cerebral cortex, mirroring its regional expression in glioma. Ectopic expression of Pax3 in normal brainstem progenitors in vitro shows that Pax3 inhibits apoptosis. Pax3-induced inhibition of apoptosis is p53-dependent, however, and in the absence of p53, Pax3 promotes proliferation of brainstem progenitors. In vivo, Pax3 enhances PDGF-B-driven gliomagenesis by shortening tumor latency and increasing tumor penetrance and grade, in a region-specific manner, while loss of Pax3 function extends survival of PDGF-B-driven;p53-deficient BSG-bearing mice by 33%. Importantly, Pax3 is regionally expressed in human glioma as well, with high PAX3 mRNA characterizing 40% of human BSG, revealing a subset of tumors that significantly associates with PDGFRA alterations, amplifications of cell cycle regulatory genes, and is exclusive of ACVR1 mutations. Collectively, these data suggest that regional Pax3 expression not only marks a novel subset of BSG but also contributes to PDGF-B-induced brainstem gliomagenesis.

Analysis of the genome and transcriptome of Cryptococcus neoformans var. grubii reveals complex RNA expression and microevolution leading to virulence attenuation.

Cryptococcus neoformans is a pathogenic basidiomycetous yeast responsible for more than 600,000 deaths each year. It occurs as two serotypes (A and D) representing two varieties (i.e. grubii and neoformans, respectively). Here, we sequenced the genome and performed an RNA-Seq-based analysis of the C. neoformans var. grubii transcriptome structure. We determined the chromosomal locations, analyzed the sequence/structural features of the centromeres, and identified origins of replication. The genome was annotated based on automated and manual curation. More than 40,000 introns populating more than 99% of the expressed genes were identified. Although most of these introns are located in the coding DNA sequences (CDS), over 2,000 introns in the untranslated regions (UTRs) were also identified. Poly(A)-containing reads were employed to locate the polyadenylation sites of more than 80% of the genes. Examination of the sequences around these sites revealed a new poly(A)-site-associated motif (AUGHAH). In addition, 1,197 miscRNAs were identified. These miscRNAs can be spliced and/or polyadenylated, but do not appear to have obvious coding capacities. Finally, this genome sequence enabled a comparative analysis of strain H99 variants obtained after laboratory passage. The spectrum of mutations identified provides insights into the genetics underlying the micro-evolution of a laboratory strain, and identifies mutations involved in stress responses, mating efficiency, and virulence.

Quantitative genetics of CTCF binding reveal local sequence effects and different modes of X-chromosome association.

Associating genetic variation with quantitative measures of gene regulation offers a way to bridge the gap between genotype and complex phenotypes. In order to identify quantitative trait loci (QTLs) that influence the binding of a transcription factor in humans, we measured binding of the multifunctional transcription and chromatin factor CTCF in 51 HapMap cell lines. We identified thousands of QTLs in which genotype differences were associated with differences in CTCF binding strength, hundreds of them confirmed by directly observable allele-specific binding bias. The majority of QTLs were either within 1 kb of the CTCF binding motif, or in linkage disequilibrium with a variant within 1 kb of the motif. On the X chromosome we observed three classes of binding sites: a minority class bound only to the active copy of the X chromosome, the majority class bound to both the active and inactive X, and a small set of female-specific CTCF sites associated with two non-coding RNA genes. In sum, our data reveal extensive genetic effects on CTCF binding, both direct and indirect, and identify a diversity of patterns of CTCF binding on the X chromosome.

Cell Fate Specification and the Regulation of RNA-dependent DNA Methylation in the Arabidopsis Root Meristem

The Arabidopsis root apical meristem (RAM) is a complex tissue capable of generating all the cell types that ultimately make up the root. The work presented in this thesis takes advantage of the versatility of high-throughput sequencing to address two independent questions about the root meristem. Although a lot of information is known regarding the cell fate decisions that occur at the RAM, cortex specification and differentiation remain poorly understood. In the first part of this thesis, I used an ethylmethanesulfonate (EMS) mutagenized marker line to perform a forward genetics screen. The goal of this screen was to identify novel genes involved in the specification and differentiation of the cortex tissue. Mapping analysis from the results obtained in this screen revealed a new allele of BRASSINOSTEROID4 with abnormal marker expression in the cortex tissue. Although this allele proved to be non-cortex specific, this project highlights new technology that allows mapping of EMS-generated mutations without the need to map-cross or back-cross. In the second part of this thesis, using fluorescence activated cell sorting (FACS) coupled with high throughput sequencing, my collaborators and I generated single-base resolution whole genome DNA methylomes, mRNA transcriptomes, and smallRNA transcriptomes for six different populations of cell types in the Arabidopsis root meristem. We were able to discover that the columella is hypermethylated in the CHH context within transposable elements. This hypermethylation is accompanied by upregulation of the RNA-dependent DNA methylation pathway (RdDM), including higher levels of 24-nt silencing RNAs (siRNAs). In summary, our studies demonstrate the versatility of high-throughput sequencing as a method for identifying single mutations or to perform complex comparative genomic analyses.

permGPU: Using graphics processing units in RNA microarray association studies.

BACKGROUND: Many analyses of microarray association studies involve permutation, bootstrap resampling and cross-validation, that are ideally formulated as embarrassingly parallel computing problems. Given that these analyses are computationally intensive, scalable approaches that can take advantage of multi-core processor systems need to be developed. RESULTS: We have developed a CUDA based implementation, permGPU, that employs graphics processing units in microarray association studies. We illustrate the performance and applicability of permGPU within the context of permutation resampling for a number of test statistics. An extensive simulation study demonstrates a dramatic increase in performance when using permGPU on an NVIDIA GTX 280 card compared to an optimized C/C++ solution running on a conventional Linux server. CONCLUSIONS: permGPU is available as an open-source stand-alone application and as an extension package for the R statistical environment. It provides a dramatic increase in performance for permutation resampling analysis in the context of microarray association studies. The current version offers six test statistics for carrying out permutation resampling analyses for binary, quantitative and censored time-to-event traits.

Co-orientation of replication and transcription preserves genome integrity.

In many bacteria, there is a genome-wide bias towards co-orientation of replication and transcription, with essential and/or highly-expressed genes further enriched co-directionally. We previously found that reversing this bias in the bacterium Bacillus subtilis slows replication elongation, and we proposed that this effect contributes to the evolutionary pressure selecting the transcription-replication co-orientation bias. This selection might have been based purely on selection for speedy replication; alternatively, the slowed replication might actually represent an average of individual replication-disruption events, each of which is counter-selected independently because genome integrity is selected. To differentiate these possibilities and define the precise forces driving this aspect of genome organization, we generated new strains with inversions either over approximately 1/4 of the chromosome or at ribosomal RNA (rRNA) operons. Applying mathematical analysis to genomic microarray snapshots, we found that replication rates vary dramatically within the inverted genome. Replication is moderately impeded throughout the inverted region, which results in a small but significant competitive disadvantage in minimal medium. Importantly, replication is strongly obstructed at inverted rRNA loci in rich medium. This obstruction results in disruption of DNA replication, activation of DNA damage responses, loss of genome integrity, and cell death. Our results strongly suggest that preservation of genome integrity drives the evolution of co-orientation of replication and transcription, a conserved feature of genome organization.