Scientific highlights

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.

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).

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.

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.

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.

Differences between sexes contribute to variation in the levels of fasting glucose and insulin. Epidemiological studies established a higher prevalence of impaired fasting glucose in men and impaired glucose tolerance in women, however, the genetic component underlying this phenomenon is not established. We assess sex-dimorphic (73,089/50,404 women and 67,506/47,806 men) and sex-combined (151,188/105,056 individuals) fasting glucose/fasting insulin genetic effects via genome-wide association study meta-analyses in individuals of European descent without diabetes. Here we report sex dimorphism in allelic effects on fasting insulin at IRS1 and ZNF12 loci, the latter showing higher RNA expression in whole blood in women compared to men. We also observe sex-homogeneous effects on fasting glucose at seven novel loci. Fasting insulin in women shows stronger genetic correlations than in men with waist-to-hip ratio and anorexia nervosa. Furthermore, waist-to-hip ratio is causally related to insulin resistance in women, but not in men. These results position dissection of metabolic and glycemic health sex dimorphism as a steppingstone for understanding differences in genetic effects between women and men in related phenotypes.

Chromoanagenesis is a genomic event responsible for the formation of complex structural chromosomal rearrangements (CCRs). Germline chromoanagenesis is rare and the majority of reported cases are associated with an affected phenotype. Here, we report a healthy female carrying two de novo CCRs involving chromosomes 4, 19, 21 and X and chromosomes 7 and 11, respectively, with a total of 137 breakpoint junctions (BPJs). We characterized the CCRs using a hybrid-sequencing approach, combining short-read sequencing, nanopore sequencing, and optical mapping. The results were validated using multiple cytogenetic methods, including fluorescence in situ hybridization, spectral karyotyping, and Sanger sequencing. We identified 137 BPJs, which to our knowledge is the highest number of reported breakpoint junctions in germline chromoanagenesis. We also performed a statistical assessment of the positioning of the breakpoints, revealing a significant enrichment of BPJ-affecting genes (96 intragenic BPJs, 26 genes, p < 0.0001), indicating that the CCRs formed during active transcription of these genes. In addition, we find that the DNA fragments are unevenly and non-randomly distributed across the derivative chromosomes indicating a multistep process of scattering and re-joining of DNA fragments. In summary, we report a new maximum number of BPJs (137) in germline chromoanagenesis. We also show that a hybrid sequencing approach is necessary for the correct characterization of complex CCRs. Through in-depth statistical assessment, it was found that the CCRs most likely was formed through an event resembling chromoplexy-a catastrophic event caused by erroneous transcription factor binding.

Autonomous regulation of the intestine requires the combined activity of functionally distinct neurons of the enteric nervous system (ENS). However, the variety of enteric neuron types and how they emerge during development remain largely unknown. Here, we define a molecular taxonomy of 12 enteric neuron classes within the myenteric plexus of the mouse small intestine using single-cell RNA sequencing. We present cell-cell communication features and histochemical markers for motor neurons, sensory neurons and interneurons, together with transgenic tools for class-specific targeting. Transcriptome analysis of the embryonic ENS uncovers a novel principle of neuronal diversification, where two neuron classes arise through a binary neurogenic branching and all other identities emerge through subsequent postmitotic differentiation. We identify generic and class-specific transcriptional regulators and functionally connect Pbx3 to a postmitotic fate transition. Our results offer a conceptual and molecular resource for dissecting ENS circuits and predicting key regulators for directed differentiation of distinct enteric neuron classes.

Clear cell renal cell carcinoma (ccRCC) is the most common form of renal cancer. Due to inactivation of the von Hippel-Lindau tumour suppressor, the hypoxia-inducible transcription factors (HIFs) are constitutively activated in these tumours, resulting in a pseudo-hypoxic phenotype. The HIFs induce the expression of genes involved in angiogenesis and cell survival, but they also reset the cellular metabolism to protect cells from oxygen and nutrient deprivation. ccRCC tumours are highly vascularized and the cytoplasm of the cancer cells is filled with lipid droplets and glycogen, resulting in the histologically distinctive pale (clear) cytoplasm. Intratumoural heterogeneity may occur, and in some tumours, areas with granular, eosinophilic cytoplasm are found. Little is known regarding these traits and how they relate to the coexistent clear cell component, yet eosinophilic ccRCC is associated with higher grade and clinically more aggressive tumours. In this study, we have for the first time performed RNA sequencing comparing histologically verified clear cell and eosinophilic areas from ccRCC tissue, aiming to analyse the characteristics of these cell types. Findings from RNA sequencing were confirmed by immunohistochemical staining of biphasic ccRCC. We found that the eosinophilic phenotype displayed a higher proliferative drive and lower differentiation, and we confirmed a correlation to tumours of higher stage. We further identified mutations of the tumour suppressor p53 (TP53) exclusively in the eosinophilic ccRCC component, where mTORC1 activity was also elevated. Also, eosinophilic areas were less vascularized, yet harboured more abundant infiltrating immune cells. The cytoplasm of clear cell ccRCC cells was filled with lipids but had very low mitochondrial content, while the reverse was found in eosinophilic tissue. We herein suggest possible transcriptional mechanisms behind these phenomena. © 2020 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd. on behalf of The Pathological Society of Great Britain and Ireland.

Chrysanthemum virus B encodes a multifunctional p12 protein that acts as a transcriptional activator in the nucleus and as a suppressor of RNA silencing in the cytoplasm. Here, we investigated the impact of p12 on accumulation of major classes of small RNAs (sRNAs). The results show dramatic changes in the sRNA profiles characterised by an overall reduction in sRNA accumulation, changes in the pattern of size distribution of canonical siRNAs and in the ratio between sense and antisense strands, lower abundance of siRNAs with a U residue at the 5'-terminus, and changes in the expression of certain miRNAs, most of which were downregulated.

Genomic selection (GS) or genomic prediction is a promising approach for tree breeding to obtain higher genetic gains by shortening time of progeny testing in breeding programs. As proof-of-concept for Scots pine (Pinus sylvestris L.), a genomic prediction study was conducted with 694 individuals representing 183 full-sib families that were genotyped with genotyping-by-sequencing (GBS) and phenotyped for growth and wood quality traits. 8719 SNPs were used to compare different genomic with pedigree prediction models. Additionally, four prediction efficiency methods were used to evaluate the impact of genomic breeding value estimations by assigning diverse ratios of training and validation sets, as well as several subsets of SNP markers. Genomic Best Linear Unbiased Prediction (GBLUP) and Bayesian Ridge Regression (BRR) combined with expectation maximization (EM) imputation algorithm showed slightly higher prediction efficiencies than Pedigree Best Linear Unbiased Prediction (PBLUP) and Bayesian LASSO, with some exceptions. A subset of approximately 6000 SNP markers, was enough to provide similar prediction efficiencies as the full set of 8719 markers. Additionally, prediction efficiencies of genomic models were enough to achieve a higher selection response, that varied between 50-143% higher than the traditional pedigree-based selection. Although prediction efficiencies were similar for genomic and pedigree models, the relative selection response was doubled for genomic models by assuming that earlier selections can be done at the seedling stage, reducing the progeny testing time, thus shortening the breeding cycle length roughly by 50%.

We previously observed an unexpected fivefold (35 vs. 200 days) difference in the survival of respiratory chain complex III (CIII) deficient Bcs1l p.S78G mice between two congenic backgrounds. Here, we identify a spontaneous homoplasmic mtDNA variant (m.G14904A, mt-Cybp.D254N), affecting the CIII subunit cytochrome b (MT-CYB), in the background with short survival. We utilize maternal inheritance of mtDNA to confirm this as the causative variant and show that it further decreases the low CIII activity in Bcs1lp.S78G tissues to below survival threshold by 35 days of age. Molecular dynamics simulations predict D254N to restrict the flexibility of MT-CYB ef loop, potentially affecting RISP dynamics. In Rhodobacter cytochrome bc1 complex the equivalent substitution causes a kinetics defect with longer occupancy of RISP head domain towards the quinol oxidation site. These findings represent a unique case of spontaneous mitonuclear epistasis and highlight the role of mtDNA variation as modifier of mitochondrial disease phenotypes.

After over a decade of developments in field collection, laboratory methods and advances in high-throughput sequencing, contamination remains a key issue in ancient DNA research. Currently, human and microbial contaminant DNA still impose challenges on cost-effective sequencing and accurate interpretation of ancient DNA data. Here we investigate whether human contaminating DNA can be found in ancient faunal sequencing datasets. We identify variable levels of human contamination, which persists even after the sequence reads have been mapped to the faunal reference genomes. This contamination has the potential to affect a range of downstream analyses. We propose a fast and simple method, based on competitive mapping, which allows identifying and removing human contamination from ancient faunal DNA datasets with limited losses of true ancient data. This method could represent an important tool for the ancient DNA field.

Biological sulfate and thiosulfate reduction under haloalkaline conditions can be applied to treat waste streams from biodesulfurization systems. However, the lack of microbial aggregation under haloalkaline conditions limits the volumetric rates of sulfate and thiosulfate reducing bioreactors. As biomass retention in haloalkaline bioreactors has not been studied before, sand was chosen as a biomass carrier material to increase cell retention and consequently raise the volumetric rates. The results showed that ~10 fold higher biomass concentrations could be achieved with sand, compared to previous studies without carrier addition. The volumetric rates of sulfate/thiosulfate reduction increased approximately 4.5 times. Biomass attachment to the sand was restricted to cavities within the sand particles. Acetate produced by acetogenic bacteria from H 2 and CO2 was used as carbon source for biomass growth, while formate that was also produced from H2 and CO2 enhanced sulfate reduction. The microbial community composition was analyzed by 16S rRNA gene amplicon sequencing, and Tindallia related bacteria were probably responsible for formate formation from hydrogen. The community attached to the sand particles was similar to the suspended fraction, but the relative abundance of sequences most closely related to Desulfohalobiaceae was much higher in the attached fraction compared to the suspended fraction (30% and 13%, respectively). The results indicated that even though the biomass attachment to sand was poor, it still increased the biomass concentration and consequently the sulfate and thiosulfate reduction volumetric rates.