Dynamic integration of splicing within gene regulatory pathways.

Precursor mRNA splicing is one of the most highly regulated processes in metazoan species. In addition to generating vast repertoires of RNAs and proteins, splicing has a profound impact on other gene regulatory layers, including mRNA transcription, turnover, transport, and translation. Conversely, factors regulating chromatin and transcription complexes impact the splicing process. This extensive crosstalk between gene regulatory layers takes advantage of dynamic spatial, physical, and temporal organizational properties of the cell nucleus, and further emphasizes the importance of developing a multidimensional understanding of splicing control.

Early Detection of Severe Acute Respiratory Syndrome Coronavirus 2 Variants Using Traveler-based Genomic Surveillance at 4 US Airports, September 2021-January 2022.

We enrolled arriving international air travelers in a severe acute respiratory syndrome coronavirus 2 genomic surveillance program. We used molecular testing of pooled nasal swabs and sequenced positive samples for sublineage. Traveler-based surveillance provided early-warning variant detection, reporting the first US Omicron BA.2 and BA.3 in North America.

Genome-wide identification of zero nucleotide recursive splicing in Drosophila.

Recursive splicing is a process in which large introns are removed in multiple steps by re-splicing at ratchet points--5' splice sites recreated after splicing. Recursive splicing was first identified in the Drosophila Ultrabithorax (Ubx) gene and only three additional Drosophila genes have since been experimentally shown to undergo recursive splicing. Here we identify 197 zero nucleotide exon ratchet points in 130 introns of 115 Drosophila genes from total RNA sequencing data generated from developmental time points, dissected tissues and cultured cells. The sequential nature of recursive splicing was confirmed by identification of lariat introns generated by splicing to and from the ratchet points. We also show that recursive splicing is a constitutive process, that depletion of U2AF inhibits recursive splicing, and that the sequence and function of ratchet points are evolutionarily conserved in Drosophila. Finally, we identify four recursively spliced human genes, one of which is also recursively spliced in Drosophila. Together, these results indicate that recursive splicing is commonly used in Drosophila, occurs in humans, and provides insight into the mechanisms by which some large introns are removed.

New insights from existing sequence data: generating breakthroughs without a pipette.

With the rapidly declining cost of data generation and the accumulation of massive data sets, molecular biology is entering an era in which incisive analysis of existing data will play an increasingly prominent role in the discovery of new biological phenomena and the elucidation of molecular mechanisms. Here, we discuss resources of publicly available sequencing data most useful for interrogating the mechanisms of gene expression. Existing next-generation sequence data sets, however, come with significant challenges in the form of technical and bioinformatic artifacts, which we discuss in detail. We also recount several breakthroughs made largely through the analysis of existing data, primarily in the RNA field.

Diverse forms of RPS9 splicing are part of an evolving autoregulatory circuit.

Ribosomal proteins are essential to life. While the functions of ribosomal protein-encoding genes (RPGs) are highly conserved, the evolution of their regulatory mechanisms is remarkably dynamic. In Saccharomyces cerevisiae, RPGs are unusual in that they are commonly present as two highly similar gene copies and in that they are over-represented among intron-containing genes. To investigate the role of introns in the regulation of RPG expression, we constructed 16 S. cerevisiae strains with precise deletions of RPG introns. We found that several yeast introns function to repress rather than to increase steady-state mRNA levels. Among these, the RPS9A and RPS9B introns were required for cross-regulation of the two paralogous gene copies, which is consistent with the duplication of an autoregulatory circuit. To test for similar intron function in animals, we performed an experimental test and comparative analyses for autoregulation among distantly related animal RPS9 orthologs. Overexpression of an exogenous RpS9 copy in Drosophila melanogaster S2 cells induced alternative splicing and degradation of the endogenous copy by nonsense-mediated decay (NMD). Also, analysis of expressed sequence tag data from distantly related animals, including Homo sapiens and Ciona intestinalis, revealed diverse alternatively-spliced RPS9 isoforms predicted to elicit NMD. We propose that multiple forms of splicing regulation among RPS9 orthologs from various eukaryotes operate analogously to translational repression of the alpha operon by S4, the distant prokaryotic ortholog. Thus, RPS9 orthologs appear to have independently evolved variations on a fundamental autoregulatory circuit.

Differences between homologous alleles of olfactory receptor genes require the Polycomb Group protein Eed.

A number of mammalian genes are expressed from only one of the two homologous chromosomes, selected at random in each cell. These include genes subject to X-inactivation, olfactory receptor (OR) genes, and several classes of immune system genes. The means by which monoallelic expression is established are only beginning to be understood. Using a cytological assay, we show that the two homologous alleles of autosomal random monoallelic loci differ from each other in embryonic stem (ES) cells, before establishment of monoallelic expression. The Polycomb Group gene Eed is required to establish this distinctive behavior. In addition, we found that when Eed mutant ES cells are differentiated, they fail to establish asynchronous replication timing at OR loci. These results suggest a common mechanism for random monoallelic expression on autosomes and the X chromosome, and implicate Eed in establishing differences between homologous OR loci before and after differentiation.

Lineage abundance estimation for SARS-CoV-2 in wastewater using transcriptome quantification techniques.

Effectively monitoring the spread of SARS-CoV-2 mutants is essential to efforts to counter the ongoing pandemic. Predicting lineage abundance from wastewater, however, is technically challenging. We show that by sequencing SARS-CoV-2 RNA in wastewater and applying algorithms initially used for transcriptome quantification, we can estimate lineage abundance in wastewater samples. We find high variability in signal among individual samples, but the overall trends match those observed from sequencing clinical samples. Thus, while clinical sequencing remains a more sensitive technique for population surveillance, wastewater sequencing can be used to monitor trends in mutant prevalence in situations where clinical sequencing is unavailable.

Variant abundance estimation for SARS-CoV-2 in wastewater using RNA-Seq quantification.

Effectively monitoring the spread of SARS-CoV-2 variants is essential to efforts to counter the ongoing pandemic. Wastewater monitoring of SARS-CoV-2 RNA has proven an effective and efficient technique to approximate COVID-19 case rates in the population. Predicting variant abundances from wastewater, however, is technically challenging. Here we show that by sequencing SARS-CoV-2 RNA in wastewater and applying computational techniques initially used for RNA-Seq quantification, we can estimate the abundance of variants in wastewater samples. We show by sequencing samples from wastewater and clinical isolates in Connecticut U.S.A. between January and April 2021 that the temporal dynamics of variant strains broadly correspond. We further show that this technique can be used with other wastewater sequencing techniques by expanding to samples taken across the United States in a similar timeframe. We find high variability in signal among individual samples, and limited ability to detect the presence of variants with clinical frequencies <10%; nevertheless, the overall trends match what we observed from sequencing clinical samples. Thus, while clinical sequencing remains a more sensitive technique for population surveillance, wastewater sequencing can be used to monitor trends in variant prevalence in situations where clinical sequencing is unavailable or impractical.

Transcriptome and in Vitro Differentiation Profile of Human Embryonic Stem Cell Derived NKX2.1-Positive Neural Progenitors.

The generation of inhibitory interneuron progenitors from human embryonic stem cells (ESCs) is of great interest due to their potential use in transplantation therapies designed to treat central nervous system disorders. The medial ganglionic eminence (MGE) is a transient embryonic structure in the ventral telencephalon that is a major source of cortical GABAergic inhibitory interneuron progenitors. These progenitors migrate tangentially to sites in the cortex and differentiate into a variety of interneuron subtypes, forming local synaptic connections with excitatory projection neurons to modulate activity of the cortical circuitry. The homeobox domain-containing transcription factor NKX2.1 is highly expressed in the MGE and pre-optic area of the ventral subpallium and is essential for specifying cortical interneuron fate. Using a combination of growth factor agonists and antagonists to specify ventral telencephalic fates, we previously optimized a protocol for the efficient generation of NKX2.1-positive MGE-like neural progenitors from human ESCs. To establish their identity, we now characterize the transcriptome of these MGE-like neural progenitors using RNA sequencing and demonstrate the capacity of these cells to differentiate into inhibitory interneurons in vitro using a neuron-astrocyte co-culture system. These data provide information on the potential origin of interneurons in the human brain.

COMT genetic variation confers risk for psychotic and affective disorders: a case control study.

BACKGROUND: Variation in the COMT gene has been implicated in a number of psychiatric disorders, including psychotic, affective and anxiety disorders. The majority of these studies have focused on the functional Val108/158Met polymorphism and yielded conflicting results, with limited studies examining the relationship between other polymorphisms, or haplotypes, and psychiatric illness. We hypothesized that COMT variation may confer a general risk for psychiatric disorders and have genotyped four COMT variants (Val158Met, rs737865, rs165599, and a SNP in the P2 promoter [-278A/G; rs2097603]) in 394 Caucasian cases and 467 controls. Cases included patients with schizophrenia (n = 196), schizoaffective disorder (n = 62), bipolar disorder (n = 82), major depression (n = 30), and patients diagnosed with either psychotic disorder NOS or depressive disorder NOS (n = 24). RESULTS: SNP rs2097603, the Val/Met variant and SNP rs165599 were significantly associated (p = 0.004; p = 0.05; p = 0.035) with a broad "all affected" diagnosis. Haplotype analysis revealed a potentially protective G-A-A-A haplotype haplotype (-278A/G; rs737865; Val108/158Met; rs165599), which was significantly underrepresented in this group (p = 0.0033) and contained the opposite alleles of the risk haplotype previously described by Shifman et al. Analysis of diagnostic subgroups within the "all affecteds group" showed an association of COMT in patients with psychotic disorders as well as in cases with affective illness although the associated variants differed. The protective haplotype remained significantly underrepresented in most of these subgroups. CONCLUSION: Our results support the view that COMT variation provides a weak general predisposition to neuropsychiatric disease including psychotic and affective disorders.

Gene expression analysis of human induced pluripotent stem cell-derived neurons carrying copy number variants of chromosome 15q11-q13.1.

BACKGROUND: Duplications of the chromosome 15q11-q13.1 region are associated with an estimated 1 to 3% of all autism cases, making this copy number variation (CNV) one of the most frequent chromosome abnormalities associated with autism spectrum disorder (ASD). Several genes located within the 15q11-q13.1 duplication region including ubiquitin protein ligase E3A (UBE3A), the gene disrupted in Angelman syndrome (AS), are involved in neural function and may play important roles in the neurobehavioral phenotypes associated with chromosome 15q11-q13.1 duplication (Dup15q) syndrome. METHODS: We have generated induced pluripotent stem cell (iPSC) lines from five different individuals containing CNVs of 15q11-q13.1. The iPSC lines were differentiated into mature, functional neurons. Gene expression across the 15q11-q13.1 locus was compared among the five iPSC lines and corresponding iPSC-derived neurons using quantitative reverse transcription PCR (qRT-PCR). Genome-wide gene expression was compared between neurons derived from three iPSC lines using mRNA-Seq. RESULTS: Analysis of 15q11-q13.1 gene expression in neurons derived from Dup15q iPSCs reveals that gene copy number does not consistently predict expression levels in cells with interstitial duplications of 15q11-q13.1. mRNA-Seq experiments show that there is substantial overlap in the genes differentially expressed between 15q11-q13.1 deletion and duplication neurons, Finally, we demonstrate that UBE3A transcripts can be pharmacologically rescued to normal levels in iPSC-derived neurons with a 15q11-q13.1 duplication. CONCLUSIONS: Chromatin structure may influence gene expression across the 15q11-q13.1 region in neurons. Genome-wide analyses suggest that common neuronal pathways may be disrupted in both the Angelman and Dup15q syndromes. These data demonstrate that our disease-specific stem cell models provide a new tool to decipher the underlying cellular and genetic disease mechanisms of ASD and may also offer a pathway to novel therapeutic intervention in Dup15q syndrome.

Reciprocal intronic and exonic histone modification regions in humans.

While much attention has been focused on chromatin at promoters and exons, human genes are mostly composed of intronic sequences. Analyzing published surveys of nucleosomes and 41 chromatin marks in humans, we identified histone modifications specifically associated with 5' intronic sequences, distinguishable from promoter marks and bulk nucleosomes. These intronic marks were spatially reciprocal to trimethylated histone H3 Lys36 (H3K36me3), typically transitioning near internal exons. Several marks transitioned near bona fide exons, but not near nucleosomes at exon-like sequences. Therefore, we examined whether splicing affects histone marking. Even with considerable changes in regulated alternative splicing, histone marks were stable. Notably, these findings are consistent with exon definition influencing histone marks. In summary, we show that the location of many intragenic marks in humans can be distilled into a simple organizing principle: association with 5' intronic or 3' exonic regions.

A role for the CHC22 clathrin heavy-chain isoform in human glucose metabolism.

Intracellular trafficking of the glucose transporter GLUT4 from storage compartments to the plasma membrane is triggered in muscle and fat during the body's response to insulin. Clathrin is involved in intracellular trafficking, and in humans, the clathrin heavy-chain isoform CHC22 is highly expressed in skeletal muscle. We found a role for CHC22 in the formation of insulin-responsive GLUT4 compartments in human muscle and adipocytes. CHC22 also associated with expanded GLUT4 compartments in muscle from type 2 diabetic patients. Tissue-specific introduction of CHC22 in mice, which have only a pseudogene for this protein, caused aberrant localization of GLUT4 transport pathway components in their muscle, as well as features of diabetes. Thus, CHC22-dependent membrane trafficking constitutes a species-restricted pathway in human muscle and fat with potential implications for type 2 diabetes.

Genetic diversity and genomic distribution of homologs encoding NBS-LRR disease resistance proteins in sunflower.

Three-fourths of the recognition-dependent disease resistance genes (R-genes) identified in plants encode nucleotide binding site (NBS) leucine-rich repeat (LRR) proteins. NBS-LRR homologs have only been isolated on a limited scale from sunflower (Helianthus annuus L.), and most of the previously identified homologs are members of two large NBS-LRR clusters harboring downy mildew R-genes. We mined the sunflower EST database and used comparative genomics approaches to develop a deeper understanding of the diversity and distribution of NBS-LRR homologs in the sunflower genome. Collectively, 630 NBS-LRR homologs were identified, 88 by mining a database of 284,241 sunflower ESTs and 542 by sequencing 1,248 genomic DNA amplicons isolated from common and wild sunflower species. DNA markers were developed from 196 unique NBS-LRR sequences and facilitated genetic mapping of 167 NBS-LRR loci. The latter were distributed throughout the sunflower genome in 44 clusters or singletons. Wild species ESTs were a particularly rich source of novel NBS-LRR homologs, many of which were tightly linked to previously mapped downy mildew, rust, and broomrape R-genes. The DNA sequence and mapping resources described here should facilitate the discovery and isolation of recognition-dependent R-genes guarding sunflower from a broad spectrum of economically important diseases. Sunflower nucleotide and amino acid sequences have been deposited in DDBJ/EMBL/GenBank under accession numbers EF 560168-EF 559378 and ABQ 58077-ABQ 57529.

Analysis of TBX1 variation in patients with psychotic and affective disorders.

A significant portion of patients with 22q11 deletion syndrome (22q11DS) develop psychiatric disorders, including schizophrenia and other psychotic and affective symptoms, and the responsible gene/s are assumed to also play a significant role in the etiology of nonsyndromic psychiatric disease. The most common psychiatric diagnosis among patients with 22q11DS is schizophrenia, thought to result from neurotransmitter imbalances and also from disturbed brain development. Several genes in the 22q11 region with known or suspected roles in neurotransmitter metabolism have been analyzed in patients with isolated schizophrenia; however, their contribution to the disease remains controversial. Haploinsufficiency of the TBX1 gene has been shown to be sufficient to cause the core physical malformations associated with 22q11DS in mice and humans and via abnormal brain development could contribute to 22q11DS-related and isolated psychiatric disease. 22q11DS populations also have increased rates of psychiatric conditions other than schizophrenia, including mood disorders. We therefore analyzed variations at the TBX1 locus in a cohort of 446 white patients with psychiatric disorders relevant to 22q11DS and 436 ethnically matched controls. The main diagnoses included schizophrenia (n = 226), schizoaffective disorder (n = 67), bipolar disorder (n = 82), and major depressive disorder (n = 29). We genotyped nine tag SNPs in this sample but did not observe significant differences in allele or haplotype frequencies in any of the analyzed groups (all affected, schizophrenia and schizoaffective disorder, schizophrenia alone, and bipolar disorder and major depressive disorder) compared with the control group. Based on these results we conclude that TBX1 variation does not make a strong contribution to the genetic etiology of nonsyndromic forms of psychiatric disorders commonly seen in patients with 22q11DS.

Comparative analysis of NBS domain sequences of NBS-LRR disease resistance genes from sunflower, lettuce, and chicory.

Plant resistance to many types of pathogens and pests can be achieved by the presence of disease resistance (R) genes. The nucleotide binding site-leucine rich repeat (NBS-LRR) class of R-genes is the most commonly isolated class of R-genes and makes up a super-family, which is often arranged in the genome as large multi-gene clusters. The NBS domain of these genes can be targeted by polymerase chain reaction (PCR) amplification using degenerate primers. Previous studies have used PCR derived NBS sequences to investigate both ancient R-gene evolution and recent evolution within specific plant families. However, comparative studies with the Asteraceae family have largely been ignored. In this study, we address recent evolution of NBS sequences within the Asteraceae and extend the comparison to the Arabidopsis thaliana genome. Using multiple sets of primers, NBS fragments were amplified from genomic DNA of three species from the family Asteraceae: Helianthus annuus (sunflower), Lactuca sativa (lettuce), and Cichorium intybus (chicory). Analysis suggests that Asteraceae species share distinct families of R-genes, composed of genes related to both coiled-coil (CC) and toll-interleukin-receptor homology (TIR) domain containing NBS-LRR R-genes. Between the most closely related species, (lettuce and chicory) a striking similarity of CC subfamily composition was identified, while sunflower showed less similarity in structure. These sequences were also compared to the A. thaliana genome. Asteraceae NBS gene subfamilies appear to be distinct from Arabidopsis gene clades. These data suggest that NBS families in the Asteraceae family are ancient, but also that gene duplication and gene loss events occur and change the composition of these gene subfamilies over time.

Association of the DTNBP1 locus with schizophrenia in a U.S. population.

Linkage and association studies have recently implicated dystrobrevin-binding protein 1 (DTNBP1) in the etiology of schizophrenia. We analyzed seven previously tested DTNBP1 single-nucleotide polymorphisms (SNPs) in a cohort of 524 individuals with schizophrenia or schizoaffective disorder and 573 control subjects. The minor alleles of three SNPs (P1578, P1763, and P1765) were positively associated with the diagnosis of schizophrenia or schizoaffective disorder in the white subset of the study cohort (258 cases, 467 controls), with P1578 showing the most significant association (odds ratio 1.76, P =.0026). The same three SNPs were also associated in a smaller Hispanic subset (51 cases, 32 controls). No association was observed in the African American subset (215 cases, 74 controls). A stratified analysis of the white and Hispanic subsets showed association with the minor alleles of four SNPs (P1578, P1763, P1320, and P1765). Again, the most significant association was observed for P1578 (P =.0006). Haplotype analysis supported these findings, with a single risk haplotype significantly overrepresented in the white sample (P =.005). Our study provides further evidence for a role of the DTNBP1 gene in the genetic etiology of schizophrenia.