This study aimed to characterise the RNA microbiome, including the virome of extended semen from Swedish breeding boars, with particular focus on Atypical porcine pestivirus (APPV). This neurotropic virus, associated with congenital tremor type A-II in piglets, was recently demonstrated to induce the disease through insemination with semen from infected boars. From 124 Artificial Insemination (AI) doses from Swedish breeding boars, APPV was detected in one dose in addition to a sparse seminal RNA virome, characterised by retroviruses, phages, and some fecal-associated contaminants. The detected seminal microbiome was large and characterized by Gram-negative bacteria from the phylum Proteobacteria, mainly consisting of apathogenic or opportunistic bacteria. The proportion of bacteria with a pathogenic potential was low, and no antimicrobial resistance genes (ARGs) were detected in the datasets. Overall, the results indicate a good health status among Swedish breeding boars. The detection of APPV in semen raises the question of whether routine screening for APPV in breeding boars should be instigated.

Unravelling the genomic landscape of divergence between lineages is key to understanding speciation. The naturally hybridizing collared flycatcher and pied flycatcher are important avian speciation models that show pre- as well as postzygotic isolation. We sequenced and assembled the 1.1-Gb flycatcher genome, physically mapped the assembly to chromosomes using a low-density linkage map and re-sequenced population samples of each species. Here we show that the genomic landscape of species differentiation is highly heterogeneous with approximately 50 'divergence islands' showing up to 50-fold higher sequence divergence than the genomic background. These non-randomly distributed islands, with between one and three regions of elevated divergence per chromosome irrespective of chromosome size, are characterized by reduced levels of nucleotide diversity, skewed allele-frequency spectra, elevated levels of linkage disequilibrium and reduced proportions of shared polymorphisms in both species, indicative of parallel episodes of selection. Proximity of divergence peaks to genomic regions resistant to sequence assembly, potentially including centromeres and telomeres, indicate that complex repeat structures may drive species divergence. A much higher background level of species divergence of the Z chromosome, and a lower proportion of shared polymorphisms, indicate that sex chromosomes and autosomes are at different stages of speciation. This study provides a roadmap to the emerging field of speciation genomics.

Arsenic, a carcinogen with immunotoxic effects, is a common contaminant of drinking water and certain food worldwide. We hypothesized that chronic arsenic exposure alters gene expression, potentially by altering DNA methylation of genes encoding central components of the immune system. We therefore analyzed the transcriptomes (by RNA sequencing) and methylomes (by target-enrichment next-generation sequencing) of primary CD4-positive T cells from matched groups of four women each in the Argentinean Andes, with fivefold differences in urinary arsenic concentrations (median concentrations of urinary arsenic in the lower- and high-arsenic groups: 65 and 276 μg/l, respectively). Arsenic exposure was associated with genome-wide alterations of gene expression; principal component analysis indicated that the exposure explained 53% of the variance in gene expression among the top variable genes and 19% of 28,351 genes were differentially expressed (false discovery rate <0.05) between the exposure groups. Key genes regulating the immune system, such as tumor necrosis factor alpha and interferon gamma, as well as genes related to the NF-kappa-beta complex, were significantly downregulated in the high-arsenic group. Arsenic exposure was associated with genome-wide DNA methylation; the high-arsenic group had 3% points higher genome-wide full methylation (>80% methylation) than the lower-arsenic group. Differentially methylated regions that were hyper-methylated in the high-arsenic group showed enrichment for immune-related gene ontologies that constitute the basic functions of CD4-positive T cells, such as isotype switching and lymphocyte activation and differentiation. In conclusion, chronic arsenic exposure from drinking water was related to changes in the transcriptome and methylome of CD4-positive T cells, both genome wide and in specific genes, supporting the hypothesis that arsenic causes immunotoxicity by interfering with gene expression and regulation.

Genome scans in recently separated species can inform on molecular mechanisms and evolutionary processes driving divergence. Large-scale polymorphism data from multiple species pairs are also key to investigate the repeatability of divergence-whether radiations tend to show parallel responses to similar selection pressures and/or underlying molecular forces. Here, we used whole-genome resequencing data from six wood white (Leptidea sp.) butterfly populations, representing three closely related species with karyomorph variation, to infer the species' demographic history and characterize patterns of genomic diversity and differentiation. The analyses supported previously established species relationships, and there was no evidence for postdivergence gene flow. We identified significant intraspecific genetic structure, in particular between karyomorph extremes in the wood white (L. sinapis)-a species with a remarkable chromosome number cline across the distribution range. The genomic landscapes of differentiation were erratic, and outlier regions were narrow and dispersed. Highly differentiated (F ST ) regions generally had low genetic diversity (θπ ), but increased absolute divergence (DXY ) and excess of rare frequency variants (low Tajima's D). A minority of differentiation peaks were shared across species and population comparisons. However, highly differentiated regions contained genes with overrepresented functions related to metabolism, response to stimulus and cellular processes, indicating recurrent directional selection on a specific set of traits in all comparisons. In contrast to the majority of genome scans in recently diverged lineages, our data suggest that divergence landscapes in Leptidea have been shaped by directional selection and genetic drift rather than stable recombination landscapes and/or introgression.

Pseudozyma hubeiensis is a basidiomycete yeast that has the highly desirable traits for lignocellulose valorisation of being equally efficient at utilization of glucose and xylose, and capable of their co-utilization. The species has previously mainly been studied for its capacity to produce secreted biosurfactants in the form of mannosylerythritol lipids, but it is also an oleaginous species capable of accumulating high levels of triacylglycerol storage lipids during nutrient starvation. In this study, we aimed to further characterize the oleaginous nature of P. hubeiensis by evaluating metabolism and gene expression responses during storage lipid formation conditions with glucose or xylose as a carbon source. The genome of the recently isolated P. hubeiensis BOT-O strain was sequenced using MinION long-read sequencing and resulted in the most contiguous P. hubeiensis assembly to date with 18.95 Mb in 31 contigs. Using transcriptome data as experimental support, we generated the first mRNA-supported P. hubeiensis genome annotation and identified 6540 genes. 80% of the predicted genes were assigned functional annotations based on protein homology to other yeasts. Based on the annotation, key metabolic pathways in BOT-O were reconstructed, including pathways for storage lipids, mannosylerythritol lipids and xylose assimilation. BOT-O was confirmed to consume glucose and xylose at equal rates, but during mixed glucose-xylose cultivation glucose was found to be taken up faster. Differential expression analysis revealed that only a total of 122 genes were significantly differentially expressed at a cut-off of |log2 fold change| ≥ 2 when comparing cultivation on xylose with glucose, during exponential growth and during nitrogen-starvation. Of these 122 genes, a core-set of 24 genes was identified that were differentially expressed at all time points. Nitrogen-starvation resulted in a larger transcriptional effect, with a total of 1179 genes with significant expression changes at the designated fold change cut-off compared with exponential growth on either glucose or xylose.

Cell replacement therapies for neurodegenerative disease have focused on transplantation of the cell types affected by the pathological process. Here we describe an alternative strategy for Parkinson's disease in which dopamine neurons are generated by direct conversion of astrocytes. Using three transcription factors, NEUROD1, ASCL1 and LMX1A, and the microRNA miR218, collectively designated NeAL218, we reprogram human astrocytes in vitro, and mouse astrocytes in vivo, into induced dopamine neurons (iDANs). Reprogramming efficiency in vitro is improved by small molecules that promote chromatin remodeling and activate the TGFβ, Shh and Wnt signaling pathways. The reprogramming efficiency of human astrocytes reaches up to 16%, resulting in iDANs with appropriate midbrain markers and excitability. In a mouse model of Parkinson's disease, NeAL218 alone reprograms adult striatal astrocytes into iDANs that are excitable and correct some aspects of motor behavior in vivo, including gait impairments. With further optimization, this approach may enable clinical therapies for Parkinson's disease by delivery of genes rather than cells.

Germline genetic variants are an important cause of dilated cardiomyopathy (DCM). However, recent sequencing studies have revealed rare variants in DCM-associated genes also in individuals without known heart disease. In this study, we investigate variant prevalence and genotype-phenotype correlations in Swedish DCM patients, and compare their genetic variants to those detected in reference cohorts. We sequenced the coding regions of 41 DCM-associated genes in 176 unrelated patients with idiopathic DCM and found 102 protein-altering variants with an allele frequency of <0.04% in reference cohorts; the majority were missense variants not previously described in DCM. Fifty-five (31%) patients had one variant, and 24 (14%) patients had two or more variants in the analysed genes. Detection of genetic variants in any gene, and in LMNA, MYH7 or TTN alone, was associated with early onset disease and reduced transplant-free survival. As expected, nonsense and frameshift variants were more common in DCM patients than in healthy individuals of the reference cohort 1000 Genomes Europeans. Surprisingly however, the prevalence, conservation and pathogenicity scores, and localization of missense variants were similar in DCM patients and healthy reference individuals. To our knowledge, this is the first study to identify correlations between genotype and prognosis when sequencing a large number of genes in unselected DCM patients. The similar distribution of missense variants in DCM patients and healthy reference individuals questions the pathogenic role of many variants, and suggests that results from genetic testing of DCM patients should be interpreted with caution.