Scientific highlights

Transcriptome analysis has mainly relied on analyzing RNA sequencing data from whole cells, overlooking the impact of subcellular RNA localization and its influence on our understanding of gene function, and interpretation of gene expression signatures in cells. Here, we separated cytosolic and nuclear RNA from human fetal and adult brain samples and performed a comprehensive analysis of cytosolic and nuclear transcriptomes. There are significant differences in RNA expression for protein-coding and lncRNA genes between cytosol and nucleus. We show that transcripts encoding the nuclear-encoded mitochondrial proteins are significantly enriched in the cytosol compared to the rest of protein-coding genes. Differential expression analysis between fetal and adult frontal cortex show that results obtained from the cytosolic RNA differ from results using nuclear RNA both at the level of transcript types and the number of differentially expressed genes. Our data provide a resource for the subcellular localization of thousands of RNA transcripts in the human brain and highlight differences in using the cytosolic or the nuclear transcriptomes for expression analysis.

Numerous megafauna species from northern latitudes went extinct during the Pleistocene/Holocene transition as a result of climate-induced habitat changes. However, several ungulate species managed to successfully track their habitats during this period to eventually flourish and recolonise the holarctic regions. So far, the genomic impacts of these climate fluctuations on ungulates from high latitudes have been little explored. Here, we assemble a de-novo genome for the European moose (Alces alces) and analyse it together with re-sequenced nuclear genomes and ancient and modern mitogenomes from across the moose range in Eurasia and North America. We found that moose demographic history was greatly influenced by glacial cycles, with demographic responses to the Pleistocene/Holocene transition similar to other temperate ungulates. Our results further support that modern moose lineages trace their origin back to populations that inhabited distinct glacial refugia during the Last Glacial Maximum (LGM). Finally, we found that present day moose in Europe and North America show low to moderate inbreeding levels resulting from post-glacial bottlenecks and founder effects, but no evidence for recent inbreeding resulting from human-induced population declines. Taken together, our results highlight the dynamic recent evolutionary history of the moose and provide an important resource for further genomic studies.

One ongoing concern about CRISPR-Cas9 genome editing is that unspecific guide RNA (gRNA) binding may induce off-target mutations. However, accurate prediction of CRISPR-Cas9 off-target activity is challenging. Here, we present SMRT-OTS and Nano-OTS, two novel, amplification-free, long-read sequencing protocols for detection of gRNA-driven digestion of genomic DNA by Cas9 in vitro. The methods are assessed using the human cell line HEK293, re-sequenced at 18x coverage using highly accurate HiFi SMRT reads. SMRT-OTS and Nano-OTS are first applied to three different gRNAs targeting HEK293 genomic DNA, resulting in a set of 55 high-confidence gRNA cleavage sites identified by both methods. Twenty-five of these sites are not reported by off-target prediction software, either because they contain four or more single nucleotide mismatches or insertion/deletion mismatches, as compared with the human reference. Additional experiments reveal that 85% of Cas9 cleavage sites are also found by other in vitro-based methods and that on- and off-target sites are detectable in gene bodies where short-reads fail to uniquely align. Even though SMRT-OTS and Nano-OTS identify several sites with previously validated off-target editing activity in cells, our own CRISPR-Cas9 editing experiments in human fibroblasts do not give rise to detectable off-target mutations at the in vitro-predicted sites. However, indel and structural variation events are enriched at the on-target sites. Amplification-free long-read sequencing reveals Cas9 cleavage sites in vitro that would have been difficult to predict using computational tools, including in dark genomic regions inaccessible by short-read sequencing.

Despite considerable progress in schizophrenia genetics, most findings have been for large rare structural variants and common variants in well-imputed regions with few genes implicated from exome sequencing. Whole genome sequencing (WGS) can potentially provide a more complete enumeration of etiological genetic variation apart from the exome and regions of high linkage disequilibrium. We analyze high-coverage WGS data from 1162 Swedish schizophrenia cases and 936 ancestry-matched population controls. Our main objective is to evaluate the contribution to schizophrenia etiology from a variety of genetic variants accessible to WGS but not by previous technologies. Our results suggest that ultra-rare structural variants that affect the boundaries of topologically associated domains (TADs) increase risk for schizophrenia. Alterations in TAD boundaries may lead to dysregulation of gene expression. Future mechanistic studies will be needed to determine the precise functional effects of these variants on biology.

Angioedema occurring in the head and neck region is a rare and sometimes life-threatening adverse reaction to angiotensin-converting enzyme inhibitors (ACEIs) and angiotensin receptor blockers (ARBs). Few studies have investigated the association of common variants with this extreme reaction, but none have explored the combined influence of rare variants yet. Adjudicated cases of ACEI-induced angioedema (ACEI-AE) or ARB-induced angioedema (ARB-AE) and controls were recruited at five different centers. Sequencing of 1,066 samples (408 ACEI-AE, ARB-AE, and 658 controls) was performed using exome-enriched sequence data. A common variant of the F5 gene that causes an increase in blood clotting (rs6025, p.Arg506Gln, also called factor V Leiden), was significantly associated with both ACEI-AE and ARB-AE (odds ratio: 2.85, 95% confidence interval (CI), 1.89-4.25). A burden test analysis of five rare missense variants in F5 was also found to be associated with ACEI-AE or ARB-AE, P = 2.09 × 10 -3 . A combined gene risk score of these variants, and the common variants rs6025 and rs6020, showed that individuals carrying at least one variant had 2.21 (95% CI, 1.49-3.27, P = 6.30 × 10-9 ) times the odds of having ACEI-AE or ARB-AE. The increased risk due to the common Leiden allele was confirmed in a genome-wide association study from the United States. A high risk of angioedema was also observed for the rs6020 variant that is the main coagulation defect-causing variant in black African and Asian populations. We found that deleterious missense variants in F5 are associated with an increased risk of ACEI-AE or ARB-AE.

A meta-analysis of genome-wide association studies (GWAS) identified eight loci that are associated with heart rate variability (HRV), but candidate genes in these loci remain uncharacterized. We developed an image- and CRISPR/Cas9-based pipeline to systematically characterize candidate genes for HRV in live zebrafish embryos. Nine zebrafish orthologues of six human candidate genes were targeted simultaneously in eggs from fish that transgenically express GFP on smooth muscle cells (Tg[acta2:GFP]), to visualize the beating heart. An automated analysis of repeated 30 s recordings of beating atria in 381 live, intact zebrafish embryos at 2 and 5 days post-fertilization highlighted genes that influence HRV (hcn4 and si:dkey-65j6.2 [KIAA1755]); heart rate (rgs6 and hcn4); and the risk of sinoatrial pauses and arrests (hcn4). Exposure to 10 or 25 µM ivabradine-an open channel blocker of HCNs-for 24 h resulted in a dose-dependent higher HRV and lower heart rate at 5 days post-fertilization. Hence, our screen confirmed the role of established genes for heart rate and rhythm (RGS6 and HCN4); showed that ivabradine reduces heart rate and increases HRV in zebrafish embryos, as it does in humans; and highlighted a novel gene that plays a role in HRV (KIAA1755).

Cyclin A2 is a key regulator of the cell cycle, implicated both in DNA replication and mitotic entry. Cyclin A2 participates in feedback loops that activate mitotic kinases in G2 phase, but why active Cyclin A2-CDK2 during the S phase does not trigger mitotic kinase activation remains unclear. Here, we describe a change in localisation of Cyclin A2 from being only nuclear to both nuclear and cytoplasmic at the S/G2 border. We find that Cyclin A2-CDK2 can activate the mitotic kinase PLK1 through phosphorylation of Bora, and that only cytoplasmic Cyclin A2 interacts with Bora and PLK1. Expression of predominately cytoplasmic Cyclin A2 or phospho-mimicking PLK1 T210D can partially rescue a G2 arrest caused by Cyclin A2 depletion. Cytoplasmic presence of Cyclin A2 is restricted by p21, in particular after DNA damage. Cyclin A2 chromatin association during DNA replication and additional mechanisms contribute to Cyclin A2 localisation change in the G2 phase. We find no evidence that such mechanisms involve G2 feedback loops and suggest that cytoplasmic appearance of Cyclin A2 at the S/G2 transition functions as a trigger for mitotic kinase activation.