Search

The unstable (CTG·CAG)n trinucleotide repeat in the myotonic dystrophy type 1 (DM1) locus is bidirectionally transcribed from genes with terminal overlap. By transcription in the sense direction, the DMPK gene produces various alternatively spliced mRNAs with a (CUG)n repeat in their 3′ UTR. Expression in opposite orientation reportedly yields (CAG)n-repeat containing RNA, but both structure and biologic significance of this antisense gene (DM1-AS) are largely unknown. Via a combinatorial approach of computational and experimental analyses of RNA from unaffected individuals and DM1 patients we discovered that DM1-AS spans >6 kb, contains alternative transcription start sites and uses alternative polyadenylation sites up- and downstream of the (CAG)n repeat. Moreover, its primary transcripts undergo alternative splicing, whereby the (CAG)n segment is removed as part of an intron. Thus, in patients a mixture of DM1-AS RNAs with and without expanded (CAG)n repeat are produced. DM1-AS expression appears upregulated in patients, but transcript abundance remains very low in all tissues analyzed. Our data suggest that DM1-AS transcripts belong to the class of long non-coding RNAs. These and other biologically relevant implications for how (CAG)n-expanded transcripts may contribute to DM1 pathology can now be explored experimentally.

The scope and roles of regulated isoform gene expression during erythroid terminal development are poorly understood. We identified hundreds of differentiation-associated isoform changes during terminal erythropoiesis. Sequences surrounding cassette exons of skipped exon events are enriched for motifs bound by the Muscleblind-like (MBNL) family of splicing factors. Knockdown of Mbnl1 in cultured murine fetal liver erythroid progenitors resulted in a strong block in erythroid differentiation and disrupted the developmentally regulated exon skipping of Ndel1 mRNA, which is bound by MBNL1 and critical for erythroid terminal proliferation. These findings reveal an unanticipated scope of the alternative splicing program and the importance of Mbnl1 during erythroid terminal differentiation.

Objective To assess the performance of a directed chromosomal analysis approach in the prenatal evaluation of fetal sex chromosome aneuploidy. Methods We analyzed 432 frozen maternal plasma samples obtained from patients prior to undergoing fetal diagnostic testing. The cohort included women greater than 18 years of age with a singleton pregnancy of greater than 10 weeks gestation. Samples were analyzed using a chromosome-selective approach (DANSRTM) and a risk algorithm that incorporates fetal fraction (FORTETM). Results The cohort included 34 cases of sex chromosome aneuploidy. The assay correctly identified 26 of 27 (92.6%) cases of Monosomy X, one case of XXX, and all six cases of XXY. There were four false positive cases of sex chromosome aneuploidy among 380 euploid cases for an overall false positive rate of less than 1%. Discussion Analysis of the risk for sex chromosome aneuploidies can be accomplished with a targeted assay with high sensitivity. © 2014 John Wiley & Sons, Ltd.

Objective: To develop a microarray-based method for noninvasive prenatal testing (NIPT) and compare it with next-generation sequencing. Methods: Maternal plasma from 878 pregnant women, including 187 trisomy cases (18 trisomy 13, 37 trisomy 18, 132 trisomy 21), was evaluated for trisomy risk. Targeted chromosomes were analyzed using Digital Analysis of Selected Regions (DANSR™) assays. DANSR products were subsequently divided between two DNA quantification methods: microarrays and next-generation sequencing. For both microarray and sequencing methodologies, the Fetal-Fraction Optimized Risk of Trisomy Evaluation (FORTE™) algorithm was used to determine trisomy risk, assay variability across samples, and compute fetal fraction variability within samples. Results: NIPT using microarrays provided faster and more accurate cell-free DNA (cfDNA) measurements than sequencing. The assay variability, a measure of variance of chromosomal cfDNA counts, was lower for microarrays than for sequencing, 0.051 versus 0.099 (p < 0.0001). Analysis time using microarrays was faster, 7.5 versus 56 h for sequencing. Additionally, fetal fraction precision was improved 1.6-fold by assaying more polymorphic sites with microarrays (p < 0.0001). Microarrays correctly classified all trisomy and nontrisomy cases. Conclusions: NIPT using microarrays delivers more accurate cfDNA analysis than next-generation sequencing and can be performed in less time.

Myotonic dystrophy Type 1 (DM1) is a rare genetic disease caused by the expansion of CTG trinucleotide repeats ((CTG)exp) in the 3' untranslated region of the DMPK gene. The repeat transcripts sequester the RNA binding protein Muscleblind-like protein 1 (MBNL1) and hamper its normal function in pre-mRNA splicing. Overexpressing exogenous MBNL1 in the DM1 mouse model has been shown to rescue the splicing defects and reverse myotonia. Although a viable therapeutic strategy, pharmacological modulators of MBNL1 expression have not been identified. Here, we engineered a ZsGreen tag into the endogenous MBNL1 locus in HeLa cells and established a flow cytometry-based screening system to identify compounds that increase MBNL1 level. The initial screen of small molecule compound libraries identified more than thirty hits that increased MBNL1 expression greater than double the baseline levels. Further characterization of two hits revealed that the small molecule HDAC inhibitors, ISOX and vorinostat, increased MBNL1 expression in DM1 patient-derived fibroblasts and partially rescued the splicing defect caused by (CUG)exp repeats in these cells. These findings demonstrate the feasibility of this flow-based cytometry screen to identify both small molecule compounds and druggable targets for MBNL1 upregulation.

Skeletal muscle is a highly specialized, postmitotic tissue that must withstand chronic mechanical and physiological stress throughout life to maintain proper contractile function. Muscle damage or disease leads to progressive weakness and disability, and manifests in more than 100 different human disorders. Current therapies to treat muscle degenerative diseases are limited mostly to the amelioration of symptoms, although promising new therapeutic directions are emerging. In this review, we discuss the pathological basis for the most common muscle degenerative diseases and highlight new and encouraging experimental and clinical opportunities to prevent or reverse these afflictions.

Myotonic dystrophy (DM) is caused by the expression of mutant RNAs containing expanded CUG repeats that sequester muscleblind-like (MBNL) proteins, leading to alternative splicing changes. Cardiac alterations, characterized by conduction delays and arrhythmia, are the second most common cause of death in DM. Using RNA sequencing, here we identify novel splicing alterations in DM heart samples, including a switch from adult exon 6B towards fetal exon 6A in the cardiac sodium channel, SCN5A. We find that MBNL1 regulates alternative splicing of SCN5A mRNA and that the splicing variant of SCN5A produced in DM presents a reduced excitability compared with the control adult isoform. Importantly, reproducing splicing alteration of Scn5a in mice is sufficient to promote heart arrhythmia and cardiac-conduction delay, two predominant features of myotonic dystrophy. In conclusion, misregulation of the alternative splicing of SCN5A may contribute to a subset of the cardiac dysfunctions observed in myotonic dystrophy., Patients with myotonic dystrophy (MD) suffer from severe cardiac issues of unknown aetiology. Freyermuth et al. show that fatal changes in cardiac electrophysiological properties in humans and mice with MD may arise from misregulation of the alternative splicing of the cardiac Na+ channel SCN5A transcript, resulting in expression of its fetal form.

The RNA binding protein and alternative splicing factor Muscleblind-like 1 (MBNL1) has been a topic of intense study due to its role in myotonic dystrophy (DM) pathogenesis. MBNL1 contains four zinc finger (ZF) RNA binding domains arranged in two pairs. Through combinatorial mutagenesis of the ZF domains, we demonstrate that the pairs of ZFs have differential affinity for RNA and subsequently differential splicing activities. We evaluated splicing and binding activity for six MBNL1-mediated splicing events and found that the splicing activity profiles for the ZF mutants vary among transcripts. Clustering analysis of splicing profiles revealed that two distinct classes of MBNL1 pre-mRNA substrates exist. For some of the RNA transcripts tested, binding and splicing activity of the ZF mutants correlated. However, for some transcripts it appears that MBNL1 exerts robust splicing activity in the absence of RNA binding. We demonstrate that functionally distinct classes of MBNL1-mediated splicing events exist as defined by requirements for ZF-RNA interactions.

Myotonic dystrophy type 1 (DM1) is a triplet repeating disorder caused by expanded CTG repeats in the 3'-untranslated region of the dystrophia myotonica protein kinase (DMPK) gene. The transcribed repeats fold into an RNA hairpin with multiple copies of a 5'CUG/3'GUC motif that binds the RNA splicing regulator muscleblind-like 1 protein (MBNL1). Sequestration of MBNL1 by expanded r(CUG) repeats causes splicing defects in a subset of pre-mRNAs including the insulin receptor, the muscle-specific chloride ion channel, sarco(endo)plasmic reticulum Ca(2+) ATPase 1, and cardiac troponin T. Based on these observations, the development of small-molecule ligands that target specifically expanded DM1 repeats could be of use as therapeutics. In the present study, chemical similarity searching was employed to improve the efficacy of pentamidine and Hoechst 33258 ligands that have been shown previously to target the DM1 triplet repeat. A series of in vitro inhibitors of the RNA-protein complex were identified with low micromolar IC(50)'s, which are20-fold more potent than the query compounds. Importantly, a bis-benzimidazole identified from the Hoechst query improves DM1-associated pre-mRNA splicing defects in cell and mouse models of DM1 (when dosed with 1 mM and 100 mg/kg, respectively). Since Hoechst 33258 was identified as a DM1 binder through analysis of an RNA motif-ligand database, these studies suggest that lead ligands targeting RNA with improved biological activity can be identified by using a synergistic approach that combines analysis of known RNA-ligand interactions with chemical similarity searching.

Objective To develop a novel prenatal assay based on selective analysis of cell-free DNA in maternal blood for evaluation of fetal Trisomy 21 (T21) and Trisomy 18 (T18). Methods Two hundred ninety-eight pregnancies, including 39 T21 and seven T18 confirmed fetal aneuploidies, were analyzed using a novel, highly multiplexed assay, termed digital analysis of selected regions (DANSR ™ ). Cell-free DNA from maternal blood samples was analyzed using DANSR assays for loci on chromosomes 21 and 18. Products from 96 separate patients were pooled and sequenced together. A standard Z-test of chromosomal proportions was used to distinguish aneuploid samples from average-risk pregnancy samples. DANSR aneuploidy discrimination was evaluated at various sequence depths. Results At the lowest sequencing depth, corresponding to 204000 sequencing counts per sample, average-risk cases where distinguished from T21 and T18 cases, with Z statistics for all cases exceeding 3.6. Increasing the sequencing depth to 410000 counts per sample substantially improved separation of aneuploid and average-risk cases. A further increase to 620000 counts per sample resulted in only marginal improvement. This depth of sequencing represents less than 5% of that required by massively parallel shotgun sequencing approaches. Conclusion Digital analysis of selected regions enables highly accurate, cost efficient, and scalable noninvasive fetal aneuploidy assessment. © 2012 John Wiley & Sons, Ltd.

Bruton's Tyrosine Kinase (BTK), a TEC-family non-receptor tyrosine kinase, plays a critical role in B-cell development and function. Targeting of BTK with covalent inhibitors like ibrutinib has become a standard approach for many B-cell malignancies, including chronic lymphocytic leukemia (CLL), Waldenstrom macroglobulinemia, mantle cell lymphoma (MCL) and marginal zone lymphoma (MZL). Yet, many patients demonstrate intrinsic or acquired resistance to covalent BTK inhibition. The prognosis of patients who relapse after ibrutinib treatment for MCL or CLL is dismal, highlighting the urgent need for new approaches that overcome resistance to current BTK inhibitors. We hypothesized that degradation of BTK could be a better alternative to inhibition alone, as it would both: 1) maintain efficacy in cells harboring the ibrutinib-resistant BTK C418S mutation and 2) target non-catalytic functions of BTK. To address this, we turned to a small molecule-mediated protein degradation platform that utilizes an E3 ligase-targeting moiety linked to the ligand of a target of interest, so that the target can be marked for ubiquitination and subsequent proteasomal degradation. We previously showed that BTK is one of the most robustly degraded kinase targets using a nonspecific kinase inhibitor linked to a imide-based core followed by agnostic proteomics. We synthesized highly potent and selective degraders of BTK using imides as a base. To do so, the parent BTK inhibitor CGI1746 was linked to thalidomide using either polyethylene glycol (DD-03-007) linkers or saturated hydrocarbon chain (DD-03-171) linkers. After verifying that these degraders induce dimerization of BTK and CRBN and penetrate cells, we explored the pharmacological effects in vitro. Both DD-03-007 and DD-03-171 reduced BTK levels at concentrations as low as 40 nM within 4h of treatment. Furthermore, cellular BTK levels remained low for 24h after washout, showing that these degraders are capable of sustaining depletion of BTK for an extended period of time (Figure). DD-03-007 and DD-03-171 are both ligase and proteasome-dependent, as shown by co-treatment with bortezomib or MLN-4924 as well as competition experiments with lenalidomide or CGI-1746. Moreover, the degraders exhibited strong synergy with the HCK inhibitor A419259, suggesting that it would be possible to recapitulate ibrutinib's previously reported polypharmacology with an HCK inhibitor. We overexpressed wild-type or C481S BTK in TMD-8 cells and ran an antiproliferation assay, which showed that DD-03-007 overcame ibrutinib resistance associated with BTK C481S mutation. Next, we explored the effects of the degraders in MCL specifically. Proteomic analysis showed that DD-03-171 is a triple-degrader, as it degrades BTK but retains degradation activity on IKZF1 and IKZF3. Finally, we performed In vivo efficacy studies in cell line and patient-derived xenograft (PDX) models of MCL. The latter was obtained from a patient who had progressed on ibrutinib. In both models, DD-03-171 caused a significant reduction in tumor burden at an early timepoint. DD-03-171 also markedly extended survival compared to treatment with ibrutinib or lenalidomide alone (Figure). In conclusion, we developed highly potent and selective BTK degraders with activity in vivo against human MCL that induce the degradation of multiple factors essential for MCL survival. Figure. Figure. Disclosures Treon: Pharmacyclics: Consultancy, Other: Travel, Accommodations, Expenses, Research Funding; Janssen: Consultancy, Other: Travel, Accommodations, Expenses; Johnson & Johnson: Consultancy; BMS: Research Funding. Weinstock:Novartis, Dragonfly, Travera, DxTerity, Travera: Consultancy; Novartis: Consultancy, Research Funding; Novartis, Astra Zeneca, Abbvie, Aileron, Surface Oncology, Daiichi Sankyo: Research Funding; Genentech/Roche, Monsanto: Consultancy; Astra Zeneca, JAX, Samumed, Regeneron, Sun Pharma, Prescient: Patents & Royalties; Travera: Equity Ownership. Gray:Syros, Soltego, Petra, C4 Therapeutics: Equity Ownership.

RNA-binding proteins (RBPs) regulate many aspects of transcription and translation in a cell- and tissue-specific manner and are frequently dysregulated in malignancy. We systematically evaluated RBPs preferentially required in acute myeloid leukemia (AML) over other forms of cancer or normal hematopoietic precursors using a CRISPR/Cas9 domain-based, loss-of-function screen targeting 490 classical RBPs with 2,900 sgRNAs (Fig. A). This screen was performed in cells lines representing AML, T-cell acute lymphoblastic leukemia (T-ALL), and lung adenocarcinoma (LUAD) and revealed multiple RBPs preferentially required for AML survival, but not for T-ALL or LUAD survival. We identified genes encoding 21 RBPs that were >3-fold depleted in AML cells and significantly overexpressed in AML patient samples versus normal adult CD34+ precursors (p-value < 0.05; Fig. B). Amongst RBPs required and upregulated in AML was RBM39, an RBP described to be involved in a number of cellular processes and to interact with key splicing proteins SF3B1 and U2AF2. Genetic ablation of Rbm39 in mouse MLL-AF9 leukemia cells dramatically delayed AML development and progression (Fig. C). In parallel, it has recently been described that a class of clinically-validated anti-cancer sulfonamide compounds (including indisulam and E7820) mediate RBM39 degradation as their dominant cellular mechanism of action. This occurs via novel interactions with the DCAF15 adapter protein of the CUL4/Ddb1 ubiquitin ligase complex with RBM39 as a neo-substrate. Treatment of MOLM-13 cells xenografted into mice with indisulam conferred significant anti-leukemic effects and improved overall survival (Fig. D). To explore the mechanism of RBM39 dependence in AML, we performed proteomic analyses of RBM39 interacting proteins in MOLM-13 cells as well as transcriptome-wide analysis of RBM39 RNA binding by enhanced UV cross-linking and immunoprecipitation (eCLIP) in the same cells. RBM39 physically interacted with an entire network of RBPs identified by our CRISPR screen as crucial for AML cell survival in addition to interacting with the core SF3b splicing complex. Further, anti-RBM39 eCLIP revealed RBM39 binding to exonic regions and most enriched at exon/intron borders at 5' and 3' splice sites of pre-mRNA (Fig. E), suggesting a prominent role of RBM39 in regulating splicing. Consistent with this, RNA-sequencing of AML cells following RBM39 deletion revealed significant effects of RBM39 loss on RNA splicing, most prominently causing increased cassette exon skipping (Fig. F). Recent studies suggest that myeloid leukemias with mutations in RNA splicing factors are sensitized to pharmacologic perturbation of RNA splicing. Analysis of the effects of RBM39 degrading compounds over a panel of 18 AML cells revealed that leukemia cells bearing splicing factor mutations or with high DCAF15 expression were the most sensitive to treatment (Fig. G). Genetic introduction of SF3B1, SRSF2, or U2AF1 hotspot mutations in K562 or NALM6 cells resulted in a 20-50% reduction in IC50 in response to sulfonamides. We next performed RNA sequencing of isogenic K562 cells with or without knockin of SF3B1K700E and SRSF2P95H mutations into the endogenous loci, and treated at the IC50 of E7820 or E7107, a small molecule that inhibits the SF3b core spliceosome complex. Treatment with either drug dramatically increased cassette exon skipping events and intron retention relative to DMSO control, with greater effects in splicing mutant cells. However, at equipotent doses, E7820 markedly increased mis-splicing compared with E7107. Furthermore, E7820 treatment resulted in mis-splicing of a number of RBP targets identified in our CRISPR screen as being required for AML survival, including SUPT6H, hnRNPH, and SRSF10, as well as RBM3 and U2AF2, consistent with previous observations (Fig. H). Here through systematic evaluation of RBPs across several cancers, we identify RBPs specifically required in AML. In so doing we identify a network of functionally and physically interacting RBPs upregulated in AML over normal precursors. Genetic or pharmacologic elimination one such RBP, RBM39, led to aberrant splicing of multiple members of this RBP network as well as of transcriptional regulators required for AML survival. These data suggest important clinical potential for anti-cancer sulfonamide treatment in splicing mutant myeloid leukemias. Disclosures Uehara: Eisai: Employment. Owa:Eisai: Employment.

We set out to develop a molecular test that distinguishes benign and malignant thyroid nodules using fine-needle aspirates (FNA).We used mRNA expression analysis to measure more than 247,186 transcripts in 315 thyroid nodules, comprising multiple subtypes. The data set consisted of 178 retrospective surgical tissues and 137 prospectively collected FNA samples. Two classifiers were trained separately on surgical tissues and FNAs. The performance was evaluated using an independent set of 48 prospective FNA samples, which included 50% with indeterminate cytopathology.Performance of the tissue-trained classifier was markedly lower in FNAs than in tissue. Exploratory analysis pointed to differences in cellular heterogeneity between tissues and FNAs as the likely cause. The classifier trained on FNA samples resulted in increased performance, estimated using both 30-fold cross-validation and an independent test set. On the test set, negative predictive value and specificity were estimated to be 96 and 84%, respectively, suggesting clinical utility in the management of patients considering surgery. Using in silico and in vitro mixing experiments, we demonstrated that even in the presence of 80% dilution with benign background, the classifier can correctly recognize malignancy in the majority of FNA samples.The FNA-trained classifier was able to classify an independent set of FNAs in which substantial RNA degradation had occurred and in the presence of blood. High tolerance to dilution makes the classifier useful in routine clinical settings where sampling error may be a concern. An ongoing multicenter clinical trial will allow us to validate molecular test performance on a larger independent test set of prospectively collected thyroid FNAs.

Pre-mRNA splicing is regulated through the combinatorial activity of RNA motifs, including splice sites and splicing regulatory elements. Here we show that the activity of the G-run (polyguanine sequence) class of splicing enhancer elements is approx4-fold higher when adjacent to intermediate strength 5' splice sites (ss) than when adjacent to weak 5' ss, and approx1.3-fold higher relative to strong 5' ss. We observed this dependence on 5' ss strength in both splicing reporters and in global microarray and mRNA-Seq analyses of splicing changes following RNA interference against heterogeneous nuclear ribonucleoprotein (hnRNP) H, which cross-linked to G-runs adjacent to many regulated exons. An exon's responsiveness to changes in hnRNP H levels therefore depends in a complex way on G-run abundance and 5' ss strength. This pattern of activity enables G-runs and hnRNP H to buffer the effects of 5' ss mutations, augmenting both the frequency of 5' ss polymorphism and the evolution of new splicing patterns. Certain other splicing factors may function similarly., American Heart Association, Human Frontier Science Program (Strasbourg, France), National Institutes of Health (U.S.), National Science Foundation (U.S.) (equipment grant DBI-0821391)

Genomic surveys in humans identify a large amount of recent positive selection. Using the 3.9-million HapMap SNP dataset, we found that selection has accelerated greatly during the last 40,000 years. We tested the null hypothesis that the observed age distribution of recent positively selected linkage blocks is consistent with a constant rate of adaptive substitution during human evolution. We show that a constant rate high enough to explain the number of recently selected variants would predict ( i ) site heterozygosity at least 10-fold lower than is observed in humans, ( ii ) a strong relationship of heterozygosity and local recombination rate, which is not observed in humans, ( iii ) an implausibly high number of adaptive substitutions between humans and chimpanzees, and ( iv ) nearly 100 times the observed number of high-frequency linkage disequilibrium blocks. Larger populations generate more new selected mutations, and we show the consistency of the observed data with the historical pattern of human population growth. We consider human demographic growth to be linked with past changes in human cultures and ecologies. Both processes have contributed to the extraordinarily rapid recent genetic evolution of our species.

Local patterns of biomechanical forces experienced by endothelial cells (ECs) in different vascular geometries appear to play an essential role in regulating EC function and determining the regional susceptibility to atherosclerosis, even in the face of systemic risk factors. To study how biomechanical forces regulate EC redox homeostasis, an important pathogenic factor in atherogenesis, we have cultured human ECs under 2 prototypic arterial shear stress waveforms, "atheroprone" and "atheroprotective," which were derived from 2 distinct vascular regions in vivo that are typically "susceptible" or "resistant" to atherosclerosis. We demonstrate that atheroprotective flow decreases EC intracellular redox level and protects ECs against oxidative stress-induced injury. To identify the molecular mechanisms that control this cellular response, we examined several major oxidative/antioxidative pathways and found that atheroprotective flow upregulated certain antioxidant genes and strongly activated the transcription factor Nrf2. Using a strategy of small interfering RNA inhibition of Nrf2 expression combined with genome-wide transcriptional profiling, we determined the downstream targets of Nrf2 activation and identified Nrf2 as a critical determinant for the changes in endothelial redox balance exerted by atheroprotective flow. In addition, we showed that atheroprotective flow activates Nrf2 via the phosphoinositol 3-kinase/Akt pathway, and this activation occurs differentially in atherosclero- sis-resistant and atherosclerosis-susceptible regions of the mouse aorta. Taken together, our data demonstrate that hemodynamic forces present in atherosclerosis-resistant and -susceptible regions of the vasculature differentially regulate EC redox state and antioxidant potential. These alterations in redox homeostasis are primarily the result of the phosphoinositol 3-kinase/Akt-dependent activation of Nrf2 and its downstream transcriptional targets. (Circ Res. 2007;101:723-733.)

Spatial restriction of mRNA to distinct subcellular locations enables local regulation and synthesis of proteins. However, the organizing principles of mRNA localization remain poorly understood. Here we analyzed subcellular transcriptomes of neural projections and soma of primary mouse cortical neurons and two neuronal cell lines and found that alternative last exons (ALEs) often confer isoform-specific localization. Surprisingly, gene-distal ALE isoforms were four times more often localized to neurites than gene-proximal isoforms. Localized isoforms were induced during neuronal differentiation and enriched for motifs associated with muscleblind-like (Mbnl) family RNA-binding proteins. Depletion of Mbnl1 and/or Mbnl2 reduced localization of hundreds of transcripts, implicating Mbnls in localization of mRNAs to neurites. We provide evidence supporting a model in which the linkage between genomic position of ALEs and subcellular localization enables coordinated induction of localization-competent mRNA isoforms through a post-transcriptional regulatory program that is induced during differentiation and reversed in cellular reprogramming and cancer. Taliaferro et al. show that mRNA isoforms that contain gene-distal alternative last exons are preferentially localized to neurites. These isoforms are induced during neuronal differentiation, suggesting that a coordinated post-transcriptional program targets messages to neurites. Localization of many neuronal mRNAs depends on muscleblind proteins.

Alternative splicing is a regulated process that results in expression of specific mRNA and protein isoforms. Alternative splicing factors determine the relative abundance of each isoform. Here we focus on MBNL1, a splicing factor misregulated in the disease myotonic dystrophy. By altering the concentration of MBNL1 in cells across a broad dynamic range, we show that different splicing events require different amounts of MBNL1 for half-maximal response, and respond more or less steeply to MBNL1. Motifs around MBNL1 exon 5 were studied to assess how cis-elements mediate the MBNL1 dose-dependent splicing response. A framework was developed to estimate MBNL concentration using splicing responses alone, validated in the cell-based model, and applied to myotonic dystrophy patient muscle. Using this framework, we evaluated the ability of individual and combinations of splicing events to predict functional MBNL concentration in human biopsies, as well as their performance as biomarkers to assay mild, moderate, and severe cases of DM., Author Summary Our studies provide insight into the mechanisms of myotonic dystrophy, the most common adult form of muscular dystrophy. In this disease, a family of RNA binding proteins is sequestered by toxic RNA, which leads to mis-regulation and disease symptoms. We have created a cellular model with one of these family members to study how these RNA binding proteins function in the absence of the toxic RNA. In parallel, we analyzed transcriptomic data from over 50 individuals (44 affected by myotonic dystrophy) with a range of disease severity. The results from the transcriptomic data provide a rational approach to select biomarkers for clinical research and therapeutic trials.

RNA binding proteins of the conserved CUGBP1, Elav-like factor (CELF) family contribute to heart and skeletal muscle development and are implicated in myotonic dystrophy (DM). To understand their genome-wide functions, we analyzed the transcriptome dynamics following induction of CELF1 or CELF2 in adult mouse heart and of CELF1 in muscle by RNA-seq, complemented by crosslinking/immunoprecipitation-sequencing (CLIP-seq) analysis of mouse cells and tissues to distinguish direct from indirect regulatory targets. We identified hundreds of mRNAs bound in their 3′ UTRs by both CELF1 and the developmentally induced MBNL1 protein, a threefold greater overlap in target messages than expected, including messages involved in development and cell differentiation. The extent of 3′ UTR binding by CELF1 and MBNL1 predicted the degree of mRNA repression or stabilization, respectively, following CELF1 induction. However, CELF1's RNA binding specificity in vitro was not detectably altered by coincubation with recombinant MBNL1. These findings support a model in which CELF and MBNL proteins bind independently to mRNAs but functionally compete to specify down-regulation or localization/stabilization, respectively, of hundreds of mRNA targets. Expression of many alternative 3′ UTR isoforms was altered following CELF1 induction, with 3′ UTR binding associated with down-regulation of isoforms and genes. The splicing of hundreds of alternative exons was oppositely regulated by these proteins, confirming an additional layer of regulatory antagonism previously observed in a handful of cases. The regulatory relationships between CELFs and MBNLs in control of both mRNA abundance and splicing appear to have evolved to enhance developmental transitions in major classes of heart and muscle genes., Myotonic Dystrophy Foundation, National Institutes of Health (U.S.) (Grant OD017865-02), National Science Foundation (U.S.) (Award 0821391), National Institutes of Health (U.S.) (Grant R01GM085319), National Institutes of Health (U.S.) (Grant RC2HG005624)

In the face of systemic risk factors, certain regions of the arterial vasculature remain relatively resistant to the development of atherosclerotic lesions. The biomechanically distinct environments in these arterial geometries exert a protective influence via certain key functions of the endothelial lining; however, the mechanisms underlying the coordinated regulation of specific mechano-activated transcriptional programs leading to distinct endothelial functional phenotypes have remained elusive. Here, we show that the transcription factor Kruppel-like factor 2 (KLF2) is selectively induced in endothelial cells exposed to a biomechanical stimulus characteristic of atheroprotected regions of the human carotid and that this flow-mediated increase in expression occurs via a MEK5/ERK5/MEF2 signaling pathway. Overexpression and silencing of KLF2 in the context of flow, combined with findings from genome-wide analyses of gene expression, demonstrate that the induction of KLF2 results in the orchestrated regulation of endothelial transcriptional programs controlling inflammation, thrombosis/hemostasis, vascular tone, and blood vessel development. Our data also indicate that KLF2 expression globally modulates IL-1beta-mediated endothelial activation. KLF2 therefore serves as a mechano-activated transcription factor important in the integration of multiple endothelial functions associated with regions of the arterial vasculature that are relatively resistant to atherogenesis.

TBX5 is a member of the T-box gene family and encodes a transcription factor involved in cardiac and limb development. Mutations of TBX5 cause Holt-Oram syndrome (HOS), an autosomal-dominant condition with congenital cardiac defects and forelimb anomalies. Here, we used a GAL4-TBX5 fusion protein in a modified yeast-one hybrid system to elucidate the TBX5 transactivating domain. Using a series of deletion mutations of TBX5, we narrowed down its functional domain to amino acids 339-379 of its C-terminal half; point mutagenesis analysis then showed that the loss of amino acids 349-351 abolished transactivation. This result was confirmed in mammalian cells. Furthermore, wild-type TBX5, but not TBX5 with mutations at the amino acids 349-351, has ability to inhibit NCI-H1299 cell growth also suggesting that these amino acids are crucial for the TBX5 function in mammalian cells. In addition, to identify the nuclear localization signal of TBX5, we searched for cluster of basic amino acids. We found that the deletion of the KRK sequence at amino acids 325-327 mislocalizes TBX5 to cytoplasm, suggesting that these amino acids serve as a nuclear localization signal. These studies enhance our understanding of the structure-function relationship of TBX5 and suggest that truncation mutations of TBX5 could cause HOS through the loss of its transactivating domain and/or the nuclear localization signal.

Motivation: Analysis of RNA sequencing (RNA-Seq) data revealed that the vast majority of human genes express multiple mRNA isoforms, produced by alternative pre-mRNA splicing and other mechanisms, and that most alternative isoforms vary in expression between human tissues. As RNA-Seq datasets grow in size, it remains challenging to visualize isoform expression across multiple samples. Results: To help address this problem, we present Sashimi plots, a quantitative visualization of aligned RNA-Seq reads that enables quantitative comparison of exon usage across samples or experimental conditions. Sashimi plots can be made using the Broad Integrated Genome Viewer or with a stand-alone command line program. Availability and implementation: Software code and documentation freely available here: http://miso.readthedocs.org/en/fastmiso/sashimi.html Contact: mesirov@broadinstitute.org, airoldi@fas.harvard.edu or cburge@mit.edu Supplementary information: Supplementary data are available at Bioinformatics online.

During postnatal development the heart undergoes a rapid and dramatic transition to adult function through transcriptional and post-transcriptional mechanisms, including alternative splicing (AS). Here we perform deep RNA-sequencing on RNA from cardiomyocytes and cardiac fibroblasts to conduct a high-resolution analysis of transcriptome changes during postnatal mouse heart development. We reveal extensive changes in gene expression and AS that occur primarily between postnatal days 1 and 28. Cardiomyocytes and cardiac fibroblasts show reciprocal regulation of gene expression reflecting differences in proliferative capacity, cell adhesion functions, and mitochondrial metabolism. We further demonstrate that AS plays a role in vesicular trafficking and membrane organization, These AS transitions are enriched among targets of two RNA-binding proteins, Celf1 and Mbnl1, which undergo developmentally regulated changes in expression. Vesicular trafficking genes affected by AS during normal development (when Celf1 is down-regulated) show a reversion to neonatal splicing patterns after Celf1 re-expression in adults. Short-term Celf1 induction in adult animals results in disrupted transverse tubule organization and calcium handling. These results identify potential roles for AS in multiple aspects of postnatal heart maturation, including vesicular trafficking and intracellular membrane dynamics.

Analysis of RNA sequencing (RNA-Seq) data revealed that the vast majority of human genes express multiple mRNA isoforms, produced by alternative pre-mRNA splicing and other mechanisms, and that most alternative isoforms vary in expression between human tissues. As RNA-Seq datasets grow in size, it remains challenging to visualize isoform expression across multiple samples. We present Sashimi plots, a quantitative multi-sample visualization of RNA-Seq reads aligned to gene annotations, which enables quantitative comparison of isoform usage across samples or experimental conditions. Given an input annotation and spliced alignments of reads from a sample, a region of interest is visualized in a Sashimi plot as follows: (i) alignments in exons are represented as read densities (optionally normalized by length of genomic region and coverage), and (ii) splice junction reads are drawn as arcs connecting a pair of exons, where arc width is drawn proportional to the number of reads aligning to the junction.

This study identifies MBNL proteins as negative regulators of alternative splicing events that are differentially regulated between ES cells and other cell types; several lines of evidence show that these proteins repress an ES cell alternative splicing program and the reprogramming of somatic cells to induced pluripotent stem cells. Ben Blencowe and colleagues identify the muscleblind-like RNA binding proteins MBNL1 and MBNL2 as negative regulators of alternative splicing events that are differentially regulated between embryonic stem cells and other cell types. Several lines of evidence show that they are involved in the regulation of embryonic-stem-cell-like alternative splicing patterns. The authors also identify a regulatory role during the reprogramming of fibroblasts to induced pluripotent stem (iPS) cells. Previous investigations of the core gene regulatory circuitry that controls the pluripotency of embryonic stem (ES) cells have largely focused on the roles of transcription, chromatin and non-coding RNA regulators1,2,3. Alternative splicing represents a widely acting mode of gene regulation4,5,6,7,8, yet its role in regulating ES-cell pluripotency and differentiation is poorly understood. Here we identify the muscleblind-like RNA binding proteins, MBNL1 and MBNL2, as conserved and direct negative regulators of a large program of cassette exon alternative splicing events that are differentially regulated between ES cells and other cell types. Knockdown of MBNL proteins in differentiated cells causes switching to an ES-cell-like alternative splicing pattern for approximately half of these events, whereas overexpression of MBNL proteins in ES cells promotes differentiated-cell-like alternative splicing patterns. Among the MBNL-regulated events is an ES-cell-specific alternative splicing switch in the forkhead family transcription factor FOXP1 that controls pluripotency9. Consistent with a central and negative regulatory role for MBNL proteins in pluripotency, their knockdown significantly enhances the expression of key pluripotency genes and the formation of induced pluripotent stem cells during somatic cell reprogramming.

Longevity is influenced by genetic and environmental factors. The brain's dopamine system may be particularly relevant, since it modulates traits (e.g., sensitivity to reward, incentive motivation, sustained effort) that impact behavioral responses to the environment. In particular, the dopamine D4 receptor (DRD4) has been shown to moderate the impact of environments on behavior and health. We tested the hypothesis that theDRD4gene influences longevity and that its impact is mediated through environmental effects. Surviving participants of a 30-year-old population-based health survey (N= 310; age range, 90–109 years; the 90+ Study) were genotyped/resequenced at theDRD4gene and compared with a European ancestry-matched younger population (N= 2902; age range, 7–45 years). We found that the oldest-old population had a 66% increase in individuals carrying theDRD47R allele relative to the younger sample (p= 3.5 × 10−9), and that this genotype was strongly correlated with increased levels of physical activity. Consistent with these results,DRD4knock-out mice, when compared with wild-type and heterozygous mice, displayed a 7–9.7% decrease in lifespan, reduced spontaneous locomotor activity, and no lifespan increase when reared in an enriched environment. These results support the hypothesis thatDRD4gene variants contribute to longevity in humans and in mice, and suggest that this effect is mediated by shaping behavioral responses to the environment.

Objective: To determine the relationship between a priori risk for fetal trisomy and the fraction of fetal cell-free DNA (cfDNA) in maternal blood. Methods: A comparative analysis on fetal cfDNA amounts was performed in subjects stratified into a priori risk groups based on maternal age, prenatal screening results, or nuchal translucency measurement. Results: Across the highest and lowest deciles within each group, there were no significant differences in the fetal cfDNA fraction. Conclusions: These data support the concept that non-invasive prenatal test performance as determined by fetal cfDNA fraction is not predicted to be different based on patient risk classification.

We sought to develop a novel biochemical assay and algorithm for the prenatal evaluation of risk for fetal trisomy 21 (T21) and trisomy 18 (T18) using cell-free DNA obtained from maternal blood.We assayed cell-free DNA from a training set and a blinded validation set of pregnant women, comprising 250 disomy, 72 T21, and 16 T18 pregnancies. We used digital analysis of selected regions in combination with a novel algorithm, fetal-fraction optimized risk of trisomy evaluation (FORTE), to determine trisomy risk for each subject.In all, 163/171 subjects in the training set passed quality control criteria. Using a Z statistic, 35/35 T21 cases and 7/7 T18 cases had Z statistic3 and 120/121 disomic cases had Z statistic3. FORTE produced an individualized trisomy risk score for each subject, and correctly discriminated all T21 and T18 cases from disomic cases. All 167 subjects in the blinded validation set passed quality control and FORTE performance matched that observed in the training set correctly discriminating 36/36 T21 cases and 8/8 T18 cases from 123/123 disomic cases.Digital analysis of selected regions and FORTE enable accurate, scalable noninvasive fetal aneuploidy detection.

To develop a novel prenatal assay based on selective analysis of cell-free DNA in maternal blood for evaluation of fetal Trisomy 21 (T21) and Trisomy 18 (T18).Two hundred ninety-eight pregnancies, including 39 T21 and seven T18 confirmed fetal aneuploidies, were analyzed using a novel, highly multiplexed assay, termed digital analysis of selected regions (DANSR™). Cell-free DNA from maternal blood samples was analyzed using DANSR assays for loci on chromosomes 21 and 18. Products from 96 separate patients were pooled and sequenced together. A standard Z-test of chromosomal proportions was used to distinguish aneuploid samples from average-risk pregnancy samples. DANSR aneuploidy discrimination was evaluated at various sequence depths.At the lowest sequencing depth, corresponding to 204,000 sequencing counts per sample, average-risk cases where distinguished from T21 and T18 cases, with Z statistics for all cases exceeding 3.6. Increasing the sequencing depth to 410,000 counts per sample substantially improved separation of aneuploid and average-risk cases. A further increase to 620,000 counts per sample resulted in only marginal improvement. This depth of sequencing represents less than 5% of that required by massively parallel shotgun sequencing approaches.Digital analysis of selected regions enables highly accurate, cost efficient, and scalable noninvasive fetal aneuploidy assessment.

MicroRNAs (miRNAs) regulate physiological and pathological processes by inducing posttranscriptional repression of target messenger RNAs (mRNAs) via incompletely understood mechanisms. To discover factors required for human miRNA activity, we performed an RNAi screen using a reporter cell line of miRNA-mediated repression of translation initiation. We report that reduced expression of ribosomal protein genes (RPGs) dissociated miRNA complexes from target mRNAs, leading to increased polysome association, translation, and stability of miRNA-targeted mRNAs relative to untargeted mRNAs. RNA sequencing of polysomes indicated substantial overlap in sets of genes exhibiting increased or decreased polysomal association after Argonaute or RPG knockdowns, suggesting similarity in affected pathways. miRNA profiling of monosomes and polysomes demonstrated that miRNAs cosediment with ribosomes. RPG knockdowns decreased miRNAs in monosomes and increased their target mRNAs in polysomes. Our data show that most miRNAs repress translation and that the levels of RPGs modulate miRNA-mediated repression of translation initiation.

The immunologic consequences of prolonged infusions of rIL-2 in doses that produce physiologic serum concentrations of this cytokine were investigated. rIL-2 in doses of 0.5-6.0 x 10(6) U/m2 per d (3.3-40 micrograms/m2 per d) was administered by continuous intravenous infusion for 90 consecutive days to patients with advanced cancer. IL-2 concentrations (25 +/- 25 and 77 +/- 64 pM, respectively) that selectively saturate high-affinity IL-2 receptors (IL-2R) were achieved in the serum of patients receiving rIL-2 infusions of 10 micrograms/m2 per d and 30 micrograms/m2 per d. A gradual, progressive expansion of natural killer (NK) cells was seen in the peripheral blood of these patients with no evidence of a plateau effect during the 3 mo of therapy. A preferential expansion of CD56bright NK cells was consistently evident. NK cytotoxicity against tumor targets was only slightly enhanced at these dose levels. However, brief incubation of these expanded NK cells with IL-2 in vitro induced potent lysis of NK-sensitive, NK-resistant, and antibody-coated targets. Infusions of rIL-2 at 40 micrograms/m2 per d produced serum IL-2 levels (345 +/- 381 pM) sufficient to engage intermediate affinity IL-2R p75, which is constitutively expressed by human NK cells. This did not result in greater NK cell expansion compared to the lower dose levels, but did produce in vivo activation of NK cytotoxicity, as evidenced by lysis of NK-resistant targets. There was no consistent change in the numbers of CD56- CD3+ T cells, CD56+ CD3+ MHC-unrestricted T cells, or B cells during infusions of rIL-2 at any of the dosages used. This study demonstrates that prolonged infusions of rIL-2 in doses that saturate only high affinity IL-2R can selectively expand human NK cells for an extended period of time with only minimal toxicity. Further activation of NK cytolytic activity can also be achieved in vivo, but it requires concentrations of IL-2 that bind intermediate affinity IL-2R p75. Clinical trials are underway attempting to exploit the differing effects of various concentrations of IL-2 on human NK cells in vivo.

Through alternative splicing, most human genes express multiple isoforms that often differ in function. To infer isoform regulation from high-throughput sequencing of cDNA fragments (RNA-seq), we developed the mixture-of-isoforms (MISO) model, a statistical model that estimates expression of alternatively spliced exons and isoforms and assesses confidence in these estimates. Incorporation of mRNA fragment length distribution in paired-end RNA-seq greatly improved estimation of alternative-splicing levels. MISO also detects differentially regulated exons or isoforms. Application of MISO implicated the RNA splicing factor hnRNP H1 in the regulation of alternative cleavage and polyadenylation, a role that was supported by UV cross-linking-immunoprecipitation sequencing (CLIP-seq) analysis in human cells. Our results provide a probabilistic framework for RNA-seq analysis, give functional insights into pre-mRNA processing and yield guidelines for the optimal design of RNA-seq experiments for studies of gene and isoform expression.

Recent genome-wide analyses have elucidated the extent of alternative splicing (AS) in mammals, often focusing on comparisons of splice isoforms between differentiated tissues. However, regulated splicing changes are likely to be important in biological transitions such as cellular differentiation, or response to environmental stimuli. To assess the extent and significance of AS in myogenesis, we used splicing-sensitive microarray analysis of differentiating C2C12 myoblasts. We identified 95 AS events that undergo robust splicing transitions during C2C12 differentiation. More than half of the splicing transitions are conserved during differentiation of avian myoblasts, suggesting the products and timing of transitions are functionally significant. The majority of splicing transitions during C2C12 differentiation fall into four temporal patterns and were dependent on the myogenic program, suggesting that they are integral components of myogenic differentiation. Computational analyses revealed enrichment of many sequence motifs within the upstream and downstream intronic regions near the alternatively spliced regions corresponding to binding sites of splicing regulators. Western analyses demonstrated that several splicing regulators undergo dynamic changes in nuclear abundance during differentiation. These findings show that within a developmental context, AS is a highly regulated and conserved process, suggesting a major role for AS regulation in myogenic differentiation., National Institutes of Health (U.S.) (grant number R01GM076493), Ford Foundation (Predoctoral Diversity Fellowship), Baylor College of Medicine. Graduate School of Biomedical Sciences (Baylor Research Advocates for Student Scientists)

Australia was probably settled soon after modern humans left Africa, but details of this ancient migration are not well understood. Debate centers on whether the Pleistocene Sahul continent (composed of New Guinea, Australia, and Tasmania) was first settled by a single wave followed by regional divergence into Aboriginal Australian and New Guinean populations (common origin) or whether different parts of the continent were initially populated independently. Australia has been the subject of relatively few DNA studies even though understanding regional variation in genomic structure and diversity will be important if disease-association mapping methods are to be successfully evaluated and applied across populations. We report on a genome-wide investigation of Australian Aboriginal SNP diversity in a sample of participants from the Riverine region. The phylogenetic relationship of these Aboriginal Australians to a range of other global populations demonstrates a deep common origin with Papuan New Guineans and Melanesians, with little evidence of substantial later migration until the very recent arrival of European colonists. The study provides valuable and robust insights into an early and important phase of human colonization of the globe. A broader survey of Australia, including diverse geographic sample populations, will be required to fully appreciate the continent's unique population history and consequent genetic heritage, as well as the importance of both to the understanding of health issues.

The parts of the genome transcribed by a cell or tissue reflect the biological processes and functions it carries out. We characterized the features of mammalian tissue transcriptomes at the gene level through analysis of RNA deep sequencing (RNA-Seq) data across human and mouse tissues and cell lines. We observed that roughly 8,000 protein-coding genes were ubiquitously expressed, contributing to around 75% of all mRNAs by message copy number in most tissues. These mRNAs encoded proteins that were often intracellular, and tended to be involved in metabolism, transcription, RNA processing or translation. In contrast, genes for secreted or plasma membrane proteins were generally expressed in only a subset of tissues. The distribution of expression levels was broad but fairly continuous: no support was found for the concept of distinct expression classes of genes. Expression estimates that included reads mapping to coding exons only correlated better with qRT-PCR data than estimates which also included 3′ untranslated regions (UTRs). Muscle and liver had the least complex transcriptomes, in that they expressed predominantly ubiquitous genes and a large fraction of the transcripts came from a few highly expressed genes, whereas brain, kidney and testis expressed more complex transcriptomes with the vast majority of genes expressed and relatively small contributions from the most expressed genes. mRNAs expressed in brain had unusually long 3′UTRs, and mean 3′UTR length was higher for genes involved in development, morphogenesis and signal transduction, suggesting added complexity of UTR-based regulation for these genes. Our results support a model in which variable exterior components feed into a large, densely connected core composed of ubiquitously expressed intracellular proteins., National Institutes of Health (U.S.), Knut and Alice Wallenbergs Foundation, Ake Wibergs Foundation

Through alternative processing of pre-messenger RNAs, individual mammalian genes often produce multiple mRNA and protein isoforms that may have related, distinct or even opposing functions. Here we report an in-depth analysis of 15 diverse human tissue and cell line transcriptomes on the basis of deep sequencing of complementary DNA fragments, yielding a digital inventory of gene and mRNA isoform expression. Analyses in which sequence reads are mapped to exon-exon junctions indicated that 92-94% of human genes undergo alternative splicing, 86% with a minor isoform frequency of 15% or more. Differences in isoform-specific read densities indicated that most alternative splicing and alternative cleavage and polyadenylation events vary between tissues, whereas variation between individuals was approximately twofold to threefold less common. Extreme or 'switch-like' regulation of splicing between tissues was associated with increased sequence conservation in regulatory regions and with generation of full-length open reading frames. Patterns of alternative splicing and alternative cleavage and polyadenylation were strongly correlated across tissues, suggesting coordinated regulation of these processes, and sequence conservation of a subset of known regulatory motifs in both alternative introns and 3' untranslated regions suggested common involvement of specific factors in tissue-level regulation of both splicing and polyadenylation.

The relationship of the dopamine D4 receptor gene (DRD4) to the behavioral trait of novelty seeking has not been uniformly consistent. A methodological shortcoming in previous studies may relate to the way different DRD4 variants were categorized. Because of evolutionary and functional (e.g., diminished potency to reduce cAMP) similarities between the 2- and 7-repeat (2R, 7R) alleles of the DRD4, we suggest grouping of these two alleles together may facilitate detection of biologically meaningful and reproducible association findings with behavioral traits. We measured novelty seeking with the Tridimensional Personality Questionnaire (TPQ) in a community sample of Caucasian, Korean, and Filipino subjects (N = 171) who were subsequently characterized for the DRD4 variable number of tandem repeats (VNTR). In the Korean sample, those with a 2R and/or 7R allele scored significantly higher on novelty seeking scale (P < 0.05). By contrast, grouping the VNTR alleles by size (2, 3, 4 vs. 5, 6, 7), as has been done in similar studies of Asian subjects, was not significant. Using the extreme discordant phenotype (EDP) strategy in the pooled sample and selecting the individuals within the upper and lower decile, we observed a trend for association with higher novelty seeking in individuals who carry the 2R and/or 7R alleles (P = 0.06). We also confirmed that the 2R allele in the Korean and Filipino subjects was the result of a one-step recombination event between the 4R and 7R alleles. This study suggests that genetic association analyses can benefit by consideration of the shared functional and evolutionary attributes of the DRD4 2R and 7R alleles.

Using the 2.6 million single nucleotide polymorphism (SNP) genotype datasets from Perlegen Sciences and the Haplotype Map (HapMap) project (Phase I freeze), a probabilistic search for the landscape exhibited by positive Darwinian selection was conducted (Wang et al., 2006). By sorting each high frequency allele by homozygosity, we search for the expected decay of adjacent SNP linkage disequilibrium (LD) at recently selected alleles, eliminating the need for inferring haplotype. We designate this approach the LD decay (LDD) test. Cluster analysis indicates that approximately 3000 sites of recent inferred selection are present in human DNA, representing approximately 1800 genes. Prior simulation studies (Wang et al., 2006) indicate that this novel LDD test, at the Mb scale employed, effectively distinguishes selection from other causes of extensive LD, such as inversions, population bottlenecks and admixture. Based on over-representation analysis, these prior studies have shown that several predominant biological themes are common in inferred selected alleles, including genes involved with DNA metabolism and repair. Here, we show that three of these DNA repair genes, ERCC8, Fanconi Anemia Complementation Group C (FANCC), and RAD51C, exhibit genomic architectures consistent with ongoing balanced selection over the last 40,000-50,000 years.

The TBestDB database contains approximately 370,000 clustered expressed sequence tag (EST) sequences from 49 organisms, covering a taxonomically broad range of poorly studied, mainly unicellular eukaryotes, and includes experimental information, consensus sequences, gene annotations and metabolic pathway predictions. Most of these ESTs have been generated by the Protist EST Program, a collaboration among six Canadian research groups. EST sequences are read from trace files up to a minimum quality cut-off, vector and linker sequence is masked, and the ESTs are clustered using phrap. The resulting consensus sequences are automatically annotated by using the AutoFACT program. The datasets are automatically checked for clustering errors due to chimerism and potential cross-contamination between organisms, and suspect data are flagged in or removed from the database. Access to data deposited in TBestDB by individual users can be restricted to those users for a limited period. With this first report on TBestDB, we open the database to the research community for free processing, annotation, interspecies comparisons and GenBank submission of EST data generated in individual laboratories. For instructions on submission to TBestDB, contact tbestdb@bch.umontreal.ca. The database can be queried at http://tbestdb.bcm.umontreal.ca/.

Objective: Our purpose was to compare the costs of prenatal care and subsequent maternal and neonatal outcomes in patients with gestational diabetes cared for in an inner-city university hospital house staff clinic versus an inner-city managed care organization. Study Design: A retrospective cohort study was conducted. The groups consisted of 115 patients with gestational diabetes who were cared for in a house staff clinic and a demographically similar group of 85 patients cared for in a neighborhood managed care organization. The groups were examined regarding baseline demographics, intensity of prenatal care, maternal and neonatal outcomes, and total cost of the provision of care. Results: There was no difference between groups in the total cost of maternal-infant care. A larger percentage of patients in the house staff group saw the physician frequently. In contrast, patients cared for in the managed care organization underwent more tests of fetal well-being. There was a greater rate of neonatal macrosomia in the managed care organization group compared with the house staff group. Conclusions: Managed care does not decrease the cost of caring for patients with gestational diabetes but does lead to a greater rate of neonatal macrosomia, which may reflect poorer glucose control.

By using the 1.6 million single-nucleotide polymorphism (SNP) genotype data set from Perlegen Sciences [Hinds, D. A., Stuve, L. L., Nilsen, G. B., Halperin, E., Eskin, E., Ballinger, D. G., Frazer, K. A. & Cox, D. R. (2005) Science 307, 1072–1079], a probabilistic search for the landscape exhibited by positive Darwinian selection was conducted. By sorting each high-frequency allele by homozygosity, we search for the expected decay of adjacent SNP linkage disequilibrium (LD) at recently selected alleles, eliminating the need for inferring haplotype. We designate this approach the LD decay (LDD) test. By these criteria, 1.6% of Perlegen SNPs were found to exhibit the genetic architecture of selection. These results were confirmed on an independently generated data set of 1.0 million SNP genotypes (International Human Haplotype Map Phase I freeze). Simulation studies indicate that the LDD test, at the megabase scale used, effectively distinguishes selection from other causes of extensive LD, such as inversions, population bottlenecks, and admixture. The ≈1,800 genes identified by the LDD test were clustered according to Gene Ontology (GO) categories. Based on overrepresentation analysis, several predominant biological themes are common in these selected alleles, including host–pathogen interactions, reproduction, DNA metabolism/cell cycle, protein metabolism, and neuronal function.

Associations have been reported of the 7-repeat (7R) allele of the human dopamine receptor D4 (DRD4) gene with both the personality trait of novelty seeking and attention-deficit/hyperactivity disorder (ADHD). The increased prevalence of the 7R allele in ADHD probands is consistent with the common variant-common disorder hypothesis, which proposes that the high frequency of many complex genetic disorders is related to common DNA variants. Recently, based on the unusual DNA sequence organization and strong linkage disequilibrium surrounding the DRD4 7R allele, we proposed that this allele originated as a rare mutational event, which nevertheless increased to high prevalence in human populations by positive selection. We have now determined, by DNA resequencing of 250 DRD4 alleles obtained from 132 ADHD probands, that most ADHD 7R alleles are of the conserved haplotype found in our previous 600 allele worldwide DNA sample. Interestingly, however, half of the 24 haplotypes uncovered in ADHD probands were novel (not one of the 56 haplotypes found in our prior population studies). Over 10 percent of the ADHD probands had these novel haplotypes, most of which were 7R allele derived. The probability that this high incidence of novel alleles occurred by chance in our ADHD sample is much less than 0.0001. These results suggest that allelic heterogeneity at the DRD4 locus may also contribute to the observed association with ADHD.

The nuclear receptor PXR (pregnane X receptor) is a broad-specificity sensor that recognizes a wide variety of synthetic drugs and xenobiotic agents. On activation by these compounds, PXR coordinately induces a network of transporters, cytochrome P450 enzymes, and other genes that effectively clear xenobiotics from the liver and intestine. Like PXR, the majority of its target genes also possess a broad specificity for exogenous compounds. Thus, PXR is both a sensor and effector in a well integrated and generalized pathway for chemical immunity. Although it is clear that PXR responds to numerous foreign compounds, it is unclear whether it possesses an endogenous ligand. To address this issue, we noted that there is substantial overlap in the substrate specificities of PXR and its critical CYP3A target gene. This prompted us to ask whether endogenous CYP3A substrates also serve as PXR ligands. We demonstrate that 5β-cholestane-3α,7α,12α-triol (triol), a cholesterol-derived CYP3A substrate, is a potent PXR agonist that effectively induces cyp3a expression in mice. This defines a critical salvage pathway that can be autoinduced to minimize triol accumulation. In contrast, triol can accumulate to very high levels in humans, and unlike mice, these people develop the severe clinical manifestations of cerebrotendinous xanthomatosis. The reason for these dramatic species differences has remained unclear. We now demonstrate that triol fails to activate human PXR or induce the CYP3A -salvage pathway. This explains why humans are more susceptible to sterol accumulation and suggests that synthetic ligands for human PXR could be used to treat cerebrotendinous xanthomatosis and other disorders of cholesterol excess.

Expansions of simple repetitive DNA sequences are responsible for a large number of human hereditary diseases. The severity and age of onset of each disorder can be estimated from the length of the inherited repeat tract; however, the subsequent expansion rate is itself a heritable trait. This points to a role for trans-acting modifier genes, whose actions may be altered by single nucleotide polymorphisms. Previously, our lab has used a yeast system to study large-scale expansions of GAA trinucleotide repeats, and a number of gene knockouts have been shown to induce expansions. While this has been highly informative, gene knockouts are a rather blunt tool, which are not suitable to approximate the spectrum of SNPs that might affect repeat expansions, as has been predicted in human populations. To address this problem, we have developed a novel, alternative screening method which employs UV mutagenesis followed by the whole-genome sequencing. This approach allowed us to detect SNPs that modified the likeli...