The Expanding Landscape of Alternative Splicing Variation in Human Populations.

Alternative splicing is a tightly regulated biological process by which the number of gene products for any given gene can be greatly expanded. Genomic variants in splicing regulatory sequences can disrupt splicing and cause disease. Recent developments in sequencing technologies and computational biology have allowed researchers to investigate alternative splicing at an unprecedented scale and resolution. Population-scale transcriptome studies have revealed many naturally occurring genetic variants that modulate alternative splicing and consequently influence phenotypic variability and disease susceptibility in human populations. Innovations in experimental and computational tools such as massively parallel reporter assays and deep learning have enabled the rapid screening of genomic variants for their causal impacts on splicing. In this review, we describe technological advances that have greatly increased the speed and scale at which discoveries are made about the genetic variation of alternative splicing. We summarize major findings from population transcriptomic studies of alternative splicing and discuss the implications of these findings for human genetics and medicine

Identification of novel components of Trypanosoma brucei editosomes

The editosome is a multiprotein complex that catalyzes the insertion and deletion of uridylates that occurs during RNA editing in trypanosomatids. We report the identification of nine novel editosome proteins in Trypanosoma brucei. They were identified by mass spectrometric analysis of functional editosomes that were purified by serial ion exchange/gel permeation chromatography, immunoaffinity chromatography specific to the TbMP63 editosome protein, or tandem affinity purification based on a tagged RNA editing ligase. The newly identified proteins have ribonuclease and/or RNA binding motifs suggesting nuclease function for at least some of these. Five of the proteins are interrelated, as are two others, and one is related to four previously identified editosome proteins. The implications of these findings are discussed

Predicting variation of DNA shape preferences in protein-DNA interaction in cancer cells with a new biophysical model

DNA shape readout is an important mechanism of target site recognition by transcription factors, in addition to the sequence readout. Several models of transcription factor-DNA binding which consider DNA shape have been developed in recent years. We present a new biophysical model of protein-DNA interaction by considering the DNA shape features, which is based on a neighbour dinucleotide dependency model BayesPI2. The parameters of the new model are restricted to a subspace spanned by the 2-mer DNA shape features, which allowing a biophysical interpretation of the new parameters as position-dependent preferences towards certain values of the features. Using the new model, we explore the variation of DNA shape preferences in several transcription factors across cancer cell lines and cellular conditions. We find evidence of DNA shape variations at FOXA1 binding sites in MCF7 cells after treatment with steroids. The new model is useful for elucidating finer details of transcription factor-DNA interaction. It may be used to improve the prediction of cancer mutation effects in the future

Genomic characterization of Gli-activator targets in sonic hedgehog-mediated neural patterning

Sonic hedgehog (Shh) acts as a morphogen to mediate the specification of distinct cell identities in the ventral neural tube through a Gli-mediated (Gli1-3) transcriptional network. Identifying Gli targets in a systematic fashion is central to the understanding of the action of Shh. We examined this issue in differentiating neural progenitors in mouse. An epitope-tagged Gli-activator protein was used to directly isolate cis-regulatory sequences by chromatin immunoprecipitation (ChIP). ChIP products were then used to screen custom genomic tiling arrays of putative Hedgehog (Hh) targets predicted from transcriptional profiling studies, surveying 50-150 kb of non-transcribed sequence for each candidate. In addition to identifying expected Gli-target sites, the data predicted a number of unreported direct targets of Shh action. Transgenic analysis of binding regions in Nkx2.2, Nkx2.1 (Titf1) and Rab34 established these as direct Hh targets. These data also facilitated the generation of an algorithm that improved in silico predictions of Hh target genes. Together, these approaches provide significant new insights into both tissue-specific and general transcriptional targets in a crucial Shh-mediated patterning process

Potent CRISPR-Cas9 inhibitors from Staphylococcus genomes.

Anti-CRISPRs (Acrs) are small proteins that inhibit the RNA-guided DNA targeting activity of CRISPR-Cas enzymes. Encoded by bacteriophage and phage-derived bacterial genes, Acrs prevent CRISPR-mediated inhibition of phage infection and can also block CRISPR-Cas-mediated genome editing in eukaryotic cells. To identify Acrs capable of inhibiting Staphylococcus aureus Cas9 (SauCas9), an alternative to the most commonly used genome editing protein Streptococcus pyogenes Cas9 (SpyCas9), we used both self-targeting CRISPR screening and guilt-by-association genomic search strategies. Here we describe three potent inhibitors of SauCas9 that we name AcrIIA13, AcrIIA14, and AcrIIA15. These inhibitors share a conserved N-terminal sequence that is dispensable for DNA cleavage inhibition and have divergent C termini that are required in each case for inhibition of SauCas9-catalyzed DNA cleavage. In human cells, we observe robust inhibition of SauCas9-induced genome editing by AcrIIA13 and moderate inhibition by AcrIIA14 and AcrIIA15. We also find that the conserved N-terminal domain of AcrIIA13-AcrIIA15 binds to an inverted repeat sequence in the promoter of these Acr genes, consistent with its predicted helix-turn-helix DNA binding structure. These data demonstrate an effective strategy for Acr discovery and establish AcrIIA13-AcrIIA15 as unique bifunctional inhibitors of SauCas9

Multiple mRNA isoforms of the transcription activator protein CREB

We have characterized cDNA clones representing mouse CREB (cyclic AMP responsive element binding protein) mRNA isoforms. These include CREBA and CREBa, of which the rat and human homologues have been previously identified. Both encode proteins with CREbinding activity and identical transactivation potential. The additional CREB mRNA isoforms potentially encode CREB related proteins. From the structural organization of the mouse CREB gene we conclude that the multiple transcripts are generated by alternative splicing. Furthermore we show that specific CREB mRNA isoforms are expressed at a high level in the adult testis. Expression of these isoforms is induced after commencement of spermatogenesis. In situ hybridization suggests that this expression occurs predominantly in the primary spermatocytes. Comparison of the CREB gene with the recently isolated CREM (cAMP responsive element modulator) cDNAs illustrates that the two genes have arisen by gene duplication and have diverged to encode transcriptional activators and repressors of the cAMP signal transduction pathway

The Flagellar Motility of \u3cem\u3eChlamydomonas pf25\u3c/em\u3e Mutant Lacking an AKAP-binding Protein Is Overtly Sensitive to Medium Conditions

Radial spokes are a conserved axonemal structural complex postulated to regulate the motility of 9 + 2 cilia and flagella via a network of phosphoenzymes and regulatory proteins. Consistently, a Chlamydomonas radial spoke protein, RSP3, has been identified by RII overlays as an A-kinase anchoring protein (AKAP) that localizes the cAMP-dependent protein kinase (PKA) holoenzyme by binding to the RIIa domain of PKA RII subunit. However, the highly conserved docking domain of PKA is also found in the N termini of several AKAP-binding proteins unrelated to PKA as well as a 24-kDa novel spoke protein, RSP11. Here, we report that RSP11 binds to RSP3 directly in vitro and colocalizes with RSP3 toward the spoke base near outer doublets and dynein motors in axonemes. Importantly, RSP11 mutant pf25 displays a spectrum of motility, from paralysis with flaccid or twitching flagella as other spoke mutants to wild-typelike swimming. The wide range of motility changes reversibly depending on the condition of liquid media without replacing defective proteins. We postulate that radial spokes use the RIIa/AKAP module to regulate ciliary and flagellar beating; absence of the spoke RIIa protein exposes a medium-sensitive regulatory mechanism that is not obvious in wild-type Chlamydomonas

A flexible integrative approach based on random forest improves prediction of transcription factor binding sites

Transcription factor binding sites (TFBSs) are DNA sequences of 6-15 base pairs. Interaction of these TFBSs with transcription factors (TFs) is largely responsible for most spatiotemporal gene expression patterns. Here, we evaluate to what extent sequence-based prediction of TFBSs can be improved by taking into account the positional dependencies of nucleotides (NPDs) and the nucleotide sequence-dependent structure of DNA. We make use of the random forest algorithm to flexibly exploit both types of information. Results in this study show that both the structural method and the NPD method can be valuable for the prediction of TFBSs. Moreover, their predictive values seem to be complementary, even to the widely used position weight matrix (PWM) method. This led us to combine all three methods. Results obtained for five eukaryotic TFs with different DNA-binding domains show that our method improves classification accuracy for all five eukaryotic TFs compared with other approaches. Additionally, we contrast the results of seven smaller prokaryotic sets with high-quality data and show that with the use of high-quality data we can significantly improve prediction performance. Models developed in this study can be of great use for gaining insight into the mechanisms of TF binding