Scientific highlights

Barnacles are key marine crustaceans in several habitats, and they constitute a common practical problem by causing biofouling on man-made marine constructions and ships. Despite causing considerable ecological and economic impacts, there is a surprising void of basic genomic knowledge, and a barnacle reference genome is lacking. We here set out to characterize the genome of the bay barnacle Balanus improvisus (= Amphibalanus improvisus) based on short-read whole-genome sequencing and experimental genome size estimation. We show both experimentally (DNA staining and flow cytometry) and computationally (k-mer analysis) that B. improvisus has a haploid genome size of ~ 740 Mbp. A pilot genome assembly rendered a total assembly size of ~ 600 Mbp and was highly fragmented with an N50 of only 2.2 kbp. Further assembly-based and assembly-free analyses revealed that the very limited assembly contiguity is due to the B. improvisus genome having an extremely high nucleotide diversity (π) in coding regions (average π ≈ 5% and average π in fourfold degenerate sites ≈ 20%), and an overall high repeat content (at least 40%). We also report on high variation in the α-octopamine receptor OctA (average π = 3.6%), which might increase the risk that barnacle populations evolve resistance toward antifouling agents. The genomic features described here can help in planning for a future high-quality reference genome, which is urgently needed to properly explore and understand proteins of interest in barnacle biology and marine biotechnology and for developing better antifouling strategies.

Sex chromosomes have evolved independently multiple times in eukaryotes and are therefore considered a prime example of convergent genome evolution. Sex chromosomes are known to emerge after recombination is halted between a homologous pair of chromosomes, and this leads to a range of non-adaptive modifications causing gradual degeneration and gene loss on the sex-limited chromosome. However, the proximal causes of recombination suppression and the pace at which degeneration subsequently occurs remain unclear. Here, we use long- and short-read single-molecule sequencing approaches to assemble and annotate a draft genome of the basket willow, Salix viminalis, a species with a female heterogametic system at the earliest stages of sex chromosome emergence. Our single-molecule approach allowed us to phase the emerging Z and W haplotypes in a female, and we detected very low levels of Z/W single-nucleotide divergence in the non-recombining region. Linked-read sequencing of the same female and an additional male (ZZ) revealed the presence of two evolutionary strata supported by both divergence between the Z and W haplotypes and by haplotype phylogenetic trees. Gene order is still largely conserved between the Z and W homologs, although the W-linked region contains genes involved in cytokinin signaling regulation that are not syntenic with the Z homolog. Furthermore, we find no support across multiple lines of evidence for inversions, which have long been assumed to halt recombination between the sex chromosomes. Our data suggest that selection against recombination is a more gradual process at the earliest stages of sex chromosome formation than would be expected from an inversion and may result instead from the accumulation of transposable elements. Our results present a cohesive understanding of the earliest genomic consequences of recombination suppression as well as valuable insights into the initial stages of sex chromosome formation and regulation of sex differentiation.

Pathogen-mediated balancing selection is commonly considered to play an important role in the maintenance of genetic diversity, in particular in immune genes. However, the factors that may influence which immune genes are the targets of such selection are largely unknown. To address this, here we focus on Pattern Recognition Receptor (PRR) signalling pathways, which play a key role in innate immunity. We used whole-genome resequencing data from a population of bank voles (Myodes glareolus) to test for associations between balancing selection, pleiotropy and gene function in a set of 123 PRR signalling pathway genes. To investigate the effect of gene function, we compared genes encoding (a) receptors for microbial ligands versus downstream signalling proteins, and (b) receptors recognizing components of microbial cell walls, flagella and capsids versus receptors recognizing features of microbial nucleic acids. Analyses based on the nucleotide diversity of full coding sequences showed that balancing selection primarily targeted receptor genes with a low degree of pleiotropy. Moreover, genes encoding receptors recognizing components of microbial cell walls etc. were more important targets of balancing selection than receptors recognizing nucleic acids. Tests for localized signatures of balancing selection in coding and noncoding sequences showed that such signatures were mostly located in introns, and more evenly distributed among different functional categories of PRR pathway genes. The finding that signatures of balancing selection in full coding sequences primarily occur in receptor genes, in particular those encoding receptors for components of microbial cell walls etc., is consistent with the idea that coevolution between hosts and pathogens is an important cause of balancing selection on immune genes.

To investigate associations between a high genetic disease risk and disease severity in patients with systemic lupus erythematosus (SLE). Patients with SLE (n=1001, discovery cohort and n=5524, replication cohort) and healthy controls (n=2802 and n=9859) were genotyped using a 200K Immunochip single nucleotide polymorphism array. A genetic risk score (GRS) was assigned to each individual based on 57 SLE risk loci. SLE was more prevalent in the high, compared with the low, GRS-quartile (OR 12.32 (9.53 to 15.71), p=7.9×10-86 and OR 7.48 (6.73 to 8.32), p=2.2×10-304 for the discovery and the replication cohorts, respectively). In the discovery cohort, patients in the high GRS-quartile had a 6-year earlier mean disease onset (HR 1.47 (1.22 to 1.75), p=4.3×10-5), displayed higher prevalence of damage accrual (OR 1.47 (1.06 to 2.04), p=2.0×10-2), renal disorder (OR 2.22 (1.50 to 3.27), p=5.9×10-5), anti-dsDNA (OR 1.83 (1.19 to 2.81), p=6.1×10-3), end-stage renal disease (ESRD) (OR 5.58 (1.50 to 20.79), p=1.0×10-2), proliferative nephritis (OR 2.42 (1.30 to 4.49), p=5.1×10-3), anti-cardiolipin-IgG (OR 1.89 (1.13 to 3.18), p=1.6×10-2), anti-β2-glycoprotein-I-IgG (OR 2.29 (1.29 to 4.06), p=4.8×10-3) and positive lupus anticoagulant test (OR 2.12 (1.16 to 3.89), p=1.5×10-2) compared with patients in the low GRS-quartile. Survival analysis showed earlier onset of the first organ damage (HR 1.51 (1.04 to 2.25), p=3.7×10-2), first cardiovascular event (HR 1.65 (1.03 to 2.64), p=2.6×10-2), nephritis (HR 2.53 (1.72 to 3.71), p=9.6×10-7), ESRD (HR 6.78 (1.78 to 26.86), p=6.5×10-3) and decreased overall survival (HR 1.83 (1.02 to 3.30), p=4.3×10-2) in high to low quartile comparison. A high GRS is associated with increased risk of organ damage, renal dysfunction and all-cause mortality. Our results indicate that genetic profiling may be useful for predicting outcomes in patients with SLE.

Extracellular vesicles (EVs) mediate cell-to-cell communication by delivering or displaying macromolecules to their recipient cells. While certain broad-spectrum EV effects reflect their protein cargo composition, others have been attributed to individual EV-loaded molecules such as specific miRNAs. In this work, we have investigated the contents of vesicular cargo using small RNA sequencing of cells and EVs from HEK293T, RD4, C2C12, Neuro2a and C17.2. The majority of RNA content in EVs (49-96%) corresponded to rRNA-, coding- and tRNA fragments, corroborating with our proteomic analysis of HEK293T and C2C12 EVs which showed an enrichment of ribosome and translation-related proteins. On the other hand, the overall proportion of vesicular small RNA was relatively low and variable (2-39%) and mostly comprised of miRNAs and sequences mapping to piRNA loci. Importantly, this is one of the few studies, which systematically links vesicular RNA and protein cargo of vesicles. Our data is particularly useful for future work in unravelling the biological mechanisms underlying vesicular RNA and protein sorting and serves as an important guide in developing EVs as carriers for RNA therapeutics.

The FADS1 and FADS2 genes in the FADS cluster encode the rate-limiting enzymes in the synthesis of long-chain polyunsaturated fatty acids (LC-PUFAs). Genetic variation in this region has been associated with a large number of diseases and traits many of them correlated to differences in metabolism of PUFAs. However, the causative variants leading to these associations have not been identified. Here we find that the multiallelic rs174557 located in an AluYe5 element in intron 1 of FADS1 is functional and lies within a PATZ1 binding site. The derived allele of rs174557, which is the common variant in most populations, diminishes binding of PATZ1, a transcription factor conferring allele-specific downregulation of FADS1. The PATZ1 binding site overlaps with a SP1 site. The competitive binding between the suppressive PATZ1 and the activating complex of SP1 and SREBP1c determines the enhancer activity of this region, which regulates expression of FADS1.

Homo sapiens evolved hundreds of thousands of years ago in Africa, later spreading across the globe1, but the early evolutionary process is debated2-6. Here we present whole-genome sequencing data for 28 ancient southern African individuals, including six individuals with 25× to 7.2× genome coverage, dated to between 10,200 and 150 calibrated years before present (cal. BP). All ancient southern Africans dated to more than 1,400 cal. BP show a genetic make-up that is outside the range of genetic variation in modern-day humans (including southern African Khoe-San people, although some retain up to 80% ancient southern African ancestry), manifesting in a large fraction of Homo sapiens-specific variants that are unique to ancient southern Africans. Homo sapiens-specific variants at amino acid-altering sites fixed for all humans-which are likely to have evolved rapidly on the Homo sapiens branch-were enriched for genes associated with kidney function. Some Homo sapiens-specific variants fixed in ancient southern Africans-which are likely to have adapted rapidly on the southern African branch-were enriched for genes associated with protection against ultraviolet light. The ancient southern Africans show little spatiotemporal stratification for 9,000 years, consistent with a large, stable Holocene population transcending archaeological phases. While southern Africa served as a long-standing geographical refugium, there is outward gene flow over 8,000 years ago; however, inward gene flow manifests only after around 1,400 years ago. The ancient genomes reported here are therefore key to the evolution of Homo sapiens, and are important for advancing our understanding of human genomic variation.

B vitamins are essential cofactors for practically all living organisms on Earth and are produced by a selection of microorganisms. An imbalance between high demand and limited production, in concert with abiotic processes, may explain the low availability of these vitamins in marine systems. Natural microbial communities from surface shelf water in the productive area off northwestern Spain were enclosed in mesocosms in winter, spring, and summer 2016. In order to explore the impact of B-vitamin availability on microbial community composition (16S and 18S rRNA gene sequence analysis) and bacterial function (metatranscriptomics analysis) in different seasons, enrichment experiments were conducted with seawater from the mesocosms. Our findings revealed that significant increases in phytoplankton or prokaryote biomass associated with vitamin B12 and/or B1 amendments were not accompanied by significant changes in community composition, suggesting that most of the microbial taxa benefited from the external B-vitamin supply. Metatranscriptome analysis suggested that many bacteria were potential consumers of vitamins B12 and B1, although the relative abundance of reads related to synthesis was ca. 3.6-fold higher than that related to uptake. Alteromonadales and Oceanospirillales accounted for important portions of vitamin B1 and B12 synthesis gene transcription, despite accounting for only minor portions of the bacterial community. Flavobacteriales appeared to be involved mostly in vitamin B12 and B1 uptake, and Pelagibacterales expressed genes involved in vitamin B1 uptake. Interestingly, the relative expression of vitamin B12 and B1 synthesis genes among bacteria strongly increased upon inorganic nutrient amendment. Collectively, these findings suggest that upwelling events intermittently occurring during spring and summer in productive ecosystems may ensure an adequate production of these cofactors to sustain high levels of phytoplankton growth and biomass. IMPORTANCE B vitamins are essential growth factors for practically all living organisms on Earth that are produced by a selection of microorganisms. An imbalance between high demand and limited production may explain the low concentration of these compounds in marine systems. In order to explore the impact of B-vitamin availability on bacteria and algae in the coastal waters off northwestern Spain, six experiments were conducted with natural surface water enclosed in winter, spring, and summer. Our findings revealed that increases in phytoplankton or bacterial growth associated with B12 and/or B1 amendments were not accompanied by significant changes in community composition, suggesting that most microorganisms benefited from the B-vitamin supply. Our analyses confirmed the role of many bacteria as consumers of vitamins B12 and B1, although the relative abundance of genes related to synthesis was ca. 3.6-fold higher than that related to uptake. Interestingly, prokaryote expression of B12 and B1 synthesis genes strongly increased when inorganic nutrients were added. Collectively, these findings suggest that upwelling of cold and nutrient-rich waters occurring during spring and summer in this coastal area may ensure an adequate production of B vitamins to sustain high levels of algae growth and biomass.

The social organization of the first fully sedentary societies that emerged during the Neolithic period in Southwest Asia remains enigmatic,1 mainly because material culture studies provide limited insight into this issue. However, because Neolithic Anatolian communities often buried their dead beneath domestic buildings,2 household composition and social structure can be studied through these human remains. Here, we describe genetic relatedness among co-burials associated with domestic buildings in Neolithic Anatolia using 59 ancient genomes, including 22 new genomes from Aşıklı Höyük and Çatalhöyük. We infer pedigree relationships by simultaneously analyzing multiple types of information, including autosomal and X chromosome kinship coefficients, maternal markers, and radiocarbon dating. In two early Neolithic villages dating to the 9th and 8th millennia BCE, Aşıklı Höyük and Boncuklu, we discover that siblings and parent-offspring pairings were frequent within domestic structures, which provides the first direct indication of close genetic relationships among co-burials. In contrast, in the 7th millennium BCE sites of Çatalhöyük and Barcın, where we study subadults interred within and around houses, we find close genetic relatives to be rare. Hence, genetic relatedness may not have played a major role in the choice of burial location at these latter two sites, at least for subadults. This supports the hypothesis that in Çatalhöyük,3-5 and possibly in some other Neolithic communities, domestic structures may have served as burial location for social units incorporating biologically unrelated individuals. Our results underscore the diversity of kin structures in Neolithic communities during this important phase of sociocultural development.

Inflammation plays an important role in cardiovascular disease (CVD) development. The NOD-like receptor protein-3 (NLRP3) inflammasome contributes to the development of atherosclerosis in animal models. Components of the NLRP3 inflammasome pathway such as interleukin-1β can therapeutically be targeted. Associations of genetically determined inflammasome-mediated systemic inflammation with CVD and mortality in humans are unknown. We explored the association of genetic NLRP3 variants with prevalent CVD and cardiovascular mortality in 538 167 subjects on the individual participant level in an explorative gene-centric approach without performing multiple testing. Functional relevance of single-nucleotide polymorphisms on NLRP3 inflammasome activation has been evaluated in monocyte-enriched peripheral blood mononuclear cells (PBMCs). Genetic analyses identified the highly prevalent (minor allele frequency 39.9%) intronic NLRP3 variant rs10754555 to affect NLRP3 gene expression. rs10754555 carriers showed significantly higher C-reactive protein and serum amyloid A plasma levels. Carriers of the G allele showed higher NLRP3 inflammasome activation in isolated human PBMCs. In carriers of the rs10754555 variant, the prevalence of coronary artery disease was significantly higher as compared to non-carriers with a significant interaction between rs10754555 and age. Importantly, rs10754555 carriers had significantly higher risk for cardiovascular mortality during follow-up. Inflammasome inducers (e.g. urate, triglycerides, apolipoprotein C3) modulated the association between rs10754555 and mortality. The NLRP3 intronic variant rs10754555 is associated with increased systemic inflammation, inflammasome activation, prevalent coronary artery disease, and mortality. This study provides evidence for a substantial role of genetically driven systemic inflammation in CVD and highlights the NLRP3 inflammasome as a therapeutic target.

With the expanding ash dieback epidemic that has spread across the European continent, an improved functional understanding of the disease development in afflicted hosts is needed. The study investigated whether differences in necrosis extension between common ash (Fraxinus excelsior) trees with different levels of susceptibility to the fungus Hymenoscyphus fraxineus are associated with, and can be explained by, the differences in gene expression patterns. We inoculated seemingly healthy branches of each of two resistant and susceptible ash genotypes with H. fraxineus grown in a common garden. Ten months after the inoculation, the length of necrosis on the resistant genotypes were shorter than on the susceptible genotypes. RNA sequencing of bark samples collected at the border of necrotic lesions and from healthy tissues distal to the lesion revealed relatively limited differences in gene expression patterns between susceptible and resistant genotypes. At the necrosis front, only 138 transcripts were differentially expressed between the genotype categories while 1082 were differentially expressed in distal, non-symptomatic tissues. Among these differentially expressed genes, several genes in the mevalonate (MVA) and iridoid pathways were found to be co-regulated, possibly indicating increased fluxes through these pathways in response to H. fraxineus. Comparison of transcriptional responses of symptomatic and non-symptomatic ash in a controlled greenhouse experiment revealed a relatively small set of genes that were differentially and concordantly expressed in both studies. This gene-set included the rate-limiting enzyme in the MVA pathway and a number of transcription factors. Furthermore, several of the concordantly expressed candidate genes show significant similarity to genes encoding players in the abscisic acid- or Jasmonate-signalling pathways. A set of candidate genes, concordantly expressed between field and greenhouse experiments, was identified. The candidates are associated with hormone signalling and specialized metabolite biosynthesis pathways indicating the involvement of these pathways in the response of the host to infection by H. fraxineus.

In this work we aim to investigate the origins and genetic affinities of Bronze Age populations (2,400-1,100 BC) from the region of southern Poland and to trace maternal kinship patterns present in the burials of those populations by the use of complete mitochondrial genomes. We performed ancient DNA analyses for Bronze Age individuals from present-day Poland associated with the Strzyżow culture, the Mierzanowice culture, and the Trzciniec Cultural circle. To obtain complete mitochondrial genomes, we sequenced genomic libraries using Illumina platform. Additionally, hybridization capture was used to enrich some of the samples for mitochondrial DNA. AMS 14 C-dating was conducted for 51 individuals to verify chronological and cultural attribution of the analyzed samples. Complete ancient mitochondrial genomes were generated for 80 of the Bronze Age individuals from present-day Poland. The results of the population genetic analyses indicate close maternal genetic affinity between Mierzanowice, Trzciniec, and Corded Ware culture-associated populations. This is in contrast to the genetically more distant Strzyżów people that displayed closer maternal genetic relation to steppe populations associated with the preceding Yamnaya culture and Catacomb culture, and with later Scythians. Potential maternal kinship relations were identified in burials of Mierzanowice and Trzciniec populations analyzed in this study. Results revealed genetic continuity from the Late Neolithic Corded Ware groups to Bronze Age Mierzanowice and Trzciniec-associated populations, and possible additional genetic contribution from the steppe to the formation of the Strzyżów-associated group at the end of 3rd millennium BC. Mitochondrial patterns indicated several pairs of potentially maternally related individuals mostly in Trzciniec-associated group.

Cereal grains, including wheat (Triticum aestivum L.), are major sources of food and feed, with wheat being dominant in temperate zones. These end uses exploit the storage reserves in the starchy endosperm of the grain, with starch being the major storage component in most cereal species. However, oats (Avena sativa L.) differs in that the starchy endosperm stores significant amounts of oil. Understanding the control of carbon allocation between groups of storage compounds, such as starch and oil, is therefore important for understanding the composition and hence end use quality of cereals. WRINKLED1 is a transcription factor known to induce triacylglycerol (TAG; oil) accumulation in several plant storage tissues. An oat endosperm homolog of WRI1 (AsWRI1) expressed from the endosperm-specific HMW1Dx5 promoter resulted in drastic changes in carbon allocation in wheat grains, with reduced seed weight and a wrinkled seed phenotype. The starch content of mature grain endosperms of AsWRI1-wheat was reduced compared to controls (from 62 to 22% by dry weight (dw)), TAG was increased by up to nine-fold (from 0.7 to 6.4% oil by dw) and sucrose from 1.5 to 10% by dw. Expression of AsWRI1 in wheat grains also resulted in multiple layers of elongated peripheral aleurone cells. RNA-sequencing, lipid analyses, and pulse-chase experiments using 14C-sucrose indicated that futile cycling of fatty acids could be a limitation for oil accumulation. Our data show that expression of oat endosperm WRI1 in the wheat endosperm results in changes in metabolism which could underpin the application of biotechnology to manipulate grain composition. In particular, the striking effect on starch synthesis in the wheat endosperm indicates that an important indirect role of WRI1 is to divert carbon allocation away from starch biosynthesis in plant storage tissues that accumulate oil.

The electrocardiographic PR interval reflects atrioventricular conduction, and is associated with conduction abnormalities, pacemaker implantation, atrial fibrillation (AF), and cardiovascular mortality. Here we report a multi-ancestry (N = 293,051) genome-wide association meta-analysis for the PR interval, discovering 202 loci of which 141 have not previously been reported. Variants at identified loci increase the percentage of heritability explained, from 33.5% to 62.6%. We observe enrichment for cardiac muscle developmental/contractile and cytoskeletal genes, highlighting key regulation processes for atrioventricular conduction. Additionally, 8 loci not previously reported harbor genes underlying inherited arrhythmic syndromes and/or cardiomyopathies suggesting a role for these genes in cardiovascular pathology in the general population. We show that polygenic predisposition to PR interval duration is an endophenotype for cardiovascular disease, including distal conduction disease, AF, and atrioventricular pre-excitation. These findings advance our understanding of the polygenic basis of cardiac conduction, and the genetic relationship between PR interval duration and cardiovascular disease.

Current extinction rates are comparable to five prior mass extinctions in the earth's history, and are strongly affected by human activities that have modified more than half of the earth's terrestrial surface. Increasing human activity restricts animal movements and isolates formerly connected populations, a particular concern for the conservation of large carnivores, but no prior research has used high throughput sequencing in a standardized manner to examine genetic connectivity for multiple species of large carnivores and multiple ecosystems. Here, we used RAD SNP genotypes to test for differences in connectivity between multiple ecosystems for African wild dogs (Lycaon pictus) and lions (Panthera leo), and to test correlations between genetic distance, geographic distance and landscape resistance due to human activity. We found weaker connectivity, a stronger correlation between genetic distance and geographic distance, and a stronger correlation between genetic distance and landscape resistance for lions than for wild dogs, and propose a new hypothesis that adaptations to interspecific competition may help to explain differences in vulnerability to isolation by humans.

Mutation and recombination are key evolutionary processes governing phenotypic variation and reproductive isolation. We here demonstrate that biodiversity within all globally known strains of Schizosaccharomyces pombe arose through admixture between two divergent ancestral lineages. Initial hybridization was inferred to have occurred ∼20-60 sexual outcrossing generations ago consistent with recent, human-induced migration at the onset of intensified transcontinental trade. Species-wide heritable phenotypic variation was explained near-exclusively by strain-specific arrangements of alternating ancestry components with evidence for transgressive segregation. Reproductive compatibility between strains was likewise predicted by the degree of shared ancestry. To assess the genetic determinants of ancestry block distribution across the genome, we characterized the type, frequency, and position of structural genomic variation using nanopore and single-molecule real-time sequencing. Despite being associated with double-strand break initiation points, over 800 segregating structural variants exerted overall little influence on the introgression landscape or on reproductive compatibility between strains. In contrast, we found strong ancestry disequilibrium consistent with negative epistatic selection shaping genomic ancestry combinations during the course of hybridization. This study provides a detailed, experimentally tractable example that genomes of natural populations are mosaics reflecting different evolutionary histories. Exploiting genome-wide heterogeneity in the history of ancestral recombination and lineage-specific mutations sheds new light on the population history of S. pombe and highlights the importance of hybridization as a creative force in generating biodiversity.